SLC19A1 (Solute Carrier Family 19 Member 1), also known as the reduced folate carrier (RFC1), is the primary transporter responsible for moving folate from your blood into your cells. Even if you produce adequate methylfolate (via MTHFR) or take methylfolate supplements, this transporter determines how efficiently that folate actually reaches the inside of your cells where it is needed.

The Mechanism

The G80A variant (rs1051266) causes a histidine-to-arginine substitution ¹¹ Histidine-to-arginine substitution at position 27 of the transporter protein (p.His27Arg) at position 27 of the transporter protein, located in transmembrane domain 1 (TMD1), a region implicated in substrate binding and translocation. The T allele (arginine variant) has altered transport kinetics, resulting in reduced folate uptake into cells. This creates a situation where blood folate levels may appear normal on a standard test, but intracellular folate levels are suboptimal.

Clinical Relevance

This variant is particularly important in the context of other methylation variants. If you have reduced MTHFR activity (making less methylfolate) AND reduced RFC1 transport (getting less folate into cells), the combined effect can be more significant than either variant alone. Studies have also linked this variant to altered methotrexate response ²² Methotrexate is an antifolate drug used for cancer and autoimmune diseases — it competes with folate for the same RFC1 transporter, since methotrexate uses the same transporter. A PharmGKB summary³³ A PharmGKB summary
Gong L et al. SLC19A1 Pharmacogenomics Summary, 2010 documents the pharmacogenomic relevance of this transporter.

The Bigger Picture

The folate pathway is like a production line: MTHFR converts folate to its active form, SLC19A1 transports it into cells, and MTHFD1 helps process it further. Bottlenecks at any step can reduce overall methylation capacity ⁴⁴ This is why looking at individual SNPs in isolation can be misleading — the whole pathway matters. By understanding which steps are compromised, you can target your supplementation more effectively.

Practical Implications

If you carry the T allele, ensuring adequate (or slightly above average) folate intake becomes important. Methylfolate may have an advantage over folic acid since it is already in the active form and may be transported more efficiently. Higher doses may help compensate for reduced transport efficiency.

Interactions

SLC19A1 interacts with MTHFR (rs1801133, rs1801131) — if both folate production and transport are impaired, the combined effect is greater. It also interacts with MTHFD1 (rs2236225) for downstream folate processing.

Every heartbeat demands that cardiomyocytes — the muscle cells of the heart — transmit enormous mechanical forces across their shared boundaries. This is the job of the desmosome¹¹ desmosome
A specialized cell-cell junction found in tissues under mechanical stress. The cardiac desmosome anchors the cytoskeletal intermediate filaments of adjacent cardiomyocytes through a multi-protein complex at the intercalated disc, and its master scaffold is desmoplakin (DSP) — the only structural component that spans from the desmosome's cytoplasmic plaque to the keratin and desmin filaments inside the cell. When desmoplakin fails, the result is not immediate structural collapse but a slow, cumulative catastrophe: cardiomyocytes detach, die, and are replaced by fatty-fibrous scar tissue that creates an electrical substrate for lethal arrhythmias.

The DSP p.Arg451Gly variant replaces a positively-charged arginine at position 451 with the smallest, most conformationally flexible amino acid — glycine — within the N-terminal globular domain that anchors desmoplakin to plakoglobin and the desmosomal plaque. Unlike truncating DSP mutations (which reduce the total amount of functional desmoplakin through haploinsufficiency), this missense variant produces a full-length protein that is selectively vulnerable to enzymatic degradation. The result is a classified Pathogenic variant²² classified Pathogenic variant
ClinVar VCV000948761, reviewed by Labcorp Genetics, GeneDx, and Ambry Genetics with no conflicts, for arrhythmogenic cardiomyopathy with woolly hair and keratoderma and arrhythmogenic right ventricular dysplasia 8 associated with biventricular arrhythmogenic cardiomyopathy (ACM).

The Mechanism

Ng et al. (JCI Insight 2019)³³ Ng et al. (JCI Insight 2019)
Patient mutations linked to arrhythmogenic cardiomyopathy enhance calpain-mediated desmoplakin degradation. JCI Insight 4:e128, 2019 identified the p.Arg451Gly substitution in a family with biventricular ACM and established its pathogenic mechanism through molecular modeling and functional assays. Molecular dynamics simulations revealed that the Arg→Gly substitution causes a "significant loss of intramolecular interactions surrounding a putative calpain⁴⁴ calpain
A calcium-activated protease that cleaves proteins at specific recognition sequences. Calpain activity is elevated in stressed cardiomyocytes and normally degrades damaged contractile proteins; excessive calpain activity at the desmosome destabilizes the intercalated disc complex target site" in the DSP N-terminal domain.

The consequence is selective proteolytic vulnerability: calpain degrades the Arg451Gly mutant desmoplakin faster than wild-type protein, reducing the effective desmoplakin concentration at the intercalated disc under conditions of mechanical or calcium stress. Patient-derived induced pluripotent stem cell (iPSC) cardiomyocytes confirmed substantially reduced connexin-43 levels at intercalated discs — reflecting broader intercalated disc remodeling downstream of desmoplakin deficiency — despite normal electrical propagation at baseline. In silico screening by the same team identified three additional ACM desmoplakin mutations (S299R, S442F, S507F) sharing enhanced calpain susceptibility, suggesting this mechanism is a recurrent theme across N-terminal domain missense variants rather than unique to p.Arg451Gly.

The Evidence

The initial identification of DSP p.Arg451Gly was in a family with biventricular ACM with a LOD score of 7.65⁵⁵ family with biventricular ACM with a LOD score of 7.65
LOD (logarithm of the odds) score: a statistical measure of genetic linkage. LOD ≥ 3 is conventionally considered significant; 7.65 is exceptionally strong evidence for co-segregation of the variant with disease in the family, establishing near-certain causality. The variant is absent from gnomAD population databases (1 carrier in 805,808 alleles globally), consistent with strong negative selection against a high-penetrance cardiac disease variant.

The broader DSP cardiomyopathy literature provides robust context for clinical risk. Smith et al. (Circulation 2020)⁶⁶ Smith et al. (Circulation 2020)
Desmoplakin cardiomyopathy, a fibrotic and inflammatory form of cardiomyopathy distinct from typical dilated or arrhythmogenic right ventricular cardiomyopathy. Circulation 141:1872–1884, 2020 studied 107 patients with pathogenic DSP variants and found that left ventricular (LV) predominant disease occurred in 55% (versus 0% for PKP2 mutations; p<0.001), subepicardial late gadolinium enhancement (LGE) on MRI in 40%, and preserved LVEF in 35% of those with imaging-positive fibrosis — establishing that fibrosis precedes and is disproportionate to systolic dysfunction. This is the defining challenge of DSP cardiomyopathy: standard echocardiography misses the fibrotic substrate that underlies most arrhythmic events.

Wang et al. (Europace 2022)⁷⁷ Wang et al. (Europace 2022)
Clinical characteristics and risk stratification of desmoplakin cardiomyopathy. Europace 24:268–277, 2022 followed 91 DSP variant carriers for a median 4.3 years and found sustained ventricular arrhythmia (VA) at 5.9 per 100 person-years and heart failure at 6.7 per 100 person-years. Acute myocardial injury episodes — episodes of chest pain, troponin elevation, and new MRI changes without obstructive coronary disease — affected 22% of participants and were the strongest predictors of subsequent arrhythmia and heart failure.

