Multiple Single-cell Datasets Research Articles

SCIntRuler: Guiding the integration of multiple single-cell RNA-seq datasets with a novel statistical metric.

The growing number of single-cell RNA-seq (scRNA-seq) studies highlights the potential benefits of integrating multiple datasets, such as augmenting sample sizes and enhancing analytical robustness. Inherent diversity and batch discrepancies within samples or across studies continue to pose significant challenges for computational analyses. Questions persist in practice, lacking definitive answers: Should we use a specific integration method or opt for simply merging the datasets during joint analysis? Among all the existing data integration methods, which one is more suitable in specific scenarios? To fill the gap, we introduce SCIntRuler, a novel statistical metric for guiding the integration of multiple scRNA-seq datasets. SCIntRuler helps researchers make informed decisions regarding the necessity of data integration and the selection of an appropriate integration method. Our simulations and real data applications demonstrate that SCIntRuler streamlines decision-making processes and facilitates the analysis of diverse scRNA-seq datasets under varying contexts, thereby alleviating the complexities associated with the integration of heterogeneous scRNA-seq datasets. The implementation of our method is available on CRAN as an open-source R package with a user- friendly manual available: https://cloud.r-project.org/web/packages/SCIntRuler/index.html.

Identification of a Macrophage marker gene signature to evaluate immune infiltration and therapeutic response in hepatocellular carcinoma

BackgroundOnly a minority of hepatocellular carcinoma (HCC) patients can benefit from systemic regimens. Macrophages, which abundantly infiltrate in HCC, could mediate tumour microenvironment remodelling and immune escape, proving to be powerful weapons in combating HCC. Thus, a deeper understanding of macrophages is necessary for improving existing antitumour treatments. MethodsWith a series of bioinformatic approaches, we comprehensively explored the role of macrophage-related genes in human HCCs from multiple single-cell and bulk RNA sequencing datasets. Unsupervised clustering was performed to cluster the macrophage marker genes (MMGs). GSVA and functional enrichment analysis were used to elucidate the functional differences among the MMG-associated clusters. Subsequently, a component analysis algorithm was used to construct a Macrosig score, and the prognosis, biological characteristics, mutation profile, TME cell infiltration status and drug response of patients with different Macrosig scores were further analysed. ResultsWe identified 13 MMGs in 574 HCC samples, based on which three MMG-associated clusters were defined. Overall survival time, clinicopathological features and immune infiltration scores differed among the different clusters. On this basis, 12 hub genes were identified among these clusters; subsequently, a scoring system was constructed to determine the Macrosig score. Importantly, patients with low-Macrosig scores, characterized by increased immune infiltration, increased mutation frequency and increased immune checkpoint expression, including CTLA-4, LAG3, PDCD1 and TIGIT, exhibited enhanced efficacy of immunotherapy when validated in an external database. Moreover, a low-Macrosig score indicates increased sensitivity to AZD.2281, A.443654, ABT.263, ABT.888, AG.014699 and ATRA, while a high Macrosig score indicates increased sensitivity to AZD6482, AKT inhibitor VIII, AS601245, AZ628, AZD.0530 and AZD6244. ConclusionsA novel scoring system was constructed to guide more effective prognostic evaluation and tailoring therapeutic regimens for HCC patients.

Beaconet: A Reference-Free Method for Integrating Multiple Batches of Single-Cell Transcriptomic Data in Original Molecular Space.

Integrating multiple single-cell datasets is essential for the comprehensive understanding of cell heterogeneity. Batch effect is the undesired systematic variations among technologies or experimental laboratories that distort biological signals and hinder the integration of single-cell datasets. However, existing methods typically rely on a selected dataset as a reference, leading to inconsistent integration performance using different references, or embed cells into uninterpretable low-dimensional feature space. To overcome these limitations, a reference-free method, Beaconet, for integrating multiple single-cell transcriptomic datasets in original molecular space by aligning the global distribution of each batch using an adversarial correction network is presented. Through extensive comparisons with 13 state-of-the-art methods, it is demonstrated that Beaconet can effectively remove batch effect while preserving biological variations and is superior to existing unsupervised methods using all possible references in overall performance. Furthermore, Beaconet performs integration in the original molecular feature space, enabling the characterization of cell types and downstream differential expression analysis directly using integrated data with gene-expression features. Additionally, when applying to large-scale atlas data integration, Beaconet shows notable advantages in both time- and space-efficiencies. In summary, Beaconet serves as an effective and efficient batch effect removal tool that can facilitate the integration of single-cell datasets in a reference-free and molecular feature-preserved mode.

