Your cardiovascular system relies on a steady supply of hydrogen sulfide (H₂S), a gaseous signaling molecule produced mainly by an enzyme called cystathionine gamma-lyase (CTH/CSE)¹¹ cystathionine gamma-lyase (CTH/CSE)
also abbreviated CSE; encoded by the CTH gene on chromosome 1 in the liver and vascular wall. CTH sits at the junction of the transsulfuration pathway — the metabolic route that converts methionine into cysteine, consuming homocysteine along the way. H₂S produced by CTH acts as a vasodilator, an anti-inflammatory mediator, and a cardioprotective signal. The rs1021737 variant (c.1208G>T, p.Ser403Ile) is a common missense change in CTH's active domain that shifts the transsulfuration balance — with consequences that appear most sharply in women and in contexts of cardiovascular stress.

The Mechanism

CTH is a pyridoxal-5-phosphate (PLP)-dependent enzyme²² pyridoxal-5-phosphate (PLP)-dependent enzyme
PLP is the active form of vitamin B6, required as a cofactor for all transsulfuration enzymes that cleaves cystathionine to yield cysteine, α-ketobutyrate, and ammonia — with H₂S generated as a byproduct of cysteine catabolism. The Ser403Ile substitution replaces a serine residue (hydroxyl side chain, hydrophilic) with isoleucine (branched aliphatic, hydrophobic) near the enzyme's C-terminus.

Importantly, in vitro kinetic studies³³ in vitro kinetic studies
steady-state enzyme parameters Km and Vmax measured in purified recombinant protein show that Ser403Ile does not reduce CTH catalytic efficiency per se — Km, Vmax, and PLP cofactor loading are unchanged compared to wild-type. This distinguishes the polymorphism from pathogenic CTH mutations (e.g. T67I, Q240E) that directly impair catalysis. The mechanism by which TT homozygotes accumulate homocysteine likely involves altered protein stability, changed substrate channeling within the transsulfuration pathway, or a yet-uncharacterized regulatory effect on the enzyme's in-vivo activity under physiological conditions.

The downstream consequence is a relative insufficiency in two cardioprotective outputs of the transsulfuration pathway: (1) less homocysteine cleared through the cystathionine route, raising plasma total homocysteine (tHcy), and (2) potentially reduced H₂S bioavailability. Elevated tHcy damages vascular endothelium through oxidative stress, impairs nitric oxide signaling, promotes thrombosis, and accelerates atherosclerosis. H₂S exerts the opposite effects — it relaxes vascular smooth muscle, reduces inflammation, and protects cardiomyocytes during ischemia-reperfusion injury.

The Evidence

The foundational association was established by Wang et al. 2004⁴⁴ Wang et al. 2004
Wang J, Huff AM, Spence JD, Hegele RA. Single nucleotide polymorphism in CTH associated with variation in plasma homocysteine concentration. Clin Genet. 2004;65(6):483–486. in 496 Caucasian subjects. TT homozygotes had significantly higher mean plasma total homocysteine than GG or GT carriers, with an effect size comparable to that of the well-established MTHFR C677T variant. This places CTH Ser403Ile alongside MTHFR as a meaningful genetic determinant of plasma homocysteine — a finding that was replicated in subsequent studies.

The most striking clinical evidence comes from a 2022 Swedish case-control study Söderström et al. 2022⁵⁵ Söderström et al. 2022
CTH G1208T and MTHFR A1298C polymorphisms are associated with a higher risk of a first myocardial infarction with fatal outcome among women. Drug Metab Pers Ther. 2022. of first myocardial infarction. In women, the CTH G1208T variant (rs1021737 T allele carriers) showed a striking OR of 3.14 [95% CI 1.16–8.54] for fatal MI in heterozygotes, rising to OR 3.22 [1.22–8.51] in the dominant model. Critically, this association was entirely sex-specific: no comparable effect was found in men, and neither genotype was associated with non-fatal MI in either sex. The MTHFR A1298C variant showed a similar sex-specific pattern, suggesting that transsulfuration-pathway impairment preferentially raises fatal cardiovascular risk in women — possibly because estrogen normally supports transsulfuration enzyme activity and H₂S production, making women more vulnerable to CTH impairment.

H₂S bioavailability was further implicated by Rajpal et al. 2018⁶⁶ Rajpal et al. 2018
Total sulfane sulfur bioavailability reflects ethnic and gender disparities in cardiovascular disease. Redox Biol. 2018;15:480–489. in 324 CVD patients and controls. The CTH T allele was significantly more common in CVD patients vs controls, and Caucasian females with CVD had markedly lower acid-labile sulfide (a key H₂S pool) than controls — pointing to CTH-mediated H₂S insufficiency as a sex-stratified cardiovascular mechanism.

The variant does not appear to affect essential hypertension risk in Chinese Han populations (Li et al. 2008; n=993; p>0.05 for both allele and genotype associations), nor was it significantly associated with preeclampsia in Caucasians (Mrozikiewicz et al. 2015). Its effects may be most pronounced in the context of cumulative cardiovascular risk and the sex-hormone environment.

Practical Actions

For TT homozygotes and GT carriers — particularly women — the actionable focus is: (1) driving homocysteine levels down through the methylation pathway (B-vitamins), (2) monitoring homocysteine as a cardiovascular biomarker, and (3) supporting H₂S production through dietary sulfur amino acid intake. The CTH pathway requires vitamin B6 (as PLP) as its essential cofactor; ensuring adequate B6 status is especially important when CTH function may be compromised.

Because CTH Ser403Ile impairs homocysteine clearance through the transsulfuration route, the MTHFR re-methylation pathway becomes more important as an alternative route. This means active B12 and methylfolate availability become critical backstops — especially for individuals who also carry MTHFR C677T (rs1801133) or A1298C (rs1801131), where the re-methylation route is also impaired.

Interactions

CTH rs1021737 operates downstream of MTHFR in the one-carbon/homocysteine cycle. The transsulfuration pathway (CTH's route) and the re-methylation pathway (MTHFR's route) are the two main homocysteine disposal routes. If both are impaired — by CTH TT and MTHFR C677T TT simultaneously — homocysteine has nowhere to go, creating a severely elevated tHcy state. This dual-impairment combination is the most clinically relevant interaction for rs1021737 and warrants a compound action.

MTRR rs1801394 (methionine synthase reductase) supports the re-methylation pathway; MTRR risk genotypes further reduce the pressure on the transsulfuration route. Any combination of CTH TT + folate pathway impairment should be flagged for aggressive B-vitamin optimization and homocysteine monitoring.

Glucokinase (GCK) is the primary glucose sensor in pancreatic beta cells — the enzyme that tells the cell how much insulin to release in response to rising blood sugar. Variants in the GCK gene itself cause maturity-onset diabetes of the young type 2 (MODY2)¹¹ maturity-onset diabetes of the young type 2 (MODY2)
GCK-MODY is an autosomal dominant monogenic form of diabetes caused by heterozygous loss-of-function mutations in GCK. Individuals have a mild, stable fasting hyperglycemia throughout life, typically 5.5–8 mmol/L (99–144 mg/dL), that rarely progresses to diabetic complications and influence fasting glucose levels across the general population. rs10278336 sits in an intron of YKT6, a v-SNARE vesicle-trafficking gene immediately downstream of GCK on chromosome 7, within the same haplotype block. It is a GWAS-identified common risk allele for type 2 diabetes.

The Mechanism