Gasperetti et al. (Eur Heart J 2025)⁸⁸ Gasperetti et al. (Eur Heart J 2025)
Clinical features and outcomes in carriers of pathogenic desmoplakin variants. Eur Heart J 46:362–376, 2025 enrolled 800 DSP pathogenic variant carriers across 26 institutions. Over 3.7 years, 17.4% developed sustained VA (3.9%/year annual rate) and 9.0% required heart failure hospitalization. Independent VA predictors included female sex (aHR 1.547), prior non-sustained VT (aHR 1.721), prior sustained VA (aHR 1.923), LVEF ≤50% (aHR 1.645), and myocardial injury episodes (HR 2.394). Critically, 32.5% of carriers met no conventional diagnostic criteria for ARVC, DCM, or NDLVC — underscoring that negative standard phenotypic screening cannot exclude genotype-positive DSP cardiomyopathy.

The DSP Risk Score (Carrick et al., Eur Heart J 2024)⁹⁹ DSP Risk Score (Carrick et al., Eur Heart J 2024)
A novel tool for arrhythmic risk stratification in desmoplakin gene variant carriers. Eur Heart J 45:2968–2979, 2024 provides a validated five-predictor model (female sex HR 1.9, NSVT history HR 1.7, 24-hour PVC burden, LVEF <50% HR 1.5, moderate-to-severe RV dysfunction HR 6.0) with c-statistic 0.782 in the development cohort and 0.791 on external validation. The score stratifies patients into low (<5%), intermediate (5–20%), and high-risk (>20%) 5-year VA probability categories with 100% negative predictive value in the low-risk group.

Practical Actions

Carrying a pathogenic DSP missense variant changes clinical management fundamentally — not because disease is certain, but because the risk is high enough and the window for prevention wide enough to warrant structured surveillance. The goals are: (1) quantify fibrosis burden before systolic dysfunction develops, (2) identify arrhythmia risk predictors while time permits ICD planning, and (3) protect first-degree relatives through cascade genetic testing.

Cardiac MRI with late gadolinium enhancement is non-negotiable — echocardiography cannot detect the subepicardial fibrosis that is the defining pathological substrate of DSP cardiomyopathy. Ambulatory ECG monitoring (at least 14 days) quantifies PVC burden, a direct DSP Risk Score input. Competitive sport restriction is recommended for confirmed disease, as exercise-associated myocardial stress can trigger "hot phase" episodes that deposit additional fibrosis and elevate downstream arrhythmia risk.

Interactions

DSP p.Arg451Gly shares a pathogenic mechanism — enhanced calpain-mediated N-terminal domain degradation — with at least three other DSP missense variants (S299R, S442F, S507F) identified in the same functional screen. Compound or digenic desmosomal genotypes (DSP plus PKP2, DSG2, or JUP variants) have been reported with more severe biventricular phenotypes and earlier disease onset, though specific DSP Arg451Gly compound data are not published. SCN5A sodium channel variants may modify arrhythmia penetrance as electrophysiological modifiers in desmosomal cardiomyopathy families.

Neurocan¹¹ Neurocan
NCAN: a chondroitin sulfate proteoglycan expressed almost exclusively in the central nervous system. It is a major component of the brain's extracellular matrix — the molecular scaffolding that supports and organises neurons and their connections (NCAN) is a structural protein in the brain's extracellular matrix that shapes how neurons grow, migrate, and form connections during development and throughout life. rs1064395 is a 3' UTR variant²² 3' UTR variant
A variant in the 3' untranslated region of the mRNA. This region controls mRNA stability, localization, and translation efficiency — making 3' UTR variants functionally important even though they don't change the protein sequence in the NCAN gene — a single-letter change in the messenger RNA that likely alters how much neurocan is produced in brain tissue. It was identified in a genome-wide association study as a significant susceptibility factor for bipolar disorder and has since been independently replicated across multiple populations and extended to schizophrenia, with neuroimaging studies showing measurable effects on hippocampal structure and memory function even in healthy individuals.

The Mechanism

Neurocan is one of the dominant chondroitin sulfate proteoglycans³³ chondroitin sulfate proteoglycans
CSPGs: a family of extracellular matrix proteins built around a protein core decorated with long chains of sulphated sugars. In the brain, CSPGs form "perineuronal nets" around inhibitory neurons and regulate synaptic plasticity by controlling which connections can form or remodel (CSPGs) in the developing and adult brain. It is highly expressed in the cortex and hippocampus — precisely the regions most relevant to memory, mood regulation, and psychiatric vulnerability. Neurocan modulates axon guidance during neural development, restricts aberrant synaptic remodelling in adults, and interacts with other extracellular matrix proteins to maintain the structural integrity of perineuronal nets⁴⁴ perineuronal nets
Dense lattices of extracellular matrix molecules that wrap around the cell bodies and proximal dendrites of certain neurons, especially fast-spiking GABAergic interneurons. Perineuronal nets consolidate synaptic connectivity and regulate critical periods of brain development.

rs1064395 falls in the 3' UTR of NCAN, which does not alter the neurocan protein sequence but likely affects mRNA stability or translational efficiency. The net functional consequence is presumed to be altered neurocan protein levels in neuronal tissue — though the exact molecular mechanism has not yet been fully characterised. Expression quantitative trait locus (eQTL) evidence from brain tissue supports a regulatory effect, consistent with the observed dose-dependent changes in hippocampal and amygdala structure with each copy of the A allele.

The Evidence

The original GWAS by Cichon and colleagues⁵⁵ original GWAS by Cichon and colleagues
Cichon S et al. Genome-wide association study identifies genetic variation in neurocan as a susceptibility factor for bipolar disorder. Am J Hum Genet, 2011 identified rs1064395 as a genome-wide significant locus for bipolar disorder in a discovery sample of 2,411 patients and 3,613 controls (OR 1.31, p = 3.02×10⁻⁸), and replicated this in 6,030 patients and 31,749 controls (OR 1.12). The meta-analysis yielded OR 1.17 and p = 2.14×10⁻⁹. A subsequent meta-analysis of 15,318 cases and 91,990 controls⁶⁶ meta-analysis of 15,318 cases and 91,990 controls
Wang L et al. Further evidence of an association between NCAN rs1064395 and bipolar disorder. Mol Neuropsychiatry, 2018 confirmed genome-wide significance (A allele OR 1.126, p = 4.92×10⁻⁹).

Mühleisen et al.⁷⁷ Mühleisen et al.
Mühleisen TW et al. Association between schizophrenia and common variation in neurocan (NCAN), a genetic risk factor for bipolar disorder. Schizophr Res, 2012 extended the association to schizophrenia across 5,061 patients and 9,655 controls (A-allele OR 1.11, p = 2.28×10⁻³), establishing rs1064395 as a shared cross-disorder psychiatric risk variant. The effect size is modest — comparable to other common psychiatric GWAS hits — consistent with the polygenic architecture of both disorders.

