The TFR2 gene encodes transferrin receptor 2, a liver-expressed iron-sensing protein that plays a central role in calibrating hepcidin¹¹ hepcidin
A peptide hormone produced by the liver that acts as the master regulator of systemic iron availability — it controls how much iron the gut absorbs and how much is released from recycling macrophages. When diferric transferrin (iron-loaded transferrin in the blood) binds TFR2, it triggers a signaling cascade that increases hepcidin secretion, thereby restricting iron entry and preventing overload. rs4434553 sits approximately 2 kb upstream of the TFR2 coding sequence in a regulatory region and appears to modulate how strongly TFR2 expression responds to iron signals.

The Mechanism

As an upstream regulatory variant, rs4434553 (c.-258+123T>C in paper notation; A>G on the GRCh38 plus strand, TFR2 being minus-strand) likely influences TFR2 transcriptional activity or mRNA processing of non-coding transcripts overlapping this region. The functional consequence is inferred from the phenotypic data: carriers of the G allele (coding-strand C) show higher circulating hepcidin and lower serum ferritin²² higher circulating hepcidin and lower serum ferritin
Consistent with increased TFR2-mediated hepcidin upregulation — the G/C allele may increase TFR2 expression or shift the iron-sensing threshold compared to AA homozygotes. The A allele (coding-strand T) appears to blunt TFR2-mediated hepcidin signaling, allowing relatively greater iron accumulation under iron-replete conditions.

Because TFR2 is also expressed in the retinal pigment epithelium³³ retinal pigment epithelium
The RPE is the metabolic support layer beneath the photoreceptors; it manages the continuous renewal of the outer photoreceptor segments and is exquisitely sensitive to iron-driven oxidative damage, the same regulatory variant may influence local iron control in the eye — a separate tissue expression mechanism from its hepatic hepcidin role.

The Evidence

The clearest finding for rs4434553 comes from a Chinese Han case-control study of non-alcoholic fatty liver disease. Pan et al. 2022⁴⁴ Pan et al. 2022
Gastroenterol Rep; case-control study at two affiliated hospitals of Fujian Medical University; TFR2 promoter variants genotyped in NAFLD cases and controls; rs4434553 GA/GG genotype associated with reduced NAFLD risk vs AA (OR=0.630, 95% CI 0.504-0.788); GA/GG correlated with lower serum ferritin and higher serum hepcidin. This direction is biologically coherent: lower hepcidin in AA carriers allows more iron absorption and release from recycling macrophages, increasing hepatic iron deposition — a recognized driver of NAFLD progression through oxidative stress and mitochondrial dysfunction⁵⁵ NAFLD progression through oxidative stress and mitochondrial dysfunction
Iron catalyzes the Fenton reaction in hepatocytes, generating reactive oxygen species that promote hepatocyte injury, fibrosis, and necroinflammation.

For age-related macular degeneration, the picture is more nuanced. Wysokinski et al. 2014⁶⁶ Wysokinski et al. 2014
Dis Markers; 493 AMD cases and 171 controls; rs4434553 showed no overall AMD association but demonstrated subgroup effects — the C allele (plus-strand G) increased AMD occurrence in subjects over age 72, and the CC genotype was enriched in AMD patients with BMI below 26. These subgroup-specific signals are consistent with a modifier role: in lean, older individuals, the retinal iron regulatory environment may differ enough that TFR2 expression level becomes a relevant AMD determinant.

A pooled Parkinson's disease analysis found a haplotype containing the A allele at rs4434553 combined with the G allele at neighboring rs10247962 showed a suggestive protective association with Parkinson's disease (OR=0.87, 95% CI 0.74-1.02)⁷⁷ suggestive protective association with Parkinson's disease (OR=0.87, 95% CI 0.74-1.02)
The confidence interval crosses 1.0, making this finding non-significant at conventional thresholds — it is hypothesis-generating only. Iron dysregulation in the substantia nigra is a recognized feature of Parkinson's pathology, and TFR2 haplotype effects on brain iron regulation remain plausible but unconfirmed.

Practical Actions

For AA homozygotes — who have the least favorable iron-regulatory profile from this variant — the most actionable response is periodic monitoring of serum ferritin and iron saturation, particularly if they carry other iron-loading variants (HFE, TMPRSS6). The AMD signal specifically for the G allele in age-stratified subgroups means that GG homozygotes approaching age 70+ may benefit from proactive retinal screening. The overall evidence base remains emerging: none of these associations have been replicated in large independent cohorts or reached genome-wide significance, so this variant should be interpreted as a secondary iron-pathway signal rather than a primary disease determinant.

Interactions

rs4434553 shares the TFR2 locus with rs2075674, a synonymous coding variant that has been studied independently for AMD risk. Both variants were used together in the Wysokinski 2014 AMD study, but they appear to tag different functional effects within TFR2: rs4434553 sits upstream and likely modulates transcription, while rs2075674 is an exon 16 synonymous variant with potential splice-regulatory activity. The upstream intronic variant rs7385804 in TFR2 is the locus tag with the strongest evidence in iron-parameter GWAS. Combined carrier status across multiple TFR2 variants may capture more of the locus's functional variation than any single SNP alone. In the broader iron-sensing pathway, HFE C282Y (rs1800562) and TMPRSS6 Ala736Val (rs855791) are the highest-impact iron-loading variants; their effects compound with TFR2 regulatory variation through parallel but distinct inputs on hepcidin production.

Your immune system maintains a constant conversation with the microbial world — sensing bacteria, fungi, and other pathogens through a set of molecular detectors called pattern recognition receptors¹¹ pattern recognition receptors
PRRs — the first-responder surveillance system of innate immunity, active before the adaptive immune response can mount a targeted attack. Toll-Like Receptor 2 (TLR2) is one of the most promiscuous of these sensors, recognizing an unusually broad array of microbial signals including bacterial lipoproteins, peptidoglycan from gram-positive cell walls, lipoteichoic acid, mycobacterial components, and fungal zymosan. The T-16934A variant (rs4696480) sits in the TLR2 promoter region — not in the coding sequence, but in the regulatory stretch of DNA that governs how much TLR2 the body produces. Carriers of the A allele appear to have altered TLR2 transcriptional output, and the consequences of this expression shift ripple through several inflammatory conditions.

The A allele at rs4696480 is remarkably common — approximately 50% frequency in Europeans, 58% in East Asians, and 38% in Africans. This means the AA genotype is present in roughly 20% of the global population, making it one of the most prevalent functional variants in the entire TLR family.

The Mechanism

At position -16934 in the TLR2 promoter, the T-to-A substitution alters the local DNA sequence recognized by transcription factors that bind this regulatory region. This is an intronic/upstream regulatory variant²² intronic/upstream regulatory variant
classified as an intron variant on dbSNP, located within TLR2 genomic sequence at chr4:153,685,974 on GRCh38 that likely affects transcription factor binding affinity and chromatin accessibility at the TLR2 locus.

Unlike the nearby missense variant R753Q (rs5743708), which produces a structurally impaired TLR2 protein, T-16934A operates at the transcriptional level — modulating how much functional TLR2 protein is synthesized. The AA genotype is associated with increased TLR2-mediated inflammatory responses in the context of skin and respiratory disease, consistent with upregulated receptor expression. A neonatal cord blood study found that AA genotype carriers born to atopic mothers show significantly increased FOXP3, GITR, and LAG3 expression in regulatory T cells³³ significantly increased FOXP3, GITR, and LAG3 expression in regulatory T cells
along with elevated Th2 cytokines and TNF-α secretion when stimulated — suggesting altered immune programming from birth in genetically susceptible individuals. This early-life immune skewing may set the trajectory for later atopic disease.

The Evidence

