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

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.

Decoding changes in tumor-infiltrating leukocytes through dynamic experimental models and single-cell technologies

The ability to characterize immune cells and explore the molecular interactions that govern their functions has never been greater, fueled in recent years by the revolutionary advance of single-cell analysis platforms. However, precisely how immune cells respond to different stimuli and where differentiation processes and effector functions operate remain incompletely understood. Inferring cellular fate within single-cell transcriptomic analyses is now omnipresent, despite the assumptions typically required in such analyses. Recently developed experimental models support dynamic analyses of the immune response, providing insights into the temporal changes that occur within cells and the tissues in which such transitions occur. Here we will review these approaches and discuss how these can be combined with single-cell technologies to develop a deeper understanding of the immune responses that should support the development of better therapeutic options for patients.

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.

Ploidy inference from single-cell data: application to human and mouse cell atlases

Ploidy is relevant to numerous biological phenomena, including development, metabolism, and tissue regeneration. Single-cell RNA-seq and other omics studies are revolutionizing our understanding of biology, yet they have largely overlooked ploidy. This is likely due to the additional assay step required for ploidy measurement. Here, we developed a statistical method to infer ploidy from single-cell ATAC-seq data, addressing this gap. When applied to data from human and mouse cell atlases, our method enabled systematic detection of polyploidy across diverse cell types. This method allows for the integration of ploidy analysis into single-cell studies. Additionally, this method can be adapted to detect the proliferating stage in the cell cycle and copy number variations in cancer cells. The software is implemented as the scPloidy package of the R software and is freely available from CRAN.

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

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

SIGNIFICANCE

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

Epigenetic regulation of innate lymphoid cells

Innate lymphoid cells (ILCs) lack antigen-specific receptors and are considered the innate arm of the immune system, phenotypically and functionally mirroring CD4+ helper T cells. ILCs are categorized into groups 1, 2, and 3 based on transcription factors and cytokine expression. ILCs predominantly reside in mucosal tissues and play important roles in regional immune responses. The development and function of ILC subsets are controlled by both transcriptional and epigenetic mechanisms, which have been extensively studied in recent years. Epigenetic regulation refers to inheritable changes in gene expression that occur without affecting DNA sequences. This mainly includes chromatin status, histone modifications, and DNA methylation. In this review, we summarize recent discoveries on epigenetic mechanisms regulating ILC development and function, and how these regulations affect disease progression under pathological conditions. Although the ablation of specific epigenetic regulators can cause global changes in corresponding epigenetic modifications to the chromatin, only partial genes with altered epigenetic modifications change their mRNA expression, resulting in specific outcomes in cell differentiation and function. Therefore, elucidating epigenetic mechanisms underlying the regulation of ILCs will provide potential targets for the diagnosis and treatment of inflammatory diseases.

Deep generative models in single-cell omics

Deep Generative Models (DGMs) are becoming instrumental for inferring probability distributions inherent to complex processes, such as most questions in biomedical research. For many years, there was a lack of mathematical methods that would allow this inference in the scarce data scenario of biomedical research. The advent of single-cell omics has finally made square the so-called "skinny matrix", allowing to apply mathematical methods already extensively used in other areas. Moreover, it is now possible to integrate data at different molecular levels in thousands or even millions of samples, thanks to the number of single-cell atlases being collaboratively generated. Additionally, DGMs have proven useful in other frequent tasks in single-cell analysis pipelines, from dimensionality reduction, cell type annotation to RNA velocity inference. In spite of its promise, DGMs need to be used with caution in biomedical research, paying special attention to its use to answer the right questions and the definition of appropriate error metrics and validation check points that confirm not only its correct use but also its relevance. All in all, DGMs provide an exciting tool that opens a bright future for the integrative analysis of single-cell -omics to understand health and disease.

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.

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

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

Building and analyzing metacells in single-cell genomics data

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

Conserved Cis-Acting Range Extender Element Mediates Extreme Long-Range Enhancer Activity in Mammals

While most mammalian enhancers regulate their cognate promoters over moderate distances of tens of kilobases (kb), some enhancers act over distances in the megabase range. The sequence features enabling such extreme-distance enhancer-promoter interactions remain elusive. Here, we used in vivo enhancer replacement experiments in mice to show that short- and medium-range enhancers cannot initiate gene expression at extreme-distance range. We uncover a novel conserved cis-acting element, Range EXtender (REX), that confers extreme-distance regulatory activity and is located next to a long-range enhancer of Sall1. The REX element itself has no endogenous enhancer activity. However, addition of the REX to other short- and mid-range enhancers substantially increases their genomic interaction range. In the most extreme example observed, addition of the REX increased the range of an enhancer by an order of magnitude, from its native 71kb to 840kb. The REX element contains highly conserved [C/T]AATTA homeodomain motifs. These motifs are enriched around long-range limb enhancers genome-wide, including the ZRS, a benchmark long-range limb enhancer of Shh. Mutating the [C/T]AATTA motifs within the ZRS does not affect its limb-specific enhancer activity at short range, but selectively abolishes its long-range activity, resulting in severe limb reduction in knock-in mice. In summary, we identify a sequence signature globally associated with long-range enhancer-promoter interactions and describe a prototypical REX element that is necessary and sufficient to confer extreme-distance gene activation by remote enhancers.

