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

Single-cell multi-omics sequencing reveals the immunological disturbance underlying STAT3-V637M Hyper-IgE syndrome

Hyper-IgE syndrome (HIES) is a primary immunodeficiency characterized by, among others, the excessive production of IgE and repetitive bacterial/fungal infections. Mutations in STAT3, a transcription factor that orchestrates immune responses, may cause HIES, but the underlying mechanisms are not fully understood. Here, we used multi-omic approaches to comprehensively decipher the immune disturbance in a male HIES patient harboring STAT3-V637M. In his peripheral blood mononuclear cell (PBMC) we found significant clonal expansion of CD8 T cells (with increased CD8 subunits expression, potentially enhancing responsiveness to MHC I molecules), but not in his CD4 T cells and B cells. Although his B cells exhibited a higher potential in producing immunoglobulin, elevated SPIC binding might bias the products toward IgE isotype. Immune checkpoint inhibitors, including CTLA4, LAG3, were overexpressed in his PBMC-CD4 T cells, accompanied by reduced CD28 and IL6ST (gp130) expression. In his CD4 T cells, integrative analyses predicted upstream transcription factors (including ETV6, KLF13, and RORA) for LAG3, IL6ST, and CD28, respectively. The down-regulation of phagocytosis and nitric oxide synthesis-related genes in his PBMC-monocytes seem to be the culprit of his disseminated bacterial/fungal infection. Counterintuitively, in his PBMC we predicted increased STAT3 binding in both naïve and mature CD4 compartments, although this was not observed in most of his PBMC. In his bronchoalveolar lavage fluid (BALF), we found two macrophage subtypes with anti-bacterial properties, which were identified by CXCL8/S100A8/S100A9, or SOD2, respectively. Together, we described how the immune cell landscape was disturbed in STAT3-V637M HIES, providing a resource for further studies.

Single-cell multiomic understanding of HIV-1 reservoir at epigenetic, transcriptional, and protein levels

PURPOSE OF REVIEW

The success of HIV-1 eradication strategies relies on in-depth understanding of HIV-1-infected cells. However, HIV-1-infected cells are extremely heterogeneous and rare. Single-cell multiomic approaches resolve the heterogeneity and rarity of HIV-1-infected cells.

RECENT FINDINGS

Advancement in single-cell multiomic approaches enabled HIV-1 reservoir profiling across the epigenetic (ATAC-seq), transcriptional (RNA-seq), and protein levels (CITE-seq). Using HIV-1 RNA as a surrogate, ECCITE-seq identified enrichment of HIV-1-infected cells in clonally expanded cytotoxic CD4+ T cells. Using HIV-1 DNA PCR-activated microfluidic sorting, FIND-seq captured the bulk transcriptome of HIV-1 DNA+ cells. Using targeted HIV-1 DNA amplification, PheP-seq identified surface protein expression of intact versus defective HIV-1-infected cells. Using ATAC-seq to identify HIV-1 DNA, ASAP-seq captured transcription factor activity and surface protein expression of HIV-1 DNA+ cells. Combining HIV-1 mapping by ATAC-seq and HIV-1 RNA mapping by RNA-seq, DOGMA-seq captured the epigenetic, transcriptional, and surface protein expression of latent and transcriptionally active HIV-1-infected cells. To identify reproducible biological insights and authentic HIV-1-infected cells and avoid false-positive discovery of artifacts, we reviewed current practices of single-cell multiomic experimental design and bioinformatic analysis.

SUMMARY

Single-cell multiomic approaches may identify innovative mechanisms of HIV-1 persistence, nominate therapeutic strategies, and accelerate discoveries.

Chromatin accessibility dynamics in colorectal cancer liver metastasis: Uncovering the liver tropism at single cell resolution

Tumor metastasis causes over 90% of cancer related death and no currently available therapies target it. However, there is limited understanding regarding the epigenetic regulation of genes during this complex process. Here by integrating single-cell ATAC-seq (scATAC-seq), single-cell RNA-seq (scRNA-seq), microarray, bulk RNA-seq, immunohistochemistry (IHC) staining, as well as proteomics datasets from paired primary and liver metastatic colorectal cancer (CRC) patient-derived xenograft (PDX) model and patients, we discovered that liver metastatic CRC cells lose their colon-specific chromatin accessible sites yet gain liver-specific ones. Importantly, we observed elevated accessibility of HNF4A, a liver-specific transcription factor, in liver metastatic CRC cells. Subsequently, we performed clustering analysis of liver metastatic CRC cells together with cells involved in liver development, revealing significant heterogeneity among the liver metastatic CRC cells. Over 50% of the liver metastatic CRC cells exhibited characteristics similar to those of erythroid progenitors and hepatocytes, showing increased expression of genes involved in oxidative phosphorylation and glycolysis. Moreover, our discovery further revealed that the MHC and IFN response genes in these cells exhibit moderate epigenetic activity, which is significantly associated with the low objective response rates in checkpoint blockade immunotherapy. Our findings uncovered the critical roles of HNF4A and the cell populations within liver metastatic CRC cells might serve as crucial therapeutic targets for addressing liver metastasis and improving the immunotherapy response in patients with CRC.

KDM6A epigenetically regulates subtype plasticity in small cell lung cancer

Small cell lung cancer (SCLC) exists broadly in four molecular subtypes: ASCL1, NEUROD1, POU2F3 and Inflammatory. Initially, SCLC subtypes were thought to be mutually exclusive, but recent evidence shows intra-tumoural subtype heterogeneity and plasticity between subtypes. Here, using a CRISPR-based autochthonous SCLC genetically engineered mouse model to study the consequences of KDM6A/UTX inactivation, we show that KDM6A inactivation induced plasticity from ASCL1 to NEUROD1 resulting in SCLC tumours that express both ASCL1 and NEUROD1. Mechanistically, KDM6A normally maintains an active chromatin state that favours the ASCL1 subtype with its loss decreasing H3K4me1 and increasing H3K27me3 at enhancers of neuroendocrine genes leading to a cell state that is primed for ASCL1-to-NEUROD1 subtype switching. This work identifies KDM6A as an epigenetic regulator that controls ASCL1 to NEUROD1 subtype plasticity and provides an autochthonous SCLC genetically engineered mouse model to model ASCL1 and NEUROD1 subtype heterogeneity and plasticity, which is found in 35-40% of human SCLCs.

Chromatin accessibility memory of donor cells disrupts bovine somatic cell nuclear transfer blastocysts development

The post-transfer developmental capacity of bovine somatic cell nuclear transfer (SCNT) blastocysts is reduced, implying that abnormalities in gene expression regulation are present at blastocyst stage. Chromatin accessibility, as an indicator for transcriptional regulatory elements mediating gene transcription activity, has heretofore been largely unexplored in SCNT embryos, especially at blastocyst stage. In the present study, single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) of in vivo and SCNT blastocysts were conducted to segregate lineages and demonstrate the aberrant chromatin accessibility of transcription factors (TFs) related to inner cell mass (ICM) development in SCNT blastocysts. Pseudotime analysis of lineage segregation further reflected dysregulated chromatin accessibility dynamics of TFs in the ICM of SCNT blastocysts compared to their in vivo counterparts. ATAC- and ChIP-seq results of SCNT donor cells revealed that the aberrant chromatin accessibility in the ICM of SCNT blastocysts was due to the persistence of chromatin accessibility memory at corresponding loci in the donor cells, with strong enrichment of trimethylation of histone H3 at lysine 4 (H3K4me3) at these loci. Correction of the aberrant chromatin accessibility through demethylation of H3K4me3 by KDM5B diminished the expression of related genes (e.g., BCL11B) and significantly improved the ICM proliferation in SCNT blastocysts. This effect was confirmed by knocking down BCL11B in SCNT embryos to down-regulate p21 and alleviate the inhibition of ICM proliferation. These findings expand our understanding of the chromatin accessibility abnormalities in SCNT blastocysts and BCL11B may be a potential target to improve SCNT efficiency.