Neuroimaging provides the most direct window into how this variant affects the brain. Dannlowski et al.⁸⁸ Dannlowski et al.
Dannlowski U et al. NCAN cross-disorder risk variant is associated with limbic gray matter deficits in healthy subjects and major depression. Neuropsychopharmacology, 2015 found that A-allele carriers had reduced gray matter volume in the amygdala and hippocampus in both 512 healthy subjects and 171 depressed inpatients — a pattern that closely mirrors the structural changes seen in bipolar disorder. The Assmann et al.⁹⁹ Assmann et al.
Assmann A et al. Neurocan genome-wide psychiatric risk variant affects explicit memory performance and hippocampal function in healthy humans. Eur J Neurosci, 2021 functional MRI study documented reduced verbal recall and elevated false alarm rates on a recognition memory task in A-allele carriers across two independent cohorts (N=572 and N=302), with fMRI showing inefficiently increased left hippocampal activation in risk-allele carriers — a sign of compensatory over-recruitment in the face of reduced neural efficiency.

Practical Implications

The absolute risk increase from a single copy of the A allele is modest. This is a common variant with an odds ratio around 1.17 — comparable to many other GWAS-identified psychiatric risk variants. The majority of A-allele carriers never develop bipolar disorder or schizophrenia. What makes the variant notable is that it also produces measurable, subclinical effects on hippocampal structure and memory performance in the general population, suggesting a neurobiological mechanism that exists on a continuum rather than as a threshold effect.

The neuroimaging data argue for paying attention to hippocampal health: sleep quality, aerobic exercise, and omega-3 fatty acid intake are the lifestyle factors with the strongest evidence for maintaining hippocampal volume and function.

Interactions

NCAN rs1064395 has not been systematically studied in combination with other specific psychiatric risk SNPs in compound heterozygosity analyses. However, since both bipolar disorder and schizophrenia are highly polygenic, the effect of rs1064395 is best understood in the context of total polygenic risk — multiple small-effect variants accumulating to meaningful susceptibility. Individuals who carry several independently identified psychiatric risk variants (e.g. in CACNA1C, ANK3, DISC1, or other GWAS-significant loci) alongside the NCAN A-allele may have a more substantially elevated personal risk profile than any single variant implies.

CYP2D6 is responsible for metabolizing roughly 25% of all prescribed medications, yet the amount of CYP2D6 enzyme your liver produces is not fixed — it varies with a cluster of regulatory variants that act like volume controls on the gene. One of the most consistently studied of these is rs1080985¹¹ rs1080985
the -1584C>G polymorphism located in the CYP2D6 promoter region, approximately 1,584 base pairs upstream of the coding sequence, a nucleotide change that alters how efficiently the gene is transcribed into protein. The C allele at this position is associated with reduced promoter activity — meaning cells produce less CYP2D6 enzyme, and drug metabolism is correspondingly slower.

The Mechanism

The GRCh38 reference allele at rs1080985 is G (the most common allele in European and global populations). The minor C allele is found in roughly 17% of Europeans but is nearly absent in East Asian and African populations. This strong ancestry-stratification is itself a clue about selective pressure on CYP2D6 expression levels.

A landmark 2001 study using 76 human liver biopsies²² A landmark 2001 study using 76 human liver biopsies
Zanger et al. Comprehensive analysis of the genetic factors determining expression and function of hepatic CYP2D6. Pharmacogenetics, 2001 established that individuals carrying the G allele at position -1584 had consistently higher CYP2D6 protein levels and enzyme activity compared to C allele carriers. The G allele functions as a stronger promoter, driving higher transcription rates. The C allele, by contrast, creates a less efficient promoter, resulting in reduced CYP2D6 expression — even when the coding sequence itself is intact.

This promoter variant is closely linked — though not identical — to the CYP2D6*2A haplotype, which combines this promoter change with coding variants including rs16947 (Arg296Cys). The rs1080985 variant is also in high linkage disequilibrium with the CYP2D6 enhancer variant rs5758550, located over 100 kb downstream. The net functional outcome therefore depends on the full haplotype context, but the -1584C allele consistently tracks with reduced expression in liver expression studies.

The Evidence

The strongest clinical evidence for rs1080985 comes from studies of donepezil — a cholinesterase inhibitor used to treat Alzheimer's disease that is substantially metabolized by CYP2D6. Because donepezil requires an active enzyme for clearance, lower CYP2D6 activity (C allele) was hypothesized to result in higher drug accumulation and, counterintuitively, better therapeutic response — or alternatively, altered exposure patterns affecting efficacy.

A prospective cohort of 127 Alzheimer's patients³³ A prospective cohort of 127 Alzheimer's patients
Pilotto et al. Effect of a CYP2D6 polymorphism on the efficacy of donepezil in patients with Alzheimer disease. Neurology, 2009 found that G allele carriers (normal/higher CYP2D6 expression) had a significantly lower likelihood of poor donepezil response — OR 3.43 (95% CI 1.49–7.90, p=0.013) — meaning C allele carriers were over three times more likely to respond poorly to the drug. This landmark finding was replicated in a larger study of 415 patients (OR 1.74, p=0.04).

A 2016 meta-analysis of 1,266 donepezil-treated patients⁴⁴ A 2016 meta-analysis of 1,266 donepezil-treated patients
Xiao et al. Effect of the CYP2D6 and APOE Polymorphisms on the Efficacy of Donepezil. CNS Drugs, 2016 confirmed that rs1080985 C allele carriers have a significantly elevated risk of poor donepezil response, and found that carrying both the C allele and APOE-ε4 further increased non-response risk (OR 1.73, p=0.03).

A smaller pharmacokinetic study of 40 patients⁵⁵ A smaller pharmacokinetic study of 40 patients
Chou et al. Impact of the CYP2D6 SNP on the concentration of and therapeutic response to donepezil. J Formos Med Assoc, 2022 found that G/G homozygotes had higher donepezil plasma concentrations and a 90% response rate compared to 50% in C allele carriers (OR 9.0, p=0.015). The underlying mechanism may involve CYP2D6 channeling donepezil into specific metabolic pathways that modulate the drug's binding properties in the brain, though this remains under investigation.

For broader drug metabolism, Llerena et al.⁶⁶ Llerena et al.
Llerena et al. CYP2D6 -1584C>G promoter polymorphism and debrisoquine ultrarapid hydroxylation. Pharmacogenomics, 2013 showed in 320 healthy volunteers that G allele carriers had lower metabolic ratios (faster debrisoquine hydroxylation), confirming the G allele's promoter-boosting effect on CYP2D6 enzyme activity.

Not all studies show consistent effects — three negative replication studies in specific ethnic populations (Polish, Chinese) found no significant association with donepezil response, highlighting that rs1080985's impact interacts with background haplotypes and is not identical across all ancestries.

Practical Actions

For individuals carrying one or two C alleles, the most clinically actionable implication is awareness before initiating donepezil therapy in Alzheimer's disease, or when prescribed other CYP2D6-substrate medications. The C allele's reduced promoter efficiency means the liver produces less CYP2D6 enzyme, shifting drug metabolism toward slower clearance — which affects both prodrug activation (e.g. codeine, which needs CYP2D6 to form morphine) and active drug elimination (e.g. tricyclic antidepressants that accumulate when metabolism is slow).

