Skip to content

NotablePeer-reviewed

Efficient and accurate mixed model association tool for single-cell eQTL analysis

In 1.2 million immune cells, a single-cell model found 48.8% more eGenes than a standard pooled pipeline

ByWei Zhou, Anna S. E. Cuomo, Linfeng Hu, Angli Xue, Masahiro Kanai, Yijia Christiana Liu and 2 others

Broad Institute · University of Helsinki · Harvard University · Massachusetts General Hospital and 4 more

2nd of 36 papers scored this week on Frontier.

What it is

Genetic effects on gene expression often vary across cell types, which has motivated mapping expression quantitative trait loci (eQTLs, DNA variants linked to how much a gene is expressed) in single-cell RNA sequencing data. Existing approaches generally aggregate each person's expression across cells (the pseudobulk approach), use Gaussian models poorly suited to read counts, or do not scale. The authors present SAIGE-QTL, a Poisson mixed model that directly models single-cell read counts, enabling calibrated and computationally scalable eQTL mapping while accounting for donor-level and cell-level covariates and unequal numbers of cells per donor. In simulations and three single-cell datasets, it improved statistical efficiency and the discovery of genes with significant eQTLs (eGenes) relative to pseudobulk methods, and produced results consistent with existing single-cell models at substantially lower computational cost. In 1.2 million immune cells from the OneK1K dataset (982 individuals, 14 cell types), it identified 48.8% more eGenes than pseudobulk analysis (17,218 against 11,569, counted per cell type). It also enabled genome-wide trans-eQTL mapping, for variants far from the gene they affect, and set-based tests of rare variants.

Fig. 1 of the paper: Performance evaluation of SAIGE-QTL.
Fig. 1 from the paperPerformance evaluation of SAIGE-QTL.Wei Zhou et al., Nature Genetics (2026), under CC BY 4.0, shown unaltered.

Why it matters

In the same immune-cell data, SAIGE-QTL found 48.8% more eGenes than the standard pseudobulk pipeline it was compared with. Part of that gap comes from the gene-level test: given the same ACAT-V combination, the pseudobulk tool TensorQTL found 12,777 eGenes instead of 11,569. SAIGE-QTL also runs set-based tests of rare variants, which the authors say are better powered than single-variant tests for low-frequency variants: 5,541 rare-variant eGenes across the 14 cell types, 2,317 unique genes, of which 483 (21%) remained significant after conditioning on the same genes' common eQTLs, and 413 trans-eQTLs for 390 genes in three cell types at genome-wide significance. Against glmer, which fits the same Poisson mixed model, it cut computation time 244-fold in a small benchmark in one B-cell type, with highly concordant results.

Underlined numbers link to their source. Every metric and quoted figure is listed under Sources and data below.

Filed underGenetic Associations and Epidemiology, Single-cell and spatial transcriptomics, Genetic Mapping and Diversity in Plants and Animals

Scored 2026-11-09How scoring works

The Frontier Score, in full

How this entry's score is built, term by term. Open any pillar to see each sub-metric's raw value, how it maps to 0 to 100, and a link to its source (OpenAlex for most sub-metrics that combine several sources, NIH iCite where its ratio is part of field-normalized impact).
Evidence 87Impact 67Novelty 78
The three pillars as one shape: a result strong on every axis fills the triangle.
How the Frontier Score is calculated for this breakthrough
TermScoreWeightPoints
EvidenceHow established and verifiable the result is.87× 0.2521.8
ImpactHow much the result matters, judged field-relative and by current momentum.67× 0.2516.7
NoveltyHow genuinely new the result is, and how fast-moving.78× 0.5039.0
Weighted pillars77.4
Frontier recency0 months since publication× 1.00
IntegrityNot retracted× 1
Frontier Score77
Confidence: Early signal. Early signal: little citation history yet, so impact is an estimate and may shift as the work is taken up.Early signal: little citation history yet, so impact is an estimate and may shift as the work is taken up.

Every sub-metric, traced to its source