The Extra-Islet Pancreas Supports Autoimmunity in Human Type 1 Diabetes

In autoimmune Type 1 diabetes (T1D), immune cells infiltrate and destroy the islets of Langerhans - islands of endocrine tissue dispersed throughout the pancreas. However, the contribution of cellular programs outside islets to insulitis is unclear. Here, using CO-Detection by indEXing (CODEX) tissue imaging and cadaveric pancreas samples, we simultaneously examine islet and extra-islet inflammation in human T1D. We identify four sub-states of inflamed islets characterized by the activation profiles of CD8 + T cells enriched in islets relative to the surrounding tissue. We further find that the extra-islet space of lobules with extensive islet-infiltration differs from the extra-islet space of less infiltrated areas within the same tissue section. Finally, we identify lymphoid structures away from islets enriched in CD45RA + T cells - a population also enriched in one of the inflamed islet sub-states. Together, these data help define the coordination between islets and the extra-islet pancreas in the pathogenesis of human T1D.

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.

BET inhibition reforms the immune microenvironment and alleviates T cell dysfunction in chronic lymphocytic leukemia

Redundant tumor microenvironment (TME) immunosuppressive mechanisms and epigenetic maintenance of terminal T cell exhaustion greatly hinder functional antitumor immune responses in chronic lymphocytic leukemia (CLL). Bromodomain and extraterminal (BET) proteins regulate key pathways contributing to CLL pathogenesis and TME interactions, including T cell function and differentiation. Herein, we report that blocking BET protein function alleviates immunosuppressive networks in the CLL TME and repairs inherent CLL T cell defects. The pan-BET inhibitor OPN-51107 reduced exhaustion-associated cell signatures resulting in improved T cell proliferation and effector function in the Eμ-TCL1 splenic TME. Following BET inhibition (BET-i), TME T cells coexpressed significantly fewer inhibitory receptors (IRs) (e.g., PD-1, CD160, CD244, LAG3, VISTA). Complementary results were witnessed in primary CLL cultures, wherein OPN-51107 exerted proinflammatory effects on T cells, regardless of leukemic cell burden. BET-i additionally promotes a progenitor T cell phenotype through reduced expression of transcription factors that maintain terminal differentiation and increased expression of TCF-1, at least in part through altered chromatin accessibility. Moreover, direct T cell effects of BET-i were unmatched by common targeted therapies in CLL. This study demonstrates the immunomodulatory action of BET-i on CLL T cells and supports the inclusion of BET inhibitors in the management of CLL to alleviate terminal T cell dysfunction and potentially enhance tumoricidal T cell activity.

Microfluidic Biochips for Single-Cell Isolation and Single-Cell Analysis of Multiomics and Exosomes

Single-cell multiomic and exosome analyses are potent tools in various fields, such as cancer research, immunology, neuroscience, microbiology, and drug development. They facilitate the in-depth exploration of biological systems, providing insights into disease mechanisms and aiding in treatment. Single-cell isolation, which is crucial for single-cell analysis, ensures reliable cell isolation and quality control for further downstream analyses. Microfluidic chips are small lightweight systems that facilitate efficient and high-throughput single-cell isolation and real-time single-cell analysis on- or off-chip. Therefore, most current single-cell isolation and analysis technologies are based on the single-cell microfluidic technology. This review offers comprehensive guidance to researchers across different fields on the selection of appropriate microfluidic chip technologies for single-cell isolation and analysis. This review describes the design principles, separation mechanisms, chip characteristics, and cellular effects of various microfluidic chips available for single-cell isolation. Moreover, this review highlights the implications of using this technology for subsequent analyses, including single-cell multiomic and exosome analyses. Finally, the current challenges and future prospects of microfluidic chip technology are outlined for multiplex single-cell isolation and multiomic and exosome analyses.

Bidirectional epigenetic editing reveals hierarchies in gene regulation

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

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.

A review of single-cell transcriptomics and epigenomics studies in maternal and child health

Single-cell sequencing technologies enhance our understanding of cellular dynamics throughout pregnancy. We outlined the workflow of single-cell sequencing techniques and reviewed single-cell studies in maternal and child health. We conducted a literature review of single cell studies on maternal and child health using PubMed. We summarized the findings from 16 single-cell atlases of the human and mammalian placenta across gestational stages and 31 single-cell studies on maternal exposures and complications including infection, obesity, diet, gestational diabetes, pre-eclampsia, environmental exposure and preterm birth. Single-cell studies provides insights on novel cell types in placenta and cell type-specific marks associated with maternal exposures and complications.

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.

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.

Mapping genotypes to chromatin accessibility profiles in single cells

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

RAG suppresses group 2 innate lymphoid cells

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

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.

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.

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.