Single-cell chromatin state analysis with Signac

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

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

Transient chromatin decompaction at the start of D. melanogaster male embryonic germline development

Embryonic germ cells develop rapidly to establish the foundation for future developmental trajectories, and in this process, they make critical lineage choices including the configuration of their unique identity and a decision on sex. Here, we use single-cell genomics patterns for the entire embryonic germline in Drosophila melanogaster along with the somatic gonadal precursors after embryonic gonad coalescence to investigate molecular mechanisms involved in the setting up and regulation of the germline program. Profiling of the early germline chromatin landscape revealed sex- and stage-specific features. In the male germline immediately after zygotic activation, the chromatin structure underwent a brief remodeling phase during which nucleosome density was lower and deconcentrated from promoter regions. These findings echoed enrichment analysis results of our genomics data in which top candidates were factors with the ability to mediate large-scale chromatin reorganization. Together, they point to the importance of chromatin regulation in the early germline and raise the possibility of a conserved epigenetic reprogramming-like process required for proper initiation of germline development.

Identification of key genes and immune infiltration in peripheral blood biomarker analysis of delayed cerebral ischemia: Valproic acid as a potential therapeutic drug

BACKGROUND

Delayed cerebral ischemia (DCI) is a common and serious complication of subarachnoid hemorrhage (SAH). Its pathogenesis is not fully understood. Here, we developed a predictive model based on peripheral blood biomarkers and validated the model using several bioinformatic multi-analysis methods.

METHODS

Six datasets were obtained from the GEO database. Characteristic genes were screened using weighted correlation network analysis (WGCNA) and differentially expressed genes. Three machine learning algorithms, elastic networks-LASSO, support vector machines (SVM-RFE) and random forests (RF), were also used to construct diagnostic prediction models for key genes. To further evaluate the performance and predictive value of the diagnostic models, nomogram model were constructed, and the clinical value of the models was assessed using Decision Curve Analysis (DCA), Area Under the Check Curve (AUC), Clinical Impact Curve (CIC), and validated in the mouse single-cell RNA-seq dataset. Mendelian randomization(MR) analysis explored the causal relationship between SAH and stroke, and the intermediate influencing factors. We validated this by retrospectively analyzing the qPCR levels of the most relevant genes in SAH and SAH-DCI patients. This experiment demonstrated a statistically significant difference between SAH and SAH-DCI and normal group controls. Finally, potential small molecule compounds interacting with the selected features were screened from the Comparative Toxicogenomics Database (CTD).

RESULTS

The fGSEA results showed that activation of Toll-like receptor signaling and leukocyte transendothelial cell migration pathways were positively correlated with the DCI phenotype, whereas cytokine signaling pathways and natural killer cell-mediated cytotoxicity were negatively correlated. Consensus feature selection of DEG genes using WGCNA and three machine learning algorithms resulted in the identification of six genes (SPOCK2, TRRAP, CIB1, BCL11B, PDZD8 and LAT), which were used to predict DCI diagnosis with high accuracy. Three external datasets and the mouse single-cell dataset showed high accuracy of the diagnostic model, in addition to high performance and predictive value of the diagnostic model in DCA and CIC. MR analysis looked at stroke after SAH independent of SAH, but associated with multiple intermediate factors including Hypertensive diseases, Total triglycerides levels in medium HDL and Platelet count. qPCR confirmed that significant differences in DCI signature genes were observed between the SAH and SAH-DCI groups. Finally, valproic acid became a potential therapeutic agent for DCI based on the results of target prediction and molecular docking of the characterized genes.

CONCLUSION

This diagnostic model can identify SAH patients at high risk for DCI and may provide potential mechanisms and therapeutic targets for DCI. Valproic acid may be an important future drug for the treatment of DCI.

Chromatin plasticity predetermines neuronal eligibility for memory trace formation

Memories are encoded by sparse populations of neurons but how such sparsity arises remains largely unknown. We found that a neuron's eligibility to be recruited into the memory trace depends on its epigenetic state prior to encoding. Principal neurons in the mouse lateral amygdala display intrinsic chromatin plasticity, which when experimentally elevated favors neuronal allocation into the encoding ensemble. Such chromatin plasticity occurred at genomic regions underlying synaptic plasticity and was accompanied by increased neuronal excitability in single neurons in real time. Lastly, optogenetic silencing of the epigenetically altered neurons prevented memory expression, revealing a cell-autonomous relationship between chromatin plasticity and memory trace formation. These results identify the epigenetic state of a neuron as a key factor enabling information encoding.

In vivo interaction screening reveals liver-derived constraints to metastasis

It is estimated that only 0.02% of disseminated tumour cells are able to seed overt metastases1. While this suggests the presence of environmental constraints to metastatic seeding, the landscape of host factors controlling this process remains largely unclear. Here, combining transposon technology2 and fluorescence niche labelling3, we developed an in vivo CRISPR activation screen to systematically investigate the interactions between hepatocytes and metastatic cells. We identify plexin B2 as a critical host-derived regulator of liver colonization in colorectal and pancreatic cancer and melanoma syngeneic mouse models. We dissect a mechanism through which plexin B2 interacts with class IV semaphorins on tumour cells, leading to KLF4 upregulation and thereby promoting the acquisition of epithelial traits. Our results highlight the essential role of signals from the liver parenchyma for the seeding of disseminated tumour cells before the establishment of a growth-promoting niche. Our findings further suggest that epithelialization is required for the adaptation of CRC metastases to their new tissue environment. Blocking the plexin-B2-semaphorin axis abolishes metastatic colonization of the liver and therefore represents a therapeutic strategy for the prevention of hepatic metastases. Finally, our screening approach, which evaluates host-derived extrinsic signals rather than tumour-intrinsic factors for their ability to promote metastatic seeding, is broadly applicable and lays a framework for the screening of environmental constraints to metastasis in other organs and cancer types.

Kmt2c restricts G-CSF-driven HSC mobilization and granulocyte production in a methyltransferase-independent manner

