FUT6 rs78060698 — The Fucosylation Switch for Vitamin B12

Deep in the cells lining your small intestine, a family of enzymes called fucosyltransferases are quietly sculpting the sugar coat on your cell surfaces. These glycan structures — chains of modified sugars attached to proteins — determine which microbes can colonize your gut, how nutrients move across the intestinal wall, and how effectively your body captures vitamin B12. The FUT6 gene encodes one of these enzymes, and a single variant in its regulatory region can meaningfully shift circulating B12 levels. This variant, rs78060698, was identified in a 2017 genome-wide association study of Indian adults — a population where B12 deficiency affects an estimated 47–70% of people — and represents one of the clearest examples of how gut biology, microbial ecology, and nutrition intersect at the genetic level.

The Mechanism

FUT6 encodes alpha-1,3-fucosyltransferase 611 alpha-1,3-fucosyltransferase 6
An enzyme that transfers fucose — a six-carbon sugar — onto glycan chains on cell surfaces, creating Lewis X and sialyl-Lewis X antigens that mediate cell-cell and host-microbe recognition
. These fucosylated glycans on the intestinal epithelium act as molecular docking sites for gut bacteria and influence the local microbial ecology of the small intestine.

The rs78060698 variant sits within an intron of FUT6, not in the protein-coding sequence itself. Despite its intronic location, it has clear regulatory function. Luciferase reporter assays using human HepG2 liver cells demonstrated that the A allele produces approximately 3× higher FUT6 promoter activity and 3.5–20× higher enhancer activity compared to the G allele. Electrophoretic mobility shift assays confirmed that this difference arises from differential binding of HNF4α22 HNF4α
Hepatocyte Nuclear Factor 4-alpha — a transcription factor that regulates many genes involved in glucose, lipid, and vitamin metabolism, and is a master regulator of fucosyltransferase expression
: the A allele binds HNF4α with ~1.18-fold greater affinity.

The proposed pathway: higher FUT6 expression → more fucosylated glycans on intestinal epithelium → altered composition of gut microbiota → changes in bacterial production or competition for vitamin B12. Unlike its close relative FUT2, whose effects on B12 appear to operate through secretor status and H. pylori susceptibility, FUT6 genotype is associated with B12 levels independently of secretor status and H. pylori antibody titers — suggesting a distinct microbial or absorptive mechanism.

The Evidence

The primary evidence comes from a 2017 GWAS in 4,419 Indians33 2017 GWAS in 4,419 Indians
Nongmaithem SS et al. GWAS identifies population-specific new regulatory variants in FUT6 associated with plasma B12 concentrations in Indians. Human Molecular Genetics, 2017
. The study combined a discovery cohort from the Pune Maternal Nutrition Study with three independent Indian replication cohorts. The rs78060698 A allele was associated with higher plasma B12 (beta = 0.22 on log scale, P = 8.3×10⁻¹⁷), with consistent effects across age groups and pregnancy status.

A critical population-frequency difference shapes the clinical relevance: the A allele frequency was 0.21 in Indians versus only 0.03 in Europeans (CEU panel, 1000 Genomes). This 7-fold enrichment means the variant explains substantially more B12 variance in South Asian populations than in European ones, and was likely not detected in earlier European GWAS because of its low frequency. In silico analysis confirmed the variant's functional prediction scores were significant across populations, but population-specific LD structure and effect size differences mean extrapolation to non-Indian populations requires caution.

Partial linkage disequilibrium (r² ≈ 0.54 in Indians) with a second independent FUT6 variant, rs3760775, suggests the two SNPs tag distinct but correlated regulatory signals in the same chromosomal region. Conditional analysis in the primary study confirmed rs78060698 retains independent association after adjusting for rs3760775.

Evidence is rated moderate: the association is highly significant and biologically supported by functional assays, but the causal mechanism remains proposed rather than experimentally confirmed in vivo, replication in non-Indian populations is limited, and no clinical intervention trials exist.

Practical Actions

The actionable implication of this variant is about baseline B12 monitoring and optimizing intake to compensate for genetic variation in absorptive capacity. Those with GG genotype carry no copies of the B12-boosting A allele and may have meaningfully lower circulating B12 than AG or AA counterparts — a difference that compounds with dietary insufficiency (vegetarian or vegan diets, low dairy intake) and age-related declines in gastric acid that impair B12 absorption from food.

Monitoring serum B12 — and specifically holotranscobalamin (active B12) when available — is the most direct way to determine whether genetically lower absorptive capacity translates to functional deficiency. For supplementation, methylcobalamin and adenosylcobalamin are the bioactive forms; sublingual methylcobalamin bypasses intestinal absorption steps entirely and is particularly useful when GI function is compromised.

Interactions

rs78060698 sits in the same gene cluster as rs3760775 (FUT6), which shows a slightly stronger B12 association (beta = 0.25, P = 1.2×10⁻²³) and is partially correlated (r² = 0.54 in Indians). The two variants likely tag overlapping but non-identical regulatory elements; individuals carrying both effect alleles may experience additive benefits to B12 status.

The FUT2 variants rs601338 and rs602662 operate on a related but distinct pathway (secretor status → holo-haptocorrin glycosylation → H. pylori susceptibility). Because FUT6 genotype is independent of secretor status, carrying GG at rs78060698 alongside a non-secretor FUT2 genotype represents two separate mechanisms converging on lower B12 — a combination worth tracking with serum monitoring.

rs9594759

TNFSF11

Strong Risk Factor

RANKL Regulatory Variant — Guardian of Bone Remodeling Balance

Your TNFSF11 gene encodes RANKL (receptor activator of nuclear factor kappa-B ligand11 receptor activator of nuclear factor kappa-B ligand
a master regulator of bone remodeling
), a cytokine that tells your body when to break down old bone through osteoclast activation. This particular variant lies in a regulatory region upstream of the RANKL gene22 regulatory region upstream of the RANKL gene
about 184 kb upstream, in an area that modulates gene expression
and influences how much RANKL your bone cells produce. Too much RANKL activity tips the balance toward bone loss; too little prevents normal bone turnover. Getting this balance right is essential for maintaining bone strength throughout life, especially as you age.

The Mechanism

