Massively parallel single-cell chromatin landscapes of human immune cell development and intratumoral T cell exhaustion

CellMag-CARWash: A High Throughput Droplet Microfluidic Device for Live Cell Isolation and Single Cell Applications

The recent push toward understanding an individual cell's behavior and identifying cellular heterogeneity has created an unmet need for technologies that can probe live cells at the single-cell level. Cells within a population are known to exhibit heterogeneous responses to environmental cues. These differences can lead to varied cellular states, behavior, and responses to therapeutics. Techniques are needed that are not only capable of processing and analyzing cellular populations at the single cell level, but also have the ability to isolate specific cell populations from a complex sample at high throughputs. The new CellMag-Coalesce-Attract-Resegment Wash (CellMag-CARWash) system combines positive magnetic selection with droplet microfluidic devices to isolate cells of interest from a mixture with >93% purity and incorporate treatments within individual droplets to observe single cell biological responses. This workflow is shown to be capable of probing the single cell extracellular vesicle (EV) secretion of MCF7 GFP cells. This article reports the first measurement of β-Estradiol's effect on EV secretion from MCF7 cells at the single cell level. Single cell processing revealed that MCF7 GFP cells possess a heterogeneous response to β-Estradiol stimulation with a 1.8-fold increase relative to the control.

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

INTRODUCTION

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

OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Discrete latent embedding of single-cell chromatin accessibility sequencing data for uncovering cell heterogeneity

Single-cell epigenomic data has been growing continuously at an unprecedented pace, but their characteristics such as high dimensionality and sparsity pose substantial challenges to downstream analysis. Although deep learning models-especially variational autoencoders-have been widely used to capture low-dimensional feature embeddings, the prevalent Gaussian assumption somewhat disagrees with real data, and these models tend to struggle to incorporate reference information from abundant cell atlases. Here we propose CASTLE, a deep generative model based on the vector-quantized variational autoencoder framework to extract discrete latent embeddings that interpretably characterize single-cell chromatin accessibility sequencing data. We validate the performance and robustness of CASTLE for accurate cell-type identification and reasonable visualization compared with state-of-the-art methods. We demonstrate the advantages of CASTLE for effective incorporation of existing massive reference datasets in a weakly supervised or supervised manner. We further demonstrate CASTLE's capacity for intuitively distilling cell-type-specific feature spectra that unveil cell heterogeneity and biological implications quantitatively.

A blueprint for tumor-infiltrating B cells across human cancers

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

Chromatin accessibility during human first-trimester neurodevelopment

The human brain develops through a tightly organized cascade of patterning events, induced by transcription factor expression and changes in chromatin accessibility. Although gene expression across the developing brain has been described at single-cell resolution1, similar atlases of chromatin accessibility have been primarily focused on the forebrain2-4. Here we describe chromatin accessibility and paired gene expression across the entire developing human brain during the first trimester (6-13 weeks after conception). We defined 135 clusters and used multiomic measurements to link candidate cis-regulatory elements to gene expression. The number of accessible regions increased both with age and along neuronal differentiation. Using a convolutional neural network, we identified putative functional transcription factor-binding sites in enhancers characterizing neuronal subtypes. We applied this model to cis-regulatory elements linked to ESRRB to elucidate its activation mechanism in the Purkinje cell lineage. Finally, by linking disease-associated single nucleotide polymorphisms to cis-regulatory elements, we validated putative pathogenic mechanisms in several diseases and identified midbrain-derived GABAergic neurons as being the most vulnerable to major depressive disorder-related mutations. Our findings provide a more detailed view of key gene regulatory mechanisms underlying the emergence of brain cell types during the first trimester and a comprehensive reference for future studies related to human neurodevelopment.

Edge-relational window-attentional graph neural network for gene expression prediction in spatial transcriptomics analysis

Spatial transcriptomics (ST), containing gene expression with fine-grained (i.e., different windows) spatial location within tissue samples, has become vital in developing innovative treatments. Traditional ST technology, however, rely on costly specialized commercial equipment. Addressing this, our article aims to creates a cost-effective, virtual ST approach using standard tissue images for gene expression prediction, eliminating the need for expensive equipment. Conventional approaches in this field often overlook the long-distance spatial dependencies between different sample windows or need prior gene expression data. To overcome these limitations, we propose the Edge-Relational Window-Attentional Network (ErwaNet), enhancing gene prediction by capturing both local interactions and global structural information from tissue images, without prior gene expression data. ErwaNet innovatively constructs heterogeneous graphs to model local window interactions and incorporates an attention mechanism for global information analysis. This dual framework not only provides a cost-effective solution for gene expression predictions but also obviates the necessity of prior knowledge gene expression information, a significant advantage in the field of cancer research where it enables a more efficient and accessible analytical paradigm. ErwaNet stands out as a prior-free and easy-to-implement Graph Convolution Network (GCN) method for predicting gene expression from tissue images. Evaluation of the two public breast cancer datasets shows that ErwaNet, without additional information, outperforms the state-of-the-art (SOTA) methods. Code is available at https://github.com/biyecc/ErwaNet.

FOXO1 is a master regulator of memory programming in CAR T cells

