Abstract:Gene expression levels can vary substantially across cells, even in a seemingly homogeneous cell population. Identifying the relationships between genetic variation and gene expression is critical for understanding the mechanisms of genome regulation. However, the genetic control of gene expression variability among the cells individuals has yet to be extensively examined. This is primarily due to the statistical challenges, such as the need for sufficiently powered cohorts and adjusting mean-variance dependence. Here, we introduce MEOTIVE (Mapping genetic Effects On inTra-Individual Variability of gene Expression), a novel statistical framework to identify genetic effects on the gene expression variability (sc-veQTL) accounting for the mean-variance dependence. Using single-cell RNA-seq data of 1.2 million peripheral blood mononuclear cells from 980 human donors, we identified 14 – 3,488 genes with significant sc-veQTLs (study-wide -value < 0.05) across different blood cell types, 2,103 of which were shared across more than one cell type. We further detected 55 SNP-gene pairs (in 34 unique genes) by directly linking genetic variations with gene expression dispersion (sc-deQTL) regardless of mean-variance dependence, and these genes were enriched in biological processes relevant to immune response and viral infection. An example is rs1131017 ( <9.08×10 ), a sc-veQTL in the 5’UTR of , which shows a ubiquitous dispersion effect across cell types, with higher dispersion levels associated with lower auto-immune disease risk, including rheumatoid arthritis and type 1 diabetes. Another example is , which is associated with antibacterial activity against bacterial species and was only detected with a monocyte-specific deQTL (rs1384) located at the 3’ UTR region ( =1.48×10 ) and replicated in an independent cohort. Our results demonstrate an efficient and robust statistical method to identify genetic effects on gene expression variability and how these associations and their involved pathways confer auto-immune disease risk. This analytical framework provides a new approach to unravelling the genetic regulation of gene expression at the single-cell resolution, advancing our understanding of complex biological processes.

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

Single-cell Eqtlgen Consortium: a Personalized Understanding of Disease

The Single-Cell Eqtlgen Consortium

Research progress on single-cell expression quantitative trait loci

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

Single-cell transcriptome-wide Mendelian randomization and colocalization analyses uncover cell-specific mechanisms in atherosclerotic cardiovascular disease

Sherlock: detecting gene-disease associations by matching patterns of expression QTL and GWAS.

Discovering single-cell eQTLs from scRNA-seq data only

scPrediXcan integrates advances in deep learning and single-cell data into a powerful cell-type–specific transcriptome-wide association study framework

Genetic variants associated with cell-type-specific intra-individual gene expression variability reveal new mechanisms of genome regulation

Integrating 3D Genomic and Epigenomic Data to Enhance Target Gene Discovery and Drug Repurposing in Transcriptome-Wide Association Studies

Enhancing Disease Risk Gene Discovery by Integrating Transcription Factor-Linked Trans-located Variants into Transcriptome-Wide Association Analyses

sc2GWAS: a comprehensive platform linking single cell and GWAS traits of human

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

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

Single-cell genomics improves the discovery of risk variants and genes of cardiac traits

The Power of Single-Cell RNA Sequencing in eQTL Discovery

Integrating Predicted Transcriptome From Multiple Tissues Improves Association Detection

Identification of genetic variants that impact gene co-expression relationships using large-scale single-cell data

Single-cell eQTL analysis of activated T cell subsets reveals activation and cell type–dependent effects of disease-risk variants

Systematic identification of trans eQTLs as putative drivers of known disease associations