CDSKNNXMBD: a novel clustering framework for large-scale single-cell data based on a stable graph structure

BackgroundAccurate and efficient cell grouping is essential for analyzing single-cell transcriptome sequencing (scRNA-seq) data. However, the existing clustering techniques often struggle to provide timely and accurate cell type groupings when dealing with datasets with large-scale or imbalanced cell types. Therefore, there is a need for improved methods that can handle the increasing size of scRNA-seq datasets while maintaining high accuracy and efficiency.MethodsWe propose CDSKNNXMBD (Community Detection based on a Stable K-Nearest Neighbor Graph Structure), a novel single-cell clustering framework integrating partition clustering algorithm and community detection algorithm, which achieves accurate and fast cell type grouping by finding a stable graph structure.ResultsWe evaluated the effectiveness of our approach by analyzing 15 tissues from the human fetal atlas. Compared to existing methods, CDSKNN effectively counteracts the high imbalance in single-cell data, enabling effective clustering. Furthermore, we conducted comparisons across multiple single-cell datasets from different studies and sequencing techniques. CDSKNN is of high applicability and robustness, and capable of balancing the complexities of across diverse types of data. Most importantly, CDSKNN exhibits higher operational efficiency on datasets at the million-cell scale, requiring an average of only 6.33 min for clustering 1.46 million single cells, saving 33.3% to 99% of running time compared to those of existing methods.ConclusionsThe CDSKNN is a flexible, resilient, and promising clustering tool that is particularly suitable for clustering imbalanced data and demonstrates high efficiency on large-scale scRNA-seq datasets.

Batch-effect correction in single-cell RNA sequencing data using JIVE.

In single-cell RNA sequencing analysis, addressing batch effects-technical artifacts stemming from factors such as varying sequencing technologies, equipment, and capture times-is crucial. These factors can cause unwanted variation and obfuscate the underlying biological signal of interest. The joint and individual variation explained (JIVE) method can be used to extract shared biological patterns from multi-source sequencing data while adjusting for individual non-biological variations (i.e. batch effect). However, its current implementation is originally designed for bulk sequencing data, making it computationally infeasible for large-scale single-cell sequencing datasets. In this study, we enhance JIVE for large-scale single-cell data by boosting its computational efficiency. Additionally, we introduce a novel application of JIVE for batch-effect correction on multiple single-cell sequencing datasets. Our enhanced method aims to decompose single-cell sequencing datasets into a joint structure capturing the true biological variability and individual structures, which capture technical variability within each batch. This joint structure is then suitable for use in downstream analyses. We benchmarked the results against four popular tools, Seurat v5, Harmony, LIGER, and Combat-seq, which were developed for this purpose. JIVE performed best in terms of preserving cell-type effects and in scenarios in which the batch sizes are balanced. The JIVE implementation used for this analysis can be found at https://github.com/oconnell-statistics-lab/scJIVE.

Comprehensive Integration of Multiple Single-Cell Transcriptomic Data Sets Defines Distinct Cell Populations and Their Phenotypic Changes in Murine Atherosclerosis.

The application of single-cell transcriptomic (single-cell RNA sequencing) analysis to the study of atherosclerosis has provided unique insights into the molecular and genetic mechanisms that mediate disease risk and pathophysiology. However, nonstandardized methodologies and relatively high costs associated with the technique have limited the size and replication of existing data sets and created disparate or contradictory findings that have fostered misunderstanding and controversy. To address these uncertainties, we have performed a conservative integration of multiple published single-cell RNA sequencing data sets into a single meta-analysis, performed extended analysis of native resident vascular cells, and used in situ hybridization to map the disease anatomic location of the identified cluster cells. To investigate the transdifferentiation of smooth muscle cells to macrophage phenotype, we have developed a classifying algorithm based on the quantification of reporter transgene expression. The reporter gene expression tool indicates that within the experimental limits of the examined studies, transdifferentiation of smooth muscle cell to the macrophage lineage is extremely rare. Validated transition smooth muscle cell phenotypes were defined by clustering, and the location of these cells was mapped to lesion anatomy with in situ hybridization. We have also characterized 5 endothelial cell phenotypes and linked these cellular species to different vascular structures and functions. Finally, we have identified a transcriptomically unique cellular phenotype that constitutes the aortic valve. Taken together, these analyses resolve a number of outstanding issues related to differing results reported with vascular disease single-cell RNA sequencing studies, and significantly extend our understanding of the role of resident vascular cells in anatomy and disease.

