Abstract:Dissecting large bulk RNA-seq data into cell proportions and cell type specific expression profiles could significantly enhance our understanding of disease mechanisms at cell level and facilitate the identification of novel drug targets and development of efficient intervention strategies. In this study, we presented a scRNA-seq marker (gene) informed cell deconvolution and expression inference (CausalCellInfer) framework. CausalCellInfer employed causal inference principles to automatically identify a small set of critical marker genes from the reference scRNA-seq dataset. It then integrates deep neural networks with regularized matrix completion algorithms to deconvolute cell proportions and estimate cell type specific (CTS) expression profiles. Most importantly, we pioneered the application of the proposed framework to imputed expression data from large-scale genome-wide association studies (GWAS). We verified the efficacy of our proposed method by comparing it against existing state-of-the-art cell deconvolution methods, including CIBERSORTx, DWLS, Scaden, and MuSic across various real and pseudo-bulk samples. Furthermore, we proposed the use of a wide range of enrichment analyses to demonstrate the reliability of CausalCellInfer in estimating CTS profiles. Our proposed framework consistently outperformed existing methods with significantly higher concordance correlation coefficient (CCC), lower mean absolute error (MAE) and root mean square error (RMSE) across all testing real and pseudo bulk samples. Importantly, it also demonstrated superior computational efficiency compared to all benchmarked methods except MuSic. We also applied our trained models to deconvolute 4 tissue-specific expression cell proportions and estimate the corresponding CTS expression profiles, leveraging UK Biobank data (UKBB). We conducted a series of cellular-level analyses, including cell proportion association analysis, causal gene detection, transcriptome wide association analysis (TWAS), for 24 phenotypes in UKBB based on estimated cell composition and CTS expression profiles. Of note, the estimated proportions of various cell types were indicative of disease onset. For example, T2DM patients demonstrated a significant decrease in the proportions of alpha and beta cells in comparison to controls. We also tested their associations. Our method exhibited satisfactory positive predictive values (PPV) in uncovering differentially expressed (DE) genes for the majority of cell types. Encouragingly, most identified CTS causally relevant genes were found to be significantly enriched in target diseases or related pathophysiology. In conclusion, we presented a novel framework for inferring cell-type proportions and CTS expression, with novel applications to GWAS-imputed expression data from large-scale biobank program. Our work also shed light on how differential cell-type proportion and CTS expression may be associated with susceptibility to different diseases and their prognoses, bridging scRNA-seq and clinical phenotypes in large-scale biobank studies.

Cell-type deconvolution of bulk-blood RNA-seq reveals biological insights into neuropsychiatric disorders

Evaluating performance and applications of sample-wise cell deconvolution methods on human brain transcriptomic data

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

Bayesian estimation of cell type–specific gene expression with prior derived from single-cell data

Single-cell Eqtlgen Consortium: a Personalized Understanding of Disease

A unified framework for cell-type-specific eQTLs prioritization by integrating bulk and scRNA-seq data

Leveraging single-cell ATAC-seq and RNA-seq to identify disease-critical fetal and adult brain cell types

Cell-type specific inference from bulk RNA-sequencing data by integrating single cell reference profiles via EPIC-unmix

Clustering-independent estimation of cell abundances in bulk tissues using single-cell RNA-seq data

Gene expression imputation and cell-type deconvolution in human brain with spatiotemporal precision and its implications for brain-related disorders

Brain cell-type shifts in Alzheimer's disease, autism, and schizophrenia interrogated using methylomics and genetics

Identifying disease-critical cell types and cellular processes by integrating single-cell RNA-sequencing and human genetics

Connecting genomic results for psychiatric disorders to human brain cell types and regions reveals convergence with functional connectivity

Single-cell and spatial transcriptomics: deciphering brain complexity in health and disease

Network-based drug repurposing for psychiatric disorders using single-cell genomics

Brain expression quantitative trait locus and network analyses reveal downstream effects and putative drivers for brain-related diseases

Multi-omics analyses of cognitive traits and psychiatric disorders highlights brain-dependent mechanisms

Deconvolution from bulk gene expression by leveraging sample-wise and gene-wise similarities and single-cell RNA-Seq data

A Novel Bayesian Framework for Harmonizing Information Across Tissues and Studies to Increase Cell Type Deconvolution Accuracy

Characteristics of allelic gene expression in human brain cells from single-cell RNA-seq data analysis

Estimating cell compositions and cell-type-specific expressions from GWAS data using invariant causal prediction, deep learning and regularized matrix completion