Granulocyte colony-stimulating factor (G-CSF) is widely used to enhance myeloid recovery after chemotherapy and to mobilize hematopoietic stem cells (HSCs) for transplantation. Unfortunately, through the course of chemotherapy, cancer patients can acquire leukemogenic mutations that cause therapy-related myelodysplastic syndrome (MDS) or acute myeloid leukemia (AML). This raises the question of whether therapeutic G-CSF might potentiate therapy-related MDS/AML by disproportionately stimulating mutant HSCs and other myeloid progenitors. A common mutation in therapy-related MDS/AML involves chromosome 7 deletions that inactivate many tumor suppressor genes, including KMT2C. Here, we show that Kmt2c deletions hypersensitize murine HSCs and myeloid progenitors to G-CSF, as evidenced by increased HSC mobilization and enhanced granulocyte production from granulocyte-monocyte progenitors (GMPs). Furthermore, Kmt2c attenuates the G-CSF response independently from its SET methyltransferase function. Altogether, the data raise concerns that monosomy 7 can hypersensitize progenitors to G-CSF, such that clinical use of G-CSF may amplify the risk of therapy-related MDS/AML.

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Elucidating allergic reaction mechanisms in response to SARS-CoV-2 mRNA vaccination in adults

BACKGROUND

During the COVID-19 pandemic, novel nanoparticle-based mRNA vaccines were developed. A small number of individuals developed allergic reactions to these vaccines although the mechanisms remain undefined.

METHODS

To understand COVID-19 vaccine-mediated allergic reactions, we enrolled 19 participants who developed allergic events within 2 h of vaccination and 13 controls, nonreactors. Using standard hemolysis assays, we demonstrated that sera from allergic participants induced stronger complement activation compared to nonallergic subjects following ex vivo vaccine exposure.

RESULTS

Vaccine-mediated complement activation correlated with anti-polyethelyne glycol (PEG) IgG (but not IgM) levels while anti-PEG IgE was undetectable in all subjects. Depletion of total IgG suppressed complement activation in select individuals. To investigate the effects of vaccine excipients on basophil function, we employed a validated indirect basophil activation test that stratified the allergic populations into high and low responders. Complement C3a and C5a receptor blockade in this system suppressed basophil response, providing strong evidence for complement involvement in vaccine-mediated basophil activation. Single-cell multiome analysis revealed differential expression of genes encoding the cytokine response and Toll-like receptor (TLR) pathways within the monocyte compartment. Differential chromatin accessibility for IL-13 and IL-1B genes was found in allergic and nonallergic participants, suggesting that in vivo, epigenetic modulation of mononuclear phagocyte immunophenotypes determines their subsequent functional responsiveness, contributing to the overall physiologic manifestation of vaccine reactions.

CONCLUSION

These findings provide insights into the mechanisms underlying allergic reactions to COVID-19 mRNA vaccines, which may be used for future vaccine strategies in individuals with prior history of allergies or reactions and reduce vaccine hesitancy.

Progress in multifactorial single-cell chromatin profiling methods

Chromatin states play a key role in shaping overall cellular states and fates. Building a complete picture of the functional state of chromatin in cells requires the co-detection of several distinct biochemical aspects. These span DNA methylation, chromatin accessibility, chromosomal conformation, histone posttranslational modifications, and more. While this certainly presents a challenging task, over the past few years many new and creative methods have been developed that now enable co-assay of these different aspects of chromatin at single cell resolution. This field is entering an exciting phase, where a confluence of technological improvements, decreased sequencing costs, and computational innovation are presenting new opportunities to dissect the diversity of chromatin states present in tissues, and how these states may influence gene regulation. In this review, I discuss the spectrum of current experimental approaches for multifactorial chromatin profiling, highlight some of the experimental and analytical challenges, as well as some areas for further innovation.

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.

A multiomic atlas of the aging hippocampus reveals molecular changes in response to environmental enrichment

Aging involves the deterioration of organismal function, leading to the emergence of multiple pathologies. Environmental stimuli, including lifestyle, can influence the trajectory of this process and may be used as tools in the pursuit of healthy aging. To evaluate the role of epigenetic mechanisms in this context, we have generated bulk tissue and single cell multi-omic maps of the male mouse dorsal hippocampus in young and old animals exposed to environmental stimulation in the form of enriched environments. We present a molecular atlas of the aging process, highlighting two distinct axes, related to inflammation and to the dysregulation of mRNA metabolism, at the functional RNA and protein level. Additionally, we report the alteration of heterochromatin domains, including the loss of bivalent chromatin and the uncovering of a heterochromatin-switch phenomenon whereby constitutive heterochromatin loss is partially mitigated through gains in facultative heterochromatin. Notably, we observed the multi-omic reversal of a great number of aging-associated alterations in the context of environmental enrichment, which was particularly linked to glial and oligodendrocyte pathways. In conclusion, our work describes the epigenomic landscape of environmental stimulation in the context of aging and reveals how lifestyle intervention can lead to the multi-layered reversal of aging-associated decline.

Pan-cancer single-cell dissection reveals phenotypically distinct B cell subtypes

Characterizing the compositional and phenotypic characteristics of tumor-infiltrating B cells (TIBs) is important for advancing our understanding of their role in cancer development. Here, we establish a comprehensive resource of human B cells by integrating single-cell RNA sequencing data of B cells from 649 patients across 19 major cancer types. We demonstrate substantial heterogeneity in their total abundance and subtype composition and observe immunoglobulin G (IgG)-skewness of antibody-secreting cell isotypes. Moreover, we identify stress-response memory B cells and tumor-associated atypical B cells (TAABs), two tumor-enriched subpopulations with prognostic potential, shared in a pan-cancer manner. In particular, TAABs, characterized by a high clonal expansion level and proliferative capacity as well as by close interactions with activated CD4 T cells in tumors, are predictive of immunotherapy response. Our integrative resource depicts distinct clinically relevant TIB subsets, laying a foundation for further exploration of functional commonality and diversity of B cells in cancer.

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

How disruptions to normal cell differentiation link to tumorigenesis remains incompletely understood. Wilms tumor, an embryonal tumor associated with disrupted organogenesis, often harbors mutations in epigenetic regulators, but their role in kidney development remains unexplored. Here, we show at single-cell resolution that a Wilms tumor-associated mutation in the histone acetylation reader ENL disrupts kidney differentiation in mice by rewiring the gene regulatory landscape. Mutant ENL promotes nephron progenitor commitment while restricting their differentiation by dysregulating transcription factors such as Hox clusters. It also induces abnormal progenitors that lose kidney-associated chromatin identity. Furthermore, mutant ENL alters the transcriptome and chromatin accessibility of stromal progenitors, resulting in hyperactivation of Wnt signaling. The impacts of mutant ENL on both nephron and stroma lineages lead to profound kidney developmental defects and postnatal mortality in mice. Notably, a small molecule inhibiting mutant ENL's histone acetylation binding activity largely reverses these defects. This study provides insights into how mutations in epigenetic regulators disrupt kidney development and suggests a potential therapeutic approach.

