The immune system's response to tissue injury, infection, or self-antigens depends on a system of activating and restraining signals. Among the activating signals, the OX40–OX40L axis¹¹ OX40–OX40L axis
OX40 (CD134) is expressed on activated T cells; OX40L (CD252), encoded by TNFSF4, is expressed on antigen-presenting cells and delivers a survival and proliferation signal that amplifies the T cell response is one of the most potent co-stimulatory pathways in adaptive immunity. rs1234314 is an intronic variant in TNFSF4 that reduces the gene's promoter activity, altering the quantity of OX40L produced. It provides an independent association signal at the TNFSF4 locus, separate from the well-characterized upstream risk haplotype tagged by rs7514229 and rs2205960, and has been associated with susceptibility to systemic sclerosis (scleroderma) and systemic lupus erythematosus across multiple cohorts.

The Mechanism

rs1234314 sits within an intron of TNFSF4 on chromosome 1 (GRCh38 chr1:173208253), a gene transcribed from the minus strand. The plus-strand notation is C>G, with C as the GRCh38 reference allele. A 2023 promoter-reporter study²² 2023 promoter-reporter study
Chen et al. J Clin Med 2023 (PMID 36983159) — dual-luciferase assay; >10 independent experiments per variant; one-way ANOVA p=0.003 demonstrated that the G allele drives approximately 68% less transcriptional activity than the C allele in promoter-reporter constructs — the G-allele construct produced only 0.32 ± 0.09 times the relative light units of the C-allele reference. This functional effect establishes rs1234314 as a regulatory variant that modulates TNFSF4 expression, not merely a passive tag for a nearby causal variant.

The consequence of reduced OX40L expression is not straightforwardly protective. OX40L is required not only for effector T cell activation but also for the differentiation and persistence of regulatory T cells (Tregs³³ Tregs
Regulatory T cells that suppress immune activation and maintain tolerance to self-antigens; their function depends partly on OX40 signaling for maintenance in inflamed tissues). Insufficient OX40 co-stimulation may selectively impair Treg maintenance in inflamed tissue, while early effector T cell responses proceed through other co-stimulatory pathways — a net dysregulation that promotes autoimmune chronicity.

The Evidence

The first identification of rs1234314 as a risk locus came from a 2010 case-control study of systemic sclerosis⁴⁴ 2010 case-control study of systemic sclerosis
Gourh et al. Ann Rheum Dis 2010 (PMID 19778912) — 1,059 SSc patients, 698 controls; nine TNFSF4 SNPs tested; FDR-corrected p=0.019 in North American patients. The minor (G) allele at rs1234314 conferred an odds ratio of 1.20 (95% CI 1.04–1.40) for SSc susceptibility. This was replicated⁵⁵ replicated
Bossini-Castillo et al. Ann Rheum Dis 2011 (PMID 21187296) — 3,014 SSc cases and 3,125 controls across Germany, Spain, The Netherlands, Italy, France, UK, Czech Republic, and Norway; multi-cohort meta-analysis in the largest European SSc genetics cohort assembled at that time (n=6,139), yielding an overall OR of 1.15 (95% CI 1.02–1.31), with stronger effects in limited cutaneous SSc (lcSSc; OR 1.22) and anti-centromere antibody (ACA)-positive patients (OR 1.23).

For systemic lupus erythematosus, a trans-ancestral fine-mapping study⁶⁶ trans-ancestral fine-mapping study
Manku et al. Ann Rheum Dis 2013 (PMID 23874208) — 17,900 SLE and control subjects spanning Amerindian/Hispanic, African-American, European, and East Asian ancestries; conditional regression analyses identified rs1234314-C as an independent non-risk signal at the TNFSF4 locus — meaning the C allele is associated with protection, and the G allele carries SLE risk as a signal independent of the primary rs2205960-T risk haplotype. A Taiwanese case-control study⁷⁷ Taiwanese case-control study
Chen et al. Front Immunol 2022 (PMID 36189300) corroborated this direction: rs1234314 genotype distribution differed significantly between SLE cases and controls (CC vs. CG vs. GG: p=0.005), with GG vs. CC comparison reaching p=0.004.

Practical Implications

The TNFSF4 locus, via multiple independent signals including rs1234314, is one of the most consistently replicated genetic associations for both scleroderma and lupus. Carriers of the G allele have a modestly elevated background risk for these conditions — approximately 15–20% higher odds per the SSc replication data. The clinical translation is not treatment of genetic risk, but informed monitoring: recognizing early features of these conditions that warrant specialist evaluation.

Systemic sclerosis (scleroderma) often begins with Raynaud's phenomenon⁸⁸ Raynaud's phenomenon
Episodic vasospasm causing fingers and toes to turn white or blue in response to cold or stress; present in >95% of SSc patients, often preceding skin and organ involvement by years — abnormal cold sensitivity — and progresses to skin thickening, internal organ fibrosis, and vascular disease. Early detection matters because immunosuppressive treatment is most effective before irreversible fibrosis occurs.

Interactions

rs1234314 provides an independent genetic signal at the TNFSF4 locus that is separate from the upstream risk haplotype tagged by rs7514229 and rs2205960. The trans-ancestral fine-mapping study (PMID 23874208) demonstrated that conditional regression analysis delineates the primary risk signal to rs2205960-T and the independent signal to rs1234314-G, confirming these two signals coexist without full mutual explanation. A carrier of both the rs7514229-T (upstream haplotype) and rs1234314-G alleles would carry two independent TNFSF4 risk signals, a combination worth noting in the context of autoimmune risk assessment.

Endometriosis — in which tissue resembling the uterine lining grows and implants outside the uterus — affects an estimated 10% of women of reproductive age and accounts for a major share of chronic pelvic pain, dyspareunia, and infertility. Its causes are multifactorial, but genetic factors account for roughly half of the susceptibility variance. FN1, encoding fibronectin 1¹¹ fibronectin 1
a large glycoprotein of the extracellular matrix (ECM) that provides scaffolding for cell adhesion, migration, proliferation, and tissue remodeling, is one of the first confirmed common genetic loci for this disease.

The rs1250248 variant sits in an intronic region of FN1 between exons 10 and 11. It does not change the fibronectin protein sequence, but the position overlaps a predicted transcription factor-binding site, raising the possibility that the A allele alters the regulation of FN1 expression in endometrial tissue.

The Mechanism

Fibronectin is a critical mediator of extracellular matrix remodeling²² extracellular matrix remodeling
ECM remodeling refers to the continuous process of synthesis and degradation of matrix proteins that governs cell behavior; dysregulated ECM remodeling is a hallmark of fibrosis and invasion in endometriotic lesions. In peritoneal endometriotic lesions, fibronectin is overexpressed compared with eutopic endometrium, and FN1-integrin signaling at the interface of mesothelial cells and ectopic endometrial stromal cells may promote progesterone resistance and lesion persistence.

A 2025 single-cell and spatial transcriptomic study identified a specific CXCR4⁺ fibroblast subpopulation³³ CXCR4⁺ fibroblast subpopulation
Fibroblasts are the main ECM-producing cells; the CXCR4⁺ subset identified here had high stemness and proliferative capacity, acting as a hub for FN1 signaling in immune and fibrotic responses within ectopic lesions as a central mediator of FN1-driven immune remodeling and fibrosis. This subpopulation coordinates both ECM deposition and immune suppression, creating an environment that may sustain ectopic implant viability. The intronic rs1250248 A allele may subtly shift FN1 expression in this context.

Additionally, plasma fibronectin concentrations are significantly elevated in women with endometriosis⁴⁴ significantly elevated in women with endometriosis
Fibronectin 292.6 ± 96.2 mg/L in endometriosis vs 226.6 ± 91.9 mg/L in controls; high-molecular-mass fibronectin-fibrin complexes absent in healthy women but present in endometriosis patients, and novel fibronectin-fibrin complexes have been identified exclusively in affected women, suggesting altered fibronectin molecular biology is a downstream consequence of the disease process that rs1250248 may predispose toward.

The Evidence