rs10278336 is located at chr7:44,205,764 (GRCh38), approximately 16 kb downstream of the GCK gene body in an intron of YKT6. Its functional annotation is intronic — it does not change an amino acid or a splice site. The most likely mechanism is regulatory: variants in this region may sit in enhancer elements that influence GCK expression in pancreatic beta cells or hepatocytes, or may tag a broader GCK haplotype through linkage disequilibrium with functional variants upstream. The A allele at rs10278336 co-segregates with alleles at nearby GCK-region SNPs (including the well-studied rs4607517²² rs4607517
GCK intronic fasting-glucose GWAS variant; A allele at ~17% global frequency associates with elevated fasting glucose by reducing GCK expression) that collectively reduce GCK-mediated glucose sensing efficiency in beta cells. The net result is a small but detectable rightward shift in the fasting glucose set-point — the blood sugar level at which the pancreas judges it is time to release more insulin.

Unlike rare MODY2 mutations (which cause large, clinically obvious glucose elevations), common GWAS variants in this region each contribute only a few milligrams per deciliter to fasting glucose. Their cumulative effect, combined with other T2D risk variants, can meaningfully increase lifetime diabetes risk.

The Evidence

The primary evidence comes from a large-scale GWAS of type 2 diabetes³³ large-scale GWAS of type 2 diabetes
Morris AP et al. Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nature Genetics 2012 (GWAS Catalog GCST005047), which identified rs10278336-A as a T2D risk allele with OR=1.07 (95% CI 1.04–1.10) at p=1×10⁻⁶. The effect size is modest and consistent with the polygenic architecture of common T2D — no single common variant contributes large risk, but the aggregate of dozens of variants is substantial.

A case-control study in 835 Chinese women with gestational diabetes⁴⁴ case-control study in 835 Chinese women with gestational diabetes
She L et al. Association of glucokinase gene and glucokinase regulatory protein gene polymorphisms with gestational diabetes mellitus. Gene 2022 (835 GDM cases, 870 controls) examined rs10278336 as part of a GCK-region panel and found no significant GDM association (P>0.05). However, the nearby GCK variant rs1799884 (-30G>A promoter) was significantly associated with GDM in this population, reinforcing that functional GCK-region variants do contribute to gestational glucose dysregulation. The null result for rs10278336 in GDM may reflect population-specific LD patterns in Han Chinese or insufficient statistical power.

The EPIC-InterAct study⁵⁵ EPIC-InterAct study
Langenberg C et al. Gene-lifestyle interaction and type 2 diabetes: the EPIC interact case-cohort study. PLoS Med 2014 (12,403 T2D cases from 9 European countries) found that the collective effect of T2D risk alleles — including GCK-region variants — was amplified in leaner individuals (p=7.49×10⁻⁹ for waist circumference interaction) and younger adults, but was not significantly modified by diet or physical activity at the population level. This suggests that for carriers of multiple T2D risk alleles, maintaining lean body composition may partially offset genetic risk.

Practical Actions

rs10278336 contributes a modest per-allele effect on T2D risk (OR~1.07). Its practical relevance is primarily in the context of cumulative polygenic risk: carriers of AA who also carry risk alleles at TCF7L2, KCNJ11, PPARG, and other major T2D loci face meaningfully elevated lifetime risk that justifies proactive monitoring. For AA homozygotes, fasting glucose tracking and attention to postprandial glucose patterns provide early warning before overt dysglycemia.

Dietary carbohydrate quality — specifically favoring low-glycemic-index sources — reduces the demand on GCK-mediated beta-cell glucose sensing: when carbohydrate absorption is slower, the acute glucose spike requiring precise incretin and GCK-mediated detection is lower. This is a mechanism-aligned (not generic) dietary strategy for this variant.

Interactions

rs10278336 acts within the GCK haplotype block. The neighboring rs4607517 (GCK intron) is a well-established fasting-glucose variant — both variants tag similar GCK regulatory architecture, and they are often in LD. A compound action covering both rs10278336-AA and rs4607517-AA (homozygous risk at both GCK-region loci) would be appropriate if both are confirmed in the same individual: the combined effect on fasting glucose would be additive.

The GCKR gene (glucokinase regulatory protein, chromosome 2) encodes a protein that sequesters GCK in the liver nucleus, buffering its activity. Variants in GCKR (rs780094, rs1260326) interact with GCK-region variants by modulating how much GCK activity is "released" to process incoming glucose. Carriers of GCK-region risk alleles who also carry GCKR high-expression variants face compounded impairment of hepatic glucose disposal.

Methionine synthase reductase (MTRR) is an enzyme whose sole job is to keep another enzyme — methionine synthase (MTR) — running. MTR converts homocysteine to methionine using methylcobalamin¹¹ methylcobalamin
The methyl-carrying, active form of vitamin B12 (active B12) as a cofactor, but during each catalytic cycle the B12 becomes oxidized to an inactive form. MTRR reduces it back to active methylcobalamin so MTR can continue. Without efficient MTRR, the methylation cycle slows, homocysteine accumulates, and the production of SAM — the universal methyl donor for DNA, RNA, and protein methylation — declines.

The His595Tyr missense variant (rs10380, c.1783C>T) swaps histidine for tyrosine at position 595 of the MTRR protein. This amino acid change sits in a functionally important region and was co-identified with the intronic variant rs162049 in the same pancreatic cancer risk haplotype — a haplotype shown in functional cell-line studies to produce less MTRR protein, elevated homocysteine in culture medium, and reduced LINE-1 methylation²² LINE-1 methylation
LINE-1 (long interspersed element-1): repetitive DNA sequences whose methylation status is a proxy for genome-wide DNA methylation. The T allele is the minor allele globally, occurring in about 11% of chromosomes in gnomAD exomes, with substantial variation by ancestry (approximately 33–34% in African and Latino populations, versus ~9–10% in Europeans).

The Mechanism

His595 is a conserved residue in MTRR's methyltransferase domain. The histidine-to-tyrosine substitution is predicted to alter the local protein fold, reducing the enzyme's ability to perform the reductive methylation of cob(II)alamin³³ reductive methylation of cob(II)alamin
The chemical step by which MTRR converts oxidized, inactive cobalamin back to methylcobalamin for MTR back to active methylcobalamin. Functional data from risk-haplotype transfectants (containing both rs10380 and rs162049) confirmed that cells carrying the risk haplotype produced less MTRR protein, higher extracellular homocysteine, and lower LINE-1 methylation compared to wild-type cells — consistent with impaired B12 cycling and downstream methylation deficit.

The Evidence

The discovery study by Ohnami et al.⁴⁴ Ohnami et al.
Ohnami S et al. His595Tyr polymorphism in the methionine synthase reductase gene is associated with pancreatic cancer risk. Gastroenterology, 2008 found rs10380 associated with pancreatic cancer risk (OR 1.45, 95% CI 1.11–1.88; P=0.0063) in a multicenter Japanese case-control study of 317 cases and 1,232 controls, surviving permutation testing for multiple comparisons (P=0.023 recessive model). The association was attributed to impaired methylation-dependent regulation of tumor-suppressor genes.

A California population-based study by Shaw et al.⁵⁵ Shaw et al.
Shaw GM et al. 118 SNPs in folate-related genes and spina bifida risk. BMC Med Genet, 2009 found heterozygous or homozygous rs10380 carriers had OR 3.4 (95% CI 1.6–7.1) for spina bifida risk among 259 cases and 359 controls — a striking result that implicates impaired MTRR function in neural tube closure, where adequate methylation is essential for proper gene regulation during embryogenesis. No association was found with conotruncal heart defects in the same study.

A more recent Han Chinese case-control study of 595 children with congenital heart disease⁶⁶ congenital heart disease
CHD: structural heart abnormalities present from birth, often linked to disrupted epigenetic regulation during cardiac development and 605 controls found the TT genotype associated with CHD risk (OR 2.27, 95% CI 1.20–4.31). Critically, maternal folic acid supplementation during pregnancy reduced CHD risk (OR 0.55), suggesting that adequate one-carbon supply during cardiac development can partially compensate for impaired MTRR function.