Atopic dermatitis is the condition most extensively studied in relation to rs4696480. A Ukrainian pediatric cohort of 103 atopic dermatitis patients and 84 healthy controls found the AA genotype associated with severe AD phenotype⁴⁴ AA genotype associated with severe AD phenotype
OR 6.395 (95% CI 1.240–32.991) — a striking effect size from a relatively small study, indicating the AA genotype substantially increases risk of severe rather than mild disease. This severity association, rather than simply disease presence, mirrors what was found for the coding variant R753Q in earlier work. A systematic review and meta-analysis⁵⁵ systematic review and meta-analysis
confirming significant association across allelic, homozygous, heterozygous, and dominant models of inheritance confirmed the overall association between rs4696480 and atopic dermatitis.

Psoriasis shows a similar pattern. A Turkish case-control study of 140 psoriasis patients and 250 controls found the AA genotype associated with elevated psoriasis risk⁶⁶ AA genotype associated with elevated psoriasis risk
adjusted OR 2.41 (95% CI 1.349–4.292), p=0.003 — robust after covariate adjustment. The involvement of TLR2 in psoriasis is biologically coherent: keratinocytes express TLR2, and its signaling contributes to the IL-17/IL-23 axis that drives psoriatic plaques.

Asthma susceptibility is also elevated. A meta-analysis of 13 studies⁷⁷ meta-analysis of 13 studies
finding OR 2.455 (95% CI 1.235–4.88) under the dominant model — meaning AA+AT carriers have roughly 2.5-fold higher asthma risk than TT homozygotes found significant association between rs4696480 and asthma susceptibility, making it the only TLR2 variant among rs4696480, rs5743708, rs3804099, and rs3804100 to show this association. This implicates airway TLR2 activity in asthma pathogenesis — possibly through enhanced responses to microbial triggers that promote type 2 inflammation in susceptible airways.

Skin tumors associated with HPV are also linked to this variant. A Croatian case-control study of 161 keratoacanthoma cases, 152 common wart cases, and 469 controls found TLR2 rs4696480 A allele and AA genotype significantly overrepresented⁸⁸ TLR2 rs4696480 A allele and AA genotype significantly overrepresented
P<0.001 — with stronger association in common warts, which are caused by HPV, suggesting TLR2 expression level affects antiviral innate responses in skin.

Inflammatory bowel disease susceptibility is also modulated by this variant. A Danish cohort of 624 Crohn's disease patients, 411 ulcerative colitis patients, and 795 controls found rs4696480 associated with IBD risk in combined patient analysis⁹⁹ associated with IBD risk in combined patient analysis
alongside other inflammatory pathway polymorphisms in TLR4, TLR9, TNFRSF1A, IL6R, IL10, IL23R, and PTPN22, consistent with TLR2's role in mucosal immune surveillance of gut bacteria.

Epistatic interactions add another layer of complexity. A study examining interactions between FCER1A (the high-affinity IgE receptor alpha chain) and TLR2 found that TT homozygotes carrying the FCER1A rs2252226 minor allele had higher SCORAD than all other combined genotype groups¹⁰¹⁰ TT homozygotes carrying the FCER1A rs2252226 minor allele had higher SCORAD than all other combined genotype groups
demonstrating that TLR2 expression level interacts with IgE receptor genetics to modify atopic dermatitis severity in a non-additive fashion. This epistatic finding shows that the "normal" TT genotype at rs4696480 is not uniformly protective — its effect depends on co-occurring immune receptor variants.

Practical Implications

The T-16934A variant affects the gain, not the function, of TLR2 signaling. Unlike R753Q (which produces a broken receptor), the A allele at rs4696480 may produce more TLR2, creating a system that amplifies inflammatory responses to microbial ligands. This creates a somewhat paradoxical pattern: higher TLR2 expression can drive more intense inflammatory skin responses (atopic dermatitis, psoriasis) while also potentially affecting mucosal immune homeostasis in the gut.

For AA genotype carriers, the actionable implications center on inflammatory skin conditions and atopic disease management. The variant does not impair pathogen recognition — the TLR2 receptors produced are functional, just potentially more numerous — so infection susceptibility differs fundamentally from the pattern seen with R753Q. The clinical priority is managing the downstream consequences of enhanced TLR2 activity: atopic inflammation, skin barrier integrity, and monitoring for inflammatory conditions in which TLR2 overactivation plays a role.

Interactions

Within the TLR2 locus itself, rs4696480 and rs5743708 (R753Q) have been studied together in several atopic dermatitis investigations. These two variants operate through opposite mechanisms — T-16934A modulates expression, R753Q impairs function — yet both associate with severe AD. Carrying both a promoter upregulation variant and a loss-of-function coding variant in the same gene creates a complex phenotype where TLR2 may be abundantly present but functionally impaired, potentially confounding the downstream inflammatory pattern.

The TLR2/FCER1A epistatic interaction (rs4696480 × rs2252226) documented in atopic dermatitis suggests that TLR2 expression level modulates the IgE-mediated arm of allergic inflammation, connecting innate pattern recognition to the adaptive allergic response. This cross-pathway interaction implies that carriers of both variants warrant heightened atopic monitoring.

For gut mucosal immunity, TLR2 operates alongside TLR4 (rs4986790) and CD14 in recognizing the microbial content of the intestinal lumen. Enhanced TLR2 expression from the A allele, combined with other inflammatory pathway variants in NOD2 or IL23R, could amplify mucosal inflammatory responses relevant to IBD phenotype and severity.

A region of chromosome 10 sitting roughly 80 kilobases downstream of the CXCL12 gene¹¹ CXCL12 gene
C-X-C motif chemokine ligand 12, also known as SDF-1 (stromal cell-derived factor 1) harbors one of the most consistently replicated non-obvious GWAS signals for coronary artery disease. rs501120 does not change any protein. What it appears to do is turn up the volume on CXCL12 production — and chronically elevated CXCL12 promotes the vascular inflammation and plaque progression that underlie most heart attacks.

The Mechanism

CXCL12 signals through two receptors, CXCR4 and ACKR3 (CXCR7), orchestrating the recruitment of bone-marrow progenitor cells to sites of vascular injury and regulating inflammatory cell trafficking in arterial walls. At physiological levels it supports vascular repair. When chronically elevated, however, it drives endothelial dysfunction, promotes neointimal hyperplasia, and accelerates atherosclerotic plaque formation²² atherosclerotic plaque formation
CXCL12 derived from endothelial cells promotes atherosclerosis in animal models and human data.

rs501120 sits in a non-coding regulatory region. The T (risk) allele is associated with measurably higher circulating CXCL12 protein — suggesting the locus influences gene expression or mRNA stability, though the precise regulatory element has not been pinpointed. In strong linkage disequilibrium³³ linkage disequilibrium
LD, correlated inheritance, with the nearby rs1746048 (r² ≈ 1.0), both variants tag the same underlying haplotype.

The Evidence

The original signal emerged in the Wellcome Trust Case Control Consortium (WTCCC) and German MI Family Study⁴⁴ Wellcome Trust Case Control Consortium (WTCCC) and German MI Family Study
Samani et al. 2007 NEJM, 1,926 CAD cases, 2,938 controls, where combined analysis reached OR 1.33 (p = 9.46 × 10⁻⁸) — genome-wide significance. This placed the 10q11 locus among the earliest replicated novel GWAS hits for coronary artery disease.

The association held up in a large European replication: Coronary Artery Disease Consortium, 2009⁵⁵ Coronary Artery Disease Consortium, 2009
11,550 cases, 11,205 controls across 9 European studies confirmed OR 1.11 (95% CI 1.05–1.18, p = 4.34 × 10⁻⁴). A notable finding was a sex-specific effect: the risk was significant in women (OR 1.29) but not in men (OR 1.03, p = 0.39), suggesting the variant's cardiovascular impact may be modulated by sex hormones or other sex-linked biology.