This SNP sits in a regulatory enhancer region33 regulatory enhancer region
a DNA sequence that controls gene expression from a distance
that responds to vitamin D and parathyroid hormone signals44 vitamin D and parathyroid hormone signals
1,25-dihydroxyvitamin D3 and PTH bind to vitamin D receptor (VDR) and CREB at this enhancer
. The T allele appears to alter the binding efficiency of these regulatory factors55 alter the binding efficiency of these regulatory factors
functional experiments show differential promoter inhibition
, potentially leading to increased RANKL expression in bone tissue. When RANKL levels rise, more osteoclasts differentiate and activate66 osteoclasts differentiate and activate
through RANK-RANKL signaling and downstream NF-κB activation
, accelerating the breakdown of bone matrix. Over time, this shifts the bone remodeling equilibrium toward net bone loss, particularly in contexts where other factors (low dietary calcium, vitamin D deficiency, hormonal changes) also promote resorption.

The Evidence

A validation study in 700 elderly Chinese subjects77 A validation study in 700 elderly Chinese subjects
350 with hip osteoporotic fractures, 350 controls
found significant association between TNFSF11 variants including rs9594759 and hip fracture risk (p=0.018). T allele carriers showed lower bone mineral density88 lower bone mineral density
particularly at the lumbar spine
in multiple cohort studies. Genome-wide association studies99 Genome-wide association studies
including the landmark 2008 GWAS
have consistently identified the TNFSF11 region at chromosome 13q14 as one of the most robust loci associated with bone mineral density variation and osteoporotic fracture risk.

The functional relevance was confirmed through enhancer deletion studies in mice1010 enhancer deletion studies in mice
deletion of RL-D2 enhancer led to high bone mass phenotype
, which demonstrated that regulatory variants in this region directly control RANKL expression and bone remodeling rates. Importantly, this regulatory region responds to vitamin D1111 this regulatory region responds to vitamin D
inhibition significantly reduced in presence of vitamin D
, suggesting that adequate vitamin D status may partially compensate for genetic risk.

Practical Implications

If you carry the T allele, your bone cells may produce more RANKL in response to normal physiological signals, increasing your baseline rate of bone turnover. This becomes particularly important after age 50, during menopause (when estrogen loss further elevates RANKL), or if your diet is low in calcium. The good news: bone health is highly modifiable through nutrition and lifestyle. Adequate calcium and vitamin D intake1212 Adequate calcium and vitamin D intake
shown to reduce RANKL levels and bone loss
can help offset genetic predisposition. Weight-bearing exercise stimulates bone formation and may help maintain the remodeling balance. Regular bone density screening becomes more important if you have two copies of the T allele, as early detection allows for targeted interventions before fractures occur.

Interactions

This variant interacts with other genes in the RANK/RANKL/OPG pathway1313 RANK/RANKL/OPG pathway
the trio that regulates bone remodeling
, including TNFRSF11A (RANK receptor) and TNFRSF11B (osteoprotegerin). Variants in the vitamin D receptor (VDR) gene also modulate risk, as VDR polymorphisms affect how bone cells respond to vitamin D1414 VDR polymorphisms affect how bone cells respond to vitamin D
combined VDR and TNFSF11 variants show gene-gene interactions
. Additionally, calcium intake directly influences RANKL expression1515 calcium intake directly influences RANKL expression
low calcium triggers secondary hyperparathyroidism and RANKL upregulation
, meaning dietary habits interact with this genetic variant to determine actual bone health outcomes.

rs2069705

IFNG IFNG Promoter -1616C/T

Moderate Risk Factor

IFNG Promoter -1616C/T — The Th1 Output Dial and Atopic Susceptibility

IFN-gamma11 IFN-gamma
interferon-gamma (IFNG) — the master cytokine of Th1 immunity, produced primarily by CD4+ Th1 cells, CD8+ cytotoxic T cells, and NK cells. It activates macrophages to kill intracellular pathogens, promotes Th1 differentiation, and critically suppresses Th2 cytokines (IL-4, IL-5, IL-13) that drive allergic inflammation
is the cornerstone of type 1 adaptive immunity. When IFN-gamma production is robust, T cells preferentially differentiate into the Th1 lineage, suppressing the IgE-producing, eosinophil-recruiting Th2 responses that underlie atopic dermatitis, allergic rhinitis, and asthma. When IFN-gamma output is reduced, the Th1/Th2 balance tips toward Th2 dominance — the molecular underpinning of atopic disease. rs2069705 sits approximately 1,616 base pairs upstream of the IFNG transcription start site (GRCh38 chr12:68161231), in the promoter region that governs how strongly the gene is switched on in response to immune stimuli.

The Mechanism

rs2069705 is annotated as a regulatory variant in the IFNG upstream promoter region. The IFNG gene sits on the minus strand of chromosome 12; papers using coding-strand notation call this position -1616C/T (where C on the coding strand corresponds to the G GRCh38 reference allele on the plus strand, and T corresponds to the A alternate allele). The A allele (coding T) creates or strengthens a binding site for STAT422 STAT4
Signal Transducer and Activator of Transcription 4 — a transcription factor activated by IL-12 and IL-18 signaling that drives Th1 differentiation and promotes IFN-gamma transcription
. A 2024 functional study by Chen et al. demonstrated that rs2069705 "boosts IFNγ transcription by promoting interaction between its promoter and STAT4," activating the downstream JAK/STAT1 pathway. The G allele (coding C) lacks this STAT4 binding enhancement, resulting in lower baseline IFN-gamma transcriptional output — translating to a reduced Th1 tone and a permissive environment for Th2-mediated allergic inflammation.

This is part of a broader circuit: T-bet (TBX21)33 T-bet (TBX21)
the master Th1 transcription factor that drives IFN-gamma expression; T-bet promotes IFNG transcription and represses GATA-3-driven Th2 differentiation
drives IFN-gamma production, and rs2069705 modulates how responsive the IFNG promoter is to T-bet-upstream signals including STAT4. Disruption at either node — reduced T-bet expression (rs4794067, TBX21 promoter variant) or reduced IFNG promoter responsiveness (rs2069705) — can depress Th1 output and increase Th2-driven atopic susceptibility.

The Evidence

The most direct evidence for the G allele's reduced IFN-gamma output comes from its association pattern across immune diseases. A 2023 Russian pediatric asthma study44 2023 Russian pediatric asthma study
Smolnikova et al. Vavilovskii Zhurnal Genet Selektsii 2023; 263 Russian children with asthma
found that the coding TT genotype (plus-strand AA) was specifically associated with mild and controlled asthma phenotypes (p<0.05) — the Th1-competent A allele homozygotes showing the most favourable disease trajectory.

The clearest functional evidence comes from a 2024 mechanistic study in primary Sjögren's syndrome55 2024 mechanistic study in primary Sjögren's syndrome
Chen et al. Am J Physiol Cell Physiol 2024; luciferase reporter assays and chromatin immunoprecipitation demonstrating STAT4-IFNG promoter interaction
which established the A allele as a transcriptional activator at this position.

