scholarly journals Bayesian cell-type deconvolution and gene expression inference reveals tumor-microenvironment interactions

2020 ◽  
Author(s):  
Tinyi Chu ◽  
Charles G. Danko
Keyword(s):  
Gene Expression ◽  
Cell Types ◽  
Malignant Cell ◽  
Malignant Cells ◽  
Cellular Composition ◽  
Rna Seq ◽  
Cell Type ◽  
Primary Glioblastoma ◽  
Multiple Tumor ◽  
Tumor Types

AbstractUnderstanding the complicated interactions between cells in their environment is a major challenge in genomics. Here we developed BayesPrism, a Bayesian method to jointly predict cellular composition and gene expression in each cell type, including heterogeneous malignant cells, from bulk RNA-seq using scRNA-seq as prior information. We conducted an integrative analysis of 1,412 bulk RNA-seq samples in primary glioblastoma, head and neck squamous cell carcinoma, and melanoma using single-cell datasets of 85 patients. We identified cell types correlated with clinical outcomes and explored spatial heterogeneity in tumor state and stromal composition. We refined subtypes using gene expression in malignant cells, after excluding confounding non-malignant cell types. Finally, we identified genes whose expression in malignant cells correlated with infiltration of macrophages, T cells, fibroblasts, and endothelial cells across multiple tumor types. Our work introduces a new lens that uses scRNA-seq to accurately infer cellular composition and expression in large cohorts of bulk data.

scholarly journals Integration of GWAS Summary Statistics and Gene Expression Reveals Target Cell Types Underlying Kidney Function Traits

2020 ◽  
Vol 31 (10) ◽  
pp. 2326-2340 ◽  
Author(s):  
Yong Li ◽  
Stefan Haug ◽  
Pascal Schlosser ◽  
Alexander Teumer ◽  
Adrienne Tin ◽  
...  
Keyword(s):  
Gene Expression ◽  
Kidney Function ◽  
Cell Types ◽  
Gene Set Enrichment Analysis ◽  
Genome Wide Association Studies ◽  
Summary Statistics ◽  
Rna Seq ◽  
Cell Type ◽  
Gene Set Enrichment ◽  
Gene Set

BackgroundGenetic variants identified in genome-wide association studies (GWAS) are often not specific enough to reveal complex underlying physiology. By integrating RNA-seq data and GWAS summary statistics, novel computational methods allow unbiased identification of trait-relevant tissues and cell types.MethodsThe CKDGen consortium provided GWAS summary data for eGFR, urinary albumin-creatinine ratio (UACR), BUN, and serum urate. Genotype-Tissue Expression Project (GTEx) RNA-seq data were used to construct the top 10% specifically expressed genes for each of 53 tissues followed by linkage disequilibrium (LD) score–based enrichment testing for each trait. Similar procedures were performed for five kidney single-cell RNA-seq datasets from humans and mice and for a microdissected tubule RNA-seq dataset from rat. Gene set enrichment analyses were also conducted for genes implicated in Mendelian kidney diseases.ResultsAcross 53 tissues, genes in kidney function–associated GWAS loci were enriched in kidney (P=9.1E-8 for eGFR; P=1.2E-5 for urate) and liver (P=6.8·10-5 for eGFR). In the kidney, proximal tubule was enriched in humans (P=8.5E-5 for eGFR; P=7.8E-6 for urate) and mice (P=0.0003 for eGFR; P=0.0002 for urate) and confirmed as the primary cell type in microdissected tubules and organoids. Gene set enrichment analysis supported this and showed enrichment of genes implicated in monogenic glomerular diseases in podocytes. A systematic approach generated a comprehensive list of GWAS genes prioritized by cell type–specific expression.ConclusionsIntegration of GWAS statistics of kidney function traits and gene expression data identified relevant tissues and cell types, as a basis for further mechanistic studies to understand GWAS loci.

scholarly journals Identifying gene expression programs of cell-type identity and cellular activity with single-cell RNA-Seq

