Follicle-stimulating hormone (FSH) is the master regulator of follicle development in women and spermatogenesis in men. It is produced by the pituitary gland when the beta-subunit gene FSHB¹¹ beta-subunit gene FSHB
Located on chromosome 11p14.1, encoding the hormone-specific subunit that confers biological activity is transcribed and translated. rs11031006 is a G-to-A variant located approximately 26 kilobases upstream of the FSHB transcription start site, sitting within a conserved regulatory enhancer rather than the coding sequence of FSHB itself. This locus has emerged as one of the most robustly replicated genetic determinants of circulating FSH levels, female reproductive timing, dizygotic twinning propensity, and PCOS susceptibility — and it also influences male spermatogenic function through its effect on FSH production.

Note on variant identity: rs11031006 is a GWAS lead SNP at the FSHB locus; it is distinct from rs10835638, the -211G>T proximal promoter variant studied in many clinical male infertility trials. Both variants affect FSHB regulation, but through different mechanisms and at different distances from the gene. The two are not in strong linkage disequilibrium. Studies of male infertility citing "FSHB c.-211G>T" refer to rs10835638, while studies citing the 11p14.1 GWAS locus and twinning associations predominantly discuss rs11031006.

The Mechanism

The rs11031006 variant sits within a ~450 base-pair region that is highly conserved across placental mammals, a hallmark of functional regulatory elements. In vitro luciferase assays demonstrate that this region acts as a transcriptional enhancer of FSHB²² In vitro luciferase assays demonstrate that this region acts as a transcriptional enhancer of FSHB
The enhancer augments activin- and GnRH-stimulated FSHB transcription in gonadotrope cell models. The minor A allele creates a stronger binding site for Steroidogenic Factor 1 (SF1)³³ Steroidogenic Factor 1 (SF1)
A nuclear receptor transcription factor essential for gonadotrope cell identity and FSH gene expression, increasing enhancer activity approximately 1.5-fold compared to the major G allele in cell culture experiments.

This in vitro finding presents a mechanistic paradox: the A allele increases FSHB transcription experimentally, yet population data consistently show that individuals carrying the A allele have lower circulating FSH levels, higher LH/FSH ratios, and altered reproductive phenotypes compared to G-allele carriers. The discrepancy likely reflects the complexity of pituitary negative feedback regulation in vivo — the G allele may be associated with higher FSH in part because its carriers have faster hypothalamic-pituitary-gonadal axis dynamics overall. Mouse models partially reconcile this: female mice homozygous for the A-equivalent mutation show fewer litters and abnormal estrous cycling⁴⁴ fewer litters and abnormal estrous cycling
Despite no reduction in baseline FSH measured in the deletion model, the point mutation itself disrupts reproductive cycling, suggesting the locus affects reproductive cycling through mechanisms beyond steady-state FSH levels.

The Evidence

The most robust evidence comes from multiple GWAS. A 2016 study of mothers of spontaneous dizygotic twins (n~95,000 births in Iceland plus replication cohorts)⁵⁵ A 2016 study of mothers of spontaneous dizygotic twins (n~95,000 births in Iceland plus replication cohorts)
Mbarek et al., American Journal of Human Genetics identified rs11031006-G as a genome-wide significant twinning variant (p=1.54×10⁻⁹), with each copy of the G allele increasing the likelihood of a mother delivering fraternal twins by approximately 18% (OR 1.18 per copy). The same G allele was associated with higher serum FSH levels, earlier age at menarche, earlier age at first child, higher lifetime parity, lower PCOS risk, and earlier age at natural menopause — a constellation that collectively points to a more "fast-cycling" reproductive phenotype.

In polycystic ovary syndrome genetics, the 11p14.1 locus containing rs11031006 was identified in European PCOS GWAS as associated with altered LH levels and LH/FSH ratio. Women carrying copies of the A allele show higher LH/FSH ratios⁶⁶ Women carrying copies of the A allele show higher LH/FSH ratios
Consistent with the gonadotropin imbalance characteristic of PCOS, a pattern distinct from G-allele carriers who tend to have higher FSH relative to LH.

The male fertility significance was established in a 2022 GWAS of 760 idiopathic infertile men (validated in 1,140)⁷⁷ 2022 GWAS of 760 idiopathic infertile men (validated in 1,140)
Schubert et al., Journal of Clinical Endocrinology & Metabolism. The 11p14.1 locus (represented by rs11031005, in high LD with rs11031006) was the top genome-wide significant hit for serum FSH levels, explaining 4.65% of FSH variance overall and 6.95% of variance in the oligozoospermic subgroup specifically — a larger effect than the well-studied proximal promoter variant rs10835638 (which explains ~3.6% of FSH variance). Lower FSH in men impairs Sertoli cell function and reduces sperm production; the FSHB locus was identified as a potential etiologic factor in approximately 28% of men with idiopathic infertility.

Practical Implications

The practical implications of this variant differ by sex. In women, the A allele (associated with lower FSH and higher LH/FSH ratio) may contribute to longer menstrual cycles, slightly delayed folliculogenesis, and a reproductive axis phenotype that overlaps with some PCOS features — although rs11031006 alone is not diagnostic of PCOS. Women carrying the AA genotype (approximately 2% of European-ancestry individuals) may wish to discuss FSH and LH panel interpretation with a reproductive endocrinologist if experiencing irregular cycles, delayed conception, or unexplained subfertility.

In men, the A allele is associated with measurably lower FSH levels at the population level, which may impair spermatogenesis. Men with low-normal FSH and idiopathic infertility who carry variants at this locus represent a distinct etiologic subgroup⁸⁸ Men with low-normal FSH and idiopathic infertility who carry variants at this locus represent a distinct etiologic subgroup
Defined as functional secondary hypogonadism with isolated FSH deficiency, this group responds to exogenous FSH treatment with improved sperm parameters. Semen analysis combined with FSH measurement is the key initial investigation for male carriers, particularly for the AA homozygous genotype.

Interactions

rs10835638 (FSHB -211G>T, proximal promoter): This is a separate variant located 211 bp upstream of the FSHB mRNA transcription start site. Both rs11031006 and rs10835638 affect FSHB expression but at different positions and through distinct mechanisms. In men, rs10835638 T allele has been extensively studied and reduces FSH by ~0.51 IU/L per allele and testicular volume by ~3.2 ml. These two variants are not in strong LD and may have partially independent effects; their combined impact on FSH levels in men with idiopathic infertility is additive and warrants separate genotyping.

rs6166 (FSHR N680S): The FSH receptor sensitivity variant interacts functionally with FSHB variants. In men, the effect of FSHB locus variants on FSH-driven spermatogenesis is modulated by FSHR genotype — men with lower FSH production (FSHB A allele) and reduced FSH receptor sensitivity (FSHR GG) have a compounded spermatogenic disadvantage. This interaction has been documented for the proximal FSHB variant, and the same pathway logic applies to rs11031006. A compound action may be warranted when both unfavorable genotypes co-occur.

RORA (RAR-related Orphan Receptor Alpha) occupies a dual role in human biology: it is a core activator of the circadian clock¹¹ circadian clock
RORA binds ROR-response elements in the BMAL1 promoter, driving its transcription and sustaining the ~24-hour oscillation in the suprachiasmatic nucleus and a key regulator of immune tolerance in the airways. The rs11071559 variant sits in the first intron of the RORA-1 transcript on chromosome 15q22.2 — a regulatory region, not a coding sequence — which means it likely influences how much RORA protein is produced rather than altering the protein itself. The best-replicated association for this specific variant is with asthma susceptibility: the common C allele carries slightly elevated risk, while the rarer T allele is protective.

The Mechanism