The practical implication depends on the specific drug: for prodrugs (codeine, tramadol), lower CYP2D6 activity means less activation and potentially reduced efficacy. For active drugs metabolized by CYP2D6 (many antidepressants, tamoxifen, beta-blockers), lower activity means slower elimination and higher plasma concentrations. Discuss your rs1080985 genotype with prescribing clinicians before starting any high-priority CYP2D6-substrate drug.

Interactions

rs1080985 interacts closely with rs16947 (CYP2D6*2 coding variant, Arg296Cys) and the enhancer SNP rs5758550. The net CYP2D6 activity score for any individual depends on the complete haplotype structure across these loci, not on any single SNP in isolation. In particular, a CYP2D6*2A haplotype that carries the -1584C allele plus the rs16947 A allele and lacks the rs5758550 enhancer variant may show meaningfully reduced activity compared to an individual with the G allele at all three positions. Full clinical CYP2D6 pharmacogenomic testing that resolves star alleles and copy number provides a more complete picture than individual SNPs. The CYP2D6*4 variant (rs3892097) is the major no-function allele in Europeans; if combined with a reduced-function *2A haplotype, the diplotype may approach intermediate metabolizer territory.

LRP1 (Low-Density Lipoprotein Receptor-Related Protein 1) is a giant scavenger receptor expressed throughout the body but especially important in the brain and blood vessels. While it was originally studied for its role in clearing lipoproteins from the bloodstream, research over the past decade has revealed that LRP1 serves as a critical hub connecting leptin signaling¹¹ leptin signaling
Leptin is the "satiety hormone" produced by fat cells that tells the brain to stop eating, glutamate neurotransmission, and vascular integrity. The variant rs11172113 sits in an intronic enhancer region that controls how much LRP1 protein your cells produce.

The Mechanism

The rs11172113 variant lies within intron 1 of LRP1 at an enhancer element that regulates gene expression. The C allele creates a binding site for the transcriptional repressor SNAIL²² transcriptional repressor SNAIL
SNAIL (encoded by SNAI1) is a zinc-finger transcription factor that silences gene expression by binding to E-box-like sequences, which reduces LRP1 expression. The T allele does not bind SNAIL and allows higher LRP1 expression. This allele-specific repression has been confirmed experimentally: SNAIL knockdown in C/C cells significantly increases LRP1 levels, while having no effect in T/T cells.

In the brain, LRP1 directly binds both leptin and the leptin receptor complex³³ leptin and the leptin receptor complex
Liu et al. showed LRP1 is required for leptin receptor phosphorylation and downstream STAT3 activation in hypothalamic neurons. When neuronal LRP1 is deleted in mice, the result is obesity driven by increased food intake and decreased energy expenditure — essentially a state of leptin resistance. Even targeted deletion of LRP1 specifically in the hypothalamus is sufficient to trigger accelerated weight gain. This means reduced LRP1 expression (as occurs with the C allele) could weaken the brain's ability to respond to leptin's satiety signals.

LRP1 also modulates NMDA receptor trafficking⁴⁴ NMDA receptor trafficking
LRP1 controls the surface distribution and internalization of the NR2B subunit of NMDA receptors, which are glutamate-gated ion channels involved in pain and migraine pathophysiology at neuronal synapses. Altered NMDA receptor dynamics may contribute to the cortical hyperexcitability that underlies migraine with aura.

The Evidence

The original GWAS discovery⁵⁵ original GWAS discovery
Chasman et al. Genome-wide Association Study Reveals Three Susceptibility Loci for Common Migraine in the General Population. Nat Genet, 2011 identified rs11172113 in a study of 5,122 migraineurs and 18,108 controls, with the T allele conferring protection (OR 0.90, 95% CI 0.87-0.93, p = 4.3 x 10^-9). This has been replicated in European populations⁶⁶ replicated in European populations
Esserlind et al. Replication and meta-analysis of common variants identifies a genome-wide significant locus in migraine. Eur J Neurol, 2013 and confirmed as the likely causal variant with posterior probability of 1.0 in fine-mapping studies.

A landmark mouse study⁷⁷ mouse study
Liu et al. Lipoprotein receptor LRP1 regulates leptin signaling and energy homeostasis in the adult central nervous system. PLoS Biol, 2011 demonstrated that brain-specific LRP1 knockout mice develop obesity with impaired hypothalamic STAT3 phosphorylation — the key downstream step in leptin signaling. LRP1 overexpression rescued the phenotype, confirming a causal role.

The most mechanistically detailed work came from Liu et al. 2024⁸⁸ Liu et al. 2024
LRP1 Repression by SNAIL Results in ECM Remodeling in Genetic Risk for Vascular Diseases. Circ Res, 2024, which used iPSC-derived smooth muscle cells to show that rs11172113 lies in an enhancer region and that the C allele permits SNAIL-mediated repression of LRP1. This variant was identified as the causal SNP for multiple traits: migraine, fibromuscular dysplasia (OR 1.34, p = 2 x 10^-10), pulse pressure, and spontaneous coronary artery dissection.

The effect size for migraine is modest (OR ~1.11 per C allele), consistent with a common regulatory variant contributing to a complex trait. However, the convergence of GWAS evidence across multiple vascular and neurological phenotypes, combined with robust functional characterization, places this variant on solid mechanistic ground.

Practical Actions

Carriers of the C allele have reduced LRP1 expression, which may impair central leptin signaling and increase migraine susceptibility. Magnesium supplementation is well-established for migraine prevention and also modulates NMDA receptor activity — the same pathway influenced by LRP1. Riboflavin (vitamin B2) at 400 mg/day has strong evidence for migraine prophylaxis and supports mitochondrial function in neurons. Omega-3 fatty acids reduce neuroinflammation and support vascular health, both relevant to the LRP1 pathway.

Monitoring leptin and metabolic markers is relevant because reduced LRP1 function in the brain may contribute to leptin resistance even in the absence of frank obesity.

Interactions

LRP1's role in leptin signaling creates a potential interaction with variants in the leptin pathway. Carriers of rs11172113 CC who also carry the LEPR rs1137101 GG (reduced leptin receptor sensitivity) may experience compounded leptin resistance — both the receptor itself and the LRP1 co-receptor that facilitates its signaling are impaired. Similarly, LEP rs7799039 AA carriers (elevated leptin production) with reduced LRP1 expression may develop more pronounced leptin resistance, as the brain's capacity to transduce the leptin signal is diminished despite high circulating levels.

The FTO rs9939609 AA genotype (increased appetite drive) may compound with reduced LRP1 signaling to further weaken central satiety regulation, though this interaction has not been directly studied in humans.

Regulatory T cells (Tregs) are the immune system's peacekeepers — specialized lymphocytes that suppress excessive inflammation and prevent the immune system from attacking harmless environmental antigens like pollen, pet dander, and food proteins. Their ability to do this depends critically on a surface protein called GARP¹¹ GARP
Glycoprotein A repetitions predominant, encoded by LRRC32. Also known as LRRC32 (leucine-rich repeat containing 32). Essential for anchoring latent TGF-beta on the Treg surface.. rs11236797 sits in a distal enhancer element at chromosome 11q13.5 that controls how much GARP your Tregs produce. The risk A allele weakens this enhancer, reducing GARP expression and — with it — the Tregs' capacity to enforce immune tolerance. The result is a measurably elevated risk for asthma, allergic rhinitis, hay fever, and inflammatory bowel disease.