The Deep Generative Decoder: MAP estimation of representations improves modelling of single-cell RNA data.

Learning low-dimensional representations of single-cell transcriptomics has become instrumental to its downstream analysis. The state of the art is currently represented by neural network models, such as variational autoencoders, which use a variational approximation of the likelihood for inference. We here present the Deep Generative Decoder (DGD), a simple generative model that computes model parameters and representations directly via maximum a posteriori estimation. The DGD handles complex parameterized latent distributions naturally unlike variational autoencoders, which typically use a fixed Gaussian distribution, because of the complexity of adding other types. We first show its general functionality on a commonly used benchmark set, Fashion-MNIST. Secondly, we apply the model to multiple single-cell datasets. Here, the DGD learns low-dimensional, meaningful, and well-structured latent representations with sub-clustering beyond the provided labels. The advantages of this approach are its simplicity and its capability to provide representations of much smaller dimensionality than a comparable variational autoencoder. scDGD is available as a python package at https://github.com/Center-for-Health-Data-Science/scDGD. The remaining code is made available here: https://github.com/Center-for-Health-Data-Science/dgd.

Limited expression of Nrf2 in neurons across the central nervous system

Nrf2, encoded by the gene Nfe2l2, is a broadly expressed transcription factor that regulates gene expression in response to reactive oxygen species (ROS) and oxidative stress. It is commonly referred to as a ubiquitous pathway, but this generalization overlooks work indicating that Nrf2 is essentially unexpressed in some neuronal populations. To explore whether this pattern extends throughout the central nervous system (CNS), we quantified Nfe2l2 expression and chromatin accessibility at the Nfe2l2 locus across multiple single cell datasets. In both the mouse and human CNS, Nfe2l2 was repressed in almost all mature neurons, but highly expressed in non-neuronal support cells, and this pattern was robust across multiple human CNS diseases. A subset of key Nrf2 target genes, like Slc7a11, also remained low in neurons. Thus, these data suggest that while most cells express Nfe2l2, with activity determined by ROS levels, neurons actively avoid Nrf2 activity by keeping Nfe2l2 expression low.

Inferring regulators of cell identity in the human adult pancreas.

Cellular identity during development is under the control of transcription factors that form gene regulatory networks. However, the transcription factors and gene regulatory networks underlying cellular identity in the human adult pancreas remain largely unexplored. Here, we integrate multiple single-cell RNA-sequencing datasets of the human adult pancreas, totaling 7393 cells, and comprehensively reconstruct gene regulatory networks. We show that a network of 142 transcription factors forms distinct regulatory modules that characterize pancreatic cell types. We present evidence that our approach identifies regulators of cell identity and cell states in the human adult pancreas. We predict that HEYL, BHLHE41 and JUND are active in acinar, beta and alpha cells, respectively, and show that these proteins are present in the human adult pancreas as well as in human induced pluripotent stem cell (hiPSC)-derived islet cells. Using single-cell transcriptomics, we found that JUND represses beta cell genes in hiPSC-alpha cells. BHLHE41 depletion induced apoptosis in primary pancreatic islets. The comprehensive gene regulatory network atlas can be explored interactively online. We anticipate our analysis to be the starting point for a more sophisticated dissection of how transcription factors regulate cell identity and cell states in the human adult pancreas.

Divergent expression of Neurl3 from hemogenic endothelial cells to hematopoietic stem progenitor cells during development

