A Single-Cell Atlas of In Vivo Mammalian Chromatin Accessibility

scCorrect: Cross-modality label transfer from scRNA-seq to scATAC-seq using domain adaptation

Cell type annotation in single-cell chromatin accessibility sequencing (scATAC-seq) is crucial for enabling researchers to identify subpopulations of cells associated with specific diseases, elucidate gene regulatory networks, and discover markers indicative of disease states. The prevailing approach for cell type annotation in single-cell research involves transferring well-delineated cell types from single-cell RNA sequencing (scRNA-seq) data to scATAC-seq data using a label propagation algorithm. However, the inherent modal discrepancies (i.e.biological interpretation) between scRNA-seq and scATAC-seq data, coupled with the intrinsic sparsity and high dimensionality of scATAC-seq data, pose significant challenges to the efficacy of this strategy. To address these challenges, we introduce a novel neural network framework, scCorrect, which operates in two distinct phases. In the first phase, scCorrect aligns the scRNA-seq and scATAC-seq datasets, generating initial annotation results. The second phase involves training a corrective network specifically designed to amend any erroneous annotations produced during the first phase. Empirical tests across multiple datasets have demonstrated that scCorrect consistently achieves superior recognition accuracy, underscoring its significant potential to enhance disease-related research in humans.

Preeclampsia risk prediction from prenatal cell-free DNA screening

Preeclampsia is characterized by placental dysfunction and results in significant morbidity, but reliable early prediction remains challenging. We investigated whether clinically obtained prenatal cell-free DNA (cfDNA) screening (PDNAS) using whole-genome sequencing (WGS) data can be leveraged to predict preeclampsia risk early in pregnancy (≤16 weeks). Using 1,854 routinely collected clinical PDNAS samples (median, 12.1 weeks) with low-coverage (0.5×) WGS data, we developed a framework to quantify maternal and fetal tissue signatures using nucleosome accessibility, revealing early placental and endothelial dysfunction. These signatures informed a prediction model for preeclampsia risk, which achieved a validation performance of 0.85 area under the receiver operating characteristic curve (AUC) (81% sensitivity at 80% specificity) for preterm phenotypes several months prior to disease onset in a separate cohort of 831 consecutively collected samples, and subsequently confirmed in an external cohort of 141 samples (AUC 0.84, 79% sensitivity). We demonstrate that assessment of cfDNA nucleosome accessibility from early-pregnancy cfDNA sequence data enables the detection of early placental and endothelial-tissue aberrations and may aid in the determination of preeclampsia risk.

SCITUNA: single-cell data integration tool using network alignment

BACKGROUND

As single-cell genomics experiments increase in complexity and scale, the need to integrate multiple datasets has grown. Such integration enhances cellular feature identification by leveraging larger data volumes. However, batch effects-technical variations arising from differences in labs, times, or protocols-pose a significant challenge. Despite numerous proposed batch correction methods, many still have limitations, such as outputting only dimension-reduced data, relying on computationally intensive models, or resulting in overcorrection for batches with diverse cell type composition.

RESULTS

We introduce a novel method for batch effect correction named SCITUNA, a Single-Cell data Integration Tool Using Network Alignment. We perform evaluations on 39 individual batches from four real datasets and a simulated dataset, which include both scRNA-seq and scATAC-seq datasets, spanning multiple organisms and tissues. A thorough comparison of existing batch correction methods using 13 metrics reveals that SCITUNA outperforms current approaches and is successful at preserving biological signals present in the original data. In particular, SCITUNA shows a better performance than the current methods in all the comparisons except for the multiple batch integration of the lung dataset where the difference is 0.004.

CONCLUSION

SCITUNA effectively removes batch effects while retaining the biological signals present in the data. Our extensive experiments reveal that SCITUNA will be a valuable tool for diverse integration tasks.

scMDCL: A Deep Collaborative Contrastive Learning Framework for Matched Single-Cell Multiomics Data Clustering

Single-cell multiomics clustering integrates multiple omics data to analyze cellular heterogeneity and is crucial for uncovering complex biological processes and disease mechanisms. However, existing matched single-cell multiomics clustering methods often neglect the full utilization of intercellular relationships and the interactions and synergy between features from different omics, leading to suboptimal clustering performance. In this paper, we propose a deep collaborative contrastive learning framework for matched single-cell multiomics data clustering, named scMDCL. This framework fully leverages intercell relationships while enhancing feature interactions among identical cells across different omics data, thereby facilitating efficient clustering of multiomics data. Specifically, to fully utilize the topological information between cells, a graph autoencoder and a feature information enhancement module are designed for different omics, enabling the extraction and augmentation of cell features. Additionally, contrastive learning techniques are employed to strengthen the interactions among the different omics features of the same cell. Ultimately, multiomics deep collaborative clustering modules are utilized to achieve single-cell multiomics clustering. Extensive experiments conducted on nine publicly available single-cell multiomics datasets demonstrate the superior performance of the proposed framework in integrating multiomics data for clustering tasks.

Identifying cross-lineage dependencies of cell-type-specific regulators in mouse gastruloids

Correct gene expression levels are crucial for normal development. Advances in genomics enable the inference of gene regulatory programs active during development but cannot capture the complex multicellular interactions occurring during mammalian embryogenesis in utero. In vitro models of mammalian development, like gastruloids, can overcome this limitation. Using time-resolved single-cell chromatin accessibility analysis, we delineated the regulatory profile during mouse gastruloid development, identifying critical drivers of developmental transitions. Gastruloids develop from bipotent progenitor cells driven by the transcription factors (TFs) OCT4, SOX2, and TBXT, differentiating into the mesoderm (characterized by the mesogenin 1 [MSGN1]) and spinal cord (characterized by CDX2). ΔCDX gastruloids fail to form spinal cord, while Msgn1 ablation inhibits paraxial mesoderm and spinal cord development. Chimeric gastruloids with ΔMSGN1 and wild-type cells formed both tissues, indicating that inter-tissue communication is necessary for spinal cord formation. Our work has important implications for studying inter-tissue communication and gene regulatory programs in development.

Inferring the Selective History of CNVs Using a Maximum Likelihood Model