Functional support came from Mehta et al. 2011 Eur Heart J⁶⁶ Mehta et al. 2011 Eur Heart J
2,939 participants, two independent cohorts: T allele carriers had higher plasma CXCL12 concentrations (TT 2.34 ± 0.49 vs CC 2.23 ± 0.53 ng/mL, replication p = 0.007, meta-analysis p = 6 × 10⁻⁴), connecting the genetic signal directly to a measurable biochemical intermediate.

Elevated plasma CXCL12 independently predicts adverse outcomes in people already diagnosed with CAD: Ghasemzadeh et al. 2015 Atherosclerosis⁷⁷ Ghasemzadeh et al. 2015 Atherosclerosis
785 CAD patients, 2.6-year follow-up found a 4.8-fold increased risk of cardiovascular death or MI in those with high CXCL12 levels (HR 4.81, p = 1 × 10⁻⁶), improving risk reclassification by 40% beyond traditional risk factors.

Replication in non-European populations adds breadth: a Chinese Han study⁸⁸ Chinese Han study
368 ischemic stroke cases, 381 controls found the C (protective) allele and CT/CC genotypes associated with lower ischemic stroke risk, particularly in men, consistent with the 10q11 locus influencing cerebrovascular as well as coronary endpoints.

Practical Actions

For TT homozygotes (~74% of Europeans), the signal justifies earlier and more intensive cardiovascular monitoring — specifically requesting a high-sensitivity CRP and lipoprotein(a) panel to contextualize the genotype-based risk. Statin therapy has been shown to reduce circulating CXCL12 levels⁹⁹ reduce circulating CXCL12 levels
Dose-dependent statin effect on CXCL12, PMC3560987, providing a mechanistic rationale, in addition to their established LDL-lowering benefit, for prioritizing statin therapy discussions in TT individuals at intermediate cardiovascular risk who might otherwise be in a "watchful waiting" category.

Nitrate-rich vegetables (beetroot, leafy greens) support endothelial nitric oxide production via a CXCL12-independent pathway, providing a dietary complement to pharmacological approaches for people with the TT genotype. Avoid high-dose supplement regimens marketed as "boosting stem cell mobilization" — several promote CXCL12/CXCR4 activation, which would be counterproductive in the context of this variant.

Interactions

rs501120 is in near-perfect LD (r² ≈ 1.0) with rs1746048, another CXCL12 downstream variant with an independent GWAS signal, meaning the two variants effectively tag the same haplotype and compound actions between them are not meaningful — they will nearly always be inherited together.

rs1801157 in the 3′UTR of CXCL12 has been investigated separately; a meta-analysis found it is not independently associated with CAD risk, though it may influence CXCL12 mRNA stability by a different mechanism.

Your FUT2 gene determines secretor status — one of the most consequential genetic traits in human biology, influencing not just your own gut immunity but the microbial world you pass on to your children. Secretors express ABO blood group antigens¹¹ Secretors express ABO blood group antigens
H-type antigens built by alpha(1,2)-fucosyltransferase — the FUT2 enzyme — on mucosal surfaces and in saliva, tears, and breast milk. These fucosylated glycans are absent in non-secretors, who carry two loss-of-function copies of FUT2 in their intestinal mucus, saliva, and breast milk. rs516246 is an intronic proxy variant that tags this biological divide in European and African populations through strong linkage disequilibrium with the functional W143X nonsense mutation (rs601338). The T allele travels with the secretor phenotype; the C allele travels with the non-secretor phenotype.

This platform already catalogues rs601338 (the functional variant, focused on vitamin B12 metabolism) and rs492602 (a synonymous proxy focused on Crohn's disease and the gut-skin axis). rs516246 adds a distinct clinical dimension: the maternal-infant axis²² maternal-infant axis
The biological channel through which a mother's FUT2 genotype shapes the gut microbiome colonisation of her breastfed infant via the fucosylated oligosaccharides in her breast milk. This angle has direct implications for infant microbiome development and long-term immune programming.

The Mechanism

The FUT2 enzyme adds fucose residues to glycan chains on the intestinal epithelium and into secreted fluids, creating H-type 1 antigens³³ H-type 1 antigens
The carbohydrate structure H type 1 is the intestinal and secretory form of the H antigen, built on Type 1 precursor chains; H type 2 is expressed on red blood cells. Only FUT2 generates H type 1 on mucosal surfaces and in milk. In breast milk, secretor mothers produce abundant 2'-fucosyllactose (2'-FL) and lacto-N-fucopentaose I (LNFP I)⁴⁴ 2'-fucosyllactose (2'-FL) and lacto-N-fucopentaose I (LNFP I)
The two most abundant human milk oligosaccharides (HMOs) produced exclusively by secretor mothers; 2'-fucosyllactose is the most prevalent HMO in secretor breast milk at concentrations of 1.5–4 g/L, collectively termed 2'-fucosylated human milk oligosaccharides (HMOs). Non-secretor mothers produce no 2'-fucosylated HMOs.

These fucosylated HMOs are the primary carbon source for Bifidobacterium species in the breastfed infant's gut — particularly B. longum subsp. infantis, which encodes a complete suite of fucosidase enzymes to harvest fucose from HMOs. Without 2'-FL and LNFP I in the milk, the infant's developing gut receives a fundamentally different microbial substrate, shifting which Bifidobacterium species flourish and which metabolic pathways dominate early colonisation.

The Evidence

The maternal secretor-infant microbiome connection was characterised by Lewis et al. 2015⁵⁵ Lewis et al. 2015
Lewis ZT et al. Maternal fucosyltransferase 2 status affects the gut bifidobacterial communities of breastfed infants. Microbiome, 2015: infants of secretor mothers established Bifidobacterium colonisation earlier and at significantly higher levels (p<0.001), with B. longum dominant in secretor-fed infants and B. breve more prevalent in non-secretor-fed infants — a difference in species composition with distinct metabolic and immune consequences.

Durham et al. 2021⁶⁶ Durham et al. 2021
Durham SD et al. A one-year study of human milk oligosaccharide profiles in the milk of healthy UK mothers and their relationship to maternal FUT2 genotype. Glycobiology, 2021 used rs516246 directly to genotype maternal secretor status in a longitudinal cohort, tracking HMO profiles across twelve months of lactation. The study confirmed that maternal FUT2 genotype at rs516246 reliably predicts 2'-fucosylated HMO composition, though it also noted occasional phenotype-genotype discordance, suggesting that enzyme expression variation can modulate predicted HMO output.

From the infant's own perspective, Thorman et al. 2023⁷⁷ Thorman et al. 2023
Thorman AW et al. Gut Microbiome Composition and Metabolic Capacity Differ by FUT2 Secretor Status in Exclusively Breastfed Infants. Nutrients, 2023 found that infant FUT2 secretor status was an even stronger driver of early microbiome composition than maternal secretor status. Full-secretor infants showed greater alpha diversity (p=0.049) and distinct metabolic profiles, with non-secretor infants displaying enriched sucrose-degradation pathways and different enterotypes.

Beyond infant biology, the broader FUT2 literature documents substantial effects on adult immunity. Non-secretors face elevated Crohn's disease risk⁸⁸ Non-secretors face elevated Crohn's disease risk
McGovern DPB et al. Fucosyltransferase 2 (FUT2) non-secretor status is associated with Crohn's disease. Hum Mol Genet, 2010 at genome-wide significance (P=4.90×10⁻⁸) and increased type 1 diabetes susceptibility OR 1.29, 95% CI 1.20–1.37⁹⁹ OR 1.29, 95% CI 1.20–1.37
Smyth DJ et al. FUT2 nonsecretor status links type 1 diabetes susceptibility and resistance to infection. Diabetes, 2011. The trade-off is near-complete protection against symptomatic norovirus GII infection — not a single non-secretor developed symptomatic norovirus¹⁰¹⁰ not a single non-secretor developed symptomatic norovirus
Thorven M et al. A homozygous nonsense mutation (428G→A) in FUT2 provides resistance to symptomatic norovirus (GGII) infections. J Virol, 2005 among outbreak-exposed adults in the Thorven 2005 study.

