Background

Expression Quantitative Trait Loci (eQTL) analysis, which identifies the impact of genomic variants on gene expression, has become a standard approach for interpreting the functions of disease-associated variants. Recently, advancements in single-cell transcriptome analysis (scRNA-seq) technology have fueled research into cell-type-specific eQTLs that act only on specific cell subtypes. This is because traditional bulk transcriptome analysis, which involves homogenizing entire tissues, struggles to capture the expression regulatory signals of rare immune cells or specific cell populations within the microenvironment.

However, significant statistical barriers existed for accurately calculating the correlation between genetic variants and gene expression at the single-cell level. Read depth varies across cells, and the proportion of zero counts (where expression is not detected) is extremely high. Researchers have primarily adopted a 'pseudobulk' approach, summing expression values of cells per individual to apply linear regression models. This approach significantly dilutes the heterogeneity and intracellular expression variation of individual cells, making it difficult to precisely correct for inter-individual relatedness or the genetic structure of multi-ethnic cohorts. Due to computational complexity, there were also clear limitations in testing the combined effects of low-frequency rare variants.

Key Findings

Researchers in the United States have reported SAIGE-QTL, a statistical genomics framework that directly handles single-cell read counts without post-processing, in the international journal Nature Genetics. This methodology is designed based on a Poisson Mixed Model to reflect the count data distribution unique to single-cell data without distortion.

By modeling the correlations among thousands of cells derived from individual donors and the genetic relatedness between samples as random effects, Type I errors (false positives) are suppressed. This structure directly controls covariates within the model, such as the total transcript count per cell or batch effects.

Compared to the pseudobulk approach, statistical power is significantly improved for rare subtypes with low cell counts or genes with low expression levels. The research team achieved computational efficiency capable of handling genomic data from tens of thousands of individuals and millions of single-cell data points by applying approximation estimation techniques validated in large-scale biobank-linked analyses.

The most prominent achievement is the implementation of rare variant set-based association testing in a single-cell environment. Because individual rare variants are low in frequency, it is nearly impossible to prove significance at the single-variant level. SAIGE-QTL directly linked burden tests and variance component tests—which group rare variants by gene or regulatory region units—to single-cell expression levels.

Significance and Outlook

SAIGE-QTL opens the path for genomic discovery spanning from common variants to rare variants at single-cell resolution. While previous studies were limited to analyzing common variants among single-nucleotide polymorphisms (SNPs), it is now possible to discern the effects of protein loss-of-function variants or rare non-coding variants in promoter regions on specific cell types.

It is expected to become a foundational tool for functional genomics research aimed at tracking the cellular pathways in which disease-susceptibility genes cease to function. It is particularly likely to be widely used in the analysis of complex diseases with highly restricted pathogenic cell types, such as neurodegenerative and autoimmune disorders.

However, the requirement for massive memory and high-performance computing resources remains a challenge. To comprehensively model spatial information between cells and the physical signals of neighboring cells in spatial transcriptomics data, further algorithmic expansion is required.