Copy number variations (CNVs)-structural variations generated by deletion and/or duplication that result in a change in DNA dosage-are prevalent in nature. CNVs can drastically affect the phenotype of an organism and have been shown to be both involved in genetic disorders and be used as raw material in adaptive evolution. Unlike single-nucleotide variations, the often large and varied effects of CNVs on phenotype hinders our ability to infer their selective advantage based on the population genetics data. Here, we present a likelihood-based approach, dubbed PoMoCNV (POlymorphism-aware phylogenetic MOdel for CNVs), that estimates the evolutionary parameters such as mutation rates among different copy numbers and relative fitness loss per copy deletion at a genomic locus based on population genetics data. As a case study, we analyze the genomics data of 40 strains of Caenorhabditis elegans, representing four different populations. We take advantage of the data on chromatin accessibility to interpret the mutation rate and fitness of copy numbers, as inferred by PoMoCNV, specifically in open or closed chromatin loci. We further test the reliability of PoMoCNV by estimating the evolutionary parameters of CNVs for mutation-accumulation experiments in C. elegans with varying levels of genetic drift.

Multiomic single-cell profiling identifies critical regulators of postnatal brain

Human brain development spans from embryogenesis to adulthood, with dynamic gene expression controlled by cell-type-specific cis-regulatory element activity and three-dimensional genome organization. To advance our understanding of postnatal brain development, we simultaneously profiled gene expression and chromatin accessibility in 101,924 single nuclei from four brain regions across ten donors, covering five key postnatal stages from infancy to late adulthood. Using this dataset and chromosome conformation capture data, we constructed enhancer-based gene regulatory networks to identify cell-type-specific regulators of brain development and interpret genome-wide association study loci for ten main brain disorders. Our analysis connected 2,318 cell-specific loci to 1,149 unique genes, representing 41% of loci linked to the investigated traits, and highlighted 55 genes influencing several disease phenotypes. Pseudotime analysis revealed distinct stages of postnatal oligodendrogenesis and their regulatory programs. These findings provide a comprehensive dataset of cell-type-specific gene regulation at critical timepoints in postnatal brain development.

Reducing batch effects in single cell chromatin accessibility measurements by pooled transposition with MULTI-ATAC

Large-scale scATAC-seq experiments are challenging because of their costs, lengthy protocols, and confounding batch effects. Several sample multiplexing technologies aim to address these challenges, but do not remove batch effects introduced when performing transposition reactions in parallel. We demonstrate that sample-to-sample variability in nuclei-to-Tn5 ratios is a major cause of batch effects and develop MULTI-ATAC, a multiplexing method that pools samples prior to transposition, as a solution. MULTI-ATAC provides high accuracy in sample classification and doublet detection while eliminating batch effects associated with variable nucleus-to-Tn5 ratio. We illustrate the power of MULTI-ATAC by performing a 96-plex multiomic drug assay targeting epigenetic remodelers in a model of primary immune cell activation, uncovering tens of thousands of drug-responsive chromatin regions, cell-type specific effects, and potent differences between matched inhibitors and degraders. MULTI-ATAC therefore enables batch-free and scalable scATAC-seq workflows, providing deeper insights into complex biological processes and potential therapeutic targets.

EpiFoundation: A Foundation Model for Single-Cell ATAC-seq via Peak-to-Gene Alignment

Foundation models exhibit strong capabilities for downstream tasks by learning generalized representations through self-supervised pre-training on large datasets. While several foundation models have been developed for single-cell RNA-seq (scRNA-seq) data, there is still a lack of models specifically tailored for single-cell ATAC-seq (scATAC-seq), which measures epigenetic information in individual cells. The principal challenge in developing such a model lies in the vast number of scATAC peaks and the significant sparsity of the data, which complicates the formulation of peak-to-peak correlations. To address this challenge, we introduce EpiFoundation, a foundation model for learning cell representations from the high-dimensional and sparse space of peaks. EpiFoundation relies on an innovative cross-modality pre-training procedure with two key technical innovations. First, EpiFoundation exclusively processes the non-zero peak set, thereby enhancing the density of cell-specific information within the input data. Second, EpiFoundation utilizes dense gene expression information to supervise the pre-training process, aligning peak-to-gene correlations. EpiFoundation can handle various types of downstream tasks, including cell-type annotation, batch correction, and gene expression prediction. To train and validate EpiFoundation, we curated MiniAtlas, a dataset of 100,000+ single cells with paired scRNA-seq and scATAC-seq data, along with diverse test sets spanning various tissues and cell types for robust evaluation. EpiFoundation demonstrates state-of-the-art performance across multiple tissues and diverse downstream tasks.

Principled PCA separates signal from noise in omics count data

Principal component analysis (PCA) is indispensable for processing high-throughput omics datasets, as it can extract meaningful biological variability while minimizing the influence of noise. However, the suitability of PCA is contingent on appropriate normalization and transformation of count data, and accurate selection of the number of principal components; improper choices can result in the loss of biological information or corruption of the signal due to excessive noise. Typical approaches to these challenges rely on heuristics that lack theoretical foundations. In this work, we present Biwhitened PCA (BiPCA), a theoretically grounded framework for rank estimation and data denoising across a wide range of omics modalities. BiPCA overcomes a fundamental difficulty with handling count noise in omics data by adaptively rescaling the rows and columns - a rigorous procedure that standardizes the noise variances across both dimensions. Through simulations and analysis of over 100 datasets spanning seven omics modalities, we demonstrate that BiPCA reliably recovers the data rank and enhances the biological interpretability of count data. In particular, BiPCA enhances marker gene expression, preserves cell neighborhoods, and mitigates batch effects. Our results establish BiPCA as a robust and versatile framework for high-throughput count data analysis.

Cross-species comparative single-cell transcriptomics highlights the molecular evolution and genetic basis of male infertility

In male animals, spermatogonia in testes differentiate into sperm, one of the most diverse cell types across species. Despite the evolutionary retention of key genes essential for spermatogenesis, the extent of their conservation remains unclear. To explore the genetic basis of spermatogenesis under strong selective pressure, we compare single-cell RNA sequencing (scRNA-seq) datasets from the testes of humans, mice, and fruit flies. Our analysis identifies conserved genes involved in key molecular programs, such as post-transcriptional regulation, meiosis, and energy metabolism. We perform gene knockout experiments of 20 candidate genes, three of which, when mutated in fruit flies, result in reduced male fertility, emphasizing the conservation of sperm centriole and steroid lipid processes across mammals and Drosophila. Additionally, deep-learning analysis uncovers potential transcriptional mechanisms driving gene-expression evolution. These findings establish a core genetic foundation for spermatogenesis, offering insights into sperm-phenotype evolution and the underlying mechanisms of male infertility.

Modular organization of enhancer network provides transcriptional robustness in mammalian development

