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.

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

Single-cell Eqtlgen Consortium: a Personalized Understanding of Disease

The Single-Cell Eqtlgen Consortium

A computational method for cell type-specific expression quantitative trait loci mapping using bulk RNA-seq data

Research progress on single-cell expression quantitative trait loci

Integrating single cell expression quantitative trait loci summary statistics to understand complex trait risk genes

Integrative modeling of eQTLs and cis-regulatory elements suggest mechanisms underlying cell type specificity of eQTLs

Detecting cell-type-specific allelic expression imbalance by integrative analysis of bulk and single-cell RNA sequencing data

Discovering single-cell eQTLs from scRNA-seq data only

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

A statistical framework for joint eQTL analysis in multiple tissues

Disentangling associations between complex traits and cell types with

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

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

Efficient and accurate mixed model association tool for single-cell eQTL analysis

The Power of Single-Cell RNA Sequencing in eQTL Discovery

Cepip: Context-Dependent Epigenomic Weighting for Prioritization of Regulatory Variants and Disease-Associated Genes

SingleQ: a comprehensive database of single-cell expression quantitative trait loci (sc-eQTLs) cross human tissues

Integration of expression QTLs with fine mapping via SuSiE

Structural Basis of Differential Gene Expression at eQTLs Loci from High-Resolution Ensemble Models of 3D Single-Cell Chromatin Conformations

Joint eQTL mapping and inference of gene regulatory network improves power of detecting both cis - and trans -eQTLs