Prior to the generation of hematopoietic stem cells (HSCs) from the hemogenic endothelial cells (HECs) mainly in the dorsal aorta in midgestational mouse embryos, multiple hematopoietic progenitors including erythro-myeloid progenitors and lymphoid progenitors are generated from yolk sac HECs. These HSC-independent hematopoietic progenitors have recently been identified as major contributors to functional blood cell production until birth. However, little is known about yolk sac HECs. Here, combining integrative analyses of multiple single-cell RNA-sequencing datasets and functional assays, we reveal that Neurl3-EGFP, in addition to marking the continuum throughout the ontogeny of HSCs from HECs, can also serve as a single enrichment marker for yolk sac HECs. Moreover, while yolk sac HECs have much weaker arterial characteristics than either arterial endothelial cells in the yolk sac or HECs within the embryo proper, the lymphoid potential of yolk sac HECs is largely confined to the arterial-biased subpopulation featured by the Unc5b expression. Interestingly, the B lymphoid potential of hematopoietic progenitors, but not for myeloid potentials, is exclusively detected in Neurl3-negative subpopulations in midgestational embryos. Taken together, these findings enhance our understanding of blood birth from yolk sac HECs and provide theoretical basis and candidate reporters for monitoring step-wise hematopoietic differentiation.

Quantifying Cell-Type-Specific Differences of Single-Cell Datasets Using Uniform Manifold Approximation and Projection for Dimension Reduction and Shapley Additive exPlanations.

With rapid advances in single-cell profiling technologies, larger-scale investigations that require comparisons of multiple single-cell datasets can lead to novel findings. Specifically, quantifying cell-type-specific responses to different conditions across single-cell datasets could be useful in understanding how the difference in conditions is induced at a cellular level. In this study, we present a computational pipeline that quantifies cell-type-specific differences and identifies genes responsible for the differences. We quantify differences observed in a low-dimensional uniform manifold approximation and projection for dimension reduction space as a proxy for the difference present in the high-dimensional space and use SHapley Additive exPlanations to quantify genes driving the differences. In this study, we applied our algorithm to the Iris flower dataset, single-cell RNA sequencing dataset, and mass cytometry dataset and demonstrate that it can robustly quantify cell-type-specific differences and it can also identify genes that are responsible for the differences.

Illustrating the biological functions and diagnostic value of transmembrane protein family members in glioma.

It is well-established that patients with glioma have a poor prognosis. Although the past few decades have witnessed unprecedented medical advances, the 5-year survival remains dismally low. This study aims to investigate the role of transmembrane protein-related genes in the development and prognosis of glioma and provide new insights into the pathogenesis of the disease. The datasets of glioma patients, including RNA sequencing data and relative clinical information, were obtained from The Cancer Genome Atlas (TCGA), Chinese Glioma Genome Atlas (CGGA) and Gene Expression Omnibus (GEO) databases. Prognostic transmembrane protein-related genes were identified by univariate Cox analysis. New disease subtypes were recognized based on the consensus clustering method, and their biological uniqueness was verified via various algorithms. The prognosis signature was constructed using the LASSO-Cox regression model, and its predictive power was validated in external datasets by receiver operating characteristic (ROC) curve analysis. An independent prognostic analysis was conducted to evaluate whether the signature could be considered a prognostic factor independent of other variables. A nomogram was constructed in conjunction with traditional clinical variables. The concordance index (C-index) and Decision Curve Analysis (DCA) were used to assess the net clinical benefit of the signature over traditional clinical variables. Seven different softwares were used to compare the differences in immune infiltration between the high- and low-risk groups to explore potential mechanisms of glioma development and prognosis. Hub genes were found using the random forest method, and their expression was based on multiple single-cell datasets. Four molecular subtypes were identified, among which the C1 group had the worst prognosis. Principal Component Analysis (PCA) results and heatmaps indicated that prognosis-related transmembrane protein genes exhibited differential expression in all four groups. Besides, the microenvironment of the four groups exhibited significant heterogeneity. The 6 gene-based signatures could predict the 1-, 2-, and 3-year overall survival (OS) of glioma patients. The signature could be used as an independent prognosis factor of glioma OS and was superior to traditional clinical variables. More immune cells were infiltrated in the high-risk group, suggesting immune escape. According to our signature, many genes were associated with the content of immune cells, which revealed that transmembrane protein-related genes might influence the development and prognosis of glioma by regulating the immune microenvironment. TMEM158 was identified as the most important gene using the random forest method. The single-cell datasets consistently showed that TMEM158 was expressed in multiple malignant cells. The expression of transmembrane protein-related genes is closely related to the immune status and prognosis of glioma patients by regulating tumor progression in various ways. The interaction between transmembrane protein-related genes and immunity during glioma development lays the groundwork for future studies on the molecular mechanism and targeted therapy of glioma.