The Mechanism

GARP functions as a docking receptor for latent TGF-beta²² docking receptor for latent TGF-beta
TGF-beta (transforming growth factor-beta) is produced as an inactive precursor tethered to a latency-associated peptide. GARP anchors this complex to the Treg surface, positioning it for integrin-mediated activation on the Treg cell surface. When a Treg contacts an inflammatory cell through T-cell receptor signaling, surface-bound latent TGF-beta is converted to its active form. The activated TGF-beta then suppresses nearby effector T cells, mast cells, and innate immune cells — the cellular machinery behind allergic reactions.

The enhancer containing rs11236797 recruits the transcription factors STAT5 and NF-κB³³ STAT5 and NF-κB
STAT5 responds to cytokines like IL-2 that maintain Treg survival; NF-κB responds to inflammatory signals, ensuring GARP is upregulated precisely when Tregs need to suppress inflammation to drive signal-dependent LRRC32 expression. In a landmark 2020 Nature study, Nasrallah et al. demonstrated that human Treg cells carrying risk variants at this locus show reduced histone acetylation at the enhancer and lower GARP protein levels⁴⁴ reduced histone acetylation at the enhancer and lower GARP protein levels
Histone acetylation is an epigenetic mark of active gene regulatory elements; reduced acetylation means the enhancer is less open and less active, producing less GARP mRNA and protein. Mice lacking this enhancer entirely retained viable Treg cells — but those Tregs were unable to control colitis in cell-transfer experiments, directly demonstrating that GARP loss impairs Treg suppressor function without abolishing Tregs altogether.

The Evidence

The earliest genome-wide signal at this locus came from a 2011 Lancet study by Ferreira et al.⁵⁵ 2011 Lancet study by Ferreira et al.
Australian Asthma Genetics Consortium, n=57,800 combined that identified 11q13.5/LRRC32 as a genome-wide significant asthma risk locus (OR=1.09, p=1.8×10⁻⁸), with a stronger signal specifically for atopic asthma (OR=1.33, p=7×10⁻⁴), consistent with the locus acting through allergic sensitization rather than non-allergic airway inflammation.

The most comprehensive allergic disease GWAS to date — Ferreira et al. Nature Genetics 2017⁶⁶ Ferreira et al. Nature Genetics 2017
360,838 participants from the UK Biobank and international cohorts; 136 independent genome-wide significant loci identified, 73 novel; 23andMe was a major contributing cohort — confirmed 11q13.5 as a shared risk locus for asthma, hay fever, and eczema. Most loci in this study act across all three allergic conditions, underscoring that allergic diseases share more genetic architecture than they differ. The GWAS Catalog records associations for rs11236797-A with asthma (OR=1.12, p=6×10⁻⁶²), childhood-onset asthma (OR=1.16, p=1×10⁻⁹³), allergic rhinitis (OR=1.14, p=5×10⁻³²), and inflammatory bowel disease (p=7×10⁻³³), pointing to a broadly immunoregulatory locus rather than a disease-specific one.

The mechanistic connection was established by the 2020 Nature study from Nasrallah et al.⁷⁷ 2020 Nature study from Nasrallah et al.
Functional genomics in 91–123 healthy human donors, confirmed in mouse enhancer-knockout models, which linked the disease-associated enhancer directly to GARP expression in human Tregs — bridging the gap between GWAS signals and biological function.

Practical Actions

The A allele does not eliminate Treg function — it reduces the efficiency of Treg-mediated suppression at the molecular level. This has practical implications: environmental and lifestyle factors that support Treg numbers and function can partially offset the genetic disadvantage.

Vitamin D promotes Treg differentiation and upregulates FoxP3, the master transcription factor for Treg identity. Maintaining adequate 25-OH vitamin D levels supports Treg biology across multiple pathways. Probiotic strains with documented Treg-supporting effects — particularly those used in clinical allergy studies — can expand peripheral Treg populations in the gut mucosa, the site where much of the GARP-dependent immunosuppression operates.

The IBD association for this locus is particularly meaningful: the gut is where GARP-expressing Tregs are most densely deployed for mucosal tolerance. Carriers of the A allele who have digestive symptoms, a family history of IBD, or significant allergic disease burden may benefit from earlier gastroenterological evaluation.

Interactions

rs11236797 (LRRC32/GARP) and rs17293632 (SMAD3) operate in the same TGF-beta tolerance pathway: GARP activates latent TGF-beta on the Treg surface; SMAD3 is the primary intracellular signal transducer that TGF-beta activates inside target cells. Reduced GARP (this variant) and reduced SMAD3 (rs17293632) would compound to impair TGF-beta signal initiation and propagation. Both variants are associated with asthma and IBD. The combination of risk alleles at both loci represents a compound Treg tolerance defect that no single-genotype analysis captures.

The immune system uses the interleukin-23 (IL-23) signalling axis as a master switch for chronic inflammatory responses. When IL-23 binds the IL-23 receptor (IL23R) on Th17 cells¹¹ Th17 cells
T helper 17 cells — a subset of CD4+ T cells that produce interleukin-17A and interleukin-17F; their chronic activation drives gut and joint inflammation in IBD, ankylosing spondylitis, and psoriasis, it sustains production of the pro-inflammatory cytokine IL-17A and drives the kind of persistent mucosal inflammation that characterises Crohn's disease (CD) and ulcerative colitis (UC). The rs11465770 variant in IL23R sits within an intronic region on chromosome 1, position 67,168,280 (GRCh38)²² chromosome 1, position 67,168,280 (GRCh38)
The IL23R gene spans chromosome 1p31.3; rs11465770 is located in the 5' portion of the gene near the block containing the well-studied protective variants. The T allele at this position marks a protective haplotype that co-segregates with other IL23R variants known to dampen receptor signalling, and it is this haplotype that reduces susceptibility to inflammatory bowel disease.

The Mechanism

rs11465770 is itself an intronic substitution (C>T on the plus strand) with no direct coding consequence. Its biological relevance lies in linkage disequilibrium (LD)³³ linkage disequilibrium (LD)
Two variants are in LD when they are inherited together on the same chromosomal stretch far more often than chance would predict; knowing one allele reliably predicts the other within a population with functional IL23R variants in the same haplotype block. The IL23R gene contains two major LD blocks: a centromeric block (containing exons 5–11) harbouring the most significant GWAS signals, and a 5' region block. The T allele at rs11465770 tags a haplotype architecture where reduced receptor activity is the net effect.