As an intronic variant with no known amino acid change, rs11071559 is thought to act as a regulatory tag SNP²² regulatory tag SNP
A tag SNP marks a block of correlated variants; rs11071559 may be in linkage disequilibrium with a causal variant that alters a transcription factor binding site or splicing regulatory element in RORA intron 1. RORA itself regulates the promoter activity of NPSR1 (neuropeptide S receptor 1), a gene independently associated with asthma and panic disorder. When RORA expression is altered, NPSR1 promoter activity changes reciprocally³³ NPSR1 promoter activity changes reciprocally
Overexpression of RORA in cell models decreased NPSR1 promoter activity, suggesting a negative regulatory loop. NPSR1 signaling in turn activates a pathway that includes circadian clock genes, creating a feedback loop linking the circadian system to airway inflammation. This may partly explain why asthma symptoms — particularly nocturnal asthma — follow a circadian pattern⁴⁴ circadian pattern
Airway inflammation and bronchoconstriction peak between 2–4 AM in susceptible individuals, mirroring diurnal rhythms in cortisol, melatonin, and mast cell activity.

The Evidence

The asthma association for rs11071559 is among the better-replicated findings in the RORA locus. The APCAT consortium meta-analysis⁵⁵ APCAT consortium meta-analysis
Ramasamy et al. PLOS One 2012 brought this variant to genome-wide significance (p = 2.4 × 10⁻⁹) across six European population-based cohorts totaling approximately 5,751 asthmatics and 28,139 controls. Within the BAMSE and PARSIFAL cohorts, Acevedo et al. 2013⁶⁶ Acevedo et al. 2013
PLOS One found the T allele protective against physician-diagnosed childhood asthma with an odds ratio of 0.71 (95% CI: 0.55–0.92, p = 0.007) in the combined dataset, a modest but replicable effect. Importantly, the RORA association was primarily with asthma diagnosis rather than atopic traits (IgE levels, eczema), suggesting RORA's role in airway inflammation specifically rather than general atopic tendency.

A note on sleep and hormonal phenotypes: RORA is a validated circadian clock gene and BMAL1 activator, and other RORA variants (notably rs75981965) have been associated with sleep duration in a Taiwanese biobank study⁷⁷ Taiwanese biobank study
10,112 subjects, p = 9.93 × 10⁻⁵ after Bonferroni correction. However, rs11071559 itself has not been associated with sleep duration, chronotype, or hormonal traits in any published GWAS. The association for this specific SNP is asthma only. The batch placed this variant in the hormones-sleep category because RORA is a circadian gene, but the actionable evidence for rs11071559 is in airway immunity. Users with the CC genotype should focus on the respiratory health implications described below.

Practical Actions

For the common C-allele carriers (CC or CT), the slightly elevated asthma risk from the RORA locus is modest and modifiable. The circadian-immune link means that circadian disruption amplifies airway inflammation⁸⁸ circadian disruption amplifies airway inflammation
Shift workers and those with social jet lag have higher rates of asthma exacerbation and more severe symptom patterns. Maintaining consistent sleep-wake timing is therefore a genotype-relevant action — not as a generic sleep hygiene recommendation, but because circadian misalignment specifically impairs RORA function in airway epithelial and immune cells, which is the pathway through which this variant acts. Avoiding known asthma triggers and monitoring for nocturnal symptoms are warranted.

Interactions

RORA interacts genetically and biologically with NPSR1 (neuropeptide S receptor 1 gene, chromosome 7p15). The interaction is strongest for nocturnal asthma: rs7164773 in RORA and NPSR1 variants jointly predict asthma phenotype better than either alone. If you have both RORA and NPSR1 risk genotypes and experience primarily night-time or early-morning asthma symptoms, the RORA-NPSR1 interaction may be the relevant mechanism.

The circadian clock gene REV-ERBα (NR1D1) is RORA's functional antagonist at the BMAL1 promoter — both bind the same response element, with opposing effects. Variants in NR1D1 that reduce its repressive activity can partially compensate for reduced RORA activity, and vice versa. Pathway-level interactions between RORA, NR1D1, and RORB on sleep and circadian phenotypes have been reported but not yet linked to this specific rsid.

The ataxia telangiectasia mutated (ATM) gene encodes one of the master regulators of the cellular DNA damage response — a kinase that springs into action when double-strand DNA breaks are detected, coordinating cell cycle arrest, DNA repair, and, when damage is irreparable, programmed cell death. What makes ATM central to longevity biology is what it does downstream: ATM activates AMPK¹¹ AMPK
AMP-activated protein kinase — the cell's master energy sensor, activated when ADP/AMP ratios rise, and the primary target of metformin's longevity effects, which in turn inhibits mTOR and triggers autophagy. This positions ATM not just as a DNA repair enzyme but as a node connecting genomic stability to the metabolic hallmarks of aging.

The rs11212617 variant sits in an intron within the ATM locus on chromosome 11q22, inside a 340 kb linkage disequilibrium block. The C allele is associated with enhanced metformin response, suggesting it may modify how efficiently ATM activates AMPK under metabolic stress — the same pathway that metformin engages when it inhibits mitochondrial complex I and raises the AMP:ATP ratio.

The Mechanism

Metformin's primary action is inhibition of mitochondrial complex I in the electron transport chain, raising the cellular AMP:ATP ratio and directly activating AMPK. But metformin also appears to engage the ATM-AMPK axis: experimental work using ATM inhibitor KU-55933 in rat hepatoma cells²² experimental work using ATM inhibitor KU-55933 in rat hepatoma cells
Zhou K et al. Common variants near ATM are associated with glycemic response to metformin in type 2 diabetes. Nature Genetics, 2011 showed attenuated AMPK phosphorylation in response to metformin when ATM was blocked. The mechanistic interpretation was complicated by a subsequent finding that KU-55933 also inhibits OCT1 (a metformin transporter), potentially reducing intracellular metformin rather than directly blocking an ATM-AMPK signal. Whether ATM acts directly on AMPK, indirectly through modulating DNA damage-sensing cascades that converge on AMPK, or primarily via effects on metformin transport remains an active question.

What is established is that activated ATM phosphorylates and stabilizes SIRT6, and that ATM activity also restrains mTORC1 via the AMPK-TSC2 pathway after DNA damage. Boosting ATM activity extended lifespan in mouse models of progeria³³ Boosting ATM activity extended lifespan in mouse models of progeria
Qian M et al. Boosting ATM activity alleviates aging and extends lifespan in a mouse model of progeria. eLife, 2018, and ATM-deficient mice show accelerated metabolic dysfunction and premature aging — establishing the gene's role in organismal longevity beyond its classic function in cancer suppression.

The Evidence

The landmark GWAS by Zhou et al. 2011⁴⁴ The landmark GWAS by Zhou et al. 2011
Common variants near ATM are associated with glycemic response to metformin in type 2 diabetes. Nature Genetics discovered rs11212617 by scanning 1,024 Scottish type 2 diabetes patients on metformin. In the combined meta-analysis of 3,920 patients, the C allele reached genome-wide significance (P=2.9×10⁻⁹) for metformin treatment success (achieving HbA1c below 7%), with an odds ratio of 1.35 (95% CI 1.22–1.49). Each additional C allele correlated with 0.11% lower HbA1c on metformin treatment (P=6.6×10⁻⁷).

Van Leeuwen et al. 2012⁵⁵ Van Leeuwen et al. 2012
A gene variant near ATM is significantly associated with metformin treatment response in type 2 diabetes: a replication and meta-analysis of five cohorts. Diabetologia replicated the association across three new cohorts and declared rs11212617 the first robustly replicated common pharmacogenetic variant for metformin (combined five-cohort OR 1.25, P=7.8×10⁻⁶).

However, Florez et al. 2012⁶⁶ Florez et al. 2012
The C allele of ATM rs11212617 does not associate with metformin response in the Diabetes Prevention Program. Diabetes Care failed to confirm the association in 2,994 participants treated with metformin for diabetes prevention (HR 1.17, P=0.13). The authors note an important distinction: the DPP enrolled prediabetic individuals while the discovery GWAS enrolled established type 2 diabetics — the biological context of metformin's action may differ substantially between these populations. Multiple other studies in European, South Asian, and East Asian populations have found inconsistent replication, suggesting the association may be population-specific or context-dependent.