Epigenetic regulation of HBV-specific tumor-infiltrating T cells in HBV-related HCC

BACKGROUND AND AIMS

HBV shapes the T-cell immune responses in HBV-related HCC. T cells can be recruited to the nidus, but limited T cells participate specifically in response to the HBV-related tumor microenvironment and HBV antigens. How epigenomic programs regulate T-cell compartments in virus-specific immune processes is unclear.

APPROACH AND RESULTS

We developed Ti-ATAC-seq. 2 to map the T-cell receptor repertoire, epigenomic, and transcriptomic landscape of αβ T cells at both the bulk-cell and single-cell levels in 54 patients with HCC. We deeply investigated HBV-specific T cells and HBV-related T-cell subsets that specifically responded to HBV antigens and the HBV + tumor microenvironment, respectively, characterizing their T-cell receptor clonality and specificity and performing epigenomic profiling. A shared program comprising NFKB1/2-, Proto-Oncogene, NF-KB Sub unit, NFATC2-, and NR4A1-associated unique T-cell receptor-downstream core epigenomic and transcriptomic regulome commonly regulated the differentiation of HBV-specific regulatory T-cell (Treg) cells and CD8 + exhausted T cells; this program was also selectively enriched in the HBV-related Treg-CTLA4 and CD8-exhausted T cell-thymocyte selection associated high mobility subsets and drove greater clonal expansion in HBV-related Treg-CTLA4 subset. Overall, 54% of the effector and memory HBV-specific T cells are governed by transcription factor motifs of activator protein 1, NFE2, and BACH1/2, which have been reported to be associated with prolonged patient relapse-free survival. Moreover, HBV-related tumor-infiltrating Tregs correlated with both increased viral titer and poor prognosis in patients.

CONCLUSIONS

This study provides insight into the cellular and molecular basis of the epigenomic programs that regulate the differentiation and generation of HBV-related T cells from viral infection and HBV + HCC unique immune exhaustion.

Exploring the impact of dexamethasone on gene regulation in myeloma cells

Among glucocorticoids (GCs), dexamethasone (Dex) is widely used in treatment of multiple myelomas. However, despite a definite benefit, all patients relapse. Moreover, the molecular basis of glucocorticoid efficacy remains elusive. To determine genomic response to Dex in myeloma cells, we generated bulk and single-cell multi-omics data and high-resolution contact maps of active enhancers and target genes. We show that a minority of glucocorticoid receptor-binding sites are associated with enhancer activity gains, increased interaction loops, and transcriptional activity. We identified and characterized a predominant enhancer enriched in cohesin (RAD21) and more accessible upon Dex exposure. Analysis of four gene-specific networks revealed the importance of the CTCF-cohesin couple and the synchronization of regulatory sequence openings for efficient transcription in response to Dex. Notably, these epigenomic changes are associated with cell-to-cell transcriptional heterogeneity, in particular, lineage-specific genes. As consequences, BCL2L11-encoding BIM critical for Dex-induced apoptosis and CXCR4 protective from chemotherapy-induced apoptosis are rather up-regulated in different cells. In summary, our work provides new insights into the molecular mechanisms involved in Dex escape.

A Roadmap for a Consensus Human Skin Cell Atlas and Single-Cell Data Standardization

Single-cell technologies have become essential to driving discovery in both basic and translational investigative dermatology. Despite the multitude of available datasets, a central reference atlas of normal human skin, which can serve as a reference resource for skin cell types, cell states, and their molecular signatures, is still lacking. For any such atlas to receive broad acceptance, participation by many investigators during atlas construction is an essential prerequisite. As part of the Human Cell Atlas project, we have assembled a Skin Biological Network to build a consensus Human Skin Cell Atlas and outline a roadmap toward that goal. We define the drivers of skin diversity to be considered when selecting sequencing datasets for the atlas and list practical hurdles during skin sampling that can result in data gaps and impede comprehensive representation and technical considerations for tissue processing and computational analysis, the accounting for which should minimize biases in cell type enrichments and exclusions and decrease batch effects. By outlining our goals for Atlas 1.0, we discuss how it will uncover new aspects of skin biology.

Ancestry-related differences in chromatin accessibility and gene expression of APOE ε4 are associated with Alzheimer's disease risk

INTRODUCTION

European local ancestry (ELA) surrounding apolipoprotein E (APOE) ε4 confers higher risk for Alzheimer's disease (AD) compared to African local ancestry (ALA). We demonstrated significantly higher APOE ε4  expression in ELA versus ALA in AD brains from APOE ε4/ε4 carriers. Chromatin accessibility differences could contribute to these expression changes.

METHODS

We performed single nuclei assays for transposase accessible chromatin sequencing from the frontal cortex of six ALA and six ELA AD brains, homozygous for local ancestry and APOE ε4.

RESULTS

Our results showed an increased chromatin accessibility at the APOE ε4  promoter area in ELA versus ALA astrocytes. This increased accessibility in ELA astrocytes extended genome wide. Genes with increased accessibility in ELA in astrocytes were enriched for synapsis, cholesterol processing, and astrocyte reactivity.

DISCUSSION

Our results suggest that increased chromatin accessibility of APOE ε4  in ELA astrocytes contributes to the observed elevated APOE ε4  expression, corresponding to the increased AD risk in ELA versus ALA APOE ε4/ε4 carriers.

Epstein-Barr virus evades restrictive host chromatin closure by subverting B cell activation and germinal center regulatory loci

Chromatin accessibility fundamentally governs gene expression and biological response programs that can be manipulated by pathogens. Here we capture dynamic chromatin landscapes of individual B cells during Epstein-Barr virus (EBV) infection. EBV+ cells that exhibit arrest via antiviral sensing and proliferation-linked DNA damage experience global accessibility reduction. Proliferative EBV+ cells develop expression-linked architectures and motif accessibility profiles resembling in vivo germinal center (GC) phenotypes. Remarkably, EBV elicits dark zone (DZ), light zone (LZ), and post-GC B cell chromatin features despite BCL6 downregulation. Integration of single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq), single-cell RNA sequencing (scRNA-seq), and chromatin immunoprecipitation sequencing (ChIP-seq) data enables genome-wide cis-regulatory predictions implicating EBV nuclear antigens (EBNAs) in phenotype-specific control of GC B cell activation, survival, and immune evasion. Knockouts validate bioinformatically identified regulators (MEF2C and NFE2L2) of EBV-induced GC phenotypes and EBNA-associated loci that regulate gene expression (CD274/PD-L1). These data and methods can inform high-resolution investigations of EBV-host interactions, B cell fates, and virus-mediated lymphomagenesis.

An introduction to representation learning for single-cell data analysis