Histone lactylation-related genes correlate with the molecular patterns and functions of cancer-associated fibroblasts and have significant clinical implications in clear cell renal cell carcinoma

Recent research emphasised the indispensable role of histone lactylation in the activation of hepatic stellate cells. The VHL mutation is extremely common in clear cell renal cell carcinoma, which normally causes a metabolic shift in cancer cells and increases lactate production, eventually creating a lactate-enriched tumour microenvironment. Cancer-associated fibroblasts (CAFs) promote tumour progression, which is also vital in clear cell renal cell carcinoma. Therefore, this study investigated histone lactylation in CAFs and its impact on patient survival. Multiomics technology was employed to determine the role of histone lactylation-related genes in the evolution of CAFs which correlated with the function and molecular signatures of CAFs. The results suggested that TIMP1 was the hub gene of histone lactylation-related genes in clear cell renal cell carcinoma.

Robust reprogramming of glia into neurons by inhibition of Notch signaling and nuclear factor I (NFI) factors in adult mammalian retina

Generation of neurons through direct reprogramming has emerged as a promising therapeutic approach for treating neurodegenerative diseases. In this study, we present an efficient method for reprogramming retinal glial cells into neurons. By suppressing Notch signaling by disrupting either Rbpj or Notch1/2, we induced mature Müller glial cells to reprogram into bipolar- and amacrine-like neurons. We demonstrate that Rbpj directly activates both Notch effector genes and genes specific to mature Müller glia while indirectly repressing expression of neurogenic basic helix-loop-helix (bHLH) factors. Combined loss of function of Rbpj and Nfia/b/x resulted in conversion of nearly all Müller glia to neurons. Last, inducing Müller glial proliferation by overexpression of dominant-active Yap promotes neurogenesis in both Rbpj- and Nfia/b/x/Rbpj-deficient Müller glia. These findings demonstrate that Notch signaling and NFI factors act in parallel to inhibit neurogenic competence in mammalian Müller glia and help clarify potential strategies for regenerative therapies aimed at treating retinal dystrophies.

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.

Progenitors of distinct lineages shape the diversity of mature type 2 conventional dendritic cells

Conventional dendritic cells (cDC) are antigen-presenting cells comprising cDC1 and cDC2, responsible for priming naive CD8+ and CD4+ T cells, respectively. Recent studies have unveiled cDC2 heterogeneity and identified various cDC2 progenitors beyond the common DC progenitor (CDP), hinting at distinct cDC2 lineages. By generating Cd300ciCre-hCD2R26tdTomato reporter mice, we identified a bone marrow pro-cDC2 progenitor exclusively generating cDC2 in vitro and in vivo. Single-cell analyses and multiparametric flow cytometry demonstrated that pro-cDC2 encompasses myeloid-derived pre-cDC2 and lymphoid-derived plasmacytoid DC (pDC)-like precursors differentiating into a transcriptionally convergent cDC2 phenotype. Cd300c-traced cDC2 had distinct transcriptomic profiles, phenotypes, and tissue distributions compared with Ms4a3CreR26tdTomato lineage-traced DC3, a monocyte-DC progenitor (MDP)-derived subset that bypasses CDP. Mice with reduced Cd300c-traced cDC2 showed impaired humoral responses to T cell-dependent antigens. We conclude that progenitors of distinct lineages shape the diversity of mature cDC2 across tissues. Thus, ontogenesis may impact tissue immune responses.

scMultiome analysis identifies a single caudal hindbrain compartment in the developing zebrafish nervous system

BACKGROUND

A key step in nervous system development involves the coordinated control of neural progenitor specification and positioning. A long-standing model for the vertebrate CNS postulates that transient anatomical compartments - known as neuromeres - function to position neural progenitors along the embryonic anteroposterior neuraxis. Such neuromeres are apparent in the embryonic hindbrain - that contains six rhombomeres with morphologically apparent boundaries - but other neuromeres lack clear morphological boundaries and have instead been defined by different criteria, such as differences in gene expression patterns and the outcomes of transplantation experiments. Accordingly, the caudal hindbrain (CHB) posterior to rhombomere (r) 6 has been variably proposed to contain from two to five 'pseudo-rhombomeres', but the lack of comprehensive molecular data has precluded a detailed definition of such structures.

METHODS

We used single-cell Multiome analysis, which allows simultaneous characterization of gene expression and chromatin state of individual cell nuclei, to identify and characterize CHB progenitors in the developing zebrafish CNS.

RESULTS

We identified CHB progenitors as a transcriptionally distinct population, that also possesses a unique profile of accessible transcription factor binding motifs, relative to both r6 and the spinal cord. This CHB population can be subdivided along its dorsoventral axis based on molecular characteristics, but we do not find any molecular evidence that it contains multiple pseudo-rhombomeres. We further observe that the CHB is closely related to r6 at the earliest embryonic stages, but becomes more divergent over time, and that it is defined by a unique gene regulatory network.

CONCLUSIONS

We conclude that the early CHB represents a single neuromere compartment that cannot be molecularly subdivided into pseudo-rhombomeres and that it may share an embryonic origin with r6.

An atlas of transcribed enhancers across helper T cell diversity for decoding human diseases

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4+ T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

Transcription factor dynamics, oscillation, and functions in human enteroendocrine cell differentiation