The association between rs1250248 and endometriosis was first identified in a genome-wide association study of 3,194 surgically confirmed cases and 7,060 controls from Australia and the UK⁵⁵ genome-wide association study of 3,194 surgically confirmed cases and 7,060 controls from Australia and the UK
Painter et al. Nature Genetics, 2011. The rs1250248 locus in FN1 reached P = 3.2 × 10⁻⁸ in the Stage III/IV-restricted analysis of the discovery dataset, though it did not replicate independently in a US validation cohort, leaving the signal sub-threshold at the time.

A subsequent meta-analysis of eight GWAS datasets in 11,506 cases and 32,678 controls⁶⁶ meta-analysis of eight GWAS datasets in 11,506 cases and 32,678 controls
Rahmioglu et al. Human Reproduction Update, 2014 confirmed the FN1 locus signal with greater precision. Across all endometriosis cases the A allele showed OR = 1.11 (P = 1.1 × 10⁻⁴), rising to OR = 1.26 (95% CI 1.16–1.38, P = 8.0 × 10⁻⁸) when restricted to Stage III/IV disease. The lead SNP in that analysis was rs1250241 (r² = 0.95 with rs1250248), representing the same genetic signal.

An Italian study in 305 laparoscopically confirmed cases and 2,710 controls⁷⁷ Italian study in 305 laparoscopically confirmed cases and 2,710 controls
Pagliardini et al. J Med Genet, 2013 confirmed genome-wide significance specifically for severe disease (P = 3.89 × 10⁻⁹) and identified an epistatic interaction with rs7521902 (WNT4): the combined effect of risk alleles at both loci reached OR = 1.56 overall and OR = 2.15 specifically for ovarian endometriosis, suggesting these two ECM-related pathways converge in promoting cyst formation.

A Greek case-control study⁸⁸ Greek case-control study
Matalliotaki et al. Mol Med Rep, 2019 found the A allele conferred OR = 1.87 (95% CI 1.34–2.61, P = 0.002); the AA genotype was present in 16.9% of patients versus 9.5% of controls (OR = 2.59). Notably, this cohort found the strongest effect in early-stage (I/II) disease — suggesting rs1250248 may have distinct stage-specific effects across populations.

Practical Implications

Carrying the A allele at rs1250248 raises the likelihood of endometriosis, with the clearest signal for moderate-to-severe disease in European populations. The absolute risk conferred by a single common intronic variant remains modest on an individual basis, but the biological plausibility is strong: fibronectin overexpression in ectopic tissue and elevated plasma fibronectin are both documented features of endometriosis, and FN1 signaling coordinates the fibrosis and immune evasion that sustain lesion persistence.

Women carrying the A allele — particularly AA homozygotes — may benefit from heightened vigilance around endometriosis symptoms, proactive specialist referral, and awareness of the potential for fibronectin-related ECM involvement when discussing treatment options with a clinician. N-acetylcysteine (NAC) has shown benefit in multiple endometriosis studies⁹⁹ benefit in multiple endometriosis studies
Porpora et al. 2013 (PMID 23737821): NAC-treated patients showed reduced cyst size after 3 months vs increase in controls; reduces oxidative stress and ECM-promoting pathways and may work in part through reduction of the oxidative ECM signaling environment that fibronectin overexpression supports.

Interactions

rs7521902 (WNT4 locus): An epistatic interaction between rs1250248 (FN1) and rs7521902 (WNT4) has been formally demonstrated in the Italian replication dataset. For women carrying the A risk allele at both loci, the combined OR for ovarian endometriosis reaches 2.15 (P = 3.12 × 10⁻⁴). WNT4 regulates Müllerian duct development and suppresses androgen production; its intersection with fibronectin-mediated ECM biology in the ovarian microenvironment represents a plausible convergent pathway for endometrioma formation.

Supervisor compound action proposal: women carrying the A risk allele at rs1250248 (FN1) AND the risk allele at rs7521902 (WNT4 locus) represent a subgroup with substantially elevated ovarian endometriosis risk (OR ~2.15). Combined recommendation: early transvaginal ultrasound surveillance targeting ovarian endometrioma formation, proactive AMH testing to establish baseline ovarian reserve, and expedited fertility consultation if conception is desired. Evidence level: moderate (single formal interaction study, biologically plausible).

rs12700667 (7p15.2 locus, near HOXA10/HOXA11): This locus, which influences HOXA10/11 regulation and endometrial receptivity, has been identified as another major endometriosis susceptibility locus. Both the 7p15.2 and FN1 loci are among the most replicated GWAS signals for endometriosis. Additive effects of risk alleles across multiple loci are expected under a polygenic architecture, though formal interaction testing between rs12700667 and rs1250248 has not been published.

Complement Factor B (CFB) is the central amplification enzyme of the alternative complement pathway¹¹ alternative complement pathway
One of three complement activation pathways; the alternative pathway provides continuous low-level surveillance and amplifies the classical and lectin pathways, the arm of the innate immune system responsible for continuous immune surveillance and the amplification of inflammation. The rs1270942 variant sits within an intron of the CFB gene in the major histocompatibility complex (MHC) class III region on chromosome 6 and was identified in the landmark 2008 SLEGEN genome-wide association study²² identified in the landmark 2008 SLEGEN genome-wide association study
Harley et al. Nature Genetics 2008 — 720 SLE cases, 2,337 controls, >317,000 SNPs as one of the strongest non-HLA signals for systemic lupus erythematosus (SLE), with an odds ratio of 2.35 and p=1.3×10⁻⁵¹. This association has since been replicated in larger multi-ancestry meta-analyses with p-values reaching 2×10⁻¹⁶⁵.

The Mechanism

CFB encodes the complement factor B protein, which binds to activated C3b on foreign or damaged surfaces to form the alternative pathway C3 convertase (C3bBb). This convertase amplifies complement deposition on target surfaces, tags them with opsonins, activates inflammatory signaling, and triggers the terminal membrane attack complex. In healthy individuals, the alternative pathway operates at low constitutive levels to clear cellular debris, apoptotic cells, and immune complexes — functions essential for preventing autoimmunity.

The rs1270942 G allele is in tight linkage disequilibrium³³ linkage disequilibrium
LD means nearby variants are inherited together; rs1270942 is a tag SNP for the broader CFB haplotype with functional variants across the CFB locus that alter complement activity. The MHC class III haplotype bearing the G allele is associated with dysregulated alternative pathway amplification⁴⁴ dysregulated alternative pathway amplification
Impaired regulation allows C3b to accumulate on self-surfaces and within glomerular immune complexes, reducing the precision of immune complex clearance and promoting chronic inflammatory responses. When immune complexes containing nuclear antigens (anti-dsDNA, anti-Smith antibodies) accumulate in glomeruli, the alternative pathway amplification loop intensifies local complement deposition, leading to glomerular inflammation and lupus nephritis.

The Evidence

The association between rs1270942 and SLE is among the strongest in the MHC class III region outside of classical HLA alleles. The original SLEGEN GWAS⁵⁵ SLEGEN GWAS
Harley JB et al. 2008 — this was the first large-scale GWAS specifically in women with SLE found per-allele OR=2.35 (p=1.3×10⁻⁵¹) for SLE risk in European ancestry women, with the G allele present at ~13% frequency in European controls rising to ~25% in SLE cases. Subsequent meta-analyses across multi-ancestry cohorts showed p-values approaching 2×10⁻¹⁶⁵, making this one of the most replicated non-HLA associations in lupus genetics.

Functional evidence for CFB's causal role in lupus nephritis is compelling. CFB-deficient MRL/lpr mice⁶⁶ CFB-deficient MRL/lpr mice
Watanabe et al. J Immunol 2000 develop significantly less proteinuria, less glomerular IgG deposition, and lower renal damage scores compared to CFB-sufficient lupus mice, establishing that alternative pathway amplification through factor B drives kidney damage. Mechanistically, review of complement pathway involvement in LN⁷⁷ review of complement pathway involvement in LN
Satyam et al. Transl Res 2022 shows intra-renal upregulation of CFB and complement factor D during renal flares, and elevated alternative pathway biomarkers in renal tissue predict worse renal outcomes⁸⁸ predict worse renal outcomes
Activation at biopsy associated with faster progression to end-stage renal disease independent of classical pathway markers. Furthermore, antisense oligonucleotides that suppress factor B production⁹⁹ suppress factor B production
Grossman et al. 2016 — significant reduction in CFB levels achieved in mice, with dose-dependent improvement in renal scores substantially improve renal pathology in lupus nephritis models, directly validating factor B as a therapeutic target and mechanistic driver.