In an independent direction, Cuyàs et al. 2019⁷⁷ Cuyàs et al. 2019
METTEN trial — Frontiers in Oncology found that C allele carriers among HER2-positive breast cancer patients had a 7.94-fold higher probability of pathological complete response when treated with neoadjuvant metformin (p=0.011), while no association existed in the control arm — extending the variant's pharmacogenetic relevance to cancer treatment.

The mixed replication record is reflected in the moderate evidence level: the initial GWAS signal is robust within European type 2 diabetes cohorts, but context-dependence, population variation, and an unresolved molecular mechanism prevent elevation to strong.

Practical Actions

The core implication of rs11212617 is pharmacogenetic: C allele carriers — both AC heterozygotes and CC homozygotes — appear to show improved glycemic response to metformin in the context of established type 2 diabetes. For individuals with AA genotype who are prescribed metformin, awareness that this variant may confer reduced metformin efficacy is worth discussing with a clinician, particularly if glycemic targets are not met at standard doses.

Beyond pharmacogenetics, the ATM pathway connects to the same AMPK-mTOR axis that underlies caloric restriction and intermittent fasting biology. The TAME trial (Targeting Aging with Metformin)⁸⁸ TAME trial (Targeting Aging with Metformin) is currently testing whether metformin can extend healthspan in non-diabetic older adults specifically through these longevity pathways. Individual rs11212617 genotype may ultimately predict differential benefit even in that prevention context.

Interactions

rs11212617 operates within the same AMPK-mTOR axis as rs2295080 (MTOR promoter variant). Individuals carrying both the ATM rs11212617 A allele (reduced ATM-AMPK coupling) and MTOR rs2295080 TT genotype (highest mTOR expression) would face a compound disadvantage: reduced capacity to activate AMPK through the ATM route combined with elevated constitutive mTOR activity. This biologically plausible interaction has not been formally tested in a published combined-genotype study. rs2802292 (FOXO3) is also a longevity pathway partner: FOXO3 activity is downstream of AMPK and upstream mTOR signaling, making all three variants part of the same regulatory circuit linking DNA damage response to cellular aging.

Your brain constantly decides how strongly to respond to its own growth signals. Brain-derived neurotrophic factor (BDNF)¹¹ Brain-derived neurotrophic factor (BDNF)
BDNF is a protein that promotes the survival, growth, and differentiation of neurons and synapses. It is one of the most important regulators of synaptic plasticity, learning, and emotional resilience. Low BDNF has been repeatedly linked to depression, anxiety, and cognitive decline is among the most potent of these signals — yet unbridled BDNF activity would be disruptive. The SORCS3 gene encodes a molecular gatekeeper that sorts the BDNF receptor TrkB between active surface positions and intracellular compartments, calibrating how much BDNF signaling actually reaches the synapse. Genetic variation at this locus has now emerged as one of the more robust molecular links between synaptic neurotrophin trafficking and psychiatric risk in the human genome.

The Mechanism

SORCS3 belongs to the VPS10 domain-containing receptor family²² VPS10 domain-containing receptor family
A family of type-I transmembrane proteins that act as intracellular sorting receptors, routing cargo proteins between cellular compartments. The family includes SORTILIN, SORCS1, SORCS2, and SORCS3 — all critical for neuronal protein trafficking. The VPS10 (vacuolar protein sorting 10) domain is the cargo-binding domain first characterized in yeast. In neurons, SORCS3 and its close relative SORCS1 function as [intracellular trafficking receptors for TrkB | Tropomyosin-related kinase B (TrkB) is the primary high-affinity receptor for BDNF. When BDNF binds TrkB at the cell surface, it triggers downstream signaling cascades that support synaptic plasticity, neuronal survival, and mood regulation], routing the receptor away from the synaptic surface and attenuating BDNF signaling.

A 2018 mouse knockout study³³ 2018 mouse knockout study
Subkhangulova A et al. SORCS1 and SORCS3 control energy balance and orexigenic peptide production. EMBO Rep, 2018 showed that neurons lacking both SORCS1 and SORCS3 display elevated TrkB phosphorylation following BDNF application — confirming that these receptors normally brake TrkB activation. The rs11599236 variant sits within an intron of SORCS3, in a position likely to influence gene expression rather than protein structure. Critically, a systematic analysis of 46 independent SORCS3 SNPs found that alleles linked to better psychiatric outcomes were consistently associated with higher SORCS3 expression⁴⁴ alleles linked to better psychiatric outcomes were consistently associated with higher SORCS3 expression
Kamran M et al. Independent Associated SNPs at SORCS3 and Its Protein Interactors for Multiple Brain-Related Disorders and Traits. Genes (Basel), 2023. This direction of effect — more SORCS3 protein equals better outcomes — points to a loss-of-function mechanism, where variants that reduce SORCS3 expression allow TrkB to accumulate at inappropriate synaptic locations, dysregulating the pro-neurotrophin / mature-neurotrophin balance at the synapse.

Beyond BDNF trafficking, SORCS3 interacts with p75NTR (the pan-neurotrophin receptor), promoting its internalization and lysosomal degradation. p75NTR is the primary receptor for the immature precursor form of BDNF (proBDNF), which has opposite effects to mature BDNF — promoting apoptosis and long-term depression at synapses rather than survival and potentiation. Dysregulated p75NTR activity is implicated in synaptic weakening associated with chronic stress and depression.

The Evidence

The association between rs11599236 and mood-related traits is among the most replicated SORCS3 findings in the literature. In the GWAS Catalog, this single variant is associated with at least six psychiatric and affective phenotypes spanning independent large-scale studies:

Mood instability and negative affect. The C allele at rs11599236 is associated with reduced subjective wellbeing (beta = -0.0074, p = 1×10⁻¹⁵, N > 300,000) and mood instability (p = 5×10⁻¹⁰) in population-based GWAS. The T allele is associated with increased ratings of feeling miserable (z-score beta = 7.05, p = 2×10⁻¹²) and higher neuroticism scores (p = 7×10⁻⁸ to p = 3×10⁻¹⁹ across independent cohorts).

Psychiatric cross-disorder risk. In the Psychiatric Genomics Consortium cross-disorder analysis spanning ADHD, autism spectrum disorder, bipolar disorder, major depressive disorder, and schizophrenia simultaneously, rs11599236 emerged as genome-wide significant (z-score beta = 6.19, p = 6×10⁻¹⁰), highlighting SORCS3 as a pleiotropic risk locus shared across neuropsychiatric conditions.

ADHD. A 2023 mega-GWAS⁵⁵ 2023 mega-GWAS
Demontis D et al. Genome-wide analyses of ADHD identify 27 risk loci, refine the genetic architecture and implicate several cognitive domains. Nat Genet, 2023 of 38,691 ADHD cases and 186,843 controls identified SORCS3 as one of 27 genome-wide significant loci, with convergent evidence from both common variants and an elevated burden of rare protein-truncating variants in SORCS3 among ADHD cases.

Autism. Tag SNPs within SORCS3 (rs9787523 and rs3750261) showed nominal association with autism in a Han Chinese replication study⁶⁶ Han Chinese replication study
Heliyon 2024 of 757 trios, and the gene has appeared in multiple independent autism GWAS datasets.

Neurodegeneration. SORCS3 expression is significantly reduced in Alzheimer's disease brain tissue compared to controls (p = 5.1×10⁻⁵⁷⁷ p = 5.1×10⁻⁵
Reitz C et al. Independent and epistatic effects of variants in VPS10-d receptors on Alzheimer disease risk and processing of the amyloid precursor protein (APP). Transl Psychiatry, 2013), and SORCS3 knockdown increases amyloid precursor protein (APP) processing by threefold — a direct functional link to the amyloid cascade. Non-coding SORCS3 variants were associated with dementia in women in a whole-genome sequencing study of the Women's Health Initiative cohort.