A major limitation of chimeric antigen receptor (CAR) T cell therapies is the poor persistence of these cells in vivo1. The expression of memory-associated genes in CAR T cells is linked to their long-term persistence in patients and clinical efficacy2-6, suggesting that memory programs may underpin durable CAR T cell function. Here we show that the transcription factor FOXO1 is responsible for promoting memory and restraining exhaustion in human CAR T cells. Pharmacological inhibition or gene editing of endogenous FOXO1 diminished the expression of memory-associated genes, promoted an exhaustion-like phenotype and impaired the antitumour activity of CAR T cells. Overexpression of FOXO1 induced a gene-expression program consistent with T cell memory and increased chromatin accessibility at FOXO1-binding motifs. CAR T cells that overexpressed FOXO1 retained their function, memory potential and metabolic fitness in settings of chronic stimulation, and exhibited enhanced persistence and tumour control in vivo. By contrast, overexpression of TCF1 (encoded by TCF7) did not enforce canonical memory programs or enhance the potency of CAR T cells. Notably, FOXO1 activity correlated with positive clinical outcomes of patients treated with CAR T cells or tumour-infiltrating lymphocytes, underscoring the clinical relevance of FOXO1 in cancer immunotherapy. Our results show that overexpressing FOXO1 can increase the antitumour activity of human CAR T cells, and highlight memory reprogramming as a broadly applicable approach for optimizing therapeutic T cell states.

Facts and Hopes in Using Omics to Advance Combined Immunotherapy Strategies

The field of oncology has been transformed by immune checkpoint inhibitors (ICI) and other immune-based agents; however, many patients do not receive a durable benefit. While biomarker assessments from pivotal ICI trials have uncovered certain mechanisms of resistance, results thus far have only scraped the surface. Mechanisms of resistance are as complex as the tumor microenvironment (TME) itself, and the development of effective therapeutic strategies will only be possible by building accurate models of the tumor-immune interface. With advancement of multi-omic technologies, high-resolution characterization of the TME is now possible. In addition to sequencing of bulk tumor, single-cell transcriptomic, proteomic, and epigenomic data as well as T-cell receptor profiling can now be simultaneously measured and compared between responders and nonresponders to ICI. Spatial sequencing and imaging platforms have further expanded the dimensionality of existing technologies. Rapid advancements in computation and data sharing strategies enable development of biologically interpretable machine learning models to integrate data from high-resolution, multi-omic platforms. These models catalyze the identification of resistance mechanisms and predictors of benefit in ICI-treated patients, providing scientific foundation for novel clinical trials. Moving forward, we propose a framework by which in silico screening, functional validation, and clinical trial biomarker assessment can be used for the advancement of combined immunotherapy strategies.

RAG suppresses group 2 innate lymphoid cells

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

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

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

Improved multi-ancestry fine-mapping identifies cis -regulatory variants underlying molecular traits and disease risk

Multi-ancestry statistical fine-mapping of cis -molecular quantitative trait loci ( cis -molQTL) aims to improve the precision of distinguishing causal cis -molQTLs from tagging variants. However, existing approaches fail to reflect shared genetic architectures. To solve this limitation, we present the Sum of Shared Single Effects (SuShiE) model, which leverages LD heterogeneity to improve fine-mapping precision, infer cross-ancestry effect size correlations, and estimate ancestry-specific expression prediction weights. We apply SuShiE to mRNA expression measured in PBMCs (n=956) and LCLs (n=814) together with plasma protein levels (n=854) from individuals of diverse ancestries in the TOPMed MESA and GENOA studies. We find SuShiE fine-maps cis -molQTLs for 16 % more genes compared with baselines while prioritizing fewer variants with greater functional enrichment. SuShiE infers highly consistent cis -molQTL architectures across ancestries on average; however, we also find evidence of heterogeneity at genes with predicted loss-of-function intolerance, suggesting that environmental interactions may partially explain differences in cis -molQTL effect sizes across ancestries. Lastly, we leverage estimated cis -molQTL effect-sizes to perform individual-level TWAS and PWAS on six white blood cell-related traits in AOU Biobank individuals (n=86k), and identify 44 more genes compared with baselines, further highlighting its benefits in identifying genes relevant for complex disease risk. Overall, SuShiE provides new insights into the cis -genetic architecture of molecular traits.

Pharmacological induction of chromatin remodeling drives chemosensitization in triple-negative breast cancer

Targeted therapies have improved outcomes for certain cancer subtypes, but cytotoxic chemotherapy remains a mainstay for triple-negative breast cancer (TNBC). The epithelial-to-mesenchymal transition (EMT) is a developmental program co-opted by cancer cells that promotes metastasis and chemoresistance. There are no therapeutic strategies specifically targeting mesenchymal-like cancer cells. We report that the US Food and Drug Administration (FDA)-approved chemotherapeutic eribulin induces ZEB1-SWI/SNF-directed chromatin remodeling to reverse EMT that curtails the metastatic propensity of TNBC preclinical models. Eribulin induces mesenchymal-to-epithelial transition (MET) in primary TNBC in patients, but conventional chemotherapy does not. In the treatment-naive setting, but not after acquired resistance to other agents, eribulin sensitizes TNBC cells to subsequent treatment with other chemotherapeutics. These findings provide an epigenetic mechanism of action of eribulin, supporting its use early in the disease process for MET induction to prevent metastatic progression and chemoresistance. These findings warrant prospective clinical evaluation of the chemosensitizing effects of eribulin in the treatment-naive setting.

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

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

Single-cell multiomics reveals the interplay of clonal evolution and cellular plasticity in hepatoblastoma