Enteroendocrine cells (EECs) secrete serotonin (enterochromaffin [EC] cells) or specific peptide hormones (non-EC cells) that serve vital metabolic functions. The basis for terminal EEC diversity remains obscure. By forcing activity of the transcription factor (TF) NEUROG3 in 2D cultures of human intestinal stem cells, we replicated physiologic EEC differentiation and examined transcriptional and cis-regulatory dynamics that culminate in discrete cell types. Abundant EEC precursors expressed stage-specific genes and TFs. Before expressing pre-terminal NEUROD1, post-mitotic precursors oscillated between transcriptionally distinct ASCL1+ and HES6hi cell states. Loss of either factor accelerated EEC differentiation substantially and disrupted EEC individuality; ASCL1 or NEUROD1 deficiency had opposing consequences on EC and non-EC cell features. These TFs mainly bind cis-elements that are accessible in undifferentiated stem cells, and they tailor subsequent expression of TF combinations that underlie discrete EEC identities. Thus, early TF oscillations retard EEC maturation to enable accurate diversity within a medically important cell lineage.

ezSingleCell: an integrated one-stop single-cell and spatial omics analysis platform for bench scientists

ezSingleCell is an interactive and easy-to-use application for analysing various single-cell and spatial omics data types without requiring prior programing knowledge. It combines the best-performing publicly available methods for in-depth data analysis, integration, and interactive data visualization. ezSingleCell consists of five modules, each designed to be a comprehensive workflow for one data type or task. In addition, ezSingleCell allows crosstalk between different modules within a unified interface. Acceptable input data can be in a variety of formats while the output consists of publication ready figures and tables. In-depth manuals and video tutorials are available to guide users on the analysis workflows and parameter adjustments to suit their study aims. ezSingleCell's streamlined interface can analyse a standard scRNA-seq dataset of 3000 cells in less than five minutes. ezSingleCell is available in two forms: an installation-free web application ( https://immunesinglecell.org/ezsc/ ) or a software package with a shinyApp interface ( https://github.com/JinmiaoChenLab/ezSingleCell2 ) for offline analysis.

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

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

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

Type 2 and type 1 diabetes (T2D, T1D) exhibit sex differences in insulin secretion, the mechanisms of which are unknown. We examined sex differences in human pancreatic islets from 52 donors with and without T2D combining single cell RNA-seq (scRNA-seq), single nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq), hormone secretion, and bioenergetics. In nondiabetic (ND) donors, sex differences in islet cells gene accessibility and expression predominantly involved sex chromosomes. Islets from T2D donors exhibited similar sex differences in sex chromosomes differentially expressed genes (DEGs), but also exhibited sex differences in autosomal genes. Comparing β cells from T2D vs. ND donors, gene enrichment of female β cells showed suppression in mitochondrial respiration, while male β cells exhibited suppressed insulin secretion. Thus, although sex differences in gene accessibility and expression of ND β cells predominantly affect sex chromosomes, the transition to T2D reveals sex differences in autosomes highlighting mitochondrial failure in females.

Single-nucleus multi-omics of Parkinson's disease reveals a glutamatergic neuronal subtype susceptible to gene dysregulation via alteration of transcriptional networks

The genetic architecture of Parkinson's disease (PD) is complex and multiple brain cell subtypes are involved in the neuropathological progression of the disease. Here we aimed to advance our understanding of PD genetic complexity at a cell subtype precision level. Using parallel single-nucleus (sn)RNA-seq and snATAC-seq analyses we simultaneously profiled the transcriptomic and chromatin accessibility landscapes in temporal cortex tissues from 12 PD compared to 12 control subjects at a granular single cell resolution. An integrative bioinformatic pipeline was developed and applied for the analyses of these snMulti-omics datasets. The results identified a subpopulation of cortical glutamatergic excitatory neurons with remarkably altered gene expression in PD, including differentially-expressed genes within PD risk loci identified in genome-wide association studies (GWAS). This was the only neuronal subtype showing significant and robust overexpression of SNCA. Further characterization of this neuronal-subpopulation showed upregulation of specific pathways related to axon guidance, neurite outgrowth and post-synaptic structure, and downregulated pathways involved in presynaptic organization and calcium response. Additionally, we characterized the roles of three molecular mechanisms in governing PD-associated cell subtype-specific dysregulation of gene expression: (1) changes in cis-regulatory element accessibility to transcriptional machinery; (2) changes in the abundance of master transcriptional regulators, including YY1, SP3, and KLF16; (3) candidate regulatory variants in high linkage disequilibrium with PD-GWAS genomic variants impacting transcription factor binding affinities. To our knowledge, this study is the first and the most comprehensive interrogation of the multi-omics landscape of PD at a cell-subtype resolution. Our findings provide new insights into a precise glutamatergic neuronal cell subtype, causal genes, and non-coding regulatory variants underlying the neuropathological progression of PD, paving the way for the development of cell- and gene-targeted therapeutics to halt disease progression as well as genetic biomarkers for early preclinical diagnosis.

A unified model-based framework for doublet or multiplet detection in single-cell multiomics data

Droplet-based single-cell sequencing techniques rely on the fundamental assumption that each droplet encapsulates a single cell, enabling individual cell omics profiling. However, the inevitable issue of multiplets, where two or more cells are encapsulated within a single droplet, can lead to spurious cell type annotations and obscure true biological findings. The issue of multiplets is exacerbated in single-cell multiomics settings, where integrating cross-modality information for clustering can inadvertently promote the aggregation of multiplet clusters and increase the risk of erroneous cell type annotations. Here, we propose a compound Poisson model-based framework for multiplet detection in single-cell multiomics data. Leveraging experimental cell hashing results as the ground truth for multiplet status, we conducted trimodal DOGMA-seq experiments and generated 17 benchmarking datasets from two tissues, involving a total of 280,123 droplets. We demonstrated that the proposed method is an essential tool for integrating cross-modality multiplet signals, effectively eliminating multiplet clusters in single-cell multiomics data-a task at which the benchmarked single-omics methods proved inadequate.

Direct neuronal reprogramming of mouse astrocytes is associated with multiscale epigenome remodeling and requires Yy1

Direct neuronal reprogramming is a promising approach to regenerate neurons from local glial cells. However, mechanisms of epigenome remodeling and co-factors facilitating this process are unclear. In this study, we combined single-cell multiomics with genome-wide profiling of three-dimensional nuclear architecture and DNA methylation in mouse astrocyte-to-neuron reprogramming mediated by Neurogenin2 (Ngn2) and its phosphorylation-resistant form (PmutNgn2), respectively. We show that Ngn2 drives multilayered chromatin remodeling at dynamic enhancer-gene interaction sites. PmutNgn2 leads to higher reprogramming efficiency and enhances epigenetic remodeling associated with neuronal maturation. However, the differences in binding sites or downstream gene activation cannot fully explain this effect. Instead, we identified Yy1, a transcriptional co-factor recruited by direct interaction with Ngn2 to its target sites. Upon deletion of Yy1, activation of neuronal enhancers, genes and ultimately reprogramming are impaired without affecting Ngn2 binding. Thus, our work highlights the key role of interactors of proneural factors in direct neuronal reprogramming.

