scholarly journals Deep autoencoder enables interpretable tissue-adaptive deconvolution and cell-type-specific gene analysis

2021 ◽  
Author(s):  
Yanshuo Chen ◽  
Yixuan Wang ◽  
Yuelong Chen ◽  
Yumeng Wei ◽  
Yunxiang Li ◽  
...  
Keyword(s):  
Gene Expression ◽  
Single Cell ◽  
Clinical Data ◽  
Gene Set Enrichment Analysis ◽  
Specific Gene ◽  
Rna Seq ◽  
Cell Type ◽  
Specific Gene Expression ◽  
Wide Range ◽  
Cell Type Specific

AbstractSingle-cell RNA-seq has become a powerful tool for researchers to study biologically significant characteristics at explicitly high resolution, but its application on emerging data is currently limited by its intrinsic techniques. Here, we introduce TAPE, a deep learning method that connects bulk RNA-seq and single-cell RNA-seq to balance the demands of big data and precision. By taking advantage of constructing an interpretable decoder and training under a unique scheme, TAPE can predict cell-type fractions and cell-type-specific gene expression tissue-adaptively. Compared with existing methods on several benchmarking datasets, TAPE is more accurate (up to 40% performnace improvement on the real bulk data) and faster than the previous methods. It is sensitive enough to provide biologically meaningful predictions. For example, only TAPE can predict the tendency of increasing monocytes-to-lymphocytes (MLR) ratio in COVID-19 patients from mild to serious symptoms, whose estimated indices are consistent with laboratory data. More importantly, through the analysis of clinical data, TAPE shows its ability to predict cell-type-specific gene expression profiles with biological significance. Combining with single-sample gene set enrichment analysis (ssGSEA), TAPE also provides valuable clues for people to investigate the immune response in different virus-infected patients. We believe that TAPE will enable and accelerate the precise analysis of high-throughput clinical data in a wide range.

scholarly journals Sources of Variation in Cell-Type RNA-Seq Profiles

2020 ◽  
Author(s):  
Johan Gustafsson ◽  
Felix Held ◽  
Jonathan Robinson ◽  
Elias Björnson ◽  
Rebecka Jörnsten ◽  
...  
Keyword(s):  
Gene Expression ◽  
Single Cell ◽  
Expression Profiles ◽  
Gene Expression Profiles ◽  
Specific Gene ◽  
Rna Seq ◽  
Cell Type ◽  
Specific Gene Expression ◽  
Cell Type Specific ◽  
Technical Factors

Abstract Background Cell-type specific gene expression profiles are needed for many computational methods operating on bulk RNA-Seq samples, such as deconvolution of cell-type fractions and digital cytometry. However, the gene expression profile of a cell type can vary substantially due to both technical factors and biological differences in cell state and surroundings, reducing the efficacy of such methods. Here, we investigated which factors contribute most to this variation. Results We evaluated different normalization methods, quantified the magnitude of variation introduced by different sources, and examined the differences between UMI-based single-cell RNA-Seq and bulk RNA-Seq. We applied methods such as random forest regression to a collection of publicly available bulk and single-cell RNA-Seq datasets containing B and T cells, and found that the technical variation across laboratories is of the same magnitude as the biological variation across cell types. Tissue of origin and cell subtype are less important but still substantial factors, while the difference between individuals is relatively small. We also show that much of the differences between UMI-based single-cell and bulk RNA-Seq methods can be explained by the number of read duplicates per mRNA molecule in the single-cell sample.Conclusions Our work shows the importance of either matching or correcting for technical factors when creating cell-type specific gene expression profiles that are to be used together with bulk samples.

scholarly journals Sources of Variation in Cell-Type RNA-Seq Profiles

2020 ◽  
Author(s):  
Johan Gustafsson ◽  
Felix Held ◽  
Jonathan Robinson ◽  
Elias Björnson ◽  
Rebecka Jörnsten ◽  
...  
Keyword(s):  
Gene Expression ◽  
Single Cell ◽  
Expression Profiles ◽  
Gene Expression Profiles ◽  
Specific Gene ◽  
Rna Seq ◽  
Cell Type ◽  
Specific Gene Expression ◽  
Cell Type Specific ◽  
Technical Factors

Abstract Cell-type specific gene expression profiles are needed for many computational methods operating on bulk RNA-Seq samples, such as deconvolution of cell-type fractions and digital cytometry. However, the gene expression profile of a cell type can vary substantially due to both technical factors and biological differences in cell state and surroundings, reducing the efficacy of such methods. Here, we investigated which factors contribute most to this variation. We evaluated different normalization methods, quantified the variance explained by different factors, evaluated the effect on deconvolution of cell type fractions, and examined the differences between UMI-based single-cell RNA-Seq and bulk RNA-Seq. We investigated a collection of publicly available bulk and single-cell RNA-Seq datasets containing B and T cells, and found that the technical variation across laboratories is substantial, even for genes specifically selected for deconvolution, and has a confounding effect on deconvolution. Tissue of origin is also a substantial factor, highlighting the challenge of applying cell type profiles derived from blood on mixtures from other tissues. We also show that much of the differences between UMI-based single-cell and bulk RNA-Seq methods can be explained by the number of read duplicates per mRNA molecule in the single-cell sample. Our work shows the importance of either matching or correcting for technical factors when creating cell-type specific gene expression profiles that are to be used together with bulk samples.