The biological mechanism of protection is best understood through the functional IL23R variants this haplotype accompanies. The key coding variant rs11209026 (R381Q) disrupts the cytoplasmic signalling domain of IL23R. When IL-23 binds the receptor, it normally phosphorylates the transcription factor STAT3, driving IL-17A transcription in Th17 cells. T cells carrying the protective haplotype produce only 5.5 pg/mL IL-17A after IL-23 stimulation, compared with 36.0 pg/mL in non-carriers⁴⁴ 36.0 pg/mL in non-carriers
6.5-fold reduction measured in Th17 effector cells; Th17 differentiation itself is unaffected, meaning pathogen defence is preserved while the chronic inflammatory overdrive is attenuated. More broadly, protective IL23R variants act through impaired protein stability and intracellular trafficking⁵⁵ impaired protein stability and intracellular trafficking
The R381Q, G149R, and V362I protective variants all display reduced cell-surface receptor expression due to ER retention or accelerated degradation, limiting the number of receptors available to respond to IL-23 stimulation. rs11465770 T, as a haplotype tag, identifies individuals who carry this reduced-receptor architecture on at least one chromosomal copy.

This mechanism is clinically validated by the pharmaceutical success of anti-IL-23 biologics. Drugs such as risankizumab, guselkumab, and mirikizumab — all targeting the IL-23 p19 subunit that binds IL23R — are now first-line or second-line treatments for moderate-to-severe Crohn's disease and ulcerative colitis, achieving remission rates of 40–60% in clinical trials. Carriers of protective IL23R haplotypes are, in effect, born with a partial pharmacological blockade of the same pathway these drugs target.

The Evidence

The IL23R locus was identified as an IBD susceptibility region in the landmark 2006 GWAS published in Science⁶⁶ published in Science
Duerr et al. — genome-wide discovery in 547 CD cases and 548 controls, with replication in Jewish and non-Jewish cohorts; the IL23R centromeric haplotype block showed P values as low as 10⁻¹³ for Crohn's disease. Multiple intronic IL23R variants in the same haplotype block as rs11465770 achieved protective odds ratios of 0.30–0.70 for CD, with protective allele frequencies in European controls reaching 10% or higher.

A meta-analysis of 60 CD case-control studies⁷⁷ meta-analysis of 60 CD case-control studies
Xu et al. Scientific Reports 2015; 22,820 CD cases and 27,401 controls across multiple populations confirmed that IL23R intronic variants — particularly rs7517847, a member of the same haplotype architecture — reduce CD risk with OR = 0.70 (95% CI 0.66–0.74, P<0.001) in Caucasian populations. The protective signal is not present in East Asian populations, consistent with the near-absence of the T allele at rs11465770 in East Asian gnomAD data (allele frequency <0.04%).

For ulcerative colitis, a meta-analysis of 33 studies⁸⁸ meta-analysis of 33 studies
Zhong et al. Oncotarget 2016; 10,527 UC cases, 15,142 controls; OR=0.76 (95% CI 0.64–0.90, P=0.002) for protective IL23R alleles vs risk alleles in Caucasians confirmed a consistent 24% reduction in UC risk among carriers of the protective haplotype. The effect size is smaller for UC than for CD, mirroring the pattern seen across all IL23R variants in the gene — the locus has a stronger and more consistently replicated association with CD than with UC across populations.

The T allele at rs11465770 is found at approximately 10% allele frequency in Europeans and South Asians (~8%), is exceedingly rare in East Asians (<0.1%), and uncommon in Africans (~1.4%). This population structure — high frequency in populations with high IBD prevalence, very low frequency in East Asian populations where IBD incidence is rising but was historically low — is consistent with a variant under modest balancing selection in environments where the IL-23/Th17 pathway plays a critical host-defence role.

Practical Actions

Carrying the T allele at rs11465770 represents a genuinely favourable finding with respect to gut inflammatory disease risk. Heterozygous CT carriers have partial haplotype protection on one chromosomal copy; TT homozygotes carry the maximum protection available at this locus. This does not eliminate IBD risk — other genetic contributors (NOD2, ATG16L1, CARD9) and environmental factors are also important — but it meaningfully shifts the baseline.

For individuals who develop IBD despite carrying the protective T allele, or who have a family history of IBD alongside this protective genotype, the IL-23/Th17 biology is still directly relevant to treatment selection. Anti-IL-23 biologics (risankizumab, guselkumab, ustekinumab for the IL-12/23 p40 subunit) target the exact pathway in which this haplotype confers its protection, and carriers of protective IL23R haplotypes may show differential pharmacological responses — though this remains an area of active research rather than established clinical guidance.

Interactions

rs11465770 is in linkage disequilibrium with the principal IL23R protective variants: rs11209026 (R381Q, the functional coding change), rs11465804, rs7517847, rs10489629, and rs1343151. Individuals carrying the T allele at rs11465770 are likely to also carry protective alleles at one or more of these correlated positions, though LD is not perfect and each provides independent information. Combining haplotype data across all typed IL23R markers gives a more complete picture of IL23R protective coverage than any single variant.

Beyond IL23R itself, the IL-23/Th17 pathway involves upstream (IL12B, STAT3, JAK2) and downstream (IL17A, IL17RA, IL22) genes whose variants interact epistatically with IL23R haplotypes. Studies of the IL23/IL17 pathway in CD have documented significant gene-gene interactions between IL23R risk haplotypes and IL17A, IL17RA, and IL12RB2 haplotypes: combined carrier status of multiple pathway risk alleles pushes CD odds ratios to 4.3 or higher, while protective IL23R haplotypes partially buffer against these epistatic risks.

Narcolepsy type 1 — the form marked by sudden muscle weakness (cataplexy) and daytime sleep attacks — is not simply a sleep disorder. It is an autoimmune disease where the immune system destroys a small cluster of neurons in the hypothalamus that produce hypocretin¹¹ hypocretin
also called orexin; a neuropeptide that stabilizes the switch between waking and sleep states. Once lost, these roughly 70,000 neurons do not regenerate, and the wakefulness-sleep boundary becomes unstable for life. The genetic variant rs1154155 sits inside the T-cell receptor alpha (TRA) locus — the gene region that encodes the alpha chain of the T-cell receptor, the molecule T cells use to recognize peptides presented by HLA class II molecules.

The Mechanism

The TRA locus spans a large stretch of chromosome 14q11.2 and contains approximately 70 J-segment genes that are somatically recombined during T-cell development to generate diverse receptor specificities. rs1154155 falls within an 18-kb region encompassing J-segment genes TRAJ24 and TRAJ28. It is in near-perfect linkage disequilibrium with a functional coding variant (rs1483979) that encodes a leucine-to-phenylalanine substitution in the complementarity-determining region 3 (CDR3) of the TRAJ24 segment — the very region that contacts the peptide-HLA complex directly.

The leading biological model is that specific TRA chain sequences encoded near the G-risk haplotype generate T-cell receptors with higher affinity for hypothalamic peptides presented by HLA-DQB1*06:02²² HLA-DQB1*06:02
the HLA allele present in >95% of narcolepsy type 1 patients. When the immune system encounters a cross-reactive antigen — most likely a peptide from influenza H1N1 hemagglutinin or related pathogens — these high-affinity T cells mount an immune response that simultaneously targets hypocretin neurons. The 2009–2010 H1N1 pandemic and the AS03-adjuvanted Pandemrix vaccine both dramatically elevated narcolepsy incidence in children carrying HLA-DQB1*06:02 and TRA risk genotypes, providing the strongest epidemiological confirmation of this autoimmune trigger model.

The Evidence