eLife ◽  
2019 ◽  
Vol 8 ◽  
Author(s):  
Dylan Kotliar ◽  
Adrian Veres ◽  
M Aurel Nagy ◽  
Shervin Tabrizi ◽  
Eran Hodis ◽  
...  
Keyword(s):  
Gene Expression ◽  
Single Cell ◽  
Matrix Factorization ◽  
Cell Types ◽  
Environmental Cues ◽  
Rna Seq ◽  
Cell Type ◽  
Type Identity ◽  
Brain Organoid ◽  
Non Negative Matrix Factorization

Identifying gene expression programs underlying both cell-type identity and cellular activities (e.g. life-cycle processes, responses to environmental cues) is crucial for understanding the organization of cells and tissues. Although single-cell RNA-Seq (scRNA-Seq) can quantify transcripts in individual cells, each cell’s expression profile may be a mixture of both types of programs, making them difficult to disentangle. Here, we benchmark and enhance the use of matrix factorization to solve this problem. We show with simulations that a method we call consensus non-negative matrix factorization (cNMF) accurately infers identity and activity programs, including their relative contributions in each cell. To illustrate the insights this approach enables, we apply it to published brain organoid and visual cortex scRNA-Seq datasets; cNMF refines cell types and identifies both expected (e.g. cell cycle and hypoxia) and novel activity programs, including programs that may underlie a neurosecretory phenotype and synaptogenesis.

scholarly journals MarkerCount: A stable, count-based cell type identifier for single cell RNA-Seq experiments

2021 ◽  
Author(s):  
Hanbyeol Kim ◽  
Joongho Lee ◽  
Keunsoo Kang ◽  
Seokhyun Yoon
Keyword(s):  
Gene Expression ◽  
Single Cell ◽  
Cell Types ◽  
Batch Effect ◽  
Expression Level ◽  
Rna Seq ◽  
Cell Type ◽  
Stable Performance ◽  
Downstream Analysis

Abstract Cell type identification is a key step to downstream analysis of single cell RNA-seq experiments. Indispensible information for this is gene expression, which is used to cluster cells, train the model and set rejection thresholds. Problem is they are subject to batch effect arising from different platforms and preprocessing. We present MarkerCount, which uses the number of markers expressed regardless of their expression level to initially identify cell types and, then, reassign cell type in cluster-basis. MarkerCount works both in reference and marker-based mode, where the latter utilizes only the existing lists of markers, while the former required pre-annotated dataset to train the model. The performance was evaluated and compared with the existing identifiers, both marker and reference-based, that can be customized with publicly available datasets and marker DB. The results show that MarkerCount provides a stable performance when comparing with other reference-based and marker-based cell type identifiers.

scholarly journals Identifying Gene Expression Programs of Cell-type Identity and Cellular Activity with Single-Cell RNA-Seq

2018 ◽  
Author(s):  
Dylan Kotliar ◽  
Adrian Veres ◽  
M. Aurel Nagy ◽  
Shervin Tabrizi ◽  
Eran Hodis ◽  
...  
Keyword(s):  
Gene Expression ◽  
Single Cell ◽  
Matrix Factorization ◽  
Cell Types ◽  
Rna Seq ◽  
Cell Type ◽  
Relative Contribution ◽  
Neuronal Synapses ◽  
Type Identity ◽  
Brain Organoid

AbstractIdentifying gene expression programs underlying both cell-type identity and cellular activities (e.g. life-cycle processes, responses to environmental cues) is crucial for understanding the organization of cells and tissues. Although single-cell RNA-Seq (scRNA-Seq) can quantify transcripts in individual cells, each cell’s expression profile may be a mixture of both types of programs, making them difficult to disentangle. Here we illustrate and enhance the use of matrix factorization as a solution to this problem. We show with simulations that a method that we call consensus non-negative matrix factorization (cNMF) accurately infers identity and activity programs, including the relative contribution of programs in each cell. Applied to published brain organoid and visual cortex scRNA-Seq datasets, cNMF refines the hierarchy of cell types and identifies both expected (e.g. cell cycle and hypoxia) and intriguing novel activity programs. We propose that one of the novel programs may reflect a neurosecretory phenotype and a second may underlie the formation of neuronal synapses. We make cNMF available to the community and illustrate how this approach can provide key insights into gene expression variation within and between cell types.