Other studies examining rs10380 in colorectal cancer, gastric cancer, and DNA methylation have reported null or inconsistent results — consistent with a variant whose effect is tissue-context-dependent and modified by folate/B12 nutritional status.

Practical Implications

The His595Tyr substitution impairs MTRR's ability to recycle B12 for MTR. The primary intervention strategy is to maintain high active-B12 supply to compensate for reduced recycling efficiency: more incoming methylcobalamin means MTR is less dependent on MTRR-mediated recycling to stay active. Hydroxocobalamin is an especially useful form because it enters both methylcobalamin and adenosylcobalamin pathways. Methylfolate (5-MTHF) upstream of MTR keeps the methyl-donor pool full. Homocysteine measurement is a direct functional readout — elevated homocysteine (above 10 µmol/L) signals that the remethylation pathway is running below capacity. The CHD data further suggest that periconceptional folic acid use is especially important for T allele carriers planning a pregnancy.

Interactions

rs10380 and rs162049 (intronic MTRR) were co-identified in the same functional haplotype; individuals carrying both likely have a compounded impairment of MTRR expression and enzymatic efficiency. Combined with MTRR A66G (rs1801394), which reduces enzyme efficiency at the protein level via p.Ile22Met, or with MTR A2756G (rs1805087), which reduces MTR activity directly, the overall B12 cycling capacity is further curtailed. The most clinically significant interaction is with MTHFR C677T (rs1801133): reduced methylfolate supply from MTHFR combined with impaired B12 recycling from MTRR creates dual pressure on homocysteine remethylation that neither variant produces alone.

The CHI3L1 gene encodes YKL-40, one of the most clinically informative biomarkers of tissue inflammation — elevated in asthma, COPD, rheumatoid arthritis, and coronary artery disease. How much YKL-40 your body produces is largely genetically predetermined, and the CHI3L1 locus on chromosome 1q32.1 contains multiple independent regulatory variants that each contribute to this set-point. rs10399931 is one of those variants — a 2 kb upstream regulatory SNP that modulates CHI3L1 mRNA levels independently of the primary promoter variant (rs4950928), adding a distinct layer of genetic control over the inflammatory YKL-40 axis.

The Mechanism

rs10399931 lies ~2,000 base pairs upstream of the CHI3L1 transcription start site, in a region with regulatory potential. Unlike the promoter SNP rs4950928, which alters a transcription factor binding site, rs10399931 appears to operate through a post-transcriptional mechanism¹¹ post-transcriptional mechanism
The T allele reduces CHI3L1 mRNA abundance without altering promoter-level reporter activity, implying effects on mRNA stability, splicing-associated regulation, or upstream chromatin context rather than direct promoter binding. The common C allele (plus strand, ~73% globally) is associated with higher CHI3L1 mRNA levels and consequently higher circulating YKL-40. The rarer T allele (~27%) reduces mRNA expression and lowers the YKL-40 set-point.

The Guerra et al. birth cohort study²² Guerra et al. birth cohort study
Guerra S et al. Genetic and epigenetic regulation of YKL-40 in childhood. J Allergy Clin Immunol, 2018 studied 68 CHI3L1 SNPs in up to 2,405 children and identified rs10399931 among 7 variants independently associated with YKL-40. Critically, alleles linked to lower YKL-40 (including the T allele at rs10399931) were associated with higher DNA methylation at five CHI3L1 CpG sites, suggesting the genetic effect is partly mediated through epigenetic remodeling — an upstream regulatory SNP shaping chromatin accessibility at the CHI3L1 locus.

The Evidence