The discovery study by Hallmayer et al.³³ Hallmayer et al.
Narcolepsy is strongly associated with the T-cell receptor alpha locus. Nature Genetics, 2009 genotyped 1,830 narcolepsy cases and 2,164 controls across European, Asian, and African-American populations. rs1154155 showed the strongest signal (p < 10⁻²¹) with an average allelic odds ratio of 1.69. Genotypic odds ratios were 1.94 for heterozygotes (GT) and 2.55 for homozygotes (GG) relative to TT — a classic additive dose-response pattern. This was the first documented association between the TRA locus and any human disease.

Ollila et al.⁴⁴ Ollila et al.
Narcolepsy risk loci outline role of T cell autoimmunity and infectious triggers in narcolepsy. Nature Communications, 2023 expanded this to 6,073 cases and 84,856 controls, confirming the TRA signal and mechanistically linking rs1154155 to altered TRAJ24 and TRAJ28 chain usage in T-cell repertoire sequencing data (posterior probability 0.958 for the TRAJ28 co-localization signal).

A Chinese case-control study (Ouyang et al. 2020⁵⁵ Ouyang et al. 2020, 903 cases, 1,981 controls) found that among HLA-DQB1*06:02-positive individuals, carrying GG or GT genotypes conferred an OR of 9.33 compared to TT (p = .017), illustrating how HLA and TRA risk genotypes multiply each other's effect. In HLA-negative individuals, the TRA association was not statistically significant — confirming that TRA genotype modifies rather than independently causes narcolepsy risk.

Practical Actions

The absolute risk of developing narcolepsy remains low even in those carrying both HLA-DQB1*06:02 and the TRA G-risk genotype (lifetime prevalence is approximately 1 in 2,000). The clinical value of knowing this genotype lies in:

• Earlier recognition: Symptom clusters (excessive daytime sleepiness, sleep paralysis, hypnagogic hallucinations, cataplexy) often precede diagnosis by 5–10 years. Carriers with emerging symptoms should seek polysomnography and hypocretin CSF measurement promptly rather than attributing symptoms to lifestyle.
• Infection and vaccination vigilance: G-risk carriers — especially those known to be HLA-DQB1*06:02 positive — should discuss vaccine timing and neurological symptoms with their physician when novel H1N1-type antigens are in circulation.
• Sleep architecture monitoring: Narcolepsy type 1 disrupts REM sleep regulation distinctively. Sleep studies (overnight polysomnography + MSLT) are the diagnostic standard and can detect REM-onset abnormalities before full symptom development.

Interactions

The TRA locus does not act alone. Narcolepsy type 1 requires HLA-DQB1*06:02 as the essential prerequisite — present in >95% of cases. rs1154155 TRA risk genotypes interact multiplicatively with HLA-DQB1*06:02 status: the OR for narcolepsy among HLA-positive individuals carrying GG or GT genotypes is dramatically higher (OR ~9.33) than in the overall population. Separately, the TCR beta locus variant rs9648789 (TRB locus) also associates with narcolepsy and may compound with TRA risk genotypes, since both alpha and beta chains must heterodimerize to form the functional TCR that recognizes peptide-HLA complexes. The precise combined effect of rs1154155 + rs9648789 has not been quantified in a single study, but pathway logic suggests additive or super-additive risk among carriers of both loci.

The FTO gene — officially the fat mass and obesity-associated gene — contains one of the most replicated loci in obesity genetics. Most research has focused on a cluster of intron-1 variants¹¹ intron-1 variants
The first non-coding interval of FTO, approximately 53,738,000–53,800,000 bp on chromosome 16 (GRCh38), contains the best-studied GWAS hits for obesity including rs9939609 and the causal variant rs1421085 led by rs9939609 and the mechanistic variant rs1421085. But the FTO gene spans over 400 kilobases across chromosome 16, harboring multiple intronic regions, and rs11644943 sits at a distinct chromosomal position²² distinct chromosomal position
chr16:53,961,672 (GRCh38) — approximately 200 kilobases downstream of the primary intron-1 obesity cluster, in a separate linkage disequilibrium block from the well-characterized FTO obesity haplotype.

The Mechanism

rs11644943 is an intronic FTO variant. Like all FTO intronic variants, it does not alter the amino acid sequence of the FTO protein — instead, it sits within a non-coding regulatory region where it may influence local chromatin accessibility, splicing efficiency, or transcriptional output of neighboring exons. The FTO protein itself functions as an N6-methyladenosine (m6A) RNA demethylase³³ N6-methyladenosine (m6A) RNA demethylase
m6A is the most abundant chemical modification on messenger RNA. FTO erases it, altering mRNA stability and translation of target genes including those involved in appetite regulation and adipocyte differentiation. Elevated FTO expression raises ghrelin mRNA stability (increasing hunger drive) and promotes white adipocyte differentiation over thermogenic beige adipocytes (reducing resting energy expenditure). The T allele at rs11644943 may act through a similar upregulatory mechanism on FTO transcript levels in relevant tissues, though the specific regulatory element disrupted by this variant has not been characterized.

Critically, the rs11644943 T allele and nearby rs7206790 are not in linkage disequilibrium with previously reported obesity-associated FTO variants⁴⁴ the rs11644943 T allele and nearby rs7206790 are not in linkage disequilibrium with previously reported obesity-associated FTO variants
Xu et al. 2014 (PMID:25251416) explicitly confirmed that rs11644943 and rs7206790 are independent of rs9939609 and other previously catalogued FTO GWAS hits in the Chinese school-age population studied. This makes rs11644943 a potentially independent FTO obesity signal — tagging genetic variation in a different part of the gene rather than serving as a proxy for the primary intron-1 cluster.

The Evidence

Published evidence for rs11644943 is limited. The primary study is a case-control analysis in Chinese school-age children⁵⁵ case-control analysis in Chinese school-age children
Xu et al. Rs7206790 and rs11644943 in FTO gene are associated with risk of obesity in Chinese school-age population. PLoS One 2014. PMID:25251416 (500 obese cases, 500 matched controls). The investigators identified the A allele as protective: carriers of the AA genotype had a substantially lower risk of obesity compared to T-allele carriers (adjusted OR 0.16, 95% CI 0.04–0.72). The association survived adjustment for age, sex, and residential location. When combined with rs9939609 in a three-SNP genetic risk score model, dose-response relationships with obesity were observed, though rs11644943's contribution was confirmed as independent of the primary FTO locus.

Direct evidence for rs11644943 outside this single study is not yet available in the peer-reviewed literature, and the variant is not catalogued in ClinVar or the GWAS Catalog. The evidence level is therefore emerging — the association requires replication in larger and ethnically diverse cohorts before it can be considered a confirmed independent FTO signal. For context, the broader FTO locus has strong, replicated evidence for influencing obesity risk across populations⁶⁶ strong, replicated evidence for influencing obesity risk across populations
Frayling et al. Science 2007. PMID:17434869, with the primary intron-1 signal explaining an additive ~3 kg difference in body weight per risk allele in Europeans.

Practical Actions