Practical Actions

For carriers of one or two C alleles, the actionable insight centers on supporting BDNF-TrkB signaling and maintaining synaptic neurotrophin balance. The evidence base for specific interventions is mostly mechanistic and derived, but the biological logic is clear: if SORCS3 hypofunction allows pro-neurotrophin signaling to dominate, strategies that boost mature BDNF availability and reduce chronic stress-mediated BDNF depletion are rational targets.

The single most evidence-backed modulator of BDNF expression is aerobic exercise — specifically continuous, moderate-intensity effort that reliably increases serum BDNF and hippocampal BDNF mRNA. This is distinct from generic exercise advice: the target is BDNF amplification in the context of a receptor that has reduced trafficking fidelity. Sleep quality independently regulates BDNF protein levels, with sleep deprivation robustly suppressing BDNF expression. Omega-3 DHA is incorporated into neuronal membranes and supports TrkB receptor signaling efficiency. Magnesium deficiency is associated with dysregulated neurotrophin signaling and increased anxiety-like behavior in animal models.

Interactions

SORCS3 sits in a gene family alongside SORCS1 and SORCS2, which interact with overlapping but distinct cargo sets. Variants in SORCS1 (rs10790256, rs600879) and SORCS2 have been independently associated with Alzheimer's risk and psychiatric phenotypes — raising the possibility that multiple VPS10-receptor variants compound their effects on neurotrophin trafficking in the same individual. SORCS3 gene-sets are enriched for synaptic biology⁸⁸ synaptic biology
Specifically, 15 biological process terms related to synapse organization, postsynaptic processes, and synaptic signaling were over-represented in SORCS3-associated genes in the Kamran et al. 2023 analysis, overlapping substantially with the biological pathways implicated by schizophrenia and bipolar disorder GWAS — suggesting that rs11599236 may act additively with variants in DRD2 (rs1800497), COMT (rs4680), and BDNF (rs6265) that alter the same dopamine-neurotrophin signaling axis.

Your blood pressure is governed in part by a family of hormones called natriuretic peptides¹¹ natriuretic peptides
ANP (atrial natriuretic peptide) and BNP (brain natriuretic peptide) — released by heart muscle cells in response to stretch when blood volume rises. These peptides act on blood vessel walls and kidneys to promote vasodilation and sodium excretion, pulling blood pressure back down. But natriuretic peptides need a mechanism to be cleared from the bloodstream — otherwise their blood-pressure-lowering signal would never switch off.

That clearance role belongs to NPR3, the natriuretic peptide receptor C gene. NPR3 encodes a decoy receptor on vascular smooth muscle and endothelial cells that binds ANP and BNP and internalizes them for degradation — functioning as a physiological "off switch" for natriuretic peptide signaling. When NPR3 is highly expressed, natriuretic peptides are cleared rapidly, reducing their protective vasodilatory signal. When NPR3 is lower, peptides linger in circulation longer, sustaining blood pressure reduction.

The variant rs1173771 sits near the NPR3 transcription start site²² transcription start site
the region of DNA from which gene expression is initiated and reached genome-wide significance for systolic blood pressure, diastolic blood pressure, pulse pressure, and hypertension in the ICBP meta-analysis³³ ICBP meta-analysis
International Consortium for Blood Pressure, combining data from ~200,000 individuals of European descent. The G allele is the BP-elevating allele.

The Mechanism

Functional studies show that the BP-elevating haplotype tagged by rs1173771 is associated with lower endogenous NPR3 mRNA and protein in vascular smooth muscle cells (VSMCs)⁴⁴ lower endogenous NPR3 mRNA and protein in vascular smooth muscle cells (VSMCs)
detected in primary human VSMCs from carriers of the BP-elevating allele vs. the protective allele, along with reduced open chromatin and nuclear protein binding at the locus. In other words, the G allele reduces NPR3 expression — but paradoxically, lower NPR3 means less clearance of vasodilatory natriuretic peptides. The key is that the downstream consequences go beyond simple clearance.

Lower NPR3 expression in VSMCs produces three compounding effects: increased cell proliferation⁵⁵ increased cell proliferation
partially reversible by natriuretic peptide treatment, linking genotype to atherogenic smooth muscle behavior, enhanced angiotensin II-induced intracellular calcium flux⁶⁶ angiotensin II-induced intracellular calcium flux
a molecular signature of heightened vasoconstriction responsiveness, and increased VSMC contraction. The net result is a vascular phenotype primed for higher resting blood pressure and greater pressor responses.

A 2024 experimental model using Dahl salt-sensitive rats⁷⁷ 2024 experimental model using Dahl salt-sensitive rats
Dahl SS rats are a classic model of salt-sensitive hypertension reinforced this picture: when researchers deleted the orthologous noncoding haplotype region, NPR3 was upregulated in arteries and salt-induced systolic blood pressure rose by approximately 10 mmHg less than in unedited controls. The 2025 renal chromatin study⁸⁸ 2025 renal chromatin study
Mapping of ATAC-seq accessible chromatin in human proximal tubule and medullary thick ascending limb segments found rs1173771 enriched within regulatory chromatin in kidney segments that control sodium-water balance, extending the relevant tissue context from vessels to kidney.

The Evidence

The core blood pressure association is established and multiply replicated⁹⁹ established and multiply replicated
GWAS Catalog accessions GCST001227, GCST001228, GCST001235, GCST001236, GCST001238. In the original ICBP meta-analysis, the G allele reached P=2×10⁻¹⁶ for systolic BP (beta +0.50 mmHg per allele) and P=9×10⁻¹² for diastolic BP (beta +0.26 mmHg), P=3×10⁻¹⁰ for hypertension (OR 1.06 per allele).

The effect is small per allele but has been replicated transethnic across European and East Asian populations. A 2013 transethnic meta-analysis¹⁰¹⁰ 2013 transethnic meta-analysis
East Asian stage 1 n=26,600, stage 2 up to 28,783 participants confirmed rs1173771 shows suggestive replication (P=0.018) for mean arterial pressure with consistent effect direction across ancestries. The Genetics of Postural Hemodynamics (GPH) Consortium study¹¹¹¹ Genetics of Postural Hemodynamics (GPH) Consortium study
Published Eur Heart J 2012 found that the A allele was modestly protective against orthostatic hypotension (OR 0.92, 95% CI 0.87–0.98, P=0.009), consistent with the G allele being the BP-elevating variant.

The functional work by Ren et al. (2018) bridges the statistical association to mechanism: using primary human VSMCs genotyped at rs1173771, they demonstrated allele-specific differences in NPR3 expression and downstream vascular cell behavior, elevating this from a statistical GWAS hit to a biologically understood variant.

Practical Actions

The per-allele blood pressure effect (~0.5 mmHg SBP per G allele) is modest in isolation. GG homozygotes carry the full additive effect (~1 mmHg SBP higher than AA), which over a lifetime can compound with other risk factors. The critical question is whether blood pressure lands in the normal, elevated, or hypertensive range in the context of your full genetic and lifestyle profile.

The most important practical implication is monitoring: home blood pressure measurement provides personalized data that genetic risk estimates cannot. If blood pressure trends upward — even into the high-normal range (120–129/80 mmHg) — early lifestyle and dietary intervention is more effective than waiting for overt hypertension to develop.

Dietary sodium is particularly relevant: the NPR3 pathway is central to sodium excretion and fluid balance. Excess dietary sodium blunts natriuretic peptide effectiveness, compounding the reduced clearance capacity created by lower NPR3 expression.

Interactions