The G allele shows marked ancestral frequency variation: ~13% in Europeans, dropping to ~5% in Africans and <0.1% in East Asians, mirroring the pattern of SLE susceptibility across these groups. In European women of reproductive age — the demographic at highest absolute SLE risk — the G allele translates into approximately 2-fold increased population risk for SLE and a particularly elevated risk for lupus nephritis when SLE develops.

Practical Implications

Carriers of the G allele — especially GG homozygotes — face substantially elevated SLE risk. SLE predominantly affects women aged 15–45 (female-to-male ratio ~9:1), and the genetics of the MHC class III region including CFB predominantly contribute to this susceptibility. If you carry this variant and have personal or family history of SLE, early rheumatologic evaluation and kidney function monitoring are warranted.

For those with established SLE, the rs1270942 G allele is a marker of alternative pathway-mediated disease — these patients are candidates for complement monitoring (serum C3, C4, alternative pathway activation products like Bb) and may benefit from emerging complement-targeted therapies. The alternative pathway's role in lupus nephritis means that standard immunosuppressants may incompletely suppress renal damage driven by this pathway; CFB inhibitors are under active clinical investigation.

Hydroxychloroquine, the backbone of SLE management, does not directly inhibit complement but reduces flare frequency and may limit the triggers that activate complement cascades. Patients with the G allele who develop lupus nephritis should discuss aggressive management with a nephrologist, including regular urinalysis, 24-hour urine protein quantification, and renal biopsy timing decisions.

Interactions

The rs1270942 G allele operates within the broader SLE genetic architecture centered on the MHC region, where multiple complement genes (C2, C4A, C4B, CFB) form an ancestral MHC haplotype that dramatically increases lupus susceptibility. The extended MHC class III haplotype harboring the CFB G allele is in partial LD with HLA-DR3/DQ2 alleles, making disentanglement of independent effects complex.

For interactions with rs7574865 (STAT4) and rs10516487 (BANK1) — genes in the JAK-STAT and B-cell signaling pathways — the biological pathways are distinct from complement but converge on SLE pathogenesis. Epidemiological data show that individuals carrying risk alleles at multiple SLE loci have additive or supra-additive risk, consistent with the polygenic architecture of SLE. The combination of complement pathway dysregulation (CFB) with T-cell signaling variants (STAT4) or B-cell activation variants (BANK1) may particularly predispose to severe nephritis.

The C3 variant rs2230199 (C3 R102G) is also in this complement pathway family, affecting the downstream amplification substrate. When both CFB and C3 carry risk alleles, the alternative pathway amplification loop may be further dysregulated, though direct interaction data in SLE are limited.

The TYK2 gene encodes a Janus kinase that sits at the junction of the IL-12, IL-23, and type I interferon signaling pathways — the three cytokine cascades most consistently implicated in autoimmune disease. Most protective TYK2 variants studied in depth are coding-sequence changes that impair kinase domain activity directly. rs12720270 operates differently: it is an intronic variant that appears to regulate pre-mRNA splicing¹¹ pre-mRNA splicing
Pre-mRNA splicing is the process by which introns are removed and exons joined to produce mature mRNA; regulated alternative splicing can include or exclude specific exons depending on sequence signals in the surrounding intron/exon context, affecting how TYK2 protein is assembled before it ever reaches the cytoplasm.

The Mechanism

rs12720270 lies in intron 7 of TYK2, positioned 36 nucleotides upstream of the intron 7/exon 8 splice acceptor boundary. Li et al. (2020)²² Li et al. (2020) demonstrated that the minor A allele of rs12720270 (together with rs2304256 in exon 8) promotes the inclusion of exon 8³³ exon 8
Exon 8 of TYK2 encodes part of the FERM domain, which mediates TYK2 binding to cytokine receptor subunits including IFNAR1 and IL-12Rβ1; without exon 8, TYK2 cannot couple properly to these receptors in the mature TYK2 transcript. The G-to-A substitution is predicted in silico to disrupt a potential branch point⁴⁴ branch point
A branch point is an adenosine nucleotide near the 3' end of an intron that initiates the lariat-formation step of RNA splicing; disrupting the branch point sequence can alter the efficiency with which the intron is removed, affecting whether the adjacent exon is included or skipped consensus sequence in intron 7, which would alter the efficiency of exon 8 recognition by the spliceosome.

Exon 8 is not dispensable: it contributes critical residues to the TYK2 FERM domain that mediate receptor binding. Without exon 8, TYK2 cannot efficiently couple to the IL-12 receptor beta-1 subunit (IL-12Rβ1) or the type I interferon receptor chain IFNAR1. Paradoxically, enhancing exon 8 inclusion — and therefore optimizing TYK2 receptor coupling — appears to reduce net autoimmune signaling output, possibly through enhanced regulatory feedback mechanisms or more precisely calibrated cytokine threshold responses.

The rs12720270 A allele and rs2304256 A allele are in strong linkage disequilibrium⁵⁵ linkage disequilibrium
Linkage disequilibrium (LD) is the non-random co-inheritance of alleles at nearby genomic positions; when two variants are in strong LD (high r²), they tend to be co-inherited as a haplotype and it can be difficult to determine which variant drives an observed association in Europeans (r² ≈ 0.83), meaning they are commonly co-inherited. The Li et al. study found that the independent effect of rs12720270 on TYK2 expression disappears after statistically accounting for rs2304256, suggesting the two variants tag the same functional haplotype in European populations. In East Asian populations, allele frequencies and LD patterns differ, which may explain why disease associations are less consistent in Asian cohorts.

The Evidence

The foundational SLE association for rs12720270 comes from two European population studies. Cunninghame Graham et al. (2007)⁶⁶ Cunninghame Graham et al. (2007) studied 380 UK SLE families and found significant under-transmission of the A allele to affected offspring (P = 0.004), placing rs12720270 among the original TYK2 variants implicated in SLE genetics. The location near an intron/exon boundary was noted at the time as mechanistically suggestive of a splicing effect, though the molecular confirmation waited for Li et al. (2020).

The Finnish SLE study⁷⁷ Finnish SLE study (277 patients, 356 controls) by Voss et al. (2009) found significant association at rs12720270 (p = 0.0031), with the G allele conferring risk (OR 1.57, 95% CI 1.16–2.21). This G allele risk framing is the reciprocal of the A allele's protection: carriers of the A allele are protected; carriers of GG lack this protection.

The most comprehensive assessment comes from the Pellenz et al. (2021) systematic review and meta-analysis⁸⁸ Pellenz et al. (2021) systematic review and meta-analysis of 34 studies, which pooled data from 2,792 SLE cases and 5,184 controls specifically for rs12720270. The A allele conferred SLE protection under the allele contrast model in Caucasian populations (FEM OR 0.84, 95% CI 0.72–0.98, P = 0.022), but the effect was not significant in Asian populations and was not evaluated for other autoimmune conditions beyond SLE. This contrasts with the co-associated rs2304256, which has been assessed across eight autoimmune conditions including RA, T1D, and MS.

Practical Implications

For users of European ancestry, the A allele frequency is approximately 18%, meaning about 33% of Europeans carry at least one A allele (AG or AA). The protective signal observed for rs12720270 in SLE is generally attributed to the same exon 8 splicing haplotype as rs2304256, and the two variants should be interpreted together where possible. When a user carries the A allele at both rs12720270 and rs2304256, the evidence is strongest — they carry the complete protective haplotype across both the intronic splicing signal and the exonic Val362Phe variant.