Because rs11644943 is an FTO intronic variant, the genotype-specific advice derives from what is known about the FTO pathway. Physical activity consistently attenuates FTO-related obesity risk: a meta-analysis of 218,166 adults⁷⁷ meta-analysis of 218,166 adults
Kilpeläinen et al. 2011. PMID:22069379 found that physically active individuals showed a 27% lower FTO allele-associated obesity effect than sedentary individuals — one of the largest known gene-environment interactions in common disease genetics. Dietary protein composition also matters: higher-protein meals improve satiety signaling in FTO risk allele carriers by compensating for impaired gut hormone (GLP-1, PYY) responses. These effects have been demonstrated across FTO intron variants generally, though not yet specifically for rs11644943.

Interactions

rs11644943 is not in LD with the primary FTO intron-1 obesity cluster (rs9939609, rs1421085, rs8050136), making it a potentially additive second FTO signal. Users who carry both the T allele at rs11644943 AND the A risk allele at rs9939609 may carry compounded FTO pathway burden through independent intronic regulatory mechanisms. The adjacent variant rs7206790 (chr16:53,763,996) was studied alongside rs11644943 in the Chinese school-age cohort but appears to be in a different intron and LD block; the two variants together may index distinct regulatory regions across the FTO gene body.

Your veins are not passive conduits. They are living structures whose walls depend on a continuous interplay of growth signals, smooth muscle tone, and extracellular matrix remodeling. VEGFA — vascular endothelial growth factor A¹¹ VEGFA — vascular endothelial growth factor A
the master regulator of angiogenesis; controls endothelial proliferation, vessel permeability, and vascular tone through binding to VEGFR1 and VEGFR2 — sits at the center of this balance. The rs11967262 variant lies approximately 7 kb upstream of the VEGFA gene on chromosome 6p21.1, in a region that likely influences VEGFA expression in vascular tissue. Carriers of the G allele show a small but reproducibly elevated risk of developing varicose veins.

Varicose veins affect an estimated 25% of women and 15% of men in Western populations, with heritability estimates ranging from 17–49%. They are not merely cosmetic — dilated, tortuous superficial veins reflect underlying venous hypertension and valve incompetence that, if unchecked, can progress to chronic venous insufficiency, venous ulceration, and superficial thrombophlebitis.

The Mechanism

VEGFA is significantly overexpressed in the walls of varicose veins compared to normal veins. This upregulation is now understood as both a cause and consequence of venous dysfunction. VEGFA is a potent inducer of vascular permeability — it loosens endothelial junctions, allows plasma proteins to leak into the vessel wall, and drives edema. In the venous context, chronically elevated VEGFA signaling through VEGFR2 promotes endothelial cell proliferation and smooth muscle cell phenotypic switching from contractile to synthetic state, weakening the structural integrity of the vein wall.

VEGFA also induces vasodilation through nitric oxide–dependent pathways, reducing venous tone. In the deep venous system this vasodilation is tightly regulated; but in the superficial venous system — which lacks surrounding muscle support — sustained VEGFA overactivity can tip the balance toward stasis, reduced return flow, and progressive dilation. Hypoxia and venous hypertension, which develop as veins dilate and valves fail, further upregulate VEGFA expression in a feed-forward loop.

The rs11967262 variant is intergenic at 6p21.1, ~7 kb upstream of the VEGFA transcription start site²² intergenic at 6p21.1, ~7 kb upstream of the VEGFA transcription start site
the exact regulatory element affected has not been functionally characterized; the association may act through eQTL effects on VEGFA expression in venous endothelial or smooth muscle cells, or through linkage disequilibrium with a causal variant in the VEGFA promoter region. The VEGFA gene has several known regulatory polymorphisms in its 5′ UTR (notably rs2010963) that alter expression levels, and rs11967262 may tag a similar effect.

The Evidence

The primary association evidence comes from the largest varicose vein GWAS published to date. Ahmed et al. 2022 in Nature Communications³³ Ahmed et al. 2022 in Nature Communications
two-stage analysis: 22,473 cases and 379,183 controls from UK Biobank; replication in 113,041 cases and 295,928 controls from 23andMe identified 49 genome-wide significant signals across 46 loci. The rs11967262 G allele reached a meta-analysis p-value of 1.45×10⁻¹⁹ with an odds ratio of 1.09 — a modest per-allele effect that is nonetheless very well-powered given the exceptional sample size.

The biological plausibility of the VEGFA locus is strong. The varicose vein wall consistently shows elevated VEGFR1, VEGFR2, and VEGFR3 expression compared to normal venous tissue, particularly in cases complicated by thrombophlebitis. Plasma VEGFA levels are significantly higher in varicose vein patients. A Mendelian randomization study⁴⁴ Mendelian randomization study
genetically instrumented VEGF levels examined as causal instrument in ~16,000 Europeans found that higher circulating VEGF levels associate with increased venous thromboembolism risk (OR 1.064, 95% CI 1.009–1.122), extending the causal inference beyond simple association.

A separate candidate gene study found the VEGFA promoter variant rs2010963 C allele to be protective against varicose veins⁵⁵ protective against varicose veins
case-control study of 448 patients and 609 controls in ethnic Russians; C allele OR 0.73, 95% CI 0.59–0.91, p=0.004, directly demonstrating that VEGFA regulatory variants influence varicose vein susceptibility in an allele-specific manner.

The OR of 1.09 per G allele is modest by clinical standards. The absolute risk increase for GG homozygotes compared to CC homozygotes is approximately 18–19% relative risk elevation — meaningful in context of a common, heritable condition but not deterministic. Most GG carriers will not develop clinically significant varicose veins if they manage modifiable risk factors aggressively.

Practical Actions

The VEGFA locus variant acts through venous wall biology — specifically vascular permeability, smooth muscle tone, and extracellular matrix remodeling. Interventions that reduce venous hydrostatic pressure, improve wall tone, and reduce conditions that upregulate VEGFA (hypoxia, venous stasis, inflammation) are directly relevant to this genotype.

Graduated compression garments (20–30 mmHg for prevention; 30–40 mmHg for established insufficiency) mechanically reduce the venous hypertension that drives VEGFA upregulation. Occupations requiring prolonged standing or sitting — which are the strongest modifiable risk factors for varicose veins independent of genetics — are particularly relevant for G allele carriers to mitigate.

Micronized purified flavonoid fractions (MPFF, e.g., diosmin-hesperidin combinations) have evidence from randomized trials for reducing venous edema and improving microcirculation in chronic venous disease. They act partly by reducing vascular permeability — a VEGFA-mediated mechanism — making them particularly relevant for VEGFA-pathway genotypes.

Interactions

The VEGFA locus interacts with other varicose vein risk loci identified in the same GWAS — including CASZ1 (vascular development), PIEZO1 (mechanosensitive ion channel), and extracellular matrix genes (ELN, LTBP3, COL3A1). Individuals carrying multiple risk alleles across these independent loci will have additive polygenic risk that compounds the modest VEGFA effect. The polygenic risk score derived in the Ahmed et al. 2022 cohort has predictive utility for varicose vein surgery, suggesting that cumulative genetic burden at VEGFA and these other loci has clinical prognostic value.

VEGFA also interacts with the coagulation cascade — elevated VEGFA increases vascular permeability that can promote local thrombus formation in venous stasis conditions. Carriers of thrombophilia variants (Factor V Leiden rs6025, Prothrombin rs1799963, MTHFR rs1801133) who also carry GG at rs11967262 may have compounded venous disease risk through convergent mechanisms.