A team of Russian researchers has produced one of the most detailed maps to date of how human genetic variation controls the behavior of immune cells, linking DNA differences to gene activity one cell type at a time and connecting those regulatory effects to COVID-19. The study, published in the journal Immunogenetics, combined whole-genome sequencing with single-cell RNA sequencing of more than 230,000 peripheral blood mononuclear cells from 30 individuals, identifying over 1.2 million expression quantitative trait loci, or cis-eQTLs, across 18 distinct immune cell types. The work represents a significant step forward in the effort to understand why people respond so differently to infectious diseases, and it demonstrates how deep learning can be used to interpret the biological meaning of variants that would otherwise remain statistical abstractions.
Expression quantitative trait loci are genomic positions, typically single nucleotide variants, that correlate with differences in the activity of nearby genes. In classical eQTL studies, which aggregate gene expression measurements across whole tissues, these associations reveal only the average effect of a variant across a mixture of cell types. That averaging is a serious limitation for immunology, because peripheral blood contains dozens of specialized cell populations—T cells, natural killer cells, monocytes, B cells, and many rarer intermediates—each with its own regulatory grammar. A variant that strongly boosts a gene in natural killer cells may be invisible in a bulk measurement dominated by monocytes. By pairing each donor’s genome with single-cell transcriptomes, the researchers could test the association between genotype and expression separately within each cell population, resolving effects that bulk approaches blur away.
The technical pipeline behind the study was substantial. Whole-genome sequencing data were processed through established variant-calling workflows, with germline small variants identified using the Strelka2 caller, following best practices endorsed by the Global Alliance for Genomics and Health for benchmarking variant calls. Coverage calculations were performed with Mosdepth, and sequencing quality was assessed against community standards to ensure that the genetic data underpinning the eQTL analysis were reliable. On the transcriptomic side, single-cell RNA-seq data from the 230,000 PBMCs were integrated across donors using the Harmony method, a widely adopted algorithm that corrects for batch effects and donor-specific technical noise while preserving genuine biological differences between cell types. Cells were then classified into 18 populations spanning the major branches of the immune system, providing the cellular resolution required for cell-type-specific association testing.
With genotypes and cell-type-resolved expression in hand, the team carried out cis-eQTL mapping using Matrix eQTL, a computational framework designed to perform the enormous matrix operations required for fast association testing across millions of variants and thousands of genes. The result was a catalog of 1,233,644 cis-eQTLs distributed across the 18 cell types. Importantly, this was not simply a numbers exercise. The researchers subjected their findings to a series of validation and interpretive analyses designed to ask whether the variants they detected showed the hallmarks of genuine regulatory elements, and, more provocatively, what evolutionary forces have shaped them.
One of the most intriguing findings to emerge from these secondary analyses concerns evolutionary conservation. When the team examined the genomic regions harboring their strongest eQTLs, they found that the most statistically significant associations tended to sit in less conserved regions of the genome—stretches of DNA that have diverged relatively rapidly between species. These variants were also concentrated in the regulatory regions of more divergent genes. This pattern suggests that immune gene regulation is an evolutionary hotspot, consistent with the well-documented observation that genes involved in host defense are frequent targets of positive selection. Rapid turnover of regulatory elements may allow populations to adapt to changing pathogen landscapes, but it may also help explain why immune-related variants are a rich source of susceptibility to chronic inflammatory and autoimmune disease in modern humans. The finding carries a double edge: the same regulatory flexibility that enabled adaptation to ancient pathogens may predispose contemporary genomes to misfire.
To move from statistical association to mechanistic understanding, the researchers turned to deep learning models of cis-regulatory sequence. Neural networks trained on genomic data can learn the relationship between DNA sequence and regulatory function, predicting how transcription factors bind to specific sequence contexts and how single-base changes alter those interactions. Drawing on approaches pioneered by tools such as DeepSEA, Basset, and the Enformer-style sequence models developed in recent years, the team applied these computational models to their eQTL catalog to ask, for each variant, which transcription factor binding sites are disrupted and in which cellular context that disruption matters. This step transformed the analysis from a list of correlated positions into a functional hypothesis-generating resource: each eQTL could now be annotated with a predicted mechanism of action grounded in sequence-level regulatory biology.
The functional analysis focused on genes with well-established roles in immunity, including NKG7, members of the HLA family, MIF, and MS4A1. NKG7 encodes a protein essential for the cytotoxic function of natural killer cells and CD8 T cells, involved in the trafficking of lytic granules that deliver the killing blow to infected or malignant target cells; variants affecting its expression could plausibly modulate antiviral and antitumor immunity. HLA genes, which encode the human leukocyte antigen molecules that present viral peptides to T cells, are among the most polymorphic loci in the human genome and have been repeatedly implicated in COVID-19 susceptibility and severity. MIF, the macrophage migration inhibitory factor, is a potent inflammatory mediator whose circulating levels correlate with severe COVID-19 pneumonia. MS4A1, better known as CD20, defines B cells and is the target of widely used monoclonal antibody therapies. By tracing how variants influence the expression of these genes through disrupted transcription factor binding—including factors such as those involved in myeloid and lymphoid differentiation—the study connected genetic variation to plausible cellular mechanisms.
The disease context of the study is explicit in its framing. Prior work by overlapping teams had used single-cell transcriptomics to identify immune cell signatures associated with severe Delta-variant COVID-19, and large consortia had shown that severe disease is marked by a dysregulated myeloid cell compartment. The GWAS Catalog lists numerous loci associated with COVID-19 outcomes, but translating those associations into functional biology has remained a central challenge for the field. The new eQTL resource provides a bridge: because it assigns regulatory variants to specific immune cell types and annotates their predicted effects on gene expression, it can be used to prioritize candidate causal variants among the many statistical associations emerging from disease-genetic studies. The authors also point toward broader applications, noting parallels with single-cell eQTL studies in brain, autoimmune disease, and other contexts that have revealed how cell-type-specific regulation shapes genetic risk.
The scale of the cellular census deserves emphasis. Thirty donors may sound modest compared with the million-person cohorts of cardiovascular genetics, but single-cell eQTL studies trade sample size for resolution: every donor contributes tens of thousands of individual cells, each measured across the whole transcriptome. The 230,000-cell dataset allowed the researchers to detect eQTLs not only in abundant populations like CD4 T cells and monocytes but in rarer states such as intermediate monocytes and specific lymphocyte subsets, where regulatory effects would be hopelessly diluted in bulk analysis. The identification of 18 cell-type-resolved regulatory landscapes from this relatively small cohort illustrates the power of the approach, and the authors’ analytical strategy—conservation analysis, transcription factor modeling, and pathway interrogation—provides a template that larger consortia are likely to follow as single-cell eQTL mapping matures.
Methodologically, the study also underscores how much of modern genomics depends on careful assembly of open computational tools. Beyond the core machinery of variant calling and eQTL mapping, the researchers drew on annotation frameworks for regulatory elements derived from more than a thousand epigenomic datasets, packages for gene ontology and pathway enrichment, tools for predicting DNA shape features that influence transcription factor binding, and simulation frameworks for testing the significance of overlaps between genomic intervals. This layered infrastructure allowed a single research group to integrate population genetics, transcriptomics, chromatin biology, and machine learning into a coherent narrative about how sequence variation becomes functional variation in the immune system.
The limitations of the study are those inherent to its design. A cohort of 30 individuals restricts statistical power for detecting rarer variants and weaker regulatory effects, and the donor population limits generalizability across ancestries—an important consideration given that eQTL effects and linkage patterns differ among populations. The deep learning predictions, while mechanistically informative, remain computational hypotheses that would need experimental validation through reporter assays or CRISPR-based perturbation of individual variants. Nevertheless, the study’s value lies in the framework it establishes: a complete chain from genome sequence, through cell-type-resolved expression, to predicted transcription factor mechanism, anchored to a disease of global significance.
As single-cell sequencing costs continue to fall and paired genotyping-transcriptomics cohorts grow, resources of this kind are expected to expand rapidly in scale and diversity. What this study demonstrates is that the payoff of such investment is not merely a longer list of associations, but a progressively sharper picture of the regulatory code that governs human immunity—and with it, new opportunities to understand, predict, and ultimately intervene in diseases where the immune system holds the balance between recovery and catastrophe. For COVID-19, whose genetic architecture continues to be dissected years after the pandemic’s peak, that sharper picture may help explain at last why the same virus produces a mild illness in one person and a life-threatening one in another.
Subject of Research: Cell-type-specific cis-eQTL mapping of immune cell function in COVID-19 using paired whole-genome sequencing and single-cell RNA sequencing of human peripheral blood mononuclear cells
Subject of Research: Biology
Article Title: Deciphering the genetic control of immune cell function at single-cell resolution: Disease-Specific Cis-eQTLs analysis of COVID-19
Article References: Romanova, E. I., Tychinin, D. I., Shaymardanov, A. M., Akimov, V. E., Korobeinikova, A. V., Shiryagin, V. V., Guskova, N. I., Astafieva, V. A., Shingaliev, A. S., Antonova, O. A., Golubnikova, L. A., Mitrofanov, S. I., Grammatikati, K. S., Yudin, V. S., Yudin, S. M., Makhotenko, A. V., Keskinov, A. A., Kraevoy, S. A., Snigir, E. A., … Skvortsova, V. I. (2026). Deciphering the genetic control of immune cell function at single-cell resolution: Disease-Specific Cis-eQTLs analysis of COVID-19. Immunogenetics, 78(1), Article 4. https://doi.org/10.1007/s00251-026-01396-0
Image Credits: AI Generated
DOI: 10.1007/s00251-026-01396-0
Keywords: cis-eQTL, single-cell RNA sequencing, COVID-19, immune cells, PBMCs, whole-genome sequencing, deep learning, gene regulation, transcription factors, HLA, NKG7, SNP
Cite Scienmag News
APA MLA Chicago
Juliet Wilcox. (September 9, 2026). Single-cell eQTL analysis reveals genetic control of immune cells in COVID-19. Scienmag. https://scienmag.com/single-cell-eqtl-analysis-reveals-genetic-control-of-immune-cells-in-covid-19/
Juliet Wilcox. “Single-cell eQTL analysis reveals genetic control of immune cells in COVID-19.” Scienmag, 9 September 2026, https://scienmag.com/single-cell-eqtl-analysis-reveals-genetic-control-of-immune-cells-in-covid-19/. Accessed 9 September 2026.
Juliet Wilcox. “Single-cell eQTL analysis reveals genetic control of immune cells in COVID-19.” Scienmag. September 9, 2026. https://scienmag.com/single-cell-eqtl-analysis-reveals-genetic-control-of-immune-cells-in-covid-19/
Copy citation Download RIS
Tags: cell-type-specific gene regulationCOVID-19 host genetic factorsCOVID-19 immune responseCOVID-19 immune response geneticsdeep learning for genetic variant interpretationdeep learning in genetic variant interpretationexpression quantitative trait loci in immune cellsexpression quantitative trait loci in immunitygene activity in immune cellsgenetic basis of differential COVID-19 responsesgenetic control of immune cell behaviorgenetic regulation of immune cellshuman genetic variation and infectious diseasesimmune cell genetic variationimmune cell type mapping in genetic studiesimmune cell type-specific gene regulationimpact of DNA variants on immune cell behaviorperipheral blood mononuclear cellsperipheral blood mononuclear cells genetic mappingsingle-cell eQTL analysissingle-cell eQTL analysis in COVID-19single-cell RNA sequencing in immunogeneticssingle-cell RNA sequencing of immune cells