scholarly journals A computational method for direct imputation of cell type-specific expression profiles and cellular compositions from bulk-tissue RNA-Seq in brain disorders

2020 ◽  
Author(s):  
Abolfazl Doostparast Torshizi ◽  
Jubao Duan ◽  
Kai Wang
Keyword(s):  
Gene Expression ◽  
Expression Profiles ◽  
Complex Diseases ◽  
Specific Gene ◽  
Cellular Composition ◽  
Rna Seq ◽  
Cell Type ◽  
Specific Expression ◽  
Cell Type Specific Expression ◽  
Cell Type Specific

AbstractThe importance of cell type-specific gene expression in disease-relevant tissues is increasingly recognized in genetic studies of complex diseases. However, the vast majority of gene expression studies are conducted on bulk tissues, necessitating computational approaches to infer biological insights on cell type-specific contribution to diseases. Several computational methods are available for cell type deconvolution (that is, inference of cellular composition) from bulk RNA-Seq data, but cannot impute cell type-specific expression profiles. We hypothesize that with external prior information such as single cell RNA-seq (scRNA-seq) and population-wide expression profiles, it can be a computationally tractable and identifiable to estimate both cellular composition and cell type-specific expression from bulk RNA-Seq data. Here we introduce CellR, which addresses cross-individual gene expression variations by employing genome-wide tissue-wise expression signatures from GTEx to adjust the weights of cell-specific gene markers. It then transforms the deconvolution problem into a linear programming model while taking into account inter/intra cellular correlations, and uses a multi-variate stochastic search algorithm to estimate the expression level of each gene in each cell type. Extensive analyses on several complex diseases such as schizophrenia, Alzheimer’s disease, Huntington’s disease, and type 2 diabetes validated efficiency of CellR, while revealing how specific cell types contribute to different diseases. We conducted numerical simulations on human cerebellum to generate pseudo-bulk RNA-seq data and demonstrated its efficiency in inferring cell-specific expression profiles. Moreover, we inferred cell-specific expression levels from bulk RNA-seq data on schizophrenia and computed differentially expressed genes within certain cell types. Using predicted gene expression profile on excitatory neurons, we were able to reproduce our recently published findings on TCF4 being a master regulator in schizophrenia and showed how this gene and its targets are enriched in excitatory neurons. In summary, CellR compares favorably (both accuracy and stability of inference) against competing approaches on inferring cellular composition from bulk RNA-seq data, but also allows direct imputation of cell type-specific gene expression, opening new doors to re-analyze gene expression data on bulk tissues in complex diseases.

scholarly journals Comprehensive characterization of tissue-specific chromatin accessibility in L2 Caenorhabditis elegans nematodes

2020 ◽  
Author(s):  
Timothy J. Durham ◽  
Riza M. Daza ◽  
Louis Gevirtzman ◽  
Darren A. Cusanovich ◽  
William Stafford Noble ◽  
...  
Keyword(s):  
Gene Expression ◽  
Single Cell ◽  
Expression Patterns ◽  
Cell Types ◽  
Chromatin Accessibility ◽  
Gene Expression Patterns ◽  
Rna Seq ◽  
Cell Type ◽  
Tissue Specific ◽  
C Elegans

AbstractRecently developed single cell technologies allow researchers to characterize cell states at ever greater resolution and scale. C. elegans is a particularly tractable system for studying development, and recent single cell RNA-seq studies characterized the gene expression patterns for nearly every cell type in the embryo and at the second larval stage (L2). Gene expression patterns are useful for learning about gene function and give insight into the biochemical state of different cell types; however, in order to understand these cell types, we must also determine how these gene expression levels are regulated. We present the first single cell ATAC-seq study in C. elegans. We collected data in L2 larvae to match the available single cell RNA-seq data set, and we identify tissue-specific chromatin accessibility patterns that align well with existing data, including the L2 single cell RNA-seq results. Using a novel implementation of the latent Dirichlet allocation algorithm, we leverage the single-cell resolution of the sci-ATAC-seq data to identify accessible loci at the level of individual cell types, providing new maps of putative cell type-specific gene regulatory sites, with promise for better understanding of cellular differentiation and gene regulation in the worm.

