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.

Measurement of selective constraint on human gene expression

Evidence for Weak Selective Constraint on Human Gene Expression.

Epistasis And Destabilizing Mutations Shape Gene Expression Variability In Humans Via Distinct Modes Of Action

Evolution of dosage-sensitive genes by tissue-restricted expression changes

Aberrant Gene Expression in Humans

Haplotype-aware modeling of cis-regulatory effects highlights the gaps remaining in eQTL data

Strength of purifying selection on the amino-acid sequence is associated with the amount of non-additive variance in gene expression.

Allele-specific transcription factor binding across human brain regions offers mechanistic insight into eQTLs

Stabilizing selection and adaptation shape cis and trans gene expression variation in C. elegans

Probing the limits of cis-acting gene regulation using a model of allelic imbalance quantitative trait loci

Genetic and selective constraints on the optimization of gene product diversity

Cis- and Trans-regulatory Effects on Gene Expression in a Natural Population of Drosophila melanogaster

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

Complex Genetic Interactions Underlying Expression Differences Between Drosophila Races: Analysis Of Chromosome Substitutions

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

Additive, epistatic, and environmental effects through the lens of expression variability QTL in a twin cohort.

Evolutionary Dynamics of Regulatory Changes Underlying Gene Expression Divergence among Saccharomyces Species

The evolution of complex gene regulation by low specificity binding sites

Functional dynamic genetic effects on gene regulation are specific to particular cell types and environmental conditions

Analyses of Allele-Specific Gene Expression in Highly Divergent Mouse Crosses Identifies Pervasive Allelic Imbalance

Unique Signatures of Highly Constrained Genes Across Publicly Available Genomic Databases