Single-cell-resolved systems biology methods, including omics- and imaging-based measurement modalities, generate a wealth of high-dimensional data characterizing the heterogeneity of cell populations. Representation learning methods are routinely used to analyze these complex, high-dimensional data by projecting them into lower-dimensional embeddings. This facilitates the interpretation and interrogation of the structures, dynamics, and regulation of cell heterogeneity. Reflecting their central role in analyzing diverse single-cell data types, a myriad of representation learning methods exist, with new approaches continually emerging. Here, we contrast general features of representation learning methods spanning statistical, manifold learning, and neural network approaches. We consider key steps involved in representation learning with single-cell data, including data pre-processing, hyperparameter optimization, downstream analysis, and biological validation. Interdependencies and contingencies linking these steps are also highlighted. This overview is intended to guide researchers in the selection, application, and optimization of representation learning strategies for current and future single-cell research applications.

Single-Cell Transcriptional and Epigenetic Profiles of Male Breast Cancer Nominate Salient Cancer-Specific Enhancers

Male breast cancer represents about 1% of all breast cancer diagnoses and, although there are some similarities between male and female breast cancer, the paucity of data available on male breast cancer makes it difficult to establish targeted therapies. To date, most male breast cancers (MBCs) are treated according to protocols established for female breast cancer (FBC). Thus, defining the transcriptional and epigenetic landscape of MBC with improved resolution is critical for developing better avenues for therapeutic intervention. In this study, we present matched transcriptional (scRNA-seq) and epigenetic (scATAC-seq) profiles at single-cell resolution of two treatment naïve MBC tumors processed immediately after surgical resection. These data enable the detection of differentially expressed genes between male and female breast tumors across immune, stromal, and malignant cell types, to highlight several genes that may have therapeutic implications. Notably, MYC target genes and mTORC1 signaling genes were significantly upregulated in the malignant cells of MBC compared to the female counterparts. To understand how the regulatory landscape of MBC gives rise to these male-specific gene expression patterns, we leveraged the scATAC-seq data to systematically link changes in chromatin accessibility to changes in gene expression within each cell type. We observed cancer-specific rewiring of several salient enhancers and posit that these enhancers have a higher regulatory load than lineage-specific enhancers. We highlight two examples of previously unannotated cancer-cell-specific enhancers of ANXA2 and PRDX4 gene expression and show evidence for super-enhancer regulation of LAMB3 and CD47 in male breast cancer cells. Overall, this dataset annotates clinically relevant regulatory networks in male breast tumors, providing a useful resource that expands our current understanding of the gene expression programs that underlie the biology of MBC.

Temporally specific gene expression and chromatin remodeling programs regulate a conserved Pdyn enhancer

Neuronal and behavioral adaptations to novel stimuli are regulated by temporally dynamic waves of transcriptional activity, which shape neuronal function and guide enduring plasticity. Neuronal activation promotes expression of an immediate early gene (IEG) program comprised primarily of activity-dependent transcription factors, which are thought to regulate a second set of late response genes (LRGs). However, while the mechanisms governing IEG activation have been well studied, the molecular interplay between IEGs and LRGs remain poorly characterized. Here, we used transcriptomic and chromatin accessibility profiling to define activity-driven responses in rat striatal neurons. As expected, neuronal depolarization generated robust changes in gene expression, with early changes (1 h) enriched for inducible transcription factors and later changes (4 h) enriched for neuropeptides, synaptic proteins, and ion channels. Remarkably, while depolarization did not induce chromatin remodeling after 1 h, we found broad increases in chromatin accessibility at thousands of sites in the genome at 4 h after neuronal stimulation. These putative regulatory elements were found almost exclusively at non-coding regions of the genome, and harbored consensus motifs for numerous activity-dependent transcription factors such as AP-1. Furthermore, blocking protein synthesis prevented activity-dependent chromatin remodeling, suggesting that IEG proteins are required for this process. Targeted analysis of LRG loci identified a putative enhancer upstream of Pdyn (prodynorphin), a gene encoding an opioid neuropeptide implicated in motivated behavior and neuropsychiatric disease states. CRISPR-based functional assays demonstrated that this enhancer is both necessary and sufficient for Pdyn transcription. This regulatory element is also conserved at the human PDYN locus, where its activation is sufficient to drive PDYN transcription in human cells. These results suggest that IEGs participate in chromatin remodeling at enhancers and identify a conserved enhancer that may act as a therapeutic target for brain disorders involving dysregulation of Pdyn.

Resolving the spatial architecture of myeloma and its microenvironment at the single-cell level

In multiple myeloma spatial differences in the subclonal architecture, molecular signatures and composition of the microenvironment remain poorly characterized. To address this shortcoming, we perform multi-region sequencing on paired random bone marrow and focal lesion samples from 17 newly diagnosed patients. Using single-cell RNA- and ATAC-seq we find a median of 6 tumor subclones per patient and unique subclones in focal lesions. Genetically identical subclones display different levels of spatial transcriptional plasticity, including nearly identical profiles and pronounced heterogeneity at different sites, which can include differential expression of immunotherapy targets, such as CD20 and CD38. Macrophages are significantly depleted in the microenvironment of focal lesions. We observe proportional changes in the T-cell repertoire but no site-specific expansion of T-cell clones in intramedullary lesions. In conclusion, our results demonstrate the relevance of considering spatial heterogeneity in multiple myeloma with potential implications for models of cell-cell interactions and disease progression.

Single-cell genomics improves the discovery of risk variants and genes of atrial fibrillation

Genome-wide association studies (GWAS) have linked hundreds of loci to cardiac diseases. However, in most loci the causal variants and their target genes remain unknown. We developed a combined experimental and analytical approach that integrates single cell epigenomics with GWAS to prioritize risk variants and genes. We profiled accessible chromatin in single cells obtained from human hearts and leveraged the data to study genetics of Atrial Fibrillation (AF), the most common cardiac arrhythmia. Enrichment analysis of AF risk variants using cell-type-resolved open chromatin regions (OCRs) implicated cardiomyocytes as the main mediator of AF risk. We then performed statistical fine-mapping, leveraging the information in OCRs, and identified putative causal variants in 122 AF-associated loci. Taking advantage of the fine-mapping results, our novel statistical procedure for gene discovery prioritized 46 high-confidence risk genes, highlighting transcription factors and signal transduction pathways important for heart development. In summary, our analysis provides a comprehensive map of AF risk variants and genes, and a general framework to integrate single-cell genomics with genetic studies of complex traits.

Whole genome deconvolution unveils Alzheimer's resilient epigenetic signature

Assay for Transposase Accessible Chromatin by sequencing (ATAC-seq) accurately depicts the chromatin regulatory state and altered mechanisms guiding gene expression in disease. However, bulk sequencing entangles information from different cell types and obscures cellular heterogeneity. To address this, we developed Cellformer, a deep learning method that deconvolutes bulk ATAC-seq into cell type-specific expression across the whole genome. Cellformer enables cost-effective cell type-specific open chromatin profiling in large cohorts. Applied to 191 bulk samples from 3 brain regions, Cellformer identifies cell type-specific gene regulatory mechanisms involved in resilience to Alzheimer's disease, an uncommon group of cognitively healthy individuals that harbor a high pathological load of Alzheimer's disease. Cell type-resolved chromatin profiling unveils cell type-specific pathways and nominates potential epigenetic mediators underlying resilience that may illuminate therapeutic opportunities to limit the cognitive impact of the disease. Cellformer is freely available to facilitate future investigations using high-throughput bulk ATAC-seq data.