Enhancer clusters, pivotal in mammalian development and diseases, can organize as enhancer networks to control cell identity and disease genes; however, the underlying mechanism remains largely unexplored. Here, we introduce eNet 2.0, a comprehensive tool for enhancer networks analysis during development and diseases based on single-cell chromatin accessibility data. eNet 2.0 extends our previous work eNet 1.0 by adding network topology, comparison and dynamics analyses to its network construction function. We reveal modularly organized enhancer networks, where inter-module interactions synergistically affect gene expression. Moreover, network alterations correlate with abnormal and dynamic gene expression in disease and development. eNet 2.0 is robust across diverse datasets. To facilitate application, we introduce eNetDB (https://enetdb.huanglabxmu.com), an enhancer network database leveraging extensive scATAC-seq (single-cell assay for transposase-accessible chromatin sequencing) datasets from human and mouse tissues. Together, our work provides a powerful computational tool and reveals that modularly organized enhancer networks contribute to gene expression robustness in mammalian development and diseases.

Leveraging Single-Cell Multi-Omics to Decode Tumor Microenvironment Diversity and Therapeutic Resistance

Recent developments in single-cell multi-omics technologies have provided the ability to identify diverse cell types and decipher key components of the tumor microenvironment (TME), leading to important advancements toward a much deeper understanding of how tumor microenvironment heterogeneity contributes to cancer progression and therapeutic resistance. These technologies are able to integrate data from molecular genomic, transcriptomic, proteomics, and metabolomics studies of cells at a single-cell resolution scale that give rise to the full cellular and molecular complexity in the TME. Understanding the complex and sometimes reciprocal relationships among cancer cells, CAFs, immune cells, and ECs has led to novel insights into their immense heterogeneity in functions, which can have important consequences on tumor behavior. In-depth studies have uncovered immune evasion mechanisms, including the exhaustion of T cells and metabolic reprogramming in response to hypoxia from cancer cells. Single-cell multi-omics also revealed resistance mechanisms, such as stromal cell-secreted factors and physical barriers in the extracellular matrix. Future studies examining specific metabolic pathways and targeting approaches to reduce the heterogeneity in the TME will likely lead to better outcomes with immunotherapies, drug delivery, etc., for cancer treatments. Future studies will incorporate multi-omics data, spatial relationships in tumor micro-environments, and their translation into personalized cancer therapies. This review emphasizes how single-cell multi-omics can provide insights into the cellular and molecular heterogeneity of the TME, revealing immune evasion mechanisms, metabolic reprogramming, and stromal cell influences. These insights aim to guide the development of personalized and targeted cancer therapies, highlighting the role of TME diversity in shaping tumor behavior and treatment outcomes.

Atlas-scale single-cell DNA methylation profiling with sciMETv3

Single-cell methods to assess DNA methylation have not achieved the same level of cell throughput per experiment compared to other modalities, with large-scale datasets requiring extensive automation, time, and other resources. Here, we describe sciMETv3, a combinatorial indexing-based technique that enables atlas-scale libraries to be produced in a single experiment. To reduce the sequencing burden, we demonstrate the compatibility of sciMETv3 with capture techniques to enrich regulatory regions, as well as the ability to leverage enzymatic conversion, which can yield higher library diversity. We showcase the throughput of sciMETv3 by producing a >140,000 cell library from human middle frontal gyrus split across four multiplexed individuals using both Illumina and Ultima sequencing instrumentation. Finally, we introduce sciMET+ATAC to enable high-throughput exploration of the interplay between chromatin accessibility and DNA methylation within the same cell.

Depth-corrected multi-factor dissection of chromatin accessibility for scATAC-seq data with PACS

Single cell ATAC-seq (scATAC-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 scATAC-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 allows complex hypothesis testing of accessibility-modulating factors while accounting for sparse and incomplete data. For differential accessibility analysis, PACS controls the false positive rate and achieves a 17% to 122% higher power on average 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 datasets from various tissues and show its ability to reveal previously undiscovered insights in scATAC-seq data.

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.

Integration of unpaired single cell omics data by deep transfer graph convolutional network

The rapid advance of large-scale atlas-level single cell RNA sequences and single-cell chromatin accessibility data provide extraordinary avenues to broad and deep insight into complex biological mechanism. Leveraging the datasets and transfering labels from scRNA-seq to scATAC-seq will empower the exploration of single-cell omics data. However, the current label transfer methods have limited performance, largely due to the lower capable of preserving fine-grained cell populations and intrinsic or extrinsic heterogeneity between datasets. Here, we present a robust deep transfer model based graph convolutional network, scTGCN, which achieves versatile performance in preserving biological variation, while achieving integration hundreds of thousands cells in minutes with low memory consumption. We show that scTGCN is powerful to the integration of mouse atlas data and multimodal data generated from APSA-seq and CITE-seq. Thus, scTGCN shows high label transfer accuracy and effectively knowledge transfer across different modalities.

Single-Cell Analysis Provides New Insights into the Roles of Tertiary Lymphoid Structures and Immune Cell Infiltration in Kidney Injury and Chronic Kidney Disease

Chronic kidney disease (CKD) is a global health concern with high morbidity and mortality. Acute kidney injury (AKI) is a pivotal risk factor for the progression of CKD, and the rate of AKI-to-CKD progression increases with aging. Intrarenal inflammation is a fundamental mechanism underlying AKI-to-CKD progression. Tertiary lymphoid structures (TLSs), ectopic lymphoid aggregates formed in nonlymphoid organs, develop in aged injured kidneys, but not in young kidneys, with prolonged inflammation and maladaptive repair, which potentially exacerbates AKI-to-CKD progression in aged individuals. Dysregulated immune responses are involved in the pathogenesis of various kidney diseases, such as IgA nephropathy, lupus nephritis, and diabetic kidney diseases, thereby deteriorating kidney function. TLSs also develop in several kidney diseases, including transplanted kidneys and renal cell carcinoma. However, the precise immunologic mechanisms driving AKI-to-CKD progression and development of these kidney diseases remain unclear, which hinders the development of novel therapeutic approaches. This review aims to describe recent findings from single-cell analysis of cellular heterogeneity and complex interactions among immune and renal parenchymal cells, which potentially contribute to the pathogenesis of AKI-to-CKD progression and other kidney diseases, highlighting the mechanisms of formation and pathogenic roles of TLSs in aged injured kidneys.

Spatiotemporal modeling of molecular holograms

Quantifying spatiotemporal dynamics during embryogenesis is crucial for understanding congenital diseases. We developed Spateo (https://github.com/aristoteleo/spateo-release), a 3D spatiotemporal modeling framework, and applied it to a 3D mouse embryogenesis atlas at E9.5 and E11.5, capturing eight million cells. Spateo enables scalable, partial, non-rigid alignment, multi-slice refinement, and mesh correction to create molecular holograms of whole embryos. It introduces digitization methods to uncover multi-level biology from subcellular to whole organ, identifying expression gradients along orthogonal axes of emergent 3D structures, e.g., secondary organizers such as midbrain-hindbrain boundary (MHB). Spateo further jointly models intercellular and intracellular interaction to dissect signaling landscapes in 3D structures, including the zona limitans intrathalamica (ZLI). Lastly, Spateo introduces "morphometric vector fields" of cell migration and integrates spatial differential geometry to unveil molecular programs underlying asymmetrical murine heart organogenesis and others, bridging macroscopic changes with molecular dynamics. Thus, Spateo enables the study of organ ecology at a molecular level in 3D space over time.

Advances and applications in single-cell and spatial genomics

The applications of single-cell and spatial technologies in recent times have revolutionized the present understanding of cellular states and the cellular heterogeneity inherent in complex biological systems. These advancements offer unprecedented resolution in the examination of the functional genomics of individual cells and their spatial context within tissues. In this review, we have comprehensively discussed the historical development and recent progress in the field of single-cell and spatial genomics. We have reviewed the breakthroughs in single-cell multi-omics technologies, spatial genomics methods, and the computational strategies employed toward the analyses of single-cell atlas data. Furthermore, we have highlighted the advances made in constructing cellular atlases and their clinical applications, particularly in the context of disease. Finally, we have discussed the emerging trends, challenges, and opportunities in this rapidly evolving field.

Should Artificial Intelligence Play a Durable Role in Biomedical Research and Practice?

During the last decade, artificial intelligence (AI) was applied to nearly all domains of human activity, including scientific research. It is thus warranted to ask whether AI thinking should be durably involved in biomedical research. This problem was addressed by examining three complementary questions (i) What are the major barriers currently met by biomedical investigators? It is suggested that during the last 2 decades there was a shift towards a growing need to elucidate complex systems, and that this was not sufficiently fulfilled by previously successful methods such as theoretical modeling or computer simulation (ii) What is the potential of AI to meet the aforementioned need? it is suggested that recent AI methods are well-suited to perform classification and prediction tasks on multivariate systems, and possibly help in data interpretation, provided their efficiency is properly validated. (iii) Recent representative results obtained with machine learning suggest that AI efficiency may be comparable to that displayed by human operators. It is concluded that AI should durably play an important role in biomedical practice. Also, as already suggested in other scientific domains such as physics, combining AI with conventional methods might generate further progress and new applications, involving heuristic and data interpretation.

Going beyond cell clustering and feature aggregation: Is there single cell level information in single-cell ATAC-seq data?

Single-cell Assay for Transposase Accessible Chromatin with sequencing (scATAC-seq) has become a widely used method for investigating chromatin accessibility at single-cell resolution. However, the resulting data is highly sparse with most data entries being zeros. As such, currently available computational methods for scATAC-seq feature a range of transformation procedures to extract meaningful information from the sparse data. Most notably, these transformations can be categorized into: 1) feature aggregation with known biological associations, 2) pseudo-bulking cells of similar biology, and 3) binarisation of count data. These strategies beg the question of whether or not scATAC-seq data actually has usable single-cell and single-region information as intended from the assay. If we can go beyond aggregated features and pooled cells, it opens up the possibility of more complex statistical tasks that require that degree of granularity. To reach the finest possible resolution of single-cell, single-region information there are inevitably many computational challenges to overcome. Here, we review the major data analysis challenges lying between raw data readout and biological discovery, and discuss the limitations of current data analysis approaches. Lastly, we conclude that chromatin accessibility profiling at true single-cell resolution is not yet achieved with current technology, but that it may be achieved with promising developments in optimising the efficiency of scATAC-seq assays.