At the related downstream SNP rs2430561 (+874T/A)66 rs2430561 (+874T/A)
A nearby IFNG variant in partial LD with rs2069705, located in an NF-kappaB binding site at position +874 relative to the IFNG TSS
, a 2009 Egyptian atopy study77 2009 Egyptian atopy study
Hussein et al. J Investig Allergol Clin Immunol 2009; Egyptian atopic patients vs healthy controls
directly linked the low-producer IFNG allele to atopic disease: atopic patients showed a significantly higher frequency of the A allele at +874 (the low-producer allele at that position), with AA homozygotes showing decreased serum IFN-gamma, elevated total IgE, and increased eosinophil counts compared to TT homozygotes — a direct quantitative readout of how reduced IFN-gamma production enables Th2 immune activation.

At the gene interaction level, a 2016 SLE study88 2016 SLE study
Leng et al. Sci Rep 2016; 3,732 Chinese Han subjects
identified a significant genetic interaction between rs2069705 (IFNG) and rs4794067 (TBX21 promoter) in SLE susceptibility — illustrating that this IFNG promoter variant does not act in isolation but as part of a coordinated T-bet/IFN-gamma regulatory axis.

Practical Implications

For GG homozygotes (coding CC, ~23% globally, ~11% of Europeans), the reduced IFNG promoter responsiveness creates the lowest baseline IFN-gamma output of the three genotypes. In the allergy-atopic context this translates to the most permissive environment for Th2-driven sensitisation: higher IgE class-switching potential, less Th1-mediated suppression of eosinophil recruitment, and increased susceptibility to sensitisation to environmental allergens. Strategies that actively support Th1 immune balance are most relevant for this group.

For AG heterozygotes (~50% globally), one G allele moderately reduces IFNG promoter responsiveness. The practical implications are milder, but the atopic tendency is real, and early attention to Th1-supporting exposures during immune development is appropriate.

AA homozygotes carry the population-common, Th1-competent genotype in most European and South Asian populations and are not at elevated atopic risk from this variant.

Interactions

The most clinically relevant interaction is between rs2069705 and rs4794067 (TBX21 promoter). TBX21 encodes T-bet, which drives IFN-gamma transcription; rs4794067-C reduces T-bet expression. Leng et al. (2016) found that while neither variant independently reached significance for SLE in their cohort, the combination was significant — suggesting that disruption at both the upstream driver (T-bet) and the downstream promoter (IFNG) produces compounded immune dysregulation exceeding either variant's individual effect. For the allergy-atopic category, carriers of both GG at rs2069705 and CC at rs4794067 face reduced IFN-gamma through two parallel mechanisms.

The nearby IFNG variant rs243056199 rs2430561
IFNG +874T/A, an NF-kappaB binding site variant extensively studied in atopic disease and autoimmunity
is in partial LD with rs2069705 and likely captures overlapping variance in IFNG promoter activity.

rs2236225

MTHFD1 G1958A

Strong Risk Factor

MTHFD1 — The Choline Connection

MTHFD1 (methylenetetrahydrofolate dehydrogenase 1) is a trifunctional enzyme that processes dietary folates through three sequential reactions. It plays a central role in one-carbon metabolism 11 One-carbon metabolism: a network of folate-dependent reactions that shuttle single carbon units for DNA synthesis and methylation, feeding into both nucleotide synthesis 22 For DNA repair and cell division — rapidly dividing cells like gut lining and blood cells are especially dependent and the methylation cycle.

The Mechanism

The G1958A variant (rs2236225) causes an arginine-to-glutamine substitution 33 Arginine-to-glutamine substitution at position 653 of the protein (p.Arg653Gln) at position 653 of the MTHFD1 protein. The A allele produces a less thermostable enzyme that loses activity more readily at body temperature. This reduces the efficiency of folate processing, particularly the 10-formyltetrahydrofolate synthetase activity that is important for purine synthesis. While the enzyme retains normal substrate affinity, its reduced stability diminishes overall metabolic activity.

The Choline Compensation

What makes MTHFD1 especially interesting is its connection to choline. When MTHFD1 activity is reduced, your body compensates by drawing more heavily on choline as an alternative methyl donor 44 The betaine pathway: choline is oxidized to betaine, which donates a methyl group directly to homocysteine, bypassing the folate cycle. This increases your dietary choline requirements significantly. Studies have shown that individuals with the AA genotype who consume low-choline diets are more likely to develop signs of choline deficiency, including fatty liver.

The Evidence

A landmark study by Kohlmeier et al.55 Kohlmeier et al.
Kohlmeier M et al. PNAS 2005 — genetic variation in folate-mediated one-carbon transfer predicts susceptibility to choline deficiency
demonstrated that the A allele is a risk factor for neural tube defects, independent of MTHFR status. A meta-analysis of nine studies66 meta-analysis of nine studies
Shen H et al. MTHFD1 polymorphisms and neural tube defect susceptibility, 2014
with 4,302 NTD patients and 4,238 controls confirmed an increased risk of neural tube defects with the AA genotype (OR=2.63). Subsequent research confirmed that this variant increases choline requirements and that adequate choline intake can compensate for the reduced MTHFD1 activity.

Practical Implications

Egg yolks are the richest common dietary source of choline, providing about 150mg per yolk. Liver is even richer. If you carry the A allele, eating 2-3 egg yolks daily provides meaningful choline support. This is one of the most actionable nutrigenomics findings — a simple dietary change (eating more eggs) can compensate for a clear genetic limitation.

Interactions

MTHFD1 interacts with MTHFR (rs1801133, rs1801131) for overall folate pathway efficiency. It also interacts with PEMT (rs7946) — both variants increase choline requirements, and the combined effect can be substantial.

DDX39B — The RNA Helicase Gating Autoimmune Risk

Deep in the MHC class III region on chromosome 6, the gene DDX39B (also known as BAT1) encodes an RNA helicase with a central role in mRNA export, pre-mRNA splicing, and nuclear export of immune transcripts. A 2017 landmark study in Cell established that a regulatory variant at rs2523506 reduces DDX39B protein production, and that this reduction propagates through a chain of molecular events to increase the risk of multiple sclerosis (MS) — the first rigorous demonstration of biological epistasis in humans. The discovery matters not just for MS genetics but as a proof of concept: two common variants in separate genes, each with modest individual effects, combine to produce a dramatically elevated risk that neither alone explains.

The Mechanism