Hepatoblastomas (HB) display heterogeneous cellular phenotypes that influence the clinical outcome, but the underlying mechanisms are poorly understood. Here, we use a single-cell multiomic strategy to unravel the molecular determinants of this plasticity. We identify a continuum of HB cell states between hepatocytic (scH), liver progenitor (scLP) and mesenchymal (scM) differentiation poles, with an intermediate scH/LP population bordering scLP and scH areas in spatial transcriptomics. Chromatin accessibility landscapes reveal the gene regulatory networks of each differentiation pole, and the sequence of transcription factor activations underlying cell state transitions. Single-cell mapping of somatic alterations reveals the clonal architecture of each tumor, showing that each genetic subclone displays its own range of cellular plasticity across differentiation states. The most scLP subclones, overexpressing stem cell and DNA repair genes, proliferate faster after neo-adjuvant chemotherapy. These results highlight how the interplay of clonal evolution and epigenetic plasticity shapes the potential of HB subclones to respond to chemotherapy.

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

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

Deciphering cell types by integrating scATAC-seq data with genome sequences

The single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) technology provides insight into gene regulation and epigenetic heterogeneity at single-cell resolution, but cell annotation from scATAC-seq remains challenging due to high dimensionality and extreme sparsity within the data. Existing cell annotation methods mostly focus on the cell peak matrix without fully utilizing the underlying genomic sequence. Here we propose a method, SANGO, for accurate single-cell annotation by integrating genome sequences around the accessibility peaks within scATAC data. The genome sequences of peaks are encoded into low-dimensional embeddings, and then iteratively used to reconstruct the peak statistics of cells through a fully connected network. The learned weights are considered as regulatory modes to represent cells, and utilized to align the query cells and the annotated cells in the reference data through a graph transformer network for cell annotations. SANGO was demonstrated to consistently outperform competing methods on 55 paired scATAC-seq datasets across samples, platforms and tissues. SANGO was also shown to be able to detect unknown tumor cells through attention edge weights learned by the graph transformer. Moreover, from the annotated cells, we found cell-type-specific peaks that provide functional insights/biological signals through expression enrichment analysis, cis-regulatory chromatin interaction analysis and motif enrichment analysis.

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

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

Fast and flexible profiling of chromatin accessibility and total RNA expression in single nuclei using Microwell-seq3

Single cell chromatin accessibility profiling and transcriptome sequencing are the most widely used technologies for single-cell genomics. Here, we present Microwell-seq3, a high-throughput and facile platform for high-sensitivity single-nucleus chromatin accessibility or full-length transcriptome profiling. The method combines a preindexing strategy and a penetrable chip-in-a-tube for single nucleus loading and DNA amplification and therefore does not require specialized equipment. We used Microwell-seq3 to profile chromatin accessibility in more than 200,000 single nuclei and the full-length transcriptome in ~50,000 nuclei from multiple adult mouse tissues. Compared with the existing polyadenylated transcript capture methods, integrative analysis of cell type-specific regulatory elements and total RNA expression uncovered comprehensive cell type heterogeneity in the brain. Gene regulatory networks based on chromatin accessibility profiling provided an improved cell type communication model. Finally, we demonstrated that Microwell-seq3 can identify malignant cells and their specific regulons in spontaneous lung tumors of aged mice. We envision a broad application of Microwell-seq3 in many areas of research.

Construction of single-cell cross-species chromatin accessibility landscapes with combinatorial-hybridization-based ATAC-seq

Despite recent advances in single-cell genomics, the lack of maps for single-cell candidate cis-regulatory elements (cCREs) in non-mammal species has limited our exploration of conserved regulatory programs across vertebrates and invertebrates. Here, we developed a combinatorial-hybridization-based method for single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) named CH-ATAC-seq, enabling the construction of single-cell accessible chromatin landscapes for zebrafish, Drosophila, and earthworms (Eisenia andrei). By integrating scATAC censuses of humans, monkeys, and mice, we systematically identified 152 distinct main cell types and around 0.8 million cell-type-specific cCREs. Our analysis provided insights into the conservation of neural, muscle, and immune lineages across species, while epithelial cells exhibited a higher organ-origin heterogeneity. Additionally, a large-scale gene regulatory network (GRN) was constructed in four vertebrates by integrating scRNA-seq censuses. Overall, our study provides a valuable resource for comparative epigenomics, identifying the evolutionary conservation and divergence of gene regulation across different species.

PACS allows comprehensive dissection of multiple factors governing chromatin accessibility from snATAC-seq data

Single nucleus ATAC-seq (snATAC-seq) experimental designs have become increasingly complex with multiple factors that might affect chromatin accessibility, including genotype, cell type, tissue of origin, sample location, batch, etc., whose compound effects are difficult to test by existing methods. In addition, current snATAC-seq data present statistical difficulties due to their sparsity and variations in individual sequence capture. To address these problems, we present a zero-adjusted statistical model, Probability model of Accessible Chromatin of Single cells (PACS), that can allow complex hypothesis testing of factors that affect accessibility while accounting for sparse and incomplete data. For differential accessibility analysis, PACS controls the false positive rate and achieves on average a 17% to 122% higher power than existing tools. We demonstrate the effectiveness of PACS through several analysis tasks including supervised cell type annotation, compound hypothesis testing, batch effect correction, and spatiotemporal modeling. We apply PACS to several datasets from a variety of tissues and show its ability to reveal previously undiscovered insights in snATAC-seq data.

txci-ATAC-seq: a massive-scale single-cell technique to profile chromatin accessibility