Systematic benchmarking of single-cell ATAC-sequencing protocols

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

Single-Nucleus Profiling Identifies Accelerated Oligodendrocyte Precursor Cell Senescence in a Mouse Model of Down Syndrome

Down syndrome (DS), the most common genetic cause of intellectual disability, is associated with lifelong cognitive deficits. However, the mechanisms by which triplication of chromosome 21 genes drive neuroinflammation and cognitive dysfunction are poorly understood. Here, using the Ts65Dn mouse model of DS, we performed an integrated single-nucleus ATAC and RNA-sequencing (snATAC-seq and snRNA-seq) analysis of the adult cortex. We identified cell type-specific transcriptional and chromatin-associated changes in the Ts65Dn cortex, including regulators of neuroinflammation, transcription and translation, myelination, and mitochondrial function. We discovered enrichment of a senescence-associated transcriptional signature in Ts65Dn oligodendrocyte (OL) precursor cells (OPCs) and epigenetic changes consistent with a loss of heterochromatin. We found that senescence is restricted to a subset of OPCs concentrated in deep cortical layers. Treatment of Ts65Dn mice with a senescence-reducing flavonoid rescued cortical OPC proliferation, restored microglial homeostasis, and improved contextual fear memory. Together, these findings suggest that cortical OPC senescence may be an important driver of neuropathology in DS.

Circulating immune cell dynamics as outcome predictors for immunotherapy in non-small cell lung cancer

The use of immune checkpoint inhibitors (ICIs) continues to transform the therapeutic landscape of non-small cell lung cancer (NSCLC), with these drugs now being evaluated at every stage of the disease. In contrast to these advances, little progress has been made with respect to reliable predictive biomarkers that can inform clinicians on therapeutic efficacy. All current biomarkers for outcome prediction, including PD-L1, tumor mutational burden or complex immune gene expression signatures, require access to tumor tissue. Besides the invasive nature of the sampling procedure, other disadvantages of tumor tissue biopsies are the inability to capture the complete spatial heterogeneity of the tumor and the difficulty to perform longitudinal follow-up on treatment. A concept emerges in which systemic immune events developing at a distance from the tumor reflect local response or resistance to immunotherapy. The importance of this cancer 'macroenvironment', which can be deciphered by comprehensive analysis of peripheral blood immune cell subsets, has been demonstrated in several cutting-edge preclinical reports, and is corroborated by intriguing data emerging from ICI-treated patients. In this review, we will provide the biological rationale underlying the potential of blood immune cell-based biomarkers in guiding treatment decision in immunotherapy-eligible NSCLC patients. Finally, we will describe new techniques that will facilitate the discovery of more immune cell subpopulations with potential to become predictive biomarkers, and reflect on ways and the remaining challenges to bring this type of analysis to the routine clinical care in the near future.

Solid-phase capture and profiling of open chromatin by spatial ATAC

Current methods for epigenomic profiling are limited in their ability to obtain genome-wide information with spatial resolution. We introduce spatial ATAC, a method that integrates transposase-accessible chromatin profiling in tissue sections with barcoded solid-phase capture to perform spatially resolved epigenomics. We show that spatial ATAC enables the discovery of the regulatory programs underlying spatial gene expression during mouse organogenesis, lineage differentiation and in human pathology.

Synergism Between IL21 and Anti-PD-1 Combination Therapy is Underpinned by the Coordinated Reprogramming of the Immune Cellular Network in the Tumor Microenvironment

T cell-stimulating cytokines and immune checkpoint inhibitors (ICI) are an ideal combination for increasing response rates of cancer immunotherapy. However, the results of clinical trials have not been satisfying. It is important to understand the mechanism of synergy between these two therapeutic modalities. Here, through integrated analysis of multiple single-cell RNA sequencing (scRNA-seq) datasets of human tumor-infiltrating immune cells, we demonstrate that IL21 is produced by tumor-associated T follicular helper cells and hyperactivated/exhausted CXCL13+CD4+ T cells in the human tumor microenvironment (TME). In the mouse model, the hyperactivated/exhausted CD4+ T cell-derived IL21 enhances the helper function of CD4+ T cells that boost CD8+ T cell-mediated immune responses during PD-1 blockade immunotherapy. In addition, we demonstrated that IL21's antitumor activity did not require T-cell trafficking. Using scRNA-seq analysis of the whole tumor-infiltrating immune cells, we demonstrated that IL21 treatment in combination with anti-PD-1 blockade synergistically drives tumor antigen-specific CD8+ T cells to undergo clonal expansion and differentiate toward the hyperactive/exhausted functional state in the TME. In addition, IL21 treatment and anti-PD-1 blockade synergistically promote dendritic cell (DC) activation and maturation to mature DC as well as monocyte to type 1 macrophage (M1) differentiation in the TME. Furthermore, the combined treatment reprograms the immune cellular network by reshaping cell-cell communication in the TME. Our study establishes unique mechanisms of synergy between IL21 and PD-1-based ICI in the TME through the coordinated promotion of type 1 immune responses.

Significance

This study reveals how cytokine and checkpoint inhibitor therapy can be combined to increase the efficacy of cancer immunotherapy.

A single-cell atlas of in vitro multiculture systems uncovers the in vivo lineage trajectory and cell state in the human lung

We present an in-depth single-cell atlas of in vitro multiculture systems on human primary airway epithelium derived from normal and diseased lungs of 27 individual donors. Our large-scale single-cell profiling identified new cell states and differentiation trajectories of rare airway epithelial cell types in human distal lungs. By integrating single-cell datasets of human lung tissues, we discovered immune-primed subsets enriched in lungs and organoids derived from patients with chronic respiratory disease. To demonstrate the full potential of our platform, we further illustrate transcriptomic responses to various respiratory virus infections in vitro airway models. Our work constitutes a single-cell roadmap for the cellular and molecular characteristics of human primary lung cells in vitro and their relevance to human tissues in vivo.

MultiVI: deep generative model for the integration of multimodal data

Jointly profiling the transcriptome, chromatin accessibility and other molecular properties of single cells offers a powerful way to study cellular diversity. Here we present MultiVI, a probabilistic model to analyze such multiomic data and leverage it to enhance single-modality datasets. MultiVI creates a joint representation that allows an analysis of all modalities included in the multiomic input data, even for cells for which one or more modalities are missing. It is available at scvi-tools.org .

Interpreting non-coding disease-associated human variants using single-cell epigenomics

Genome-wide association studies (GWAS) have linked hundreds of thousands of sequence variants in the human genome to common traits and diseases. However, translating this knowledge into a mechanistic understanding of disease-relevant biology remains challenging, largely because such variants are predominantly in non-protein-coding sequences that still lack functional annotation at cell-type resolution. Recent advances in single-cell epigenomics assays have enabled the generation of cell type-, subtype- and state-resolved maps of the epigenome in heterogeneous human tissues. These maps have facilitated cell type-specific annotation of candidate cis-regulatory elements and their gene targets in the human genome, enhancing our ability to interpret the genetic basis of common traits and diseases.

A fast chemical reprogramming system promotes cell identity transition through a diapause-like state