Cell type-specific endometrial transcriptome changes during initial recognition of pregnancy in the mare

2019 ◽  
Vol 31 (3) ◽  
pp. 496 ◽  
Author(s):  
Iside Scaravaggi ◽  
Nicole Borel ◽  
Rebekka Romer ◽  
Isabel Imboden ◽  
Susanne E. Ulbrich ◽  
...  
Keyword(s):  
Gene Expression ◽  
Spatial Information ◽  
Gene Expression Regulation ◽  
Cell Types ◽  
Rna Seq ◽  
Cell Type ◽  
Expression Of Genes ◽  
Prostaglandin F ◽  
Cell Type Specific ◽  
Gene Expression Studies

Previous endometrial gene expression studies during the time of conceptus migration did not provide final conclusions on the mechanisms of maternal recognition of pregnancy (MRP) in the mare. This called for a cell type-specific endometrial gene expression analysis in response to embryo signals to improve the understanding of gene expression regulation in the context of MRP. Laser capture microdissection was used to collect luminal epithelium (LE), glandular epithelium and stroma from endometrial biopsies from Day 12 of pregnancy and Day 12 of the oestrous cycle. RNA sequencing (RNA-Seq) showed greater expression differences between cell types than between pregnant and cyclic states; differences between the pregnant and cyclic states were mainly found in LE. Comparison with a previous RNA-Seq dataset for whole biopsy samples revealed the specific origin of gene expression differences. Furthermore, genes specifically differentially expressed (DE) in one cell type were found that were not detectable as DE in biopsies. Overall, this study revealed spatial information about endometrial gene expression during the phase of initial MRP. The conceptus induced changes in the expression of genes involved in blood vessel development, specific spatial regulation of the immune system, growth factors, regulation of prostaglandin synthesis, transport prostaglandin receptors, specifically prostaglandin F receptor (PTGFR) in the context of prevention of luteolysis.

scholarly journals CDSeqR: fast complete deconvolution for gene expression data from bulk tissues

2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Kai Kang ◽  
Caizhi Huang ◽  
Yuanyuan Li ◽  
David M. Umbach ◽  
Leping Li
Keyword(s):  
Gene Expression ◽  
Cell Types ◽  
Biological Tissues ◽  
Specific Gene ◽  
Specific Cell ◽  
Specific Information ◽  
Expression Data ◽  
Rna Seq ◽  
Cell Type ◽  
Cell Type Specific

Abstract Background Biological tissues consist of heterogenous populations of cells. Because gene expression patterns from bulk tissue samples reflect the contributions from all cells in the tissue, understanding the contribution of individual cell types to the overall gene expression in the tissue is fundamentally important. We recently developed a computational method, CDSeq, that can simultaneously estimate both sample-specific cell-type proportions and cell-type-specific gene expression profiles using only bulk RNA-Seq counts from multiple samples. Here we present an R implementation of CDSeq (CDSeqR) with significant performance improvement over the original implementation in MATLAB and an added new function to aid cell type annotation. The R package would be of interest for the broader R community. Result We developed a novel strategy to substantially improve computational efficiency in both speed and memory usage. In addition, we designed and implemented a new function for annotating the CDSeq estimated cell types using single-cell RNA sequencing (scRNA-seq) data. This function allows users to readily interpret and visualize the CDSeq estimated cell types. In addition, this new function further allows the users to annotate CDSeq-estimated cell types using marker genes. We carried out additional validations of the CDSeqR software using synthetic, real cell mixtures, and real bulk RNA-seq data from the Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) project. Conclusions The existing bulk RNA-seq repositories, such as TCGA and GTEx, provide enormous resources for better understanding changes in transcriptomics and human diseases. They are also potentially useful for studying cell–cell interactions in the tissue microenvironment. Bulk level analyses neglect tissue heterogeneity, however, and hinder investigation of a cell-type-specific expression. The CDSeqR package may aid in silico dissection of bulk expression data, enabling researchers to recover cell-type-specific information.

