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

PRDM16 determines specification of ventricular cardiomyocytes by suppressing alternative cell fates

PRDM16 is a transcription factor with histone methyltransferase activity expressed at the earliest stages of cardiac development. Pathogenic mutations in humans lead to cardiomyopathy, conduction abnormalities, and heart failure. PRDM16 is specifically expressed in ventricular but not atrial cardiomyocytes, and its expression declines postnatally. Because in other tissues PRDM16 is best known for its role in binary cell fate decisions, we hypothesized a similar decision-making function in cardiomyocytes. Here, we demonstrated that cardiomyocyte-specific deletion of Prdm16 during cardiac development results in contractile dysfunction and abnormal electrophysiology of the postnatal heart, resulting in premature death. By combined RNA+ATAC single-cell sequencing, we found that PRDM16 favors ventricular working cardiomyocyte identity, by opposing the activity of master regulators of ventricular conduction and atrial fate. Myocardial loss of PRDM16 during development resulted in hyperplasia of the (distal) ventricular conduction system. Hence, PRDM16 plays an indispensable role during cardiac development by driving ventricular working cardiomyocyte identity.

scBlood: A comprehensive single-cell accessible chromatin database of blood cells

The advent of single cell transposase-accessible chromatin sequencing (scATAC-seq) technology enables us to explore the genomic characteristics and chromatin accessibility of blood cells at the single-cell level. To fully make sense of the roles and regulatory complexities of blood cells, it is critical to collect and analyze these rapidly accumulating scATAC-seq datasets at a system level. Here, we present scBlood (https://bio.liclab.net/scBlood/), a comprehensive single-cell accessible chromatin database of blood cells. The current version of scBlood catalogs 770,907 blood cells and 452,247 non-blood cells from ∼400 high-quality scATAC-seq samples covering 30 tissues and 21 disease types. All data hosted on scBlood have undergone preprocessing from raw fastq files and multiple standards of quality control. Furthermore, we conducted comprehensive downstream analyses, including multi-sample integration analysis, cell clustering and annotation, differential chromatin accessibility analysis, functional enrichment analysis, co-accessibility analysis, gene activity score calculation, and transcription factor (TF) enrichment analysis. In summary, scBlood provides a user-friendly interface for searching, browsing, analyzing, visualizing, and downloading scATAC-seq data of interest. This platform facilitates insights into the functions and regulatory mechanisms of blood cells, as well as their involvement in blood-related diseases.

Antiviral innate immune memory in alveolar macrophages following SARS-CoV-2 infection ameliorates secondary influenza A virus disease

Pathogen encounter can result in epigenetic remodeling that shapes disease caused by heterologous pathogens. Here, we examined innate immune memory in the context of commonly circulating respiratory viruses. Single-cell analyses of airway-resident immune cells in a disease-relevant murine model of SARS-CoV-2 recovery revealed epigenetic reprogramming in alveolar macrophages following infection. Post-COVID-19 human monocytes exhibited similar epigenetic signatures. In airway-resident macrophages, past SARS-CoV-2 infection increased activity of type I interferon (IFN-I)-related transcription factors and epigenetic poising of antiviral genes. Viral pattern recognition and canonical IFN-I signaling were required for the establishment of this innate immune memory and augmented secondary antiviral responses. Antiviral innate immune memory mounted by airway-resident macrophages post-SARS-CoV-2 was necessary and sufficient to ameliorate secondary disease caused by influenza A virus and curtailed hyperinflammatory dysregulation and mortality. Our findings provide insights into antiviral innate immune memory in the airway that may facilitate the development of broadly effective therapeutic strategies.

scImmOmics: a manually curated resource of single-cell multi-omics immune data

Single-cell sequencing technology has enabled the discovery and characterization of subpopulations of immune cells with unique functions, which is critical for revealing immune responses under healthy or disease conditions. Efforts have been made to collect and curate single-cell RNA sequencing (scRNA-seq) data, yet an immune-specific single-cell multi-omics atlas with harmonized metadata is still lacking. Here, we present scImmOmics (https://bio.liclab.net/scImmOmics/home), a manually curated single-cell multi-omics immune database constructed based on high-quality immune cells with known immune cell labels. Currently, scImmOmics documents >2.9 million cell-type labeled immune cells derived from seven single-cell sequencing technologies, involving 131 immune cell types, 47 tissues and 4 species. To ensure data consistency, we standardized the nomenclature of immune cell types and presented them in a hierarchical tree structure to clearly describe the lineage relationships within the immune system. scImmOmics also provides comprehensive immune regulatory information, including T-cell/B-cell receptor sequencing clonotype information, cell-specific regulatory information (e.g. gene/chromatin accessibility/protein/transcription factor states within known cell types, cell-to-cell communication and co-expression networks) and immune cell responses to cytokines. Collectively, scImmOmics is a comprehensive and valuable platform for unraveling the heterogeneity and diversity of immune cells and elucidating the specific regulatory mechanisms at the single-cell level.

Comparative single-cell analyses identify shared and divergent features of human and mouse kidney development

The mammalian kidney maintains fluid homeostasis through diverse epithelial cell types generated from nephron and ureteric progenitor cells. To extend a developmental understanding of the kidney's epithelial networks, we compared chromatin organization (single-nuclear assay for transposase-accessible chromatin sequencing [ATAC-seq]; 112,864 nuclei) and gene expression (single-cell/nuclear RNA sequencing [RNA-seq]; 109,477 cells/nuclei) in the developing human (10.6-17.6 weeks; n = 10) and mouse (post-natal day [P]0; n = 10) kidney, supplementing analysis with published mouse datasets from earlier stages. Single-cell/nuclear datasets were analyzed at a species level, and then nephron and ureteric cellular lineages were extracted and integrated into a common, cross-species, multimodal dataset. Comparative computational analyses identified conserved and divergent features of chromatin organization and linked gene activity, identifying species-specific and cell-type-specific regulatory programs. In situ validation of human-enriched gene activity points to human-specific signaling interactions in kidney development. Further, human-specific enhancer regions were linked to kidney diseases through genome-wide association studies (GWASs), highlighting the potential for clinical insight from developmental modeling.

Comparative analysis of tongue cancer organoids among patients identifies the heritable nature of minimal residual disease

The relapse of tongue cancer (TC) after chemotherapy is caused by minimal residual disease (MRD), which is a few remaining cancer cells after chemotherapy. To understand the mechanism of MRD in TC, we created a library of TC organoids (TCOs) from 28 untreated TC patients at diverse ages and cancer stages. These TCOs reproduced the primary TC tissues both in vitro and in a xenograft model, and several TCO lines survived after cisplatin treatment (chemo-resistant TCOs). Of note, the chemo-resistant TCOs showed "heritable" embryonic diapause-like features before treatment and activation of the autophagy and cholesterol biosynthetic pathways. Importantly, inhibiting these pathways with specific inhibitors converted the chemo-resistant TCOs into chemo-sensitive TCOs. Conversely, autophagy activation with mTOR inhibitors conferred chemo-resistance on the chemo-sensitive TCOs. This unique model provides insights into the mechanism of MRD formation in TCs, leading to effective therapeutic approaches to reduce the recurrence of TC.

Longitudinal analysis of genetic and epigenetic changes in human pluripotent stem cells in the landscape of culture-induced abnormality

Human embryonic stem cells (hESCs) are naturally equipped to maintain genome integrity to minimize genetic mutations during early embryo development. However, genetic aberration risks and subsequent cellular changes in hESCs during in vitro culture pose a significant threat to stem cell therapy. While a few studies have reported specific somatic mutations and copy number variations (CNVs), the molecular mechanisms underlying the acquisition of 'culture-adapted phenotypes' by hESCs are largely unknown. Therefore, we conducted comprehensive genomic, single-cell transcriptomic, and single-cell ATAC-seq analyses of an isogenic hESC model displaying definitive 'culture-adapted phenotypes'. We found that hESCs lacking TP53, in which loss-of-function mutations were identified in human pluripotent stem cells (hPSCs), presented a surge in somatic mutations. Notably, hPSCs with a copy number gain of 20q11.21 during early passage did not present 'culture-adapted phenotypes' or BCL2L1 induction. Single-cell RNA-seq and ATAC-seq analyses revealed active transcriptional regulation at the 20q11.21 locus. Furthermore, the induction of BCL2L1 and TPX2 to trigger 'culture-adapted phenotypes' was associated with epigenetic changes facilitating TEA domain (TEAD) binding. These results suggest that 20q11.21 copy number gain and additional epigenetic changes are necessary for expressing 'culture-adapted phenotypes' by activating gene transcription at this specific locus.

Integrated Single-Cell Multiomic Profiling of Caudate Nucleus Suggests Key Mechanisms in Alcohol Use Disorder

Alcohol use disorder (AUD) induces complex transcriptional and regulatory changes across multiple brain regions including the caudate nucleus, which remains understudied. Using paired single-nucleus RNA-seq and ATAC-seq on caudate samples from 143 human postmortem brains, including 74 with AUD, we identified 17 distinct cell types. We found that a significant portion of the alcohol-induced changes in gene expression occurred through altered chromatin accessibility. Notably, we identified novel transcriptional and chromatin accessibility differences in medium spiny neurons, impacting pathways such as RNA metabolism and immune response. A small cluster of D1/D2 hybrid neurons showed distinct differences, suggesting a unique role in AUD. Microglia exhibited distinct activation states deviating from classical M1/M2 designations, and astrocytes entered a reactive state partially regulated by JUND , affecting glutamatergic synapse pathways. Oligodendrocyte dysregulation, driven in part by OLIG2 , was linked to demyelination and increased TGF-β1 signaling from microglia and astrocytes. We also observed increased microglia-astrocyte communication via the IL-1β pathway. Leveraging our multiomic data, we performed cell type-specific expression quantitative trait loci analysis, integrating that with public genome-wide association studies to identify AUD risk genes such as ADAL and PPP2R3C , providing a direct link between genetic variants, chromatin accessibility, and gene expression in AUD. These findings not only provide new insights into the genetic and cellular mechanisms in the caudate related to AUD but also demonstrate the broader utility of large-scale multiomic studies in uncovering complex gene regulation across diverse cell types, which has implications beyond the substance use field.

Differential chromatin accessibility and transcriptional dynamics define breast cancer subtypes and their lineages

Breast cancer (BC) is defined by distinct molecular subtypes with different cells of origin. The transcriptional networks that characterize the subtype-specific tumor-normal lineages are not established. In this work, we applied bulk, single-cell and single-nucleus multi-omic techniques as well as spatial transcriptomics and multiplex imaging on 61 samples from 37 patients with BC to show characteristic links in gene expression and chromatin accessibility between BC subtypes and their putative cells of origin. Regulatory network analysis of transcription factors underscored the importance of BHLHE40 in luminal BC and luminal mature cells and KLF5 in basal-like tumors and luminal progenitor cells. Furthermore, we identify key genes defining the basal-like (SOX6 and KCNQ3) and luminal A/B (FAM155A and LRP1B) lineages. Exhausted CTLA4-expressing CD8+ T cells were enriched in basal-like BC, suggesting an altered means of immune dysfunction. These findings demonstrate analysis of paired transcription and chromatin accessibility at the single-cell level is a powerful tool for investigating cancer lineage and highlight transcriptional networks that define basal and luminal BC lineages.

Degradation of IKZF1 prevents epigenetic progression of T cell exhaustion in an antigen-specific assay

In cancer, chronic antigen stimulation drives effector T cells to exhaustion, limiting the efficacy of T cell therapies. Recent studies have demonstrated that epigenetic rewiring governs the transition of T cells from effector to exhausted states and makes a subset of exhausted T cells non-responsive to PD1 checkpoint blockade. Here, we describe an antigen-specific assay for T cell exhaustion that generates T cells phenotypically and transcriptionally similar to those found in human tumors. We perform a screen of human epigenetic regulators, identifying IKZF1 as a driver of T cell exhaustion. We determine that the IKZF1 degrader iberdomide prevents exhaustion by blocking chromatin remodeling at T cell effector enhancers and preserving the binding of AP-1, NF-κB, and NFAT. Thus, our study uncovers a role for IKZF1 as a driver of T cell exhaustion through epigenetic modulation, providing a rationale for the use of iberdomide in solid tumors to prevent T cell exhaustion.

Multiomics 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 kidney failure. 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 catalog 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.

Restriction of innate Tγδ17 cell plasticity by an AP-1 regulatory axis

IL-17-producing γδ T (Tγδ17) cells are innate-like mediators of intestinal barrier immunity. While Th17 cell and ILC3 plasticity have been extensively studied, the mechanisms governing Tγδ17 cell effector flexibility remain undefined. Here, we combined type 3 fate-mapping with single cell ATAC/RNA-seq multiome profiling to define the cellular features and regulatory networks underlying Tγδ17 cell plasticity. During homeostasis, Tγδ17 cell effector identity was stable across tissues, including for intestinal T-bet+ Tγδ17 cells that restrained IFNγ production. However, S. typhimurium infection induced intestinal Vγ6+ Tγδ17 cell conversion into type 1 effectors, with loss of IL-17A production and partial RORγt downregulation. Multiome analysis revealed a trajectory along Vγ6+ Tγδ17 effector conversion, with TIM-3 marking ex-Tγδ17 cells with enhanced type 1 functionality. Lastly, we characterized and validated a critical AP-1 regulatory axis centered around JunB and Fosl2 that controls Vγ6+ Tγδ17 cell plasticity by stabilizing type 3 identity and restricting type 1 effector conversion.

Epigenetic and 3D genome reprogramming during the aging of human hippocampus

Age-related cognitive decline is associated with altered physiology of the hippocampus. While changes in gene expression have been observed in aging brain, the regulatory mechanisms underlying these changes remain underexplored. We generated single-nucleus gene expression, chromatin accessibility, DNA methylation, and 3D genome data from 40 human hippocampal tissues spanning adult lifespan. We observed a striking loss of astrocytes, OPC, and endothelial cells during aging, including astrocytes that play a role in regulating synapses. Microglia undergo a dramatic switch from a homeostatic state to a primed inflammatory state through DNA methylome and 3D genome reprogramming. Aged cells experience erosion of their 3D genome architecture. Our study identifies age-associated changes in cell types/states and gene regulatory features that provide insight into cognitive decline during human aging.

The molecular chronology of mammary epithelial cell fate switching

The adult mammary gland is maintained by lineage-restricted progenitor cells through pregnancy, lactation, involution, and menopause. Injury resolution and transplantation-associated mammary gland reconstitution are unique exceptions, wherein mammary basal cells gain the ability to reprogram to a luminal state. Here, we leverage newly developed cell-identity reporter mouse strains, and time-resolved single-cell epigenetic and transcriptomic analyses to decipher the molecular programs underlying basal-to-luminal fate switching in vivo. We demonstrate that basal cells rapidly reprogram toward plastic cycling intermediates that appear to hijack molecular programs we find in bipotent fetal mammary stem cells and puberty-associatiated cap cells. Loss of basal-cell specifiers early in dedifferentiation coincides with activation of Notch and BMP, among others. Pharmacologic blockade of each pathway disrupts basal-to-luminal transdifferentiation. Our studies provide a comprehensive map and resource for understanding the coordinated molecular changes enabling terminally differentiated epithelial cells to transition between cell lineages and highlights the stunning rapidity by which epigenetic reprogramming can occur in response to disruption of tissue structure.

Systems immunology approaches to study T cells in health and disease

T cells are dynamically regulated immune cells that are implicated in a variety of diseases ranging from infection, cancer and autoimmunity. Recent advancements in sequencing methods have provided valuable insights in the transcriptional and epigenetic regulation of T cells in various disease settings. In this review, we identify the key sequencing-based methods that have been applied to understand the transcriptomic and epigenomic regulation of T cells in diseases.

Loss of ARID3A perturbs intestinal epithelial proliferation-differentiation ratio and regeneration

Intestinal stem cells at the crypt divide and give rise to progenitor cells that proliferate and differentiate into various mature cell types in the transit-amplifying (TA) zone. Here, we showed that the transcription factor ARID3A regulates intestinal epithelial cell proliferation and differentiation at the TA progenitors. ARID3A forms an expression gradient from the villus tip to the upper crypt mediated by TGF-β and WNT. Intestinal-specific deletion of Arid3a reduces crypt proliferation, predominantly in TA cells. Bulk and single-cell transcriptomic analysis shows increased enterocyte and reduced secretory differentiation in the Arid3a cKO intestine, accompanied by enriched upper-villus gene signatures of both cell lineages. We find that the enhanced epithelial differentiation in the Arid3a-deficient intestine is caused by increased binding and transcription of HNF1 and HNF4. Finally, we show that loss of Arid3a impairs irradiation-induced regeneration with sustained cell death and reprogramming. Our findings imply that Arid3a functions to fine-tune the proliferation-differentiation dynamics at the TA progenitors, which are essential for injury-induced regeneration.

Evolution of plant cell-type-specific cis-regulatory elements

Cis-regulatory elements (CREs) are critical in regulating gene expression, and yet understanding of CRE evolution remains challenging. Here, we constructed a comprehensive single-cell atlas of chromatin accessibility in Oryza sativa, integrating data from 103,911 nuclei representing 126 discrete cell states across nine distinct organs. We used comparative genomics to compare cell-type resolved chromatin accessibility between O. sativa and 57,552 nuclei from four additional grass species (Zea mays, Sorghum bicolor, Panicum miliaceum, and Urochloa fusca). Accessible chromatin regions (ACRs) had different levels of conservation depending on the degree of cell-type specificity. We found a complex relationship between ACRs with conserved noncoding sequences, cell-type specificity, conservation, and tissue-specific switching. Additionally, we found that epidermal ACRs were less conserved compared to other cell types, potentially indicating that more rapid regulatory evolution has occurred in the L1-derived epidermal layer of these species. Finally, we identified and characterized a conserved subset of ACRs that overlapped the repressive histone modification H3K27me3, implicating them as potentially silencer-like CREs maintained by evolution. Collectively, this comparative genomics approach highlights the dynamics of plant cell-type-specific CRE evolution.

The G4 resolvase Dhx36 modulates cardiomyocyte differentiation and ventricular conduction system development

Extensive genetic studies have elucidated cardiomyocyte differentiation and associated gene networks using single-cell RNA-seq, yet the intricate transcriptional mechanisms governing cardiac conduction system (CCS) development and working cardiomyocyte differentiation remain largely unexplored. Here we show that mice deleted for Dhx36 (encoding the Dhx36 helicase) in the embryonic or neonatal heart develop overt dilated cardiomyopathy, surface ECG alterations related to cardiac impulse propagation, and (in the embryonic heart) a lack of a ventricular conduction system (VCS). Heart snRNA-seq and snATAC-seq reveal the role of Dhx36 in CCS development and in the differentiation of working cardiomyocytes. Dhx36 deficiency directly influences cardiomyocyte gene networks by disrupting the resolution of promoter G-quadruplexes in key cardiac genes, impacting cardiomyocyte differentiation and CCS morphogenesis, and ultimately leading to dilated cardiomyopathy and atrioventricular block. These findings further identify crucial genes and pathways that regulate the development and function of the VCS/Purkinje fiber (PF) network.

Cross-disorder and disease-specific pathways in dementia revealed by single-cell genomics

The development of successful therapeutics for dementias requires an understanding of their shared and distinct molecular features in the human brain. We performed single-nuclear RNA-seq and ATAC-seq in Alzheimer's disease (AD), frontotemporal dementia (FTD), and progressive supranuclear palsy (PSP), analyzing 41 participants and ∼1 million cells (RNA + ATAC) from three brain regions varying in vulnerability and pathological burden. We identify 32 shared, disease-associated cell types and 14 that are disease specific. Disease-specific cell states represent glial-immune mechanisms and selective neuronal vulnerability impacting layer 5 intratelencephalic neurons in AD, layer 2/3 intratelencephalic neurons in FTD, and layer 5/6 near-projection neurons in PSP. We identify disease-associated gene regulatory networks and cells impacted by causal genetic risk, which differ by disorder. These data illustrate the heterogeneous spectrum of glial and neuronal compositional and gene expression alterations in different dementias and identify therapeutic targets by revealing shared and disease-specific cell states.

On the identification of differentially-active transcription factors from ATAC-seq data

ATAC-seq has emerged as a rich epigenome profiling technique, and is commonly used to identify Transcription Factors (TFs) underlying given phenomena. A number of methods can be used to identify differentially-active TFs through the accessibility of their DNA-binding motif, however little is known on the best approaches for doing so. Here we benchmark several such methods using a combination of curated datasets with various forms of short-term perturbations on known TFs, as well as semi-simulations. We include both methods specifically designed for this type of data as well as some that can be repurposed for it. We also investigate variations to these methods, and identify three particularly promising approaches (a chromVAR-limma workflow with critical adjustments, monaLisa and a combination of GC smooth quantile normalization and multivariate modeling). We further investigate the specific use of nucleosome-free fragments, the combination of top methods, and the impact of technical variation. Finally, we illustrate the use of the top methods on a novel dataset to characterize the impact on DNA accessibility of TRAnscription Factor TArgeting Chimeras (TRAFTAC), which can deplete TFs-in our case NFkB-at the protein level.

The type 2 cytokine Fc-IL-4 revitalizes exhausted CD8+ T cells against cancer

Current cancer immunotherapy predominately focuses on eliciting type 1 immune responses fighting cancer; however, long-term complete remission remains uncommon1,2. A pivotal question arises as to whether type 2 immunity can be orchestrated alongside type 1-centric immunotherapy to achieve enduring response against cancer3,4. Here we show that an interleukin-4 fusion protein (Fc-IL-4), a typical type 2 cytokine, directly acts on CD8+ T cells and enriches functional terminally exhausted CD8+ T (CD8+ TTE) cells in the tumour. Consequently, Fc-IL-4 enhances antitumour efficacy of type 1 immunity-centric adoptive T cell transfer or immune checkpoint blockade therapies and induces durable remission across several syngeneic and xenograft tumour models. Mechanistically, we discovered that Fc-IL-4 signals through both signal transducer and activator of transcription 6 (STAT6) and mammalian target of rapamycin (mTOR) pathways, augmenting the glycolytic metabolism and the nicotinamide adenine dinucleotide (NAD) concentration of CD8+ TTE cells in a lactate dehydrogenase A-dependent manner. The metabolic modulation mediated by Fc-IL-4 is indispensable for reinvigorating intratumoural CD8+ TTE cells. These findings underscore Fc-IL-4 as a potent type 2 cytokine-based immunotherapy that synergizes effectively with type 1 immunity to elicit long-lasting responses against cancer. Our study not only sheds light on the synergy between these two types of immune responses, but also unveils an innovative strategy for advancing next-generation cancer immunotherapy by integrating type 2 immune factors.

Single-cell CAR T atlas reveals type 2 function in 8-year leukaemia remission

Despite a high response rate in chimeric antigen receptor (CAR) T cell therapy for acute lymphocytic leukaemia (ALL)1-3, approximately 50% of patients relapse within the first year4-6, representing an urgent question to address in the next stage of cellular immunotherapy. Here, to investigate the molecular determinants of ultralong CAR T cell persistence, we obtained a single-cell multi-omics atlas from 695,819 pre-infusion CAR T cells at the basal level or after CAR-specific stimulation from 82 paediatric patients with ALL enrolled in the first two CAR T ALL clinical trials and 6 healthy donors. We identified that elevated type 2 functionality in CAR T infusion products is significantly associated with patients maintaining a median B cell aplasia duration of 8.4 years. Analysis of ligand-receptor interactions revealed that type 2 cells regulate a dysfunctional subset to maintain whole-population homeostasis, and the addition of IL-4 during antigen-specific activation alleviates CAR T cell dysfunction while enhancing fitness at both transcriptomic and epigenomic levels. Serial proteomic profiling of sera after treatment revealed a higher level of circulating type 2 cytokines in 5-year or 8-year relapse-free responders. In a leukaemic mouse model, type 2high CAR T cell products demonstrated superior expansion and antitumour activity, particularly after leukaemia rechallenge. Restoring antitumour efficacy in type 2low CAR T cells was attainable by enhancing their type 2 functionality, either through incorporating IL-4 into the manufacturing process or by priming manufactured CAR T products with IL-4 before infusion. Our findings provide insights into the mediators of durable CAR T therapy response and suggest potential therapeutic strategies to sustain long-term remission by boosting type 2 functionality in CAR T cells.

DRCTdb: disease-related cell type analysis to decode cell type effect and underlying regulatory mechanisms

Understanding the molecular mechanisms underlying genetic diseases is challenging due to environmental and genetic factors. Genome-wide association studies (GWAS) have identified numerous genetic loci, but their functional implications are largely unknown. Single-cell multiomics sequencing has emerged as a powerful tool to study disease-specific cell types and their relationship with genetic variants. However, comprehensive databases for exploring these mechanisms across different tissues are lacking. We present the Disease-Related Cell Type database (DRCTdb), integrating GWAS and single-cell multiomics data to identify disease-related cell types and elucidate their regulatory mechanisms. DRCTdb contains well-processed data from 16 studies, covering 4 million cells within 28 tissues. Users can browse relationships and regulatory mechanisms between SNPs of 42 genetic diseases and cell types based on GWAS and single-cell data. DRCTdb also offers data downloads and is available at https://singlecellatlas.top/DRCTDB .

The chromatin landscape of high-grade serous ovarian cancer metastasis identifies regulatory drivers in post-chemotherapy residual tumour cells

Disease recurrence following chemotherapy is a major clinical challenge in ovarian cancer (OC), but little is known regarding how the tumour epigenome regulates transcriptional programs underpinning chemoresistance. We determine the single cell chromatin accessibility landscape of omental OC metastasis from treatment-naïve and neoadjuvant chemotherapy-treated patients and define the chromatin accessibility profiles of epithelial, fibroblast, myeloid and lymphoid cells. Epithelial tumour cells display open chromatin regions enriched with motifs for the oncogenic transcription factors MEIS and PBX. Post chemotherapy microenvironments show profound tumour heterogeneity and selection for cells with accessible chromatin enriched for TP53, TP63, TWIST1 and resistance-pathway-activating transcription factor binding motifs. An OC chemoresistant tumour subpopulation known to be present prior to treatment, and characterised by stress-associated gene expression, is enriched post chemotherapy. Nuclear receptors RORa, NR2F6 and HNF4G are uncovered as candidate transcriptional drivers of these cells whilst closure of binding sites for E2F2 and E2F4 indicate post-treated tumour having low proliferative capacity. Delineation of the gene regulatory landscape of ovarian cancer cells surviving chemotherapy treatment therefore reveals potential core transcriptional regulators of chemoresistance, suggesting novel therapeutic targets for improving clinical outcome.

Single-nucleus multiomics reveals the gene-regulatory networks underlying sex determination of murine primordial germ cells

Accurate specification of female and male germ cells during embryonic development is critical for sexual reproduction. Primordial germ cells (PGCs) are the bipotential precursors of mature gametes that commit to an oogenic or spermatogenic fate in response to sex-determining cues from the fetal gonad. The critical processes required for PGCs to integrate and respond to signals from the somatic environment in gonads are not understood. In this study, we developed the first single-nucleus multiomics map of chromatin accessibility and gene expression during murine PGC development in both XX and XY embryos. Profiling of cell-type specific transcriptomes and regions of open chromatin from the same cell captured the molecular signatures and gene networks underlying PGC sex determination. Joint RNA and ATAC data for single PGCs resolved previously unreported PGC subpopulations and cataloged a multimodal reference atlas of differentiating PGC clusters. We discovered that regulatory element accessibility precedes gene expression during PGC development, suggesting that changes in chromatin accessibility may prime PGC lineage commitment prior to differentiation. Similarly, we found that sexual dimorphism in chromatin accessibility and gene expression increased temporally in PGCs. Combining single-nucleus sequencing data, we computationally mapped the cohort of transcription factors that regulate the expression of sexually dimorphic genes in PGCs. For example, the gene regulatory networks of XX PGCs are enriched for the transcription factors, TFAP2c, TCFL5, GATA2, MGA, NR6A1, TBX4, and ZFX. Sex-specific enrichment of the forkhead-box and POU6 families of transcription factors was also observed in XY PGCs. Finally, we determined the temporal expression patterns of WNT, BMP, and RA signaling during PGC sex determination, and our discovery analyses identified potentially new cell communication pathways between supporting cells and PGCs. Our results illustrate the diversity of factors involved in programming PGCs towards a sex-specific fate.

Acute and persistent responses after H5N1 vaccination in humans

To gain insight into how an adjuvant impacts vaccination responses, we use systems immunology to study human H5N1 influenza vaccination with or without the adjuvant AS03, longitudinally assessing 14 time points including multiple time points within the first day after prime and boost. We develop an unsupervised computational framework to discover high-dimensional response patterns, which uncover adjuvant- and immunogenicity-associated early response dynamics, including some that differ post prime versus boost. With or without adjuvant, some vaccine-induced transcriptional patterns persist to at least 100 days after initial vaccination. Single-cell profiling of surface proteins, transcriptomes, and chromatin accessibility implicates transcription factors in the erythroblast-transformation-specific (ETS) family as shaping these long-lasting signatures, primarily in classical monocytes but also in CD8+ naive-like T cells. These cell-type-specific signatures are elevated at baseline in high-antibody responders in an independent vaccination cohort, suggesting that antigen-agnostic baseline immune states can be modulated by vaccine antigens alone to enhance future responses.

scooby: Modeling multi-modal genomic profiles from DNA sequence at single-cell resolution

Understanding how regulatory DNA elements shape gene expression across individual cells is a fundamental challenge in genomics. Joint RNA-seq and epigenomic profiling provides opportunities to build unifying models of gene regulation capturing sequence determinants across steps of gene expression. However, current models, developed primarily for bulk omics data, fail to capture the cellular heterogeneity and dynamic processes revealed by single-cell multi-modal technologies. Here, we introduce scooby, the first model to predict scRNA-seq coverage and scATAC-seq insertion profiles along the genome from sequence at single-cell resolution. For this, we leverage the pre-trained multi-omics profile predictor Borzoi as a foundation model, equip it with a cell-specific decoder, and fine-tune its sequence embeddings. Specifically, we condition the decoder on the cell position in a precomputed single-cell embedding resulting in strong generalization capability. Applied to a hematopoiesis dataset, scooby recapitulates cell-specific expression levels of held-out genes and cells, and identifies regulators and their putative target genes through in silico motif deletion. Moreover, accurate variant effect prediction with scooby allows for breaking down bulk eQTL effects into single-cell effects and delineating their impact on chromatin accessibility and gene expression. We anticipate scooby to aid unraveling the complexities of gene regulation at the resolution of individual cells.

Protocol to identify defined reprogramming factor expression using a factor-indexing single-nuclei multiome sequencing approach

Ectopic expression of lineage-specific transcription factors (TFs) of another cell type can induce cell fate reprogramming. However, the heterogeneity of reprogramming cells has been a challenge for data interpretation and model evaluation. Here, we present a protocol to characterize cells expressing defined factors during direct cell reprogramming using a factor-indexing approach based on single-nuclei multiome sequencing (FI-snMultiome-seq). We describe the steps for barcoding TFs, converting human fibroblasts to pancreatic ductal-like cells using defined TFs, and preparing library for FI-snMultiome-seq analysis. For complete details on the use and execution of this protocol, please refer to Fei et al.1.

Transformation of brain myeloid cell populations by SIV in rhesus macaques revealed by multiomics

The primary immune constituents in the brain, microglia and macrophages, are the target for HIV in people and simian immunodeficiency virus (SIV) in nonhuman primates. This infection can lead to neurological dysfunction, known as HIV-associated neurocognitive disorder (HAND). Given the gaps in our knowledge on how these cells respond in vivo to CNS infection, we performed single-cell multiomic sequencing, including gene expression and ATAC-seq, on myeloid cells from the brains of rhesus macaques with SIV-induced encephalitis (SIVE) as well as uninfected controls. We found that the myeloid cell populations were significantly changed by SIVE. In SIVE microglia-like cells express high levels of chemoattractants capable of recruiting highly activated CAM-like cells to the site of infection/inflammation. A unique population of microglia-like cells was found in which the chromatin accessibility of genes diverged from their RNA expression. Additionally, we observed a dramatic shift of upstream gene regulators and their targets in brain myeloid cells during SIVE. In summary, this study further uncovers the transcriptome, gene regulatory events and potential roles of different brain myeloid phenotypes in SIVE.

Organoid culture promotes dedifferentiation of mouse myoblasts into stem cells capable of complete muscle regeneration

Experimental cell therapies for skeletal muscle conditions have shown little success, primarily because they use committed myogenic progenitors rather than true muscle stem cells, known as satellite cells. Here we present a method to generate in vitro-derived satellite cells (idSCs) from skeletal muscle tissue. When transplanted in small numbers into mouse muscle, mouse idSCs fuse into myofibers, repopulate the satellite cell niche, self-renew, support multiple rounds of muscle regeneration and improve force production on par with freshly isolated satellite cells in damaged skeletal muscle. We compared the epigenomic and transcriptional signatures between idSCs, myoblasts and satellite cells and used these signatures to identify core signaling pathways and genes that confer idSC functionality. Finally, from human muscle biopsies, we successfully generated satellite cell-like cells in vitro. After further development, idSCs may provide a scalable source of cells for the treatment of genetic muscle disorders, trauma-induced muscle damage and age-related muscle weakness.

Epigenetic regulation of cell state by H2AFY governs immunogenicity in high-risk neuroblastoma

Childhood neuroblastoma with MYCN amplification is classified as high risk and often relapses after intensive treatments. Immune checkpoint blockade therapy against the PD-1/L1 axis shows limited efficacy in patients with neuroblastoma, and the cancer intrinsic immune regulatory network is poorly understood. Here, we leverage genome-wide CRISPR/Cas9 screens and identify H2AFY as a resistance gene to the clinically approved PD-1 blocking antibody nivolumab. Analysis of single-cell RNA-Seq datasets reveals that H2AFY mRNA is enriched in adrenergic cancer cells and is associated with worse patient survival. Genetic deletion of H2afy in MYCN-driven neuroblastoma cells reverts in vivo resistance to PD-1 blockade by eliciting activation of the adaptive and innate immunity. Mapping of the epigenetic and translational landscape demonstrates that H2afy deletion promotes cell transition to a mesenchymal-like state. With a multiomics approach, we uncovered H2AFY-associated genes that are functionally relevant and prognostic in patients. Altogether, our study elucidates the role of H2AFY as an epigenetic gatekeeper for cell states and immunogenicity in high-risk neuroblastoma.

Joint trajectory inference for single-cell genomics using deep learning with a mixture prior

Trajectory inference methods are essential for analyzing the developmental paths of cells in single-cell sequencing datasets. It provides insights into cellular differentiation, transitions, and lineage hierarchies, helping unravel the dynamic processes underlying development and disease progression. However, many existing tools lack a coherent statistical model and reliable uncertainty quantification, limiting their utility and robustness. In this paper, we introduce VITAE (Variational Inference for Trajectory by AutoEncoder), a statistical approach that integrates a latent hierarchical mixture model with variational autoencoders to infer trajectories. The statistical hierarchical model enhances the interpretability of our framework, while the posterior approximations generated by our variational autoencoder ensure computational efficiency and provide uncertainty quantification of cell projections along trajectories. Specifically, VITAE enables simultaneous trajectory inference and data integration, improving the accuracy of learning a joint trajectory structure in the presence of biological and technical heterogeneity across datasets. We show that VITAE outperforms other state-of-the-art trajectory inference methods on both real and synthetic data under various trajectory topologies. Furthermore, we apply VITAE to jointly analyze three distinct single-cell RNA sequencing datasets of the mouse neocortex, unveiling comprehensive developmental lineages of projection neurons. VITAE effectively reduces batch effects within and across datasets and uncovers finer structures that might be overlooked in individual datasets. Additionally, we showcase VITAE's efficacy in integrative analyses of multiomic datasets with continuous cell population structures.

Mechano-osmotic signals control chromatin state and fate transitions in pluripotent stem cells

Acquisition of specific cell shapes and morphologies is a central component of cell fate transitions. Although signaling circuits and gene regulatory networks that regulate pluripotent stem cell differentiation have been intensely studied, how these networks are integrated in space and time with morphological transitions and mechanical deformations to control state transitions remains a fundamental open question. Here, we focus on two distinct models of pluripotency, primed pluripotent stem cells and pre-implantation inner cell mass cells of human embryos to discover that cell fate transitions associate with rapid changes in nuclear shape and volume which collectively alter the nuclear mechanophenotype. Mechanistic studies in human induced pluripotent stem cells further reveal that these phenotypical changes and the associated active fluctuations of the nuclear envelope arise from growth factor signaling-controlled changes in chromatin mechanics and cytoskeletal confinement. These collective mechano-osmotic changes trigger global transcriptional repression and a condensation-prone environment that primes chromatin for a cell fate transition by attenuating repression of differentiation genes. However, while this mechano-osmotic chromatin priming has the potential to accelerate fate transitions and differentiation, sustained biochemical signals are required for robust induction of specific lineages. Our findings uncover a critical mechanochemical feedback mechanism that integrates nuclear mechanics, shape and volume with biochemical signaling and chromatin state to control cell fate transition dynamics.

Flexible use of conserved motif vocabularies constrains genome access in cell type evolution

Cell types evolve into a hierarchy with related types grouped into families. How cell type diversification is constrained by the stable separation between families over vast evolutionary times remains unknown. Here, integrating single-nucleus multiomic sequencing and deep learning, we show that hundreds of sequence features (motifs) divide into distinct sets associated with accessible genomes of specific cell type families. This division is conserved across highly divergent, early-branching animals including flatworms and cnidarians. While specific interactions between motifs delineate cell type relationships within families, surprisingly, these interactions are not conserved between species. Consistently, while deep learning models trained on one species can predict accessibility of other species' sequences, their predictions frequently rely on distinct, but synonymous, motif combinations. We propose that long-term stability of cell type families is maintained through genome access specified by conserved motif sets, or 'vocabularies', whereas cell types diversify through flexible use of motifs within each set.

Mir221/222 drive synovial hyperplasia and arthritis by targeting cell cycle inhibitors and chromatin remodeling components

miRNAs constitute fine-tuners of gene expression and are implicated in a variety of diseases spanning from inflammation to cancer. miRNA expression is deregulated in rheumatoid arthritis (RA); however, their specific role in key arthritogenic cells such as the synovial fibroblast (SF) remains elusive. Previous studies have shown that Mir221/222 expression is upregulated in RA SFs. Here, we demonstrate that TNF and IL-1β but not IFN-γ activated Mir221/222 gene expression in murine SFs. SF-specific overexpression of Mir221/222 in huTNFtg mice led to further expansion of SFs and disease exacerbation, while its total ablation led to reduced SF expansion and attenuated disease. Mir221/222 overexpression altered the SF transcriptional profile igniting pathways involved in cell cycle and ECM (extracellular matrix) regulation. Validation of targets of Mir221/222 revealed cell cycle inhibitors Cdkn1b and Cdkn1c, as well as the epigenetic regulator Smarca1. Single-cell ATAC-seq data analysis revealed increased Mir221/222 gene activity in pathogenic SF subclusters and transcriptional regulation by Rela, Relb, Junb, Bach1, and Nfe2l2. Our results establish an SF-specific pathogenic role of Mir221/222 in arthritis and suggest that its therapeutic targeting in specific subpopulations could lead to novel fibroblast-targeted therapies.

An integrated transcription factor framework for Treg identity and diversity

Vertebrate cell identity depends on the combined activity of scores of transcription factors (TF). While TFs have often been studied in isolation, a systematic perspective on their integration has been missing. Focusing on FoxP3+ regulatory T cells (Tregs), key guardians of immune tolerance, we combined single-cell chromatin accessibility, machine learning, and high-density genetic variation, to resolve a validated framework of diverse Treg chromatin programs, each shaped by multi-TF inputs. This framework identified previously unrecognized Treg controllers (Smarcc1) and illuminated the mechanism of action of FoxP3, which amplified a pre-existing Treg identity, diversely activating or repressing distinct programs, dependent on different regulatory partners. Treg subpopulations in the colon relied variably on FoxP3, Helios+ Tregs being completely dependent, but RORγ+ Tregs largely independent. These differences were rooted in intrinsic biases decoded by the integrated framework. Moving beyond master regulators, this work unravels how overlapping TF activities coalesce into Treg identity and diversity.

Single-cell EpiChem jointly measures drug-chromatin binding and multimodal epigenome

Studies of molecular and cellular functions of small-molecule inhibitors in cancer treatment, eliciting effects by targeting genome and epigenome associated proteins, requires measurement of drug-target engagement in single-cell resolution. Here we present EpiChem for in situ single-cell joint mapping of small molecules and multimodal epigenomic landscape. We demonstrate single-cell co-assays of three small molecules together with histone modifications, chromatin accessibility or target proteins in human colorectal cancer (CRC) organoids. Integrated multimodal analysis reveals diverse drug interactions in the context of chromatin states within heterogeneous CRC organoids. We further reveal drug genomic binding dynamics and adaptive epigenome across cell types after small-molecule drug treatment in CRC organoids. This method provides a unique tool to exploit the mechanisms of cell type-specific drug actions.

Deciphering the Transcriptional Regulatory Network Governing Starch and Storage Protein Biosynthesis in Wheat for Breeding Improvement

Starch and seed storage protein (SSP) composition profoundly impact wheat grain yield and quality. To unveil regulatory mechanisms governing their biosynthesis, transcriptome, and epigenome profiling is conducted across key endosperm developmental stages, revealing that chromatin accessibility, H3K27ac, and H3K27me3 collectively regulate SSP and starch genes with varying impact. Population transcriptome and phenotype analyses highlight accessible promoter regions' crucial role as a genetic variation resource, influencing grain yield and quality in a core collection of wheat accessions. Integration of time-serial RNA-seq and ATAC-seq enables the construction of a hierarchical transcriptional regulatory network governing starch and SSP biosynthesis, identifying 42 high-confidence novel candidates. These candidates exhibit overlap with genetic regions associated with grain size and quality traits, and their functional significance is validated through expression-phenotype association analysis among wheat accessions and loss-of-function mutants. Functional analysis of wheat abscisic acid insensitive 3-A1 (TaABI3-A1) with genome editing knock-out lines demonstrates its role in promoting SSP accumulation while repressing starch biosynthesis through transcriptional regulation. Excellent TaABI3-A1Hap1 with enhanced grain weight is selected during the breeding process in China, linked to altered expression levels. This study unveils key regulators, advancing understanding of SSP and starch biosynthesis regulation and contributing to breeding enhancement.

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

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

Dual function of PHF16 in reinstating homeostasis of murine intestinal epithelium after crypt regeneration

Intestinal stem cells (ISCs) are highly vulnerable to damage, being in a constant state of proliferation. Reserve stem cells repair the intestinal epithelium following damage-induced ablation of ISCs. Here, we report that the epigenetic regulator plant homology domain (PHD) finger protein 16 (PHF16) restores homeostasis of the intestinal epithelium after initial damage-induced repair. In Phf16-/Y mice, revival stem cells (revSCs) showed defects in exiting the regenerative state, and intestinal crypt regeneration failed even though revSCs were still induced in response to tissue damage, as observed by single-cell RNA sequencing (scRNA-seq). Analysis of Phf16-/Y intestinal organoids by RNA sequencing (RNA-seq) and ATAC sequencing identified that PHF16 restores homeostasis of the intestinal epithelium by inducing retinoic acid receptor (RAR)/retinoic X receptor (RXR) target genes through HBO1-mediated histone H3K14 acetylation, while at the same time counteracting YAP/TAZ activity by ubiquitination of CDC73. Together, our findings demonstrate the importance of timely suppression of regenerative activity by PHF16 for the restoration of gut homeostasis after acute tissue injury.

Dissecting the regulatory logic of specification and differentiation during vertebrate embryogenesis

The interplay between transcription factors and chromatin accessibility regulates cell type diversification during vertebrate embryogenesis. To systematically decipher the gene regulatory logic guiding this process, we generated a single-cell multi-omics atlas of RNA expression and chromatin accessibility during early zebrafish embryogenesis. We developed a deep learning model to predict chromatin accessibility based on DNA sequence and found that a small number of transcription factors underlie cell-type-specific chromatin landscapes. While Nanog is well-established in promoting pluripotency, we discovered a new function in priming the enhancer accessibility of mesendodermal genes. In addition to the classical stepwise mode of differentiation, we describe instant differentiation, where pluripotent cells skip intermediate fate transitions and terminally differentiate. Reconstruction of gene regulatory interactions reveals that this process is driven by a shallow network in which maternally deposited regulators activate a small set of transcription factors that co-regulate hundreds of differentiation genes. Notably, misexpression of these transcription factors in pluripotent cells is sufficient to ectopically activate their targets. This study provides a rich resource for analyzing embryonic gene regulation and reveals the regulatory logic of instant differentiation.

Single-cell multiomic analysis identifies macrophage subpopulations in promoting cardiac repair

Cardiac mononuclear phagocytic cells (Cardiac MPCs) participate in maintaining homeostasis and orchestrating cardiac responses upon injury. However, the function of specific MPC subtypes and the related cell fate commitment mechanisms remain elusive in regenerative and nonregenerative hearts due to their cellular heterogeneities. Using spatiotemporal single-cell epigenomic analysis of cardiac MPCs in regenerative (P1) and nonregenerative (P10) mouse hearts after injury, we found that P1 hearts accumulate reparative Arg1+ macrophages, while proinflammatory S100a9+Ly6c+ monocytes are uniquely abundant during nonregenerative remodeling. Moreover, blocking chemokine CXCR2 to inhibit the specification of the S100a9+Ly6c+-biased inflammatory fate in P10 hearts resulted in elevated wound repair responses and marked improvements in cardiac function after injury. Single-cell RNA-Seq further confirmed an increased Arg1+ macrophage subpopulation after CXCR2 blockade, which was accomplished by increased expression of wound repair-related genes and reduced expression of proinflammatory genes. Collectively, our findings provide instructive insights into the molecular mechanisms underlying the function and fate specification of heterogeneous MPCs during cardiac repair and identify potential therapeutic targets for myocardial infarction.

IL-4-induced SOX9 confers lineage plasticity to aged adult lung stem cells

Wound healing in response to acute injury is mediated by the coordinated and transient activation of parenchymal, stromal, and immune cells that resolves to homeostasis. Environmental, genetic, and epigenetic factors associated with inflammation and aging can lead to persistent activation of the microenvironment and fibrosis. Here, we identify opposing roles of interleukin-4 (IL-4) cytokine signaling in interstitial macrophages and type II alveolar epithelial cells (ATIIs). We show that IL4Ra signaling in macrophages promotes regeneration of the alveolar epithelium after bleomycin-induced lung injury. Using organoids and mouse models, we show that IL-4 directly acts on a subset of ATIIs to induce the expression of the transcription factor SOX9 and reprograms them toward a progenitor-like state with both airway and alveolar lineage potential. In the contexts of aging and bleomycin-induced lung injury, this leads to aberrant epithelial cell differentiation and bronchiolization, consistent with cellular and histological changes observed in interstitial lung disease.

Baf155 controls hematopoietic differentiation and regeneration through chromatin priming

Chromatin priming promotes cell-type-specific gene expression, lineage differentiation, and development. The mechanism of chromatin priming has not been fully understood. Here, we report that mouse hematopoietic stem and progenitor cells (HSPCs) lacking the Baf155 subunit of the BAF (BRG1/BRM-associated factor) chromatin remodeling complex produce a significantly reduced number of mature blood cells, leading to a failure of hematopoietic regeneration upon transplantation and 5-fluorouracil (5-FU) injury. Baf155-deficient HSPCs generate particularly fewer neutrophils, B cells, and CD8+ T cells at homeostasis, supporting a more immune-suppressive tumor microenvironment and enhanced tumor growth. Single-nucleus multiomics analysis reveals that Baf155-deficient HSPCs fail to establish accessible chromatin in selected regions that are enriched for putative enhancers and binding motifs of hematopoietic lineage transcription factors. Our study provides a fundamental mechanistic understanding of the role of Baf155 in hematopoietic lineage chromatin priming and the functional consequences of Baf155 deficiency in regeneration and tumor immunity.

Epigenetic insights into GABAergic development in Dravet Syndrome iPSC and therapeutic implications

Dravet syndrome (DS) is a devastating early-onset refractory epilepsy syndrome caused by variants in the SCN1A gene. A disturbed GABAergic interneuron function is implicated in the progression to DS but the underlying developmental and pathophysiological mechanisms remain elusive, in particularly at the chromatin level. Induced pluripotent stem cells (iPSCs) derived from DS cases and healthy donors were used to model disease-associated epigenetic abnormalities of GABAergic development. Chromatin accessibility was assessed at multiple time points (Day 0, Day 19, Day 35, and Day 65) of GABAergic differentiation. Additionally, the effects of the commonly used anti-seizure drug valproic acid (VPA) on chromatin accessibility were elucidated in GABAergic cells. The distinct dynamics in the chromatin profile of DS iPSC predicted accelerated early GABAergic development, evident at D19, and diverged further from the pattern in control iPSC with continued differentiation, indicating a disrupted GABAergic maturation. Exposure to VPA at D65 reshaped the chromatin landscape at a variable extent in different iPSC-lines and rescued the observed dysfunctional development of some DS iPSC-GABA. The comprehensive investigation on the chromatin landscape of GABAergic differentiation in DS-patient iPSC offers valuable insights into the epigenetic dysregulations associated with interneuronal dysfunction in DS. Moreover, the detailed analysis of the chromatin changes induced by VPA in iPSC-GABA holds the potential to improve the development of personalized and targeted anti-epileptic therapies.

Innate immune memory after brain injury drives inflammatory cardiac dysfunction

The medical burden of stroke extends beyond the brain injury itself and is largely determined by chronic comorbidities that develop secondarily. We hypothesized that these comorbidities might share a common immunological cause, yet chronic effects post-stroke on systemic immunity are underexplored. Here, we identify myeloid innate immune memory as a cause of remote organ dysfunction after stroke. Single-cell sequencing revealed persistent pro-inflammatory changes in monocytes/macrophages in multiple organs up to 3 months after brain injury, notably in the heart, leading to cardiac fibrosis and dysfunction in both mice and stroke patients. IL-1β was identified as a key driver of epigenetic changes in innate immune memory. These changes could be transplanted to naive mice, inducing cardiac dysfunction. By neutralizing post-stroke IL-1β or blocking pro-inflammatory monocyte trafficking with a CCR2/5 inhibitor, we prevented post-stroke cardiac dysfunction. Such immune-targeted therapies could potentially prevent various IL-1β-mediated comorbidities, offering a framework for secondary prevention immunotherapy.

Transcriptome data are insufficient to control false discoveries in regulatory network inference

Inference of causal transcriptional regulatory networks (TRNs) from transcriptomic data suffers notoriously from false positives. Approaches to control the false discovery rate (FDR), for example, via permutation, bootstrapping, or multivariate Gaussian distributions, suffer from several complications: difficulty in distinguishing direct from indirect regulation, nonlinear effects, and causal structure inference requiring "causal sufficiency," meaning experiments that are free of any unmeasured, confounding variables. Here, we use a recently developed statistical framework, model-X knockoffs, to control the FDR while accounting for indirect effects, nonlinear dose-response, and user-provided covariates. We adjust the procedure to estimate the FDR correctly even when measured against incomplete gold standards. However, benchmarking against chromatin immunoprecipitation (ChIP) and other gold standards reveals higher observed than reported FDR. This indicates that unmeasured confounding is a major driver of FDR in TRN inference. A record of this paper's transparent peer review process is included in the supplemental information.

On the identification of differentially-active transcription factors from ATAC-seq data

ATAC-seq has emerged as a rich epigenome profiling technique, and is commonly used to identify Transcription Factors (TFs) underlying given phenomena. A number of methods can be used to identify differentially-active TFs through the accessibility of their DNA-binding motif, however little is known on the best approaches for doing so. Here we benchmark several such methods using a combination of curated datasets with various forms of short-term perturbations on known TFs, as well as semi-simulations. We include both methods specifically designed for this type of data as well as some that can be repurposed for it. We also investigate variations to these methods, and identify three particularly promising approaches (a chromVAR-limma workflow with critical adjustments, monaLisa and a combination of GC smooth quantile normalization and multivariate modeling). We further investigate the specific use of nucleosome-free fragments, the combination of top methods, and the impact of technical variation. Finally, we illustrate the use of the top methods on a novel dataset to characterize the impact on DNA accessibility of TRAnscription Factor TArgeting Chimeras (TRAFTAC), which can deplete TFs - in our case NFkB - at the protein level.

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.