Cellular reprogramming by only small molecules holds enormous potentials for regenerative medicine. However, chemical reprogramming remains a slow process and labour intensive, hindering its broad applications and the investigation of underlying molecular mechanisms. Here, through screening of over 21,000 conditions, we develop a fast chemical reprogramming (FCR) system, which significantly improves the kinetics of cell identity rewiring. We find that FCR rapidly goes through an interesting route for pluripotent reprogramming, uniquely transitioning through a developmentally diapause-like state. Furthermore, FCR critically enables comprehensive characterizations using multi-omics technologies, and has revealed unexpected important features including key regulatory factors and epigenetic dynamics. Particularly, activation of pluripotency-related endogenous retroviruses via inhibition of heterochromatin significantly enhances reprogramming. Our studies provide critical insights into how only environmental cues are sufficient to rapidly reinstate pluripotency in somatic cells, and make notable technical and conceptual advances for solving the puzzle of regeneration.

PD-1/PD-L1 blockade restores tumor-induced COVID-19 vaccine bluntness

The COVID-19 vaccinations are crucial in protecting against the global pandemic. However, accumulating studies revealed the severely blunted COVID-19 vaccine effectiveness in cancer patients. The PD-1/PD-L1 immune checkpoint blockade (ICB) therapy leads to durable therapeutic responses in a subset of cancer patients and has been approved to treat a wide spectrum of cancers in the clinic. In this regard, it is pivotal to explore the potential impact of PD-1/PD-L1 ICB therapy on COVID-19 vaccine effectiveness during ongoing malignancy. In this study, using preclinical models, we found that the tumor-suppressed COVID-19 vaccine responses are largely reverted in the setting of PD-1/PD-L1 ICB therapy. We also identified that the PD-1/PD-L1 blockade-directed restoration of COVID-19 vaccine effectiveness is irrelevant to anti-tumor therapeutic outcomes. Mechanistically, the restored COVID-19 vaccine effectiveness is entwined with the PD-1/PD-L1 blockade-driven preponderance of follicular helper T cell and germinal center responses during ongoing malignancy. Thus, our findings indicate that PD-1/PD-L1 blockade will greatly normalize the responses of cancer patients to COVID-19 vaccination, while regardless of its anti-tumor efficacies on these patients.

A Review of Single-Cell RNA-Seq Annotation, Integration, and Cell-Cell Communication

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for investigating cellular biology at an unprecedented resolution, enabling the characterization of cellular heterogeneity, identification of rare but significant cell types, and exploration of cell-cell communications and interactions. Its broad applications span both basic and clinical research domains. In this comprehensive review, we survey the current landscape of scRNA-seq analysis methods and tools, focusing on count modeling, cell-type annotation, data integration, including spatial transcriptomics, and the inference of cell-cell communication. We review the challenges encountered in scRNA-seq analysis, including issues of sparsity or low expression, reliability of cell annotation, and assumptions in data integration, and discuss the potential impact of suboptimal clustering and differential expression analysis tools on downstream analyses, particularly in identifying cell subpopulations. Finally, we discuss recent advancements and future directions for enhancing scRNA-seq analysis. Specifically, we highlight the development of novel tools for annotating single-cell data, integrating and interpreting multimodal datasets covering transcriptomics, epigenomics, and proteomics, and inferring cellular communication networks. By elucidating the latest progress and innovation, we provide a comprehensive overview of the rapidly advancing field of scRNA-seq analysis.

Application of single-cell omics in inflammatory bowel disease

Over the past decade, the advent of single cell RNA-sequencing has revolutionized the approach in cellular transcriptomics research. The current technology offers an unbiased platform to understand how genotype correlates to phenotype. Single-cell omics applications in gastrointestinal (GI) research namely inflammatory bowel disease (IBD) has become popular in the last few years with multiple publications as single-cell omics techniques can be applied directly to the target organ, the GI tract at the tissue level. Through examination of mucosal tissue and peripheral blood in IBD, the recent boom in single cell research has identified a myriad of key immune players from enterocytes to tissue resident memory T cells, and explored functional heterogeneity within cellular subsets previously unreported. As we begin to unravel the complex mucosal immune system in states of health and disease like IBD, the power of exploration through single-cell omics can change our approach to translational research. As novel techniques evolve through multiplexing single-cell omics and spatial transcriptomics come to the forefront, we can begin to fully comprehend the disease IBD and better design targets of treatment. In addition, hopefully these techniques can ultimately begin to identify biomarkers of therapeutic response and answer clinically relevant questions in how to tailor individual therapy to patients through personalized medicine.

ChromaFold predicts the 3D contact map from single-cell chromatin accessibility

The identification of cell-type-specific 3D chromatin interactions between regulatory elements can help to decipher gene regulation and to interpret the function of disease-associated non-coding variants. However, current chromosome conformation capture (3C) technologies are unable to resolve interactions at this resolution when only small numbers of cells are available as input. We therefore present ChromaFold, a deep learning model that predicts 3D contact maps and regulatory interactions from single-cell ATAC sequencing (scATAC-seq) data alone. ChromaFold uses pseudobulk chromatin accessibility, co-accessibility profiles across metacells, and predicted CTCF motif tracks as input features and employs a lightweight architecture to enable training on standard GPUs. Once trained on paired scATAC-seq and Hi-C data in human cell lines and tissues, ChromaFold can accurately predict both the 3D contact map and peak-level interactions across diverse human and mouse test cell types. In benchmarking against a recent deep learning method that uses bulk ATAC-seq, DNA sequence, and CTCF ChIP-seq to make cell-type-specific predictions, ChromaFold yields superior prediction performance when including CTCF ChIP-seq data as an input and comparable performance without. Finally, fine-tuning ChromaFold on paired scATAC-seq and Hi-C in a complex tissue enables deconvolution of chromatin interactions across cell subpopulations. ChromaFold thus achieves state-of-the-art prediction of 3D contact maps and regulatory interactions using scATAC-seq alone as input data, enabling accurate inference of cell-type-specific interactions in settings where 3C-based assays are infeasible.

Single-cell chromatin accessibility profiling of acute myeloid leukemia reveals heterogeneous lineage composition upon therapy-resistance

Acute myeloid leukemia (AML) is a heterogeneous disease characterized by high rate of therapy resistance. Since the cell of origin can impact response to therapy, it is crucial to understand the lineage composition of AML cells at time of therapy resistance. Here we leverage single-cell chromatin accessibility profiling of 22 AML bone marrow aspirates from eight patients at time of therapy resistance and following subsequent therapy to characterize their lineage landscape. Our findings reveal a complex lineage architecture of therapy-resistant AML cells that are primed for stem and progenitor lineages and spanning quiescent, activated and late stem cell/progenitor states. Remarkably, therapy-resistant AML cells are also composed of cells primed for differentiated myeloid, erythroid and even lymphoid lineages. The heterogeneous lineage composition persists following subsequent therapy, with early progenitor-driven features marking unfavorable prognosis in The Cancer Genome Atlas AML cohort. Pseudotime analysis further confirms the vast degree of heterogeneity driven by the dynamic changes in chromatin accessibility. Our findings suggest that therapy-resistant AML cells are characterized not only by stem and progenitor states, but also by a continuum of differentiated cellular lineages. The heterogeneity in lineages likely contributes to their therapy resistance by harboring different degrees of lineage-specific susceptibilities to therapy.

Coordinated inheritance of extrachromosomal DNA species in human cancer cells