The variant is relevant in the same clinical contexts as other protective TYK2 alleles: autoimmune disease workup interpretation, biologic therapy prescribing context, and family history counseling for SLE and related conditions. As with rs2304256, if a TYK2 inhibitor such as deucravacitinib is prescribed, the A allele at this locus contributes to a genetic background of partially modulated TYK2 signaling.

Interactions

rs12720270 and rs2304256 (Val362Phe) are in strong linkage disequilibrium in Europeans and were identified together in the same functional splicing study. They promote exon 8 inclusion through complementary intronic and exonic splicing signals, respectively — the intronic branch point effect of rs12720270 and the exonic splicing enhancer effect encoded by the Val362Phe substitution of rs2304256. Carriers of the A allele at rs12720270 are very likely to also carry the A allele at rs2304256.

The splicing/expression arm of TYK2 protection (rs12720270 and rs2304256) is mechanistically independent from the kinase-domain-impairment arm (rs34536443 P1104A, rs12720356 I684S, rs35018800 A928V). Individuals who carry the exon 8 splicing haplotype (rs12720270/rs2304256 A alleles) together with coding protective alleles at the other TYK2 loci have multiple independent layers of TYK2 attenuation through distinct structural mechanisms.

Beyond TYK2, this variant acts within the same type I interferon pathway as IRF5 (rs10954213), where functional interaction between TYK2 and IRF5 variants on SLE risk has been documented in Finnish and other European population studies.

Psoriasis affects roughly 2-3% of the global population, and one of the most consequential decisions in its management is determining which patients will go on to develop psoriatic arthritis¹¹ psoriatic arthritis
A chronic inflammatory arthritis affecting peripheral joints, the spine, entheses (tendon/ligament insertions), and nails; occurs in 25-30% of people with psoriasis and can cause irreversible joint damage if untreated. Early PsA rarely announces itself dramatically — joint stiffness, tendon insertion pain, and subtle finger or toe swelling are easily attributed to overuse or aging. By the time the diagnosis is established through X-ray changes, structural damage is often already present. rs12883343 in the regulatory region of NFKBIA is one of the few genetic markers that specifically distinguishes people at elevated risk of PsA from those whose disease will stay skin-limited, making it a valuable early stratification signal.

The Mechanism

NFKBIA (Nuclear Factor Kappa B Inhibitor Alpha) encodes IκB-alpha²² IκB-alpha
The primary cytoplasmic inhibitor of NF-κB; it binds the p65/p50 NF-κB dimer and sequesters it in the cytoplasm, preventing it from entering the nucleus to drive inflammatory gene transcription. Pro-inflammatory signals trigger IκB-alpha phosphorylation, ubiquitination, and proteasomal degradation, releasing NF-κB to activate cytokine genes, the primary cytoplasmic brake on NF-κB inflammatory signaling. Without adequate IκB-alpha, NF-κB dimers enter the nucleus constitutively and drive sustained transcription of TNF-alpha, IL-1beta, IL-6, IL-8, and other cytokines central to joint inflammation and synovial hyperplasia.

rs12883343 sits approximately 18 kb from the NFKBIA gene body in a downstream regulatory region. The variant is not a coding change — it does not alter the IκB-alpha protein sequence — but lies in a region consistent with transcriptional regulatory activity, likely influencing NFKBIA enhancer function or chromatin accessibility. A 2025 single-cell RNA sequencing study³³ 2025 single-cell RNA sequencing study
Garrido et al. scRNAseq of circulating immune cells in PsA vs cutaneous psoriasis and healthy controls found that NFKBIA is overexpressed at the mRNA level in PsA immune cells relative to cutaneous-only psoriasis, but IκB-alpha protein is paradoxically reduced in PsA CD8+ T cells — suggesting translational suppression as an additional regulatory layer. The G allele at rs12883343 likely tags a haplotype that impairs NFKBIA regulatory function in the synovial microenvironment, allowing NF-κB to drive joint-specific inflammatory pathways more effectively than skin-focused inflammatory cascades.

The Evidence

The primary evidence for rs12883343 as a PsA-specific marker comes from a Chinese case-control study⁴⁴ Chinese case-control study
Zhao Q et al. Identification of a Single Nucleotide Polymorphism in NFKBIA with Different Effects on Psoriatic Arthritis and Cutaneous Psoriasis in China. Acta Derm Venereol 2019 enrolling 379 PsA patients, 376 cutaneous psoriasis patients, and 760 healthy controls. The G allele showed a significantly different association profile between the two conditions — its effect on PsA risk (OR=2.371, p=4.93×10⁻¹⁰) was substantially larger than its association with cutaneous-only psoriasis. The extreme statistical significance (10 orders of magnitude below the conventional threshold) despite a modest sample size indicates a robust biological distinction, not a chance finding.

That NFKBIA is a PsA-stratifying locus is further supported by a European cohort study⁵⁵ European cohort study
Coto-Segura P et al. Gene Variant in the NF-kappaB Pathway Inhibitor NFKBIA Distinguishes Patients with Psoriatic Arthritis within the Spectrum of Psoriatic Disease. Acta Derm Venereol 2019 of 690 psoriatic disease patients and 550 controls. Studying a different NFKBIA variant (rs7152376), the rare C allele was more frequent in PsA compared to both controls (OR=2.03, 95% CI 1.3-3.1, p<0.01) and compared to pure cutaneous psoriasis patients (OR=3.2, 95% CI 2.1-5.1, p<0.001). Both Chinese and European cohorts independently converge on NFKBIA as encoding the molecular threshold that separates joint-involving from skin-limited disease.

At a population level, the Stuart et al. GWAS⁶⁶ Stuart et al. GWAS
Stuart PE et al. Genome-wide Association Analysis of Psoriatic Arthritis and Cutaneous Psoriasis Reveals Differences in Their Genetic Architecture. Am J Hum Genet 2015 of over 9,000 European psoriasis cases and 13,670 controls confirmed that NFKBIA achieves genome-wide significance independently for PsA and cutaneous psoriasis, with the two forms of disease showing partially distinct genetic architectures. This establishes NFKBIA not merely as a psoriasis gene, but as a locus where regulatory variation specifically determines the trajectory toward joint involvement.

The clinical stakes are high: PsA is detectable and treatable, but joint damage — erosions, osteolysis, and ankylosis — can occur within the first two years of active disease. Approximately 20% of PsA patients develop severe, disabling joint destruction even with modern treatment. Genetic stratification of which psoriasis patients carry elevated PsA risk enables proactive rheumatological surveillance before that window for prevention closes.

Practical Actions

For carriers of the G allele — particularly GG homozygotes — the clinical imperative is awareness and early symptom recognition, not alarm. The G allele raises relative risk substantially, but most psoriasis patients, regardless of genotype, will not develop severe PsA. What the genetic finding changes is the threshold for seeking rheumatological evaluation. Joint symptoms that a skin-only psoriasis patient might dismiss as overuse deserve prompt attention in a G-allele carrier.

Monitoring for the earliest PsA features — dactylitis (whole-finger swelling, "sausage digit"), enthesitis (pain at tendon insertions, especially Achilles and plantar fascia), and new-onset asymmetric peripheral arthritis — is the most evidence-based response to this genetic signal. Early initiation of DMARDs (methotrexate, sulfasalazine) or biologics (TNF inhibitors, IL-17 inhibitors, IL-23 inhibitors) in confirmed early PsA demonstrably reduces radiographic progression.

NF-κB pathway activity can also be modulated through documented nutritional interventions: high-dose omega-3 fatty acids suppress NF-κB signaling through GPR120 and PPARγ pathways, and vitamin D receptor activation directly induces NFKBIA transcription in immune cells, offering two evidence-based strategies to support the NF-κB brake that this variant may partially impair.

Interactions

NFKBIA rs12883343 sits within the broader NF-κB regulatory network that governs psoriatic disease progression. TNFAIP3/A20 (tagged by rs9321623 and rs5029937) is the other major NF-κB negative regulator in psoriatic disease — it degrades the ubiquitin chains that activate NF-κB upstream of IκB-alpha degradation. Individuals carrying risk alleles at both NFKBIA and TNFAIP3 regulatory loci would have impaired NF-κB suppression through two independent mechanisms, potentially compounding PsA risk.