Practical Actions

For non-secretors (CC), the most targeted actions address the deficit in gut microbial diversity. Without host-derived fucosylated glycans in the gut lumen, dietary prebiotic fibers serve as an alternative carbon source for Bifidobacterium species. Targeted supplementation with multiple Bifidobacterium strains can partially compensate for reduced colonisation.

For secretor mothers (CT or TT) breastfeeding non-secretor infants, or non-secretor mothers who cannot produce 2'-fucosylated HMOs, awareness of the downstream infant microbiome implications is clinically relevant. Some infant formula products now supplement with 2'-fucosyllactose, and this genotype-driven difference in milk composition may influence breastfeeding supplementation decisions in consultation with a paediatrician.

Vitamin B12 monitoring is also relevant for non-secretors: standard total serum B12 may overestimate functional status because elevated haptocorrin-bound B12 inflates the total measurement without increasing bioavailable transcobalamin-bound B12.

Interactions

rs516246 is an intronic proxy in strong LD with rs601338 (FUT2 W143X, the functional variant) and rs492602 (another synonymous proxy) in European and African populations. Carrying the C allele at rs516246 very likely means also carrying the A (non-secretor) allele at rs601338. The three entries address distinct clinical dimensions of the same biology: rs601338 covers B12 metabolism, rs492602 covers Crohn's disease and the gut-skin axis, and this entry covers maternal-infant HMO biology.

In East Asian populations, the primary non-secretor determinant is rs1047781 (A385T), since the W143X-tagged haplotype is nearly absent in that ancestry.

The OXTR gene encodes the oxytocin receptor¹¹ oxytocin receptor
A G-protein coupled receptor expressed throughout the brain, uterus, and cardiovascular system that mediates the effects of the neuropeptide oxytocin, the protein through which the neuropeptide oxytocin exerts its wide-ranging effects on social bonding, empathy, trust, and stress regulation. Oxytocin is sometimes called the "love hormone," but its biology is far more nuanced than that label suggests — it modulates social salience, making social cues more prominent, for better or worse.

The rs53576 variant is a common A-to-G polymorphism in intron 3²² intron 3
An intron is a non-coding region within a gene. While it doesn't change the protein sequence, intronic variants can affect gene expression by altering regulatory elements, mRNA splicing, or chromatin structure of the OXTR gene on chromosome 3. Despite not directly altering the receptor protein, it is the single most studied variant in the oxytocin system, with over 245 published studies linking it to differences in empathy, stress reactivity, social behavior, and mental health outcomes. The G allele is generally associated with enhanced social sensitivity and greater benefit from social support, while the A allele is associated with reduced empathy scores, lower parental sensitivity, and diminished stress buffering from social connections.

The Mechanism

As an intronic variant, rs53576 does not change the amino acid sequence of the oxytocin receptor itself. Its functional effects are thought to arise through regulatory mechanisms — potentially influencing OXTR gene expression levels, mRNA stability, or epigenetic modification³³ epigenetic modification
DNA methylation at the OXTR locus has been shown to affect receptor expression; rs53576 genotype may influence susceptibility to methylation changes that alter how much receptor protein is produced. The variant is in complete linkage disequilibrium⁴⁴ linkage disequilibrium
LD means two genetic variants are inherited together so frequently that knowing one genotype effectively predicts the other with rs4686302, a missense variant in OXTR that causes a Thr-to-Met amino acid change — raising the possibility that rs53576 is a marker for a functional change at this nearby site.

Neuroimaging studies⁵⁵ Neuroimaging studies
Tost H et al. A common allele in the oxytocin receptor gene impacts prosocial temperament and human hypothalamic-limbic structure and function. PNAS, 2010 have shown that A-allele carriers have altered hypothalamic and amygdala structure and function. Male A-allele carriers show reduced hypothalamic volume and increased amygdala volume compared to GG carriers, and these structural differences predict lower scores on prosocial temperament measures. These findings suggest the variant shapes the neural architecture underlying social cognition.

The Evidence

The landmark Rodrigues et al. 2009 study⁶⁶ Rodrigues et al. 2009 study
Rodrigues SM et al. Oxytocin receptor genetic variation relates to empathy and stress reactivity in humans. PNAS, 2009 first established the behavioral significance of rs53576 in 192 participants. GG homozygotes were 22.7% less likely to make errors on the Reading the Mind in the Eyes Test (a measure of empathic accuracy) and showed lower heart-rate reactivity during a startle anticipation task compared to A-allele carriers.

Saphire-Bernstein et al. (2011)⁷⁷ Saphire-Bernstein et al. (2011)
Saphire-Bernstein S et al. Oxytocin receptor gene is related to psychological resources. PNAS, 2011 extended these findings to psychological resources in 348 participants: A-allele carriers had lower optimism, self-esteem, and mastery, along with higher depressive symptomatology. The effect on depression appeared to be mediated by reduced psychological resources.

A pivotal study by Chen et al. (2011)⁸⁸ Chen et al. (2011)
Chen FS et al. Common oxytocin receptor gene polymorphism and social support interact to reduce stress in humans. PNAS, 2011 demonstrated genotype-dependent social buffering in 194 men. G-allele carriers who received social support before a psychosocial stress test showed significantly lower cortisol and subjective stress responses, while AA homozygotes derived much less benefit from the same social support. This is one of the clearest demonstrations that rs53576 modulates the stress- protective effects of social connection.

