Joint analysis of heterogeneous single-cell RNA-seq dataset collections.
Ontology highlight
ABSTRACT:
SUBMITTER: Barkas N
PROVIDER: S-EPMC6684315 | biostudies-literature | 2019 Aug
REPOSITORIES: biostudies-literature
Publications
Joint analysis of heterogeneous single-cell RNA-seq dataset collections.
Barkas Nikolas N Petukhov Viktor V Nikolaeva Daria D Lozinsky Yaroslav Y Demharter Samuel S Khodosevich Konstantin K Kharchenko Peter V PV
Nature methods 20190715 8
Single-cell RNA sequencing is often applied in study designs that include multiple individuals, conditions or tissues. To identify recurrent cell subpopulations in such heterogeneous collections, we developed Conos, an approach that relies on multiple plausible inter-sample mappings to construct a global graph connecting all measured cells. The graph enables identification of recurrent cell clusters and propagation of information between datasets in multi-sample or atlas-scale collections. ...[more]
Similar Datasets
Project description:BackgroundGlioblastoma (GBM) is the most common primary malignant brain tumor in adults. For patients with GBM, the median overall survival (OS) is 14.6 months and the 5-year survival rate is 7.2%. It is imperative to develop a reliable model to predict the survival probability in new GBM patients. To date, most prognostic models for predicting survival in GBM were constructed based on bulk RNA-seq dataset, which failed to accurately reflect the difference between tumor cores and peripheral regions, and thus show low predictive capability. An effective prognostic model is desperately needed in clinical practice.MethodsWe studied single-cell RNA-seq dataset and The Cancer Genome Atlas-glioblastoma multiforme (TCGA-GBM) dataset to identify differentially expressed genes (DEGs) that impact the OS of GBM patients. We then applied the least absolute shrinkage and selection operator (LASSO) Cox penalized regression analysis to determine the optimal genes to be included in our risk score prognostic model. Then, we used another dataset to test the accuracy of our risk score prognostic model.ResultsWe identified 2128 DEGs from the single-cell RNA-seq dataset and 6461 DEGs from the bulk RNA-seq dataset. In addition, 896 DEGs associated with the OS of GBM patients were obtained. Five of these genes (LITAF, MTHFD2, NRXN3, OSMR, and RUFY2) were selected to generate a risk score prognostic model. Using training and validation datasets, we found that patients in the low-risk group showed better OS than those in the high-risk group. We validated our risk score model with the training and validating datasets and demonstrated that it can effectively predict the OS of GBM patients.ConclusionWe constructed a novel prognostic model to predict survival in GBM patients by integrating a scRNA-seq dataset and a bulk RNA-seq dataset. Our findings may advance the development of new therapeutic targets and improve clinical outcomes for GBM patients.
| S-EPMC9120724 | biostudies-literature
Project description:Dimension reduction and (spatial) clustering is usually performed sequentially; however, the low-dimensional embeddings estimated in the dimension-reduction step may not be relevant to the class labels inferred in the clustering step. We therefore developed a computation method, Dimension-Reduction Spatial-Clustering (DR-SC), that can simultaneously perform dimension reduction and (spatial) clustering within a unified framework. Joint analysis by DR-SC produces accurate (spatial) clustering results and ensures the effective extraction of biologically informative low-dimensional features. DR-SC is applicable to spatial clustering in spatial transcriptomics that characterizes the spatial organization of the tissue by segregating it into multiple tissue structures. Here, DR-SC relies on a latent hidden Markov random field model to encourage the spatial smoothness of the detected spatial cluster boundaries. Underlying DR-SC is an efficient expectation-maximization algorithm based on an iterative conditional mode. As such, DR-SC is scalable to large sample sizes and can optimize the spatial smoothness parameter in a data-driven manner. With comprehensive simulations and real data applications, we show that DR-SC outperforms existing clustering and spatial clustering methods: it extracts more biologically relevant features than conventional dimension reduction methods, improves clustering performance, and offers improved trajectory inference and visualization for downstream trajectory inference analyses.
| S-EPMC9262606 | biostudies-literature
Project description:SummaryA common bioinformatics task in single-cell data analysis is to purify a cell type or cell population of interest from heterogeneous datasets. Here, we present scGate, an algorithm that automatizes marker-based purification of specific cell populations, without requiring training data or reference gene expression profiles. scGate purifies a cell population of interest using a set of markers organized in a hierarchical structure, akin to gating strategies employed in flow cytometry. scGate outperforms state-of-the-art single-cell classifiers and it can be applied to multiple modalities of single-cell data (e.g. RNA-seq, ATAC-seq, CITE-seq). scGate is implemented as an R package and integrated with the Seurat framework, providing an intuitive tool to isolate cell populations of interest from heterogeneous single-cell datasets.Availability and implementationscGate is available as an R package at https://github.com/carmonalab/scGate (https://doi.org/10.5281/zenodo.6202614). Several reproducible workflows describing the main functions and usage of the package on different single-cell modalities, as well as the code to reproduce the benchmark, can be found at https://github.com/carmonalab/scGate.demo (https://doi.org/10.5281/zenodo.6202585) and https://github.com/carmonalab/scGate.benchmark. Test data are available at https://doi.org/10.6084/m9.figshare.16826071.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC9048671 | biostudies-literature
Project description:We introduce BacDrop, a highly scalable technology for bacterial single-cell RNA sequencing that has overcome many challenges hindering the development of scRNA-seq in bacteria. BacDrop can be applied to thousands to millions of cells from both gram-negative and gram-positive species. It features universal ribosomal RNA depletion and combinatorial barcodes that enable multiplexing and massively parallel sequencing. We applied BacDrop to study Klebsiella pneumoniae clinical isolates and to elucidate their heterogeneous responses to antibiotic stress. In an unperturbed population presumed to be homogeneous, we found within-population heterogeneity largely driven by the expression of mobile genetic elements that promote the evolution of antibiotic resistance. Under antibiotic perturbation, BacDrop revealed transcriptionally distinct subpopulations associated with different phenotypic outcomes including antibiotic persistence. BacDrop thus can capture cellular states that cannot be detected by bulk RNA-seq, which will unlock new microbiological insights into bacterial responses to perturbations and larger bacterial communities such as the microbiome.
| S-EPMC10014032 | biostudies-literature
Project description:BackgroundSingle-cell RNA-sequencing (scRNA-seq) measures gene expression in single cells, while single-nucleus ATAC-sequencing (snATAC-seq) quantifies chromatin accessibility in single nuclei. These two data types provide complementary information for deciphering cell types and states. However, when analyzed individually, they sometimes produce conflicting results regarding cell type/state assignment. The power is compromised since the two modalities reflect the same underlying biology. Recently, it has become possible to measure both gene expression and chromatin accessibility from the same nucleus. Such paired data enable the direct modeling of the relationships between the two modalities. Given the availability of the vast amount of single-modality data, it is desirable to integrate the paired and unpaired single-modality datasets to gain a comprehensive view of the cellular complexity.ResultsWe benchmark nine existing single-cell multi-omic data integration methods. Specifically, we evaluate to what extent the multiome data provide additional guidance for analyzing the existing single-modality data, and whether these methods uncover peak-gene associations from single-modality data. Our results indicate that multiome data are helpful for annotating single-modality data. However, we emphasize that the availability of an adequate number of nuclei in the multiome dataset is crucial for achieving accurate cell type annotation. Insufficient representation of nuclei may compromise the reliability of the annotations. Additionally, when generating a multiome dataset, the number of cells is more important than sequencing depth for cell type annotation.ConclusionsSeurat v4 is the best currently available platform for integrating scRNA-seq, snATAC-seq, and multiome data even in the presence of complex batch effects.
| S-EPMC10594700 | biostudies-literature
Project description:The ubiquity of RNA-seq has led to many methods that use RNA-seq data to analyze variations in RNA splicing. However, available methods are not well suited for handling heterogeneous and large datasets. Such datasets scale to thousands of samples across dozens of experimental conditions, exhibit increased variability compared to biological replicates, and involve thousands of unannotated splice variants resulting in increased transcriptome complexity. We describe here a suite of algorithms and tools implemented in the MAJIQ v2 package to address challenges in detection, quantification, and visualization of splicing variations from such datasets. Using both large scale synthetic data and GTEx v8 as benchmark datasets, we assess the advantages of MAJIQ v2 compared to existing methods. We then apply MAJIQ v2 package to analyze differential splicing across 2,335 samples from 13 brain subregions, demonstrating its ability to offer insights into brain subregion-specific splicing regulation.
| S-EPMC9984406 | biostudies-literature
Project description:Single-cell RNA-seq is a powerful tool for revealing heterogeneity in cancer cells. However, each of the current single-cell RNA-seq platforms has inherent advantages and disadvantages. Here, we show that combining the different single-cell RNA-seq platforms can be an effective approach to obtaining complete information about expression differences and a sufficient cellular population to understand transcriptional heterogeneity in cancers. We demonstrate that it is possible to estimate missing expression information. We further demonstrate that even in the cases where precise information for an individual gene cannot be inferred, the activity of given transcriptional modules can be analyzed. Interestingly, we found that two distinct transcriptional modules, one associated with the Aurora kinase gene and the other with the DUSP gene, are aberrantly regulated in a minor population of cells and may thus contribute to the possible emergence of dormancy or eventual drug resistance within the population.
| S-EPMC5823859 | biostudies-literature
Project description:Breast cancer is the most commonly diagnosed cancer in the world, with triple-negative breast cancer (TNBC) making up 12% of these diagnoses. TNBC tumours are highly heterogeneous in both inter-tumour and intra-tumour gene expression profiles, where they form subclonal populations of varying levels of aggressiveness. These aspects make it difficult to study and treat TNBC, requiring further research into tumour heterogeneity as well as potential therapeutic targets and biomarkers. Recently, it was discovered that the majority of the transcribed genome comprises non-coding RNAs, in particular long non-coding RNAs (lncRNAs). LncRNAs are transcripts of >200 nucleotides in length that do not encode a protein. They have been characterised as regulatory molecules and their expression can be associated with a malignant phenotype. We set out to explore TNBC tumour heterogeneity in vivo at a single cell level to investigate whether lncRNA expression varies across different cells within the tumour, even if cells are coming from the same cell line, and whether lncRNA expression is sufficient to define cellular subpopulations. We applied single-cell expression profiling due to its ability to capture expression signals of lncRNAs expressed in small subpopulations of cells. Overall, we observed most lncRNAs to be expressed at low, but detectable levels in TNBC xenografts, with a median of 25 lncRNAs detected per cell. LncRNA expression alone was insufficient to define a subpopulation of cells, and lncRNAs showed highly heterogeneous expression patterns, including ubiquitous expression, subpopulation-specific expression, and a hybrid pattern of lncRNAs expressed in several, but not all subpopulations. These findings reinforce that transcriptionally defined tumour cell subpopulations can be identified in cell-line derived xenografts, and uses single-cell RNA-seq (scRNA-seq) to detect and characterise lncRNA expression across these subpopulations in xenografted tumours. Future studies will aim to investigate the spatial distribution of lncRNAs within xenografts and patient tissues, and study the potential of subclone-specific lncRNAs as new therapeutic targets and/or biomarkers.
| S-EPMC8533545 | biostudies-literature
Project description:One primary reason that makes single-cell RNA-seq analysis challenging is dropouts, where the data only captures a small fraction of the transcriptome of each cell. Almost all computational algorithms developed for single-cell RNA-seq adopted gene selection, dimension reduction or imputation to address the dropouts. Here, an opposite view is explored. Instead of treating dropouts as a problem to be fixed, we embrace it as a useful signal. We represent the dropout pattern by binarizing single-cell RNA-seq count data, and present a co-occurrence clustering algorithm to cluster cells based on the dropout pattern. We demonstrate in multiple published datasets that the binary dropout pattern is as informative as the quantitative expression of highly variable genes for the purpose of identifying cell types. We expect that recognizing the utility of dropouts provides an alternative direction for developing computational algorithms for single-cell RNA-seq analysis.
| S-EPMC7054558 | biostudies-literature
Project description:Human red blood cells (RBCs), or erythrocytes, are the most abundant blood cells responsible for gas exchange. RBC diseases affect hundreds of millions of people and impose enormous financial and personal burdens. One well-recognized, but poorly understood feature of RBC populations within the same individual are their phenotypic heterogeneity. The granular characterization of phenotypic RBC variation in normative and disease states may allow us to identify the genetic determinants of red cell diseases and reveal novel therapeutic approaches for their treatment. Previously, we discovered diverse RNA transcripts in RBCs that has allowed us to dissect the phenotypic heterogeneity and malaria resistance of sickle red cells. However, these analyses failed to capture the heterogeneity found in RBC sub-populations. To overcome this limitation, we have performed single cell RNA-Seq to analyze the transcriptional heterogeneity of RBCs from three adult healthy donors which have been stored in the blood bank conditions and assayed at day 1 and day 15. The expression pattern clearly separated RBCs into seven distinct clusters that include one RBC cluster that expresses HBG2 and a small population of RBCs that express fetal hemoglobin (HbF) that we annotated as F cells. Almost all HBG2-expessing cells also express HBB, suggesting bi-allelic expression in single RBC from the HBG2/HBB loci, and we annotated another cluster as reticulocytes based on canonical gene expression. Additional RBC clusters were also annotated based on the enriched expression of NIX, ACVR2B and HEMGN, previously shown to be involved in erythropoiesis. Finally, we found the storage of RBC was associated with an increase in the ACVR2B and F-cell clusters. Collectively, these data indicate the power of single RBC RNA-Seq to capture and discover known and unexpected heterogeneity of RBC population.
| S-EPMC9065680 | biostudies-literature