We develop a large-scale single-cell ATAC-seq method by combining Tn5-based pre-indexing with 10× Genomics barcoding, enabling the indexing of up to 200,000 nuclei across multiple samples in a single reaction. We profile 449,953 nuclei across diverse tissues, including the human cortex, mouse brain, human lung, mouse lung, mouse liver, and lung tissue from a club cell secretory protein knockout (CC16-/-) model. Our study of CC16-/- nuclei uncovers previously underappreciated technical artifacts derived from remnant 129 mouse strain genetic material, which cause profound cell-type-specific changes in regulatory elements near many genes, thereby confounding the interpretation of this commonly referenced mouse model.

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

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

The transcription factor NF-κB orchestrates nucleosome remodeling during the primary response to Toll-like receptor 4 signaling

Inducible nucleosome remodeling at hundreds of latent enhancers and several promoters shapes the transcriptional response to Toll-like receptor 4 (TLR4) signaling in macrophages. We aimed to define the identities of the transcription factors that promote TLR-induced remodeling. An analysis strategy based on ATAC-seq and single-cell ATAC-seq that enriched for genomic regions most likely to undergo remodeling revealed that the transcription factor nuclear factor κB (NF-κB) bound to all high-confidence peaks marking remodeling during the primary response to the TLR4 ligand, lipid A. Deletion of NF-κB subunits RelA and c-Rel resulted in the loss of remodeling at high-confidence ATAC-seq peaks, and CRISPR-Cas9 mutagenesis of NF-κB-binding motifs impaired remodeling. Remodeling selectivity at defined regions was conferred by collaboration with other inducible factors, including IRF3- and MAP-kinase-induced factors. Thus, NF-κB is unique among TLR4-activated transcription factors in its broad contribution to inducible nucleosome remodeling, alongside its ability to activate poised enhancers and promoters assembled into open chromatin.

scPerturb: harmonized single-cell perturbation data

Analysis across a growing number of single-cell perturbation datasets is hampered by poor data interoperability. To facilitate development and benchmarking of computational methods, we collect a set of 44 publicly available single-cell perturbation-response datasets with molecular readouts, including transcriptomics, proteomics and epigenomics. We apply uniform quality control pipelines and harmonize feature annotations. The resulting information resource, scPerturb, enables development and testing of computational methods, and facilitates comparison and integration across datasets. We describe energy statistics (E-statistics) for quantification of perturbation effects and significance testing, and demonstrate E-distance as a general distance measure between sets of single-cell expression profiles. We illustrate the application of E-statistics for quantifying similarity and efficacy of perturbations. The perturbation-response datasets and E-statistics computation software are publicly available at scperturb.org. This work provides an information resource for researchers working with single-cell perturbation data and recommendations for experimental design, including optimal cell counts and read depth.

High-resolution omics of vascular ageing and inflammatory pathways in neurodegeneration

High-resolution omics, particularly single-cell and spatial transcriptomic profiling, are rapidly enhancing our comprehension of the normal molecular diversity of gliovascular cells, as well as their age-related changes that contribute to neurodegeneration. With more omic profiling studies being conducted, it is becoming increasingly essential to synthesise valuable information from the rapidly accumulating findings. In this review, we present an overview of the molecular features of neurovascular and glial cells that have been recently discovered through omic profiling, with a focus on those that have potentially significant functional implications and/or show cross-species differences between human and mouse, and that are linked to vascular deficits and inflammatory pathways in ageing and neurodegenerative disorders. Additionally, we highlight the translational applications of omic profiling, and discuss omic-based strategies to accelerate biomarker discovery and facilitate disease course-modifying therapeutics development for neurodegenerative conditions.

Applications of single-cell omics for chimeric antigen receptor T cell therapy

Chimeric antigen receptor (CAR) T cell therapy is a promising cancer treatment modality. The breakthroughs in CAR T cell therapy were, in part, possible with the help of cell analysis methods, such as single-cell analysis. Bulk analyses have provided invaluable information regarding the complex molecular dynamics of CAR T cells, but their results are an average of thousands of signals in CAR T or tumour cells. Since cancer is a heterogeneous disease where each minute detail of a subclone could change the outcome of the treatment, single-cell analysis could prove to be a powerful instrument in deciphering the secrets of tumour microenvironment for cancer immunotherapy. With the recent studies in all aspects of adoptive cell therapy making use of single-cell analysis, a comprehensive review of the recent preclinical and clinical findings in CAR T cell therapy was needed. Here, we categorized and summarized the key points of the studies in which single-cell analysis provided insights into the genomics, epigenomics, transcriptomics and proteomics as well as their respective multi-omics of CAR T cell therapy.

RNA structure profiling at single-cell resolution reveals new determinants of cell identity

RNA structure is critical for multiple steps in gene regulation. However, how the structures of transcripts differ both within and between individual cells is unknown. Here we develop a SHAPE-inspired method called single-cell structure probing of RNA transcripts that enables simultaneous determination of transcript secondary structure and abundance at single-cell resolution. We apply single-cell structure probing of RNA transcripts to human embryonic stem cells and differentiating neurons. Remarkably, RNA structure is more homogeneous in human embryonic stem cells compared with neurons, with the greatest homogeneity found in coding regions. More extensive heterogeneity is found within 3' untranslated regions and is determined by specific RNA-binding proteins. Overall RNA structure profiles better discriminate cell type identity and differentiation stage than gene expression profiles alone. We further discover a cell-type variable region of 18S ribosomal RNA that is associated with cell cycle and translation control. Our method opens the door to the systematic characterization of RNA structure-function relationships at single-cell resolution.

Harnessing HetHydrogel: A Universal Platform to Dropletize Single-Cell Multiomics