DDX39B is a DEAD-box RNA helicase11 DEAD-box RNA helicase
A family of enzymes that use ATP hydrolysis to unwind RNA secondary structures and remodel RNA-protein complexes, enabling downstream RNA processing steps
. One of its critical functions is promoting the inclusion of IL7R exon 622 IL7R exon 6
Exon 6 of the IL7R gene encodes a transmembrane anchor; when included, IL7R is expressed on the T cell surface as a membrane-bound receptor. When skipped, the resulting mRNA produces a soluble, secreted form that cannot signal properly
during pre-mRNA splicing. When exon 6 is properly included, the interleukin-7 receptor (IL7R) anchors to the T cell surface and signals for T cell survival and differentiation. When exon 6 is skipped, the mRNA encodes a soluble IL7R (sIL7R)33 soluble IL7R (sIL7R)
The secreted form acts as a decoy receptor, sequestering IL-7 away from surface IL7R and dysregulating T cell homeostasis
isoform that is secreted rather than membrane-bound, dysregulating IL-7 signaling in ways that predispose to autoimmunity.

The rs2523506 T allele (reported as "A" on the coding/minus strand in the original paper) lies in the 5' untranslated region of DDX39B. This regulatory change reduces the translational efficiency44 translational efficiency
How efficiently a cell's ribosomes convert DDX39B mRNA into DDX39B protein; a less efficient 5' UTR means fewer protein molecules are produced from the same amount of mRNA
of DDX39B mRNA — cells with the T allele produce less DDX39B protein. With less DDX39B helicase available, the spliceosome is less able to promote IL7R exon 6 inclusion, and exon 6 skipping increases. This cascade — T allele → less DDX39B protein → more exon 6 skipping → more soluble IL7R → impaired IL-7 signaling → dysregulated T cell homeostasis — provides the mechanistic basis for the variant's MS association.

Subsequent research has broadened the picture considerably. A 2023 eLife study showed that DDX39B is also essential for proper splicing of FOXP3, the master transcription factor of regulatory T cells (Tregs). When DDX39B is depleted, FOXP3 introns are retained, FOXP3 protein is lost, and Treg function collapses — providing a second immune tolerance mechanism under DDX39B control. A 2024 mechanistic study confirmed that DDX39B's ATPase activity (not its helicase activity per se) is required for efficient pre-spliceosome assembly at FOXP3 introns, and demonstrated that MS susceptibility genes are disproportionately enriched among genes affected by DDX39B depletion (p = 0.00013). DDX39B thus emerges as a broad guardian of immune gene splicing, with the T allele at rs2523506 chronically reducing this protection.

The Evidence

Galarza-Muñoz et al. 201755 Galarza-Muñoz et al. 2017
Human Epistatic Interaction Controls IL7R Splicing and Increases Multiple Sclerosis Risk. Cell 169(1):72-84
is the primary study. The team conducted genetic association analysis, reporter assays for translational efficiency, and splicing experiments in primary CD4+ T cells and lymphoblastoid cell lines. Key findings: (1) the DDX39B 5' UTR T allele reduces translation in reporter assays, (2) DDX39B depletion causes increased IL7R exon 6 skipping preferentially in the context of the IL7R risk allele (rs6897932 C allele), and (3) carriers of both risk alleles — the DDX39B T allele and the IL7R C allele — show a combined OR of approximately 2.75 for MS, substantially exceeding what either variant contributes alone.

The epistatic architecture is critical to understand. The IL7R rs6897932 C allele reduces exon 6 splicing efficiency intrinsically (the exon splice site is weaker). The DDX39B T allele reduces the level of the helicase that compensates for this weakness. Together, they remove both the intrinsic and compensatory mechanisms for exon 6 inclusion — a synthetic depletion that explains why the combined genotype confers dramatically elevated risk whereas either alone produces more modest effects.

The 2023 FOXP3 finding (Hirano et al. 202366 Hirano et al. 2023
The RNA helicase DDX39B activates FOXP3 RNA splicing to control T regulatory cell fate. eLife 12
) extends this model: MS susceptibility genes are substantially enriched among transcripts sensitive to DDX39B levels (empirical p = 0.00013), suggesting the T allele impairs a broad immune-regulatory splicing program rather than a single target.

The variant maps to chromosome 6p21.3, within the MHC class III region — a genomic area with some of the highest density of immune-related genes in the human genome. Its position in this region means it may co-segregate with other MHC haplotype-specific effects, which complicates precise effect-size estimation but is consistent with strong evolutionary selection pressure on this locus.

Practical Implications

No pharmacological intervention exists that specifically compensates for reduced DDX39B expression. The actionable implications for T allele carriers center on early awareness and monitoring for MS and related autoimmune conditions, and on avoiding environmental triggers that compound autoimmune risk.

MS is a complex disease requiring multiple hits — genetic, environmental (low vitamin D, Epstein-Barr virus infection, smoking, obesity, shift work disrupting circadian rhythms), and stochastic. The T allele is common enough (~16% frequency in Europeans) that most carriers will never develop MS. Nevertheless, the combined double-risk genotype at rs2523506 and rs6897932 confers an OR of approximately 2.75, placing double-risk carriers in a higher-surveillance population.

Modifiable risk factors for MS that are well-supported by evidence include: maintaining adequate serum 25(OH)D (≥50 nmol/L; high-dose D supplementation in MS-risk populations has been studied in clinical trials), avoiding smoking (one of the strongest MS environmental risk factors), maintaining a healthy BMI in adolescence and early adulthood, and — given the EBV connection — awareness that EBV seroconversion in adolescence substantially increases MS risk in genetically susceptible individuals.

Interactions

The defining interaction is between rs2523506 (DDX39B) and rs6897932 (IL7R). These two variants act in the same molecular pathway: DDX39B promotes IL7R exon 6 inclusion; the IL7R rs6897932 C allele weakens the exon 6 splice site. When DDX39B levels are low (T allele at rs2523506) and the splice site is already weak (C allele at rs6897932), exon 6 skipping becomes nearly complete, maximizing sIL7R production and MS risk. This is the epistatic pair described in the 2017 Cell paper with combined OR ≈ 2.75.

This compound interaction is the most important clinical finding — the individual SNP results for rs2523506 are substantially amplified when rs6897932 genotype is known. See the compound action in consolidated_actions.yml for the combined recommendation.

CYP2A6*7 — The Near-Silent Nicotine Enzyme

CYP2A611 CYP2A6
Cytochrome P450 2A6, the main enzyme converting nicotine to cotinine in the liver
is responsible for approximately 80–90% of nicotine's metabolic clearance. The CYP2A6*7 allele carries a missense mutation — isoleucine to threonine at codon 471 — that produces a dramatically unstable enzyme with near-zero nicotine-metabolizing activity. While this allele is rare in most populations (<1%), it reaches 6–15% frequency in East Asian populations, making it clinically relevant for a substantial portion of people with Japanese, Korean, or Chinese ancestry.