rs1173771 is part of a 17.4 kb haplotype block at the NPR3 locus containing 11 SNPs. The sentinel variant rs1173771 tags this haplotype, and the entire block functions as a regulatory unit. Other SNPs in the NPR3 gene region — including rs2270915¹²¹² rs2270915
a previously studied NPR3 promoter polymorphism associated with blood pressure in type 2 diabetes — may have partially overlapping or independent effects. Other blood-pressure GWAS loci in the natriuretic peptide pathway include NPPA/NPPB (ANP and BNP precursor genes) and NPPB itself; collectively, variation across the natriuretic peptide system accumulates to influence resting blood pressure trajectory.

Most genetic diseases work by a familiar logic: if you change a gene's sequence, you change its protein. But the PYGM Trp798Arg mutation (also written W798R or c.2392T>C) tells a stranger story. Patients who carry this variant produce normal amounts of PYGM messenger RNA — the gene is being read and copied faithfully — yet their muscle cells contain virtually no myophosphorylase¹¹ myophosphorylase
muscle-specific glycogen phosphorylase, the enzyme encoded by PYGM that cleaves glucose-1-phosphate from glycogen at the start of exercise protein. The mutation seems to hijack a post-translational quality-control system, condemning the newly made protein to rapid degradation before it can fold and function. The result is a complete block in muscle glycogenolysis indistinguishable from a frameshift or nonsense variant.

This is McArdle disease²² McArdle disease
glycogen storage disease type V (GSD-V), an autosomal recessive metabolic myopathy first described by Brian McArdle in 1951 and caused by loss of functional myophosphorylase in skeletal muscle: exercise intolerance, painful cramps in the first minutes of exertion, and the characteristic "second wind" as blood-glucose and fatty-acid delivery eventually compensate for the absent glycogenolytic pathway.

ClinVar classifies Trp798Arg as Pathogenic³³ ClinVar classifies Trp798Arg as Pathogenic
VCV000002312, criteria provided, multiple submitters, no conflicts for glycogen storage disease type V, with eleven contributing submissions.

The Mechanism

Tryptophan 798 sits in the C-terminal catalytic domain of myophosphorylase, near the region involved in dimer interface contacts and allosteric effector binding. The substitution of tryptophan (neutral, aromatic, bulky) with arginine (positively charged, flexible) creates a physicochemical mismatch that the cell's protein quality control system apparently cannot tolerate.

The critical evidence comes from an iPSC-based skeletal muscle model developed by Lucía et al.⁴⁴ Lucía et al.
Lucía A et al. Creation of an iPSC-based skeletal muscle model of McArdle disease harbouring the mutation c.2392T>C (p.Trp798Arg) in the PYGM gene. Biomedicines, 2023. Patient-derived myofibers showed: (1) PYGM mRNA expression statistically indistinguishable from healthy controls; (2) complete absence of myophosphorylase protein by Western blot and immunofluorescence; and (3) PAS-positive glycogen accumulation around cell nuclei — the morphological signature of McArdle disease. The mRNA-protein discordance points to a post-translational mechanism: the mutant protein is produced but immediately recognized as misfolded and degraded, likely by the proteasomal or chaperone-mediated quality-control machinery.

This "unexpected consequence" of a missense mutation has been observed broadly across PYGM pathogenic missense variants — many of which cause complete myophosphorylase absence despite their theoretical missense status. The clinical implication is that Trp798Arg homozygotes and compound heterozygotes have the same functional deficit as patients carrying nonsense or frameshift mutations: zero functional enzyme.

The Evidence

The variant was identified as a recurrent allele in the Spanish McArdle population by Rubio et al.⁵⁵ Rubio et al.
Rubio JC et al. A proposed molecular diagnostic flowchart for myophosphorylase deficiency (McArdle disease) in blood samples from Spanish patients. Hum Mutat, 2007, who found it in 9 of 55 Spanish patients — making it the third most frequent PYGM mutation in Spain alongside p.R50X and p.Gly205Ser, and characterizing it as "a virtually Spanish-private mutation" not identified in British, American, German, French, or Italian cohorts.

A 2015 comprehensive update of all documented PYGM mutations by Nogales-Gadea et al.⁶⁶ Nogales-Gadea et al.
Nogales-Gadea G et al. McArdle disease: update of reported mutations and polymorphisms in the PYGM gene. Hum Mutat, 2015 catalogued 147 pathogenic mutations. Trp798Arg is documented as the second most frequent pathogenic missense allele in the Spanish population, with no genotype-phenotype correlation found across PYGM variants — the specific mutation does not predict disease severity.

Vieitez et al.⁷⁷ Vieitez et al.
Vieitez I et al. Molecular and clinical study of McArdle's disease in a cohort of 123 European patients. Neuromuscul Disord, 2011 confirmed the dominance of p.R50X (61.7% allelic frequency) across European populations and the population-specific distribution of secondary alleles, with no cross-population case of Trp798Arg in non-Spanish Europeans.

McArdle disease affects approximately 1 in 100,000–167,000 people globally. Life expectancy is normal; about 11% of patients develop permanent proximal weakness after age 40, and rhabdomyolysis episodes carry acute kidney injury risk.

Practical Actions

The management of McArdle disease is identical regardless of which specific PYGM pathogenic mutation is present — phenotype tracks with enzyme absence, not genotype. The most validated, genotype-specific interventions are:

Pre-exercise sucrose: 25–40 g of sucrose 5 minutes before exercise significantly improves tolerance by raising blood glucose before glycolysis demand peaks. This is the most validated nutritional intervention specific to myophosphorylase deficiency — it bypasses the glycogen block by providing circulating glucose directly.

Second-wind strategy: Starting exercise at very low intensity for 8–10 minutes allows blood-glucose delivery and fatty acid mobilization to compensate for the absent glycogenolytic pathway before effort is increased. Patients who exploit this phenomenon achieve substantially better exercise capacity.

Myoglobinuria recognition: Dark urine after exertion signals rhabdomyolysis and requires immediate rest, aggressive hydration, and emergency evaluation for acute kidney injury.

Interactions

This is an autosomal recessive variant. Compound heterozygosity — inheriting Trp798Arg on one chromosome and any other pathogenic PYGM allele (R50X, Gly205Ser, or any of the ~147 documented variants) on the other — produces full McArdle disease with zero functional myophosphorylase. The most common pairings in Spanish patients involve compound heterozygosity between Trp798Arg and p.R50X (rs116855232).

The "virtually Spanish-private" epidemiology of Trp798Arg means this allele is of greatest clinical relevance in patients of Iberian or Latin American descent. Genetic panels for McArdle disease in Spanish-ancestry individuals should explicitly include this variant alongside R50X and Gly205Ser.

Your serum uric acid level is regulated in large part by transporters in your kidneys, and the strongest genetic region for that regulation sits in the SLC2A9 gene. Most people have heard of the coding variant rs3733591 (Arg265His), but the SLC2A9 locus contains a second, genetically independent signal — rs11942223 — that contributes its own effect on urate transport. Understanding which variant you carry matters because their effects are additive, and the intronic signal captured by rs11942223 is the one with the most pronounced sex-specific effect.

SLC2A9 encodes GLUT9¹¹ GLUT9
Glucose Transporter 9, a high-capacity urate transporter in the kidney proximal tubule that mediates urate reabsorption; despite its name, it transports urate far more efficiently than glucose. The rs11942223 variant lies within an intron and does not change the protein directly, but it tags a regulatory or structural haplotype that modulates the efficiency of renal urate clearance. This intronic signal is in strong LD with rs12498742, another well-characterised urate-GWAS proxy for the same haplotype block.

The Mechanism

rs11942223 is an intron variant that does not alter the amino acid sequence of GLUT9. Its effect on urate transport is regulatory in nature — it tags a haplotype block within the SLC2A9 gene that influences transporter expression or splicing. The variant explains 1.2% of variance in serum urate in men and up to 6% in women (Döring et al., 2008)²² (Döring et al., 2008), a difference attributed to an interaction with estrogen signalling: estrogen independently stimulates renal urate excretion, which amplifies the apparent genetic effect in women compared to men.

