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.

Gene expression variability across cells and species shapes innate immunity

Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells

Epigenetic Variation Impacts Individual Differences in the Transcriptional Response to Influenza Infection

Dissecting human population variation in single-cell responses to SARS-CoV-2

Single-cell Sequencing Reveals the Evolution of Immune Molecules Across Multiple Vertebrate Species.

Coordinated variation in the immune systems and microbiomes of healthy humans is linked to tonic interferon states

Quantifying dynamic pro-inflammatory gene expression and heterogeneity in single macrophage cells

Epigenetic programming of monocyte-to-macrophage differentiation and trained innate immunity

Gene expression variability across cells and species shapes the relationship between renal resident macrophages and infiltrated macrophages

Divergent cytokine and transcriptional signatures control functional T follicular helper cell heterogeneity

Conservation and divergence in Toll-like receptor 4-regulated gene expression in primary human versus mouse macrophages

Molecular mechanisms governing circulating immune cell heterogeneity across different species revealed by single‐cell sequencing

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

Interactions between immune cell types facilitate the evolution of immune traits

Multi-organ single-cell transcriptomics of immune cells uncovered organ-specific gene expression and functions

Behavioural immune landscapes of inflammation

Single-cell RNA sequencing reveals induction of distinct trained-immunity programs in human monocytes

Impact of Genetic Polymorphisms on Human Immune Cell Gene Expression

Demographic and genetic factors influence the abundance of infiltrating immune cells in human tissues

Shaping Variation in the Human Immune System

Single-cell RNA-sequencing of peripheral blood mononuclear cells reveals widespread, context-specific gene expression regulation upon pathogenic exposure