The chromosomal theory of inheritance has dominated human genetics, including cancer genetics. Genes on the same chromosome segregate together while genes on different chromosomes assort independently, providing a fundamental tenet of Mendelian inheritance. Extrachromosomal DNA (ecDNA) is a frequent event in cancer that drives oncogene amplification, dysregulated gene expression and intratumoral heterogeneity, including through random segregation during cell division. Distinct ecDNA sequences, herein termed ecDNA species, can co-exist to facilitate intermolecular cooperation in cancer cells. However, how multiple ecDNA species within a tumor cell are assorted and maintained across somatic cell generations to drive cancer cell evolution is not known. Here we show that cooperative ecDNA species can be coordinately inherited through mitotic co-segregation. Imaging and single-cell analyses show that multiple ecDNAs encoding distinct oncogenes co-occur and are correlated in copy number in human cancer cells. EcDNA species are coordinately segregated asymmetrically during mitosis, resulting in daughter cells with simultaneous copy number gains in multiple ecDNA species prior to any selection. Computational modeling reveals the quantitative principles of ecDNA co-segregation and co-selection, predicting their observed distributions in cancer cells. Finally, we show that coordinated inheritance of ecDNAs enables co-amplification of specialized ecDNAs containing only enhancer elements and guides therapeutic strategies to jointly deplete cooperating ecDNA oncogenes. Coordinated inheritance of ecDNAs confers stability to oncogene cooperation and novel gene regulatory circuits, allowing winning combinations of epigenetic states to be transmitted across cell generations.

The current state and future of T-cell exhaustion research

'Exhaustion' is a term used to describe a state of native and redirected T-cell hypo-responsiveness resulting from persistent antigen exposure during chronic viral infections or cancer. Although a well-established phenotype across mice and humans, exhaustion at the molecular level remains poorly defined and inconsistent across the literature. This is, in part, due to an overreliance on surface receptors to define these cells and explain exhaustive behaviours, an incomplete understanding of how exhaustion arises, and a lack of clarity over whether exhaustion is the same across contexts, e.g. chronic viral infections versus cancer. With the development of systems-based genetic approaches such as single-cell RNA-seq and CRISPR screens applied to in vivo data, we are moving closer to a consensus view of exhaustion, although understanding how it arises remains challenging given the difficulty in manipulating the in vivo setting. Accordingly, producing and studying exhausted T-cells ex vivo are burgeoning, allowing experiments to be conducted at scale up and with high throughput. Here, we first review what is currently known about T-cell exhaustion and how it's being studied. We then discuss how improvements in their method of isolation/production and examining the impact of different microenvironmental signals and cell interactions have now become an active area of research. Finally, we discuss what the future holds for the analysis of this physiological condition and, given the diversity of ways in which exhausted cells are now being generated, propose the adoption of a unified approach to clearly defining exhaustion using a set of metabolic-, epigenetic-, transcriptional-, and activation-based phenotypic markers, that we call 'M.E.T.A'.

In Synergy: Optimizing CAR T Development and Personalizing Patient Care Using Single-Cell Technologies

Chimeric antigen receptor (CAR) T therapies hold immense promise to revolutionize cancer treatment. Nevertheless, key challenges, primarily in solid tumor settings, continue to hinder the application of this technology. Understanding CAR T-cell mechanism of action, in vivo activity, and clinical implications is essential for harnessing its full therapeutic potential. Single-cell genomics and cell engineering tools are becoming increasingly effective for the comprehensive research of complex biological systems. The convergence of these two technologies can accelerate CAR T-cell development. Here, we examine the potential of applying single-cell multiomics for the development of next-generation CAR T-cell therapies.

SIGNIFICANCE

Although CAR T-cell therapies have demonstrated remarkable clinical results in treating cancer, their effectiveness in most patients and tumor types remains limited. Single-cell technologies, which are transforming our understanding of molecular biology, provide new opportunities to overcome the challenges of CAR T-cell therapies. Given the potential of CAR T-cell therapy to tip the balance in the fight against cancer, it is important to understand how single-cell multiomic approaches can be leveraged to develop the next generations of more effective and less toxic CAR T-cell products and to provide powerful decision-making tools for clinicians to optimize treatment and improve patient outcomes.

Gene regulatory network study of rheumatoid arthritis in single-cell chromatin landscapes of peripheral blood mononuclear cells

OBJECTIVES

Assays for transposase-accessible chromatin with single-cell sequencing (scATAC-seq) contribute to the progress in epigenetic studies. The purpose of our project was to discover the transcription factors (TFs) that were involved in the pathogenesis of rheumatoid arthritis (RA) at a single-cell resolution using epigenetic technology.

METHODS

Peripheral blood mononuclear cells of seven RA patients and seven natural controls were extracted nuclei suspensions for library construction. Subsequently, scATAC-seq was performed to generate a high-resolution map of active regulatory DNA for bioinformatics analysis.

RESULTS

We obtained 22 accessible chromatin patterns. Then, 10 key TFs were involved in RA pathogenesis by regulating the activity of mitogen-activated protein kinase. Consequently, two genes (PTPRC and SPAG9) regulated by 10 key TFs were found, which may be associated with RA disease pathogenesis, and these TFs were obviously enriched in RA patients (P < .05, fold change value > 1.2). With further quantitative polymerase chain reaction validation on PTPRC and SPAG9 in monocytes, we found differential expression of these two genes, which were regulated by eight TFs [ZNF384, HNF1B, DMRTA2, MEF2A, NFE2L1, CREB3L4 (var. 2), FOSL2::JUNB (var. 2), and MEF2B], showing highly accessible binding sites in RA patients.

CONCLUSIONS

These findings demonstrate the value of using scATAC-seq to reveal transcriptional regulatory variation in RA-derived peripheral blood mononuclear cells, providing insights into therapy from an epigenetic perspective.

The Use of Single-Cell RNA-Sequencing and Spatial Transcriptomics in Understanding the Pathogenesis and Treatment of Skin Diseases

The development of multiomic profiling tools has rapidly expanded in recent years, along with their use in profiling skin tissues in various contexts, including dermatologic diseases. Among these tools, single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics (ST) have emerged as widely adopted and powerful assays for elucidating key cellular components and their spatial arrangement within skin disease. In this paper, we review the recent biological insights gained from the use of scRNA-seq and ST and the advantages of combining both for profiling skin diseases, including aberrant wound healing, inflammatory skin diseases, and cancer. We discuss the role of scRNA-seq and ST in improving skin disease treatments and moving toward the goal of achieving precision medicine in dermatology, whereby patients can be optimally matched to treatments that maximize therapeutic response.

Clonal hematopoiesis is associated with protection from Alzheimer's disease

Clonal hematopoiesis of indeterminate potential (CHIP) is a premalignant expansion of mutated hematopoietic stem cells. As CHIP-associated mutations are known to alter the development and function of myeloid cells, we hypothesized that CHIP may also be associated with the risk of Alzheimer's disease (AD), a disease in which brain-resident myeloid cells are thought to have a major role. To perform association tests between CHIP and AD dementia, we analyzed blood DNA sequencing data from 1,362 individuals with AD and 4,368 individuals without AD. Individuals with CHIP had a lower risk of AD dementia (meta-analysis odds ratio (OR) = 0.64, P = 3.8 × 10-5), and Mendelian randomization analyses supported a potential causal association. We observed that the same mutations found in blood were also detected in microglia-enriched fraction of the brain in seven of eight CHIP carriers. Single-nucleus chromatin accessibility profiling of brain-derived nuclei in six CHIP carriers revealed that the mutated cells comprised a large proportion of the microglial pool in the samples examined. While additional studies are required to validate the mechanistic findings, these results suggest that CHIP may have a role in attenuating the risk of AD.