A universal platform is developed for dropletizing single cell plate-based multiomic assays, consisting of three main pillars: a miniaturized open Heterogeneous Hydrogel reactor (abbreviated HetHydrogel) for multi-step biochemistry, its tunable permeability that allows Tn5 tagmentation, and single cell droplet barcoding. Through optimizing the HetHydrogel manufacturing procedure, the chemical composition, and cell permeation conditions, simultaneous high-throughput mitochondrial DNA genotyping and chromatin profiling at the single-cell level are demonstrated using a mixed-species experiment. This platform offers a powerful way to investigate the genotype-phenotype relationships of various mtDNA mutations in biological processes. The HetHydrogel platform is believed to have the potential to democratize droplet technologies, upgrading a whole range of plate-based single cell assays to high throughput format.

Accelerating Single-Cell Sequencing Data Analysis with SciDAP: A User-Friendly Approach

Single-cell (sc) RNA, ATAC and Multiome sequencing became powerful tools for uncovering biological and disease mechanisms. Unfortunately, manual analysis of sc data presents multiple challenges due to large data volumes and complexity of configuration parameters. This complexity, as well as not being able to reproduce a computational environment, affects the reproducibility of analysis results. The Scientific Data Analysis Platform (https://SciDAP.com) allows biologists without computational expertise to analyze sequencing-based data using portable and reproducible pipelines written in Common Workflow Language (CWL). Our suite of computational pipelines addresses the most common needs in scRNA-Seq, scATAC-Seq and scMultiome data analysis. When executed on SciDAP, it offers a user-friendly alternative to manual data processing, eliminating the need for coding expertise. In this protocol, we describe the use of SciDAP to analyze scMultiome data. Similar approaches can be used for analysis of scRNA-Seq, scATAC-Seq and scVDJ-Seq datasets.

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

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

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

Cis -regulatory elements (CREs) are critical in regulating gene expression, and yet our understanding of CRE evolution remains a challenge. Here, we constructed a comprehensive single-cell atlas of chromatin accessibility in Oryza sativa , integrating data from 104,029 nuclei representing 128 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 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 critical silencer CREs maintained by evolution. Collectively, this comparative genomics approach highlights the dynamics of cell-type-specific CRE evolution in plants.

Accessible chromatin maps of inflammatory bowel disease intestine nominate cell-type mediators of genetic disease risk

Inflammatory Bowel Disease ( IBD ) is a chronic and often debilitating autoinflammatory condition, with an increasing incidence in children. Standard-of-care therapies lead to sustained transmural healing and clinical remission in fewer than one-third of patients. For children, TNFα inhibition remains the only FDA-approved biologic therapy, providing an even greater urgency to understanding mechanisms of response. Genome-wide association studies ( GWAS ) have identified 418 independent genetic risk loci contributing to IBD, yet the majority are noncoding and their mechanisms of action are difficult to decipher. If causal, they likely alter transcription factor ( TF ) binding and downstream gene expression in particular cell types and contexts. To bridge this knowledge gap, we built a novel resource: multiome-seq (tandem single-nuclei ( sn )RNA-seq and chromatin accessibility ( snATAC )-seq) of intestinal tissue from pediatric IBD patients, where anti-TNF response was defined by endoscopic healing. From the snATAC-seq data, we generated a first-time atlas of chromatin accessibility (putative regulatory elements) for diverse intestinal cell types in the context of IBD. For cell types/contexts mediating genetic risk, we reasoned that accessible chromatin will co-localize with genetic disease risk loci. We systematically tested for significant co-localization of our chromatin accessibility maps and risk variants for 758 GWAS traits. Globally, genetic risk variants for IBD, autoimmune and inflammatory diseases are enriched in accessible chromatin of immune populations, while other traits (e.g., colorectal cancer, metabolic) are enriched in epithelial and stromal populations. This resource opens new avenues to uncover the complex molecular and cellular mechanisms mediating genetic disease risk.

Cell-type-directed design of synthetic enhancers

Transcriptional enhancers act as docking stations for combinations of transcription factors and thereby regulate spatiotemporal activation of their target genes1. It has been a long-standing goal in the field to decode the regulatory logic of an enhancer and to understand the details of how spatiotemporal gene expression is encoded in an enhancer sequence. Here we show that deep learning models2-6, can be used to efficiently design synthetic, cell-type-specific enhancers, starting from random sequences, and that this optimization process allows detailed tracing of enhancer features at single-nucleotide resolution. We evaluate the function of fully synthetic enhancers to specifically target Kenyon cells or glial cells in the fruit fly brain using transgenic animals. We further exploit enhancer design to create 'dual-code' enhancers that target two cell types and minimal enhancers smaller than 50 base pairs that are fully functional. By examining the state space searches towards local optima, we characterize enhancer codes through the strength, combination and arrangement of transcription factor activator and transcription factor repressor motifs. Finally, we apply the same strategies to successfully design human enhancers, which adhere to enhancer rules similar to those of Drosophila enhancers. Enhancer design guided by deep learning leads to better understanding of how enhancers work and shows that their code can be exploited to manipulate cell states.

Application of Single-Cell Assay for Transposase-Accessible Chromatin with High Throughput Sequencing in Plant Science: Advances, Technical Challenges, and Prospects

The Single-cell Assay for Transposase-Accessible Chromatin with high throughput sequencing (scATAC-seq) has gained increasing popularity in recent years, allowing for chromatin accessibility to be deciphered and gene regulatory networks (GRNs) to be inferred at single-cell resolution. This cutting-edge technology now enables the genome-wide profiling of chromatin accessibility at the cellular level and the capturing of cell-type-specific cis-regulatory elements (CREs) that are masked by cellular heterogeneity in bulk assays. Additionally, it can also facilitate the identification of rare and new cell types based on differences in chromatin accessibility and the charting of cellular developmental trajectories within lineage-related cell clusters. Due to technical challenges and limitations, the data generated from scATAC-seq exhibit unique features, often characterized by high sparsity and noise, even within the same cell type. To address these challenges, various bioinformatic tools have been developed. Furthermore, the application of scATAC-seq in plant science is still in its infancy, with most research focusing on root tissues and model plant species. In this review, we provide an overview of recent progress in scATAC-seq and its application across various fields. We first conduct scATAC-seq in plant science. Next, we highlight the current challenges of scATAC-seq in plant science and major strategies for cell type annotation. Finally, we outline several future directions to exploit scATAC-seq technologies to address critical challenges in plant science, ranging from plant ENCODE(The Encyclopedia of DNA Elements) project construction to GRN inference, to deepen our understanding of the roles of CREs in plant biology.

scaDA: A Novel Statistical Method for Differential Analysis of Single-Cell Chromatin Accessibility Sequencing Data

Single-cell ATAC-seq sequencing data (scATAC-seq) has been widely used to investigate chromatin accessibility on the single-cell level. One important application of scATAC-seq data analysis is differential chromatin accessibility analysis. However, the data characteristics of scATAC-seq such as excessive zeros and large variability of chromatin accessibility across cells impose a unique challenge for DA analysis. Existing statistical methods focus on detecting the mean difference of the chromatin accessible regions while overlooking the distribution difference. Motivated by real data exploration that distribution difference exists among cell types, we introduce a novel composite statistical test named "scaDA", which is based on zero-inflated negative binomial model (ZINB), for performing differential distribution analysis of chromatin accessibility by jointly testing the abundance, prevalence and dispersion simultaneously. Benefiting from both dispersion shrinkage and iterative refinement of mean and prevalence parameter estimates, scaDA demonstrates its superiority to both ZINB-based likelihood ratio tests and published methods by achieving the highest power and best FDR control in a comprehensive simulation study. In addition to demonstrating the highest power in three real sc-multiome data analyses, scaDA successfully identifies differentially accessible regions in microglia from sc-multiome data for an Alzheimer's disease (AD) study, regions which are most enriched in GO terms related to neurogenesis, the clinical phenotype of AD, and SNPs identified in AD-associated GWAS.

scENCORE: leveraging single-cell epigenetic data to predict chromatin conformation using graph embedding

Dynamic compartmentalization of eukaryotic DNA into active and repressed states enables diverse transcriptional programs to arise from a single genetic blueprint, whereas its dysregulation can be strongly linked to a broad spectrum of diseases. While single-cell Hi-C experiments allow for chromosome conformation profiling across many cells, they are still expensive and not widely available for most labs. Here, we propose an alternate approach, scENCORE, to computationally reconstruct chromatin compartments from the more affordable and widely accessible single-cell epigenetic data. First, scENCORE constructs a long-range epigenetic correlation graph to mimic chromatin interaction frequencies, where nodes and edges represent genome bins and their correlations. Then, it learns the node embeddings to cluster genome regions into A/B compartments and aligns different graphs to quantify chromatin conformation changes across conditions. Benchmarking using cell-type-matched Hi-C experiments demonstrates that scENCORE can robustly reconstruct A/B compartments in a cell-type-specific manner. Furthermore, our chromatin confirmation switching studies highlight substantial compartment-switching events that may introduce substantial regulatory and transcriptional changes in psychiatric disease. In summary, scENCORE allows accurate and cost-effective A/B compartment reconstruction to delineate higher-order chromatin structure heterogeneity in complex tissues.

scDecouple: decoupling cellular response from infected proportion bias in scCRISPR-seq

Single-cell clustered regularly interspaced short palindromic repeats-sequencing (scCRISPR-seq) is an emerging high-throughput CRISPR screening technology where the true cellular response to perturbation is coupled with infected proportion bias of guide RNAs (gRNAs) across different cell clusters. The mixing of these effects introduces noise into scCRISPR-seq data analysis and thus obstacles to relevant studies. We developed scDecouple to decouple true cellular response of perturbation from the influence of infected proportion bias. scDecouple first models the distribution of gene expression profiles in perturbed cells and then iteratively finds the maximum likelihood of cell cluster proportions as well as the cellular response for each gRNA. We demonstrated its performance in a series of simulation experiments. By applying scDecouple to real scCRISPR-seq data, we found that scDecouple enhances the identification of biologically perturbation-related genes. scDecouple can benefit scCRISPR-seq data analysis, especially in the case of heterogeneous samples or complex gRNA libraries.

Single-cell sequencing analysis within biologically relevant dimensions

The currently predominant approach to transcriptomic and epigenomic single-cell analysis depends on a rigid perspective constrained by reduced dimensions and algorithmically derived and annotated clusters. Here, we developed Seqtometry (sequencing-to-measurement), a single-cell analytical strategy based on biologically relevant dimensions enabled by advanced scoring with multiple gene sets (signatures) for examination of gene expression and accessibility across various organ systems. By utilizing information only in the form of specific signatures, Seqtometry bypasses unsupervised clustering and individual annotations of clusters. Instead, Seqtometry combines qualitative and quantitative cell-type identification with specific characterization of diverse biological processes under experimental or disease conditions. Comprehensive analysis by Seqtometry of various immune cells as well as other cells from different organs and disease-induced states, including multiple myeloma and Alzheimer's disease, surpasses corresponding cluster-based analytical output. We propose Seqtometry as a single-cell sequencing analysis approach applicable for both basic and clinical research.

scATAnno: Automated Cell Type Annotation for single-cell ATAC Sequencing Data

The recent advances in single-cell epigenomic techniques have created a growing demand for scATAC-seq analysis. One key task is to determine cell types based on epigenetic profiling. We introduce scATAnno, a workflow designed to automatically annotate scATAC-seq data using large-scale scATAC-seq reference atlases. This workflow can generate scATAC-seq reference atlases from publicly available datasets, and enable accurate cell type annotation by integrating query data with reference atlases, without the aid of scRNA-seq profiling. To enhance annotation accuracy, we have incorporated KNN-based and weighted distance-based uncertainty scores to effectively detect unknown cell populations within the query data. We showcase the utility of scATAnno across multiple datasets, including peripheral blood mononuclear cell (PBMC), basal cell carcinoma (BCC) and Triple Negative Breast Cancer (TNBC), and demonstrate that scATAnno accurately annotates cell types across conditions. Overall, scATAnno is a powerful tool for cell type annotation in scATAC-seq data and can aid in the interpretation of new scATAC-seq datasets in complex biological systems.

Single-cell joint profiling of multiple epigenetic proteins and gene transcription

Sculpting the epigenome with a combination of histone modifications and transcription factor occupancy determines gene transcription and cell fate specification. Here, we first develop uCoTarget, utilizing a split-pool barcoding strategy for realizing ultrahigh-throughput single-cell joint profiling of multiple epigenetic proteins. Through extensive optimization for sensitivity and multimodality resolution, we demonstrate that uCoTarget enables simultaneous detection of five histone modifications (H3K27ac, H3K4me3, H3K4me1, H3K36me3, and H3K27me3) in 19,860 single cells. We applied uCoTarget to the in vitro generation of hematopoietic stem/progenitor cells (HSPCs) from human embryonic stem cells, presenting multimodal epigenomic profiles in 26,418 single cells. uCoTarget reveals establishment of pairing of HSPC enhancers (H3K27ac) and promoters (H3K4me3) and RUNX1 engagement priming for H3K27ac activation along the HSPC path. We then develop uCoTargetX, an expansion of uCoTarget to simultaneously measure transcriptome and multiple epigenome targets. Together, our methods enable generalizable, versatile multimodal profiles for reconstructing comprehensive epigenome and transcriptome landscapes and analyzing the regulatory interplay at single-cell level.

scATAC-Ref: a reference of scATAC-seq with known cell labels in multiple species

Chromatin accessibility profiles at single cell resolution can reveal cell type-specific regulatory programs, help dissect highly specialized cell functions and trace cell origin and evolution. Accurate cell type assignment is critical for effectively gaining biological and pathological insights, but is difficult in scATAC-seq. Hence, by extensively reviewing the literature, we designed scATAC-Ref (https://bio.liclab.net/scATAC-Ref/), a manually curated scATAC-seq database aimed at providing a comprehensive, high-quality source of chromatin accessibility profiles with known cell labels across broad cell types. Currently, scATAC-Ref comprises 1 694 372 cells with known cell labels, across various biological conditions, >400 cell/tissue types and five species. We used uniform system environment and software parameters to perform comprehensive downstream analysis on these chromatin accessibility profiles with known labels, including gene activity score, TF enrichment score, differential chromatin accessibility regions, pathway/GO term enrichment analysis and co-accessibility interactions. The scATAC-Ref also provided a user-friendly interface to query, browse and visualize cell types of interest, thereby providing a valuable resource for exploring epigenetic regulation in different tissues and cell types.

Histone lactylation couples cellular metabolism with developmental gene regulatory networks

Embryonic cells exhibit diverse metabolic states. Recent studies have demonstrated that metabolic reprogramming drives changes in cell identity by affecting gene expression. However, the connection between cellular metabolism and gene expression remains poorly understood. Here we report that glycolysis-regulated histone lactylation couples the metabolic state of embryonic cells with chromatin organization and gene regulatory network (GRN) activation. We found that lactylation marks genomic regions of glycolytic embryonic tissues, like the neural crest (NC) and pre-somitic mesoderm. Histone lactylation occurs in the loci of NC genes as these cells upregulate glycolysis. This process promotes the accessibility of active enhancers and the deployment of the NC GRN. Reducing the deposition of the mark by targeting LDHA/B leads to the downregulation of NC genes and the impairment of cell migration. The deposition of lactyl-CoA on histones at NC enhancers is supported by a mechanism that involves transcription factors SOX9 and YAP/TEAD. These findings define an epigenetic mechanism that integrates cellular metabolism with the GRNs that orchestrate embryonic development.

Repression of latent NF-κB enhancers by PDX1 regulates β cell functional heterogeneity

Interactions between lineage-determining and activity-dependent transcription factors determine single-cell identity and function within multicellular tissues through incompletely known mechanisms. By assembling a single-cell atlas of chromatin state within human islets, we identified β cell subtypes governed by either high or low activity of the lineage-determining factor pancreatic duodenal homeobox-1 (PDX1). β cells with reduced PDX1 activity displayed increased chromatin accessibility at latent nuclear factor κB (NF-κB) enhancers. Pdx1 hypomorphic mice exhibited de-repression of NF-κB and impaired glucose tolerance at night. Three-dimensional analyses in tandem with chromatin immunoprecipitation (ChIP) sequencing revealed that PDX1 silences NF-κB at circadian and inflammatory enhancers through long-range chromatin contacts involving SIN3A. Conversely, Bmal1 ablation in β cells disrupted genome-wide PDX1 and NF-κB DNA binding. Finally, antagonizing the interleukin (IL)-1β receptor, an NF-κB target, improved insulin secretion in Pdx1 hypomorphic islets. Our studies reveal functional subtypes of single β cells defined by a gradient in PDX1 activity and identify NF-κB as a target for insulinotropic therapy.

Single-Nucleus ATAC-seq for Mapping Chromatin Accessibility in Individual Cells of Murine Hearts

During the last decade a wide range of single-cell and single-nucleus next-generation sequencing techniques have been developed, which revolutionized detection of rare cell populations, enabling creation of comprehensive cell atlases of complex organs and tissues. State-of-the-art methods do not only allow classical transcriptomics of individual cells but also comprise a number of epigenetic approaches, including assessment of chromatin accessibility by single-nucleus Assay for Transposase Accessible Chromatin ATAC-seq (snATAC-seq). The snATAC-seq assay detects "open chromatin," a term for low nucleosome occupancy of genomic regions, which is a prerequisite for effective transcription factor binding. Information about open chromatin at the single-nucleus level helps to recognize epigenetic changes, sometimes before transcription of respective genes occurs. snATAC-seq detects cellular heterogeneity in otherwise still transcriptionally and/or morphologically homogeneous cell populations. Chromatin accessibility assays may be used to detect epigenetic changes in cardiac lineages during heart development, chromatin landscape changes during aging, and epigenetic alterations in heart diseases. Here, we provide an optimized protocol for snATAC-seq of murine hearts. We describe isolation of single nuclei from snap-frozen hearts, provide hints for preparation of libraries suitable for snATAC-seq next-generation sequencing (NGS) using the Chromium 10× platform, and give general recommendations for downstream analysis using conventional bioinformatic pipelines and packages. The protocol should serve as a beginner's guide to generate high-quality snATAC-seq datasets and to perform chromatin accessibility analysis of individual heart-derived cell nuclei.

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

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

Modeling fragment counts improves single-cell ATAC-seq analysis

Single-cell ATAC sequencing coverage in regulatory regions is typically binarized as an indicator of open chromatin. Here we show that binarization is an unnecessary step that neither improves goodness of fit, clustering, cell type identification nor batch integration. Fragment counts, but not read counts, should instead be modeled, which preserves quantitative regulatory information. These results have immediate implications for single-cell ATAC sequencing analysis.

Dynamic chromatin remodeling in cycling human endometrium at single-cell level

Estrogen-dependent proliferation followed by progesterone-dependent differentiation of the endometrium culminates in a short implantation window. We performed single-cell assay for transposase-accessible chromatin with sequencing on endometrial samples obtained across the menstrual cycle to investigate the regulation of temporal gene networks that control embryo implantation. We identify uniquely accessible chromatin regions in all major cellular constituents of the endometrium, delineate temporal patterns of coordinated chromatin remodeling in epithelial and stromal cells, and gain mechanistic insights into the emergence of a receptive state through integrated analysis of enriched transcription factor (TF) binding sites in dynamic chromatin regions, chromatin immunoprecipitation sequencing analyses, and gene expression data. We demonstrate that the implantation window coincides with pervasive cooption of transposable elements (TEs) into the regulatory chromatin landscape of decidualizing cells and expression of TE-derived transcripts in a spatially defined manner. Our data constitute a comprehensive map of the chromatin changes that control TF activities in a cycling endometrium at cellular resolution.

Single-cell multi-omics analysis of human testicular germ cell tumor reveals its molecular features and microenvironment

Seminoma is the most common malignant solid tumor in 14 to 44 year-old men. However, its molecular features and tumor microenvironment (TME) is largely unexplored. Here, we perform a series of studies via genomics profiling (single cell multi-omics and spatial transcriptomics) and functional examination using seminoma samples and a seminoma cell line. We identify key gene expression programs share between seminoma and primordial germ cells, and further characterize the functions of TFAP2C in promoting tumor invasion and migration. We also identify 15 immune cell subtypes in TME, and find that subtypes with exhaustion features were located closer to the tumor region through combined spatial transcriptome analysis. Furthermore, we identify key pathways and genes that may facilitate seminoma disseminating beyond the seminiferous tubules. These findings advance our knowledge of seminoma tumorigenesis and produce a multi-omics atlas of in situ human seminoma microenvironment, which could help discover potential therapy targets for seminoma.

Characterization of human pluripotent stem cell differentiation by single-cell dual-omics analyses

In vivo differentiation of human pluripotent stem cells (hPSCs) has unique advantages, such as multilineage differentiation, angiogenesis, and close cell-cell interactions. To systematically investigate multilineage differentiation mechanisms of hPSCs, we constructed the in vivo hPSC differentiation landscape containing 239,670 cells using teratoma models. We identified 43 cell types, inferred 18 cell differentiation trajectories, and characterized common and specific gene regulation patterns during hPSC differentiation at both transcriptional and epigenetic levels. Additionally, we developed the developmental single-cell Basic Local Alignment Search Tool (dscBLAST), an R-based cell identification tool, to simplify the identification processes of developmental cells. Using dscBLAST, we aligned cells in multiple differentiation models to normally developing cells to further understand their differentiation states. Overall, our study offers new insights into stem cell differentiation and human embryonic development; dscBLAST shows favorable cell identification performance, providing a powerful identification tool for developmental cells.