scholarly journals Dynamic changes in cis-regulatory occupancy by Six1 and its cooperative interactions with distinct cofactors drive lineage-specific gene expression programs during progressive differentiation of the auditory sensory epithelium

2020 ◽  
Vol 48 (6) ◽  
pp. 2880-2896 ◽  
Author(s):  
Jun Li ◽  
Ting Zhang ◽  
Aarthi Ramakrishnan ◽  
Bernd Fritzsch ◽  
Jinshu Xu ◽  
...  
Keyword(s):  
Gene Expression ◽  
Regulatory Elements ◽  
Sensory Epithelium ◽  
Hair Bundle ◽  
Specific Gene ◽  
Cell Type ◽  
Specific Gene Expression ◽  
Wide Range ◽  
Cell Type Specific ◽  
Lineage Specific Gene

Abstract The transcription factor Six1 is essential for induction of sensory cell fate and formation of auditory sensory epithelium, but how it activates gene expression programs to generate distinct cell-types remains unknown. Here, we perform genome-wide characterization of Six1 binding at different stages of auditory sensory epithelium development and find that Six1-binding to cis-regulatory elements changes dramatically at cell-state transitions. Intriguingly, Six1 pre-occupies enhancers of cell-type-specific regulators and effectors before their expression. We demonstrate in-vivo cell-type-specific activity of Six1-bound novel enhancers of Pbx1, Fgf8, Dusp6, Vangl2, the hair-cell master regulator Atoh1 and a cascade of Atoh1’s downstream factors, including Pou4f3 and Gfi1. A subset of Six1-bound sites carry consensus-sequences for its downstream factors, including Atoh1, Gfi1, Pou4f3, Gata3 and Pbx1, all of which physically interact with Six1. Motif analysis identifies RFX/X-box as one of the most significantly enriched motifs in Six1-bound sites, and we demonstrate that Six1-RFX proteins cooperatively regulate gene expression through binding to SIX:RFX-motifs. Six1 targets a wide range of hair-bundle regulators and late Six1 deletion disrupts hair-bundle polarity. This study provides a mechanistic understanding of how Six1 cooperates with distinct cofactors in feedforward loops to control lineage-specific gene expression programs during progressive differentiation of the auditory sensory epithelium.

scholarly journals Improved estimation of cell type-specific gene expression through deconvolution of bulk tissues with matrix completion

2021 ◽  
Author(s):  
Weixu Wang ◽  
Jun Yao ◽  
Yi Wang ◽  
Chao Zhang ◽  
Wei Tao ◽  
...  
Keyword(s):  
Gene Expression ◽  
Single Cells ◽  
Matrix Completion ◽  
Superior Performance ◽  
Specific Gene ◽  
Rna Seq ◽  
Cell Type ◽  
Specific Gene Expression ◽  
Large Patient ◽  
Cell Type Specific

Cell type-specific gene expression (CSE) brings novel biological insights into both physiological and pathological processes compared with bulk tissue gene expression. Although fluorescence-activated cell sorting (FACS) and single-cell RNA sequencing (scRNA-seq) are two widely used techniques to detect gene expression in a cell type-specific manner, the constraints of cost and labor force make it impractical as a routine on large patient cohorts. Here, we present ENIGMA, an algorithm that deconvolutes bulk RNA-seq into cell type-specific expression matrices and cell type fraction matrices without the need of physical sorting or sequencing of single cells. ENIGMA used cell type signature matrix generated from either FACS RNA-seq or scRNA-seq as reference, and applied matrix completion technique to achieve fast and accurate deconvolution. We demonstrated the superior performance of ENIGMA to previously published algorithms (TCA, bMIND and CIBERSORTx) while requiring much less running time on both simulated and realistic datasets. To prove its value in biological discovery, we applied ENIGMA to bulk RNA-seq from arthritis patients and revealed a pseudo-differentiation trajectory that could reflect monocyte to macrophage transition. We also applied ENIGMA to bulk RNA-seq data of pancreatic islet tissue from type 2 diabetes (T2D) patients and discovered a beta cell-specific gene co-expression module related to senescence and apoptosis that possibly contributed to the pathogenesis of T2D. Together, ENIGMA provides a new framework to improve the CSE estimation by integrating FACS RNA-seq and scRNA-seq with tissue bulk RNA-seq data, and will extend our understandings about cell heterogeneity on population level with no need for experimental tissue disaggregation.

scholarly journals An Efficient and Flexible Method for Deconvoluting Bulk RNA-Seq Data with Single-Cell RNA-Seq Data

Cells ◽  
2019 ◽  
Vol 8 (10) ◽  
pp. 1161 ◽  
Author(s):  
Xifang Sun ◽  
Shiquan Sun ◽  
Sheng Yang
Keyword(s):  
Gene Expression ◽  
Single Cell ◽  
Rna Sequencing ◽  
Specific Gene ◽  
Rna Seq ◽  
Cell Type ◽  
Disease Etiology ◽  
Expression Levels ◽  
Cell Type Specific ◽  
Gene Expression Levels

