1
|
Nguyen H, Nguyen H, Tran D, Draghici S, Nguyen T. Fourteen years of cellular deconvolution: methodology, applications, technical evaluation and outstanding challenges. Nucleic Acids Res 2024; 52:4761-4783. [PMID: 38619038 PMCID: PMC11109966 DOI: 10.1093/nar/gkae267] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2023] [Revised: 03/01/2024] [Accepted: 04/02/2024] [Indexed: 04/16/2024] Open
Abstract
Single-cell RNA sequencing (scRNA-Seq) is a recent technology that allows for the measurement of the expression of all genes in each individual cell contained in a sample. Information at the single-cell level has been shown to be extremely useful in many areas. However, performing single-cell experiments is expensive. Although cellular deconvolution cannot provide the same comprehensive information as single-cell experiments, it can extract cell-type information from bulk RNA data, and therefore it allows researchers to conduct studies at cell-type resolution from existing bulk datasets. For these reasons, a great effort has been made to develop such methods for cellular deconvolution. The large number of methods available, the requirement of coding skills, inadequate documentation, and lack of performance assessment all make it extremely difficult for life scientists to choose a suitable method for their experiment. This paper aims to fill this gap by providing a comprehensive review of 53 deconvolution methods regarding their methodology, applications, performance, and outstanding challenges. More importantly, the article presents a benchmarking of all these 53 methods using 283 cell types from 30 tissues of 63 individuals. We also provide an R package named DeconBenchmark that allows readers to execute and benchmark the reviewed methods (https://github.com/tinnlab/DeconBenchmark).
Collapse
Affiliation(s)
- Hung Nguyen
- Department of Computer Science and Software Engineering, Auburn University, Auburn, AL, USA
| | - Ha Nguyen
- Department of Computer Science and Software Engineering, Auburn University, Auburn, AL, USA
| | - Duc Tran
- Department of Medicine, Washington University School of Medicine, St. Louis, MO, USA
| | - Sorin Draghici
- Department of Computer Science, Wayne State University, Detroit, MI, USA
- Advaita Bioinformatics, Ann Arbor, MI, USA
| | - Tin Nguyen
- Department of Computer Science and Software Engineering, Auburn University, Auburn, AL, USA
| |
Collapse
|
2
|
Ouyang Z, Bourgeois-Tchir N, Lyashenko E, Cundiff PE, Cullen PF, Challa R, Li K, Zhang X, Casey F, Engle SJ, Zhang B, Zavodszky MI. Characterizing the composition of iPSC derived cells from bulk transcriptomics data with CellMap. Sci Rep 2022; 12:17394. [PMID: 36253414 PMCID: PMC9576729 DOI: 10.1038/s41598-022-22115-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2021] [Accepted: 10/10/2022] [Indexed: 01/10/2023] Open
Abstract
Induced pluripotent stem cell (iPSC) derived cell types are increasingly employed as in vitro model systems for drug discovery. For these studies to be meaningful, it is important to understand the reproducibility of the iPSC-derived cultures and their similarity to equivalent endogenous cell types. Single-cell and single-nucleus RNA sequencing (RNA-seq) are useful to gain such understanding, but they are expensive and time consuming, while bulk RNA-seq data can be generated quicker and at lower cost. In silico cell type decomposition is an efficient, inexpensive, and convenient alternative that can leverage bulk RNA-seq to derive more fine-grained information about these cultures. We developed CellMap, a computational tool that derives cell type profiles from publicly available single-cell and single-nucleus datasets to infer cell types in bulk RNA-seq data from iPSC-derived cell lines.
Collapse
Affiliation(s)
- Zhengyu Ouyang
- BioInfoRx, Inc., 510 Charmany Dr, Suite 275A, Madison, WI 53719 USA
| | - Nathanael Bourgeois-Tchir
- grid.417832.b0000 0004 0384 8146Translational Biology, Research and Development, Biogen, Inc., 225 Binney St, Cambridge, MA 02142 USA
| | - Eugenia Lyashenko
- grid.417832.b0000 0004 0384 8146Translational Biology, Research and Development, Biogen, Inc., 225 Binney St, Cambridge, MA 02142 USA ,grid.417555.70000 0000 8814 392XPresent Address: Genomic Medicine Unit, Sanofi, 225 2nd Ave, Waltham, MA 02451 USA
| | - Paige E. Cundiff
- grid.417832.b0000 0004 0384 8146Translational Biology, Research and Development, Biogen, Inc., 225 Binney St, Cambridge, MA 02142 USA
| | - Patrick F. Cullen
- grid.417832.b0000 0004 0384 8146Translational Biology, Research and Development, Biogen, Inc., 225 Binney St, Cambridge, MA 02142 USA
| | - Ravi Challa
- grid.417832.b0000 0004 0384 8146Translational Biology, Research and Development, Biogen, Inc., 225 Binney St, Cambridge, MA 02142 USA
| | - Kejie Li
- grid.417832.b0000 0004 0384 8146Translational Biology, Research and Development, Biogen, Inc., 225 Binney St, Cambridge, MA 02142 USA
| | - Xinmin Zhang
- BioInfoRx, Inc., 510 Charmany Dr, Suite 275A, Madison, WI 53719 USA
| | - Fergal Casey
- grid.417832.b0000 0004 0384 8146Translational Biology, Research and Development, Biogen, Inc., 225 Binney St, Cambridge, MA 02142 USA
| | - Sandra J. Engle
- grid.417832.b0000 0004 0384 8146Translational Biology, Research and Development, Biogen, Inc., 225 Binney St, Cambridge, MA 02142 USA
| | - Baohong Zhang
- grid.417832.b0000 0004 0384 8146Translational Biology, Research and Development, Biogen, Inc., 225 Binney St, Cambridge, MA 02142 USA
| | - Maria I. Zavodszky
- grid.417832.b0000 0004 0384 8146Translational Biology, Research and Development, Biogen, Inc., 225 Binney St, Cambridge, MA 02142 USA
| |
Collapse
|
3
|
Bhattacharya A, Hamilton AM, Troester MA, Love MI. DeCompress: tissue compartment deconvolution of targeted mRNA expression panels using compressed sensing. Nucleic Acids Res 2021; 49:e48. [PMID: 33524140 PMCID: PMC8096278 DOI: 10.1093/nar/gkab031] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2020] [Revised: 12/21/2020] [Accepted: 01/12/2021] [Indexed: 12/13/2022] Open
Abstract
Targeted mRNA expression panels, measuring up to 800 genes, are used in academic and clinical settings due to low cost and high sensitivity for archived samples. Most samples assayed on targeted panels originate from bulk tissue comprised of many cell types, and cell-type heterogeneity confounds biological signals. Reference-free methods are used when cell-type-specific expression references are unavailable, but limited feature spaces render implementation challenging in targeted panels. Here, we present DeCompress, a semi-reference-free deconvolution method for targeted panels. DeCompress leverages a reference RNA-seq or microarray dataset from similar tissue to expand the feature space of targeted panels using compressed sensing. Ensemble reference-free deconvolution is performed on this artificially expanded dataset to estimate cell-type proportions and gene signatures. In simulated mixtures, four public cell line mixtures, and a targeted panel (1199 samples; 406 genes) from the Carolina Breast Cancer Study, DeCompress recapitulates cell-type proportions with less error than reference-free methods and finds biologically relevant compartments. We integrate compartment estimates into cis-eQTL mapping in breast cancer, identifying a tumor-specific cis-eQTL for CCR3 (C-C Motif Chemokine Receptor 3) at a risk locus. DeCompress improves upon reference-free methods without requiring expression profiles from pure cell populations, with applications in genomic analyses and clinical settings.
Collapse
Affiliation(s)
- Arjun Bhattacharya
- Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California-Los Angeles, Los Angeles, CA 90095, USA
| | - Alina M Hamilton
- Department of Pathology and Laboratory Medicine, University of North Carolina-Chapel Hill, Chapel Hill, NC 27516, USA
| | - Melissa A Troester
- Department of Pathology and Laboratory Medicine, University of North Carolina-Chapel Hill, Chapel Hill, NC 27516, USA
- Department of Epidemiology, University of North Carolina-Chapel Hill, Chapel Hill, NC 27516, USA
| | - Michael I Love
- Department of Biostatistics, University of North Carolina-Chapel Hill, Chapel Hill, NC 27516, USA
- Department of Genetics, University of North Carolina-Chapel Hill, Chapel Hill, NC 27516, USA
| |
Collapse
|
4
|
Wang X, Allen M, Li S, Quicksall ZS, Patel TA, Carnwath TP, Reddy JS, Carrasquillo MM, Lincoln SJ, Nguyen TT, Malphrus KG, Dickson DW, Crook JE, Asmann YW, Ertekin-Taner N. Deciphering cellular transcriptional alterations in Alzheimer's disease brains. Mol Neurodegener 2020; 15:38. [PMID: 32660529 PMCID: PMC7359236 DOI: 10.1186/s13024-020-00392-6] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2020] [Accepted: 06/27/2020] [Indexed: 02/06/2023] Open
Abstract
Large-scale brain bulk-RNAseq studies identified molecular pathways implicated in Alzheimer's disease (AD), however these findings can be confounded by cellular composition changes in bulk-tissue. To identify cell intrinsic gene expression alterations of individual cell types, we designed a bioinformatics pipeline and analyzed three AD and control bulk-RNAseq datasets of temporal and dorsolateral prefrontal cortex from 685 brain samples. We detected cell-proportion changes in AD brains that are robustly replicable across the three independently assessed cohorts. We applied three different algorithms including our in-house algorithm to identify cell intrinsic differentially expressed genes in individual cell types (CI-DEGs). We assessed the performance of all algorithms by comparison to single nucleus RNAseq data. We identified consensus CI-DEGs that are common to multiple brain regions. Despite significant overlap between consensus CI-DEGs and bulk-DEGs, many CI-DEGs were absent from bulk-DEGs. Consensus CI-DEGs and their enriched GO terms include genes and pathways previously implicated in AD or neurodegeneration, as well as novel ones. We demonstrated that the detection of CI-DEGs through computational deconvolution methods is promising and highlight remaining challenges. These findings provide novel insights into cell-intrinsic transcriptional changes of individual cell types in AD and may refine discovery and modeling of molecular targets that drive this complex disease.
Collapse
Affiliation(s)
- Xue Wang
- Department of Health Sciences Research, Mayo Clinic Florida, Jacksonville, FL, USA.
| | - Mariet Allen
- Department of Neuroscience, Mayo Clinic Florida, Jacksonville, FL, USA
| | - Shaoyu Li
- Department of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC, USA
| | - Zachary S Quicksall
- Department of Health Sciences Research, Mayo Clinic Florida, Jacksonville, FL, USA
| | - Tulsi A Patel
- Department of Neuroscience, Mayo Clinic Florida, Jacksonville, FL, USA
| | - Troy P Carnwath
- Department of Neuroscience, Mayo Clinic Florida, Jacksonville, FL, USA
| | - Joseph S Reddy
- Department of Health Sciences Research, Mayo Clinic Florida, Jacksonville, FL, USA
| | | | - Sarah J Lincoln
- Department of Neuroscience, Mayo Clinic Florida, Jacksonville, FL, USA
| | - Thuy T Nguyen
- Department of Neuroscience, Mayo Clinic Florida, Jacksonville, FL, USA
| | | | - Dennis W Dickson
- Department of Neuroscience, Mayo Clinic Florida, Jacksonville, FL, USA
| | - Julia E Crook
- Department of Health Sciences Research, Mayo Clinic Florida, Jacksonville, FL, USA
| | - Yan W Asmann
- Department of Health Sciences Research, Mayo Clinic Florida, Jacksonville, FL, USA
| | - Nilüfer Ertekin-Taner
- Department of Neuroscience, Mayo Clinic Florida, Jacksonville, FL, USA.
- Department of Neurology, Mayo Clinic Florida, Jacksonville, FL, USA.
| |
Collapse
|
5
|
Ng'oma E, Williams-Simon PA, Rahman A, King EG. Diverse biological processes coordinate the transcriptional response to nutritional changes in a Drosophila melanogaster multiparent population. BMC Genomics 2020; 21:84. [PMID: 31992183 PMCID: PMC6988245 DOI: 10.1186/s12864-020-6467-6] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2019] [Accepted: 01/08/2020] [Indexed: 12/19/2022] Open
Abstract
Background Environmental variation in the amount of resources available to populations challenge individuals to optimize the allocation of those resources to key fitness functions. This coordination of resource allocation relative to resource availability is commonly attributed to key nutrient sensing gene pathways in laboratory model organisms, chiefly the insulin/TOR signaling pathway. However, the genetic basis of diet-induced variation in gene expression is less clear. Results To describe the natural genetic variation underlying nutrient-dependent differences, we used an outbred panel derived from a multiparental population, the Drosophila Synthetic Population Resource. We analyzed RNA sequence data from multiple female tissue samples dissected from flies reared in three nutritional conditions: high sugar (HS), dietary restriction (DR), and control (C) diets. A large proportion of genes in the experiment (19.6% or 2471 genes) were significantly differentially expressed for the effect of diet, and 7.8% (978 genes) for the effect of the interaction between diet and tissue type (LRT, Padj. < 0.05). Interestingly, we observed similar patterns of gene expression relative to the C diet, in the DR and HS treated flies, a response likely reflecting diet component ratios. Hierarchical clustering identified 21 robust gene modules showing intra-modularly similar patterns of expression across diets, all of which were highly significant for diet or diet-tissue interaction effects (FDR Padj. < 0.05). Gene set enrichment analysis for different diet-tissue combinations revealed a diverse set of pathways and gene ontology (GO) terms (two-sample t-test, FDR < 0.05). GO analysis on individual co-expressed modules likewise showed a large number of terms encompassing many cellular and nuclear processes (Fisher exact test, Padj. < 0.01). Although a handful of genes in the IIS/TOR pathway including Ilp5, Rheb, and Sirt2 showed significant elevation in expression, many key genes such as InR, chico, most insulin peptide genes, and the nutrient-sensing pathways were not observed. Conclusions Our results suggest that a more diverse network of pathways and gene networks mediate the diet response in our population. These results have important implications for future studies focusing on diet responses in natural populations.
Collapse
Affiliation(s)
- E Ng'oma
- University of Missouri, 401 Tucker Hall, Columbia, MO, 65211, USA.
| | | | - A Rahman
- University of Missouri, 401 Tucker Hall, Columbia, MO, 65211, USA
| | - E G King
- University of Missouri, 401 Tucker Hall, Columbia, MO, 65211, USA
| |
Collapse
|
6
|
Berens ME, Sood A, Barnholtz-Sloan JS, Graf JF, Cho S, Kim S, Kiefer J, Byron SA, Halperin RF, Nasser S, Adkins J, Cuyugan L, Devine K, Ostrom Q, Couce M, Wolansky L, McDonough E, Schyberg S, Dinn S, Sloan AE, Prados M, Phillips JJ, Nelson SJ, Liang WS, Al-Kofahi Y, Rusu M, Zavodszky MI, Ginty F. Multiscale, multimodal analysis of tumor heterogeneity in IDH1 mutant vs wild-type diffuse gliomas. PLoS One 2019; 14:e0219724. [PMID: 31881020 PMCID: PMC6934292 DOI: 10.1371/journal.pone.0219724] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2019] [Accepted: 11/12/2019] [Indexed: 12/31/2022] Open
Abstract
Glioma is recognized to be a highly heterogeneous CNS malignancy, whose diverse cellular composition and cellular interactions have not been well characterized. To gain new clinical- and biological-insights into the genetically-bifurcated IDH1 mutant (mt) vs wildtype (wt) forms of glioma, we integrated data from protein, genomic and MR imaging from 20 treatment-naïve glioma cases and 16 recurrent GBM cases. Multiplexed immunofluorescence (MxIF) was used to generate single cell data for 43 protein markers representing all cancer hallmarks, Genomic sequencing (exome and RNA (normal and tumor) and magnetic resonance imaging (MRI) quantitative features (protocols were T1-post, FLAIR and ADC) from whole tumor, peritumoral edema and enhancing core vs equivalent normal region were also collected from patients. Based on MxIF analysis, 85,767 cells (glioma cases) and 56,304 cells (GBM cases) were used to generate cell-level data for 24 biomarkers. K-means clustering was used to generate 7 distinct groups of cells with divergent biomarker profiles and deconvolution was used to assign RNA data into three classes. Spatial and molecular heterogeneity metrics were generated for the cell data. All features were compared between IDH mt and IDHwt patients and were finally combined to provide a holistic/integrated comparison. Protein expression by hallmark was generally lower in the IDHmt vs wt patients. Molecular and spatial heterogeneity scores for angiogenesis and cell invasion also differed between IDHmt and wt gliomas irrespective of prior treatment and tumor grade; these differences also persisted in the MR imaging features of peritumoral edema and contrast enhancement volumes. A coherent picture of enhanced angiogenesis in IDHwt tumors was derived from multiple platforms (genomic, proteomic and imaging) and scales from individual proteins to cell clusters and heterogeneity, as well as bulk tumor RNA and imaging features. Longer overall survival for IDH1mt glioma patients may reflect mutation-driven alterations in cellular, molecular, and spatial heterogeneity which manifest in discernable radiological manifestations.
Collapse
Affiliation(s)
- Michael E. Berens
- Translational Genomics Research Institute, Phoenix, AZ, United States of America
- * E-mail: (MEB); (AS); (FG)
| | - Anup Sood
- GE Research Center, Niskayuna, NY, United States of America
- * E-mail: (MEB); (AS); (FG)
| | - Jill S. Barnholtz-Sloan
- Department of Population and Quantitative Health Sciences and Case Comprehensive Cancer Center, Case Western Reserve University School of Medicine, Cleveland, OH, United States of America
| | - John F. Graf
- GE Research Center, Niskayuna, NY, United States of America
| | - Sanghee Cho
- GE Research Center, Niskayuna, NY, United States of America
| | - Seungchan Kim
- Department of Electrical and Computer Engineering, Roy G. Perry College of Engineering, Prairie View A&M University, Prairie View, TX, United States of America
| | - Jeffrey Kiefer
- Translational Genomics Research Institute, Phoenix, AZ, United States of America
| | - Sara A. Byron
- Translational Genomics Research Institute, Phoenix, AZ, United States of America
| | - Rebecca F. Halperin
- Translational Genomics Research Institute, Phoenix, AZ, United States of America
| | - Sara Nasser
- Translational Genomics Research Institute, Phoenix, AZ, United States of America
| | - Jonathan Adkins
- Translational Genomics Research Institute, Phoenix, AZ, United States of America
| | - Lori Cuyugan
- Translational Genomics Research Institute, Phoenix, AZ, United States of America
| | - Karen Devine
- Department of Population and Quantitative Health Sciences and Case Comprehensive Cancer Center, Case Western Reserve University School of Medicine, Cleveland, OH, United States of America
| | - Quinn Ostrom
- Department of Population and Quantitative Health Sciences and Case Comprehensive Cancer Center, Case Western Reserve University School of Medicine, Cleveland, OH, United States of America
| | - Marta Couce
- Department of Population and Quantitative Health Sciences and Case Comprehensive Cancer Center, Case Western Reserve University School of Medicine, Cleveland, OH, United States of America
| | - Leo Wolansky
- Department of Population and Quantitative Health Sciences and Case Comprehensive Cancer Center, Case Western Reserve University School of Medicine, Cleveland, OH, United States of America
| | | | | | - Sean Dinn
- GE Research Center, Niskayuna, NY, United States of America
| | - Andrew E. Sloan
- Department of Neurosurgery, University Hospitals-Seidman Cancer Center, Cleveland, OH, United States of America
| | - Michael Prados
- Department of Neurological Surgery, Helen Diller Cancer Center, University of California San Francisco, San Francisco, CA, United States of America
| | - Joanna J. Phillips
- Department of Neurological Surgery, Helen Diller Cancer Center, University of California San Francisco, San Francisco, CA, United States of America
| | - Sarah J. Nelson
- Department of Radiology and Biomedical Imaging, University of California, San Francisco, CA, United States of America
| | - Winnie S. Liang
- Translational Genomics Research Institute, Phoenix, AZ, United States of America
| | | | - Mirabela Rusu
- GE Research Center, Niskayuna, NY, United States of America
| | | | - Fiona Ginty
- GE Research Center, Niskayuna, NY, United States of America
- * E-mail: (MEB); (AS); (FG)
| |
Collapse
|