The Taiwanese population study by Tsai et al.³³ Tsai et al.
Tsai Y et al. CHI3L1 polymorphisms associate with asthma in a Taiwanese population. BMC Med Genet, 2014 genotyped rs10399931 in 628 asthma patients and controls. The CC genotype (plus-strand, equivalent to the paper's "GG" in minus-strand gene notation) was significantly associated with asthma risk (aOR=1.77, 95% CI 1.13–2.77) and the highest serum YKL-40 levels (133 ng/mL vs 91 ng/mL for TT carriers). YKL-40 levels in CC carriers correlated with FEV1 decline (P=0.004) and FVC reduction (P=0.011), and this effect persisted even among patients on inhaled corticosteroids — suggesting rs10399931-linked YKL-40 elevation is refractory to standard anti-inflammatory therapy.

A meta-analysis across 16 publications involving 5,005 cases and 9,725 controls⁴⁴ 16 publications involving 5,005 cases and 9,725 controls
Huang QP et al. Assessment of the Association between Genetic Polymorphisms in the CHI3L1 Gene and Asthma Risk. Int Arch Allergy Immunol, 2022 confirmed that the TT genotype is protective against asthma in East Asian populations (TT vs CC: OR=0.77, 95% CI 0.61–0.98, P=0.030), with population allele frequencies making this signal particularly detectable in Asian cohorts.

A functional study by Chen et al.⁵⁵ Chen et al.
Chen G et al. Functional study of the association of CHI3L1 polymorphisms with asthma susceptibility in the Southwest Chinese Han population. Biosci Rep, 2019 directly measured mRNA expression in human subjects — CT/TT genotype carriers showed significantly reduced CHI3L1 mRNA (P=0.002), establishing molecular causality.

In the coronary artery disease context, Chou et al.⁶⁶ Chou et al.
Chou HH et al. Circulating YKL-40 levels but not CHI3L1 or TRIB1 gene variants predict long-term outcomes in multivessel CAD. Sci Rep, 2024 used conditional analysis to confirm rs10399931 as an independent YKL-40 determinant at the CHI3L1 locus, distinct from rs4950928 — meaning the two variants capture separate genetic contributions to YKL-40 variation.

Practical Actions

For CC carriers, the practical implications mirror those of the rs4950928 CC genotype: an elevated YKL-40 baseline reflects a higher inflammatory setpoint. However, because rs10399931 acts at a different regulatory level than rs4950928, a CC carrier at this locus who also carries CG or GG at rs4950928 has a compound picture — neither variant alone describes the full genetic contribution to their YKL-40 level. Serum YKL-40 testing is the direct readout of both variants together.

For TT carriers, the low-YKL-40 genotype confers partial protection against asthma and airway inflammation, but because TT is uncommon (~7%), clinical reference ranges are calibrated to the CC majority and may flag TT-typical low values as anomalous.

Interactions

rs10399931 operates in the same YKL-40 QTL cluster as rs4950928 (the primary CHI3L1 promoter variant, also in related_snps here) and rs872129, which each contribute independent signals to circulating YKL-40 levels confirmed by sequential conditional analyses (Chou et al. 2024, PMID 39592699). The intronic variant rs12141494 independently affects airway tissue YKL-40 expression and lung function, adding a fourth distinct CHI3L1 regulatory layer. Together, these variants can substantially compound or offset one another's contributions to the inflammatory setpoint.

The CYP2B6 enzyme accounts for only 1–4% of total hepatic cytochrome P450 content yet handles metabolism of several high-stakes medications: the HIV antiretroviral efavirenz, the addiction-treatment drug methadone, the antidepressant and smoking-cessation aid bupropion, the cancer drug cyclophosphamide, and the anesthetic and antidepressant ketamine. CYP2B6 is the most polymorphic human CYP enzyme¹¹ CYP2B6 is the most polymorphic human CYP enzyme
more than 100 defined alleles make CYP2B6 activity highly variable between individuals, producing up to 100-fold inter-individual differences in plasma drug concentrations at standard doses.

rs10403955 is an intronic variant (c.172-468T>G) located deep in intron 1 of CYP2B6 on chromosome 19. It does not alter the enzyme's amino acid sequence but marks a haplotype block in the gene's regulatory and intronic architecture that tags reduced CYP2B6 metabolic activity at the population level.

The Mechanism

Unlike the well-characterized missense variant rs3745274 (516G>T, p.Q172H), rs10403955 sits 468 base pairs upstream of exon 2 within an intronic region. Its G allele acts as a tagging SNP²² tagging SNP
a variant in strong linkage disequilibrium with one or more functional variants in the same haplotype block, allowing it to predict activity even without itself being causal. The T allele belongs to haplotype blocks associated with higher CYP2B6 expression and enzyme activity; the G allele marks haplotypes where expression or splicing is impaired.

Wang et al. (2011)³³ Wang et al. (2011)
CYP2B6 polymorphisms influence the plasma concentration and clearance of the methadone S-enantiomer identified a trinucleotide haplotype block (rs8100458–rs10500282–rs10403955) in intron 1 of CYP2B6 where the T allele of rs10403955 was significantly associated with higher S-methadone clearance and lower plasma concentrations (P < 0.0017). Conversely, carriers of the G allele have lower clearance — resulting in higher drug accumulation at standard doses.

The Evidence

In a prospective study of 366 Taiwanese patients receiving methadone maintenance therapy, Wang et al. demonstrated⁴⁴ Wang et al. demonstrated
CYP2B6 polymorphisms influence the plasma concentration and clearance of the methadone S-enantiomer that the rs10403955 T/G haplotype block explained a meaningful fraction of inter-individual variability in S-methadone plasma concentration-to-dose ratios. The S-enantiomer is the pharmacologically active form responsible for mu-opioid receptor binding and opioid maintenance effects, and also the form that prolongs cardiac QT intervals, making its accumulation clinically significant.

For efavirenz, Carr et al. (2010)⁵⁵ Carr et al. (2010)
Haplotype structure of CYP2B6 and association with plasma efavirenz concentrations in a Chilean HIV cohort included rs10403955 as one of three CYP2B6 tagging SNPs representing 11 associated variants. Among 219 HIV-infected Chilean patients, a composite genetic model containing these tagging SNPs predicted efavirenz concentrations exceeding the CNS-toxicity threshold of 4 μg/mL with an odds ratio of 48.1 (95% CI 13.5–207.7; P < 0.001). The model had 97% specificity and 92% negative predictive value, demonstrating that CYP2B6 haplotype tagging — including this locus — can meaningfully stratify toxicity risk before treatment.

The G allele frequency varies by ancestry: approximately 25% in Europeans, 37% in Africans, 21% in East Asians, 39% in South Asians, and 35% in Latino populations, reflecting the same broad ancestry stratification seen in the primary CYP2B6 functional variant rs3745274.

Practical Implications

The haplotype context of rs10403955 makes it clinically informative for anyone prescribed a CYP2B6-metabolized drug. Homozygous GG carriers represent approximately 9% of the global population and are most likely to experience reduced enzymatic clearance, leading to higher plasma concentrations of efavirenz (with CNS toxicity risk), methadone (with QT-prolongation and accumulation risk), and other CYP2B6 substrates. Heterozygous GT carriers (~42% of the population) show intermediate accumulation.

Importantly, rs10403955 should be interpreted alongside primary functional variants in CYP2B6 (especially rs3745274 and rs28399499) for the most accurate metabolizer phenotype assignment. CPIC guidelines define CYP2B6 metabolizer phenotypes primarily through the star-allele system: poor metabolizers should receive efavirenz dose reductions to 400 mg or 200 mg daily instead of the standard 600 mg.

Interactions

rs10403955 lies within a haplotype block with rs8100458 and rs10500282, forming an intronic triad that modulates CYP2B6 expression as a unit. The primary functional variants in CYP2B6 — rs3745274 (516G>T, defining CYP2B6*6/*9) and rs28399499 (983T>C, defining CYP2B6*16/*18) — are the classical anchors of the CYP2B6 star-allele system, and rs10403955 captures independent haplotype information not fully explained by those coding variants alone. CYP2B6 activity is also strongly inducible by rifampin and by efavirenz itself, meaning that drug-drug interactions can partially overcome or exaggerate the genetic effect depending on co-medications.

Cortisol, the body's primary stress hormone, communicates with cells through the glucocorticoid receptor (GR) encoded by NR3C1. When cortisol binds the GR, it triggers gene expression programs that regulate inflammation, blood glucose, immune function, and fat distribution — all processes central to both metabolic health and the biology of aging. The rs10482605 variant sits in the promoter region of NR3C1, upstream of the coding sequence, where it influences how much GR protein the cell produces in the first place. Less GR means weaker cortisol signaling — a subtle but persistent dampening of the hormone's metabolic and anti-inflammatory effects that compounds over a lifetime.

Two findings from independent research programs define this SNP's clinical significance. First, Kumsta et al. (2009)¹¹ Kumsta et al. (2009)
Kumsta R et al. Characterization of a glucocorticoid receptor gene (GR, NR3C1) promoter polymorphism reveals functionality and extends a haplotype with putative clinical relevance. Am J Med Genet B Neuropsychiatr Genet. 2009;150B(4):476-82 demonstrated in reporter gene assays that the risk allele reduces GR transcriptional activity in brain-derived cell lines under both basal and stimulated conditions — making this a functionally characterized regulatory variant, not merely an association signal. Second, Kolb et al. (2023)²² Kolb et al. (2023)
Kolb KL et al. Glucocorticoid Receptor Gene (NR3C1) Polymorphisms and Metabolic Syndrome: Insights from the Mennonite Population. Genes (Basel). 2023;14(9):1805 found that homozygotes for the risk allele had a 4.74-fold increased odds of metabolic syndrome — a finding that places this promoter variant in direct conversation with the global epidemic of insulin resistance, central obesity, hypertension, and dyslipidemia.

The Mechanism

rs10482605 is located at chromosome 5, position 143,403,956 (GRCh38), within the promoter region of NR3C1. Because NR3C1 spans the minus strand of chromosome 5, the alleles described in published papers use coding-strand notation: the T allele (coding strand) corresponds to the A allele on the plus strand (reference, major), and the C allele (coding strand) corresponds to the G allele on the plus strand (alternate, minor). Genome files report plus-strand alleles, so the risk allele in this database is G (plus strand, ~32% global frequency).

The NR3C1 gene uses multiple alternative first exons — nine non-coding exons (1A, 1B, 1C, 1D, 1E, 1F, 1H, 1I, 1J) that permit tissue-specific and stimulus-specific control of GR expression. The rs10482605 position maps to a CpG island³³ CpG island
A stretch of DNA where CpG dinucleotides are present at higher than expected frequency — CpG islands in gene promoters often regulate transcriptional activity and are sensitive to methylation-based silencing in the 5'UTR region. The G (risk) allele alters the sequence within this regulatory context, reducing the transcriptional drive on GR expression. In reporter assays, this reduction in activity was observed under both resting and stimulated conditions, suggesting a constitutive dampening of GR production rather than a context-specific effect.

The variant sits in high linkage disequilibrium⁴⁴ linkage disequilibrium
LD means two variants are inherited together more often than chance would predict — when one is present, the other tends to be too with rs6198, the 9β variant in NR3C1's 3'UTR that stabilizes the glucocorticoid-resistant GRβ mRNA isoform. This creates a double-hit haplotype: the G allele at rs10482605 reduces GR transcription at the promoter, while the co-occurring C allele at rs6198 shifts the expressed mRNA toward the dominant-negative GRβ isoform. The net result is less GR protein being produced, and a larger fraction of what is produced being the cortisol-resistant form — a compounding attenuation of glucocorticoid signaling.

The Evidence

Kumsta et al. (2009) genotyped 219 subjects and performed in vitro reporter gene assays to establish that rs10482605 is a functional variant. The functional evidence — reduced transcriptional activity in the risk allele — is the foundation for understanding why this SNP matters biologically. The observation of high LD with rs6198 extended the known NR3C1 haplotype architecture and proposed a mechanism for depression risk: blunted GR expression reduces negative feedback on the HPA axis, permitting prolonged cortisol elevation under stress.

The metabolic syndrome connection was established by Kolb et al. (2023) in a genetically isolated Brazilian Mennonite community — a founder population with reduced genetic background noise. The study genotyped 74 MetS cases and 138 unaffected controls (212 total), plus a replication set of 236 individuals. Homozygosity for the risk allele (G/G on plus strand; C/C on coding strand) was associated with OR = 4.74 (95% CI 1.10–20.28, pcorr = 0.024). An independent haplotype analysis confirmed this signal (TTCGTTGATT haplotype, OR = 4.74, pcorr = 0.048). Critically, the association was independent of age, physical activity, and family environment — pointing to a direct genetic contribution to metabolic risk rather than a lifestyle confounder.

The evidence level is rated moderate: the functional characterization (Kumsta 2009) is solid, and the metabolic syndrome OR of 4.74 is large. However, the metabolic finding comes from a single study in a founder population (which both boosts power and limits generalizability), and cross-population replication is not yet established. The mechanistic link between reduced GR expression and metabolic syndrome — while biologically plausible — is not directly demonstrated in the Mennonite study.

The biological pathway is well-supported: chronic GR insufficiency impairs the normal glucocorticoid suppression of inflammatory cytokines, and chronic low-grade inflammation is a driver of insulin resistance, central fat deposition, and dyslipidemia — the cardinal features of metabolic syndrome. Furthermore, GR-mediated transcriptional regulation of hepatic lipid metabolism genes⁵⁵ GR-mediated transcriptional regulation of hepatic lipid metabolism genes
The GR directly regulates PCSK9, BHLHE40, and SREBP-2 pathway genes in liver cells, meaning altered GR activity can independently shift cholesterol and triglyceride metabolism provides an additional route by which reduced GR transcription could produce dyslipidemia.

Practical Implications

For homozygous GG carriers (~10% of the population globally), the combination of reduced GR transcriptional activity and the metabolic syndrome risk signal warrants proactive metabolic monitoring. The actionable targets are the five components of metabolic syndrome: waist circumference, fasting glucose, blood pressure, triglycerides, and HDL cholesterol. GG carriers who develop metabolic syndrome may have a component that is driven by impaired glucocorticoid signaling — meaning that lifestyle interventions need to account for the possibility that standard inflammatory and metabolic setpoints are shifted at the receptor level.

For AG heterozygotes (~44% of the population), the functional reduction in GR expression is partial. The metabolic syndrome OR in heterozygotes is not separately reported in the available literature, but given the additive inheritance pattern of most GR variants, a graded effect is likely. Heterozygotes benefit from metabolic awareness without the same urgency as GG homozygotes.

For AA homozygotes (~46%), full GR promoter activity is maintained. This genotype represents the ancestral configuration with no identified metabolic risk from this specific locus.

Interactions

rs10482605 is in high linkage disequilibrium with rs6198 (9β)⁶⁶ rs6198 (9β), the NR3C1 3'UTR variant that stabilizes GRβ mRNA and blunts glucocorticoid sensitivity. When both risk alleles are present on the same haplotype, the individual faces both reduced GR transcription (rs10482605 G) and a shift toward the dominant-negative GRβ isoform (rs6198 C). This compound haplotype was the original focus of Kumsta et al. (2009) and likely represents the maximum NR3C1-driven glucocorticoid resistance achievable from coding and promoter variation in this gene.

Within the same NR3C1 gene, two longevity-associated intronic variants are already catalogued: rs2963154⁷⁷ rs2963154 and rs10515522⁸⁸ rs10515522, both from the Polish centenarian cohort. Those variants associate with survival to extreme old age and altered cholesterol metabolism through GR-driven hepatic lipid regulation. rs10482605 adds the metabolic syndrome dimension — showing that the same gene's promoter regulation influences metabolic risk decades before extreme longevity outcomes become observable.

The BclI variant rs41423247⁹⁹ BclI variant rs41423247 sits on the same gene and increases GR sensitivity. Combined NR3C1 haplotype analysis — incorporating rs10482605 (promoter activity), rs6198 (GRβ isoform balance), and rs41423247 (receptor sensitivity) — represents the most comprehensive picture of an individual's glucocorticoid receptor biology. No single-study haplotype analysis covering all four variants has been published, but the mechanistic logic supports composite profiling.

A single molecular accident — a C-to-T transition in exon 4 of CYP17A1 — replaces arginine at position 239 with a premature stop codon. The resulting truncated protein retains only the first 238 of 508 amino acid residues, lacking both the substrate-binding pocket and the heme coordination site that give CYP17A1¹¹ CYP17A1
Cytochrome P450 17α-hydroxylase/17,20-lyase, a bifunctional enzyme encoded on chromosome 10q24.32 that catalyzes two sequential reactions required for cortisol and sex steroid biosynthesis its catalytic power. The protein is non-functional before it is even folded. The result is the most complete form of 17α-hydroxylase/17,20-lyase deficiency (17OHD)²² 17α-hydroxylase/17,20-lyase deficiency (17OHD)
A rare form of congenital adrenal hyperplasia accounting for approximately 1% of all CAH cases, caused by biallelic loss-of-function variants in CYP17A1.

The Mechanism

CYP17A1 sits at the branch point of steroidogenesis in both the adrenal cortex and the gonads. Its first reaction, 17α-hydroxylation³³ 17α-hydroxylation
Addition of a hydroxyl group to the 17-position of pregnenolone and progesterone, generating the precursors required for cortisol synthesis, converts pregnenolone to 17α-hydroxypregnenolone and progesterone to 17α-hydroxyprogesterone. Its second reaction, 17,20-lyase activity⁴⁴ 17,20-lyase activity
Cleavage of the C17–C20 bond to produce DHEA and androstenedione, the entry point into sex steroid synthesis, is what allows the adrenal gland and gonads to produce sex hormone precursors at all.

Without CYP17A1 function, three simultaneous hormonal disasters unfold. First, cortisol is absent; ACTH rises without feedback restraint, driving bilateral adrenal hyperplasia. Second, the pregnenolone pool that cannot enter the cortisol or sex steroid branches is shunted instead to mineralocorticoids — specifically 11-deoxycorticosterone (DOC) and corticosterone — causing sodium retention, hypertension, and hypokalemia that can persist for years before diagnosis. Third, neither the adrenal gland nor the gonads can synthesize DHEA or androstenedione, so no sex steroids are produced from any source.

The Arg239* truncation is classified as Pathogenic with multiple-submitter concordance in ClinVar⁵⁵ classified as Pathogenic with multiple-submitter concordance in ClinVar
ClinVar VCV000001782.18, reviewed by GeneDx, Baylor Genetics, and Labcorp Genetics/Invitae with no conflicts. OMIM documents it as allelic variant 609300.0006.

The Evidence

The c.715C>T (Arg239*) variant was identified as the causative mutation in a 2018 case report of a 46,XY adolescent with 17OHD⁶⁶ 2018 case report of a 46,XY adolescent with 17OHD
Zhang et al. 2018, Gynecological Endocrinology, PMID 29345162. The patient developed a rare complication — an adrenal crisis on the first postoperative day after gonadectomy — attributed to insufficient glucocorticoid coverage during surgery, underscoring that 17OHD confers genuine cortisol insufficiency despite the compensatory high corticosterone.

Population-level evidence comes from Willemsen et al. 2025⁷⁷ Willemsen et al. 2025
Meta-analysis of 465 patients across 178 studies, J Clin Endocrinol Metab, PMID 39500362. Across this global cohort, hypertension was present in 57% of patients, hypokalemia in 45%, and primary amenorrhea in 38% of females. Male patients (46,XY) were typically diagnosed earlier because genital dysplasia is apparent at birth or infancy; females were diagnosed later through investigation of amenorrhea and hypertension. Complete loss-of-function variants — including nonsense mutations like Arg239* — were associated with the most severe phenotypes: complete hypocortisolism and absent sexual development.

A literature review of 198 reported 46,XY 17OHD cases⁸⁸ literature review of 198 reported 46,XY 17OHD cases
Kawashima et al. 2025, Endocrine Journal, PMID 40545346 found that 7% of patients with typical female external genitalia developed spontaneous breast tissue, likely from peripheral aromatization of elevated corticosterone. This clinical variability creates diagnostic challenges: 17OHD can be misidentified as complete androgen insensitivity syndrome or primary ovarian insufficiency without genetic testing.

The variant is absent from gnomAD population databases (all ancestries), consistent with strong negative selection against complete CYP17A1 loss-of-function. TOPMed estimates an allele frequency of approximately 1 in 100,000 chromosomes.

Practical Actions

For heterozygous carriers, no hormonal dysfunction occurs — a single functional copy of CYP17A1 provides sufficient enzyme activity. The clinically meaningful action is reproductive genetic counseling: if both partners carry any CYP17A1 loss-of-function variant, each pregnancy has a 25% chance of producing an affected child with complete 17OHD.

For homozygous or compound heterozygous individuals, management follows a two-pronged approach established across multiple case series. First, glucocorticoid replacement — typically hydrocortisone — suppresses ACTH, halts DOC accumulation, and thereby resolves the mineralocorticoid-mediated hypertension and hypokalemia. The hypertension of 17OHD is not essential hypertension; it is driven by DOC excess and resolves when ACTH is adequately suppressed. Starting antihypertensive therapy without glucocorticoid replacement addresses the symptom but not the cause.

Second, sex hormone replacement is added to establish and maintain secondary sexual development appropriate to the individual's gender identity. In 46,XX individuals: estrogen followed by combined estrogen-progestogen replacement. In 46,XY individuals with female phenotype: estrogen replacement after gonadectomy (gonads are typically dysgenetic or undescended and carry a small malignancy risk), with careful perioperative glucocorticoid stress dosing to prevent adrenal crisis.

Interactions

The Arg239* allele most commonly causes disease in the compound heterozygous state paired with a second CYP17A1 pathogenic variant on the other chromosome. Documented partners include Asian founder variants p.H373L, p.W406R, and p.Y329Kfs, as well as other exon 4–6 mutations. Because the Arg239* variant truncates the protein at codon 239 and eliminates all catalytic function, compound heterozygosity with any other loss-of-function allele produces the same complete-deficiency phenotype as homozygosity.

Related variants in the same gene that are tracked in published case series include the nearby missense p.Arg239Gln (rs2439628), which causes incomplete enzyme impairment and a milder phenotype, and the Ser106Pro variant (rs104894135), another complete loss-of-function allele documented in Chinese and Middle Eastern populations.

Your heart contracts and relaxes roughly 100,000 times a day. Each cycle depends on a molecular off-switch: a protein called tropomyosin¹¹ tropomyosin
Tropomyosin is a long, rod-shaped protein that wraps around actin filaments in muscle cells. At rest it physically blocks the sites where myosin (the motor protein) can grab actin, preventing contraction until calcium signals it to step aside that sits like a lid on actin filaments, blocking the cardiac motor machinery until calcium says "go." The TPM1 E180G variant — a single amino-acid substitution replacing glutamic acid with glycine at position 180 of cardiac alpha-tropomyosin — weakens that lid. The result is a heart that can't fully disengage its own throttle, pushing toward hypertrophic cardiomyopathy (HCM): progressive thickening of the left ventricular wall that stiffens the pump and, in some carriers, triggers dangerous arrhythmias.

E180G was first identified in the landmark 1994 Cell paper by Thierfelder and colleagues²² landmark 1994 Cell paper by Thierfelder and colleagues
Thierfelder L et al., Cell 77:701–712, 1994 — the paper that established familial HCM as "a disease of the sarcomere" by showing mutations in multiple distinct sarcomeric proteins produce the same cardiac phenotype, which catalogued the first alpha-tropomyosin mutations causing familial HCM on chromosome 15q2. It is classified as ClinVar Pathogenic (VCV000012455), supported by multiple independent submitters without conflicts, and listed as the first allelic variant in OMIM entry 191010.

The Mechanism

Normally, tropomyosin exists in a tightly coiled double-helix that shifts between three positions — blocked, closed, and open — depending on calcium levels. The Glu180Gly substitution removes a charged glutamic acid residue and replaces it with tiny glycine, which has almost no side chain. This dramatically increases the local and global flexibility³³ local and global flexibility
Measured as persistence length — a quantitative index of a polymer's resistance to bending. Wild-type tropomyosin has a persistence length of approximately 150 nm; E180G reduces this, making the filament more prone to curving and bending of the tropomyosin filament by approximately 35% compared to wild-type protein.

The consequence is mechanical: excess flexibility impedes normal propagation of the "blocked → open" activation signal along the thin filament, and requires a smaller calcium signal to trigger contraction. Atomic force microscopy⁴⁴ Atomic force microscopy
A technique that traces the physical contour of a single protein molecule along a surface, allowing direct measurement of bending at nanometer scale and molecular dynamics simulations⁵⁵ molecular dynamics simulations
Computer models that simulate how individual atoms in a protein move over time, revealing conformational changes too fast to observe experimentally both confirm this increased flexibility, which destabilizes the low-calcium "off" state of cardiac muscle. The result: higher resting actin-myosin interaction, enhanced contractility, and — critically — failure of the muscle to fully relax between beats (diastolic dysfunction).

The Evidence

The molecular evidence for E180G's pathogenic mechanism is strong and consistent across independent methods:

Structural mechanics: Li et al., 2012⁶⁶ Li et al., 2012
Li XE et al. Biochem Biophys Res Commun 2012 — examined both D175N and E180G using electron microscopy of isolated tropomyosin molecules; persistence length reductions were statistically significant across >200 molecules per condition showed increased bending flexibility leads to excess Ca²⁺-activation and shifts regulatory equilibrium toward the "on" state even at diastolic calcium concentrations.

Calcium sensitivity and kinetics: Sewanan et al., 2016⁷⁷ Sewanan et al., 2016
Sewanan LR et al. Front Physiol 2016 — computational myofilament model incorporating E180G-specific stiffness and duty-cycle changes; validated against published in vitro motility data predicted E180G increases both maximum and diastolic force generation, and slows the time from peak tension to 50% relaxation — a signature of HCM diastolic dysfunction.

Signaling cascades: Robinson et al., 2018⁸⁸ Robinson et al., 2018
Robinson P et al. J Biol Chem 2018 — guinea pig cardiomyocyte model expressing thin-filament HCM mutations including the closely related D175N demonstrated that increased myofilament calcium buffering from these mutations elevates diastolic calcium and activates CaMKII and calcineurin/NFAT⁹⁹ CaMKII and calcineurin/NFAT
Two calcium-sensitive kinases that, when chronically activated, trigger gene programs causing hypertrophy, fibrosis, and arrhythmia — the cardinal features of HCM signaling cascades — providing a mechanism linking the sarcomeric defect to macroscopic cardiac remodeling.

Actin-myosin interaction: Kopylova et al., 2019¹⁰¹⁰ Kopylova et al., 2019
Kopylova GV et al. J Muscle Res Cell Motil 2019 — combined single-molecule optical trap and ensemble in vitro motility assay; both HCM mutations E180G and D175N increased calcium sensitivity in ensemble measurements, whereas dilated cardiomyopathy mutations showed opposite effects confirmed E180G increases thin-filament sliding velocity and calcium sensitivity in reconstituted assays, in the same direction as D175N — and in the opposite direction from dilated cardiomyopathy TPM1 mutations, validating the disease-specific mechanism.

E180G is ultra-rare globally (absent from gnomAD population databases), consistent with its being a disease-causing variant with strong negative selection pressure. The variant is classified Pathogenic with a 3-star ClinVar review status reflecting consistent classification across multiple independent submitters.

Practical Actions

Identifying an E180G carrier changes clinical management in ways that directly reduce morbidity and mortality. HCM is one of the most common causes of sudden cardiac death in individuals under 35. For carriers, the priorities are: (1) confirm diagnosis and baseline LV morphology, (2) identify high-risk features that warrant ICD implantation or septal reduction therapy, and (3) extend family testing to first-degree relatives. The autosomal dominant inheritance means each biological child of a carrier has a 50% chance of inheriting the variant.

Avoidance of competitive athletics and extreme exertion is recommended pending full clinical evaluation — sudden cardiac death events in HCM are disproportionately exercise-associated.

Interactions

E180G's closely related neighbor on chromosome 15, TPM1 D175N (rs104894503), shares nearly identical functional properties: both increase tropomyosin flexibility and calcium sensitivity, both were identified in the same 1994 paper, and both are classified Pathogenic for HCM. D175N is a Finnish founder mutation (accounting for 6.5% of Finnish HCM cases in a cohort of 306 patients) but is not known to interact with E180G as a compound heterozygote — these are independent dominant mutations in the same gene affecting adjacent codons.

Other sarcomeric protein HCM genes — MYBPC3 (rs36211723), MYH7, TNNT2 — can produce overlapping clinical phenotypes. Patients with multiple sarcomeric variants ("double positive" HCM) tend to have more severe hypertrophy and earlier onset, though compound heterozygosity specifically for E180G has not been studied in published literature.

The NOD2 gene encodes a cytosolic pattern recognition receptor¹¹ pattern recognition receptor
NOD2 senses muramyl dipeptide (MDP), a fragment of bacterial cell-wall peptidoglycan that detects bacterial cell-wall fragments and mounts an immune defense in the intestinal lining. It was the first gene identified for Crohn's disease²² the first gene identified for Crohn's disease
NOD2 was discovered as a susceptibility locus in the IBD1 linkage region on chromosome 16q12 susceptibility and remains the most studied single gene in inflammatory bowel disease genetics. The rs104895444 variant changes valine to methionine at protein position 793 (Val793Met), a site within the leucine-rich repeat domain that participates in bacterial ligand sensing.

The Mechanism

NOD2 protein consists of two N-terminal CARD (caspase-recruitment) domains³³ CARD (caspase-recruitment) domains
CARDs mediate protein-protein interactions that activate the NF-κB signaling cascade, a central NOD domain, and C-terminal leucine-rich repeats (LRRs). The LRR domain directly recognizes muramyl dipeptide (MDP), the minimal bioactive peptidoglycan fragment from both Gram-positive and Gram-negative bacteria. Successful MDP binding triggers NOD2 oligomerization, RIPK2 recruitment, and downstream NF-κB activation⁴⁴ NF-κB activation
NF-κB activates genes encoding antimicrobial peptides (defensins), pro-inflammatory cytokines, and barrier-strengthening proteins — a coordinated response that clears bacteria while reinforcing the epithelial barrier.

The Val793Met change lies within the LRR domain at a position conserved across mammals. While ClinVar currently classifies this variant as likely benign based on existing submissions, it was identified in deep resequencing⁵⁵ identified in deep resequencing
Rivas et al. 2011 whole-exome resequencing of 16,054 CD and 12,153 UC cases identified V793M among additional rare NOD2 risk variants beyond the three classic mutations as one of several rare NOD2 missense variants independently associated with IBD in a genome-wide study of over 45,000 individuals. The substitution of valine (a small non-polar amino acid) for methionine (a larger sulfur-containing amino acid) could alter LRR domain folding or MDP-binding geometry, potentially reducing the sensitivity of bacterial sensing.

The Evidence

The genetic case for Val793Met rests on its identification in a landmark deep-resequencing study⁶⁶ landmark deep-resequencing study
Rivas et al. sequenced 163 genes in 9 GWAS loci in 16,054 CD cases, 12,153 UC cases, and 17,575 controls, identifying V793M among four additional independent NOD2 risk factors that interrogated over 45,000 IBD cases and controls. The study identified Val793Met (V793M) as one of four additional independent NOD2 risk factors beyond the three classic variants (R702W, G908R, L1007fs), each of which carries odds ratios of 2–4 for heterozygotes and 17–40 for compound heterozygotes or homozygotes.

The broader NOD2-Crohn's connection is established beyond any doubt⁷⁷ established beyond any doubt
More than 145 published studies confirm NOD2 as the strongest non-HLA locus for Crohn's disease. The three classic variants together account for ~80% of NOD2-attributable CD risk in European populations; 27 rare mutations account for the remaining ~20%⁸⁸ 27 rare mutations account for the remaining ~20%
Lesage et al. 2002 found that in 612 IBD patients, the three major variants contributed 81% of disease-associated alleles while 27 rare variants contributed 19%. Val793Met occupies this second tier — rare individually, but clinically meaningful in the context of compound heterozygosity.

The evidence for Val793Met as a standalone risk factor is moderate rather than established: it is supported by one large resequencing study but lacks the replication depth of the three major NOD2 variants. Its individual contribution to Crohn's disease risk is likely modest. The variant's significance is amplified when a second NOD2 hit is present on the opposite chromosome.

Practical Implications

Carriers of one copy of the A allele face a mildly elevated risk for Crohn's disease, primarily through two mechanisms: direct LRR domain disruption and, more importantly, the possibility of compound heterozygosity with another NOD2 variant. Because Val793Met is rare (~0.2% allele frequency in Europeans), most carriers will not have a second NOD2 hit — but those who do face substantially elevated risk.

NOD2-associated Crohn's disease has a characteristic phenotype: ileal or ileocolonic location, fibrostenotic (stricturing) complications, and a tendency toward surgical resection. Patients with two NOD2 mutations⁹⁹ Patients with two NOD2 mutations
Double-dose NOD2 carriers show earlier onset, 53% stricturing rate vs 28%, and more frequent ileal-only disease show earlier disease onset, more stricturing behavior, and higher surgical rates than single-copy carriers or non-carriers.

NOD2 variants are also being used as biomarkers to guide therapy¹⁰¹⁰ biomarkers to guide therapy
NOD2 genotype informs risk of steroid refractoriness, surgical necessity, and consideration of early biologic use decisions in Crohn's disease management. Carriers with established IBD may benefit from early specialist involvement to monitor for fibrostenotic complications.

Interactions

The clinical significance of Val793Met is substantially magnified by compound heterozygosity. When this variant occurs on one chromosome alongside a different NOD2 variant — such as R702W (rs2066844), G908R (rs2066845), or L1007fs (rs2066847) — on the opposite chromosome, both alleles carry impaired function and the net NOD2 activity is severely reduced. Compound heterozygosity for two NOD2 mutations¹¹¹¹ Compound heterozygosity for two NOD2 mutations
In the large NOD2 genotype-phenotype study, double-mutation carriers showed OR ~40 for CD compared to OR ~3-4 for single-copy carriers dramatically amplifies Crohn's disease risk compared to either mutation alone.

The rs104895431 variant in the same batch belongs to the same NOD2 compound heterozygosity risk cluster. Individuals carrying both rs104895444 and rs104895431 on opposite chromosomes would have effectively zero functional NOD2 alleles, creating a high-risk state for Crohn's disease with the ileal-stricturing phenotype typical of double-NOD2-mutation carriers. This interaction is the central clinical significance of both variants — individually rare and modest, combined potentially severe.

NOD2 genotype also interacts with smoking status: NOD2 mutation carriers who smoke¹²¹² NOD2 mutation carriers who smoke
Kuenzig et al. 2017 found NOD2-smoking interaction modified stricturing phenotype risk in CD show modified disease outcomes, with smoking amplifying the stricturing tendency already present with NOD2 variants. The pathobiological pathway converges on the gut microbiome — NOD2-deficient mice develop dysbiosis and intestinal inflammation that mirrors CD, confirming the bacterial-sensing deficit hypothesis¹³¹³ confirming the bacterial-sensing deficit hypothesis
NOD2-knockout mouse models develop altered gut microbiome composition and show heightened intestinal inflammation.

Every cell in your body runs a continuous protein quality-control operation called the ubiquitin-proteasome system (UPS)¹¹ ubiquitin-proteasome system (UPS)
The UPS tags damaged or misfolded proteins with ubiquitin chains and feeds them to the proteasome — a barrel-shaped multi-subunit protease — for destruction. It also controls the half-lives of regulatory proteins, including the inhibitors of inflammatory signalling. The 20S core proteasome is assembled from alpha and beta subunits; PSMA6 encodes the alpha-6 subunit, which forms part of the outer rings that gate substrate entry. The rs1048990 variant sits just eight bases upstream of the translation start in the 5'UTR of PSMA6 exon 1 — a location that has a quiet but measurable effect on how much PSMA6 protein the cell produces.

The G allele enhances PSMA6 transcription²² enhances PSMA6 transcription
Demonstrated by reporter assays and confirmed by the correlation between genotype and mRNA expression levels across multiple population cohorts. More PSMA6 means more assembled proteasomes, and more proteasome activity accelerates the degradation of IκB — the inhibitor of NF-κB³³ IκB — the inhibitor of NF-κB
IκB normally sequesters NF-κB in the cytoplasm, preventing it from driving expression of inflammatory cytokines (TNF-α, IL-1β, IL-6). When IκB is phosphorylated after an immune stimulus, it is tagged for proteasomal degradation. More proteasomes degrade phospho-IκB faster, releasing NF-κB to the nucleus more readily. The result is a lower threshold for inflammatory gene expression — not a dramatic immune deficiency, but a subtle amplification of the inflammatory signal across the lifetime.

The Mechanism

The -8C>G substitution does not change any protein but alters the regulatory architecture of the PSMA6 5'UTR. Functional studies demonstrated that PSMA6 knockdown using siRNA reduced NF-κB activation by stabilising phosphorylated IκB in cultured cells⁴⁴ Functional studies demonstrated that PSMA6 knockdown using siRNA reduced NF-κB activation by stabilising phosphorylated IκB in cultured cells
The experiment directly connects proteasome abundance to NF-κB activity: less proteasome → slower IκB degradation → more IκB retained → less NF-κB in the nucleus → less inflammatory gene expression. Conversely, the G allele's enhanced transcription runs this pathway in the opposite direction — faster IκB turnover, lower threshold for NF-κB activation, and greater cytokine output in response to inflammatory stimuli.

This mechanism is relevant across multiple disease contexts: chronic low-grade inflammation drives atherosclerosis and plaque instability in cardiovascular disease, and overactive NF-κB signalling is a central feature of autoimmune conditions including rheumatoid arthritis, SLE, and psoriasis.

The Evidence

The original discovery study⁵⁵ The original discovery study
Ozaki K et al. A functional SNP in PSMA6 confers risk of myocardial infarction in the Japanese population. Nat Genet. 2006 found the G allele in 2,592 MI cases and 2,851 controls from Japan, with striking significance (chi²=21.1, p=4.4×10⁻⁶) and a minor allele frequency of ~0.35 in this population — much higher than in Europeans (~0.17). Replication was confirmed in an independent Japanese cohort.

A meta-analysis pooling 15,991 cases and 16,784 controls⁶⁶ meta-analysis pooling 15,991 cases and 16,784 controls
Wang M et al. Quantitative assessment of the influence of PSMA6 variant (rs1048990) on coronary artery disease risk. Mol Biol Rep. 2013 found a per-G-allele OR of 1.09 (95% CI 1.02–1.16, p=0.006) globally, with the GG homozygous recessive model yielding OR 1.38 (1.22–1.57). The effect was significant in East Asians but did not reach significance in Europeans when analysed separately — likely because GG homozygotes are far rarer in European populations (~2% vs ~9% in Japanese), reducing statistical power.

A UK replication study in 6,946 MI cases and 2,720 controls⁷⁷ A UK replication study in 6,946 MI cases and 2,720 controls
Balmforth AJ et al. The exon 1-8C/G SNP in the PSMA6 gene contributes only a small amount to the burden of myocardial infarction. Atherosclerosis. 2008 confirmed the GG homozygote frequency is far lower in Europeans (2.1% vs 8.9% in Japanese) and that the individual-study result was non-significant, but the combined meta-analysis consistently returned OR ~1.15 per G allele and OR ~1.38 for GG in the recessive model.

Beyond cardiovascular disease, rs1048990 has been associated with type 1 diabetes⁸⁸ type 1 diabetes
OR 2.04, 95% CI 1.38–3.03, p<0.001 in a Latvian cohort, with expression correlates across 14 UPS-related genes, juvenile idiopathic arthritis⁹⁹ juvenile idiopathic arthritis
Polish children study with PSMA6 rs1048990 significant for JIA specifically, and psoriasis¹⁰¹⁰ psoriasis
PSMA6 -8C>G among significant variants in Polish adult psoriasis patients. A male-specific association with asthma was also found in Lithuanian patients, consistent with sex-differential NF-κB signalling.

The breadth of associations across inflammatory and autoimmune phenotypes is expected from a variant that modulates a master inflammatory transcription factor, though the effect size at any single locus is modest (OR 1.1–1.4), fitting a complex polygenic architecture.

Practical Implications

The G allele's elevated NF-κB tone is not manageable through generic anti-inflammatory strategies (these are excluded per platform rules), but there are genotype-specific consequences worth acting on. GG homozygotes carry the highest cardiovascular inflammatory burden. CG heterozygotes carry intermediate risk. The most actionable areas are cardiovascular biomarker monitoring — since the mechanism is inflammatory plaque instability — and awareness of the autoimmune associations for early detection.

Interactions

PSMA6 encodes one of 14 alpha-subunit isoforms of the 20S proteasome. Other proteasome subunit genes have been studied alongside rs1048990: rs2277460 in PSMA6 (a second variant in the same gene), rs2295826 and rs2295827 in PSMC6 (a 19S regulatory particle ATPase), and rs2348071 in PSMA3. Haplotypes spanning PSMA6 and the neighbouring KIAA0391 gene cluster on chromosome 14q13 show stronger associations with coronary artery disease than rs1048990 alone. Individuals carrying risk alleles across multiple proteasome subunit loci may have compounded NF-κB dysregulation, though published compound heterozygosity data remain limited.