Long Noncoding RNA TRIBAL Links the 8q24.13 Locus to Hepatic Lipid Metabolism and Coronary Artery Disease

BACKGROUND

Genome-wide association studies identified a 20-Kb region of chromosome 8 (8q24.13) associated with plasma lipids, hepatic steatosis, and risk for coronary artery disease. The region is proximal to TRIB1, and given its well-established role in lipid regulation in animal models, TRIB1 has been proposed to mediate the contribution of the 8q24.13 locus to these traits. This region overlaps a gene encoding the primate-specific long noncoding RNA transcript TRIBAL/TRIB1AL (TRIB1-associated locus), but the contribution of TRIBAL to coronary artery disease risk remains untested.

METHODS

Using recently available expression quantitative trait loci data and hepatocyte models, we further investigated this locus by Mendelian randomization analysis. Following antisense oligonucleotide targeting of TRIBAL, transcription array, quantitative reverse transcription polymerase chain reaction, and enrichment analyses were performed and effects on apoB and triglyceride secretion were determined.

RESULTS

Mendelian randomization analysis supports a causal relationship between genetically determined hepatic TRIBAL expression and markers of hepatic steatosis and coronary artery disease risk. By contrast, expression data sets did not support expression quantitative trait loci relationships between coronary artery disease-associated variants and TRIB1. TRIBAL suppression reduced the expression of key regulators of triglyceride metabolism and bile acid synthesis. Enrichment analyses identified patterns consistent with impaired metabolic functions, including reduced triglyceride and cholesterol handling ability. Furthermore, TRIBAL suppression was associated with reduced hepatocyte secretion of triglycerides.

CONCLUSIONS

This work identifies TRIBAL as a gene bridging the genotype-phenotype relationship at the 8q24.13 locus with effects on genes regulating hepatocyte lipid metabolism and triglyceride secretion.

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

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

Recruitment of homodimeric proneural factors by conserved CAT-CAT E-boxes drives major epigenetic reconfiguration in cortical neurogenesis

Proneural factors of the basic helix-loop-helix family coordinate neurogenesis and neurodifferentiation. Among them, NEUROG2 and NEUROD2 subsequently act to specify neurons of the glutamatergic lineage. Disruption of these factors, their target genes and binding DNA motifs has been linked to various neuropsychiatric disorders. Proneural factors bind to specific DNA motifs called E-boxes (hexanucleotides of the form CANNTG, composed of two CAN half sites on opposed strands). While corticogenesis heavily relies on E-box activity, the collaboration of proneural factors on different E-box types and their chromatin remodeling mechanisms remain largely unknown. Here, we conducted a comprehensive analysis using chromatin immunoprecipitation followed by sequencing (ChIP-seq) data for NEUROG2 and NEUROD2, along with time-matched single-cell RNA-seq, ATAC-seq and DNA methylation data from the developing mouse cortex. Our findings show that these factors are highly enriched in transiently active genomic regions during intermediate stages of neuronal differentiation. Although they primarily bind CAG-containing E-boxes, their binding in dynamic regions is notably enriched in CAT-CAT E-boxes (i.e. CATATG, denoted as 5'3' half sites for dimers), which undergo significant DNA demethylation and exhibit the highest levels of evolutionary constraint. Aided by HT-SELEX data reanalysis, structural modeling and DNA footprinting, we propose that these proneural factors exert maximal chromatin remodeling influence during intermediate stages of neurogenesis by binding as homodimers to CAT-CAT motifs. This study provides an in-depth integrative analysis of the dynamic regulation of E-boxes during neuronal development, enhancing our understanding of the mechanisms underlying the binding specificity of critical proneural factors.