IL-23R (tagged by rs12044149) contributes to PsA-specific risk through the Th17 axis, which is itself partially NF-κB-dependent. The convergence of NFKBIA regulatory impairment with IL-23R susceptibility alleles may define a high-risk PsA subgroup most likely to benefit from early IL-17 or IL-23 inhibitor therapy.

The color of your eyes is determined primarily by a single nucleotide change on chromosome 15, not in a pigmentation gene itself, but in a regulatory enhancer¹¹ regulatory enhancer
A DNA sequence that controls when and where genes are turned on located deep within intron 86 of the HERC2 gene. This variant, rs12913832, functions as a dimmer switch for the nearby OCA2 gene, which encodes a protein essential for melanin²² melanin
The pigment responsible for eye, skin, and hair color production in the iris. The ancestral A allele permits full OCA2 expression and results in brown eyes, while the derived G allele—which emerged 6,000-10,000 years ago in Europe³³ emerged 6,000-10,000 years ago in Europe
Likely selected during the agricultural transition when lighter pigmentation became advantageous in low-UV environments—reduces OCA2 transcription and produces blue eyes.

The Mechanism

This variant sits within a melanocyte-specific enhancer⁴⁴ melanocyte-specific enhancer
Active only in pigment-producing cells that physically contacts the OCA2 promoter via a long-range chromatin loop spanning 21 kilobases. When you carry the A allele, transcription factors including MITF, LEF1, and HLTF⁵⁵ MITF, LEF1, and HLTF
Master regulators of melanocyte development and function bind efficiently to the enhancer region, pulling it into close proximity with the OCA2 promoter through three-dimensional DNA folding. This chromatin looping⁶⁶ chromatin looping
Physical interaction between distant DNA regions brought together in 3D space dramatically increases OCA2 transcription, leading to robust melanin synthesis and darker eye colors ranging from brown to hazel. The G allele disrupts this process by reducing the binding efficiency of these transcription factors, weakening the chromatin loop formation and decreasing OCA2 expression by approximately 60-70%⁷⁷ 60-70%
Measured in melanocyte cell culture experiments comparing A vs G alleles. With less OCA2 protein available, melanosomes cannot maintain the optimal conditions for melanin production, resulting in the blue structural color that emerges when light scatters through a relatively pigment-free iris stroma⁸⁸ stroma
The fibrous middle layer of the iris.

The Evidence

The association between rs12913832 and eye color is among the strongest genotype-phenotype correlations⁹⁹ genotype-phenotype correlations
Statistical relationships between genetic variants and observable traits in human genetics. A 2008 Danish family study¹⁰¹⁰ A 2008 Danish family study
Eiberg et al. Blue eye color in humans may be caused by a perfectly associated founder mutation in a regulatory element located within the HERC2 gene inhibiting OCA2 expression. Human Genetics, 2008 identified rs12913832 as perfectly associated with blue versus brown eye color across multiple pedigrees, with the GG genotype predicting blue eyes in 99% of cases among Europeans. Visser et al. 2012¹¹¹¹ Visser et al. 2012
HERC2 rs12913832 modulates human pigmentation by attenuating chromatin-loop formation between a long-range enhancer and the OCA2 promoter. Genome Research, 2012 used chromosome conformation capture (3C) technology to directly demonstrate that the G allele weakens the physical interaction between the HERC2 enhancer and OCA2 promoter in melanocytes. This single SNP explains approximately 68% of eye color variance¹²¹² 68% of eye color variance
R²=0.68 for blue-brown eye color differences attributable to this variant in European populations, far exceeding the effect of any other genetic variant. A 2010 Danish population study¹³¹³ 2010 Danish population study
Mengel-From et al. Human eye colour and HERC2, OCA2 and MATP. Forensic Science International: Genetics, 2010 of 395 individuals found that diplotype analysis combining three HERC2 sequence variations with one OCA2 variation yielded a likelihood ratio of 29.3 for predicting light versus dark eye color.

The G allele also has measurable effects beyond the iris. The same chromatin loop mechanism affects melanin production in skin, with GG individuals showing lighter constitutive skin pigmentation¹⁴¹⁴ GG individuals showing lighter constitutive skin pigmentation
Measured by reflectance spectrophotometry on sun-protected skin compared to AA individuals (mean difference of 2-3 units on the melanin index). Gelmi et al. 2025¹⁵¹⁵ Gelmi et al. 2025
Survival in patients with uveal melanoma is linked to genetic variation at HERC2 single nucleotide polymorphism rs12913832. Ophthalmology, 2025 analyzed 392 uveal melanoma patients and found that GG genotype carriers showed significantly worse survival (p=0.017) and higher rates of high-risk tumors with monosomy 3 (p=0.04), though this prognostic effect was mediated through tumor genetics rather than representing an independent factor beyond the chromosomal abnormality. Barón et al. 2014¹⁶¹⁶ Barón et al. 2014
Interactions between ultraviolet light and MC1R and OCA2 variants are determinants of childhood nevus and freckle phenotypes. Cancer Epidemiology, Biomarkers & Prevention, 2014 documented a significant gene-environment interaction where the GG genotype combined with waterside vacations predicted higher total body nevus counts in children ages 6-10, suggesting that blue-eyed children with this genotype may be particularly susceptible to UV-induced melanocytic proliferation.

Practical Implications

If you have the GG genotype (blue eyes), you likely have lower baseline melanin production not just in your irises but also in your skin, which has direct implications for photoprotection¹⁷¹⁷ photoprotection
Natural defense against UV radiation damage. Melanin functions as nature's sunscreen, absorbing UV photons before they can damage DNA in keratinocytes and melanocytes. With reduced melanin, GG individuals face higher melanoma risk compared to AA individuals, with epidemiological studies consistently reporting elevated odds in blue-eyed populations, particularly when combined with intermittent high-intensity sun exposure patterns like beach vacations. The interaction with UV is not simply additive—blue-eyed children show a steeper dose-response curve for nevus development per unit of sun exposure, suggesting a qualitative difference in how their skin responds to UV stress.

The AG genotype produces intermediate effects. A substantial minority of Europeans with AG genotype have intermediate eye colors including green, grey, and hazel, though the majority still have brown eyes. This dosage effect¹⁸¹⁸ dosage effect
One functional copy of the A allele partially rescues OCA2 expression is consistent with the partial restoration of the chromatin loop observed in cell culture studies. Your melanin production is intermediate, and so is your UV sensitivity—higher than AA individuals but lower than GG.

Eye color is polygenic¹⁹¹⁹ polygenic
Controlled by multiple genes with additive effects, and rs12913832 does not explain all variation. Approximately 3% of Europeans with GG genotype have brown eyes due to variants in other pigmentation genes including TYR, TYRP1, SLC24A4, and IRF4. Similarly, some AA or AG individuals have blue eyes due to rare variants discovered through massively parallel sequencing²⁰²⁰ massively parallel sequencing
Next-generation DNA sequencing technology of the OCA2-HERC2 region, including rs191109490 and several others at very low frequencies (0.2-8%). If your eye color doesn't match your rs12913832 genotype, you likely carry one of these modifying variants.

Interactions

This variant is in near-perfect linkage disequilibrium²¹²¹ near-perfect linkage disequilibrium
Two genetic variants inherited together >95% of the time with rs1129038 (r²>0.95), another intronic SNP 29.8 kb away in HERC2, forming a stable haplotype that defines the "blue eye" chromosome in Europeans. The entire 166 kb region spanning from intron 86 of HERC2 through the OCA2 gene shows remarkably low recombination, suggesting strong positive selection²²²² strong positive selection
Evolutionary pressure favoring the blue-eye haplotype in European populations over the past 6,000-10,000 years.