ID factors regulate the ability of Müller glia to become proliferating neurogenic progenitor-like cells

The purpose of this study was to investigate how ID factors regulate the ability of Müller glia (MG) to reprogram into proliferating MG-derived progenitor cells (MGPCs) in the chick retina. We found that ID1 is transiently expressed by maturing MG (mMG), whereas ID4 is maintained in mMG in embryonic retinas. In mature retinas, ID4 was prominently expressed by resting MG, but following retinal damage ID4 was rapidly upregulated and then downregulated in MGPCs. By contrast, ID1, ID2, and ID3 were low in resting MG and then upregulated in MGPCs. Inhibition of ID factors following retinal damage decreased numbers of proliferating MGPCs. Inhibition of IDs, after MGPC proliferation, significantly increased numbers of progeny that differentiated as neurons. In damaged or undamaged retinas inhibition of IDs increased levels of p21Cip1 in MG. In response to damage or insulin+FGF2 levels of CDKN1A message and p21Cip1 protein were decreased, absent in proliferating MGPCs, and elevated in MG returning to a resting phenotype. Inhibition of notch- or gp130/Jak/Stat-signaling in damaged retinas increased levels of ID4 but not p21Cip1 in MG. Although ID4 is the predominant isoform expressed by MG in the chick retina, id1 and id2a are predominantly expressed by resting MG and downregulated in activated MG and MGPCs in zebrafish retinas. We conclude that ID factors have a significant impact on regulating the responses of MG to retinal damage, controlling the ability of MG to proliferate by regulating levels of p21Cip1, and suppressing the neurogenic potential of MGPCs.

Alcohol-induced epigenetic changes prevent fibrosis resolution after alcohol cessation in miceresolution

BACKGROUND AND AIMS

Alcohol-associated liver disease is a major cause of alcohol-associated mortality. Recently, we identified hepatic demethylases lysine demethylase (KDM)5B and KDM5C as important epigenetic regulators of alcohol response in the liver. In this study, we aimed to investigate the role of KDM5 demethylases in alcohol-associated liver disease resolution.

APPROACH AND RESULTS

We showed that alcohol-induced liver steatosis rapidly resolved after alcohol cessation. In contrast, fibrosis persisted in the liver for up to 8 weeks after the end of alcohol exposure. Defects in fibrosis resolution were in part due to alcohol-induced KDM5B and KDM5C-dependent epigenetic changes in hepatocytes. Using cell-type-specific knockout mice, we found that adeno-associated virus-mediated knockout of KDM5B and KDM5C demethylases in hepatocytes at the time of alcohol withdrawal promoted fibrosis resolution. Single-cell ATAC sequencing analysis showed that during alcohol-associated liver disease resolution epigenetic cell states largely reverted to control conditions. In addition, we found unique epigenetic cell states distinct from both control and alcohol states and identified associated transcriptional regulators, including liver X receptor (LXR) alpha (α). In vitro and in vivo analysis confirmed that knockout of KDM5B and KDM5C demethylases promoted LXRα activity, likely through regulation of oxysterol biosynthesis, and this activity was critical for the fibrosis resolution process. Reduced LXR activity by small molecule inhibitors prevented fibrosis resolution in KDM5-deficient mice.

CONCLUSIONS

In summary, KDM5B and KDM5C demethylases prevent liver fibrosis resolution after alcohol cessation in part through suppression of LXR activity.

Non-stem cell lineages as an alternative origin of intestinal tumorigenesis in the context of inflammation

According to conventional views, colon cancer originates from stem cells. However, inflammation, a key risk factor for colon cancer, has been shown to suppress intestinal stemness. Here, we used Paneth cells as a model to assess the capacity of differentiated lineages to trigger tumorigenesis in the context of inflammation in mice. Upon inflammation, Paneth cell-specific Apc mutations led to intestinal tumors reminiscent not only of those arising in patients with inflammatory bowel disease, but also of a larger fraction of human sporadic colon cancers. The latter is possibly because of the inflammatory consequences of western-style dietary habits, a major colon cancer risk factor. Machine learning methods designed to predict the cell-of-origin of cancer from patient-derived tumor samples confirmed that, in a substantial fraction of sporadic cases, the origins of colon cancer reside in secretory lineages and not in stem cells.

Single-cell epigenomic reconstruction of developmental trajectories from pluripotency in human neural organoid systems

Cell fate progression of pluripotent progenitors is strictly regulated, resulting in high human cell diversity. Epigenetic modifications also orchestrate cell fate restriction. Unveiling the epigenetic mechanisms underlying human cell diversity has been difficult. In this study, we use human brain and retina organoid models and present single-cell profiling of H3K27ac, H3K27me3 and H3K4me3 histone modifications from progenitor to differentiated neural fates to reconstruct the epigenomic trajectories regulating cell identity acquisition. We capture transitions from pluripotency through neuroepithelium to retinal and brain region and cell type specification. Switching of repressive and activating epigenetic modifications can precede and predict cell fate decisions at each stage, providing a temporal census of gene regulatory elements and transcription factors. Removing H3K27me3 at the neuroectoderm stage disrupts fate restriction, resulting in aberrant cell identity acquisition. Our single-cell epigenome-wide map of human neural organoid development serves as a blueprint to explore human cell fate determination.

Similarity Network Analysis of the Adaptive Immune Response in the Proximal Airway

OBJECTIVES

Recent immunologic study of the adaptive immune repertoire in the subglottic airway demonstrated high-frequency T cell clones that do not overlap between individuals. However, the anatomic distribution and antigenic target of the T cell repertoire in the proximal airway mucosa remain unresolved.

METHODS

Single-cell RNA sequencing of matched scar and unaffected mucosa from idiopathic subglottic stenosis patients (iSGS, n = 32) was performed and compared with airway mucosa from healthy controls (n = 10). T cell receptor (TCR) sequences were interrogated via similarity network analysis to explore antigenic targets using the published algorithm: Grouping of Lymphocyte Interactions by Paratope Hotspots (GLIPH2).