The rs11942223 T allele (major allele, ~74–78% frequency in Europeans) is associated with reduced renal urate excretion and higher steady-state serum uric acid. The C allele (minor allele, ~22–26% in Europeans, rare in Polynesian populations) is protective — each copy reducing serum urate by approximately 0.23–0.36 mg/dL in men and 0.36–0.46 mg/dL in women.

Critically, this signal is independent of rs3733591 (Arg265His). Linkage disequilibrium between rs11942223 and rs3733591 is very low (r² = 0.03–0.05 across Māori, Eastern Polynesian, Western Polynesian, and Caucasian populations), confirming they tag distinct biological mechanisms at the same locus. Carriers of T alleles at both variants carry additive genetic risk for elevated uric acid.

The Evidence

Discovery of the SLC2A9 locus with sex-specific effects: Döring et al. (2008) mapped intronic SLC2A9 variants explaining up to 6% of urate variance in women but only 1.2% in men, making this the first demonstration of a major sex-specific quantitative trait locus for serum urate. The study used KORA German cohort data with independent replication (Döring et al., 2008)³³ (Döring et al., 2008).

Multi-ethnic gout association: Hollis-Moffatt et al. (2009) genotyped rs11942223 alongside rs3733591 and two other SLC2A9 variants in Māori, Pacific Island, and Caucasian gout case-control sets. The major (T) allele of rs11942223 was associated with gout across all populations (P = 1.6 × 10⁻⁶ for the combined sample), with odds ratio > 2.0 in Māori and Pacific Island subjects. Crucially, LD between rs11942223 and rs3733591 was r² = 0.03–0.05 — confirming these are independent loci acting in parallel (Hollis-Moffatt et al., 2009)⁴⁴ (Hollis-Moffatt et al., 2009).

Gene-environment interaction with dietary fructose: A fructose-challenge study in 76 volunteers (25 Māori, 26 Pacific, 25 European Caucasian) found that the C allele (protective) was associated with an attenuated hyperuricemic response and increased fractional excretion of uric acid after a fructose load, but only in the Caucasian subgroup — not in Māori or Pacific participants (Batt et al., 2013)⁵⁵ (Batt et al., 2013). This suggests the protective haplotype tagged by the C allele has functionally different penetrance across ancestries.

Sugar-sweetened beverages reverse the protective effect: A follow-up study found that the normally urate-lowering C allele reversed its effect under regular sugar-sweetened beverage (SSB) exposure — C allele carriers showed a 15% increase in gout risk per daily SSB serving, compared to 12% in non-carriers (Dalbeth et al., 2014)⁶⁶ (Dalbeth et al., 2014). This gene-environment interaction is clinically actionable: the C allele's urate-lowering benefit disappears under high fructose load, making SSB avoidance particularly important for carriers of either allele.

Practical Actions

The rs11942223 T allele (major allele) elevates serum urate by reducing renal clearance, and the effect is amplified in women (especially post-menopausal) and by dietary fructose. Since this is an independent signal from rs3733591, carrying T alleles at both variants compounds the urate burden — each locus adds independently to the genetic elevation.

The most effective dietary intervention for T/T homozygotes is eliminating sugar-sweetened beverages (sodas, fruit juices, energy drinks), since fructose both raises urate directly and — for this specific variant — can override the dose-response that would normally moderate the effect. Purine restriction (organ meats, shellfish, beer) remains important. Low-fat dairy, coffee, vitamin C, and adequate hydration all contribute to modest but evidence-supported reductions in serum urate.

Given the pronounced sex-specific effect (~5× larger variance explained in women), women who carry T alleles should have serum uric acid checked at menopause, when the loss of estrogen's uricosuric protection unmasks genetic risk that was previously buffered.

Interactions

Independent from rs3733591 (SLC2A9 Arg265His): The r² between rs11942223 and rs3733591 is only 0.03–0.05 — these variants are effectively unlinked and tag distinct regulatory and coding effects in the same gene. Carrying risk alleles at both compounds the urate elevation additively. Individuals with TT at rs11942223 and CC at rs3733591 carry two independent SLC2A9 risk signals.

SLC2A9 and ABCG2 (rs2231142): ABCG2 Q141K reduces intestinal urate secretion through a completely different pathway. Carrying risk alleles at rs11942223 (renal reabsorption) and ABCG2 rs2231142 (intestinal secretion) is additive and can produce mean serum urate above 7 mg/dL even in healthy individuals without dietary provocation.

Fructose and sugar-sweetened beverages: A documented gene-environment interaction reverses the C allele's protective effect under SSB exposure. This is unique among known urate GWAS variants and means that SSB avoidance is the single most important lifestyle action for this specific genotype.

Sex and menopausal status: The intronic SLC2A9 signal has a substantially larger effect in women (6% of urate variance vs. 1.2% in men), mediated through estrogen's interaction with renal urate handling. Post-menopausal women carrying T alleles lose this hormonal buffer and represent the highest-risk subgroup for this variant.

The heart is among the most metabolically active organs in the body, and triiodothyronine (T3)¹¹ triiodothyronine (T3)
the active form of thyroid hormone that directly enters cardiomyocyte nuclei and regulates gene transcription is one of its principal regulators. T3 governs cardiac contractility, heart rate, coronary vasodilation, and — critically — the regenerative response to ischemic injury. The DIO1 gene encodes type 1 iodothyronine deiodinase²² type 1 iodothyronine deiodinase
a selenoprotein enzyme expressed in liver, kidney, thyroid, and to a lesser extent cardiac tissue, one of three deiodinase enzymes responsible for converting the inactive prohormone T4 into active T3 in peripheral tissues.

The rs12095080 variant sits in the 3' untranslated region of DIO1, a regulatory segment that influences how efficiently the mRNA is translated into functional enzyme. Carriers of the G allele appear to have altered DIO1 activity that reduces T4→T3 conversion, creating a tissue-level hypothyroid state³³ tissue-level hypothyroid state
where serum TSH can be normal while cardiomyocytes receive less T3 than required for optimal function. This genetic predisposition becomes acutely dangerous in the context of myocardial infarction, when the cardiac thyroid axis is already under physiological stress.

The Mechanism

DIO1 is a selenocysteine-containing enzyme⁴⁴ selenocysteine-containing enzyme
requiring dietary selenium for its active site that catalyzes the outer-ring deiodination of T4 to T3. The rs12095080 G allele falls in the 3' UTR — a region that controls mRNA stability, translational efficiency, and microRNA binding. 3' UTR variants can reduce the half-life of messenger RNA or impair ribosomal loading⁵⁵ reduce the half-life of messenger RNA or impair ribosomal loading
resulting in lower steady-state protein levels without any change in the protein sequence. Lower DIO1 protein levels mean less T4→T3 conversion in the liver and periphery, reducing the systemic pool of active T3 available to the heart.

During and after acute myocardial infarction, non-thyroidal illness syndrome⁶⁶ non-thyroidal illness syndrome
also called low-T3 syndrome, a transient suppression of T4 deiodination driven by ischemic stress, inflammatory cytokines, and cortisol compounds this baseline genetic deficiency. Cardiomyocytes deprived of T3 show impaired sarco-endoplasmic reticulum calcium ATPase (SERCA2a)⁷⁷ sarco-endoplasmic reticulum calcium ATPase (SERCA2a)
the pump that cycles calcium between cytosol and SR to enable muscle relaxation expression, slower myosin heavy chain isoform switching toward the efficient β-MHC form, and reduced capacity for post-ischemic repair.

The Evidence