scholarly journals DSAVE: Detection of misclassified cells in single-cell RNA-Seq data

PLoS ONE ◽  
2020 ◽  
Vol 15 (12) ◽  
pp. e0243360
Author(s):  
Johan Gustafsson ◽  
Jonathan Robinson ◽  
Juan S. Inda-Díaz ◽  
Elias Björnson ◽  
Rebecka Jörnsten ◽  
...  
Keyword(s):  
Gene Expression ◽  
Single Cell ◽  
Cell Types ◽  
Rna Seq ◽  
Cell Type ◽  
Log Likelihood ◽  
Single Cell Rna Sequencing ◽  
Cell Transcriptome ◽  
Average Gene ◽  
Single Cell Transcriptome

Single-cell RNA sequencing has become a valuable tool for investigating cell types in complex tissues, where clustering of cells enables the identification and comparison of cell populations. Although many studies have sought to develop and compare different clustering approaches, a deeper investigation into the properties of the resulting populations is lacking. Specifically, the presence of misclassified cells can influence downstream analyses, highlighting the need to assess subpopulation purity and to detect such cells. We developed DSAVE (Down-SAmpling based Variation Estimation), a method to evaluate the purity of single-cell transcriptome clusters and to identify misclassified cells. The method utilizes down-sampling to eliminate differences in sampling noise and uses a log-likelihood based metric to help identify misclassified cells. In addition, DSAVE estimates the number of cells needed in a population to achieve a stable average gene expression profile within a certain gene expression range. We show that DSAVE can be used to find potentially misclassified cells that are not detectable by similar tools and reveal the cause of their divergence from the other cells, such as differing cell state or cell type. With the growing use of single-cell RNA-seq, we foresee that DSAVE will be an increasingly useful tool for comparing and purifying subpopulations in single-cell RNA-Seq datasets.

scholarly journals Single-nucleus and spatial transcriptomics of archival pancreatic cancer reveals multi-compartment reprogramming after neoadjuvant treatment

2020 ◽  
Author(s):  
William L. Hwang ◽  
Karthik A. Jagadeesh ◽  
Jimmy A. Guo ◽  
Hannah I. Hoffman ◽  
Payman Yadollahpour ◽  
...  
Keyword(s):  
Molecular Taxonomy ◽  
Neoadjuvant Treatment ◽  
Cell Types ◽  
Malignant Cell ◽  
Ductal Adenocarcinoma ◽  
Precision Oncology ◽  
Malignant Cells ◽  
Rna Seq ◽  
Cytotoxic Treatment ◽  
Single Nucleus

ABSTRACTPancreatic ductal adenocarcinoma (PDAC) remains a treatment-refractory disease. Characterizing PDAC by mRNA profiling remains particularly challenging. Previously identified bulk expression subtypes were influenced by contaminating stroma and have not yet informed clinical management, whereas single cell RNA-seq (scRNA-seq) of fresh tumors under-represented key cell types. Here, we developed a robust single-nucleus RNA-seq (snRNA-seq) technique for frozen archival PDAC specimens and used it to study both untreated tumors and those that received neoadjuvant chemotherapy and radiotherapy (CRT). Gene expression programs learned across untreated malignant cell and fibroblast profiles uncovered a clinically relevant molecular taxonomy with improved prognostic stratification compared to prior classifications. Moreover, in the increasingly-adopted neoadjuvant treatment context, there was a depletion of classical-like phenotypes in malignant cells in favor of basal-like phenotypes associated with TNF-NFkB and interferon signaling as well as the presence of novel acinar and neuroendocrine classical-like states, which may be more resilient to cytotoxic treatment. Spatially-resolved transcriptomics revealed an association between malignant cells expressing these basal-like programs and higher immune infiltration with increased lymphocytic content, whereas those exhibiting classical-like programs were linked to sparser macrophage-predominant microniches, perhaps pointing to susceptibility to distinct therapeutic strategies. Our refined molecular taxonomy and spatial resolution can help advance precision oncology in PDAC through informative stratification in clinical trials and insights into differential therapeutic targeting leveraging the immune system.

Sign in / Sign up

Export Citation Format

Share Document