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.

Quantification of escape from X chromosome inactivation with single-cell omics data reveals heterogeneity across cell types and tissues

Quantification of the escape from X chromosome inactivation with the million cell-scale human single-cell omics datasets reveals heterogeneity of escape across cell types and tissues

Inferring genes that escape X-Chromosome inactivation reveals important contribution of variable escape genes to sex-biased diseases

Escape from X Inactivation Varies in Mouse Tissues

Global Survey of Escape from X Inactivation by RNA-sequencing in Mouse.

Identification of Genes Escaping X Inactivation by Allelic Expression Analysis in a Novel Hybrid Mouse Model

Functional Inference of Gene Regulation Using Single-Cell Multi-Omics

Clustered transcripts that escape X inactivation at mouse XqD

Integrated analysis of Xist upregulation and X-chromosome inactivation with single-cell and single-allele resolution

Random X-chromosome inactivation dynamics in vivo by single-cell RNA sequencing

Escape from X Inactivation in Mice and Humans.

Escape from X inactivation is directly modulated by levels of Xist non-coding RNA

Single-cell RNA sequencing of peripheral blood links cell-type-specific regulation of splicing to autoimmune and inflammatory diseases

B cell-specific XIST complex enforces X-inactivation and restrains atypical B cells

Integrated analysis of Xist upregulation and gene silencing at the onset of random X-chromosome inactivation at high temporal and allelic resolution

Massively parallel single-cell chromatin landscapes of human immune cell development and intratumoral T cell exhaustion

A whole-organism landscape of X-inactivation in humans

Deconvolution Of Single-Cell Multi-Omics Layers Reveals Regulatory Heterogeneity

Genes that escape from X inactivation

Characterization of X Chromosome Inactivation Using Integrated Analysis of Whole-Exome and mRNA Sequencing

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