The primary evidence for rs12095080 comes from a 2020 longitudinal study by Kazukauskiene et al.⁸⁸ 2020 longitudinal study by Kazukauskiene et al.
published in Scientific Reports following 290 patients with acute myocardial infarction over a multi-year period. After controlling for NT-pro-BNP, ejection fraction, and other standard prognostic variables, the AG genotype at rs12095080 independently predicted cardiac-related mortality with an odds ratio of 3.97 (95% CI 1.45–10.89; p = 0.005)⁹⁹ odds ratio of 3.97 (95% CI 1.45–10.89; p = 0.005)
a nearly 4-fold increase in mortality risk. The overall variant (any G allele) carried a hazard ratio of 1.74 (95% CI 1.04–2.91; p = 0.034)¹⁰¹⁰ hazard ratio of 1.74 (95% CI 1.04–2.91; p = 0.034) as an independent predictor. The same group previously found rs12095080 associated with hypertension prevalence¹¹¹¹ rs12095080 associated with hypertension prevalence
in the same AMI cohort, suggesting broader cardiovascular vulnerability.

The genetic finding is reinforced by extensive evidence linking low T3 syndrome to cardiac outcomes independent of the variant. A 2017 meta-analysis of 41 studies by Wang et al.¹²¹² 2017 meta-analysis of 41 studies by Wang et al.
International Journal of Cardiology found that low-T3 syndrome in cardiovascular patients was associated with cardiac mortality HR 2.06 (95% CI 1.58–2.69)¹³¹³ cardiac mortality HR 2.06 (95% CI 1.58–2.69)
and all-cause mortality HR 2.52. In a large 2018 propensity-matched analysis of 2,459 AMI patients¹⁴¹⁴ 2018 propensity-matched analysis of 2,459 AMI patients
J Cardiol, low-T3 syndrome (prevalence 23.3%) was associated with in-hospital cardiovascular death rates of 4.7% vs. 1.7% in euthyroid controls, and adding the T3 status to the TIMI risk score meaningfully improved mortality prediction.

Important caveats: the rs12095080 cardiac mortality association comes from a single study with 290 patients from a Lithuanian AMI cohort. Replication in other populations and larger samples is needed before this variant reaches a higher evidence grade. The G allele is also notably rare in East Asian populations (essentially absent), so these findings may not generalize across all ancestries.

Practical Actions

If you carry the G allele, the central actionable priority is thyroid hormone monitoring — specifically free T3 (fT3), not just TSH and fT4. TSH is normal in low-T3 syndrome, making standard thyroid panels uninformative. Requesting an fT3 level as part of routine cardiovascular assessment establishes your baseline and provides meaningful prognostic context.

Selenium is the critical dietary cofactor for all three deiodinase enzymes. DIO1 is a selenoprotein — its catalytic activity depends on adequate selenium availability. Brazil nuts (one or two per day) provide the daily requirement¹⁵¹⁵ Brazil nuts (one or two per day) provide the daily requirement
each nut contains ~70–90 µg selenium; excess selenium is toxic, so do not exceed 400 µg/day. Selenium-containing foods and supplementation at 100–200 µg/day support baseline deiodinase capacity.

In any cardiac event or major illness, low-T3 syndrome develops rapidly and dramatically worsens outcomes. Ensuring your cardiologist is aware of this genetic finding before any planned cardiovascular procedures allows pre-emptive thyroid monitoring during hospitalization.

Interactions