The Li et al. 2015 meta-analysis⁹⁹ Li et al. 2015 meta-analysis
Li J et al. Association of OXTR rs53576 polymorphism with sociality: a meta-analysis. PLoS ONE, 2015 pooled 24 samples (n=4,955) and confirmed that GG homozygotes show greater general sociality than A-allele carriers (Cohen's d=0.11). However, the effect was specific to general social behavior and did not extend to close relationships, suggesting rs53576 primarily affects broader social orientation rather than intimate bonding.

A 2021 systematic review by Chander et al.¹⁰¹⁰ 2021 systematic review by Chander et al.
Chander RJ et al. The influence of rs53576 polymorphism in the OXTR gene on empathy in healthy adults by subtype and ethnicity. Psychoneuroendocrinology, 2021 found that the GG-empathy association was significant primarily in young to middle-aged adults and showed differential effects by ethnicity, with stronger cognitive empathy differences in Asian cohorts.

Practical Implications

This is fundamentally a gene-environment variant. The rs53576 genotype does not operate in isolation — its effects are consistently modulated by the social environment. G-allele carriers appear to be more socially sensitive in both positive and negative directions: they benefit more from social support but are also more affected by social adversity. Bradley et al. (2013)¹¹¹¹ Bradley et al. (2013)
Bradley B et al. Association between childhood maltreatment and adult emotional dysregulation: moderation by oxytocin receptor gene. Dev Psychopathol, 2013 found that GG carriers exposed to severe childhood maltreatment showed greater emotional dysregulation than A carriers — consistent with a differential susceptibility model where the G allele amplifies environmental influence rather than simply being "better."

For GG and AG individuals, the practical takeaway is that social connection is not just pleasant but physiologically protective. Investing in close relationships, seeking support during stress, and maintaining social engagement may be especially important for stress management. For AA individuals, the biology suggests that solitary stress-management strategies (exercise, mindfulness, structured routines) may be relatively more effective than relying primarily on social support.

Population frequencies vary dramatically across ancestries. The A allele predominates in East Asian populations (~65%), while the G allele predominates in European (~68%) and African (~77%) populations. Cultural factors interact with these genetic differences: Kim et al. (2010)¹²¹² Kim et al. (2010)
Kim HS et al. Culture, distress, and oxytocin receptor polymorphism interact to influence emotional support seeking. PNAS, 2010 showed that the GG-genotype association with emotional support seeking under distress appeared in American but not Korean participants, suggesting that cultural norms modulate how genetic sensitivity is expressed behaviorally.

Interactions

OXTR rs53576 likely interacts with COMT rs4680 (Val158Met) in shaping social-emotional phenotypes. COMT determines dopamine clearance speed in the prefrontal cortex — slow COMT (Met/Met) increases baseline dopamine and emotional sensitivity, while fast COMT (Val/Val) clears dopamine rapidly. An individual carrying both OXTR GG (high social sensitivity) and COMT Met/Met (high emotional reactivity) may experience amplified responses to social environments, both positive and negative. Conversely, OXTR AA combined with COMT Val/Val could produce a profile of relative emotional and social resilience. While this interaction has theoretical grounding in overlapping neurocircuitry, direct gene-gene interaction studies at the rs53576-by-rs4680 level are preliminary.

Von Willebrand factor (VWF) is the molecular bridge between injured blood vessel walls and circulating platelets. Without adequate VWF, even a minor cut can lead to prolonged bleeding. The D141G variant — a single amino acid substitution replacing aspartate with glycine at position 141 of the VWF protein — was catalogued by the International Society on Thrombosis and Haemostasis VWF mutation database and identified in a UK cohort study of type 1 VWD families¹¹ UK cohort study of type 1 VWD families
Cumming et al., Thromb Haemost 2006 — VWF gene sequenced in 32 confirmed UK type 1 VWD families as likely causative of von Willebrand disease in index cases. The variant is extremely rare in the general population and is not well-represented in large population databases such as gnomAD, consistent with a disease-causing rather than neutral polymorphism.

The Mechanism

Position 141 of VWF falls within the D1 propeptide domain²² D1 propeptide domain
The VWF propeptide (D1-D2) is cleaved after secretion and is essential for correct disulfide-bond formation and multimerisation of VWF. Aspartate at position 141 participates in the structural organisation of this domain; substitution with glycine — a conformationally flexible but electrically neutral residue — disrupts the local protein fold. Both SIFT (score 0, deleterious) and PolyPhen-2 (score 0.996, probably damaging)³³ SIFT (score 0, deleterious) and PolyPhen-2 (score 0.996, probably damaging)
Ensembl VEP predictions for NM_000552.5:c.422A>G classify D141G as likely deleterious. The consequence may be impaired multimerisation, reduced secretion efficiency, or accelerated clearance of VWF — mechanisms common to type 1 VWD missense variants and collectively resulting in lower circulating VWF:Ag and VWF activity (VWF:RCo).

The VWF gene sits on chromosome 12 (12p13.31) and is transcribed off the minus strand. On a genome file, this variant is reported as T→C on the plus strand (the C allele encoding the Gly141 substitution when translated from the minus-strand coding sequence).

The Evidence

The MCMDM-1VWD European cohort⁴⁴ MCMDM-1VWD European cohort
Goodeve et al. Blood 2007 — 150 type 1 VWD index cases across 9 European countries found that 70% of type 1 VWD index cases carried identifiable VWF mutations, with missense variants concentrated in the D and C domains responsible for multimerisation and platelet binding. Penetrance is incomplete: some D141G carriers may have VWF levels within the low-normal range (30–50 IU/dL) and minimal symptoms, while others present with VWF:Ag well below 30 IU/dL and frank bleeding diathesis.

The ASH/ISTH/NHF/WFH 2021 joint guidelines⁵⁵ ASH/ISTH/NHF/WFH 2021 joint guidelines
James et al. Blood Adv 2021 — 11 evidence-based recommendations using GRADE methodology define type 1 VWD as VWF:Ag <30 IU/dL with a positive bleeding history, and recommend structured bleeding assessment tools (ISTH-BAT) for all suspected VWD cases. Carriers of D141G who have never had laboratory testing should have VWF:Ag, VWF:RCo, and FVIII measured to determine their functional phenotype.

Atiq et al. (Blood Adv 2019)⁶⁶ Atiq et al. (Blood Adv 2019)
122 type 1 VWD patients followed prospectively; desmopressin challenge correlated with subsequent bleeding phenotype demonstrated that the magnitude of VWF and factor VIII rise after desmopressin challenge predicts clinical bleeding severity: patients achieving the highest FVIII quartile at 3 hours post-dose had substantially lower Tosetto bleeding scores. This means the desmopressin challenge test is both diagnostic and prognostic for type 1 VWD carriers.

Practical Actions

Any carrier of D141G who has not had haematological evaluation should request VWF panel testing (VWF:Ag, VWF:RCo, FVIII:C) from their physician and, if results are borderline or low, a referral to a haematologist for desmopressin challenge testing. Before surgical or dental procedures, excessive bleeding should be anticipated and preventive measures discussed with the treating team. Tranexamic acid (an anti-fibrinolytic agent) is effective for mucosal bleeding episodes (nosebleeds, heavy menstrual bleeding) in mild-moderate VWD and can often be managed without haematology input once the diagnosis is confirmed. Desmopressin (DDAVP) is first-line for type 1 VWD when a documented response is confirmed on prior challenge testing.

Interactions

VWF interacts genetically with ABO blood group: O blood group individuals naturally have approximately 25% lower VWF:Ag than non-O groups, which can unmask borderline VWF variants. Carriers of D141G with blood group O may have a more symptomatic phenotype than those with A, B, or AB blood groups. Additionally, high VWF clearance (measurable as elevated VWFpp/VWF:Ag ratio) compounds the functional deficit in type 1 VWD regardless of the causal variant; some D141G carriers may have both a synthesis defect and enhanced clearance contributing to their low VWF levels.

The kidney's glomerular filter depends on an intricate architecture of specialized cells called podocytes¹¹ podocytes
Podocytes: highly differentiated epithelial cells that wrap around glomerular capillaries with interdigitating foot processes connected by slit diaphragms, forming the final barrier against protein leakage into urine. rs62408925 sits in the intergenic region at chromosome 3q22.3, in the vicinity of NCK1 — a molecular adaptor that physically connects the podocyte's slit diaphragm protein nephrin to its actin cytoskeleton. Genetic variation at this locus has been linked to elevated diabetic nephropathy (DN) risk, the most common serious complication of type 1 diabetes and a leading cause of end-stage renal disease worldwide.

The Mechanism

rs62408925 is not itself a coding variant — it sits in an intergenic region and does not change any amino acid sequence. It marks one boundary of an 11 kilobase segment that contains three highly conserved regulatory elements, and it is in strong linkage disequilibrium²² strong linkage disequilibrium
Linkage disequilibrium (LD): the non-random co-inheritance of nearby alleles on the same chromosome. r²=0.95 means the two variants are almost always co-inherited — knowing one nearly perfectly predicts the other (r²=0.95) with rs1866813, the functional variant at this locus.

Mechanistic studies³³ Mechanistic studies
He B et al. A remote cis-acting variant at 3q links glomerular NCK1 to diabetic nephropathy. PLoS One, 2013 show that rs1866813 operates as a remote cis-regulatory element approximately 70 kb upstream of the NCK1 transcription start site. The risk allele drives higher NCK1 expression specifically in glomerular tissue. NCK1 encodes an adaptor protein with SH2 and SH3 domains that binds phosphorylated nephrin at the slit diaphragm and recruits actin polymerization machinery — it is the molecular bridge⁴⁴ it is the molecular bridge
Jones N et al. Nck adaptor proteins link nephrin to the actin cytoskeleton of kidney podocytes. Nature, 2006 between the filtration barrier's membrane architecture and the cytoskeletal scaffolding that holds foot processes in shape.

Under hyperglycemia, podocytes are already under metabolic and oxidative stress. Elevated NCK1 expression may dysregulate the tight balance of nephrin-actin signaling needed to maintain foot process architecture. Adult podocyte-specific NCK deletion in mice rapidly produced proteinuria, glomerulosclerosis, and foot process effacement⁵⁵ rapidly produced proteinuria, glomerulosclerosis, and foot process effacement
Jones N et al. Nck proteins maintain the adult glomerular filtration barrier. J Am Soc Nephrol, 2009 within 1–2 weeks — demonstrating that NCK dosage dysregulation in either direction can compromise the filtration barrier.

The Evidence

The initial genetic association came from a multistage case-control study⁶⁶ multistage case-control study
He B et al. Association of genetic variants at 3q22 with nephropathy in patients with type 1 diabetes mellitus. AJHG, 2009 of 1,822 diabetic nephropathy cases and 1,874 controls from Finland, Iceland, and the British Isles. rs1866813 (r²=0.95 with rs62408925) showed a combined odds ratio of 1.33 for DN, with the Finnish replication cohort yielding OR=1.38 (95% CI 1.18–1.62, p=4.7×10⁻⁵). The association was codominant, meaning each risk allele copy incrementally increased risk.

The 2013 mechanistic paper provided supporting biological plausibility — allele-specific luciferase reporter assays in lymphocytes and transgenic zebrafish with podocyte-specific GFP reporters confirmed differential expression driven by the two alleles. This molecular evidence substantially strengthens the interpretation beyond a statistical association alone.

The evidence level is rated moderate: the original association was replicated within the same 2009 study across populations, and mechanistic validation exists, but independent large-scale replication by separate groups is not established in the literature. The locus has not appeared in the largest diabetic nephropathy GWAS meta-analyses to date, possibly reflecting heterogeneity across type 1 versus type 2 diabetes cohorts or statistical power constraints.

Practical Actions

For carriers of the CT or TT genotype with type 1 diabetes, the primary actionable implication is earlier and more frequent kidney function monitoring. Diabetic nephropathy progresses through stages detectable years before clinical kidney impairment — microalbuminuria (UACR 30–300 mg/g) is the earliest measurable sign. Early detection enables RAAS blockade (ACE inhibitors or ARBs) that substantially slows DN progression.

For the roughly 86% of people who carry the common CC genotype, standard nephropathy screening guidelines for type 1 diabetes apply (annual UACR from 5 years after diagnosis).

For people without type 1 diabetes, this locus has no identified clinical relevance — the association is specific to the diabetic hyperglycemic milieu that places podocyte stress on NCK1 dysregulation.

Interactions

The neighboring gene IL20RB at 3q22.3 encodes the IL-20 receptor beta subunit, which participates in cytokine signaling relevant to kidney inflammation. IL20RB protein has been detected in biopsies from diabetic nephropathy, IgA nephropathy, and lupus nephritis patients, raising the possibility that the 3q22 locus contains multiple independent regulatory signals. The LD structure of the region means rs62408925 and rs1866813 may co-tag both NCK1 and IL20RB regulatory elements, though this has not been mechanistically resolved.

Diabetic nephropathy risk is substantially modified by glycemic control — the hyperglycemic environment is the necessary cofactor for NCK1-related podocyte injury to manifest. Variants affecting insulin secretion (TCF7L2 rs7903146) or insulin sensitivity (PPARG rs1801282) are independent DN risk modifiers through glucose control pathways.

Approximately 185 kilobases upstream of the TNFAIP3 gene sits a regulatory region that controls how much A20 protein your immune cells produce. A20, encoded by TNFAIP3¹¹ A20, encoded by TNFAIP3
TNFAIP3 stands for TNF Alpha Induced Protein 3; A20 is its common protein name and a key negative regulator of NF-kB signaling is the primary brake on NF-kB-driven inflammation — the pathway that amplifies immune responses after infection or injury. The rs6920220 variant resides in an intergenic region between OLIG3 and TNFAIP3 at chromosome 6q23.3, and the A allele reduces TNFAIP3 transcription²² reduces TNFAIP3 transcription
CRISPR-Cas9 editing demonstrated that rs6920220 A-allele cells show significant downregulation of TNFAIP3 mRNA compared to G-allele controls, allowing inflammatory signals to persist longer and reach higher intensities than they should.

The Mechanism

The 6q23 intergenic region around rs6920220 contains putative transcriptional regulatory elements. Luciferase reporter assays³³ Luciferase reporter assays
In vitro functional studies using reporter constructs in T lymphoblastoid cell lines demonstrated repressor activity at rs6920220 and two neighboring SNPs in high linkage disequilibrium, confirming that this region actively modulates TNFAIP3 gene expression. When the A allele is present, this repressor activity is altered, and TNFAIP3 mRNA levels fall.

The downstream consequences are measurable: CRISPR-engineered salivary gland epithelial cells⁴⁴ CRISPR-engineered salivary gland epithelial cells
Cells with the A allele introduced by CRISPR showed markedly increased NF-κB mRNA levels, elevated IL-6 and IL-8 expression, and increased IL-1β in co-culture with immune cells. Since A20 normally terminates NF-kB signaling after immune activation, reduced A20 expression means inflammatory cascades run longer and produce more cytokines before shutting down. This is mechanistically distinct from the missense variant rs2230926 (F127C), which reduces A20 enzymatic activity — rs6920220 reduces how much A20 is made in the first place.

The Evidence

The rs6920220 A allele was first linked to rheumatoid arthritis through the Wellcome Trust Case Control Consortium GWAS⁵⁵ Wellcome Trust Case Control Consortium GWAS
Discovery was in a genome-wide association study subsequently replicated in independent cohorts and achieved unequivocal replication with P=1.1×10⁻⁸ and OR=1.22 (95% CI 1.15-1.33). A meta-analysis of 21 case-control studies⁶⁶ meta-analysis of 21 case-control studies
Included stratified analysis by ethnicity confirming Caucasian-specific effect confirmed OR 1.36 (95% CI 1.24–1.50, P<0.001) for homozygous AA carriers versus GA+GG in the overall analysis, with stratified analysis showing Caucasian-specific risk (OR 1.37, 95% CI 1.24–1.51 in the recessive model). There is no significant RA association in East Asian populations for this SNP, in contrast to other 6q23 variants.

At 6q23, rs6920220 is one of three independent RA risk signals⁷⁷ three independent RA risk signals
Conditional logistic regression identified three independent associations at 6q23: rs6920220 (risk), rs13207033 (protective), and rs5029937 (risk). Carrying both A alleles of rs6920220 and rs5029937 while lacking the protective allele of rs13207033 raises the combined OR to 1.86, substantially amplifying the individual effects.

The variant is associated with faster radiological joint destruction⁸⁸ faster radiological joint destruction
Median Larsen radiological damage scores 31 (GG) vs 36 (GA/AA), P=0.02, in autoantibody-positive RA, specifically in autoantibody-positive patients — suggesting this variant exacerbates the inflammatory joint damage characteristic of seropositive disease.

Associations extend beyond RA: rs6920220 was associated with SLE⁹⁹ rs6920220 was associated with SLE
Two independent signals near TNFAIP3 included rs6920220 with P=0.03, confirming shared autoimmune susceptibility across conditions in an Italian cohort with OR 1.53 for SLE and OR 1.69 for primary Sjogren's syndrome. The variant also conferred JIA risk¹⁰¹⁰ conferred JIA risk
OR 1.30 (95% CI 1.05-1.61), P=0.015, particularly for oligoarticular JIA, confirming its shared autoimmune susceptibility across conditions with pediatric onset.

The variant shows striking population stratification¹¹¹¹ population stratification
A allele frequency 21% in Europeans, 0.2% in East Asians, 11% in Africans by dbSNP/ALFA consortium data: high frequency in Europeans, low in East Asians. RA associations with rs6920220 are confirmed in European populations but not Asian populations, which carry the A allele at less than 0.3% frequency.

Practical Implications

If you carry one or two A alleles, your cells produce somewhat less A20 protein — the primary brake on NF-kB inflammatory signaling. This creates a genetically primed state for inflammatory amplification when immune triggers arise. For most carriers, this manifests as modestly elevated risk for autoimmune diseases rather than active disease.

The most clinically actionable implication involves anti-TNF biologic therapy. For carriers who develop autoimmune conditions — particularly rheumatoid arthritis — preliminary evidence from psoriatic arthritis cohorts¹²¹² preliminary evidence from psoriatic arthritis cohorts
Spanish PsA cohort (n=20) found rs6920220 significantly associated with quality of life improvement at 3 and 6 months with TNF inhibitor treatment suggests this variant may influence how well TNF inhibitors work. The biological logic is coherent: since reduced A20 drives excess NF-kB activity and the TNF signaling pathway feeds directly into NF-kB, blocking TNF upstream could be particularly relevant for carriers. However, this evidence is from a small cohort and should be treated as emerging rather than established guidance. Sharing your genotype with your rheumatologist provides useful context when selecting between biologics with different mechanisms.

The variant's association with more severe radiological damage in autoantibody-positive RA (higher Larsen scores) underscores the value of early and aggressive treatment for A-allele carriers who develop seropositive disease.

Interactions

rs6920220 is one of three independent association signals at the 6q23 locus. The other two variants — rs5029937 (intronic in TNFAIP3, also risk-conferring) and rs13207033 (protective) — are independently inherited and their combined effect can substantially raise or lower RA risk beyond rs6920220 alone.

The TNFAIP3 missense variant rs2230926 (F127C) impairs A20 enzymatic activity through a completely independent mechanism: rs6920220 reduces how much A20 is made; rs2230926 makes A20 less effective at its job. Carriers of both the rs6920220 A allele and rs2230926 G allele could have both impaired A20 expression and impaired A20 function, representing compounded NF-kB dysregulation.

PTPN22 R620W (rs2476601) and TNF-alpha -308 (rs1800629) define the broader genetic context for autoimmune and anti-TNF response risk — each operating through distinct mechanisms that converge on T-cell activation and inflammatory amplification.

Uterine fibroids (leiomyomas) are benign smooth-muscle tumors of the uterus that affect up to 70% of women by age 50. While they are common, their underlying biology is clonal — each fibroid arises from a single cell that escapes normal growth control. A genome-wide association study of more than 400,000 women in the UK Biobank pinpointed rs72709458 in the TERT gene as one of the strongest genetic risk factors for leiomyoma, with a genome-wide significant odds ratio of 1.12 (p = 6.9 × 10⁻¹⁶). This is not the same biological story as TERT longevity variants — this variant acts through a different mechanism focused on the uterine smooth-muscle cell's capacity to maintain genome integrity¹¹ uterine smooth-muscle cell's capacity to maintain genome integrity
Telomere maintenance failure allows chromosomal instability, which can initiate clonal tumor growth.

TERT encodes the catalytic subunit of telomerase, the enzyme that rebuilds telomere caps on chromosome ends after cell division. Without adequate telomerase activity, telomeres shorten with each cell division until they trigger a DNA damage response — usually stopping the cell from dividing further, but sometimes allowing chromosomal instability and aberrant clonal expansion instead.

The rs72709458 variant sits in an intron of TERT at chromosome 5, position 1,283,640 (GRCh38). It does not change the protein sequence directly; rather, it is a regulatory tag variant in strong linkage disequilibrium with nearby functional elements that influence TERT expression levels or splicing in uterine tissue. Importantly, the three LD-independent TERT associations with leiomyoma risk (rs72709458, rs2736100, rs2853676) each contribute independently, and their combined signal points to the same biological axis: telomere maintenance capacity in the uterine myometrium²² telomere maintenance capacity in the uterine myometrium
The myometrium is estrogen-responsive, undergoes repeated cycles of growth and regression, and depends on intact genome surveillance to prevent clonal outgrowth. The TERT, TERC, and OBFC1 gene cluster together represent the telomere maintenance axis in the leiomyoma GWAS, grouping with TP53 and ATM as a "genome stability" gene cluster distinct from the "genitourinary development" gene cluster (WT1, ESR1/SYNE1, CDC42/WNT4).

Notably, rs72709458 and rs2853676 at TERT did not independently associate with blood leukocyte telomere length in this study — unlike rs2736100, which was significantly associated with shorter telomere length in tumor tissue. This suggests rs72709458 may act through a tissue-specific or context-specific mechanism rather than globally shortening telomeres. The combined effect of all three TERT SNPs plus TERC and OBFC1 variants showed a negative trend toward shorter telomeres (p = 0.055), consistent with a cumulative model.

Välimäki et al. (2018, eLife)³³ Välimäki et al. (2018, eLife)
Genetic predisposition to uterine leiomyoma is determined by loci for genitourinary development and genome stability conducted the largest uterine leiomyoma GWAS to date using the UK Biobank, identifying 22 genome-wide significant loci. The rs72709458 T allele at the TERT locus showed OR = 1.12, p = 6.9 × 10⁻¹⁶ — a highly reproducible, statistically unambiguous association. The effect allele frequency of approximately 0.21 means that a substantial fraction of women carry at least one T allele.

A subsequent Russian case-control study⁴⁴ A subsequent Russian case-control study
GWAS-Significant Loci and Uterine Fibroids Risk: Analysis of Associations, Gene-Gene and Gene-Environmental Interactions genotyped 737 hospitalized women with uterine fibroids and 451 controls for seven GWAS-validated SNPs including rs72709458. The study identified synergistic gene-gene interactions among these loci and a gene-environment interaction between the TERT variant and prior abortion history — reinforcing the view that rs72709458 acts within a broader biological context of uterine tissue stress and cell division history.

Vitamin D merits specific mention as a modifiable factor in fibroid biology. A comprehensive mechanistic review⁵⁵ A comprehensive mechanistic review
The Role of Vitamin D in Uterine Fibroid Biology, Baird & Wise 2015 documents that vitamin D (1,25-dihydroxyvitamin D3) exerts direct anti-fibroid effects: it induces G1 cell cycle arrest in leiomyoma cells, downregulates anti-apoptotic proteins (Bcl-2, Bcl-xL), reduces extracellular matrix remodeling enzymes (MMP-2/9), and suppresses estrogen and progesterone receptor expression. Three large observational studies found inverse associations between serum vitamin D levels and fibroid prevalence, with adjusted ORs ranging from 0.68 to 2.5 depending on the study. African American women — who have 2–3× higher fibroid incidence and carry the rs72709458 T allele at a frequency of ~14.5% — are also 10× more likely to be vitamin D deficient. Whether optimizing vitamin D status specifically attenuates fibroid risk in genetically predisposed women has not been formally tested in a randomized trial, but the mechanistic and epidemiological rationale is substantial.

The rs72709458 T allele also has prior associations to endometriosis, lung adenocarcinoma, and glioma — illustrating the broad neoplasia-predisposing context of TERT variation. This pleiotropy is consistent with the telomere maintenance model: cells that divide frequently in hormone-responsive or inflammatory contexts (uterine myometrium, lung epithelium, neural tissue) are particularly dependent on adequate TERT function.

For women who carry the T risk allele at rs72709458, the key message is that this genetic signal warrants increased awareness and proactive monitoring — not alarm. Uterine fibroids are common, usually benign, and often manageable. The OR of 1.12 per allele confers a modest individual risk increment; the variant matters most when combined with other risk factors (family history of fibroids, Black ethnicity, obesity, early menarche, low vitamin D).

The most evidence-supported modifiable interventions targeting fibroid biology include: - Optimizing vitamin D status: The mechanistic evidence is strong and the intervention is low-risk. Maintaining 25(OH)D levels ≥ 40 ng/mL is a reasonable target. - Early symptom reporting: Many fibroids are asymptomatic; those that become symptomatic cause heavy menstrual bleeding, pelvic pressure, or fertility impairment. Earlier diagnosis allows more conservative management options (myomectomy, UAE, progesterone IUDs). - Pelvic ultrasound awareness: Women with a genetic or family history of fibroids can discuss with their clinician whether surveillance ultrasound is appropriate, particularly in the reproductive years.

rs2736100 (TERT, intron 2): The strongest of the three LD-independent TERT signals for leiomyoma, and the one most directly associated with shorter telomeres in tumor tissue. If you carry risk alleles at both rs72709458 and rs2736100, the combined effect on TERT-related genome stability may be greater than either alone. A compound action analysis across these two variants would capture the combined telomere maintenance burden.

rs2853676 (TERT): A third independent TERT signal for leiomyoma risk, also not individually associated with telomere length. Its biological relationship to rs72709458 is not yet fully characterized.

rs12696304 (TERC): The RNA component of telomerase, which acts with TERT to maintain telomeres. Combined risk alleles at TERT and TERC loci showed a negative trend toward shorter telomeres (p = 0.055), suggesting an additive model of telomere maintenance failure.

CD226 (also called DNAM-1, DNAX Accessory Molecule-1)¹¹ CD226 (also called DNAM-1, DNAX Accessory Molecule-1)
A co-stimulatory receptor expressed on T cells, natural killer cells, and myeloid cells. CD226 binds the ligands PVR (CD155) and Nectin-2 (CD112) on antigen-presenting cells and cancer cells, delivering a secondary activation signal that lowers the threshold for T-cell and NK-cell cytotoxicity occupies a central position in both antitumour immunity and autoimmune regulation. rs763361 is a missense variant in exon 7 of CD226 on chromosome 18q22 that substitutes glycine for serine at protein position 307. Unlike loss-of-function variants, the Ser307 form of the protein alters co-stimulatory signalling in a way that has been consistently — across more than 50 independent studies — associated with increased susceptibility to a broad spectrum of autoimmune diseases.

The Mechanism

The Gly307Ser substitution lies in the cytoplasmic tail of CD226, near a tyrosine-based signalling motif²² tyrosine-based signalling motif
A short amino acid sequence (typically YxxL or similar) within a receptor's intracellular domain that, when phosphorylated, recruits signalling kinases and adaptor proteins to propagate downstream immune activation. The Gly→Ser change at position 307 is thought to alter the phosphorylation state of this domain, affecting how efficiently CD226 recruits downstream kinases such as Grb2 and PI3K after ligand engagement. The net biological result is that T cells and NK cells carrying the Ser307 form have a subtly altered activation threshold: they respond more readily to self-antigens presented in the context of inflammatory signals, a pattern consistent with the variant's observed association with multiple autoimmune diseases that have T-cell hyperactivation as a core pathogenic feature.

CD226 also functions in a competitive balance with its inhibitory counterpart TIGIT³³ TIGIT
A co-inhibitory receptor on T cells that competes with CD226 for the same ligands (PVR, Nectin-2) but delivers an inhibitory rather than activating signal, normally dampening T-cell responses. When CD226 Ser307 shifts this co-stimulatory/co-inhibitory balance toward excess activation, the self-tolerance mechanism that TIGIT normally enforces is partially undermined. This mechanistic framework positions CD226 Gly307Ser alongside PTPN22 R620W (rs2476601) and CTLA4 variants (rs3087243, rs231775) as part of a broader network of T-cell checkpoint genes where single-nucleotide variants collectively tune the threshold at which T cells become self-reactive.

The Evidence

The association was first established by Hafler et al. 2009⁴⁴ Hafler et al. 2009
Hafler JP et al. CD226 Gly307Ser association with multiple autoimmune diseases. Genes Immun, 2009 in a European cohort, finding association with type 1 diabetes (P=3.46×10⁻⁹), multiple sclerosis (P=4.2×10⁻⁴), and rheumatoid arthritis (P=0.017). This multi-disease pattern immediately distinguished CD226 from disease-specific loci and placed it among the small number of pan-autoimmune susceptibility genes.

The most comprehensive meta-analysis to date, Bai et al. 2020⁵⁵ Bai et al. 2020
Bai et al. Role of CD226 rs763361 polymorphism in susceptibility to multiple autoimmune diseases. Immunol Invest, 2020; 29 reports, 51 studies, 18,157 cases, 29,904 controls, confirmed associations with rheumatoid arthritis, SLE, type 1 diabetes, and multiple sclerosis, while finding no stable association for systemic sclerosis. An earlier meta-analysis by Song et al. 2012⁶⁶ meta-analysis by Song et al. 2012
Song G et al. Association between CD226 rs763361 polymorphism and susceptibility to autoimmune diseases. Lupus, 2012; 17 studies, 8,900 cases, 10,295 controls quantified the per-allele risk: overall OR 1.162 (95% CI 1.097–1.230, P<10⁻⁸), with T1D showing a stronger signal (OR 1.353) than SLE (OR 1.150) or systemic sclerosis (OR 1.126).

A further meta-analysis by Qiu et al. 2013⁷⁷ Qiu et al. 2013
Qiu et al. CD226 Gly307Ser association with multiple autoimmune diseases: a meta-analysis. Hum Immunol, 2013 reported an overall OR of 1.19 (95% CI 1.12–1.27), with markedly higher risk in South American populations (OR 1.72) compared to Asian (OR 1.46) or European (OR 1.29) cohorts. Both the additive and codominant models were significant, confirming that each T allele contributes incrementally to risk — TT homozygotes face approximately OR² ≈ 1.42× the risk of CT heterozygotes relative to CC homozygotes.

Sex-specific analysis in a Brazilian T1D cohort⁸⁸ Brazilian T1D cohort
Mattana et al. 2014, PMID 24891767 found the TT genotype association was driven primarily by females (P=0.0012), and TT carriers showed both higher GAD65 autoantibody frequency (31.9% vs 24.5%, OR 1.57) and lower residual C-peptide levels — markers of more aggressive autoimmune beta-cell destruction.

Practical Implications

The per-allele OR of ~1.16–1.19 places rs763361 among the moderate-effect common autoimmune risk variants. For CT heterozygotes (approximately 50% of the population), the risk elevation is modest and the absolute lifetime risk of any specific disease remains low. For TT homozygotes (~28%), the cumulative risk is more meaningful, especially when combined with other autoimmune susceptibility variants in the T-cell checkpoint network.

The diseases most consistently associated — T1D, MS, SLE, and RA — all have well-defined autoantibody and inflammatory biomarkers that can detect subclinical disease years before clinical presentation. Targeted monitoring is the primary actionable response.

Interactions

CD226 rs763361 is one node in a T-cell checkpoint network that includes PTPN22 rs2476601 (LYP phosphatase, reduces TCR signalling), CTLA4 rs3087243 and rs231775 (inhibitory co-receptor on T cells), and IL2RA rs2104286 (IL-2 receptor alpha, T-regulatory cell survival). These four genes converge on the decision of whether a T cell becomes activated or tolerized when encountering antigen. Carrying risk alleles at multiple checkpoints compounds the tolerance deficit multiplicatively. The compound effect of CD226 T + PTPN22 A (rs2476601) on T1D risk in particular deserves a compound action: both variants shift the T-cell activation threshold in the same direction (lower tolerance, higher autoimmune activation), and their combined risk is plausibly additive or super-additive.