OCA2 itself contains additional functional variants that modify eye color independently of rs12913832. The missense variant rs1800407²³²³ rs1800407
p.Arg419Gln in OCA2 reduces OCA2 protein function directly and is associated with lighter eye colors when combined with rs12913832 AG or GG genotypes. Two other nonsynonymous OCA2 variants, rs74653330 (p.Ala481Thr) and rs121918166 (p.Val443Ile)²⁴²⁴ rs74653330 (p.Ala481Thr) and rs121918166 (p.Val443Ile)
Both reduce OCA2 protein activity, produce blue eyes even in individuals with rs12913832 AA or AG genotypes, demonstrating that impaired OCA2 protein function can override high transcription levels from an intact enhancer.

For melanoma risk stratification, rs12913832 interacts with variants in MC1R, the red hair color gene. The Barón et al. cohort found that sunburns increased larger nevi (≥2mm) specifically in children with both rs12913832 blue-eye alleles and MC1R variants, demonstrating compound gene-UV interaction on melanocytic proliferation markers. This makes biological sense: MC1R variants shift melanin synthesis from protective eumelanin (brown-black) toward pheomelanin (red-yellow), which is not only less photoprotective but may actually generate reactive oxygen species²⁵²⁵ reactive oxygen species
Highly damaging molecules that attack DNA upon UV exposure, compounding the melanin deficiency caused by reduced OCA2 expression.

The variant also shows significant gene-environment interactions²⁶²⁶ gene-environment interactions
Genetic effects that vary depending on environmental exposures with UV exposure patterns. Longitudinal childhood cohort data²⁷²⁷ Longitudinal childhood cohort data
Following children from ages 6-10 with annual sun exposure questionnaires and nevus counts showed that waterside vacations strongly increased total nevus counts specifically in children with rs12913832 blue-eye alleles, while sunburns had a distinct interaction with MC1R variants for larger nevi (≥2mm). This demonstrates that the same sun exposure produces more melanocytic proliferation in blue-eyed children—a qualitatively different UV response with implications for lifelong melanoma risk, given that larger nevi are stronger melanoma precursors than small ones.

The 6q23 chromosomal region harbors one of the strongest known genetic influences on autoimmune disease risk. Most people know of this locus through its risk variants — but rs13207033 tells the other side of the story. This intergenic SNP, located near the TNFAIP3 gene on chromosome 6, tags a protective haplotype¹¹ protective haplotype
A haplotype is a set of alleles inherited together on the same chromosome; rs13207033 marks a chromosomal block associated with higher TNFAIP3 activity that appears to enhance the activity of A20, the immune system's primary brake on NF-kB-driven inflammation. Carriers of the A allele show meaningfully reduced risk for rheumatoid arthritis, particularly in the autoantibody-positive form that causes the most destructive joint disease.

TNFAIP3 encodes A20, a ubiquitin-editing enzyme²² ubiquitin-editing enzyme
A20 uses two functional domains: a deubiquitinase that removes activating ubiquitin chains and an E3 ligase that adds inhibitory chains, together terminating NF-kB signal transduction that terminates NF-kB inflammatory signaling after an immune response has served its purpose. Without sufficient A20 activity, inflammatory signals persist and self-reactive immune responses become more likely — the fundamental mechanism underlying most autoimmune diseases.

The Mechanism

rs13207033 is located in the intergenic region near TNFAIP3 on chromosome 6q23.3, in a ~60-kb stretch between the OLIG3 and TNFAIP3 genes that does not encode any transcripts. It does not alter any protein directly; instead, it tags a regulatory haplotype³³ regulatory haplotype
A haplotype block in high linkage disequilibrium with rs13207033, meaning these alleles travel together across generations that appears to increase TNFAIP3 expression or enhance A20 function. The SNP is in perfect linkage disequilibrium⁴⁴ linkage disequilibrium
LD means the alleles co-occur more frequently than chance would predict, so rs13207033 and rs10499194 are essentially interchangeable markers for the same protective haplotype (r² = 1) with rs10499194, the variant more frequently tested in Asian-population studies. Carriers of the rs10499194 T allele (which perfectly tags the rs13207033 A allele) express higher levels of A20 mRNA, providing a plausible mechanism for the observed protective effect.

The consequence of enhanced A20 expression is more efficient termination of NF-kB signaling cascades. When the immune system encounters an antigen, NF-kB drives cytokine production and T-cell activation. A20 shuts this process down once the threat is cleared. Higher A20 expression means inflammatory responses resolve more completely, reducing the probability that immune activation becomes chronic and self-directed.

The Evidence

The protective nature of rs13207033 was established through conditional logistic regression⁵⁵ conditional logistic regression
A statistical method that identifies whether a variant remains independently significant after accounting for the effects of other nearby variants fine-mapping of the 6q23 locus. After accounting for the known risk variants rs6920220 and rs5029937, rs13207033 remained independently associated with RA protection (OR 0.86, 95% CI 0.80–0.93, P=0.0001) in 3,962 RA patients and 3,531 controls from six UK centers.

The protective effect is most pronounced in ACPA-positive RA⁶⁶ ACPA-positive RA
ACPA stands for anti-citrullinated protein antibodies — autoantibodies that mark the most aggressive, erosive form of rheumatoid arthritis; ACPA-positive RA has stronger genetic determinants than seronegative RA, where the A allele confers OR 0.80 (95% CI 0.66–0.96, P=0.015). This specificity makes biological sense: ACPA-positive RA is a more genetically driven, immune-mediated disease where A20 expression capacity directly constrains disease initiation and severity.

A meta-analysis of six studies⁷⁷ meta-analysis of six studies
Including 11,166 RA cases and 11,231 controls (11,166 cases / 11,231 controls) confirmed the protective association for GA vs GG genotype (OR 0.88, 95% CI 0.79–0.99, P=0.034) with the effect concentrated in Caucasian populations. No significant association was found in Asian or African-American populations, consistent with the lower A allele frequency in East Asian populations (~9.7% vs ~28% in Europeans).

Beyond RA, the rs13207033A–rs10499194T protective haplotype has been associated with reduced risk of ankylosing spondylitis⁸⁸ ankylosing spondylitis
A chronic inflammatory arthritis primarily of the spine and sacroiliac joints, classified as a spondyloarthropathy in Eastern Chinese Han populations (OR 0.60, 95% CI 0.42–0.87, P=0.006) — extending the protective signal beyond RA to other NF-kB-driven inflammatory arthritides.

The three-SNP model at 6q23 illustrates why individual variants tell only part of the story. The most dangerous haplotype combination — homozygous risk alleles at rs6920220 and rs5029937 with complete absence of the rs13207033 protective A allele — carries a combined OR of 1.86⁹⁹ combined OR of 1.86
95% CI 1.51–2.29, from formal three-SNP conditional analysis in 3,962 RA cases for RA, compared to OR 0.86 for the protective A allele alone.

Practical Implications

Carrying one or two copies of the A allele at rs13207033 is genetically favorable — it indicates greater capacity for NF-kB self-regulation via A20. For most A-allele carriers, this is a reassuring counterweight against other inflammatory risk factors. For homozygous GG carriers (lacking the protective allele entirely), the loss of this inherent protective signal is clinically meaningful primarily in the context of the companion 6q23 risk variants; as an isolated finding, GG genotype represents the common baseline population risk.

Two large-scale RCT interventions directly target the NF-kB pathway that A20 regulates. The VITAL trial¹⁰¹⁰ VITAL trial
A 5-year randomized controlled trial of 25,871 adults assigned to vitamin D3 2000 IU/day, omega-3 1 g/day, both, or placebo found vitamin D3 supplementation reduced incident autoimmune disease by 22% (HR 0.78, P=0.05) and omega-3 fatty acids reduced it by 15%. For GG carriers who lack the protective allele, these interventions represent the most evidence-backed approach to compensating for reduced genetic NF-kB braking capacity.

Interactions

rs13207033 is one of three statistically independent signals at the 6q23 locus. The other two — rs6920220 (intergenic, reduces TNFAIP3 transcription) and rs5029937 (intronic in TNFAIP3, also risk) — are independently inherited and their combination with rs13207033 genotype substantially alters overall RA risk. Carrying both risk alleles for the other two SNPs while being GG at rs13207033 (i.e., lacking the protective A allele) yields OR 1.86, the highest risk configuration at this locus.