Single-cell multi-omics of mitochondrial DNA disorders reveals dynamics of purifying selection across human immune cells

Pathogenic mutations in mitochondrial DNA (mtDNA) compromise cellular metabolism, contributing to cellular heterogeneity and disease. Diverse mutations are associated with diverse clinical phenotypes, suggesting distinct organ- and cell-type-specific metabolic vulnerabilities. Here we establish a multi-omics approach to quantify deletions in mtDNA alongside cell state features in single cells derived from six patients across the phenotypic spectrum of single large-scale mtDNA deletions (SLSMDs). By profiling 206,663 cells, we reveal the dynamics of pathogenic mtDNA deletion heteroplasmy consistent with purifying selection and distinct metabolic vulnerabilities across T-cell states in vivo and validate these observations in vitro. By extending analyses to hematopoietic and erythroid progenitors, we reveal mtDNA dynamics and cell-type-specific gene regulatory adaptations, demonstrating the context-dependence of perturbing mitochondrial genomic integrity. Collectively, we report pathogenic mtDNA heteroplasmy dynamics of individual blood and immune cells across lineages, demonstrating the power of single-cell multi-omics for revealing fundamental properties of mitochondrial genetics.

Integrative scATAC-seq and scRNA-seq analyses map thymic iNKT cell development and identify Cbfβ for its commitment

Unlike conventional αβT cells, invariant natural killer T (iNKT) cells complete their terminal differentiation to functional iNKT1/2/17 cells in the thymus. However, underlying molecular programs that guide iNKT subset differentiation remain unclear. Here, we profiled the transcriptomes of over 17,000 iNKT cells and the chromatin accessibility states of over 39,000 iNKT cells across four thymic iNKT developmental stages using single-cell RNA sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) to define their developmental trajectories. Our study discovered novel features for iNKT precursors and different iNKT subsets and indicated that iNKT2 and iNKT17 lineage commitment may occur as early as stage 0 (ST0) by two distinct programs, while iNKT1 commitments may occur post ST0. Both iNKT1 and iNKT2 cells exhibit extensive phenotypic and functional heterogeneity, while iNKT17 cells are relatively homogenous. Furthermore, we identified that a novel transcription factor, Cbfβ, was highly expressed in iNKT progenitor commitment checkpoint, which showed a similar expression trajectory with other known transcription factors for iNKT cells development, Zbtb16 and Egr2, and could direct iNKT cells fate and drive their effector phenotype differentiation. Conditional deletion of Cbfβ blocked early iNKT cell development and led to severe impairment of iNKT1/2/17 cell differentiation. Overall, our findings uncovered distinct iNKT developmental programs as well as their cellular heterogeneity, and identified a novel transcription factor Cbfβ as a key regulator for early iNKT cell commitment.

Lipopolysaccharide increases bitter taste sensitivity via epigenetic changes in Tas2r gene clusters

T2R bitter receptors, encoded by Tas2r genes, are not only critical for bitter taste signal transduction but also important for defense against bacteria and parasites. However, little is known about whether and how Tas2r gene expression are regulated. Here, we show that in an inflammation model mimicking bacterial infection using lipopolysaccharide, the expression of many Tas2rs was significantly upregulated and mice displayed markedly increased neural and behavioral responses to bitter compounds. Using single-cell assays for transposase-accessible chromatin with sequencing (scATAC-seq), we found that the chromatin accessibility of Tas2rs was highly celltype specific and lipopolysaccharide increased the accessibility of many Tas2rs. scATAC-seq also revealed substantial chromatin remodeling in immune response genes in taste tissue stem cells, suggesting potential long-lasting effects. Together, our results suggest an epigenetic mechanism connecting inflammation, Tas2r gene regulation, and altered bitter taste, which may explain heightened bitter taste that can occur with infections and cancer treatments.

Single-Cell RNA Sequencing: Technological Progress and Biomedical Application in Cancer Research

Single-cell RNA-seq (scRNA-seq) is a revolutionary technology that allows for the genomic investigation of individual cells in a population, allowing for the discovery of unusual cells associated with cancer and metastasis. ScRNA-seq has been used to discover different types of cancers with poor prognosis and medication resistance such as lung cancer, breast cancer, ovarian cancer, and gastric cancer. Besides, scRNA-seq is a promising method that helps us comprehend the biological features and dynamics of cell development, as well as other disorders. This review gives a concise summary of current scRNA-seq technology. We also explain the main technological steps involved in implementing the technology. We highlight the present applications of scRNA-seq in cancer research, including tumor heterogeneity analysis in lung cancer, breast cancer, and ovarian cancer. In addition, this review elucidates potential applications of scRNA-seq in lineage tracing, personalized medicine, illness prediction, and disease diagnosis, which reveals that scRNA-seq facilitates these events by producing genetic variations on the single-cell level.

A multi-omics atlas of the human retina at single-cell resolution

Cell classes in the human retina are highly heterogeneous with their abundance varying by several orders of magnitude. Here, we generated and integrated a multi-omics single-cell atlas of the adult human retina, including more than 250,000 nuclei for single-nuclei RNA-seq and 137,000 nuclei for single-nuclei ATAC-seq. Cross-species comparison of the retina atlas among human, monkey, mice, and chicken revealed relatively conserved and non-conserved types. Interestingly, the overall cell heterogeneity in primate retina decreases compared with that of rodent and chicken retina. Through integrative analysis, we identified 35,000 distal cis-element-gene pairs, constructed transcription factor (TF)-target regulons for more than 200 TFs, and partitioned the TFs into distinct co-active modules. We also revealed the heterogeneity of the cis-element-gene relationships in different cell types, even from the same class. Taken together, we present a comprehensive single-cell multi-omics atlas of the human retina as a resource that enables systematic molecular characterization at individual cell-type resolution.

Picking up speed: cell cycle regulation during effector CD8+ T cell differentiation

Clonal expansion and development of immunological memory are two hallmarks of adaptive immune responses. Resolving the intricate pathways that regulate cell cycle activity and lead to the generation of diverse effector and memory T cell subsets is essential for improving our understanding of protective T cell immunity. A deeper knowledge of cell cycle regulation in T cells also has translational implications for adoptive cell therapies and vaccinations against infectious diseases. Here, we summarize recent evidence for an early diversification of effector and memory CD8⁺ T cell fates and discuss how this process is coupled to discrete changes in division speed. We further review technical advances in lineage tracing and cell cycle analysis and outline how these techniques have shed new light on the population dynamics of CD8⁺ T cell responses, thereby refining our current understanding of the developmental organization of the memory T cell pool.

Defining and using immune archetypes to classify and treat cancer

Tumours are surrounded by a host immune system that can suppress or promote tumour growth. The tumour microenvironment (TME) has often been framed as a singular entity, suggesting a single type of immune state that is defective and in need of therapeutic intervention. By contrast, the past few years have highlighted a plurality of immune states that can surround tumours. In this Perspective, we suggest that different TMEs have 'archetypal' qualities across all cancers - characteristic and repeating collections of cells and gene-expression profiles at the level of the bulk tumour. We discuss many studies that together support a view that tumours typically draw from a finite number (around 12) of 'dominant' immune archetypes. In considering the likely evolutionary origin and roles of these archetypes, their associated TMEs can be predicted to have specific vulnerabilities that can be leveraged as targets for cancer treatment with expected and addressable adverse effects for patients.