Epigenetics-targeted drugs: current paradigms and future challenges

Epigenetics governs a chromatin state regulatory system through five key mechanisms: DNA modification, histone modification, RNA modification, chromatin remodeling, and non-coding RNA regulation. These mechanisms and their associated enzymes convey genetic information independently of DNA base sequences, playing essential roles in organismal development and homeostasis. Conversely, disruptions in epigenetic landscapes critically influence the pathogenesis of various human diseases. This understanding has laid a robust theoretical groundwork for developing drugs that target epigenetics-modifying enzymes in pathological conditions. Over the past two decades, a growing array of small molecule drugs targeting epigenetic enzymes such as DNA methyltransferase, histone deacetylase, isocitrate dehydrogenase, and enhancer of zeste homolog 2, have been thoroughly investigated and implemented as therapeutic options, particularly in oncology. Additionally, numerous epigenetics-targeted drugs are undergoing clinical trials, offering promising prospects for clinical benefits. This review delineates the roles of epigenetics in physiological and pathological contexts and underscores pioneering studies on the discovery and clinical implementation of epigenetics-targeted drugs. These include inhibitors, agonists, degraders, and multitarget agents, aiming to identify practical challenges and promising avenues for future research. Ultimately, this review aims to deepen the understanding of epigenetics-oriented therapeutic strategies and their further application in clinical settings.

DECA: harnessing interpretable transformer model for cellular deconvolution of chromatin accessibility profile

The assay for transposase-accessible chromatin with sequencing (ATAC-seq) identifies chromatin accessibility across the genome, crucial for gene expression regulating. However, bulk ATAC-seq obscures cellular heterogeneity, while single-cell ATAC-seq suffers from issues such as sparsity and costliness. To this end, we introduce DECA, a sophisticated deep learning model based on vision transformer to deconvolve cell type information from bulk chromatin accessibility profiles, utilizing single-cell ATAC-seq datasets as reference for enhanced precision and resolution. Notably, patch attention generated by DECA's multi-head attention mechanism aligns with chromatin interactions detected by Hi-C. Additionally, DECA predicted lineage-specific cell composition changes due to genetic perturbation. The chromatin accessibility signatures predicted by DECA are enriched with cell-type specific genetic variations. Ultimately, we applied DECA on pan-cancer ATAC-seq datasets and demonstrated its capability to deconvolve cell type proportions with clinical significance. Taken together, DECA deconvolves cellular proportions and predicts their chromatin accessibility profiles from bulk chromatin accessibility data, which enable exploring the gene regulatory programs in development and diseases.

Exploring the Unknown: How Can We Improve Single-cell RNAseq Cell Type Annotations in Non-model Organisms?

Single-cell RNA sequencing (scRNAseq) is a powerful tool to describe cell types in multicellular organisms across the animal kingdom. In standard scRNAseq analysis pipelines, clusters of cells with similar transcriptional signatures are given cell type labels based on marker genes that infer specialized known characteristics. Since these analyses are designed for model organisms, such as humans and mice, problems arise when attempting to label cell types of distantly related, non-model species that have unique or divergent cell types. Consequently, this leads to limited discovery of novel species-specific cell types and potential mis-annotation of cell types in non-model species while using scRNAseq. To address this problem, we discuss recently published approaches that help annotate scRNAseq clusters for any non-model organism. We first suggest that annotating with an evolutionary context of cell lineages will aid in the discovery of novel cell types and provide a marker-free approach to compare cell types across distantly related species. Secondly, machine learning has greatly improved bioinformatic analyses, so we highlight some open-source programs that use reference-free approaches to annotate cell clusters. Lastly, we propose the use of unannotated genes as potential cell markers for non-model organisms, as many do not have fully annotated genomes and these data are often disregarded. Improving single-cell annotations will aid the discovery of novel cell types and enhance our understanding of non-model organisms at a cellular level. By unifying approaches to annotate cell types in non-model organisms, we can increase the confidence of cell annotation label transfer and the flexibility to discover novel cell types.

Single-cell multi-modal integrative analyses highlight functional dynamic gene regulatory networks directing human cardiac development

Illuminating the precise stepwise genetic programs directing cardiac development provides insights into the mechanisms of congenital heart disease and strategies for cardiac regenerative therapies. Here, we integrate in vitro and in vivo human single-cell multi-omic studies with high-throughput functional genomic screening to reveal dynamic, cardiac-specific gene regulatory networks (GRNs) and transcriptional regulators during human cardiomyocyte development. Interrogating developmental trajectories reconstructed from single-cell data unexpectedly reveal divergent cardiomyocyte lineages with distinct gene programs based on developmental signaling pathways. High-throughput functional genomic screens identify key transcription factors from inferred GRNs that are functionally relevant for cardiomyocyte lineages derived from each pathway. Notably, we discover a critical heat shock transcription factor 1 (HSF1)-mediated cardiometabolic GRN controlling cardiac mitochondrial/metabolic function and cell survival, also observed in fetal human cardiomyocytes. Overall, these multi-modal genomic studies enable the systematic discovery and validation of coordinated GRNs and transcriptional regulators controlling the development of distinct human cardiomyocyte populations.

Human vascularized macrophage-islet organoids to model immune-mediated pancreatic β cell pyroptosis upon viral infection

There is a paucity of human models to study immune-mediated host damage. Here, we utilized the GeoMx spatial multi-omics platform to analyze immune cell changes in COVID-19 pancreatic autopsy samples, revealing an accumulation of proinflammatory macrophages. Single-cell RNA sequencing (scRNA-seq) analysis of human islets exposed to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) or coxsackievirus B4 (CVB4) viruses identified activation of proinflammatory macrophages and β cell pyroptosis. To distinguish viral versus proinflammatory-macrophage-mediated β cell pyroptosis, we developed human pluripotent stem cell (hPSC)-derived vascularized macrophage-islet (VMI) organoids. VMI organoids exhibited enhanced marker expression and function in both β cells and endothelial cells compared with separately cultured cells. Notably, proinflammatory macrophages within VMI organoids induced β cell pyroptosis. Mechanistic investigations highlighted TNFSF12-TNFRSF12A involvement in proinflammatory-macrophage-mediated β cell pyroptosis. This study established hPSC-derived VMI organoids as a valuable tool for studying immune-cell-mediated host damage and uncovered the mechanism of β cell damage during viral exposure.

