Tandem repeat variation shapes immune cell type-specific gene expression

By integrating multiomic data from over 5.4 million immune cells across nearly 2,000 individuals, this study reveals that tandem repeat variation is a major, cell type-specific regulator of gene expression and a key driver of complex immune and hematological traits.

The Gist

Imagine your genome (your body's instruction manual) is a massive library. For decades, scientists have been obsessed with reading the "main text"—the specific letters (A, C, T, G) that make up our genes. They've been very good at spotting typos where a single letter is wrong.

But there's a whole other section of the library that has been largely ignored: the Tandem Repeats.

Think of Tandem Repeats like a chorus in a song or a stutter in a sentence. Instead of a single letter, it's a short phrase that gets repeated over and over.

• Normal: "The cat sat on the mat."
• Tandem Repeat: "The cat sat sat sat sat on the mat."

In the past, trying to count exactly how many times "sat" was repeated was like trying to count the rungs on a ladder while it's shaking. It was too blurry and difficult to get right. Because of this, scientists mostly ignored these repeats, assuming they were just "junk" noise.

This paper is a game-changer because it finally learned how to count the rungs on that shaky ladder.

Here is what the researchers discovered, broken down into simple concepts:

1. The "Cell Type" Detective Work

Previously, scientists looked at blood samples as a big, blended smoothie. They couldn't tell if a genetic instruction was coming from a T-cell, a B-cell, or a macrophage. It was like trying to figure out who in a crowded room is speaking by listening to the whole room at once.

This study used single-cell technology, which is like giving every person in that crowded room a microphone. They looked at over 5.4 million individual immune cells from nearly 2,000 people. They found that these "repeating phrases" (Tandem Repeats) act as volume knobs for genes, but only for specific types of cells.

• Analogy: A repeat might turn up the volume on a "fight infection" gene in a T-cell, but leave it silent in a B-cell.

2. The "Volume Knob" Effect

The researchers found over 69,000 instances where the length of a repeat directly changed how much a gene was expressed.

• The Metaphor: Imagine a gene is a lightbulb. A Tandem Repeat is the dimmer switch.
  • If you have a short repeat (e.g., 10 repeats), the light is dim.
  • If you have a long repeat (e.g., 20 repeats), the light is bright.
  • Crucially, they found that for many genes, the "dimmer switch" is made of these repeats, not the single-letter typos (SNPs) scientists usually study.

3. The "Cellular State" Surprise

Cells aren't static; they change their mood and job. A T-cell might be "resting" or "angry" (activated). The study found that some of these repeat switches only work when the cell is in a specific "mood."

• Analogy: It's like a car's cruise control. The switch (the repeat) only turns on the speed control when you are on the highway (activated cell), but does nothing when you are in a parking lot (resting cell).

4. The "Architectural" Connection

To understand how these repeats work, the team also looked at the "scaffolding" of the cell (chromatin accessibility).

• The Discovery: They found that when a repeat gets longer, it often physically pulls the DNA open, making it easier for the cell's machinery to read the gene.
• The Metaphor: Imagine the DNA is a tightly wound ball of yarn. These repeats act like a hook that, when they get longer, pull a loop of the yarn loose, making it easy to grab and read the instructions.

5. Why This Matters for Disease

The most exciting part is that these repeats are linked to real-world diseases.

• The team found that these "volume knobs" are likely the reason why some people are more prone to asthma, inflammatory bowel disease, or blood disorders.
• The "Missing Link": For years, genome studies found a "suspected criminal" (a location in the DNA) associated with a disease, but they couldn't find the "weapon" (the specific mutation). This paper suggests that for many diseases, the weapon is actually a Tandem Repeat that was previously too blurry to see.

The Big Picture

This paper is like upgrading from a black-and-white, blurry map to a high-definition, 3D GPS.

• Old View: We thought our genes were controlled mostly by single-letter typos.
• New View: We now know that "repeating phrases" (Tandem Repeats) are massive, dynamic regulators that act as specialized volume knobs for different immune cells.

By learning to read these repeats, we are finally unlocking a hidden layer of human biology that explains why our immune systems work the way they do, and why some people get sick while others stay healthy. It's a reminder that in the library of life, sometimes the most important instructions aren't the unique words, but the rhythm of the repetition.

Technical Summary

1. Problem Statement

Tandem repeats (TRs), including short tandem repeats (STRs) and variable number tandem repeats (VNTRs), constitute a highly polymorphic and mutable class of human genetic variation. Despite their known potential to regulate gene expression by altering transcription factor binding, nucleosome occupancy, and promoter/enhancer activity, they remain significantly underexplored in population genetics.

• Limitations of Current Approaches: Historically, short-read sequencing has struggled to accurately genotype repetitive regions, leading to reliance on imputation from surrounding SNVs or limited catalog definitions.
• Resolution Gap: Most existing expression quantitative trait locus (eQTL) studies rely on bulk RNA-seq, which averages signals across heterogeneous cell populations, thereby masking cell type-specific regulatory effects.
• Knowledge Gap: There is a lack of comprehensive, genome-wide maps linking TR variation to gene expression at single-cell resolution, particularly regarding how TRs influence chromatin accessibility and complex trait architecture in the immune system.

2. Methodology

The study leverages TenK10K Phase 1, the largest paired human whole-genome sequencing (WGS) and single-cell RNA sequencing (scRNA-seq) resource to date.

• Cohorts: 1,925 individuals (European ancestry) from the Tasmanian Ophthalmic Biobank (TOB) and the BioHEART cohort, providing matched WGS and scRNA-seq data.
• Multi-omic Integration:
  • scRNA-seq: >5.4 million peripheral blood mononuclear cells (PBMCs) across 28 immune cell types.
  • scATAC-seq: Matched single-cell ATAC sequencing from >3.4 million nuclei in 922 individuals.
  • Long-read Validation: PacBio HiFi sequencing in 25 individuals and capillary electrophoresis for truth set generation.
• TR Catalog Curation:
  • Merged multiple reference-based and population-derived catalogs (including 1000 Genomes and Human Pangenome Reference Consortium) to create a unified catalog of ~4.9 million loci.
  • Genotyped ~2.6 million polymorphic TR loci using ExpansionHunter (v5) after rigorous quality filtering, resulting in a high-confidence callset of ~1.7 million autosomal TRs.
• Association Analysis:
  • sc-eTR Mapping: Used associaTR to perform pseudobulk eQTL mapping within ±100 kb of genes, testing associations between TR length and gene expression across 28 cell types. Meta-analysis combined results from TOB and BioHEART cohorts.
  • Cell State Inference: Applied scDeepID (a deep learning framework) to infer continuous cell states (e.g., NK cell activation) and test for state-dependent TR effects.
  • caQTL Mapping: Tested associations between TR length and chromatin accessibility (ATAC peaks) in matched cell types.
  • Fine-mapping: Used SuSiE to identify candidate causal TRs (Posterior Inclusion Probability, PIP ≥ 0.7) and distinguish them from linked SNVs.
  • Colocalization: Performed Bayesian colocalization (using coloc) between sc-eTRs and GWAS signals for immune/hematological traits and diseases.

3. Key Contributions

• First Large-Scale sc-eTR Atlas: Generated a genome-wide catalog of >69,000 single-cell expression TR loci (sc-eTRs) across 28 immune cell types.
• Multi-omic Mechanistic Insight: Integrated scRNA-seq, scATAC-seq, and DNA methylation data to demonstrate that TRs drive coordinated changes in chromatin accessibility and gene expression.
• Cell Type Specificity Framework: Developed a classification system to distinguish between shared and cell type-specific TR regulatory effects, revealing that ~30.7% of sc-eTRs are specific to a single cell type.
• Causal Inference: Identified 1,490 fine-mapped candidate causal TRs, demonstrating that direct TR genotyping captures regulatory variation often missed by SNV-based approaches.

4. Key Results

A. Discovery and Specificity of sc-eTRs

• Scale: Identified 69,210 unique sc-eTRs associated with 15,889 genes (FDR < 5%).
• Cell Type Specificity: 30.7% of sc-eTRs are specific to one of the 28 immune cell types. Cell type-specific TRs are enriched in intergenic regions and di/tetranucleotide motifs, whereas shared TRs are enriched in proximal promoters and homopolymers.
• Dynamic Regulation: Cell state inference revealed 579 "cell state eTRs" where the regulatory effect of a TR depends on the functional state of the cell (e.g., NK cell activation levels modulating SYNGR1 expression).
• Comparison to Bulk: 91.5% of sc-eTRs overlapped with bulk blood eTRs, but single-cell resolution resolved many bulk signals to specific cell types (e.g., FCLR5 in naive B cells) and uncovered unique cell type-specific drivers not visible in bulk data.

B. Mechanistic Link to Chromatin and Methylation

• caQTLs: Nearly one-third (31.8%) of expression-associated TRs were also associated with local chromatin accessibility (caQTLs).
• Directionality: 78% of caQTL+ sc-eTRs showed concordant effects: longer repeat lengths generally increased chromatin accessibility and gene expression, particularly in promoters and 5' UTRs.
• Multi-omic Examples:
  • POMC: A poly(GCCGCTGCT) repeat in CD8TEM cells associated with increased expression, increased accessibility, and decreased methylation.
  • INF2: An intronic repeat in CD4TCM cells associated with decreased expression and accessibility.

C. Causal Drivers and Fine-mapping

• Causal TRs: Fine-mapping identified 1,490 TRs as candidate causal drivers (PIP ≥ 0.7) for 1,354 genes.
• Independence from SNVs: Many candidate causal TRs showed moderate linkage disequilibrium (LD) with lead SNVs (median $r^2$ = 0.57). Over 50% remained significant after conditioning on the lead SNV, indicating they capture independent regulatory variation.
• Genomic Context: Causal TRs were strongly enriched in proximal regulatory elements (promoters, 5' UTRs) and GC-rich regions, often overlapping constrained primate-specific transcription factor binding sites.

D. Impact on Complex Traits

• Colocalization: Identified 215 genes (502 gene-trait pairs) where candidate causal sc-eTRs colocalize with GWAS signals for blood traits, serum markers, and immune diseases (e.g., asthma, IBD, Parkinson's).
• Contextual Drivers:
  • LIME1: A poly(TG) repeat in CD4 naive cells colocalizes with asthma risk.
  • KIF16B: A poly(T) repeat in NK cells colocalizes with mean platelet volume.
  • UBA7: A poly(CG) repeat in NK cells colocalizes with inflammatory bowel disease (IBD).
• Disease Relevance: 24.9% of genes linked to causal sc-eTRs overlap with clinically relevant gene panels (PanelApp), suggesting roles in variable penetrance or modifier effects for Mendelian disorders.

5. Significance

• Redefining Regulatory Architecture: This study establishes TRs as a distinct, major layer of regulatory variation that shapes cell type-specific gene expression, often in ways not captured by SNVs.
• Methodological Advancement: It demonstrates the critical necessity of direct TR genotyping and multi-omic single-cell approaches to fully resolve the genetic architecture of complex traits.
• Clinical Implications: By linking TR variation to specific immune cell contexts and complex diseases, the findings provide new candidates for therapeutic targeting and improve the interpretation of non-coding genetic variation in disease risk.
• Resource Generation: The study provides a publicly available catalog of polymorphic TRs, summary statistics, and fine-mapped causal variants, serving as a foundational resource for future population genomics and functional studies.

In conclusion, this work moves beyond the "SNV-centric" view of genetics, revealing that tandem repeat variation is a pervasive, dynamic, and cell-type-specific driver of human immune regulation and disease susceptibility.