NUCB2 rs1330 — Intronic Nesfatin-1 Variant with Sex-Specific Metabolic and Oncological Effects
Nucleobindin-2 (NUCB2) encodes the precursor protein that is proteolytically cleaved to release
nesfatin-111 nesfatin-1
An 82-amino acid neuropeptide that suppresses appetite and modulates energy balance via melanocortin MC3/MC4 receptors and CRF2 — operates independently of the leptin pathway,
a neuropeptide with broad roles in appetite suppression, glucose regulation, sleep-wake cycling,
and — more recently discovered — cancer biology. The rs1330 variant sits within an
intron of NUCB2 on chromosome 11 (GRCh38 position 17,294,482), meaning it does not
alter the nesfatin-1 amino acid sequence directly. Instead, its effects are regulatory:
the T allele is thought to influence NUCB2 splicing efficiency, transcript stability,
or expression levels, though the precise molecular mechanism has not yet been resolved.
The T allele is the minor allele in most populations (~30–39% frequency globally, with notably lower frequency in African populations at ~17%). The reference C/C genotype represents wild-type nesfatin-1 expression; T allele carriers show associations across multiple disease contexts with a consistent direction — higher risk — though with sex-specific patterns that suggest hormonal context shapes how this intronic variant exerts its effects.
The Mechanism
As an intron variant, rs1330 does not change the nesfatin-1 peptide sequence. Its regulatory impact is inferred from its population-level associations. The T allele may alter the efficiency of RNA splicing at one of NUCB2's five transcript variants, reduce mRNA stability, or affect transcription factor binding in an intronic regulatory element — any of which would reduce the amount of functional NUCB2 protein available for cleavage into nesfatin-1. Reduced circulating nesfatin-1 is the common downstream consequence observed in obesity and insulin-resistant states across multiple NUCB2 studies, and the rs1330 T allele likely compounds this deficiency via a cis-regulatory mechanism.
The sex-specific pattern observed in the
Zegers et al. obesity study22 Zegers et al. obesity study
Zegers D et al. Association between polymorphisms of the Nesfatin gene, NUCB2, and obesity in men. Mol Genet Metab, 2011
— where rs1330 associated with obesity only in males — and the reciprocal
Li et al. T2D finding33 Li et al. T2D finding
Li XS et al. NUCB2 polymorphisms are associated with an increased risk for type 2 diabetes in the Chinese population. Endocr Connect, 2020
in females (OR 1.31–1.42) may reflect estrogen and androgen modulation of
NUCB2 expression. Estrogen upregulates nesfatin-1 production in certain
hypothalamic circuits, potentially masking the effect of a reduced-expression
variant in pre-menopausal women while leaving males more exposed to its impact
on appetite and energy balance.
The Evidence
Obesity (males). The
Zegers et al. 201144 Zegers et al. 2011
Zegers D et al. Association between polymorphisms of the Nesfatin gene, NUCB2, and obesity in men. Mol Genet Metab, 2011
case-control study genotyped 1,049 obese and 315 normal-weight Caucasian subjects
across eight NUCB2 tagSNPs. Three variants — rs1330, rs214101, and rs757081 —
showed association with obesity protection, but only in the male sub-analysis.
Linear regression further linked rs1330 to BMI, body weight, and fat-free mass in
men. This was the first evidence that NUCB2 intronic variants influence energy
homeostasis in humans.
Type 2 diabetes (females). A Chinese Han study by
Li et al. (2020)55 Li et al. (2020)
Li XS et al. NUCB2 polymorphisms are associated with an increased risk for type 2 diabetes in the Chinese population. Endocr Connect, 2020
genotyped 578 T2DM patients against 1,609 healthy controls and identified rs1330
as significantly associated with T2DM risk in women (OR 1.31–1.42, P<0.05).
The variant was also correlated with BMI in the female subpopulation, suggesting
its influence spans both adiposity and glucose metabolism in women. No association
was detected in men in this cohort — the reverse sex-specificity pattern from
the Zegers obesity data.
Colorectal cancer. A Mexican cross-sectional study by
Macías-Gómez et al. (2025)66 Macías-Gómez et al. (2025)
Macías-Gómez NM et al. Variants in the neuropeptide gene NUCB2 as a possible biomarker for colorectal cancer. 2025
enrolled 397 CRC patients and 383 healthy controls. The TT genotype at rs1330
was significantly associated with colorectal cancer (OR 2.66, P<0.001), while
CT heterozygotes showed an apparent protective pattern (OR 0.61), an unusual
non-additive effect possibly reflecting over-dominant heterozygote advantage.
Hardy-Weinberg deviation in the CRC group, however, warrants cautious interpretation.
Oral cancer. A Taiwanese cohort study by
Yu et al. (2026)77 Yu et al. (2026)
Yu CC et al. Association of NUCB2 genetic variants with the clinicopathological features of oral cancer. 2026
examined four NUCB2 polymorphisms in men with oral cancer. Compared to C/C wild-type,
carriage of at least one T allele at rs1330 was associated with elevated risk of
disease progression to stage III/IV, particularly in patients aged ≥60 years.
Sleep biology. rs1330 has not been studied directly in sleep GWAS, but the
parent gene's sleep-regulatory role is established. In rat models,
Vas et al. (2013)88 Vas et al. (2013)
Vas S et al. Nesfatin-1/NUCB2 as a potential new element of sleep regulation in rats. PLoS One, 2013
demonstrated that central nesfatin-1 reduces REM sleep and increases wakefulness,
and that hypothalamic NUCB2 expression declines during REM sleep deprivation and
rebounds during recovery. Any NUCB2 variant that reduces nesfatin-1 output may
therefore modulate sleep architecture, though this remains inferred biology
at the rs1330 level.
Practical Implications
For CT and TT carriers, the most actionable implications are in metabolic monitoring and cancer awareness. Given the sex-specific patterns — obesity/BMI associations in males, T2D risk in females — sex-targeted monitoring is appropriate. Dietary interventions that support nesfatin-1 release, particularly high-protein morning meals, are theoretically beneficial but are based on functional biology rather than rs1330-specific trial data.
Interactions
rs1330 was identified in the same study (Zegers et al. 2011) as rs757081 and rs214101 — all three intronic and coding variants in NUCB2 showed male-specific obesity associations together, suggesting they may function as a regulatory haplotype influencing total NUCB2/nesfatin-1 output. Combined carriage of rs1330 T allele and rs757081 C allele (the obesity-risk form at the missense position) may additively reduce effective nesfatin-1 activity through both regulatory and protein-level mechanisms.
IL23R rs1343151 — A Haplotype Tag With Dual Immune Consequences
The interleukin-23 receptor (IL23R) gene encodes a key checkpoint in adaptive
immunity. When IL-23 binds IL23R on T helper 17 (Th17) cells, it triggers
STAT3 phosphorylation11 STAT3 phosphorylation
STAT3 is a transcription factor that, once activated,
drives production of the inflammatory cytokines IL-17A and IL-22,
sustaining chronic inflammation in the gut, spine, and skin. The rs1343151 variant
sits in an intron of IL23R at chromosome 1 position 67,253,446 (GRCh38) and carries
no direct amino acid change. Its clinical significance comes from two distinct
properties: it tags the protective IL23R haplotype that dampens Th17 responses in
ankylosing spondylitis (AS) and inflammatory bowel disease (IBD), and it carries an
independent association with rheumatoid arthritis (RA) susceptibility that operates
through a separate biological mechanism.
The Mechanism
As an intronic variant, rs1343151 does not directly alter the IL-23 receptor protein.
Instead, it acts as a tag SNP22 tag SNP
a variant that is inherited together with nearby
functional variants so frequently that knowing its allele predicts those of its
neighbours; here rs1343151 co-inherits with rs11465804 and the functional missense
variant rs11209026 (R381Q) for the AS
and IBD associations. The biological engine is the R381Q substitution at rs11209026:
replacing arginine with glutamine at position 381 in the cytoplasmic tail of IL23R
partially uncouples the receptor from JAK2/STAT3 signalling, reducing IL-17A output
in Th17 effector cells from approximately 36 pg/ml to 5.5 pg/ml under IL-23
stimulation — a 6.5-fold reduction.
The RA association, by contrast, appears to be at least partially independent of the
R381Q haplotype. Hollis-Moffatt et al. 200933 Hollis-Moffatt et al. 2009
Ann Rheum Dis; 3,000+ Caucasian RA
cases and 3,800+ controls found no
association of rs11209026 with RA (OR 1.01) while rs1343151 showed an OR of 1.14
in the same dataset — an unusual dissociation that suggests rs1343151 is tagging a
different regulatory element or haplotype block for RA than for IBD and AS.
The Evidence
For Crohn's disease and ulcerative colitis, the protective effect of the A allele
is well-replicated in European populations. A New Zealand cohort study
(Ferguson et al. 201144 (Ferguson et al. 2011
Gastroenterol Res Pract; 339 CD cases, 407 controls)
found an allelic OR of 0.68 (P=0.001), with AA homozygotes showing an OR of 0.29
(95% CI 0.16–0.53) relative to GG individuals — a striking 71% reduction in CD odds.
A subsequent meta-analysis of 15 studies55 meta-analysis of 15 studies
Ding et al. Sci Rep 2015
confirmed this with a pooled OR of 0.725 (95% CI 0.690–0.763) across Caucasian
populations. Notably, no protective signal was detected in Asian populations, consistent
with the low A allele frequency in East Asian cohorts (~8%).
For ankylosing spondylitis, a meta-analysis of 25 case-control studies66 meta-analysis of 25 case-control studies
Zhong
et al. Expert Rev Clin Immunol 2018; 8,431 AS cases, 8,972 controls
confirmed that the A allele frequency was significantly lower in AS patients than
controls (P<0.001), ranking rs1343151 among four IL23R polymorphisms with consistent
protective signals against AS. The study-level OR for rs1343151 in an earlier meta-analysis
of ankylosing spondylitis cohorts was approximately 0.68 (95% CI 0.55–0.83).
For rheumatoid arthritis, the picture reverses: the A allele is associated with a modest increase in RA risk (OR ~1.11–1.14 in Caucasians), which has been independently replicated across multiple cohorts. This disease-direction reversal at the same allele — protective for gut and spine inflammation, risky for synovial inflammation — reflects the differential roles of IL-23/Th17 signalling in different tissue compartments and disease pathologies.
Practical Implications
For individuals concerned about IBD or spondyloarthritis, the A allele at rs1343151 is a protective marker. AA homozygotes have substantially lower lifetime odds of Crohn's disease and ankylosing spondylitis than GG homozygotes, in parallel with other IL23R protective variants. For RA, the relationship inverts — carriers of the A allele have slightly elevated susceptibility. However, the absolute risk difference conferred by this single variant for RA is small (OR ~1.11) compared with the protective effect on IBD and AS (OR ~0.68–0.72). IL-23 pathway targeting (with biologics such as risankizumab or guselkumab) is highly effective for AS and CD; carriers of the GG genotype at rs1343151 — lacking the protective haplotype — may be among those most likely to benefit from such therapies if they develop these conditions.
Interactions
rs1343151 is in moderate-to-strong linkage disequilibrium with several other IL23R variants on chromosome 1p31.3, including rs11465804, rs10489629, and the functional missense variant rs11209026 (R381Q). For the IBD and AS associations, these variants are largely interchangeable as they tag the same protective haplotype block; a person carrying the protective A allele at rs1343151 will almost always also carry the protective alleles at rs11465804 and rs10489629.
The RA signal at rs1343151, being independent of rs11209026, represents a distinct haplotype association and likely reflects different regulatory architecture in synovial versus intestinal/entheseal immune environments. The IL23R susceptibility variants rs2201841 and rs1004819 tag a separate, risk-conferring haplotype block and are tracked individually in this database.
FBN1 Asp1113Gly — A Missing Anchor in the Aortic Wall
The aorta tolerates tens of millions of pressure pulses over a lifetime because its
wall is reinforced by an elastic scaffold made of
microfibrils11 microfibrils
rope-like extracellular matrix structures assembled from fibrillin-1
molecules linked end-to-end, providing elasticity and tensile strength throughout
connective tissue.
Fibrillin-1 is a 2,871-amino-acid glycoprotein encoded by FBN1 on chromosome 15.
It contains 47 epidermal growth factor-like (EGF-like) domains, 43 of which are
calcium-binding (cbEGF). Each cbEGF domain requires calcium to fold into a rigid
rod — calcium ions are not just incidental cofactors but the structural glue that
holds the entire module together. Aspartate residues at conserved positions within
these domains directly coordinate the calcium ion, and substituting any one of
them disrupts the fold.
The p.Asp1113Gly variant (coding-strand c.3338A>G on the minus-strand FBN1 transcript; plus-strand T→C at chr15:48,487,437) removes the calcium-coordinating aspartate at position 1113 and replaces it with glycine — the smallest amino acid, with no side chain at all. ClinVar classifies this variant as Pathogenic or Likely Pathogenic across three independent submissions (RCV000663632, RCV000802396, RCV002284420), with conditions listed as Marfan syndrome and familial thoracic aortic aneurysm and aortic dissection (FTAAD). The global allele frequency of the C allele is approximately 0.02% (20 of 88,940 ALFA samples), consistent with a rare, high-penetrance disease variant under negative selection.
The Mechanism
Within a cbEGF domain, a conserved sequence motif coordinates a single calcium
ion through the side chains of several residues — aspartate residues contribute
critically to this coordination shell.
Disrupting calcium coordination destabilizes the cbEGF fold22 Disrupting calcium coordination destabilizes the cbEGF fold
Dietz et al.
1993 showed that all early Marfan syndrome missense mutations at that time
affected residues with calcium-binding significance in EGF-like domains,
preventing fibrillin-1 from achieving the extended rigid conformation needed for
microfibril polymerization. The misfolded monomer acts through a
dominant-negative mechanism33 dominant-negative mechanism
the defective protein physically interferes with
normal fibrillin-1 assembly, impairing the entire microfibril network even when
the second FBN1 allele is intact —
which is why a single copy of this variant is sufficient to cause disease.
Weakened microfibrils in the aortic media reduce mechanical resilience of the aortic wall and, importantly, fail to sequester transforming growth factor-beta (TGF-β) in the extracellular matrix. Excess free TGF-β drives smooth muscle cell dysfunction, progressive aortic wall stiffening, and aneurysmal dilation at the sinuses of Valsalva — the anatomical pinch-point where nearly all Marfan-related dissections originate.
The Evidence
The broader clinical evidence base for FBN1 cbEGF domain missense variants predicting aortic disease is well-established:
Tan et al. (Human Molecular Genetics, 2017)44 Tan et al. (Human Molecular Genetics, 2017) sequenced FBN1 in 687 patients with sporadic non-syndromic aortic dissection and found pathogenic FBN1 variants in 3.9% (27 patients), the majority being missense mutations — demonstrating that FBN1 pathogenic variants cause FTAAD even in the absence of full syndromic Marfan features.
The most important prognostic data for FBN1 pathogenic variant carriers come from Milleron et al. (JACC, 2020)55 Milleron et al. (JACC, 2020): 954 patients with FBN1 pathogenic variants followed over 8,594 patient-years showed that type A aortic dissection occurred at only 0.4 events per 1,000 patient-years when maximum aortic diameter remained below 50 mm on guideline-directed treatment. This is a critically important number: risk is not inevitable but is profoundly modifiable by structured surveillance and timely intervention.
On treatment, Brooke et al. (NEJM, 2008)66 Brooke et al. (NEJM, 2008) showed that losartan (an angiotensin receptor blocker that blocks TGF-β pathway overactivation) reduced aortic root dilation rate from 3.54 mm/yr to 0.46 mm/yr in 18 pediatric Marfan patients — approximately an 87% reduction. This finding underpins current guidelines recommending ARBs when aortic root dilation progresses faster than 5 mm/year or root diameter exceeds 40 mm in children.
Practical Actions
The cornerstone of management is annual transthoracic echocardiography to track aortic root diameter at the sinuses of Valsalva. Medical therapy with a beta-blocker (to reduce heart rate and aortic wall stress) or ARB (losartan, to block excess TGF-β signaling) should be initiated at diagnosis or upon detecting progressive dilation. Prophylactic aortic root replacement is recommended when the aortic root approaches 5.0 cm in adults (4.5–5.0 cm in rapid progressors or patients with a family history of dissection).
Fluoroquinolone antibiotics (ciprofloxacin, levofloxacin, moxifloxacin) inhibit lysyl oxidase — the enzyme that cross-links collagen and elastin — and carry an FDA black-box warning for aortic aneurysm rupture and dissection risk, particularly in patients with connective tissue disorders. All prescribers must be informed of this contraindication.
Strenuous isometric exercise (heavy weightlifting, contact sports, competitive athletics) should be avoided — these modalities generate abrupt spikes in arterial pressure that place disproportionate mechanical stress on an already vulnerable aortic wall.
Cascade genetic testing of first-degree relatives is clinically essential. A single positive finding in one family member should trigger testing of all parents, siblings, and adult children, since autosomal dominant transmission gives each relative a 50% chance of carrying the variant — and several may harbor silent aortic dilation detectable only by imaging.
Interactions
FBN1 pathogenic variants, including D1113G, interact clinically with variants in TGF-β pathway genes (TGFBR1, TGFBR2, SMAD2, SMAD3) — the genetic basis of Loeys-Dietz syndrome. Carriers who also harbor pathogenic variants in ACTA2 (smooth muscle alpha-actin) or MYH11 (smooth muscle myosin heavy chain), both associated with independent familial TAAD syndromes, may warrant more aggressive surveillance thresholds, though formal compound genotype data for D1113G specifically are not published.
Other FBN1 cbEGF domain missense variants in the database — including rs193922239 (p.Gly2627Arg) — share the same dominant-negative mechanism and clinical management protocol. The diagnosis of Marfan syndrome versus FTAAD-without-full-Marfan-features depends on the presence of ocular and skeletal manifestations alongside the cardiovascular findings; either presentation requires the same aortic surveillance program.
FN1 and Alzheimer's Disease — A Blood-Brain Barrier Protective Variant
Fibronectin 1 (FN1) encodes a large glycoprotein that forms the scaffold of the extracellular matrix
around blood vessels. In the brain, fibronectin plays a structural role in the neurovascular unit11 neurovascular unit
the interface between blood vessels and neurons, comprising endothelial cells, pericytes, astrocyte
endfeet, and neurons, helping maintain blood-brain barrier
integrity and coordinating the local response to injury. The rs140926439 variant introduces an amino
acid substitution in fibronectin (p.Gly357Glu) that appears to reduce its pathological accumulation
at the blood-brain barrier — and in APOE ε4 carriers, this reduction translates into substantially
lower Alzheimer's disease risk.
The Mechanism
In APOE ε4 carriers, fibronectin abnormally accumulates in the walls of brain blood vessels22 In APOE ε4 carriers, fibronectin abnormally accumulates in the walls of brain blood vessels
This vascular fibronectin deposition correlates with impaired clearance of toxic amyloid-beta
peptides and heightened neuroinflammation. The proposed
pathway: APOE ε4 drives excess fibronectin deposition → thickened vascular walls impair gliovascular
remodeling and microglial amyloid clearance → amyloid-beta accumulates → neurodegeneration
accelerates.
The rs140926439 T allele encodes a Glu (glutamic acid) at position 357 instead of the common Gly
(glycine). This Gly→Glu substitution introduces a charged residue into a region of fibronectin involved
in matrix assembly, likely disrupting the protein's ability to multimerize or bind to vascular
components in the same pathological configuration. Functional studies including zebrafish models
support the idea that reduced FN1 function in the neurovascular unit enhances microglial responses
and gliovascular remodeling33 Functional studies including zebrafish models
support the idea that reduced FN1 function in the neurovascular unit enhances microglial responses
and gliovascular remodeling
Loss-of-function FN1 variants appear to be protective rather than
harmful in the context of APOE ε4-driven amyloid accumulation.
The Evidence
The primary evidence comes from a multi-cohort whole-genome sequencing study combining NIA-AD FBS,
WHICAP, and EFIGA cohorts, followed by validation in 7,185 APOE ε4 homozygous carriers44 The primary evidence comes from a multi-cohort whole-genome sequencing study combining NIA-AD FBS,
WHICAP, and EFIGA cohorts, followed by validation in 7,185 APOE ε4 homozygous carriers
Bhattarai et al. (Acta Neuropathologica, 2024) identified rs140926439 as a rare coding protective
variant with OR = 0.29 (95% CI 0.11–0.78, p = 0.014), and found it delayed Alzheimer's disease
onset by 3.37 years (95% CI 0.42–6.32, p = 0.025) in ε4 carriers.
The pathway analysis showed significant enrichment in extracellular matrix processes among protective
variants, consistent with FN1's role in vascular matrix biology.
An independent replication using UK Biobank data confirmed the protective signal55 An independent replication using UK Biobank data confirmed the protective signal
Lehrer & Rheinstein (Alzheimer Disease and Associated Disorders, 2025) found AD prevalence of 0.43%
among APOE ε4 carriers without the variant versus 0.10% among carriers with it — essentially
eliminating the ε4 risk increment. The authors propose
that fibronectin-driven pathological brain wound-healing may underlie a subset of Alzheimer's cases,
and that variants impairing this cascade are protective.
Evidence is classified as emerging: the finding comes from two cohort studies (one discovery + functional validation, one replication), without clinical trial data. The T allele is rare (global frequency ~0.3%), limiting statistical power and making replication in larger datasets essential. The ClinVar classification is "Likely Benign" for the variant in isolation, consistent with its rarity and the context-dependent (APOE ε4-specific) protective effect.
Practical Actions
The T allele confers meaningful protection specifically in APOE ε4 carriers. For CT carriers who are also APOE ε4 positive, this variant provides a counterbalancing protective signal that significantly modifies their overall Alzheimer's risk trajectory. Since the variant is rare and evidence emerging, the actionable value is primarily informational — it informs how aggressively to pursue standard Alzheimer's prevention strategies rather than introducing new interventions specific to this variant alone.
Monitoring for early cognitive changes remains prudent for ε4 carriers regardless of FN1 status, as the T allele reduces but does not eliminate risk.
Interactions
The protective effect of rs140926439 is specifically documented in APOE ε4 carriers (rs429358 C allele + rs7412 C allele combination). Outside of the ε4 context, the variant has not shown significant effect on Alzheimer's risk — the interaction between FN1 and APOE at the neurovascular unit appears central to the mechanism.
Other Alzheimer's risk variants modulating blood-brain barrier function include rs3851179 (PICALM), which affects amyloid-beta transcytosis independently of APOE. An individual carrying both protective rs140926439 and the PICALM protective A allele would theoretically benefit from two complementary blood-brain barrier protection mechanisms, though no direct evidence for this combination exists.
GDF5 and Joint Health — A Genetic Influence on Cartilage Longevity
The GDF5 gene encodes growth differentiation factor 511 growth differentiation factor 5
a member of the bone morphogenetic
protein (BMP) family essential for skeletal development,
particularly in forming and maintaining cartilage in synovial joints.
rs143383
is a C to T transition SNP located in the 5'untranslated region (5'UTR) of the GDF5 gene
.
This regulatory region controls how much GDF5 protein your cells produce, and
the T
allele of the SNP is associated with increased risk of osteoarthritis (OA) in Europeans and in
Asians
.
GDF5 is on the minus strand of chromosome 20, and this is an intron variant located in a critical regulatory region. While technically classified as an intron variant in some databases, it functions as a regulatory element in the gene's 5' UTR, affecting transcription.
The Mechanism
The A allele produces less GDF5 transcript relative to the G allele, a phenomenon known as differential allelic expression (DAE) . Studies show an average of 27% lower expression of the disease-associated A allele than the G allele in synovial joint tissues from OA patients . The mechanism involves transcription factor binding:
Sp1, Sp3, and DEAF-1 are repressors of GDF5 expression, with DEAF-1 modulating the differential allelic effect — the rs143383 A allele being repressed to a significantly greater extent than the G allele .
This reduced expression matters because GDF5 is essential for cartilage homeostasis. It promotes
chondrocyte differentiation, stimulates production of cartilage matrix proteins like aggrecan and
type II collagen22 cartilage matrix proteins like aggrecan and
type II collagen
key structural components of healthy cartilage,
and supports joint repair processes. Less GDF5 means less cartilage maintenance capacity over time.
The Evidence
GDF5 is the most compelling candidate association signal so far reported for OA, with the rs143383 single nucleotide polymorphism (SNP) showing association in both Europeans and Asians and at a significance level of P < 5.0 × 10⁻⁸ . The evidence spans multiple joint sites:
Knee osteoarthritis33 Knee osteoarthritis
Valdes et al. (2011). The GDF5 rs143383 polymorphism is associated with
osteoarthritis of the knee with genome-wide statistical significance. Ann Rheum Dis.:
A significant random-effects summary OR for knee OA was demonstrated for rs143383 (1.15 [95% confidence interval 1.09-1.22]) (P=9.4×10⁻⁷), with no significant between-study heterogeneity . This is the strongest and most consistent association.
Lumbar disc degeneration44 Lumbar disc degeneration
Williams et al. (2011). GDF5 single-nucleotide polymorphism rs143383 is
associated with lumbar disc degeneration in Northern European women. Arthritis Rheum.:
An association between LDD and the SNP rs143383 was identified in women, with the same risk allele as in knee and hip OA (odds ratio 1.72 [95% confidence interval 1.15–2.57], P = 0.008) . The association was specific to women and particularly evident with severe disc degeneration.
Meta-analysis across musculoskeletal conditions55 Meta-analysis across musculoskeletal conditions
Liu et al. (2018). Association between GDF5
rs143383 genetic polymorphism and musculoskeletal degenerative diseases susceptibility: a
meta-analysis. BMC Med Genet.:
Meta-analysis of GDF5 rs143383 polymorphism was statistically associated with increased risk of musculoskeletal degenerative diseases under each genetic model (allele model: OR = 1.32, 95% CI 1.19–1.48, P = 0.000; homozygote model: OR = 1.80, 95%CI 1.49–2.16, P = 0.000) , covering 5,915 cases and 12,252 controls across both osteoarthritis and intervertebral disc degeneration.
The effect sizes are modest but highly reproducible — classic for common variants affecting complex traits. The A allele doesn't guarantee joint problems, but it tips the scales toward faster cartilage degradation over decades.
Practical Implications
This variant influences your joints' capacity to maintain and repair cartilage throughout life. The A allele creates a slight ongoing deficit in GDF5 expression, which compounds with age, mechanical stress, and other risk factors.
Weight matters more for you. With reduced cartilage maintenance capacity, excess mechanical load accelerates degeneration. Each extra 5 kg of body weight increases knee OA risk, and this effect is amplified when your baseline cartilage repair is compromised.
Joint-protective nutrients and supplements may help compensate. Glucosamine and chondroitin66 Glucosamine and chondroitin
natural
components of cartilage
support cartilage structure. Studies show they can slow cartilage loss and reduce pain in OA, particularly
the glucosamine sulfate form combined with chondroitin. Omega-3 fatty acids (EPA/DHA) reduce inflammatory
responses in joints. Vitamin D and K support bone health underlying cartilage. SAM-e has shown cartilage-protective
effects and pain relief comparable to NSAIDs.
Activity patterns should favor joint preservation. Low-impact exercise (swimming, cycling, elliptical) maintains joint health without excessive wear. Strength training builds muscle support around joints, offloading cartilage. Avoid chronic high-impact activities and repetitive joint stress if possible.
Interactions
This effect is influenced by a second SNP (rs143384, C/T) in the same area . The two variants work together to regulate GDF5 expression through methylation:
The G alleles of both SNPs form CpG dinucleotides. Demethylation of both SNP's increases GDF5 expression . When you carry the A allele at rs143383 along with the A allele at rs143384, the reduction in GDF5 expression is most pronounced. This represents a compound heterozygosity scenario where the combined genotype creates a stronger effect than either variant alone.
Other genes in cartilage homeostasis pathways may also interact with GDF5 function, including COL2A1 (type II collagen), ACAN (aggrecan), and other BMP family members, though specific compound implications require individual research into those variants.
ANRIL at 9p21.3 — A Cell Cycle Gatekeeper in Endometriosis
The 9p21.3 region of chromosome 9 is one of the most functionally
complex loci in the human genome. It encodes three tumor suppressors
— p16INK4a, p14ARF (encoded by CDKN2A) and p15INK4b (encoded by CDKN2B)11 p16INK4a, p14ARF (encoded by CDKN2A) and p15INK4b (encoded by CDKN2B)
These proteins are cyclin-dependent kinase inhibitors that halt cell-cycle progression and trigger senescence
— and is also home to ANRIL (Antisense Non-coding RNA in the INK4 Locus)22 ANRIL (Antisense Non-coding RNA in the INK4 Locus)
Also designated CDKN2B-AS1 or CDKN2BAS — a long non-coding RNA transcribed on the antisense strand overlapping CDKN2A and CDKN2B,
a molecular rheostat for how readily cells stop dividing and become
senescent. rs1537377 is an independent endometriosis-risk signal at
this locus, distinct from but in partial linkage disequilibrium with
the index CDKN2BAS variants identified in Japanese GWAS.
The Mechanism
ANRIL regulates its neighboring tumor suppressors CDKN2A/B through
epigenetic silencing. It recruits the
PRC2 complex33 PRC2 complex
Polycomb Repressive Complex 2 — a histone methyltransferase that adds H3K27me3 marks, compacting chromatin and silencing genes
to deposit H3K27me3 on the CDKN2A/B promoters, and the PRC1 complex
to maintain repression through histone H2A ubiquitination. The net
effect is that ANRIL expression level determines how effectively
cells keep their p15/p16 senescence checkpoint suppressed — allowing
continued proliferation.
In endometriosis, this biology is directly relevant: ectopic
endometrial cells implant at peritoneal and ovarian sites and must
evade normal senescence and apoptosis to survive and proliferate in
an ectopic environment. Functional fine-mapping at the 9p21.3
endometriosis locus
showed that protective alleles drive stronger chromatin interaction with the ANRIL promoter44 showed that protective alleles drive stronger chromatin interaction with the ANRIL promoter
Borghese et al. 2016 — Allelic Imbalance in Regulation of ANRIL through Chromatin Interaction at 9p21 Endometriosis Risk Locus. PLOS Genetics.
via TCF7L2 and EP300 binding, increasing ANRIL transcription. More
ANRIL means more CDKN2A/B silencing, which counterintuitively appears
protective against endometriosis in this regulatory context, suggesting
that the disease-associated haplotype disrupts normal ANRIL-mediated
cell cycle regulation in endometrial stromal cells.
Consistent with this model, elevated ANRIL expression has been
independently found in ectopic endometriotic tissue relative to normal
ovarian epithelium
with expression levels increasing progressively from ASRM stage II to stage IV55 with expression levels increasing progressively from ASRM stage II to stage IV
Endometriosis Prognosis Correlates With Elevated Expression of LncRNA-ANRIL. Obstetrics & Gynecology International, 2025.,
suggesting that ANRIL dysregulation tracks disease severity rather
than simply marking disease presence.
The Evidence
The first GWAS to implicate CDKN2BAS in endometriosis was conducted in
1,907 Japanese women with endometriosis and 5,292 controls66 1,907 Japanese women with endometriosis and 5,292 controls
Uno et al. 2010 — A genome-wide association study identifies genetic variants in the CDKN2BAS locus associated with endometriosis in Japanese. Nature Genetics..
The lead SNP rs10965235 reached p=5.57×10⁻¹² (OR 1.44) and mapped to
intron 16 of CDKN2BAS. rs1537377 represents a distinct signal at this
locus identified in subsequent trans-ethnic meta-analysis.
In the
Nyholt et al. (2012) GWAS meta-analysis of 4,604 cases and 9,393 controls77 Nyholt et al. (2012) GWAS meta-analysis of 4,604 cases and 9,393 controls
Genome-wide association meta-analysis identifies new endometriosis risk loci. Nature Genetics.
across Japanese and European samples, rs1537377 at 9p21.3 emerged as
genome-wide significant in analyses restricted to moderate-to-severe
European cases (OR approximately 1.15, p<5×10⁻⁸). The larger
Rahmioglu et al. (2014) meta-analysis of 11,506 cases and 32,678 controls88 Rahmioglu et al. (2014) meta-analysis of 11,506 cases and 32,678 controls
Genetic variants underlying risk of endometriosis: insights from meta-analysis of eight GWAS and replication datasets. Human Reproduction Update.
confirmed the association at p=1.5×10⁻⁸, and noted that eight of nine
identified loci — including 9p21.3 — showed stronger effect sizes
among Stage III/IV cases, consistent with a role in facilitating
ectopic tissue invasion and persistence.
Belgian replication
in 998 cases and 783 controls99 in 998 cases and 783 controls
Sapkota et al. 2015 — Independent Replication and Meta-Analysis for Endometriosis Risk Loci. Twin Research and Human Genetics.
confirmed rs1537377 among nine loci retaining genome-wide significance
in European populations through meta-analysis, and noted that coding
variants within CDKN2B-AS1 near rs1537377 also showed nominally
significant associations.
Practical Actions
For women carrying one or two C alleles, the endometriosis risk signal at 9p21.3 is most actionable in two ways: earlier diagnostic vigilance for endometriosis (particularly moderate-to-severe presentations), and awareness that the 9p21.3 locus is shared with multiple other diseases including cardiovascular disease and several cancers, so any management of the locus's senescence pathway has broad implications.
The shared biology of 9p21.3 with ANRIL's role across diseases also means that interventions targeting cellular senescence — such as senolytic compounds that clear p16-positive cells — have theoretical relevance to the ectopic cell survival biology at this locus, though direct evidence in endometriosis specifically is not yet available.
Interactions
rs10965235 (CDKN2BAS): The index Japanese GWAS SNP at the 9p21.3 endometriosis locus. rs1537377 and rs10965235 tag the same chromosomal region but represent partially independent signals in trans-ethnic meta-analysis; combined carrier status would represent the full 9p21.3 endometriosis haplotype burden.
rs2811712 (CDKN2BAS): Another regulatory variant in ANRIL that primarily associates with physical aging and functional impairment through the same CDKN2B/p16 pathway. The two variants operate in the same ANRIL regulatory context but through different aspects of the locus.
rs1333049 (9p21.3 CAD locus): The primary coronary artery disease SNP at 9p21.3. Women with both the endometriosis-risk C allele at rs1537377 and the CAD-risk C allele at rs1333049 carry risk haplotypes at both faces of this pleiotropic locus, warranting combined cardiovascular and reproductive monitoring.
IL-1 Beta Promoter Variant — A Master Regulator of Inflammation
The IL1B gene encodes interleukin-1 beta (IL-1β), one of the most potent pro-inflammatory cytokines in the human body11 one of the most potent pro-inflammatory cytokines in the human body
IL-1β drives inflammation, activates immune cells, and plays a central role in atherosclerosis, sepsis, and autoimmune disease. The rs16944 variant sits in the promoter region at position -511, where it functions as a genetic dimmer switch controlling how much IL-1β your cells produce22 where it functions as a genetic dimmer switch controlling how much IL-1β your cells produce
The A allele is associated with higher IL-1β mRNA expression, while the G allele produces less.
This isn't just an academic curiosity — rs16944 influences your risk of sepsis, cardiovascular disease mortality, and inflammatory complications across dozens of conditions33 rs16944 influences your risk of sepsis, cardiovascular disease mortality, and inflammatory complications across dozens of conditions
From aspirin-induced asthma to coronary artery lesions in children, this variant shapes inflammatory outcomes.
The Mechanism
The -511 position in the IL1B promoter contains a binding site for transcription factors that regulate gene expression44 contains a binding site for transcription factors that regulate gene expression
The C-to-T change (G-to-A on the forward strand) alters the binding affinity of these regulatory proteins. The A allele creates a promoter configuration that permits higher transcription rates55 The A allele creates a promoter configuration that permits higher transcription rates
This results in elevated IL-1β mRNA and protein levels after immune stimulation.
IL-1β itself drives a cascade of inflammatory responses: it induces other cytokines like IL-6, activates endothelial cells to express adhesion molecules, promotes prostaglandin synthesis, and recruits immune cells to sites of inflammation66 drives a cascade of inflammatory responses: it induces other cytokines like IL-6, activates endothelial cells to express adhesion molecules, promotes prostaglandin synthesis, and recruits immune cells to sites of inflammation
This amplification loop means a small genetic change in IL-1β production gets magnified throughout the immune system.
Critically, IL-1β is the key output of the NLRP3 inflammasome77 IL-1β is the key output of the NLRP3 inflammasome
When danger signals like cholesterol crystals, uric acid, or pathogens activate this molecular complex, IL-1β is cleaved from its inactive precursor and released. Your rs16944 genotype determines how much raw material is available for this process.
The Evidence
The clearest evidence comes from sepsis studies. In 471 preterm infants, the AA genotype was significantly more common in those with early-onset sepsis (p=0.012) and was even more strongly associated with lethal outcomes (p=0.011)88 In 471 preterm infants, the AA genotype was significantly more common in those with early-onset sepsis (p=0.012) and was even more strongly associated with lethal outcomes (p=0.011). In adult sepsis, AA carriers showed higher mortality risk99 In adult sepsis, AA carriers showed higher mortality risk.
For cardiovascular disease, a 15-year follow-up of 2,010 Northern Ireland men found the A allele associated with increased all-cause mortality (HR 1.18, p=0.005)1010 a 15-year follow-up of 2,010 Northern Ireland men found the A allele associated with increased all-cause mortality (HR 1.18, p=0.005). The effect was dose-dependent: one A copy increased risk 18%, two copies increased it 43%1111 The effect was dose-dependent: one A copy increased risk 18%, two copies increased it 43%. This aligns with extensive evidence that IL-1β drives atherosclerosis progression1212 extensive evidence that IL-1β drives atherosclerosis progression
The CANTOS trial demonstrated that blocking IL-1β with canakinumab reduces cardiovascular events in high-risk patients.
In more specific inflammatory conditions, AA homozygotes show 2.98-fold increased risk of aspirin-exacerbated respiratory disease1313 AA homozygotes show 2.98-fold increased risk of aspirin-exacerbated respiratory disease, and in children under 12 months with Kawasaki disease, GG carriers (lower IL-1β producers) had significantly reduced risk of coronary artery lesions1414 in children under 12 months with Kawasaki disease, GG carriers (lower IL-1β producers) had significantly reduced risk of coronary artery lesions.
A meta-analysis found the variant associated with silent myocardial ischemia in diabetic patients1515 A meta-analysis found the variant associated with silent myocardial ischemia in diabetic patients
Under multiple inheritance models, the CC/CT genotypes (corresponding to AA/AG in forward orientation) increased risk with OR of 4.68 for homozygotes.
Practical Implications
If you carry one or two A alleles, you have a genetic predisposition to mount stronger IL-1β responses. This is a double-edged sword: potentially more effective at clearing infections initially, but prone to excessive inflammation that damages your own tissues1616 prone to excessive inflammation that damages your own tissues
Chronic low-grade inflammation accelerates atherosclerosis, increases thrombosis risk, and contributes to age-related disease.
The cardiovascular connection is particularly important. IL-1β induces IL-6 production, which drives hepatic synthesis of fibrinogen, plasminogen activator inhibitor, and C-reactive protein1717 IL-1β induces IL-6 production, which drives hepatic synthesis of fibrinogen, plasminogen activator inhibitor, and C-reactive protein
This shifts hemostasis toward a prothrombotic state while creating an inflammatory milieu that destabilizes atherosclerotic plaques. If you're AA and have existing cardiovascular risk factors, you're in a higher-risk category for events.
The sepsis association matters for surgical planning and critical illness. AA carriers may benefit from more aggressive infection monitoring and earlier intervention when signs of systemic inflammation appear.
Interactions
This variant sits within a tightly linked haplotype block with rs1143627 (-31C>T), another functional IL1B promoter SNP1818 tightly linked haplotype block with rs1143627 (-31C>T), another functional IL1B promoter SNP
The two variants are in nearly complete linkage disequilibrium, meaning they're usually inherited together. When evaluating IL-1β-related risk, consider both variants as a unit.
The IL-1 gene cluster on chromosome 2q13 also includes IL1A and IL1RN (encoding the IL-1 receptor antagonist). Variants in IL1RN can modulate the overall balance between pro- and anti-inflammatory signaling1919 Variants in IL1RN can modulate the overall balance between pro- and anti-inflammatory signaling
High IL-1Ra production may partially buffer the effects of high IL-1β.
From a pathway perspective, IL-1β functions upstream of many inflammatory cascades. Variants in downstream genes like IL6, TNF, and CRP may compound or mitigate the effects of rs16944 on disease risk.
The Nicotine Dependence Variant: How One Amino Acid Shapes Addiction Risk
The CHRNA5 gene encodes the alpha-5 subunit of the nicotinic acetylcholine receptor (nAChR), a critical component of the brain's response to nicotine. The rs16969968 variant11 The rs16969968 variant
This SNP is one of the most replicated genetic findings in addiction research, identified independently by multiple genome-wide association studies replaces aspartic acid with asparagine at position 398 of the alpha-5 subunit, fundamentally altering how your brain responds to nicotine. This single amino acid change has emerged as the strongest known genetic risk factor for heavy smoking and nicotine dependence.
The variant sits in the CHRNA5-CHRNA3-CHRNB4 gene cluster on chromosome 15q25, a region that has been consistently associated with smoking quantity22 consistently associated with smoking quantity
Meta-analyses report p-values as low as 5.57×10⁻⁷², making this one of the most significant genetic associations with any behavior across populations of European, Asian, African, and Latino ancestry. About 28% of Europeans carry at least one copy of the risk allele, though this frequency varies dramatically by population—only 2% of East Asians carry the A allele, while 15% of South Asians and Latinos do.
The Mechanism
The Asp398Asn substitution alters the structure of the alpha-5 subunit in a critical region called the second intracellular loop, which is highly conserved across species33 highly conserved across species
The aspartic acid at position 398 is preserved in mammals, suggesting strong evolutionary pressure to maintain this amino acid. When incorporated into α4β2α5 nicotinic receptors, the 398N (risk) variant reduces receptor function by approximately 50%44 reduces receptor function by approximately 50%
In vitro studies show the risk allele produces twofold lower maximal response to nicotinic agonists compared to the protective allele, measured by calcium influx and electrophysiological responses.
These α4β2α5 receptors are particularly abundant in the medial habenula and interpeduncular nucleus55 medial habenula and interpeduncular nucleus
These brain regions form a critical pathway that normally limits nicotine intake by generating aversive responses to high nicotine doses, structures that act as a natural brake on nicotine consumption. The reduced receptor function in risk allele carriers weakens this braking system. Behavioral studies confirm this: individuals with the AA genotype report significantly lower aversive effects from nicotine66 individuals with the AA genotype report significantly lower aversive effects from nicotine
In controlled intravenous nicotine administration, AA carriers rated nicotine as less unpleasant (P<5×10⁻⁸), with the effect most pronounced at higher doses, meaning they experience less nausea, dizziness, and discomfort that would normally discourage heavy smoking.
The Evidence
Multiple genome-wide association studies77 genome-wide association studies
The initial 2008 GWAS from three independent research groups all converged on the same chromosomal region identified rs16969968 as the most significant variant associated with smoking quantity. A 2010 meta-analysis pooling data from European populations88 meta-analysis pooling data from European populations
Analyzing over 140,000 individuals found rs16969968 with p=5.57×10⁻⁷² for cigarettes per day found that each A allele increases smoking quantity by approximately 1 cigarette per day, with an odds ratio of 1.3 for being a heavy smoker (≥20 cigarettes/day) versus a light smoker (≤10 cigarettes/day).
The effect extends across multiple populations99 multiple populations
A cross-ancestry meta-analysis found consistent OR=1.33 across European, Asian, and African populations, though effect sizes vary, though with varying effect sizes. In European Americans, each A allele confers OR=1.3 for nicotine dependence. In African Americans, where the allele is rarer (6% frequency), the effect is similar (OR=1.3) when present. Studies in Mexican populations1010 Mexican populations
Mexican cohort showed OR=3.12 for heavy smoking in AA carriers versus GG and Middle Eastern populations1111 Middle Eastern populations
Palestinian lung cancer cases showed 36.7% A allele frequency versus 17.5% in controls report even stronger associations, though these may reflect population-specific factors.
The variant also affects smoking cessation outcomes1212 smoking cessation outcomes
Meta-analyses show AA carriers have delayed time to cessation and lower success rates across multiple cessation interventions. Individuals with the AA genotype quit smoking later in life and have lower success rates with standard cessation interventions. One meta-analysis found the AA genotype is associated with a 7-year earlier lung cancer diagnosis1313 7-year earlier lung cancer diagnosis
Among lung cancer patients, AA carriers were diagnosed 7 years earlier on average (HR=0.68, p=4.9×10⁻¹⁰) in smokers who do develop lung cancer.
Beyond nicotine dependence, the variant shows a dose-dependent association with lung cancer risk1414 a dose-dependent association with lung cancer risk
Each A allele increases lung cancer risk with OR=1.3-1.6, even after adjusting for smoking quantity, with odds ratios ranging from 1.3 to 1.6 per risk allele in different studies. This association persists even after adjusting for smoking quantity, suggesting both a behavioral pathway (more smoking) and potentially a direct biological effect. The variant has also been linked to chronic obstructive pulmonary disease (COPD)1515 chronic obstructive pulmonary disease (COPD)
Mexican cohort showed OR=1.91 for COPD in A allele carriers, another smoking-related disease.
Interestingly, the same allele that increases nicotine dependence appears to be protective against cocaine dependence1616 protective against cocaine dependence
In two independent samples, the A allele showed OR=0.67 for cocaine dependence—the opposite direction from nicotine, with odds ratios in the opposite direction (OR=0.67). This paradox suggests the variant's effects are specific to nicotine rather than reflecting general addiction vulnerability.
Practical Implications
If you carry one or two copies of the A allele, understand that your brain's natural aversion to nicotine is blunted. You're genetically predisposed to find smoking less unpleasant than others, making it easier to escalate to heavy smoking and harder to quit. This is not a character flaw—it's neurobiology.
For smoking prevention, awareness matters most before starting. If you've never smoked and carry the AA genotype, you face approximately 4-fold increased risk1717 approximately 4-fold increased risk
Studies report ORs ranging from 1.9 for heterozygotes to 3.6 for AA homozygotes for developing nicotine dependence of developing severe nicotine dependence if you do start, compared to GG carriers. Even "social smoking" may escalate more quickly.
For current smokers with the risk genotype, cessation requires more intensive support. Standard approaches like nicotine replacement therapy1818 nicotine replacement therapy
Some evidence suggests NRT efficacy varies by genotype, though this remains under investigation and behavioral counseling may be insufficient. Evidence suggests varenicline (Chantix)1919 varenicline (Chantix)
Varenicline acts on α4β2 nAChRs and appears effective regardless of CHRNA5 genotype may be particularly effective for risk allele carriers because it works regardless of your CHRNA5 genotype, unlike some cessation aids whose efficacy varies by genetic background.
The lung cancer risk deserves serious attention. Even among smokers, those with the AA genotype develop lung cancer earlier2020 those with the AA genotype develop lung cancer earlier
Meta-analysis of 12,690 smokers found AA carriers diagnosed 7 years earlier on average and at younger ages. However, the same study found that quitting smoking reduces lung cancer risk equally2121 quitting smoking reduces lung cancer risk equally
Smoking cessation showed OR=0.48 for lung cancer across all genotypes across all genotypes—the protective effect of quitting does not depend on your CHRNA5 status. Quitting smoking cuts your lung cancer risk in half regardless of your genetic background.
Interactions
The CHRNA5 rs16969968 variant is in strong linkage disequilibrium with rs1051730 in the CHRNA3 gene (r²=1 in Europeans), meaning they are almost always inherited together and their effects are difficult to separate. The 15q25 region also contains CHRNB4 and several other variants (rs588765, rs578776, rs6495309) that show associations with smoking behaviors, though rs16969968 appears to be the primary functional variant based on in vitro studies and cross-population analyses.
The variant's effect on smoking behavior is modified by age of smoking initiation2222 modified by age of smoking initiation
Meta-analysis found the genetic effect is strongest in those who started smoking before age 16, with stronger genetic effects in those who start smoking during adolescence. Environmental factors like childhood adversity and peer smoking also interact with the genotype to influence dependence risk, though the specific mechanisms remain under investigation.
While rs16969968 is the strongest single genetic predictor of nicotine dependence, it explains only a small fraction of overall addiction risk. Smoking behavior is highly polygenic, with dozens of additional variants across the genome contributing smaller effects. Your genotype at this SNP should inform risk assessment and treatment planning, but it does not determine your destiny.
FADS1 rs174541 — Delta-5 Desaturase Depth
Your body's ability to build long-chain omega-3 and omega-6 fatty acids from
dietary precursors hinges on a single enzyme: delta-5 desaturase11 delta-5 desaturase
FADS1
(Fatty Acid Desaturase 1) — the enzyme that adds a double bond at the fifth
carbon position, converting DGLA to arachidonic acid in the omega-6 pathway and
eicosatetraenoic acid (ETA) to EPA in the omega-3 pathway.
rs174541 is an intronic variant in the FADS gene cluster on chromosome 11q12.2
that acts as an independent regulator of how much FADS1 enzyme your cells make.
The C allele dampens FADS1 expression, reducing the throughput of both the
omega-6 pathway (less arachidonic acid from dietary linoleic acid) and the
omega-3 pathway (less EPA from plant-derived ALA). Because the C allele also
affects circulating triglyceride concentrations, this variant sits at the
junction between fatty acid metabolism and broader cardiometabolic risk.
The Mechanism
rs174541 sits within FADS1's intronic regulatory architecture, in high linkage disequilibrium with the established functional cluster of FADS1 variants (rs174546, rs174547, rs174548, rs174537). Intronic variants in this region influence transcription factor binding sites and enhancer elements between FADS1 and FADS2 — the C allele at rs174541 tracks with reduced FADS1 mRNA levels and lower delta-5 desaturase enzyme activity across liver and blood cells.
The functional consequence plays out across two metabolic pathways simultaneously. In the omega-6 pathway: dietary linoleic acid (LA) → GLA → DGLA → [delta-5 desaturase] → arachidonic acid (AA). Reduced FADS1 activity means more DGLA accumulates and less AA is produced. In the omega-3 pathway: plant ALA → stearidonic acid → ETA → [delta-5 desaturase] → EPA. Again, the rate-limiting step is impaired, meaning less EPA is synthesised from the plant-sourced precursor. Neither pathway can compensate for the other — both require the same enzyme.
The Evidence
A Bayesian quantitative trait nucleotide analysis22 Bayesian quantitative trait nucleotide analysis
Voruganti et al. Variants in
CPT1A, FADS1, and FADS2 are Associated with Higher Levels of Estimated Plasma and
Erythrocyte Delta-5 Desaturases in Alaskan Eskimos. Front Genet,
2012 in 761 Alaskan Eskimos assigned
rs174541 a posterior probability >0.8 for functional effect on erythrocyte delta-5
desaturase activity — the strongest statistical evidence available in a Bayesian
framework for a variant being causally linked to its phenotype, not merely in LD
with the causal site.
The broader FADS1 locus evidence is overwhelming. A landmark GWAS in the
InCHIANTI study33 landmark GWAS in the
InCHIANTI study
Tanaka et al. Genome-wide association study of plasma
polyunsaturated fatty acids in the InCHIANTI Study. PLoS Genet,
2009 of 1,075 Italian adults found
that FADS1 cluster variants explain 18.6% of all additive variance in circulating
arachidonic acid (p=5.95×10⁻⁴⁶), by far the largest explained variance for any
common PUFA-metabolism variant. The CHARGE Consortium meta-analysis44 CHARGE Consortium meta-analysis
Lemaitre
et al. Genetic loci associated with plasma phospholipid n-3 fatty acids. PLoS Genet,
2011 across 8,866 European ancestry
participants confirmed that FADS1 minor alleles are the dominant genetic predictor
of lower circulating EPA (p=5×10⁻⁵⁸) and higher plant ALA (p=3×10⁻⁶⁴).
Beyond fatty acid ratios, FADS1 variants have direct consequences for clinical lipid
panels. In 21,004 Japanese individuals55 21,004 Japanese individuals
Nakayama et al. A single nucleotide
polymorphism in the FADS1/FADS2 gene is associated with plasma lipid profiles.
Hum Genet,
2010, the C allele at the tightly
linked rs174547 was significantly associated with higher triglycerides
(p=1.5×10⁻⁶) and lower HDL-C (p=0.03), demonstrating that FADS1 reduced activity
has effects visible on a standard fasting lipid panel — not just on specialised
PUFA measurements. A study of 8,842 Korean adults66 study of 8,842 Korean adults
Lee et al. Functional Impact
of the FADS1 rs174546 Single Nucleotide Polymorphism on Serum Lipid Levels. Mol Nutr
Food Res,
2024 quantified this: the FADS1
minor allele increases fasting serum triglycerides by 6.48 ± 1.84 mg/dL per allele,
mediated through reduced LC-PUFA production and downstream effects on VLDL
assembly and clearance.
Practical Actions
For C allele carriers, the core problem is that dietary plant-based omega-3 sources (flaxseed, chia, walnuts, canola oil) supply ALA which requires FADS1 to reach EPA. When FADS1 activity is reduced, that conversion chain slows — the ALA enters the bloodstream but stalls before reaching EPA. Direct supplementation with preformed EPA and DHA from marine or algae-based sources entirely bypasses the impaired step. The dosage depends on genotype: CT carriers benefit from 1–2 g EPA+DHA daily; CC homozygotes need 2–4 g to overcome the more complete impairment.
The triglyceride finding adds a monitoring dimension. If you carry the C allele and have other cardiovascular risk factors, a fasting lipid panel captures the clinical footprint of this genotype — elevated triglycerides and reduced HDL are the measurable downstream signal of impaired FADS1 activity.
Interactions
rs174541 is in high linkage disequilibrium (r² >0.8) with the established FADS1 functional cluster including rs174547, rs174548, rs174546, and rs174537. Users may carry risk alleles at multiple sites on the same haplotype — carrying the C allele at rs174541 in combination with risk alleles at rs174548 and rs174547 increases the probability of being on the full reduced-function FADS1 haplotype. No additional effect beyond what each individual variant predicts is needed for interpretation, as the variants tag the same underlying expression phenotype.
The ELOVL2 gene variant rs17606561 (also in the platform) encodes elongase 2, which converts EPA to DHA. A user who carries both FADS1 reduced-activity alleles and an ELOVL2 impairment faces a double block in the ALA → EPA → DHA pathway, potentially producing the most severe DHA deficiency of any single-enzyme genotype combination. For such users, a DHA-specific supplement target (≥500 mg DHA per day) matters beyond total EPA+DHA.
CXCL12 — The Chemokine That Guides Inflammatory Cells Into Coronary Arteries
Beneath the surface of every coronary artery, a biochemical signaling system
continuously directs immune cells to sites of injury and repair. At the center of
this system is CXCL1211 CXCL12
C-X-C motif chemokine ligand 12, also known as stromal
cell-derived factor 1 (SDF-1), a
chemokine that acts as a powerful attractant for progenitor cells, natural killer
cells, and immune progenitors via its receptor CXCR4. The rs1746048 variant lies
approximately 80 kilobases downstream of the CXCL12 gene on chromosome 10q11.21
and is one of the earliest and most consistently replicated genome-wide significant
loci for coronary artery disease (CAD). It was among the first loci to reach
genome-wide significance in large GWAS meta-analyses and has been replicated
across European, East Asian, and South Asian populations.
The Mechanism
The rs1746048 variant does not alter the CXCL12 protein directly. It is an
intergenic regulatory variant that modulates how much CXCL12 protein circulates
in the bloodstream. Carriers of the risk C allele have demonstrably higher plasma
CXCL12 levels than T allele carriers.
Mehta et al. (2011)22 Mehta et al. (2011)
The novel atherosclerosis locus at 10q11 regulates plasma CXCL12 levels, Eur Heart J 2011
showed a dose-dependent relationship: CC carriers had the highest CXCL12 levels,
CT intermediate, and TT the lowest (2.33 vs 2.27 vs 2.21 ng/mL, P=0.034), and
identified elevated CXCL12 transcript expression in natural killer cells and liver
tissue as the downstream effectors.
Elevated CXCL12 promotes atherosclerosis through several mechanisms: it
accelerates the migration and retention of hematopoietic progenitor cells33 hematopoietic progenitor cells
bone
marrow-derived precursor cells including monocyte precursors
into arterial walls, increases the homing of inflammatory immune cells to
developing plaques, and drives neovascularization44 neovascularization
formation of new blood vessels
within plaques, which are structurally fragile and prone to hemorrhage
of atherosclerotic lesions. Additionally, a 2017 Pakistani study found that the C
risk allele is associated with elevated IL-18 and a higher IL-18:IL-10 pro-inflammatory
ratio, pointing toward dysregulated cytokine balance as a downstream consequence of
elevated CXCL12 signaling.
The variant's effect on carotid intima-media thickness — a direct measure of subclinical atherosclerosis — has been quantified. A 2015 meta-analysis found that each C allele is associated with approximately 0.008 mm thicker carotid walls, suggesting the variant promotes vascular remodeling independent of established lipid risk factors.
The Evidence
The evidence base for rs1746048 is substantial and unusually consistent across
ethnic groups. The most comprehensive meta-analysis,
Chen et al. (2017)55 Chen et al. (2017)
Meta-analysis pooling 48,852 CHD patients and 64,386 controls from 12 studies, Medicine (Baltimore),
found per-C-allele odds ratios of 1.07 in Asian populations and 1.14 in Caucasian
populations (1.11 overall). A concurrent meta-analysis by
Huang et al. (2013)66 Huang et al. (2013)
Case-control study in Han Chinese plus meta-analysis of >107,000 individuals, Gene
confirmed an OR of 1.12 (P<0.0001), with subgroup analysis showing the effect is
amplified in older patients (OR 1.91 for age ≥65) and males under a recessive model
(OR 1.72 for CC+CT vs TT).
The longitudinal evidence is equally convincing.
Wirtwein et al. (2017)77 Wirtwein et al. (2017)
1,345 confirmed CAD patients with 8.6-year follow-up, Int J Cardiol
showed that rs1746048 predicted both the need for revascularization procedures and
major adverse cardiovascular events (MACE), indicating that the variant does not
merely associate with incident CAD but also with its clinical severity and
progression over time.
Practical Actions
The rs1746048 C allele is common — approximately 68% of people globally carry CC. Its per-allele effect is modest (OR ~1.12 per C allele), but it is one of the more robustly replicated GWAS hits for CAD, and its effect through elevated CXCL12 levels means it operates through an inflammatory pathway distinct from the standard lipid-driven mechanisms targeted by statins.
The T allele is protective and dose-dependent: CT carriers have measurably lower plasma CXCL12 levels and modestly lower CAD risk than CC carriers, and TT homozygotes have the lowest CXCL12 levels of all. For CC homozygotes, the key clinical implication is earlier subclinical atherosclerosis screening — particularly carotid intima-media thickness measurement and coronary artery calcium scoring — since the mechanistic evidence points directly to arterial wall thickening as a downstream effect.
Interactions
The strongest documented interaction is with rs501120, a neighboring variant also on chromosome 10q11.21 and also downstream of CXCL12. Mehta et al. showed that rs1746048 and rs501120 are both independently associated with CAD risk and both regulate plasma CXCL12 levels, suggesting that the 10q11.21 locus contains multiple functional regulatory elements rather than a single causal variant. Combined genotype data for these two SNPs has not been formally analyzed in a compound heterozygosity framework.
rs1746048 also co-occurs in multi-locus CAD genetic risk scores alongside rs4977574 (CDKN2B-AS1 / 9p21 locus) and rs1333049 (another 9p21 variant). The 9p21 locus and the 10q11.21 locus operate through distinct mechanisms — cell cycle regulation versus chemokine-driven inflammation respectively — meaning that carrying risk alleles at both loci compounds risk through non-overlapping biological pathways.