Epigenomic landscape exhibits interferon signaling suppression in the patient of myocarditis after BNT162b2 vaccination

After the outbreak of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic, a novel mRNA vaccine (BNT162b2) was developed at an unprecedented speed. Although most countries have achieved widespread immunity from vaccines and infections, yet people, even who have recovered from SARS-CoV-2 infection, are recommended to receive vaccination due to their effectiveness in lowering the risk of recurrent infection. However, the BNT162b2 vaccine has been reported to increase the risk of myocarditis. To our knowledge, for the first time in this study, we tracked changes in the chromatin dynamics of peripheral blood mononuclear cells (PBMCs) in the patient who underwent myocarditis after BNT162b2 vaccination. A longitudinal study of chromatin accessibility using concurrent analysis of single-cell assays for transposase-accessible chromatin with sequencing and single-cell RNA sequencing showed downregulation of interferon signaling and upregulated RUNX2/3 activity in PBMCs. Considering BNT162b2 vaccination increases the level of interferon-α/γ in serum, our data highlight the immune responses different from the conventional responses to the vaccination, which is possibly the key to understanding the side effects of BNT162b2 vaccination.

scME: a dual-modality factor model for single-cell multiomics embedding

MOTIVATION

Single-cell multiomics technologies are emerging to characterize different molecular features of cells. This gives rise to an issue of combining various kinds of molecular features to dissect cell heterogeneity. Most single-cell multiomics integration methods focus on shared information among modalities while complementary information specific to each modality is often discarded.

RESULTS

To disentangle and combine shared and complementary information across modalities, we develop a dual-modality factor model named scME by using deep factor modeling. Our results demonstrate that scME can generate a better joint representation of multiple modalities than those generated by other single-cell multiomics integration algorithms, which gives a clear elucidation of nuanced differences among cells. We also demonstrate that the joint representation of multiple modalities yielded by scME can provide salient information to improve both single-cell clustering and cell-type classification. Overall, scME will be an efficient method for combining various kinds of molecular features to facilitate the dissection of cell heterogeneity.

AVAILABILITY AND IMPLEMENTATION

The code is public for academic use and available on the GitHub site (https://github.com/bucky527/scME).

Deciphering the tumour immune microenvironment cell by cell

Immune checkpoint inhibitors (ICIs) have rejuvenated therapeutic approaches in oncology. Although responses tend to be durable, response rates vary in many cancer types. Thus, the identification and validation of predictive biomarkers is a key clinical priority, the answer to which is likely to lie in the tumour microenvironment (TME). A wealth of data demonstrates the huge impact of the TME on ICI response and resistance. However, these data also reveal the complexity of the TME composition including the spatiotemporal interactions between different cell types and their dynamic changes in response to ICIs. Here, we briefly review some of the modalities that sculpt the TME, in particular the metabolic milieu, hypoxia and the role of cancer-associated fibroblasts. We then discuss recent approaches to dissect the TME with a focus on single-cell RNA sequencing, spatial transcriptomics and spatial proteomics. We also discuss some of the clinically relevant findings these multi-modal analyses have yielded.

Integrating genetics with single-cell multiomic measurements across disease states identifies mechanisms of beta cell dysfunction in type 2 diabetes

Dysfunctional pancreatic islet beta cells are a hallmark of type 2 diabetes (T2D), but a comprehensive understanding of the underlying mechanisms, including gene dysregulation, is lacking. Here we integrate information from measurements of chromatin accessibility, gene expression and function in single beta cells with genetic association data to nominate disease-causal gene regulatory changes in T2D. Using machine learning on chromatin accessibility data from 34 nondiabetic, pre-T2D and T2D donors, we identify two transcriptionally and functionally distinct beta cell subtypes that undergo an abundance shift during T2D progression. Subtype-defining accessible chromatin is enriched for T2D risk variants, suggesting a causal contribution of subtype identity to T2D. Both beta cell subtypes exhibit activation of a stress-response transcriptional program and functional impairment in T2D, which is probably induced by the T2D-associated metabolic environment. Our findings demonstrate the power of multimodal single-cell measurements combined with machine learning for characterizing mechanisms of complex diseases.

Single-cell technologies in multiple myeloma: new insights into disease pathogenesis and translational implications

Multiple myeloma (MM) is a hematological malignancy characterized by clonal proliferation of plasma cells. Although therapeutic advances have been made to improve clinical outcomes and to prolong patients' survival in the past two decades, MM remains largely incurable. Single-cell sequencing (SCS) is a powerful method to dissect the cellular and molecular landscape at single-cell resolution, instead of providing averaged results. The application of single-cell technologies promises to address outstanding questions in myeloma biology and has revolutionized our understanding of the inter- and intra-tumor heterogeneity, tumor microenvironment, and mechanisms of therapeutic resistance in MM. In this review, we summarize the recently developed SCS methodologies and latest MM research progress achieved by single-cell profiling, including information regarding the cancer and immune cell landscapes, tumor heterogeneities, underlying mechanisms and biomarkers associated with therapeutic response and resistance. We also discuss future directions of applying transformative SCS approaches with contribution to clinical translation.

Integrated analysis of single-cell chromatin state and transcriptome identified common vulnerability despite glioblastoma heterogeneity

In 2021, the World Health Organization reclassified glioblastoma, the most common form of adult brain cancer, into isocitrate dehydrogenase (IDH)-wild-type glioblastomas and grade IV IDH mutant (G4 IDHm) astrocytomas. For both tumor types, intratumoral heterogeneity is a key contributor to therapeutic failure. To better define this heterogeneity, genome-wide chromatin accessibility and transcription profiles of clinical samples of glioblastomas and G4 IDHm astrocytomas were analyzed at single-cell resolution. These profiles afforded resolution of intratumoral genetic heterogeneity, including delineation of cell-to-cell variations in distinct cell states, focal gene amplifications, as well as extrachromosomal circular DNAs. Despite differences in IDH mutation status and significant intratumoral heterogeneity, the profiled tumor cells shared a common chromatin structure defined by open regions enriched for nuclear factor 1 transcription factors (NFIA and NFIB). Silencing of NFIA or NFIB suppressed in vitro and in vivo growths of patient-derived glioblastomas and G4 IDHm astrocytoma models. These findings suggest that despite distinct genotypes and cell states, glioblastoma/G4 astrocytoma cells share dependency on core transcriptional programs, yielding an attractive platform for addressing therapeutic challenges associated with intratumoral heterogeneity.

Epigenetic plasticity cooperates with cell-cell interactions to direct pancreatic tumorigenesis

The response to tumor-initiating inflammatory and genetic insults can vary among morphologically indistinguishable cells, suggesting as yet uncharacterized roles for epigenetic plasticity during early neoplasia. To investigate the origins and impact of such plasticity, we performed single-cell analyses on normal, inflamed, premalignant, and malignant tissues in autochthonous models of pancreatic cancer. We reproducibly identified heterogeneous cell states that are primed for diverse, late-emerging neoplastic fates and linked these to chromatin remodeling at cell-cell communication loci. Using an inference approach, we revealed signaling gene modules and tissue-level cross-talk, including a neoplasia-driving feedback loop between discrete epithelial and immune cell populations that was functionally validated in mice. Our results uncover a neoplasia-specific tissue-remodeling program that may be exploited for pancreatic cancer interception.