Multi-modal quantification of pathway activity with MAYA

Signaling pathways can be activated through various cascades of genes depending on cell identity and biological context. Single-cell atlases now provide the opportunity to inspect such complexity in health and disease. Yet, existing reference tools for pathway scoring resume activity of each pathway to one unique common metric across cell types. Here, we present MAYA, a computational method that enables the automatic detection and scoring of the diverse modes of activation of biological pathways across cell populations. MAYA improves the granularity of pathway analysis by detecting subgroups of genes within reference pathways, each characteristic of a cell population and how it activates a pathway. Using multiple single-cell datasets, we demonstrate the biological relevance of identified modes of activation, the robustness of MAYA to noisy pathway lists and batch effect. MAYA can also predict cell types starting from lists of reference markers in a cluster-free manner. Finally, we show that MAYA reveals common modes of pathway activation in tumor cells across patients, opening the perspective to discover shared therapeutic vulnerabilities.

Integrating Multiple Single-Cell RNA Sequencing Datasets Using Adversarial Autoencoders.

Single-cell RNA sequencing (RNA-seq) has been demonstrated to be a proven method for quantifying gene-expression heterogeneity and providing insight into the transcriptome at the single-cell level. When combining multiple single-cell transcriptome datasets for analysis, it is common to first correct the batch effect. Most of the state-of-the-art processing methods are unsupervised, i.e., they do not utilize single-cell cluster labeling information, which could improve the performance of batch correction methods, especially in the case of multiple cell types. To better utilize known labels for complex dataset scenarios, we propose a novel deep learning model named IMAAE (i.e., integrating multiple single-cell datasets via an adversarial autoencoder) to correct the batch effects. After conducting experiments with various dataset scenarios, the results show that IMAAE outperforms existing methods for both qualitative measures and quantitative evaluation. In addition, IMAAE is able to retain both corrected dimension reduction data and corrected gene expression data. These features make it a potential new option for large-scale single-cell gene expression data analysis.

Review of single-cell RNA-seq data clustering for cell-type identification and characterization.

In recent years, the advances in single-cell RNA-seq techniques have enabled us to perform large-scale transcriptomic profiling at single-cell resolution in a high-throughput manner. Unsupervised learning such as data clustering has become the central component to identify and characterize novel cell types and gene expression patterns. In this study, we review the existing single-cell RNA-seq data clustering methods with critical insights into the related advantages and limitations. In addition, we also review the upstream single-cell RNA-seq data processing techniques such as quality control, normalization, and dimension reduction. We conduct performance comparison experiments to evaluate several popular single-cell RNA-seq clustering approaches on simulated and multiple single-cell transcriptomic data sets.

SCIBER: a simple method for removing batch effects from single-cell RNA-sequencing data.

Integrative analysis of multiple single-cell RNA-sequencing datasets allows for more comprehensive characterizations of cell types, but systematic technical differences between datasets, known as 'batch effects', need to be removed before integration to avoid misleading interpretation of the data. Although many batch-effect-removal methods have been developed, there is still a large room for improvement: most existing methods only give dimension-reduced data instead of expression data of individual genes, are based on computationally demanding models and are black-box models and thus difficult to interpret or tune. Here, we present a new batch-effect-removal method called SCIBER (Single-Cell Integrator and Batch Effect Remover) and study its performance on real datasets. SCIBER matches cell clusters across batches according to the overlap of their differentially expressed genes. As a simple algorithm that has better scalability to data with a large number of cells and is easy to tune, SCIBER shows comparable and sometimes better accuracy in removing batch effects on real datasets compared to the state-of-the-art methods, which are much more complicated. Moreover, SCIBER outputs expression data in the original space, that is, the expression of individual genes, which can be used directly for downstream analyses. Additionally, SCIBER is a reference-based method, which assigns one of the batches as the reference batch and keeps it untouched during the process, making it especially suitable for integrating user-generated datasets with standard reference data such as the Human Cell Atlas. SCIBER is publicly available as an R package on CRAN: https://cran.r-project.org/web/packages/SCIBER/. A vignette is included in the CRAN R package. Supplementary data are available at Bioinformatics online.

