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.

Trans-eQTL hotspots shape complex traits by modulating cellular states

Cross-phenotype meta-analysis reveals large-scale trans-eQTLs mediating patterns of transcriptional co-regulation

Distant regulatory effects of genetic variation in multiple human tissues

Transcriptional Regulatory Pattern in Yeast Revealed Through Expression Quantitative Trait Locus Mapping

Protein-protein interactions shape trans-regulatory impact of genetic variation on protein expression and complex traits

The genetic architecture of cell-type-specific cis-regulation

The Importance of Regulatory Network Structure for Complex Trait Heritability and Evolution

Transcripts with high distal heritability mediate genetic effects on complex metabolic traits

Mutational sources of trans-regulatory variation affecting gene expression in Saccharomyces cerevisiae

Distant Eqtls and Non-coding Sequences Play Critical Roles in Regulating Gene Expression and Quantitative Trait Variation in Maize.

Large-scale cis- and trans-eQTL analyses identify thousands of genetic loci and polygenic scores that regulate blood gene expression

-eQTL mapping in gene sets identifies network effects of genetic variants

Stitching Together Multiple Data Dimensions Reveals Interacting Metabolomic and Transcriptomic Networks That Modulate Cell Regulation

Gene expression quantitative trait locus analysis of 16 000 barley genes reveals a complex pattern of genome-wide transcriptional regulation.

Single-cell eQTL mapping in yeast reveals a tradeoff between growth and reproduction

The impact of cell type and context-dependent regulatory variants on human immune traits

Uncovering genetic regulatory network divergence between duplicate genes using yeast eQTL landscape.

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

Modeling Regulatory Network Topology Improves Genome-Wide Analyses of Complex Human Traits

Mendelian randomization integrating GWAS and eQTL data reveals genetic determinants of complex and clinical traits

Shared and distinct molecular effects of regulatory genetic variants provide insight into mechanisms of distal enhancer-promoter communication