Evidence87
Evidence sub-metrics
Sub-metricRaw value0 to 100WeightPointsSource
Peer review & venueWhether the work has cleared peer review, and how selective its venue is. Peer review is necessary but not sufficient: before citations accrue, a fresh result in one of the five general flagships (Nature, Science, Cell, NEJM, the Lancet) is a stronger evidence signal than one in a selective specialist or society journal, which in turn outranks one in a legitimate but very high-volume mega-journal, and any of them outranks a preprint that has not been reviewed at all.Maps to 0 to 100: Peer-reviewed article graded by venue tier: flagship journal (Nature, Science, Cell, NEJM, the Lancet) 100; selective journal (every other journal the index admits: the specialist journals of those families, JAMA, PNAS, Physical Review Letters and X) 90; high-volume mega-journal (Nature Communications, Science Advances) 78; venue not recorded 88. Other records: review 90, book 70, dataset 60, preprint 45, an article that has not been peer reviewed 45, other 55.Peer-reviewed, selective journal900.3027.0OpenAlex type / primary_location / source
CorroborationIndependent corroboration that the result is on the record and reaching beyond a single source, read as the BREADTH of independent records rather than the size of the citation pile. Counts how many of the five independent citation indices (OpenAlex, Crossref, Semantic Scholar, OpenCitations, Europe PMC) report a citation count for the work at all, zero included, plus orthogonal, non-citation lines of corroboration: citation in Wikipedia, registered clinical trials or clinical citations, released datasets or software or a repository copy, public funding on record, and technical-community discussion. Citation MAGNITUDE is deliberately judged under Impact, not here, so Evidence and Impact measure genuinely different things instead of both rewarding the same citation pile twice.Maps to 0 to 100: Starts at 28 and rises by 9 for each of the five independent citation indices that returns a citation count for the work (zero counts; OpenCitations returns zero even for a DOI it does not hold) and by 6 for each orthogonal non-citation line (Wikipedia, clinical trials or clinical citations, open datasets, software or a repository copy, public funding, community discussion), up to the top of the 0-100 scale.3 of 5 citation indices report a count, 2 independent non-citation lines670.3020.1OpenAlex, Crossref, Semantic Scholar, OpenCitations, Europe PMC, DataCite, NIH RePORTER, Wikipedia, Hacker News count of indices that report a count for the work + orthogonal reach
VerifiabilityHow openly the result can be read, reused, and reproduced. Rewards open access, a Creative Commons BY or CC0 license, a green repository copy anyone can archive, openly minable full text, and released datasets or software a reader can actually run. Cross-checked across OpenAlex, Unpaywall, Europe PMC, and DataCite.Maps to 0 to 100: Closed access 45; open access 88; any Creative Commons BY license (including the NC, ND and SA variants) or CC0 100. A repository copy raises the floor to 85, openly minable full text to 90, and released data or software to 92.Open access (hybrid), CC-BY, repository copy1000.2020.0OpenAlex, Unpaywall, Europe PMC, DataCite is_oa / license / has_repository_copy / open datasets
IntegrityRetraction and post-publication correction status. A retracted result is not evidence of anything and collapses the whole Frontier Score; a correction or expression of concern is a softer flag. A correction the publisher labels as its own error (a Publisher Correction) is not a flag on the work. Read across OpenAlex and Crossref (both update-to and updated-by notices).Maps to 0 to 100: Clean 100 (also when the only notice is a Publisher Correction); a correction or expression of concern 75; retracted 0. Flagged if OpenAlex OR Crossref records the notice.No retraction or concern on record1000.2020.0OpenAlex, Crossref is_retracted / update-to / updated-by
Impact67
Impact sub-metrics
Sub-metricRaw value0 to 100WeightPointsSource
Field-normalized impactImpact relative to the world average for the same field, year, and type, so a small field and a large field are judged fairly. Triangulates three independent field-normalized measures: OpenAlex's Field-Weighted Citation Impact, its citation percentile within the exact field-and-year cohort, and the NIH iCite Relative Citation Ratio. 1.0x is average; the percentile is the share of same-field, same-year work it out-cites.Maps to 0 to 100: Each ratio (1.0 = field average) maps 1.0 to 50, 3.0 to 75, 9.0 to 90; the field+year percentile maps directly (top 1% -> ~99). The available measures are averaged. When none is available yet (very recent work), this sub-metric is left out and Impact is computed from the others.39.6x FWCI, 1.0x RCR, top 0.1% of field-year vs field830.4537.3OpenAlex, NIH iCite fwci / citation_normalized_percentile / relative_citation_ratio
Influential citationsCitations that genuinely build on the work rather than mention it in passing (Semantic Scholar's influential-citation measure). A sharper signal of real impact than a raw count.Maps to 0 to 100: Saturating: 15 influential citations maps to 50, and none maps to 0. When Semantic Scholar has no measure for the work, this sub-metric is left out and Impact is computed from the others.Not yet availablen/aredistributed0.0Semantic Scholar influentialCitationCount
Citation velocityCitations received in the trailing twelve months (estimated from yearly citation counts: this year's so far plus the share of last year's that falls inside the window), judged as a rate rather than a lifetime total. Momentum is what a frontier index cares about: a result being taken up fast right now is landing, whether its lifetime pile is large or still small. Where an expected field citation rate is known, the momentum is also read relative to it, so a fast-moving result in a quiet field is not overlooked.Maps to 0 to 100: Saturating: 30 citations in the last twelve months maps to 50. When NIH iCite provides the field's expected rate, this is blended evenly with the ratio of observed to expected momentum.18 in last 12 months (1.2x field rate)460.3516.0OpenAlex, NIH iCite counts_by_year / field_citation_rate
Novelty78
Novelty sub-metrics
Sub-metricRaw value0 to 100WeightPointsSource
Conceptual noveltyHow genuinely new the contribution is, measured at publication from the work's OWN content, not its date: how atypical its combination of research areas is versus all prior work, and whether it builds on a broad, deep base rather than only the newest papers in a fast-moving area. This is what separates a real advance from a recent-but-commodity result, and it needs no forward citations, so it is stable for brand-new work. An honestly noisy proxy (even state-of-the-art novelty measures agree only moderately with expert judgment), so it carries bounded weight and never drives the ranking alone.Maps to 0 to 100: Two publication-time signals, each 0-100. The primary one is atypical combination: how unusual the work's median pairing of research topics is versus the prior literature, by pointwise mutual information (a pairing made far more often than chance scores low). The second, foundational reach (the median age of the work's references and the share older than five years), can only lift it: when reach is higher, the score rises by 35% of the gap, and it is never lowered. When only one signal is measurable it is used alone; when neither is (a work too thin in topics or references), this sub-metric is left out and Novelty is computed from the other three. A review is damped to 60%.Median topic pairing Genetic Associations and Epidemiology x Single-cell and spatial transcriptomics (896 prior works made this pairing); builds on refs a median 6.9y deep640.4528.7OpenAlex topics + referenced_works (publication-time)
RecencyHow recently the work appeared. The frontier is now, so newer scores higher.Maps to 0 to 100: Exponential decay with a 12-month half-life from the publication date.This month990.3029.8OpenAlex publication_date
Attention accelerationWhether attention is rising: citations in the latest full year versus the year before. Accelerating interest signals an active, opening frontier.Maps to 0 to 100: Latest/prior-year ratio: 1x maps to 50, 2x to 75, 4x to 100. The prior year counts as at least 3 citations, so a jump from one citation to a few is not read as explosive growth. With no citations in the latest full year it reads as unknown, and unknown is neutral (50).2.0x rising750.107.5OpenAlex counts_by_year
Edge of reviewPreprints sit ahead of the peer-review process and brand-new papers have only just cleared it. Both earn novelty here (while a preprint is discounted under Evidence), because that is where the frontier forms first.Maps to 0 to 100: Preprint 95; any other record 80 when under 6 months old, 65 under 18 months, 45 after that.Peer-reviewed record800.1512.0OpenAlex type / publication_date
ProvenanceAll sources

Sources and data

The source figures behind this entry, grouped by the source that holds them: up to five citation indices cross-checked, field-normalized impact, then real-world reach (Wikipedia, community discussion, clinical trials, released data and public funding) where it exists. Each value links to its record.

OpenAlex

Primary record
Total citations
24
Field-weighted citation impact
39.6x field average
Citation percentile
Top 0.1% of its field-year
Citations, recent 12 months
18
Publication date
2026-11-06

Crossref

Citations
0
Works it builds on
62

OpenCitations

Citations
0

NIH iCite

Relative Citation Ratio
1.04 (NIH-funded median 1.0)
Potential to translate
25%

Unpaywall

Reuse license
CC-BY

NIH RePORTER

NIH grants funding it
3
NIH award funding
$28,713,502

Figures quoted in the write-up

FigureValueSource
Immune cells analysed (Abstract: the OneK1K single-cell RNA sequencing dataset)1.2 million immune cellsNature Genetics
Additional eGenes found compared with pseudobulk analysis (Abstract and main text: against TensorQTL on mean pseudobulk counts with the same covariates and window size, gene-level P values by beta permutation (SAIGE-QTL used ACAT-V); eGenes summed over 14 cell types. The Discussion gives 48.4%)48.8% more eGenes than pseudobulk analysisNature Genetics
People and cell types in the OneK1K dataset (Main text)982 unique individuals, 14 distinct (sub)cell typesNature Genetics
eGenes found by SAIGE-QTL and by pseudobulk TensorQTL (Main text)17,218 versus 11,569 eGenes across all 14 cell typesNature Genetics
Rare-variant eGenes from set-based tests (Main text: Burden, SKAT and ACAT-V tests of rare variants within 1 Mb of each gene)5,541 rare variant (MAF ≤ 5%) eGenes across all 14 cell types (2,317 unique genes)Nature Genetics
Rare-variant eGenes independent of common eQTL signals (Main text: still significant (FDR < 0.05) after conditioning on all independent significant common eQTLs)483 (21%)Nature Genetics
Trans-eQTLs found genome-wide (Main text: at P < 5 × 10−8, MAF ≥ 10%)413 trans-eQTLs for 390 genes across the three cell typesNature Genetics
Speed-up over glmer, a single-cell Poisson mixed model (Main text: benchmark of 100 genes in B cells; Pearson R2 = 0.991 for −log10 P values)244-foldNature Genetics
How much faster the pseudobulk tools ran (Main text: SAIGE-QTL used less memory than both)TensorQTL and jaxQTL were 4.4-fold and 4.0-fold faster than SAIGE-QTLNature Genetics
eGenes found by pseudobulk TensorQTL given the same gene-level ACAT-V test (Main text: the paper's other TensorQTL comparisons use beta permutation)12,777 versus 11,569 eGenes (10.4% more than with beta permutation)Nature Genetics
For the authors

Show this score where your work lives

A badge for your lab page, your CV or the paper's code repository. It links readers to this entry, and it always carries the entry's current score.
This paper's Frontier Score badge, as it appears where you place itLive preview
<a href="https://frontierresearch.co/breakthroughs/efficient-and-accurate-mixed-model-association-tool-for-w4397016827?utm_source=badge&amp;utm_medium=badge&amp;utm_campaign=author-kit"><img src="https://frontierresearch.co/api/badge/W4397016827.svg" alt="Frontier Score for this paper" height="20"></a>

For a lab or personal website, a CV page, or anywhere that takes HTML.

The number is live. The badge is drawn from the same signals as this page, so it never shows a score the entry does not. Every entry is re-checked weekly against the same public sources, and a result's standing eases as it ages, so the active edge of the field stays on top. How scoring works

Your other papers. The free paper check shows any paper's field-weighted citation impact, its citations in five open indexes and any retraction on record, whether or not it is in the index, and the researcher check does it for your whole record from your ORCID iD.

Something wrong or out of date? Write to hello@frontierresearch.co and the entry will be corrected.

Biomedicine & HealthThe whole field

Back to the Frontier Index, or browse every breakthrough by date.

The Frontier Brief

The week's Biomedicine & Health, scored.

Every Monday, up to six of the week's new Biomedicine & Health entries in the index, each with its score and a one-line read on why it matters, then the week's best from the rest of the frontier. What you get

Prefer a feed reader? Follow Biomedicine & Health by RSS

One email a week: up to six of the week's new Biomedicine & Health entries, scored, then the best elsewhere. Unsubscribe anytime.