Single-cell sequencing: Current applications in various tuberculosis specimen types

Tuberculosis (TB) is a chronic disease caused by Mycobacterium tuberculosis (M.tb) and responsible for millions of deaths worldwide each year. It has a complex pathogenesis that primarily affects the lungs but can also impact systemic organs. In recent years, single-cell sequencing technology has been utilized to characterize the composition and proportion of immune cell subpopulations associated with the pathogenesis of TB disease since it has a high resolution that surpasses conventional techniques. This paper reviews the current use of single-cell sequencing technologies in TB research and their application in analysing specimens from various sources of TB, primarily peripheral blood and lung specimens. The focus is on how these technologies can reveal dynamic changes in immune cell subpopulations, genes and proteins during disease progression after M.tb infection. Based on the current findings, single-cell sequencing has significant potential clinical value in the field of TB research. Next, we will focus on the real-world applications of the potential targets identified through single-cell sequencing for diagnostics, therapeutics and the development of effective vaccines.

Characterization of Single-Cell Cis-regulatory Elements Informs Implications for Cell Differentiation

Cis-regulatory elements govern the specific patterns and dynamics of gene expression in cells during development, which are the fundamental mechanisms behind cell differentiation. However, the genomic characteristics of single-cell cis-regulatory elements closely linked to cell differentiation during development remain unclear. To explore this, we systematically analyzed ∼250,000 putative single-cell cis-regulatory elements obtained from snATAC-seq analysis of the developing mouse cerebellum. We found that over 80% of these single-cell cis-regulatory elements show pleiotropic effects, being active in 2 or more cell types. The pleiotropic degrees of proximal and distal single-cell cis-regulatory elements are positively correlated with the density and diversity of transcription factor binding motifs and GC content. There is a negative correlation between the pleiotropic degrees of single-cell cis-regulatory elements and their distances to the nearest transcription start sites, and proximal single-cell cis-regulatory elements display higher relevance strengths than distal ones. Furthermore, both proximal and distal single-cell cis-regulatory elements related to cell differentiation exhibit enhanced sequence-level evolutionary conservation, increased density and diversity of transcription factor binding motifs, elevated GC content, and greater distances from their nearest genes. Together, our findings reveal the general genomic characteristics of putative single-cell cis-regulatory elements and provide insights into the genomic and evolutionary mechanisms by which single-cell cis-regulatory elements regulate cell differentiation during development.

scPerb: Predict single-cell perturbation via style transfer-based variational autoencoder

INTRODUCTION

Traditional methods for obtaining cellular responses after perturbation are usually labor-intensive and costly, especially when working with multiple different experimental conditions. Therefore, accurate prediction of cellular responses to perturbations is of great importance in computational biology. Existing methodologies, such as graph-based approaches, vector arithmetic, and neural networks, either mix perturbation-related variances with cell-type-specific patterns or implicitly distinguish them within black-box models.

OBJECTIVES

This study aims to introduce and demonstrate a novel framework, scPerb, which explicitly extracts perturbation-related variances and transfers them from unperturbed to perturbed cells to accurately predict the effect of perturbation in single-cell level.

METHODS

scPerb utilizes a style transfer strategy by incorporating a style encoder into the architecture of a variational autoencoder. The style encoder captures the differences in latent representations between unperturbed and perturbed cells, enabling accurate prediction of post-perturbation gene expression data.

RESULTS

Comprehensive comparisons with existing methods demonstrate that scPerb delivers improved performance and higher accuracy in predicting cellular responses to perturbations. Notably, scPerb outperforms other methods across multiple datasets, achieving superior R2 values of 0.98, 0.98, and 0.96 on three benchmarking datasets.

CONCLUSION

scPerb offers a significant advancement in predicting cellular responses by effectively separating and transferring perturbation-related variances. This framework not only enhances prediction accuracy but also provides a robust tool for computational biology, addressing the limitations of current methodologies.

MATES: a deep learning-based model for locus-specific quantification of transposable elements in single cell

Transposable elements (TEs) are crucial for genetic diversity and gene regulation. Current single-cell quantification methods often align multi-mapping reads to either 'best-mapped' or 'random-mapped' locations and categorize them at the subfamily levels, overlooking the biological necessity for accurate, locus-specific TE quantification. Moreover, these existing methods are primarily designed for and focused on transcriptomics data, which restricts their adaptability to single-cell data of other modalities. To address these challenges, here we introduce MATES, a deep-learning approach that accurately allocates multi-mapping reads to specific loci of TEs, utilizing context from adjacent read alignments flanking the TE locus. When applied to diverse single-cell omics datasets, MATES shows improved performance over existing methods, enhancing the accuracy of TE quantification and aiding in the identification of marker TEs for identified cell populations. This development facilitates the exploration of single-cell heterogeneity and gene regulation through the lens of TEs, offering an effective transposon quantification tool for the single-cell genomics community.

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

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

Investigating the cis-regulatory basis of C3 and C4 photosynthesis in grasses at single-cell resolution

While considerable knowledge exists about the enzymes pivotal for C4 photosynthesis, much less is known about the cis-regulation important for specifying their expression in distinct cell types. Here, we use single-cell-indexed ATAC-seq to identify cell-type-specific accessible chromatin regions (ACRs) associated with C4 enzymes for five different grass species. This study spans four C4 species, covering three distinct photosynthetic subtypes: Zea mays and Sorghum bicolor (NADP-dependent malic enzyme), Panicum miliaceum (NAD-dependent malic enzyme), Urochloa fusca (phosphoenolpyruvate carboxykinase), along with the C3 outgroup Oryza sativa. We studied the cis-regulatory landscape of enzymes essential across all C4 species and those unique to C4 subtypes, measuring cell-type-specific biases for C4 enzymes using chromatin accessibility data. Integrating these data with phylogenetics revealed diverse co-option of gene family members between species, showcasing the various paths of C4 evolution. Besides promoter proximal ACRs, we found that, on average, C4 genes have two to three distal cell-type-specific ACRs, highlighting the complexity and divergent nature of C4 evolution. Examining the evolutionary history of these cell-type-specific ACRs revealed a spectrum of conserved and novel ACRs, even among closely related species, indicating ongoing evolution of cis-regulation at these C4 loci. This study illuminates the dynamic and complex nature of cis-regulatory elements evolution in C4 photosynthesis, particularly highlighting the intricate cis-regulatory evolution of key loci. Our findings offer a valuable resource for future investigations, potentially aiding in the optimization of C3 crop performance under changing climatic conditions.