ScAB detects multiresolution cell states with clinical significance by integrating single-cell genomics and bulk sequencing data.

Although single-cell sequencing has provided a powerful tool to deconvolute cellular heterogeneity of diseases like cancer, extrapolating clinical significance or identifying clinically-relevant cells remains challenging. Here, we propose a novel computational method scAB, which integrates single-cell genomics data with clinically annotated bulk sequencing data via a knowledge- and graph-guided matrix factorization model. Once combined, scAB provides a coarse- and fine-grain multiresolution perspective of phenotype-associated cell states and prognostic signatures previously not visible by single-cell genomics. We use scAB to enhance live cancer single-cell RNA-seq data, identifying clinically-relevant previously unrecognized cancer and stromal cell subsets whose signatures show a stronger poor-survival association. The identified fine-grain cell subsets are associated with distinct cancer hallmarks and prognosis power. Furthermore, scAB demonstrates its utility as a biomarker identification tool, with the ability to predict immunotherapy, drug responses and survival when applied to melanoma single-cell RNA-seq datasets and glioma single-cell ATAC-seq datasets. Across multiple single-cell and bulk datasets from different cancer types, we also demonstrate the superior performance of scAB in generating prognosis signatures and survival predictions over existing models. Overall, scAB provides an efficient tool for prioritizing clinically-relevant cell subsets and predictive signatures, utilizing large publicly available databases to improve prognosis and treatments.

Claudin5 protects the peripheral endothelial barrier in an organ and vessel-type-specific manner.

Dysfunctional and leaky blood vessels resulting from disruption of the endothelial cell (EC) barrier accompanies numerous diseases. The EC barrier is established through endothelial cell tight and adherens junctions. However, the expression pattern and precise contribution of different junctional proteins to the EC barrier is poorly understood. Here, we focus on organs with continuous endothelium to identify structural and functional in vivo characteristics of the EC barrier. Assembly of multiple single-cell RNAseq datasets into a single integrated database revealed the variability and commonalities of EC barrier patterning. Across tissues, Claudin5 exhibited diminishing expression along the arteriovenous axis, correlating with EC barrier integrity. Functional analysis identified tissue-specific differences in leakage properties and response to the leakage agonist histamine. Loss of Claudin5 enhanced histamine-induced leakage in an organotypic and vessel type-specific manner in an inducible, EC-specific, knock-out mouse. Mechanistically, Claudin5 loss left junction ultrastructure unaffected but altered its composition, with concomitant loss of zonula occludens-1 and upregulation of VE-Cadherin expression. These findings uncover the organ-specific organisation of the EC barrier and distinct importance of Claudin5 in different vascular beds, providing insights to modify EC barrier stability in a targeted, organ-specific manner.

Identifying tumor cells at the single-cell level using machine learning

Tumors are complex tissues of cancerous cells surrounded by a heterogeneous cellular microenvironment with which they interact. Single-cell sequencing enables molecular characterization of single cells within the tumor. However, cell annotation—the assignment of cell type or cell state to each sequenced cell—is a challenge, especially identifying tumor cells within single-cell or spatial sequencing experiments. Here, we propose ikarus, a machine learning pipeline aimed at distinguishing tumor cells from normal cells at the single-cell level. We test ikarus on multiple single-cell datasets, showing that it achieves high sensitivity and specificity in multiple experimental contexts.

Cell Type Diversity Statistic: An Entropy-Based Metric to Compare Overall Cell Type Composition Across Samples.

Changes of cell type composition across samples can carry biological significance and provide insight into disease and other conditions. Single cell transcriptomics has made it possible to study cell type composition at a fine resolution. Most single cell studies investigate compositional changes between samples for each cell type independently, not accounting for the fixed number of cells per sample in sequencing data. Here, we provide a metric of the distribution of cell type proportions in a sample that can be used to compare the overall distribution of cell types across multiple samples and biological conditions. This is the first method to measure overall cell type composition at the single cell level. We use the method to assess compositional changes in peripheral blood mononuclear cells (PBMCs) related to aging and extreme old age using multiple single cell datasets from individuals of four age groups across the human lifespan.