The TNFAIP3 missense variant rs2230926 (F127C) operates through a different mechanism — impaired A20 enzymatic activity rather than altered expression — and is independent of the 6q23 haplotype block. Carriers of both the rs2230926 G risk allele and rs13207033 GG genotype lack both the quantitative (expression-level) and qualitative (enzyme-activity) components of A20 function, a combination that warrants a dedicated compound action assessment.

CD36 (also called fatty acid translocase, or FAT) is one of the body's key gatekeepers for dietary fat. It sits on the surface of cells throughout the intestine, tongue taste buds, platelets, macrophages, muscle, and adipose tissue. Its job is to bind long-chain fatty acids — the molecules that make up most of the fat in food — and shuttle them into cells. When CD36 expression is higher, cells take up more fat; when it is lower, fat stays in the bloodstream longer. rs13236689 is an intronic variant¹¹ intronic variant
A variant within a non-coding intron of the CD36 gene; these often act as expression regulators by altering transcription factor binding sites in the surrounding DNA that modulates how much CD36 protein cells produce.

The Mechanism

rs13236689 does not change the CD36 protein sequence — it lies within an intron on chromosome 7 at position 80,606,698 (GRCh38). Its biological significance lies in its status as a platelet expression quantitative trait locus (eQTL)²² platelet expression quantitative trait locus (eQTL)
An eQTL is a genomic variant that explains variation in how much of a gene's mRNA — and ultimately protein — is produced in a given tissue. Carriers of the G allele are in strong linkage disequilibrium with the functional regulatory variants rs2366739 and rs1194196, which sit approximately 13–55 kb upstream of the CD36 transcriptional start site and alter transcription factor binding. When these regulatory variants are active, platelet CD36 mRNA and surface protein levels increase. Higher platelet CD36 means stronger responses to oxidized LDL (oxLDL) — the chemically modified form of cholesterol that accumulates in atherosclerotic plaques — and greater platelet activation in lipid-rich environments.

The metabolic consequences extend beyond platelets. In the gut, enterocyte CD36 expression levels determine how efficiently long-chain fatty acids trigger satiety signals: CD36 converts dietary fatty acids into oleoylethanolamide (OEA)³³ oleoylethanolamide (OEA)
A bioactive lipid that activates PPARα in the gut wall, sending a sustained satiety signal to the brain via the vagus nerve, which in turn activates PPARα to prolong the feeling of fullness after a fat-rich meal. Variants that push CD36 expression higher also increase uptake of oxLDL by macrophages in artery walls — a key step in foam cell formation and atherosclerosis.

The Evidence

The primary genetic evidence for rs13236689 comes from the COGENT consortium meta-analysis⁴⁴ COGENT consortium meta-analysis
Cheng et al. (2012), PLoS Genetics — a meta-analysis of 7 genome-wide association studies in over 16,000 African Americans from population-based cohorts, which identified rs13236689 as a genome-wide significant locus for platelet count (p = 2.84×10⁻⁹; β = +4.18 × 10⁹ platelets/L per G allele). The association nominally replicated in European Americans. The same variant was subsequently shown by Madan et al. (2019, PLoS Genetics)⁵⁵ Madan et al. (2019, PLoS Genetics)
Using a massively parallel reporter assay to screen 81 CD36 eQTLs, the authors confirmed that rs13236689 is a bona fide eQTL for platelet CD36 mRNA, with the causal signal mapping to two nearby regulatory variants (rs2366739, rs1194196) in high LD to be a platelet CD36 eQTL: G allele carriers have higher platelet CD36 surface expression.

The downstream cardiovascular significance of CD36 expression variation is well-established: higher platelet CD36 augments oxLDL-induced platelet activation and has been associated with thromboembolism risk. At the gene level, Love-Gregory et al. (2008)⁶⁶ Love-Gregory et al. (2008)
Population study of 2,020 African Americans from the HyperGEN cohort; PMID 18305138 showed that multiple CD36 variants modulate HDL-C, triglycerides, and metabolic syndrome risk. The CD36 AGGIG haplotype in Caucasians was associated with 31% higher free fatty acids and OR 2.3 for cardiovascular events in type 2 diabetics (Corpeleijn et al. 2006, PMID 15282206)⁷⁷ (Corpeleijn et al. 2006, PMID 15282206).

Evidence for rs13236689's direct effect on dietary fat handling and plasma lipids is moderate — the lipid associations are inferred from the eQTL relationship and the well-established biology of CD36 expression, rather than from dedicated dietary intervention trials targeting this specific variant.

Practical Actions

Because rs13236689 G allele carriers have higher platelet CD36 expression, they are likely to mount stronger platelet responses to circulating oxidized LDL. This is most relevant in the context of a diet high in saturated and trans fats, which drive oxLDL production. Limiting foods that generate oxLDL — principally processed and ultra-processed foods containing oxidized vegetable oils and trans fats — is directly relevant to this genotype. Longer-chain omega-3 fatty acids (EPA/DHA) compete with pro-inflammatory fatty acids at the CD36 binding site and have been shown to reduce platelet activation. Monitoring fasting triglycerides and HDL-C provides a practical window into CD36-mediated fat metabolism; the CD36 eQTL mechanism links directly to postprandial lipemia.

Interactions

rs13236689 is in linkage disequilibrium with other CD36 expression-regulating variants. The most studied CD36 SNP, rs1761667 (promoter, −31118G>A), directly affects fat taste sensitivity — AA homozygotes have 8-fold higher detection thresholds for oleic acid, meaning they perceive fat less intensely and tend to eat more of it. If a person carries both rs13236689 G (higher CD36 in platelets/ macrophages) and rs1761667 AA (lower fat taste sensitivity from reduced tongue CD36), the combined effect — increased cardiovascular risk from platelet hyperactivation alongside blunted dietary fat sensing — warrants attention. rs3211938 is a CD36 truncation variant (stop-gain) that dramatically reduces CD36 protein; this variant is functionally distinct from rs13236689 and largely confined to African-ancestry populations.

When a macrophage migrates into a lipid-laden arterial wall and begins absorbing oxidized cholesterol, it faces a critical choice: offload that cargo to HDL particles and survive, or become overwhelmed by cholesterol accumulation and transform into a foam cell. This offloading depends heavily on the ABCG1 transporter¹¹ ABCG1 transporter
ATP-binding cassette subfamily G member 1, a membrane pump that mediates cholesterol and phospholipid efflux from macrophages to mature HDL particles — the rate-limiting step in macrophage reverse cholesterol transport. When macrophage ABCG1 function is compromised, cholesterol accumulates, foam cells form, and the necrotic core of atherosclerotic plaques expands — a process that is directly linked to plaque instability and rupture, the proximate cause of most atherothrombotic strokes.

The rs1378577 variant sits approximately 2 kilobases upstream of the ABCG1 transcription start site on chromosome 21. Located in a regulatory region that influences ABCG1 expression, this T>G substitution is part of the promoter-region haplotype that has been studied in relation to ischemic stroke in Asian populations. The G allele at rs1378577 appears to be associated with reduced stroke risk — particularly the atherothrombotic subtype, which is characterized by large-artery atherosclerosis and plaque-related events.

The Mechanism

As an upstream regulatory variant, rs1378577 likely influences the binding of transcriptional regulators to the ABCG1 locus, modulating expression levels rather than altering the protein sequence directly. The ABCG1 gene on chromosome 21 is expressed widely, with particularly important roles in macrophages and monocytes that reside within arterial walls.

ABCG1 facilitates the transfer of cholesterol from macrophage plasma membranes to mature, spherical HDL particles²² ABCG1 facilitates the transfer of cholesterol from macrophage plasma membranes to mature, spherical HDL particles
this step follows ABCA1-mediated lipidation of apolipoprotein A-I into nascent HDL; the two transporters work sequentially — ABCA1 initiates the HDL particle, ABCG1 enlarges it. Reduced ABCG1 expression impairs this efflux step, leaving macrophages with excess intracellular cholesterol that promotes foam cell differentiation. ABCG1 deficiency also amplifies Toll-like receptor (TLR4) signaling in macrophages, heightening inflammatory cytokine release — a mechanism that independently accelerates plaque vulnerability independent of cholesterol level effects.