Single-cell transcriptome and chromatin accessibility mapping of upper lip and primary palate fusion

Cleft lip and/or primary palate (CL/P) represent a prevalent congenital malformation, the aetiology of which is highly intricate. Although it is generally accepted that the condition arises from failed fusion between the upper lip and primary palate, the precise mechanism underlying this fusion process remains enigmatic. In this study, we utilized transposase-accessible chromatin sequencing (scATAC-seq) and single-cell RNA sequencing (scRNA-seq) to interrogate lambdoidal junction tissue derived from C57BL/6J mouse embryos at critical stages of embryogenesis (10.5, 11.5 and 12.5 embryonic days). We successfully identified distinct subgroups of mesenchymal and ectodermal cells involved in the fusion process and characterized their unique transcriptional profiles. Furthermore, we conducted cell differentiation trajectory analysis, revealing a dynamic repertoire of genes that are sequentially activated or repressed during pseudotime, facilitating the transition of relevant cell types. Additionally, we employed scATAC data to identify key genes associated with the fusion process and demonstrated differential chromatin accessibility across major cell types. Finally, we constructed a dynamic intercellular communication network and predicted upstream transcriptional regulators of critical genes involved in important signalling pathways. Our findings provide a valuable resource for future studies on upper lip and primary palate development, as well as congenital defects.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generate single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Epigenetics Research in Evolutionary Biology: Perspectives on Timescales and Mechanisms

Epigenetics research in evolutionary biology encompasses a variety of research areas, from regulation of gene expression to inheritance of environmentally mediated phenotypes. Such divergent research foci can occasionally render the umbrella term "epigenetics" ambiguous. Here I discuss several areas of contemporary epigenetics research in the context of evolutionary biology, aiming to provide balanced views across timescales and molecular mechanisms. The importance of epigenetics in development is now being assessed in many nonmodel species. These studies not only confirm the importance of epigenetic marks in developmental processes, but also highlight the significant diversity in epigenetic regulatory mechanisms across taxa. Further, these comparative epigenomic studies have begun to show promise toward enhancing our understanding of how regulatory programs evolve. A key property of epigenetic marks is that they can be inherited along mitotic cell lineages, and epigenetic differences that occur during early development can have lasting consequences on the organismal phenotypes. Thus, epigenetic marks may play roles in short-term (within an organism's lifetime or to the next generation) adaptation and phenotypic plasticity. However, the extent to which observed epigenetic variation occurs independently of genetic influences remains uncertain, due to the widespread impact of genetics on epigenetic variation and the limited availability of comprehensive (epi)genomic resources from most species. While epigenetic marks can be inherited independently of genetic sequences in some species, there is little evidence that such "transgenerational inheritance" is a general phenomenon. Rather, molecular mechanisms of epigenetic inheritance are highly variable between species.

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

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

Single-cell analysis of chromatin and expression reveals age- and sex-associated alterations in the human heart

Sex differences and age-related changes in the human heart at the tissue, cell, and molecular level have been well-documented and many may be relevant for cardiovascular disease. However, how molecular programs within individual cell types vary across individuals by age and sex remains poorly characterized. To better understand this variation, we performed single-nucleus combinatorial indexing (sci) ATAC- and RNA-Seq in human heart samples from nine donors. We identify hundreds of differentially expressed genes by age and sex and find epigenetic signatures of variation in ATAC-Seq data in this discovery cohort. We then scale up our single-cell RNA-Seq analysis by combining our data with five recently published single nucleus RNA-Seq datasets of healthy adult hearts. We find variation such as metabolic alterations by sex and immune changes by age in differential expression tests, as well as alterations in abundance of cardiomyocytes by sex and neurons with age. In addition, we compare our adult-derived ATAC-Seq profiles to analogous fetal cell types to identify putative developmental-stage-specific regulatory factors. Finally, we train predictive models of cell-type-specific RNA expression levels utilizing ATAC-Seq profiles to link distal regulatory sequences to promoters, quantifying the predictive value of a simple TF-to-expression regulatory grammar and identifying cell-type-specific TFs. Our analysis represents the largest single-cell analysis of cardiac variation by age and sex to date and provides a resource for further study of healthy cardiac variation and transcriptional regulation at single-cell resolution.

Benchmarking computational methods for single-cell chromatin data analysis

BACKGROUND

Single-cell chromatin accessibility assays, such as scATAC-seq, are increasingly employed in individual and joint multi-omic profiling of single cells. As the accumulation of scATAC-seq and multi-omics datasets continue, challenges in analyzing such sparse, noisy, and high-dimensional data become pressing. Specifically, one challenge relates to optimizing the processing of chromatin-level measurements and efficiently extracting information to discern cellular heterogeneity. This is of critical importance, since the identification of cell types is a fundamental step in current single-cell data analysis practices.

RESULTS

We benchmark 8 feature engineering pipelines derived from 5 recent methods to assess their ability to discover and discriminate cell types. By using 10 metrics calculated at the cell embedding, shared nearest neighbor graph, or partition levels, we evaluate the performance of each method at different data processing stages. This comprehensive approach allows us to thoroughly understand the strengths and weaknesses of each method and the influence of parameter selection.

CONCLUSIONS

Our analysis provides guidelines for choosing analysis methods for different datasets. Overall, feature aggregation, SnapATAC, and SnapATAC2 outperform latent semantic indexing-based methods. For datasets with complex cell-type structures, SnapATAC and SnapATAC2 are preferred. With large datasets, SnapATAC2 and ArchR are most scalable.

MOCHA's advanced statistical modeling of scATAC-seq data enables functional genomic inference in large human cohorts

Single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) is being increasingly used to study gene regulation. However, major analytical gaps limit its utility in studying gene regulatory programs in complex diseases. In response, MOCHA (Model-based single cell Open CHromatin Analysis) presents major advances over existing analysis tools, including: 1) improving identification of sample-specific open chromatin, 2) statistical modeling of technical drop-out with zero-inflated methods, 3) mitigation of false positives in single cell analysis, 4) identification of alternative transcription-starting-site regulation, and 5) modules for inferring temporal gene regulatory networks from longitudinal data. These advances, in addition to open chromatin analyses, provide a robust framework after quality control and cell labeling to study gene regulatory programs in human disease. We benchmark MOCHA with four state-of-the-art tools to demonstrate its advances. We also construct cross-sectional and longitudinal gene regulatory networks, identifying potential mechanisms of COVID-19 response. MOCHA provides researchers with a robust analytical tool for functional genomic inference from scATAC-seq data.