DIO1 rs12095080 acts within the broader thyroid hormone metabolism pathway, which includes DIO2 (type 2 deiodinase¹⁶¹⁶ DIO2 (type 2 deiodinase
expressed in cardiac tissue and the brain, plays a complementary T4→T3 conversion role) and DIO3 (type 3 deiodinase¹⁷¹⁷ DIO3 (type 3 deiodinase
the inactivating enzyme that converts T3 to reverse T3). Variants in DIO2 (rs225014, rs2235544) and DIO3 (rs945006) may compound the effect of rs12095080 by further disrupting T3 availability. Additionally, SLCO1C1/OATP1C1 variants¹⁸¹⁸ SLCO1C1/OATP1C1 variants
controlling thyroid hormone transport into cells interact with deiodinase activity to determine net intracellular T3 levels. Individuals with variants in multiple thyroid hormone metabolism genes face the greatest aggregate risk of tissue-level hypothyroidism under cardiac stress.

Every cell in the body depends on zinc for more than 300 enzymes and 2,000+ transcription factors, yet the human body has no dedicated zinc storage organ — it must be continuously absorbed from food. In the intestine, that absorption flows primarily through a single gateway: ZIP4¹¹ ZIP4
The Zrt/Irt-like protein 4, encoded by SLC39A4 on chromosome 8q24.3, is the primary zinc importer on the apical surface of duodenal and jejunal enterocytes. When both copies of the SLC39A4 gene carry loss-of-function variants, dietary zinc cannot cross the gut wall. The consequence — hereditary acrodermatitis enteropathica (AE) — is a severe systemic zinc deficiency that is uniformly fatal without treatment but fully manageable with lifelong oral zinc supplementation.

The rs121434291 variant (C>T on the GRCh38 plus strand) replaces glycine at position 305 of the ZIP4 protein with aspartate (p.Gly305Asp; also annotated as p.Gly330Asp in the longer transcript isoform). Glycine 305 lies in a region of ZIP4 critical for normal protein folding and zinc transport function. The introduction of aspartate — a negatively charged, polar residue — in place of the structurally neutral glycine is predicted to disrupt the protein's tertiary structure and abolish zinc transport activity. An alternative alternate allele at the same position (C>A, p.Gly305Val) exists at similarly rare global frequency (~0.000004 in gnomAD exomes) and is expected to be pathogenic by the same mechanism, though it has not been independently classified in ClinVar.

The Mechanism

ZIP4 is expressed on the apical (luminal-facing) membrane of enterocytes, with expression upregulated when zinc levels fall. Its function is to move zinc ions from the intestinal lumen into absorptive cells, from where zinc enters circulation via basolateral zinc exporters. The Gly305Asp missense substitutes a bulky, charged residue into a structurally sensitive region of the ZIP4 protein. Because AE is autosomal recessive²² autosomal recessive
Both gene copies must be non-functional for disease; one functional copy is sufficient for normal zinc absorption, a single defective copy has no measurable impact on zinc status. Homozygotes and compound heterozygotes — who lose all functional ZIP4 activity — suffer progressive systemic zinc deficiency within weeks of birth, since the body cannot synthesise or store meaningful zinc reserves.

The Evidence

Wang et al. (2002)³³ Wang et al. (2002)
Wang K et al. A novel member of a zinc transporter family is defective in acrodermatitis enteropathica. Am J Hum Genet, 2002 identified SLC39A4 (encoding hZIP4) as the gene defective in AE by positional cloning and functional characterisation, establishing that pathogenic variants throughout the gene eliminate intestinal zinc uptake. The rs121434291 T allele (p.Gly305Asp) is classified Pathogenic in ClinVar (RCV000003720) on the basis of this gene identification and subsequent case-series evidence.

Küry et al. (2003)⁴⁴ Küry et al. (2003)
Küry S et al. Mutation spectrum of human SLC39A4 in a panel of patients with acrodermatitis enteropathica. Hum Mutat, 2003 documented seven additional SLC39A4 mutations in 12 AE families from France, Tunisia, Austria, and Lithuania — including missense, nonsense, and splice-site variants — confirming that pathogenic variants are distributed throughout the ZIP4 protein and that all result in clinically indistinguishable AE phenotype.

A later mutation update by Schmitt et al. (2009)⁵⁵ Schmitt et al. (2009)
Schmitt S et al. An update on mutations of the SLC39A4 gene in acrodermatitis enteropathica. Hum Mutat, 2009 catalogued 31 pathogenic SLC39A4 variants, confirming missense mutations as the most common class. Clinically, untreated AE presents in formula-fed infants within the first 4–10 weeks of life with a triad of acral and perioral dermatitis, diarrhoea, and alopecia. Breast-fed infants are typically protected by the high zinc bioavailability of breast milk and present at weaning. Without supplementation the disease is fatal; with it, prognosis is excellent.

Practical Implications

Oral zinc supplementation fully corrects the phenotype in homozygous AE patients. Treatment is initiated at 5–10 mg/kg/day of elemental zinc during the acute phase, then reduced to a maintenance dose of 1–2 mg/kg/day for life. Doses must be adjusted upward during growth phases, illness, and pregnancy. Regular monitoring of serum zinc is essential to avoid both deficiency relapse and zinc toxicity from over-supplementation.

Carriers (heterozygotes) are clinically unaffected under normal dietary conditions, but this variant is important for family planning: two carrier parents have a 25% probability of having an affected child with each pregnancy.

Interactions

AE illustrates how completely the body's zinc economy depends on ZIP4. Variants in other SLC39A (ZIP family) and SLC30A (ZnT family) genes modulate zinc homeostasis but do not cause AE. Dietary phytates in cereals and legumes form insoluble zinc complexes that compete with ZIP4-mediated uptake; this is especially relevant for heterozygous carriers whose single functional ZIP4 copy must operate efficiently. Co-administration of oral zinc with quinolone antibiotics (ciprofloxacin) or tetracyclines (doxycycline) should be separated by at least 2 hours to avoid chelation interactions that reduce absorption of both compounds.

Compound heterozygosity — carrying one copy of the Gly305Asp allele (rs121434291) on one chromosome and a different SLC39A4 pathogenic variant on the other — causes full AE and would appear as a CT genotype at this locus. A closely related variant, rs121434288 (p.Gly501Arg), affects a different conserved residue in the ZIP4 transmembrane domain and is catalogued separately in GeneOps.

Protein S is a vitamin K-dependent plasma glycoprotein that acts as an essential cofactor for activated protein C (APC)¹¹ activated protein C (APC)
Activated protein C is the central anticoagulant enzyme; it inactivates clotting factors Va and VIIIa to brake the coagulation cascade. Without adequate protein S, the APC system loses efficiency, and the body's ability to brake clot formation is compromised. The rs121918473 variant in the PROS1 gene — encoding the p.Asn258Ser substitution — was first identified in 1995 by Formstone and colleagues²² first identified in 1995 by Formstone and colleagues
In a pedigree with autosomal dominant protein S deficiency, using RT-PCR mutation screening as a cause of familial autosomal dominant thrombophilia. It is classified as pathogenic by ClinVar (VCV000013317) and is catalogued in OMIM as allelic variant 176880.0002.

Because this variant is so rare — absent in 594,524 gnomAD exome samples — it is not a common population polymorphism but a rare disease allele³³ rare disease allele
Most functional PROS1 disease variants are private or found only in small family pedigrees; the rarity makes population frequency estimates unreliable that is identified in clinical or research sequencing of families with unexplained thrombophilia.

The Mechanism

The PROS1 gene on chromosome 3 is transcribed from the minus strand. The rs121918473 variant creates an A>G transition in the coding sequence (NM_000313.4:c.773A>G), which on the plus strand corresponds to a T>C change at position 93,898,524 (GRCh38). The resulting amino acid substitution — asparagine to serine at position 258 (p.Asn258Ser) — falls within the fourth EGF-like domain of protein S⁴⁴ fourth EGF-like domain of protein S
Protein S contains four tandem EGF-like domains (residues 178–345) that mediate binding to activated protein C and phospholipid surfaces; the fourth domain is critical for cofactor activity.

Loss of asparagine at position 258 likely disrupts local domain folding or glycosylation, impairing protein S secretion or cofactor function. Heterozygous carriers produce only one functional PROS1 allele, resulting in reduced free protein S levels consistent with type I protein S deficiency⁵⁵ type I protein S deficiency
Type I (quantitative): both total and free protein S antigen levels are reduced, typically to 40–60% of normal in heterozygotes; type II is functional deficiency with normal antigen levels; type III has selectively reduced free protein S. Approximately 95% of protein S-deficient patients have type I or type III deficiency; missense variants like p.Asn258Ser typically cause type I.

The Evidence

The clearest quantification of risk from high-impact PROS1 variants comes from a 2025 population-scale study in JAMA⁶⁶ 2025 population-scale study in JAMA
Chaudhry et al., using UK Biobank and All of Us, n > 500,000: carriers of the most damaging PROS1 variants have protein S levels at 48% of normal and face an odds ratio of 14.01 (95% CI 6.98–27.14) for venous thromboembolism. A meta-analysis of 14 observational studies⁷⁷ meta-analysis of 14 observational studies
Di Minno et al., 4,955 VTE cases and 9,267 controls estimated OR 5.37 (95% CI 2.70–10.67) for a first VTE event across protein S-deficient individuals broadly, with broader uncertainty at the low end reflecting population heterogeneity.

A Danish clinical cohort⁸⁸ Danish clinical cohort
Larsen et al. 2021, n=22 PROS1 index cases found VTE frequency of 43% among PROS1 variant carriers versus 17% in those with normal PROS1 sequence (p=0.05), with significantly lower free protein S levels (0.51 vs 0.62 × 10³ IU/L) among carriers.

Protein S deficiency is not purely a venous thrombosis risk. A meta-analysis of inherited thrombophilias and arterial stroke⁹⁹ meta-analysis of inherited thrombophilias and arterial stroke
Chiasakul et al. 2019 found OR 2.26 (95% CI 1.34–3.80) for arterial ischemic stroke in protein S-deficient individuals, though the mechanism (paradoxical clot embolism vs. direct arterial effects) is less well understood.

For women, the interaction with combined hormonal contraceptives is critical: severe thrombophilias including protein S deficiency combined with oral contraceptive use carry RR 7.15 (95% CI 2.93–17.45)¹⁰¹⁰ RR 7.15 (95% CI 2.93–17.45)
van Vlijmen et al. 2016 systematic review, absolute risk 4.3–4.6 VTE events per 100 pill-years for VTE, representing a contraindication to estrogen-containing methods. Pregnancy itself carries a postpartum absolute risk of 4.2% (95% CrI 0.7–9.4%) for a first VTE event in protein S-deficient women per a Bayesian meta-analysis¹¹¹¹ Bayesian meta-analysis
Croles et al. 2017, BMJ.

Practical Implications

Because this is a pathogenic variant causing autosomal dominant thrombophilia, the clinical implications are immediate and concrete. First, free protein S levels should be measured to confirm biochemical deficiency and guide clinical classification. Second, all estrogen-containing hormonal preparations — combined oral contraceptives, patches, rings, and menopausal hormone therapy — are contraindicated. Third, standard surgical and perioperative thromboprophylaxis applies, with disclosure to all treating providers before any procedure. Fourth, anticoagulation after a first VTE event is managed in consultation with a hematologist, with particular attention to whether indefinite anticoagulation is warranted.

Direct oral anticoagulants (DOACs: rivaroxaban, apixaban, dabigatran) are the modern standard for thromboprophylaxis and treatment in most thrombophilia patients, though an ISTH SSC communication¹²¹² ISTH SSC communication
Kovac et al. 2024 notes that treatment failure risk is elevated in severe protein S deficiency (levels below 20%), warranting specialist management in those cases.

Interactions

The most clinically important interaction is with Factor V Leiden (rs6025, F5 R506Q)¹³¹³ Factor V Leiden (rs6025, F5 R506Q)
The most common inherited thrombophilia in Europeans (~5% carrier frequency); FVL prevents APC from inactivating factor Va — a complementary defect to protein S deficiency, which impairs APC cofactor activity. Simultaneous deficiency of protein S (impairing APC cofactor activity) and Factor V Leiden (making factor Va resistant to APC inactivation) attacks the same anticoagulant pathway at two separate points, compounding VTE risk substantially — likely exceeding the risk of either defect alone.

Similarly, Prothrombin G20210A (rs1799963, F2)¹⁴¹⁴ Prothrombin G20210A (rs1799963, F2)
Second most common inherited thrombophilia; raises prothrombin levels 30%, increasing substrate for thrombin generation would compound with protein S deficiency by simultaneously elevating the pro-coagulant drive while impairing the APC anticoagulant response.