RESULTS

The mucosal T cell repertoire in healthy control airways consisted of highly expressed T cell clones conserved across anatomic subsites (trachea, bronchi, bronchioles, and lung). In iSGS, high-frequency clones were equally represented in both scar and adjacent non-scar tissue. Significant differences in repertoire structure between iSGS scar and unaffected mucosa was observed, driven by unique low-frequency clones. GLIPH2 results suggest low-frequency clones share targets between multiple iSGS patients.

CONCLUSION

Healthy airway mucosa has a highly conserved T cell repertoire across multiple anatomic subsites. Similarly, iSGS patients have highly expressed T cell clones present in both scar and unaffected mucosa. iSGS airway scar possesses an abundance of less highly expanded clones with predicted antigen targets shared between patients. Interrogation of these shared motifs suggests abundant adaptive immunity to viral targets in iSGS airway scar. These results provide insight into disease pathogenesis and illuminate new treatment strategies in iSGS.

LEVEL OF EVIDENCE

NA Laryngoscope, 134:3245-3252, 2024.

Single-nucleus multiomics reveals the disrupted regulatory programs in three brain regions of sporadic early-onset Alzheimer's disease

Sporadic early-onset Alzheimer's disease (sEOAD) represents a significant but less-studied subtype of Alzheimer's disease (AD). Here, we generated a single-nucleus multiome atlas derived from the postmortem prefrontal cortex, entorhinal cortex, and hippocampus of nine individuals with or without sEOAD. Comprehensive analyses were conducted to delineate cell type-specific transcriptomic changes and linked candidate cis- regulatory elements (cCREs) across brain regions. We prioritized seven conservative transcription factors in glial cells in multiple brain regions, including RFX4 in astrocytes and IKZF1 in microglia, which are implicated in regulating sEOAD-associated genes. Moreover, we identified the top 25 altered intercellular signaling between glial cells and neurons, highlighting their regulatory potential on gene expression in receiver cells. We reported 38 cCREs linked to sEOAD-associated genes overlapped with late-onset AD risk loci, and sEOAD cCREs enriched in neuropsychiatric disorder risk loci. This atlas helps dissect transcriptional and chromatin dynamics in sEOAD, providing a key resource for AD research.

REUNION: transcription factor binding prediction and regulatory association inference from single-cell multi-omics data

MOTIVATION

Profiling of gene expression and chromatin accessibility by single-cell multi-omics approaches can help to systematically decipher how transcription factors (TFs) regulate target gene expression via cis-region interactions. However, integrating information from different modalities to discover regulatory associations is challenging, in part because motif scanning approaches miss many likely TF binding sites.

RESULTS

We develop REUNION, a framework for predicting genome-wide TF binding and cis-region-TF-gene "triplet" regulatory associations using single-cell multi-omics data. The first component of REUNION, Unify, utilizes information theory-inspired complementary score functions that incorporate TF expression, chromatin accessibility, and target gene expression to identify regulatory associations. The second component, Rediscover, takes Unify estimates as input for pseudo semi-supervised learning to predict TF binding in accessible genomic regions that may or may not include detected TF motifs. Rediscover leverages latent chromatin accessibility and sequence feature spaces of the genomic regions, without requiring chromatin immunoprecipitation data for model training. Applied to peripheral blood mononuclear cell data, REUNION outperforms alternative methods in TF binding prediction on average performance. In particular, it recovers missing region-TF associations from regions lacking detected motifs, which circumvents the reliance on motif scanning and facilitates discovery of novel associations involving potential co-binding transcriptional regulators. Newly identified region-TF associations, even in regions lacking a detected motif, improve the prediction of target gene expression in regulatory triplets, and are thus likely to genuinely participate in the regulation.

AVAILABILITY AND IMPLEMENTATION

All source code is available at https://github.com/yangymargaret/REUNION.

The activity of early-life gene regulatory elements is hijacked in aging through pervasive AP-1-linked chromatin opening

A mechanistic connection between aging and development is largely unexplored. Through profiling age-related chromatin and transcriptional changes across 22 murine cell types, analyzed alongside previous mouse and human organismal maturation datasets, we uncovered a transcription factor binding site (TFBS) signature common to both processes. Early-life candidate cis-regulatory elements (cCREs), progressively losing accessibility during maturation and aging, are enriched for cell-type identity TFBSs. Conversely, cCREs gaining accessibility throughout life have a lower abundance of cell identity TFBSs but elevated activator protein 1 (AP-1) levels. We implicate TF redistribution toward these AP-1 TFBS-rich cCREs, in synergy with mild downregulation of cell identity TFs, as driving early-life cCRE accessibility loss and altering developmental and metabolic gene expression. Such remodeling can be triggered by elevating AP-1 or depleting repressive H3K27me3. We propose that AP-1-linked chromatin opening drives organismal maturation by disrupting cell identity TFBS-rich cCREs, thereby reprogramming transcriptome and cell function, a mechanism hijacked in aging through ongoing chromatin opening.

Systematic single-cell analysis reveals dynamic control of transposable element activity orchestrating the endothelial-to-hematopoietic transition

BACKGROUND

The endothelial-to-hematopoietic transition (EHT) process during definitive hematopoiesis is highly conserved in vertebrates. Stage-specific expression of transposable elements (TEs) has been detected during zebrafish EHT and may promote hematopoietic stem cell (HSC) formation by activating inflammatory signaling. However, little is known about how TEs contribute to the EHT process in human and mouse.

RESULTS

We reconstructed the single-cell EHT trajectories of human and mouse and resolved the dynamic expression patterns of TEs during EHT. Most TEs presented a transient co-upregulation pattern along the conserved EHT trajectories, coinciding with the temporal relaxation of epigenetic silencing systems. TE products can be sensed by multiple pattern recognition receptors, triggering inflammatory signaling to facilitate HSC emergence. Interestingly, we observed that hypoxia-related signals were enriched in cells with higher TE expression. Furthermore, we constructed the hematopoietic cis-regulatory network of accessible TEs and identified potential TE-derived enhancers that may boost the expression of specific EHT marker genes.

CONCLUSIONS