Human Vascularized Macrophage-Islet Organoids to Model Immune-Mediated Pancreatic β cell Pyroptosis upon Viral Infection

There is a paucity of human models to study immune-mediated host damage. Here, we utilized the GeoMx spatial multi-omics platform to analyze immune cell changes in COVID-19 pancreatic autopsy samples, revealing an accumulation of proinflammatory macrophages. Single cell RNA-seq analysis of human islets exposed to SARS-CoV-2 or Coxsackievirus B4 (CVB4) viruses identified activation of proinflammatory macrophages and β cell pyroptosis. To distinguish viral versus proinflammatory macrophage-mediated β cell pyroptosis, we developed human pluripotent stem cell (hPSC)-derived vascularized macrophage-islet (VMI) organoids. VMI organoids exhibited enhanced marker expression and function in both β cells and endothelial cells compared to separately cultured cells. Notably, proinflammatory macrophages within VMI organoids induced β cell pyroptosis. Mechanistic investigations highlighted TNFSF12-TNFRSF12A involvement in proinflammatory macrophage-mediated β cell pyroptosis. This study established hPSC-derived VMI organoids as a valuable tool for studying immune cell-mediated host damage and uncovered mechanism of β cell damage during viral exposure.

SMARCA5 reprograms AKR1B1-mediated fructose metabolism to control leukemogenesis

A significant variation in chromatin accessibility is an epigenetic feature of leukemia. The cause of this variation in leukemia, however, remains elusive. Here, we identify SMARCA5, a core ATPase of the imitation switch (ISWI) chromatin remodeling complex, as being responsible for aberrant chromatin accessibility in leukemia cells. We find that SMARCA5 is required to maintain aberrant chromatin accessibility for leukemogenesis and then promotes transcriptional activation of AKR1B1, an aldo/keto reductase, by recruiting transcription co-activator DDX5 and transcription factor SP1. Higher levels of AKR1B1 are associated with a poor prognosis in leukemia patients and promote leukemogenesis by reprogramming fructose metabolism. Moreover, pharmacological inhibition of AKR1B1 has been shown to have significant therapeutic effects in leukemia mice and leukemia patient cells. Thus, our findings link the aberrant chromatin state mediated by SMARCA5 to AKR1B1-mediated endogenous fructose metabolism reprogramming and shed light on the essential role of AKR1B1 in leukemogenesis, which may provide therapeutic strategies for leukemia.

Genome-Wide Screening Approaches for Biochemical Reactions Independent of Cell Growth

Genome-wide screening is a potent approach for comprehensively understanding the molecular mechanisms of biological phenomena. However, despite its widespread use in the past decades across various biological targets, its application to biochemical reactions with temporal and reversible biological outputs remains a formidable challenge. To uncover the molecular machinery underlying various biochemical reactions, we have recently developed the revival screening method, which combines flow cytometry-based cell sorting with library reconstruction from collected cells. Our refinements to the traditional genome-wide screening technique have proven successful in revealing the molecular machinery of biochemical reactions of interest. In this article, we elucidate the technical basis of revival screening, focusing on its application to CRISPR-Cas9 single guide RNA (sgRNA) library screening. Finally, we also discuss the future of genome-wide screening while describing recent achievements from in vitro and in vivo screening.

Detecting novel cell type in single-cell chromatin accessibility data via open-set domain adaptation

Recent advances in single-cell technologies enable the rapid growth of multi-omics data. Cell type annotation is one common task in analyzing single-cell data. It is a challenge that some cell types in the testing set are not present in the training set (i.e. unknown cell types). Most scATAC-seq cell type annotation methods generally assign each cell in the testing set to one known type in the training set but neglect unknown cell types. Here, we present OVAAnno, an automatic cell types annotation method which utilizes open-set domain adaptation to detect unknown cell types in scATAC-seq data. Comprehensive experiments show that OVAAnno successfully identifies known and unknown cell types. Further experiments demonstrate that OVAAnno also performs well on scRNA-seq data. Our codes are available online at https://github.com/lisaber/OVAAnno/tree/master.

Evidence of innate training in bovine γδ T cells following subcutaneous BCG administration

The Bacillus Calmette Guerin (BCG) vaccine has been shown to induce non-specific protection against diseases other than tuberculosis in vaccinated individuals, attributed to the induction of trained immunity. We have previously demonstrated that BCG administration induces innate immune training in mixed peripheral blood mononuclear cells and monocytes in calves. Gamma Delta (γδ) T cells are non-conventional T cells that exhibit innate and adaptive immune system features. They are in higher proportion in the peripheral blood of cattle than humans or rodents and play an essential role in bovine immune response to pathogens. In the current study, we determined if BCG administration induced innate immune training in bovine γδ T cells. A group of 16 pre-weaned Holstein calves (2-4 d age) were enrolled in the study and randomly assigned to vaccine and control groups (n=8/group). The vaccine group received two doses of 106 colony forming units (CFU) BCG Danish strain subcutaneously, separated by 2 weeks. The control group remained unvaccinated. Gamma delta T cells were purified from peripheral blood using magnetic cell sorting three weeks after receiving the 1st BCG dose. We observed functional changes in the γδ T cells from BCG-treated calves shown by increased IL-6 and TNF-α cytokine production in response to in vitro stimulation with Escherichia coli LPS and PAM3CSK4. ATAC-Seq analysis of 78,278 regions of open chromatin (peaks) revealed that γδ T cells from BCG-treated calves had an altered epigenetic status compared to cells from the control calves. Differentially accessible peaks (DAP) found near the promoters of innate immunity-related genes like Siglec14, Irf4, Ifna2, Lrrfip1, and Tnfrsf10d were 1 to 4-fold more accessible in cells from BCG-treated calves. MOTIF enrichment analysis of the sequences within DAPs, which explores transcription factor binding motifs (TFBM) upstream of regulatory elements, revealed TFBM for Eomes and IRF-5 were among the most enriched transcription factors. GO enrichment analysis of genes proximal to the DAPs showed enrichment of pathways such as regulation of IL-2 production, T-cell receptor signaling pathway, and other immune regulatory pathways. In conclusion, our study shows that subcutaneous BCG administration in pre-weaned calves can induce innate immune memory in the form of trained immunity in γδ T cells. This memory is associated with increased chromatin accessibility of innate immune response-related genes, thereby inducing a functional trained immune response evidenced by increased IL-6 and TNF-α cytokine production.

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

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

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

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