The Mechanism

The c.1412T>C variant22 c.1412T>C variant
nucleotide change on the coding strand; A>G on the plus strand at chr19:40843869
substitutes threonine for isoleucine at position 471 in the CYP2A6 protein. This amino acid change occurs in a region critical for heme coordination and structural stability, causing remarkable reduction of holoenzyme stability at body temperature33 remarkable reduction of holoenzyme stability at body temperature
Ariyoshi N et al. Drug Metab Dispos, 2001
. The result is an enzyme that almost entirely lacks nicotine C-oxidase activity44 almost entirely lacks nicotine C-oxidase activity
the primary reaction that converts nicotine to cotinine, measured as cotinine/nicotine ratio
while paradoxically retaining partial coumarin 7-hydroxylase activity. This substrate-selective loss of function reflects the different binding orientations of the two substrates within the enzyme active site.

The Evidence

A Japanese/Korean phenotyping study55 Japanese/Korean phenotyping study
Yoshida R et al. Br J Clin Pharmacol, 2002
in 301 subjects measured plasma nicotine and cotinine concentrations after a standardized nicotine gum dose. All five individuals with impaired metabolism carried either CYP2A6*7 or CYP2A6*4, confirming that CYP2A6*7 alone causes meaningful in vivo metabolic impairment. The allele frequency was 6.5% in Japanese and 3.6% in Korean subjects.

Population-level consequences are measurable. Schoedel et al. Pharmacogenetics, 200466 Schoedel et al. Pharmacogenetics, 2004 compared 375 smokers to 224 non-smokers, finding that slow metabolizers consumed approximately 21.3 cigarettes/day versus 28.2 for normal metabolizers (p = 0.003) and were significantly underrepresented among nicotine-dependent smokers (OR 0.52; 95% CI 0.29–0.95).

For oncology, letrozole (an aromatase inhibitor used in breast cancer) is substantially cleared by CYP2A6. A study of 284 breast cancer patients found >10-fold interpatient variability in letrozole plasma concentrations strongly associated with CYP2A6 genotype77 >10-fold interpatient variability in letrozole plasma concentrations strongly associated with CYP2A6 genotype
Desta Z et al. Clin Pharmacol Ther, 2011
(p < 0.0001). A larger 617-patient pharmacokinetic study confirmed CYP2A6 slow metabolizers have 46% lower apparent letrozole clearance88 CYP2A6 slow metabolizers have 46% lower apparent letrozole clearance
Puszkiel A et al. Eur J Pharm Sci, 2024
compared to normal metabolizers, predicting substantially elevated drug exposure.

For tegafur (a prodrug converted to 5-fluorouracil by CYP2A6), in vitro characterization of recombinant CYP2A6*7 enzyme demonstrated markedly lower Vmax values99 demonstrated markedly lower Vmax values
Yamamiya I et al. Drug Metab Dispos, 2014
for both enantiomers versus wild-type, suggesting reduced prodrug activation in carriers.

Practical Actions

Carriers of CYP2A6*7 metabolize nicotine slowly. In practical terms, this means nicotine lingers longer in the body after each cigarette, reducing the urge to smoke again quickly. Slow metabolizers naturally smoke fewer cigarettes and become less deeply dependent. However, the cessation picture is nuanced: nicotine replacement therapy may be less effective for slow metabolizers1010 nicotine replacement therapy may be less effective for slow metabolizers
Chen LS et al. Biol Psychiatry, 2014
because NRT works partly by maintaining steady nicotine levels to suppress withdrawal — but slow metabolizers already clear nicotine slowly, so the patch's advantage is attenuated. Non-nicotine cessation approaches (varenicline, bupropion) are unaffected by CYP2A6 status.

For letrozole treatment, elevated drug exposure from slow CYP2A6 metabolism may increase the risk of musculoskeletal side effects (arthralgia, myalgia) that are dose-dependent. Prescribers managing breast cancer patients with CYP2A6*7 should be aware that standard doses will produce higher circulating letrozole concentrations.

Interactions

CYP2A6 activity is additive across alleles. Carriers of one CYP2A6*7 allele plus another reduced-function allele (such as rs1801272 CYP2A6*2, rs28399463 CYP2A6*5, or the deletion allele CYP2A6*4) may have near-complete loss of CYP2A6 activity — a compound heterozygote phenotype equivalent to a null metabolizer. The nicotine metabolite ratio (3′-hydroxycotinine to cotinine) is a reliable functional biomarker of overall CYP2A6 activity and can guide cessation treatment selection independently of genotype.

DYRK1B R102C — The Rare Kinase Mutation Behind Hereditary Metabolic Syndrome

DYRK1B11 Dual-specificity tyrosine-phosphorylation-regulated kinase 1B — a serine/threonine kinase that regulates cell cycle exit, quiescence, and differentiation sits at a cellular crossroads: it puts dividing cells to sleep and shepherds them toward differentiated fates. In fat tissue, this means steering precursor cells toward full adipocyte differentiation. The R102C variant tips that balance hard in one direction, flooding carriers with fat-promoting signals they never should have received.