Estimating cell type compositions for complex diseases is an important step to investigate the cellular heterogeneity for understanding disease etiology and potentially facilitate early disease diagnosis and prevention. Here, we developed a computationally statistical method, referring to Multi-Omics Matrix Factorization (MOMF), to estimate the cell-type compositions of bulk RNA sequencing (RNA-seq) data by leveraging cell type-specific gene expression levels from single-cell RNA sequencing (scRNA-seq) data. MOMF not only directly models the count nature of gene expression data, but also effectively accounts for the uncertainty of cell type-specific mean gene expression levels. We demonstrate the benefits of MOMF through three real data applications, i.e., Glioblastomas (GBM), colorectal cancer (CRC) and type II diabetes (T2D) studies. MOMF is able to accurately estimate disease-related cell type proportions, i.e., oligodendrocyte progenitor cells and macrophage cells, which are strongly associated with the survival of GBM and CRC, respectively.

scholarly journals A novel computational complete deconvolution method using RNA-seq data

2018 ◽  
Author(s):  
Kai Kang ◽  
Qian Meng ◽  
Igor Shats ◽  
David M. Umbach ◽  
Melissa Li ◽  
...  
Keyword(s):  
Gene Expression ◽  
Expression Profiles ◽  
Gene Expression Profiles ◽  
Specific Gene ◽  
Rna Seq ◽  
Cell Type ◽  
Specific Gene Expression ◽  
Cell Type Composition ◽  
Type Composition ◽  
Cell Type Specific

AbstractThe cell type composition of many biological tissues varies widely across samples. Such sample heterogeneity hampers efforts to probe the role of each cell type in the tissue microenvironment. Current approaches that address this issue have drawbacks. Cell sorting or single-cell based experimental techniques disrupt in situ interactions and alter physiological status of cells in tissues. Computational methods are flexible and promising; but they often estimate either sample-specific proportions of each cell type or cell-type-specific gene expression profiles, not both, by requiring the other as input. We introduce a computational Complete Deconvolution method that can estimate both sample-specific proportions of each cell type and cell-type-specific gene expression profiles simultaneously using bulk RNA-Seq data only (CDSeq). We assessed our method’s performance using several synthetic and experimental mixtures of varied but known cell type composition and compared its performance to the performance of two state-of-the art deconvolution methods on the same mixtures. The results showed CDSeq can estimate both sample-specific proportions of each component cell type and cell-typespecificgene expression profiles with high accuracy. CDSeq holds promise for computationally deciphering complex mixtures of cell types, each with differing expression profiles, using RNA-seq data measured in bulk tissue (MATLAB code is available at https://github.com/kkang7/CDSeq_011).

scholarly journals A novel method for predicting cell abundance based on single-cell RNA-seq data

2021 ◽  
Vol 22 (S9) ◽  
Author(s):  
Jiajie Peng ◽  
Lu Han ◽  
Xuequn Shang
Keyword(s):  
Gene Expression ◽  
Least Squares ◽  
Single Cell ◽  
Specific Gene ◽  
Rna Seq ◽  
Cell Type ◽  
Gene Expression Matrix ◽  
Cell Type Specific ◽  
Novel Method ◽  
Signature Matrix

Abstract Background It is important to understand the composition of cell type and its proportion in intact tissues, as changes in certain cell types are the underlying cause of disease in humans. Although compositions of cell type and ratios can be obtained by single-cell sequencing, single-cell sequencing is currently expensive and cannot be applied in clinical studies involving a large number of subjects. Therefore, it is useful to apply the bulk RNA-Seq dataset and the single-cell RNA dataset to deconvolute and obtain the cell type composition in the tissue. Results By analyzing the existing cell population prediction methods, we found that most of the existing methods need the cell-type-specific gene expression profile as the input of the signature matrix. However, in real applications, it is not always possible to find an available signature matrix. To solve this problem, we proposed a novel method, named DCap, to predict cell abundance. DCap is a deconvolution method based on non-negative least squares. DCap considers the weight resulting from measurement noise of bulk RNA-seq and calculation error of single-cell RNA-seq data, during the calculation process of non-negative least squares and performs the weighted iterative calculation based on least squares. By weighting the bulk tissue gene expression matrix and single-cell gene expression matrix, DCap minimizes the measurement error of bulk RNA-Seq and also reduces errors resulting from differences in the number of expressed genes in the same type of cells in different samples. Evaluation test shows that DCap performs better in cell type abundance prediction than existing methods. Conclusion DCap solves the deconvolution problem using weighted non-negative least squares to predict cell type abundance in tissues. DCap has better prediction results and does not need to prepare a signature matrix that gives the cell-type-specific gene expression profile in advance. By using DCap, we can better study the changes in cell proportion in diseased tissues and provide more information on the follow-up treatment of diseases.

Sign in / Sign up

Export Citation Format

Share Document