Our study provides a systematic vision of how TEs are dynamically controlled to promote the hematopoietic fate decisions through transcriptional and cis-regulatory networks, and pre-train the immunity of nascent HSCs.

A cis-regulatory module underlies retinal ganglion cell genesis and axonogenesis

Atoh7 is transiently expressed in retinal progenitor cells (RPCs) and is required for retinal ganglion cell (RGC) differentiation. In humans, a deletion in a distal non-coding regulatory region upstream of ATOH7 is associated with optic nerve atrophy and blindness. Here, we functionally interrogate the significance of the Atoh7 regulatory landscape to retinogenesis in mice. Deletion of the Atoh7 enhancer structure leads to RGC deficiency, optic nerve hypoplasia, and retinal blood vascular abnormalities, phenocopying inactivation of Atoh7. Further, loss of the Atoh7 remote enhancer impacts ipsilaterally projecting RGCs and disrupts proper axonal projections to the visual thalamus. Deletion of the Atoh7 remote enhancer is also associated with the dysregulation of axonogenesis genes, including the derepression of the axon repulsive cue Robo3. Our data provide insights into how Atoh7 enhancer elements function to promote RGC development and optic nerve formation and highlight a key role of Atoh7 in the transcriptional control of axon guidance molecules.

Multi-omic analysis of SDHB-deficient pheochromocytomas and paragangliomas identifies metastasis and treatment-related molecular profiles

Hereditary SDHB-mutant pheochromocytomas (PC) and paragangliomas (PG) are rare tumours with a high propensity to metastasize although their clinical behaviour is unpredictable. To characterize the genomic landscape of these tumours and identify metastasis biomarkers, we performed multi-omic analysis on 94 tumours from 79 patients using seven molecular methods. Sympathetic (chromaffin cell) and parasympathetic (non-chromaffin cell) PCPG had distinct molecular profiles reflecting their cell-of-origin and biochemical profile. TERT and ATRX-alterations were associated with metastatic PCPG and these tumours had an increased mutation load, and distinct transcriptional and telomeric features. Most PCPG had quiet genomes with some rare co-operative driver events observed, including EPAS1/HIF-2α mutations. Two mechanisms of acquired resistance to DNA alkylating chemotherapies were also detected - MGMT overexpression and mismatch repair-deficiency causing hypermutation. Our comprehensive multi-omic analysis of SDHB-mutant PCPG therefore identified features of metastatic disease and treatment response, expanding our understanding of these rare neuroendocrine tumours.

B cells modulate lung antiviral inflammatory responses via the neurotransmitter acetylcholine

The rapid onset of innate immune defenses is critical for early control of viral replication in an infected host, yet it can also lead to irreversible tissue damage, especially in the respiratory tract. Intricate regulatory mechanisms must exist that modulate inflammation, while controlling the infection. Here, B cells expressing choline acetyl transferase (ChAT), an enzyme required for production of the metabolite and neurotransmitter acetylcholine (ACh) are identified as such regulators of the immediate early response to influenza A virus. Lung tissue ChAT + B cells are shown to interact with a7 nicotinic Ach receptor-expressing lung interstitial macrophages in mice within 24h of infection to control their production of TNFa, shifting the balance towards reduced inflammation at the cost of enhanced viral replication. Thus, innate-stimulated B cells are key participants of an immediate-early regulatory cascade that controls lung tissue damage after viral infection.

Differential impact of quiescent non-coding loci on chromatin entropy

Non-coding regions of the human genome are important for functional regulations, but their mechanisms remain elusive. We used machine learning to guide a CRISPR screening on hubs (i.e. non-coding loci forming many 3D contacts) and significantly increased the discovery rate of hubs essential for cell growth. We found no clear genetic or epigenetic differences between essential and nonessential hubs, but we observed that some neighboring hubs in the linear genome have distinct spatial contacts and opposite effects on cell growth. One such pair in an epigenetically quiescent region showed different impacts on gene expression, chromatin accessibility and chromatin organization. We also found that deleting the essential hub altered the genetic network activity and increased the entropy of chromatin accessibility, more severe than that caused by deletion of the nonessential hub, suggesting that they are critical for maintaining an ordered chromatin structure. Our study reveals new insights into the system-level roles of non-coding regions in the human genome.

SARS-CoV-2 inflammation durably imprints memory CD4 T cells

Memory CD4 T cells are critical to human immunity, yet it is unclear whether viral inflammation during memory formation has long-term consequences. Here, we compared transcriptional and epigenetic landscapes of Spike (S)-specific memory CD4 T cells in 24 individuals whose first exposure to S was via SARS-CoV-2 infection or mRNA vaccination. Nearly 2 years after memory formation, S-specific CD4 T cells established by infection remained enriched for transcripts related to cytotoxicity and for interferon-stimulated genes, likely because of a chromatin accessibility landscape altered by inflammation. Moreover, S-specific CD4 T cells primed by infection had reduced proliferative capacity in vitro relative to vaccine-primed cells. Furthermore, the transcriptional state of S-specific memory CD4 T cells was minimally altered by booster immunization and/or breakthrough infection. Thus, infection-associated inflammation durably imprints CD4 T cell memory, which affects the function of these cells and may have consequences for long-term immunity.

Single-cell multi-omics analysis reveals cooperative transcription factors for gene regulation in oligodendrocytes

Oligodendrocytes are the myelinating cells within the central nervous system. Many oligodendrocyte genes have been associated with brain disorders. However, how transcription factors (TFs) cooperate for gene regulation in oligodendrocytes remains largely uncharacterized. To address this, we integrated scRNA-seq and scATAC-seq data to identify the cooperative TFs that co-regulate the target gene (TG) expression in oligodendrocytes. First, we identified co- binding TF pairs whose binding sites overlapped in oligodendrocyte-specific regulatory regions. Second, we trained a deep learning model to predict the expression level of each TG using the expression levels of co-binding TFs. Third, using the trained models, we computed the TF importance and TF-TF interaction scores for predicting TG expression by the Shapley interaction scores. We found that the co-binding TF pairs involving known important TF pairs for oligodendrocyte differentiation, such as SOX10-TCF12, SOX10-MYRF, and SOX10-OLIG2, exhibited significantly higher Shapley scores than others (t-test, p-value < 1e-4). Furthermore, we identified 153 oligodendrocyte-associated eQTLs that reside in oligodendrocyte-specific enhancers or promoters where their eGenes (TGs) are regulated by cooperative TFs, suggesting potential novel regulatory roles from genetic variants. We also experimentally validated some identified TF pairs such as SOX10-OLIG2 and SOX10-NKX2.2 by co-enrichment analysis, using ChIP-seq data from rat peripheral nerve.

