Clonealign: statistical integration of independent single-cell RNA and DNA sequencing data from human cancers.
Ontology highlight
ABSTRACT:
SUBMITTER: Campbell KR
PROVIDER: S-EPMC6417140 | biostudies-literature | 2019 Mar
REPOSITORIES: biostudies-literature
Publications
clonealign: statistical integration of independent single-cell RNA and DNA sequencing data from human cancers.
Campbell Kieran R KR Steif Adi A Laks Emma E Zahn Hans H Lai Daniel D McPherson Andrew A Farahani Hossein H Kabeer Farhia F O'Flanagan Ciara C Biele Justina J Brimhall Jazmine J Wang Beixi B Walters Pascale P Bouchard-Côté Alexandre A Aparicio Samuel S Shah Sohrab P SP
Genome biology 20190312 1
Measuring gene expression of tumor clones at single-cell resolution links functional consequences to somatic alterations. Without scalable methods to simultaneously assay DNA and RNA from the same single cell, parallel single-cell DNA and RNA measurements from independent cell populations must be mapped for genome-transcriptome association. We present clonealign, which assigns gene expression states to cancer clones using single-cell RNA and DNA sequencing independently sampled from a heterogene ...[more]
Similar Datasets
Project description:Single-cell RNA-sequencing (scRNA-seq) is a recent high-throughput genomic technology used to study the expression dynamics of genes at single-cell level. Analyzing the scRNA-seq data in presence of biological confounding factors including dropout events is a challenging task. Thus, this article presents a novel statistical approach for various analyses of the scRNA-seq Unique Molecular Identifier (UMI) counts data. The various analyses include modeling and fitting of observed UMI data, cell type detection, estimation of cell capture rates, estimation of gene specific model parameters, estimation of the sample mean and sample variance of the genes, etc. Besides, the developed approach is able to perform differential expression, and other downstream analyses that consider the molecular capture process in scRNA-seq data modeling. Here, the external spike-ins data can also be used in the approach for better results. The unique feature of the method is that it considers the biological process that leads to severe dropout events in modeling the observed UMI counts of genes. • The differential expression analysis of observed scRNA-seq UMI counts data is performed after adjustment for cell capture rates. • The statistical approach performs downstream differential zero inflation analysis, classification of influential genes, and selection of top marker genes. • Cell auxiliaries including cell clusters and other cell variables (e.g., cell cycle, cell phase) are used to remove unwanted variation to perform statistical tests reliably.
| S-EPMC8720898 | biostudies-literature
Project description:Due to the development of next-generation RNA sequencing (NGS) technologies, there has been tremendous progress in research involving determining the role of genomics, transcriptomics and epigenomics in complex biological systems. However, scientists have realized that information obtained using earlier technology, frequently called 'bulk RNA-seq' data, provides information averaged across all the cells present in a tissue. Relatively newly developed single cell (scRNA-seq) technology allows us to provide transcriptomic information at a single-cell resolution. Nevertheless, these high-resolution data have their own complex natures and demand novel statistical data analysis methods to provide effective and highly accurate results on complex biological systems. In this review, we cover many such recently developed statistical methods for researchers wanting to pursue scRNA-seq statistical and computational research as well as scientific research about these existing methods and free software tools available for their generated data. This review is certainly not exhaustive due to page limitations. We have tried to cover the popular methods starting from quality control to the downstream analysis of finding differentially expressed genes and concluding with a brief description of network analysis.
| S-EPMC9400796 | biostudies-literature
Project description:Recent advances in technology have enabled the measurement of RNA levels for individual cells. Compared to traditional tissue-level bulk RNA-seq data, single cell sequencing yields valuable insights about gene expression profiles for different cell types, which is potentially critical for understanding many complex human diseases. However, developing quantitative tools for such data remains challenging because of high levels of technical noise, especially the "dropout" events. A "dropout" happens when the RNA for a gene fails to be amplified prior to sequencing, producing a "false" zero in the observed data. In this paper, we propose a Unified RNA-Sequencing Model (URSM) for both single cell and bulk RNA-seq data, formulated as a hierarchical model. URSM borrows the strength from both data sources and carefully models the dropouts in single cell data, leading to a more accurate estimation of cell type specific gene expression profile. In addition, URSM naturally provides inference on the dropout entries in single cell data that need to be imputed for downstream analyses, as well as the mixing proportions of different cell types in bulk samples. We adopt an empirical Bayes' approach, where parameters are estimated using the EM algorithm and approximate inference is obtained by Gibbs sampling. Simulation results illustrate that URSM outperforms existing approaches both in correcting for dropouts in single cell data, as well as in deconvolving bulk samples. We also demonstrate an application to gene expression data on fetal brains, where our model successfully imputes the dropout genes and reveals cell type specific expression patterns.
| S-EPMC6114100 | biostudies-literature
Project description:In complex tissues containing cells that are difficult to dissociate, single-nucleus RNA-sequencing (snRNA-seq) has become the preferred experimental technology over single-cell RNA-sequencing (scRNA-seq) to measure gene expression. To accurately model these data in downstream analyses, previous work has shown that droplet-based scRNA-seq data are not zero-inflated, but whether droplet-based snRNA-seq data follow the same probability distributions has not been systematically evaluated. Using pseudonegative control data from nuclei in mouse cortex sequenced with the 10x Genomics Chromium system and mouse kidney sequenced with the DropSeq system, we found that droplet-based snRNA-seq data follow a negative binomial distribution, suggesting that parametric statistical models applied to scRNA-seq are transferable to snRNA-seq. Furthermore, we found that the quantification choices in adapting quantification mapping strategies from scRNA-seq to snRNA-seq can play a significant role in downstream analyses and biological interpretation. In particular, reference transcriptomes that do not include intronic regions result in significantly smaller library sizes and incongruous cell type classifications. We also confirmed the presence of a gene length bias in snRNA-seq data, which we show is present in both exonic and intronic reads, and investigate potential causes for the bias.
| S-EPMC11247185 | biostudies-literature
Project description:The growing number of available single-cell gene expression datasets from different species creates opportunities to explore evolutionary relationships between cell types across species. Cross-species integration of single-cell RNA-sequencing data has been particularly informative in this context. However, in order to do so robustly it is essential to have rigorous benchmarking and appropriate guidelines to ensure that integration results truly reflect biology. Here, we benchmark 28 combinations of gene homology mapping methods and data integration algorithms in a variety of biological settings. We examine the capability of each strategy to perform species-mixing of known homologous cell types and to preserve biological heterogeneity using 9 established metrics. We also develop a new biology conservation metric to address the maintenance of cell type distinguishability. Overall, scANVI, scVI and SeuratV4 methods achieve a balance between species-mixing and biology conservation. For evolutionarily distant species, including in-paralogs is beneficial. SAMap outperforms when integrating whole-body atlases between species with challenging gene homology annotation. We provide our freely available cross-species integration and assessment pipeline to help analyse new data and develop new algorithms.
| S-EPMC10576752 | biostudies-literature
Project description:BackgroundSingle-cell RNA sequencing (scRNA-seq) methods have been advantageous for quantifying cell-to-cell variation by profiling the transcriptomes of individual cells. For scRNA-seq data, variability in gene expression reflects the degree of variation in gene expression from one cell to another. Analyses that focus on cell-cell variability therefore are useful for going beyond changes based on average expression and, instead, identifying genes with homogeneous expression versus those that vary widely from cell to cell.ResultsWe present a novel statistical framework, scShapes, for identifying differential distributions in single-cell RNA-sequencing data using generalized linear models. Most approaches for differential gene expression detect shifts in the mean value. However, as single-cell data are driven by overdispersion and dropouts, moving beyond means and using distributions that can handle excess zeros is critical. scShapes quantifies gene-specific cell-to-cell variability by testing for differences in the expression distribution while flexibly adjusting for covariates if required. We demonstrate that scShapes identifies subtle variations that are independent of altered mean expression and detects biologically relevant genes that were not discovered through standard approaches.ConclusionsThis analysis also draws attention to genes that switch distribution shapes from a unimodal distribution to a zero-inflated distribution and raises open questions about the plausible biological mechanisms that may give rise to this, such as transcriptional bursting. Overall, the results from scShapes help to expand our understanding of the role that gene expression plays in the transcriptional regulation of a specific perturbation or cellular phenotype. Our framework scShapes is incorporated into a Bioconductor R package (https://www.bioconductor.org/packages/release/bioc/html/scShapes.html).
| S-EPMC9871437 | biostudies-literature
Project description:The cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway has emerged as a critical innate immune pathway that could virtually impact nearly all aspects of tumorigenesis including colorectal cancer. This work aimed to develop and validate molecular subtypes related to cGAS-STING pathways for colorectal cancer using Bulk RNA-seq and single-cell RNA-seq (scRNA-seq) data. Bulk RNA-seq data were acquired from The Cancer Genome Atlas dataset (training dataset) and Gene Expression Omnibus dataset (validation dataset). Univariate COX survival analysis was utilized to identify prognostic differentially expressed genes (DEGs) from 6 immune pathways related to cGAS-STING. ConsensusClusterPlus package was used to classify different subtypes based on DEGs. scRNA-seq data were used to validate differences in immune status between different subtypes. Two clusters with distinct prognosis were identified based on 27 DEGs. The six cGAS-STING-related pathways had different levels of significance between the two clusters. Clust1 had most number of amplified CNVs and clust2 had the most number of loss CNVs. TP53 was the top mutated gene of which missense mutations contributed the most of single-nucleotide variants. Immune score of clust1 was higher than that in clust2, as reflected in macrophages, T cells, and natural killer cells. Three unfavorable genes and 31 protection factors were screened between the two clusters in three datasets. ScRNA-seq data analysis demonstrated that macrophages were more enriched in clust1, and tumor cells and immune cells had close interaction. We classified two distinct subtypes with different prognosis, mutation landscape, and immune characteristics.
| S-EPMC10290593 | biostudies-literature
Project description:Adductomics studies are used for the detection and characterization of various chemical modifications (adducts) of nucleic acids and proteins. The advancements in liquid chromatography coupled with high-resolution tandem mass spectrometry (HRMS/MS) have resulted in efficient methods for qualitative and quantitative adductomics. We developed an HRMS-based method for the simultaneous analysis of RNA and DNA adducts in a single run and demonstrated its application using Baltic amphipods, useful sentinels of environmental disturbances, as test organisms. The novelty of this method is screening for RNA and DNA adducts by a single injection on an Orbitrap HRMS instrument using full scan and data-independent acquisition. The MS raw files were processed with an open-source program, nLossFinder, to identify and distinguish RNA and DNA adducts based on the characteristic neutral loss of ribonucleosides and 2'-deoxyribonucleosides, respectively. In the amphipods, in addition to the nearly 150 putative DNA adducts characterized earlier, we detected 60 putative RNA adducts. For the structural identification of the detected RNA adducts, the MODOMICS database was used. The identified RNA adducts included simple mono- and dimethylation and other larger functional groups on different ribonucleosides and deaminated product inosine. However, 54 of these RNA adducts are not yet structurally identified, and further work on their characterization may uncover new layers of information related to the transcriptome and help understand their biological significance. Considering the susceptibility of nucleic acids to environmental factors, including pollutants, the developed multi-adductomics methodology with further advancement has the potential to provide biomarkers for diagnostics of pollution effects in biota.
| S-EPMC10523582 | biostudies-literature
Project description:MotivationRecombinant DNA technology is widely used for different applications in biology, medicine and bio-technology. Viral transduction and plasmid transfection are among the most frequently used techniques to generate recombinant cell lines. Many of these methods result in the random integration of the plasmid into the host genome. Rapid identification of the integration sites is highly desirable in order to characterize these engineered cell lines.ResultsWe developed detectIS: a pipeline specifically designed to identify genomic integration sites of exogenous DNA, either a plasmid containing one or more transgenes or a virus. The pipeline is based on a Nextflow workflow combined with a Singularity image containing all the necessary software, ensuring high reproducibility and scalability of the analysis. We tested it on simulated datasets and RNA-seq data from a human sample infected with Hepatitis B virus. Comparisons with other state of the art tools show that our method can identify the integration site in different recombinant cell lines, with accurate results, lower computational demand and shorter execution times.Availability and implementationThe Nextflow workflow, the Singularity image and a test dataset are available at https://github.com/AstraZeneca/detectIS.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC9502153 | biostudies-literature
Project description:Large-scale comprehensive single-cell experiments are often resource-intensive and require the involvement of many laboratories and/or taking measurements at various times. This inevitably leads to batch effects, and systematic variations in the data that might occur due to different technology platforms, reagent lots, or handling personnel. Such technical differences confound biological variations of interest and need to be corrected during the data integration process. Data integration is a challenging task due to the overlapping of biological and technical factors, which makes it difficult to distinguish their individual contribution to the overall observed effect. Moreover, the choice of integration method may impact the downstream analyses, including searching for differentially expressed genes. From the existing data integration methods, we selected only those that return the full expression matrix. We evaluated six methods in terms of their influence on the performance of differential gene expression analysis in two single-cell datasets with the same biological study design that differ only in the way the measurement was done: one dataset manifests strong batch effects due to the measurements of each sample at a different time. Integrated data were visualized using the UMAP method. The evaluation was done both on individual gene level using parametric and non-parametric approaches for finding differentially expressed genes and on gene set level using gene set enrichment analysis. As an evaluation metric, we used two correlation coefficients, Pearson and Spearman, of the obtained test statistics between reference, test, and corrected studies. Visual comparison of UMAP plots highlighted ComBat-seq, limma, and MNN, which reduced batch effects and preserved differences between biological conditions. Most of the tested methods changed the data distribution after integration, which negatively impacts the use of parametric methods for the analysis. Two algorithms, MNN and Scanorama, gave very poor results in terms of differential analysis on gene and gene set levels. Finally, we highlight ComBat-seq as it led to the highest correlation of test statistics between reference and corrected dataset among others. Moreover, it does not distort the original distribution of gene expression data, so it can be used in all types of downstream analyses.
| S-EPMC9663917 | biostudies-literature