Variants in the ABCG1 promoter region that reduce transcriptional activity could therefore impair reverse cholesterol transport at the arterial wall specifically, contributing to plaque growth and instability. The nearby promoter variant rs57137919 has been confirmed by luciferase assay to alter ABCG1 promoter activity; rs1378577 is part of the same upstream regulatory haplotype and is in linkage disequilibrium with the promoter region variants studied in Chinese Han populations.

The Evidence

The primary study of rs1378577 is Li et al. 2015 (Journal of Stroke and Cerebrovascular Diseases)³³ Li et al. 2015 (Journal of Stroke and Cerebrovascular Diseases), a case-control study of 389 ischemic stroke patients and 380 healthy controls in a Chinese Han population. In the overall cohort, no significant association was detected — but when stroke was subtyped, the G allele and TG+GG genotypes of rs1378577 were specifically associated with reduced risk of atherothrombotic stroke, the subtype driven by large-artery plaque rupture rather than cardioembolic or lacunar mechanisms. The association was particularly pronounced in the hypertriglyceridemic subgroup (144 cases, 115 controls), where the GG genotype was markedly less frequent among stroke patients, suggesting the G allele's protection is most evident when metabolic stress on macrophage lipid handling is highest.

Yang et al. 2022 (Gene)⁴⁴ Yang et al. 2022 (Gene) studied 10 ABCA1/G1 SNPs in 249 ischemic stroke patients and 226 controls, confirming that ABCG1 variants are associated with plasma lipid differences and stroke susceptibility in this population. The studies are consistent in pointing to ABCG1's role in HDL-mediated cholesterol transport as the biological substrate linking genotype to stroke risk.

Wang et al. 2020 (Annals of Vascular Surgery)⁵⁵ Wang et al. 2020 (Annals of Vascular Surgery) provided functional context by showing that ABCG1 promoter polymorphisms are associated with plasma HDL-C and LDL-C differences in Chinese Han individuals, consistent with the transporter's known role in systemic lipoprotein metabolism.

Important limitation: all published studies are from Chinese Han populations. Effect size estimates may differ in European, African, or South Asian ancestry groups. The G allele is notably more common in African populations (~50%) than European (~23%), which may reflect population-specific selection pressures and could affect how absolute risk differences translate across ancestries. Mechanistic evidence for ABCG1's role in macrophage cholesterol efflux and atherosclerosis is well-established across populations; the variant-level association requires replication in non-Asian cohorts.

Evidence level is moderate: replicated in multiple Chinese Han studies with a biologically plausible mechanism, but replication in diverse non-Asian populations is incomplete, and exact odds ratios from the abstract-available data are limited.

Practical Actions

For individuals carrying the T/T genotype (most common, ~58% globally): HDL-mediated cholesterol clearance from macrophages operates at baseline levels without the G-allele influence. Supporting ABCG1-dependent reverse cholesterol transport through specific dietary strategies — particularly increasing HDL particle function and reducing LDL oxidation exposure at arterial walls — is the most directly relevant action for this genotype.

For G/T heterozygotes: partial G-allele benefit may improve macrophage cholesterol efflux capacity and lipid profiles to an intermediate degree. Monitoring lipids and tracking HDL-C trends confirms whether the expected benefit is expressed.

For GG homozygotes: the G allele's protective association with atherothrombotic stroke is most pronounced in this genotype. This is particularly relevant in the context of elevated triglycerides, where ABCG1-dependent efflux faces the greatest metabolic challenge.

Interactions

rs1378577 is part of the same genomic region as rs57137919, an ABCG1 promoter variant associated with HDL-C, LDL-C, and CAD/stroke risk. Both variants are in linkage disequilibrium in Chinese Han populations, and their effects on ABCG1 expression are likely partially overlapping. Users carrying G alleles at both loci may experience compounded ABCG1 expression effects in macrophages.

ABCG1 works sequentially with ABCA1 in macrophage cholesterol efflux — ABCA1 (rs4149338) creates nascent HDL particles, ABCG1 (rs1378577) then loads cholesterol onto those particles. The Li et al. 2015 study examined both genes: ABCA1 variants showed opposing risk directions compared to ABCG1 variants, highlighting that dysfunction at either step in the efflux cascade can impair reverse cholesterol transport. The two-gene combination is a candidate for compound interaction analysis.

HMGA1 (High Mobility Group AT-Hook 1) encodes a non-histone chromosomal protein that functions as a transcriptional regulator¹¹ transcriptional regulator
HMGA1 binds to A/T-rich regions of DNA and remodels chromatin structure, enabling other transcription factors to access their target genes. It acts as an architectural transcription factor — rather than directly activating genes, it reshapes the DNA landscape to allow or prevent other regulators from doing their work. The rs1776897 variant sits in an intergenic region near HMGA1 on chromosome 6 and influences where your body preferentially stores fat.

The Mechanism

HMGA1 plays a documented role in insulin signaling and glucose metabolism. Mouse studies²² Mouse studies
HMGA1-deficient mice develop obesity, glucose intolerance, and insulin resistance, demonstrating a causal role for this gene in metabolic regulation have shown that HMGA1 deficiency leads to obesity, glucose intolerance, and insulin resistance. The protein regulates the expression of the insulin receptor itself, and reduced HMGA1 function leads to decreased insulin receptor expression on cell surfaces.

The rs1776897 G allele is associated with altered HMGA1 regulatory activity in adipose tissue, promoting central (abdominal) fat deposition rather than peripheral storage. This variant shows marked sexual dimorphism³³ sexual dimorphism
The effect on WHR is substantially stronger in women than in men, a pattern shared by many fat distribution loci, with stronger effects in women.

The Evidence

The Shungin et al. 2015 meta-analysis⁴⁴ Shungin et al. 2015 meta-analysis
Shungin et al. New genetic loci link adipose and insulin biology to body fat distribution. Nature, 2015 identified the HMGA1 locus among 49 genome-wide significant loci for waist-to-hip ratio adjusted for BMI (WHRadjBMI) across 224,459 individuals. HMGA1 was specifically highlighted as one of the transcriptional regulators at WHRadjBMI loci, and the signal showed stronger effects in women.

A Mendelian randomization study⁵⁵ Mendelian randomization study
Emdin et al. Genetic association of waist-to-hip ratio with cardiometabolic traits, type 2 diabetes, and coronary heart disease. JAMA, 2017 demonstrated that WHR-raising variants, including the HMGA1 locus, are causally linked to increased risk of type 2 diabetes (OR 1.77 per 1-SD increase in WHR) and coronary heart disease (OR 1.46). This moves the evidence beyond mere association — the fat distribution pattern itself drives metabolic disease.

The largest body fat distribution GWAS to date⁶⁶ largest body fat distribution GWAS to date
Pulit et al. Meta-analysis of genome-wide association studies for body fat distribution in 694,649 individuals of European ancestry. Hum Mol Genet, 2019 confirmed the HMGA1 locus among replicated fat distribution signals across nearly 700,000 individuals.

Practical Actions

The HMGA1 variant's effect on central fat deposition means that waist-to-hip ratio is a more meaningful health metric than BMI or total weight for carriers of the G allele. Central adiposity, even without overall obesity, carries metabolic risk. Monitoring insulin sensitivity markers and focusing on reducing visceral fat specifically are appropriate responses.

Interactions

HMGA1 sits within a broader network of fat distribution loci. The Shungin et al. 2015 pathway analysis implicated adipogenesis, angiogenesis, transcriptional regulation, and insulin resistance as interconnected processes. Carriers of WHR-raising alleles at multiple loci (including TFAP2B rs987237 and VEGFA rs6905288) may have compounded central adiposity risk, though formal gene-gene interaction studies are limited for these specific combinations.