Protocol for the isolation and single-nuclei multiomic analyses of the human fetal epicardium

The characterization of cell populations that reside in the outer layer of the heart has been hindered by difficulties in their isolation. Here, we present a protocol for isolation and single-nuclei multiomic analyses of the human fetal epicardium. We describe steps for microdissection, isolation, and enrichment of epicardial cells by mechanical dissociations and direct lysis. We then detail procedures for integrating transcriptome and chromatin accessibility datasets. This approach allows the analysis of diverse cell populations, marked by unique cis-regulatory elements. For complete details on the use and execution of this protocol, please refer to Travisano et al.1.

Multiplexed single-cell characterization of alternative polyadenylation regulators

Most mammalian genes have multiple polyA sites, representing a substantial source of transcript diversity regulated by the cleavage and polyadenylation (CPA) machinery. To better understand how these proteins govern polyA site choice, we introduce CPA-Perturb-seq, a multiplexed perturbation screen dataset of 42 CPA regulators with a 3' scRNA-seq readout that enables transcriptome-wide inference of polyA site usage. We develop a framework to detect perturbation-dependent changes in polyadenylation and characterize modules of co-regulated polyA sites. We find groups of intronic polyA sites regulated by distinct components of the nuclear RNA life cycle, including elongation, splicing, termination, and surveillance. We train and validate a deep neural network (APARENT-Perturb) for tandem polyA site usage, delineating a cis-regulatory code that predicts perturbation response and reveals interactions between regulatory complexes. Our work highlights the potential for multiplexed single-cell perturbation screens to further our understanding of post-transcriptional regulation.

Massively parallel in vivo Perturb-seq reveals cell-type-specific transcriptional networks in cortical development

Leveraging AAVs' versatile tropism and labeling capacity, we expanded the scale of in vivo CRISPR screening with single-cell transcriptomic phenotyping across embryonic to adult brains and peripheral nervous systems. Through extensive tests of 86 vectors across AAV serotypes combined with a transposon system, we substantially amplified labeling efficacy and accelerated in vivo gene delivery from weeks to days. Our proof-of-principle in utero screen identified the pleiotropic effects of Foxg1, highlighting its tight regulation of distinct networks essential for cell fate specification of Layer 6 corticothalamic neurons. Notably, our platform can label >6% of cerebral cells, surpassing the current state-of-the-art efficacy at <0.1% by lentivirus, to achieve analysis of over 30,000 cells in one experiment and enable massively parallel in vivo Perturb-seq. Compatible with various phenotypic measurements (single-cell or spatial multi-omics), it presents a flexible approach to interrogate gene function across cell types in vivo, translating gene variants to their causal function.

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.

GATA2 heterozygosity causes an epigenetic feedback mechanism resulting in myeloid and erythroid dysplasia

The transcription factor GATA2 has a pivotal role in haematopoiesis. Heterozygous germline GATA2 mutations result in a syndrome characterized by immunodeficiency, bone marrow failure and predispositions to myelodysplastic syndrome (MDS) and acute myeloid leukaemia. Clinical symptoms in these patients are diverse and mechanisms driving GATA2-related phenotypes are largely unknown. To explore the impact of GATA2 haploinsufficiency on haematopoiesis, we generated a zebrafish model carrying a heterozygous mutation of gata2b (gata2b+/-), an orthologue of GATA2. Morphological analysis revealed myeloid and erythroid dysplasia in gata2b+/- kidney marrow. Because Gata2b could affect both transcription and chromatin accessibility during lineage differentiation, this was assessed by single-cell (sc) RNA-seq and single-nucleus (sn) ATAC-seq. Sn-ATAC-seq showed that the co-accessibility between the transcription start site (TSS) and a -3.5-4.1 kb putative enhancer was more robust in gata2b+/- zebrafish HSPCs compared to wild type, increasing gata2b expression and resulting in higher genome-wide Gata2b motif use in HSPCs. As a result of increased accessibility of the gata2b locus, gata2b+/- chromatin was also more accessible during lineage differentiation. scRNA-seq data revealed myeloid differentiation defects, that is, impaired cell cycle progression, reduced expression of cebpa and cebpb and increased signatures of ribosome biogenesis. These data also revealed a differentiation delay in erythroid progenitors, aberrant proliferative signatures and down-regulation of Gata1a, a master regulator of erythropoiesis, which worsened with age. These findings suggest that cell-intrinsic compensatory mechanisms, needed to obtain normal levels of Gata2b in heterozygous HSPCs to maintain their integrity, result in aberrant lineage differentiation, thereby representing a critical step in the predisposition to MDS.

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

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

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

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

CAraCAl: CAMML with the integration of chromatin accessibility

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Leveraging single cell multiomic analyses to identify factors that drive human chondrocyte cell fate

Cartilage plays a crucial role in skeletal development and function, and abnormal development contributes to genetic and age-related skeletal disease. To better understand how human cartilage develops in vivo , we jointly profiled the transcriptome and open chromatin regions in individual nuclei recovered from distal femurs at 2 fetal timepoints. We used these multiomic data to identify transcription factors expressed in distinct chondrocyte subtypes, link accessible regulatory elements with gene expression, and predict transcription factor-based regulatory networks that are important for growth plate or epiphyseal chondrocyte differentiation. We developed a human pluripotent stem cell platform for interrogating the function of predicted transcription factors during chondrocyte differentiation and used it to test NFATC2 . We expect new regulatory networks we uncovered using multiomic data to be important for promoting cartilage health and treating disease, and our platform to be a useful tool for studying cartilage development in vitro .

Statement of Significance

The identity and integrity of the articular cartilage lining our joints are crucial to pain-free activities of daily living. Here we identified a gene regulatory landscape of human chondrogenesis at single cell resolution, which is expected to open new avenues of research aimed at mitigating cartilage diseases that affect hundreds of millions of individuals world-wide.