This variant is not a common risk factor — it is a rare, high-penetrance disease mutation. First described in a landmark 2014 New England Journal of Medicine study22 2014 New England Journal of Medicine study
Keramati et al. A form of the metabolic syndrome associated with mutations in DYRK1B. NEJM, 2014
, the R102C substitution was found to co-segregate perfectly with a devastating metabolic syndrome phenotype — dubbed AOMS3 (Abdominal Obesity-Metabolic Syndrome 3; OMIM #615812) — across three unrelated multigenerational families. Every carrier developed the syndrome; every non-carrier did not.

The Mechanism

The R102C mutation swaps arginine for cysteine at position 102, a residue within the kinase-like domain33 kinase-like domain
The catalytic core of the DYRK1B protein that phosphorylates target proteins
that is conserved across vertebrates. At the molecular level, the mutation impairs the HSP90/CDC37 chaperone-assisted folding44 impairs the HSP90/CDC37 chaperone-assisted folding
The R102C protein shows enhanced CDC37 co-chaperone binding, indicating conformational instability of the kinase domain
of the kinase domain, causing more than half of the mutant protein to accumulate in detergent-insoluble aggregates55 detergent-insoluble aggregates
Approximately 50% of mutant DYRK1B vs ~10% of wild-type accumulates in insoluble cytoplasmic fractions
within cells. Paradoxically — and critically — the soluble mutant fraction retains and potentiates gain-of-function effects on downstream targets.

The cellular consequence unfolds through two main channels. First, DYRK1B normally inhibits the SHH (Sonic Hedgehog) and Wnt signaling pathways66 SHH (Sonic Hedgehog) and Wnt signaling pathways
Two developmental pathways that suppress fat cell differentiation when active
. R102C potentiates this inhibition, removing the brake on adipogenesis and causing precursor cells to differentiate into fat cells faster and more completely. Cells expressing R102C show elevated PPARγ, C/EBPα, and C/EBPβ expression77 PPARγ, C/EBPα, and C/EBPβ expression
Master transcription factors that drive adipocyte differentiation; elevated levels in patient-derived adipose stem cells carrying R102C
— the master drivers of fat cell identity. Second, DYRK1B promotes expression of glucose-6-phosphatase (G6Pase)88 glucose-6-phosphatase (G6Pase)
The enzyme that releases glucose from the liver into the bloodstream, a key target for diabetes drugs
, the enzyme central to hepatic glucose output; the R102C variant further amplifies G6Pase expression, driving excess glucose release into the bloodstream and compounding insulin resistance.

The Evidence

The original Keramati et al. NEJM 201499 Keramati et al. NEJM 2014
Keramati AR et al. A form of the metabolic syndrome associated with mutations in DYRK1B. N Engl J Med, 2014
study identified the R102C mutation in three large multigenerational Iranian families. The mutation exhibited complete cosegregation1010 complete cosegregation
Every family member carrying the mutation had the full metabolic syndrome phenotype; none without it did
with the disease phenotype and was absent in 2,000 ethnically-matched Iranian controls, 3,600 Caucasian US controls, and multiple large exome databases. The clinical phenotype was striking: heterozygous carriers developed central obesity, hypertriglyceridemia, low HDL cholesterol, type 2 diabetes, hypertension, and early-onset coronary artery disease — the full constellation of metabolic syndrome — typically before age 40.

A 2021 Orphanet study1111 2021 Orphanet study
Mendoza-Caamal et al. Two novel variants in DYRK1B causative of AOMS3. Orphanet J Rare Dis, 2021
expanded the AOMS3 spectrum by identifying two further DYRK1B mutations in unrelated families, with the same autosomal dominant inheritance and full penetrance. Strikingly, affected members showed age-dependent progression1212 age-dependent progression
Central obesity beginning in childhood, morbid obesity by the third decade, overt type 2 diabetes before 40, hypertension emerging in the fifth decade
: central obesity begins in childhood, morbid obesity and hypertriglyceridemia develop before age 40, and hypertension emerges in the fifth decade. The 2022 adipose stem cell study1313 2022 adipose stem cell study
Armanmehr et al. DYRK1B, PPARG, and CEBPB Expression in Adipose-Derived Stem Cells from Patients Carrying DYRK1B R102C. Metab Syndr Relat Disord, 2022
directly compared fat precursor cells from R102C patients versus healthy controls, confirming accelerated adipogenic differentiation at the transcriptional level.

Given the pathogenic classification in ClinVar, autosomal dominant inheritance, and complete penetrance, this variant warrants aggressive metabolic monitoring and intervention in carriers.

Practical Actions

Because DYRK1B R102C causes early-onset central adiposity and metabolic dysregulation through a defined genetic mechanism, carriers benefit from earlier and more intensive monitoring than the general population. Annual fasting lipid panels, glucose, HbA1c, and blood pressure checks from childhood or early adulthood are warranted. Pharmaceutical strategies that reduce hepatic glucose output (metformin) and improve insulin sensitivity (GLP-1 receptor agonists, SGLT-2 inhibitors) address the specific pathways disrupted by this variant. Dietary approaches that limit glycaemic load and saturated fat directly oppose the adipogenic and gluconeogenic overdrive caused by R102C.

Interactions

DYRK1B interacts with PPARγ (PPARG) biology — the R102C mutation drives elevated PPARγ expression, making the transcriptional target of thiazolidinedione drugs highly active. This raises the question of whether pioglitazone or rosiglitazone might paradoxically worsen adipogenesis in carriers. Clinicians managing AOMS3 carriers should be aware that standard insulin sensitizers designed to activate PPARγ may not behave as expected in this context. The interaction between DYRK1B gain-of-function and PPARG rs1801282 (Pro12Ala) genotype has not been studied but could modulate phenotype severity.

rs3816183

HAAO Ile37Val

Moderate Risk Factor

HAAO Ile37Val — When the Sleep Pathway Tips Toward Neurotoxicity

Your body uses the amino acid tryptophan three ways: to make serotonin and melatonin (the sleep hormones), to fuel energy metabolism as NAD+, and — when inflammation is present — to produce quinolinic acid, a potent excitotoxin that keeps the brain in a state of heightened arousal. The enzyme HAAO (3-hydroxyanthranilate 3,4-dioxygenase) sits at the critical branch point that decides how much of your tryptophan ends up as quinolinic acid. The rs3816183 Ile37Val variant alters this enzyme's function, skewing the pathway in ways that raise insomnia risk.

The Mechanism

In the kynurenine pathway, most dietary tryptophan (~95%) is metabolized through kynurenine11 through kynurenine
rather than the 1-2% that becomes serotonin
. At the HAAO step, 3-hydroxyanthranilic acid (3-HAA) is converted by HAAO into a semialdehyde intermediate that spontaneously cyclizes into quinolinic acid (QUIN)22 quinolinic acid (QUIN)
a potent NMDA receptor agonist with excitotoxic properties
. QUIN is then metabolized to NAD+, but this conversion is easily overwhelmed: the neuronal enzyme that degrades QUIN (QPRT) becomes saturated at concentrations around 300 nM, allowing excess QUIN to accumulate and continuously stimulate NMDA receptors.

The Ile37Val substitution (p.Ile37Val) changes a bulky isoleucine to the smaller valine at position 37, near the active site of this iron-dependent dioxygenase. While detailed kinetic studies on this specific variant are limited, population genetics clearly signals functional impact: the T allele (encoding Val37) has been independently identified in two types of studies — as an insomnia risk locus in the landmark 2019 insomnia GWAS, and as a hypospadias risk factor, consistent with altered quinolinic acid synthesis affecting embryonic development pathways.

Critically, inflammation amplifies this genetic susceptibility. Pro-inflammatory cytokines upregulate the upstream enzyme IDO1, dramatically increasing flux through the entire kynurenine pathway. With an already-altered HAAO, more tryptophan is shunted toward quinolinic acid production instead of toward serotonin and its downstream conversion to melatonin.

The Evidence

The primary genetic evidence comes from a massive insomnia GWAS: Jansen et al. (2019) in Nature Genetics33 Jansen et al. (2019) in Nature Genetics
Genome-wide analysis of insomnia in 1,331,010 individuals identifies new risk loci and functional pathways
identified 202 genome-wide significant loci and performed pathway enrichment analysis that highlighted kynurenine pathway genes — HAAO, KYNU, QPRT, and ACMSD — as a convergent insomnia risk cluster. This is mechanistically coherent: all four encode enzymes in the same metabolic branch responsible for determining the kynurenic acid to quinolinic acid ratio.

The metabolic link to sleep quality is directly demonstrated in human data: Cho et al. (2017)44 Cho et al. (2017)
Sleep disturbance and kynurenine metabolism in depression. Journal of Psychosomatic Research
found that sleep disturbance was significantly associated with a reduced kynurenic acid/quinolinic acid (KynA/QA) ratio — meaning a shift toward neurotoxic quinolinic acid dominance.

The mechanistic pathway from QUIN elevation to disrupted sleep is established in experimental models: Pocivavsek et al. (2018)55 Pocivavsek et al. (2018)
Acute kynurenine challenge disrupts sleep-wake architecture in rats
demonstrated that elevated kynurenine (which feeds HAAO to produce QUIN) reduced total REM duration, delayed REM onset, and increased wakefulness, with EEG evidence of impaired theta power during REM — a signature of hippocampal arousal. The mechanism is NMDA receptor hyperactivation and antagonism of α7 nicotinic acetylcholine receptors.

Finally, the NAD+ connection provides a second route to intervention: Weiss (2026)66 Weiss (2026)
Vitamin B3 ameliorates sleep duration and quality
reviewed clinical evidence that nicotinamide riboside (NR) supplementation — which bypasses HAAO to replenish NAD+ — improves sleep efficiency in individuals with insomnia and supports circadian clock gene function (BMAL1, PER2).

Practical Actions

For T-allele carriers, two complementary strategies address the underlying mechanism: (1) reducing inflammation to limit IDO1 induction and the resulting kynurenine flood through HAAO, and (2) supporting NAD+ levels via a pathway that bypasses HAAO entirely. Timing of tryptophan intake also matters — consuming tryptophan-rich foods in the evening rather than splitting intake evenly across the day preferentially supports serotonin/melatonin synthesis while the liver's kynurenine pathway activity is lower.

Interactions

HAAO Ile37Val acts within a broader kynurenine pathway context. The upstream enzyme KMO (kynurenine 3-monooxygenase) and downstream enzyme QPRT both have functional variants that modulate the same neurotoxic/neuroprotective balance. TT carriers may benefit from looking at their KMO and QPRT variants to understand the full picture of their kynurenine pathway function.

The COMT gene (rs4680) interacts with this pathway indirectly: COMT metabolizes catecholamines that modulate hypothalamic arousal, and individuals with both HAAO TT and COMT AA (low dopamine clearance) may experience compounded sleep-onset difficulty via both NMDA hyperactivation and elevated dopaminergic arousal.

rs4148102

ABCG1 ABCG1 intronic PUFA-interaction variant

Emerging Risk Factor

ABCG1 — When Dietary PUFAs Raise, Not Lower, Cholesterol

ABCG1 (ATP-binding cassette transporter G111 ATP-binding cassette transporter G1
a membrane-spanning pump that moves cholesterol and phospholipids from the inner leaflet of cell membranes onto maturing HDL particles, completing the second step of reverse cholesterol transport after ABCA1 initiates lipid loading onto nascent HDL
) encodes one of the body's principal cholesterol efflux transporters. Located on chromosome 21q22.3, ABCG1 is expressed in macrophages, liver, and many other tissues, where it loads mature HDL particles with surplus cellular cholesterol for transport back to the liver — a central step in preventing foam cell formation and atherosclerotic plaque development.

rs4148102 is an intronic variant in ABCG1 that has no known effect on the protein sequence itself, but its location within an intron has raised the possibility of subtle effects on ABCG1 splicing, expression, or its response to dietary lipid signals. What distinguishes this variant is not its baseline effect on cholesterol but the way it conditions the body's response to high polyunsaturated fatty acid (PUFA) intake — producing a counterintuitive pattern where a seemingly heart-healthy dietary pattern raises LDL-cholesterol in AA homozygotes.

The Mechanism

The biological basis for this gene-diet interaction is not fully characterized at the molecular level. However, ABCG1 transcription is under LXR22 LXR
liver X receptor — a nuclear receptor activated by oxysterols and certain fatty acid derivatives; once activated, LXR drives expression of ABCG1, ABCA1, and other cholesterol homeostasis genes
control. Polyunsaturated fatty acids and their metabolites serve as LXR ligands and can modulate its activity. An intronic variant could influence how effectively the ABCG1 gene responds to these dietary lipid signals — for example, by altering a regulatory element that fine-tunes LXR responsiveness.

When ABCG1 function is subtly impaired or mis-regulated, cholesterol efflux from cells to HDL is reduced. Under high-PUFA dietary conditions, which normally promote LDL receptor upregulation and LDL clearance, the AA variant appears to interfere with this adaptive response, resulting in paradoxically elevated LDL and total cholesterol rather than the reduction seen in GG carriers on the same diet.

The Evidence

The primary evidence comes from a study by Abellán et al.33 study by Abellán et al.
Abellán R et al. Dietary polyunsaturated fatty acids may increase plasma LDL-cholesterol and plasma cholesterol concentrations in carriers of an ABCG1 gene single nucleotide polymorphism: study in two Spanish populations. Atherosclerosis, 2011
examining 1,941 participants across two independent Spanish cohorts (the Hortega study, n=1,178, and the Pizarra study, n=763). In the Hortega population, AA homozygotes consuming high-PUFA diets had LDL-cholesterol of 149.8 ± 37.9 mg/dL compared to 111.4 ± 32.1 mg/dL in G allele carriers (p=0.005). Total cholesterol diverged similarly (242.1 vs 198.0 mg/dL, p=0.003). Pooling both cohorts strengthened the signal: gene-diet interaction p=0.006 for total cholesterol, p=0.003 for LDL-cholesterol.

The A allele is relatively uncommon globally (~14.5%), making AA homozygosity rare (~2% of the population). This limits replication opportunities and explains why no dedicated follow-up trials exist yet. The evidence level is therefore emerging — a well-designed two-cohort study, but not yet replicated in independent populations or mechanistically confirmed in cell-based assays.

Practical Actions

For AA homozygotes, the clinical implication is specific: high-PUFA diets — which most guidelines recommend as beneficial — may actually raise LDL-cholesterol in this genotype. Rather than increasing total PUFA from vegetable oils (sunflower, corn, safflower), focus on pre-formed long-chain omega-3s (EPA and DHA from fatty fish or supplements) and monitor fasting LDL when making significant dietary changes. Saturated fat should still be limited — this is not a signal to increase saturated fat intake. The issue is specifically with high omega-6 PUFA loads from vegetable oils, not with EPA/DHA.

Carriers of one A allele (AG) represent roughly 25% of people. The gene-diet interaction was strongest in AA homozygotes; AG heterozygotes showed intermediate patterns in some analyses, but the effect was not consistent across both study populations. Standard dietary PUFA guidance applies.

Interactions

ABCG1's role in cholesterol efflux depends on the upstream ABCA1 step (rs1044317 is another ABCG1 variant in the cholesterol_lipoproteins category). If both ABCG1 steps in the efflux pathway are impaired — the initial ABCA1-mediated lipid loading onto nascent HDL, and the ABCG1-mediated maturation step — the combined cholesterol efflux deficit could be larger than either variant alone. FADS1 variants (rs174547, rs174537) also interact with this locus biologically: impaired FADS1 conversion of plant omega-3s to EPA/DHA shifts the circulating PUFA ratio toward omega-6, which may amplify the LDL-raising signal seen with rs4148102 AA on high-PUFA diets.

rs4810424

R3HDML R3HDML rs4810424

Moderate Risk Factor

HNF4A P2 Promoter: A Beta-Cell Transcription Signal

The R3HDML gene on chromosome 20q13 sits immediately adjacent to a well-characterized regulatory landmark: the P2 promoter of HNF4A (| Hepatocyte Nuclear Factor 4-alpha, a transcription factor critical for pancreatic islet development and insulin secretion). rs4810424 falls within an R3HDML intron, but its biological significance comes from being a reliable tag for the HNF4A P2 promoter risk haplotype — a cluster of SNPs (rs4810424, rs1884613, rs1884614, rs2144908) in nearly complete linkage disequilibrium (D' > 0.97, r² > 0.95) that were originally identified in Finnish and Ashkenazi Jewish populations.

The Mechanism

HNF4A is an orphan nuclear receptor | A transcription factor that binds DNA to regulate gene expression; called "orphan" because its natural ligand was originally unknown that controls genes governing glucose-stimulated insulin secretion, gluconeogenesis, and beta-cell maintenance. Critically, while both the P1 and P2 promoters drive HNF4A expression in liver cells, pancreatic islets rely almost exclusively on the P2 promoter. Common variants tagging the P2 haplotype are therefore expected to alter HNF4A expression disproportionately in the tissue most relevant to type 2 diabetes: the insulin-secreting beta cell.

The P2 haplotype risk alleles are associated with reduced HNF4A transcript levels in islets, impairing the transcriptional cascade that drives expression of insulin, glucose transporter-2 (GLUT2), L-pyruvate kinase, and other metabolic coupling genes. The functional consequence is a blunted insulin secretory response to glucose — a mechanism consistent with the OGTT-based associations observed in Japanese and Thai cohorts | Tokunaga et al., Endocrine J, 2008, PMID 18654034; Jongjaroenprasert et al., Acta Diabetol, 2007, PMID 17805472.

The Evidence

The original discovery came from the FUSION study in Finland, where Silander et al.11 Silander et al.
Silander K et al. Genetic variation near the hepatocyte nuclear factor-4 alpha gene predicts susceptibility to type 2 diabetes. Diabetes, 2004
identified the P2 haplotype association with T2D (OR 1.33, 95% CI 1.06-1.65, p = 0.011) in 495 Finnish affected sibling pairs, with independent replication in Ashkenazi Jewish samples.

A large UK replication study Weedon et al.22 Weedon et al.
Weedon MN et al. Common variants of the hepatocyte nuclear factor-4alpha P2 promoter are associated with type 2 diabetes in the UK population. Diabetes, 2004
confirmed the association in 5,256 Caucasians (OR 1.15, 95% CI 1.02-1.33, p = 0.02). The smaller effect size in the UK compared with Ashkenazi populations (OR ~1.7) suggests the causal variant within the haplotype block differs in frequency or is modified by population background.

In the prospective STOP-NIDDM trial33 STOP-NIDDM trial
Andrulionyte L et al. SNPs of the HNF4alpha gene are associated with the conversion to type 2 diabetes mellitus: the STOP-NIDDM trial. J Mol Med, 2006
, female carriers of the C allele at rs4810424 had a 1.7-fold elevated risk of converting from impaired glucose tolerance to diabetes (OR 1.7, 95% CI 1.09-2.66, p = 0.020), while no association was detected in men. A combined Swedish-Finnish cohort Holmkvist et al.44 Holmkvist et al.
Holmkvist J et al. Common variants in maturity-onset diabetes of the young genes and future risk of type 2 diabetes. Diabetes, 2008
found rs4810424 predicted future T2D (HR 1.3, 95% CI 1.0-1.6, p = 0.03) across 17,831 participants.

The association is notably population-specific: the risk effect is strongest in Ashkenazi Jewish individuals, moderate in Scandinavians and Finns, and not consistently replicated in UK Caucasians. The C allele is also far more common in East Asian populations (~44% vs ~17% in Europeans), though its T2D association in Asian cohorts is primarily mediated by rs1884614 and rs2144908 rather than rs4810424 itself.

Practical Actions

Carriers of the risk haplotype — particularly CC homozygotes — have a modestly elevated risk of beta-cell secretory insufficiency that is amenable to dietary and monitoring strategies. A diet that reduces the post-meal glucose load (lower glycemic index, reduced refined carbohydrate) directly offsets the beta-cell demand amplified by impaired HNF4A expression in islets. Prospective data from the STOP-NIDDM trial also showed that acarbose — an alpha-glucosidase inhibitor — partially offset the progression to T2D in high-risk individuals with impaired glucose tolerance, consistent with a glucose-spike mechanism.

Interactions

rs4810424 is in high LD with rs2144908 (r² ≈ 0.97), rs1884614, and rs1884613 within a ~15 kb P2 haplotype block. Carrying risk alleles at multiple haplotype members compounds the beta-cell expression deficit but does not add independent risk beyond the haplotype signal. Separately, published data show an interaction between HNF4A rs2144908 and KCNJ11 E23K (rs5219) — carriers of both risk variants show markedly greater impairment of insulin secretion than either alone, consistent with additive disruption of the ATP- sensitive potassium channel pathway that HNF4A transcriptionally regulates.