MT-ND5 m.13042G>A — A Mitochondrial Throttle Stuck in Low Gear

Every cell in your body runs on ATP — the molecular currency of energy — manufactured by the electron transport chain11 electron transport chain
a series of five protein complexes embedded in the inner mitochondrial membrane
. Complex I (NADH:ubiquinone oxidoreductase) is the first and largest of these complexes, accepting electrons from cellular metabolism and using their energy to pump protons that ultimately drive ATP synthesis. MT-ND5 encodes one of the seven core subunits of complex I that are encoded in the mitochondrial genome itself rather than in the cell nucleus. The m.13042G>A variant substitutes alanine for threonine at a position (Ala236) that is conserved across vertebrates and sits within a functionally critical domain of the ND5 subunit.

Unlike nuclear DNA, where you have two copies (one from each parent), you carry hundreds to thousands of copies of mitochondrial DNA in every cell. This variant is heteroplasmic — meaning some copies carry the G>A mutation while others are normal. The ratio of mutant to normal copies (the heteroplasmy level) is the principal determinant of whether, and how severely, the variant causes disease. At low heteroplasmy levels (<50-60%), most individuals are asymptomatic or only mildly affected. Above 70-80%, complex I activity falls sufficiently to produce clinical disease. Heteroplasmy levels can differ substantially between tissues and between generations.

The Mechanism

The p.Ala236Thr substitution replaces a non-polar alanine with a larger, polar threonine at a position under strong evolutionary constraint in the ND5 protein. The consequence is reduced complex I assembly and/or catalytic efficiency. When mutant mtDNA exceeds the biochemical threshold in a tissue, that tissue cannot meet its ATP demand — and the most energy-intensive tissues (brain, skeletal muscle, heart, retina) fail first.

The phenotypic range is broad because heteroplasmy levels vary between individuals and tissues. In the first reported family, affected members presented with LHON-like optic neuropathy, retinopathy, cataracts, strokes, and early deaths — all tracing through the maternal line22 LHON-like optic neuropathy, retinopathy, cataracts, strokes, and early deaths — all tracing through the maternal line
Valentino ML et al. The 13042G>A/ND5 mutation in mtDNA is pathogenic and can be associated also with a prevalent ocular phenotype. J Med Genet, 2006
. The original case report described a 25-year-old man with MELAS/MERRF overlap (strokes, seizures, myoclonus) and confirmed isolated complex I deficiency in muscle biopsy.

The Evidence

The pathogenicity of m.13042G>A is supported by multiple independent case series and by expert panel review. ClinVar VCV00000970333 ClinVar VCV000009703
ClinVar expert panel review — Likely Pathogenic, ★★★★ review status
classifies the variant as Likely Pathogenic for mitochondrial disease, based on six unrelated individuals with disease onset from the first year of life to adulthood, variable phenotypes including CPEO, LHON, and Leigh-like syndrome, and demonstrated disease segregation within families. The criteria met include PS4_moderate (case enrichment), PP1_moderate (co-segregation), PP3 (computational evidence), and PM2_supporting (rare in population databases).

Naini et al.44 Naini et al.
Naini AB et al. A novel heteroplasmic point mutation in the mitochondrial tRNA gene. Arch Neurol, 2005
established the first clinical description and biochemically confirmed complex I deficiency in muscle. Blok et al.55 Blok et al.
Blok MJ et al. Mutations in the ND5 subunit of complex I of the mitochondrial DNA are a frequent cause of oxidative phosphorylation disease. J Med Genet, 2007
screened 116 complex I–deficient patients and identified m.13042G>A in a patient with Leigh-like syndrome, concluding that ND5 represents a diagnostic hot spot. Danhelovska et al.66 Danhelovska et al.
Danhelovska T et al. Multisystem mitochondrial diseases due to mutations in mtDNA-encoded subunits of complex I. BMC Pediatrics, 2020
showed in a 13-patient MT-ND cohort that higher heteroplasmy correlates with worse clinical outcomes, though biochemical complex I activity does not reliably track mutation load — making heteroplasmy quantification more diagnostically informative than enzyme assay alone.

Practical Actions

Management of mitochondrial complex I disease focuses on supporting residual mitochondrial function, avoiding metabolic stressors, and monitoring affected organs. Supplementation with mitochondrial cofactors is standard practice in specialist centres, though robust RCT evidence is limited by disease rarity. The GeneReviews MELAS protocol recommends CoQ10 (adults 200–400 mg/day in 3 divided doses, children 5–10 mg/kg/day), L-carnitine (adults 3 g/day, children 100 mg/kg/day), and creatine for individuals with this class of mitochondrial disease. Energy demand from intercurrent illness, fasting, extreme exertion, and valproate can precipitate acute metabolic decompensation.

Heteroplasmy measurement — ideally in muscle or urine (more informative than blood) — guides prognosis and should be repeated at clinical reassessments. Organ-specific surveillance includes annual cardiac, ophthalmological, and audiological evaluation, given the multi-organ tropism of ND5 mutations.

Interactions

All mitochondrial disease variants share the same metabolic pathway (oxidative phosphorylation) and interact through the common mechanism of ATP depletion and reactive oxygen species excess. Combination of m.13042G>A with nuclear-encoded complex I subunit variants (NDUFS1, NDUFS2, NDUFV1, etc.) in the same individual is theoretically additive, though clinical documentation of such combinations is very limited given the rarity of the variant. The maternal inheritance pattern means that risk to relatives flows exclusively through the maternal lineage — all children of an affected or carrier mother are at risk.

SH2B3 R262W — The Pleiotropic Inflammation Dial

SH2B3 (also known as LNK) encodes an adaptor protein that acts as a negative regulator of cytokine and growth-factor signalling11 negative regulator of cytokine and growth-factor signalling
SH2B3 contains a pleckstrin homology domain and an SH2 domain; it suppresses JAK2, c-kit, MPL, and IL-7R signalling by competing with downstream effectors at phosphotyrosine docking sites
. The rs3184504 variant swaps arginine for tryptophan at position 262 in the pleckstrin homology domain, partially disabling this brake. The result is a hypomorphic allele22 hypomorphic allele
a variant that reduces but does not abolish protein function
— the protein is still made, but it represses downstream signalling less effectively.

This makes rs3184504 one of the most studied pleiotropic SNPs in the human genome: a single amino acid change33 single amino acid change
R262W located in exon 3 of SH2B3 on chromosome 12q24
that influences blood pressure, platelet count, eosinophil count, coronary artery disease, type 1 diabetes, and celiac disease — each replicated across multiple large GWAS and validated mechanistically.

The Mechanism

SH2B3 normally dampens the proliferation and activation of hematopoietic progenitor cells and T cells. The R262W change reduces SH2B3's ability to bind membrane phosphoinositides and to suppress JAK/STAT and ERK1/2 cascades44 JAK/STAT and ERK1/2 cascades
downstream of cytokine receptors including thrombopoietin receptor MPL, stem-cell factor receptor c-kit, and IL-12 receptor
. In blood cells, reduced SH2B3 function leads to expansion of hematopoietic stem cells and enhanced megakaryopoiesis — more platelets and leukocytes. In T cells, the Trp variant is less repressive of IL-12 signalling55 less repressive of IL-12 signalling
IL-12 triggers Stat4 phosphorylation and IFNγ production; Trp-SH2B3 suppresses this less effectively than Arg-SH2B3
, leading to greater IFNγ output from CD8+ T cells when stimulated.

Knockin mouse studies66 Knockin mouse studies
Mice engineered to carry the equivalent Trp substitution exhibit a phenotype similar to full Sh2b3 knockout, confirming R262W is at least a partial loss of function
confirm the mechanism: Trp/Trp animals develop approximately 10 mmHg higher systolic blood pressure under chronic angiotensin II challenge, alongside greater renal infiltration by CD8+ T cells, perivascular fibrosis, and albuminuria — a cardiovascular–immune axis linking the variant to hypertension and end-organ damage.

The Evidence

Blood pressure. The Global BPgen consortium77 Global BPgen consortium
34,433 Europeans tested for 2.5 million SNPs; SH2B3 reached p = 3×10−18 for diastolic blood pressure
identified SH2B3 rs3184504 as one of eight genome-wide-significant blood pressure loci in 2009. The T allele is associated with modestly elevated diastolic and systolic blood pressure across large European cohorts.

Coronary artery disease. A PRISMA-compliant meta-analysis88 PRISMA-compliant meta-analysis
12 studies, 25,845 CHD cases and 68,910 controls
pooling 25,845 cases and 68,910 controls found OR = 1.12 (95% CI 1.09–1.15) for the T allele overall; OR = 1.13 (1.08–1.18) for myocardial infarction. The association was not significant in Asian populations, consistent with the T allele being extremely rare outside European ancestry.

Cardiac inflammation and fibrosis. Sh2b3-knockout rats develop 2.2-fold greater post-MI fibrosis99 2.2-fold greater post-MI fibrosis
measured by collagen staining and echocardiographic fractional shortening in Sh2b3-KO vs wild-type rats after myocardial infarction
and significantly more leukocyte infiltration than wild-type controls, with impaired left-ventricular function — evidence that SH2B3 actively restrains cardiac inflammation after injury.

Atherosclerosis and thrombosis. Using human cord blood, the common TT genotype showed expansion of hematopoietic stem cells1010 expansion of hematopoietic stem cells
increased MPL/thrombopoietin signalling, enhanced platelet production and activation, leukocytosis
and enhanced megakaryopoiesis; in hypercholesterolaemic mice, LNK deficiency exacerbated both atherosclerotic plaque burden and thrombosis. Eosinophil-specific LNK deficiency promotes arterial thrombosis independently, mediated by eosinophil degranulation onto vessel walls1111 mediated by eosinophil degranulation onto vessel walls
Blood 2024
.

Autoimmune disease. In 2008 rs3184504 was identified as a celiac disease risk locus1212 identified as a celiac disease risk locus
Nature Genetics study of 2,410 cases; rs3184504 T allele OR ~1.18
in the first post-HLA celiac GWAS. A 2008 NEJM study confirmed it as one of seven loci shared between type 1 diabetes and celiac disease1313 seven loci shared between type 1 diabetes and celiac disease
Smyth et al., NEJM 2008, 8,064 T1D patients and 2,828 trios
. The same T allele that raises cardiovascular risk also raises autoimmune risk — a genuine trade-off driven by the same underlying mechanism: less immune suppression, more immune activation.

Longevity trade-off. In a UK Biobank PheWAS of 379,758 Europeans1414 UK Biobank PheWAS of 379,758 Europeans
phenome-wide association study of 52 aging traits, genotype in Hardy–Weinberg equilibrium (p=0.642)
, CC homozygotes were 18% more likely to have extremely long-lived parents (OR 1.18, 95% CI 1.07–1.29), with lower blood pressure, fewer cardiovascular events, and lower T1D and hypothyroidism rates. TT homozygotes had modest cancer protection (any cancer OR 0.97) — a biological trade-off between cardiovascular longevity and immune surveillance.

Practical Actions

The T allele increases cardiovascular risk primarily through two channels: elevated blood pressure and a mildly prothrombotic haematological profile (higher platelet counts, higher eosinophils). For TT carriers especially, blood pressure monitoring is the most directly actionable intervention — even a 10 mmHg difference in systolic BP is clinically meaningful for lifetime cardiovascular risk.

On the autoimmune side, the same T allele slightly elevates risk for celiac disease and type 1 diabetes. This is worth knowing if you have unexplained gastrointestinal symptoms or a family history of autoimmune disease.

The dietary pattern most aligned with the biology of this variant is one that specifically targets platelet activation and blood pressure: omega-3 fatty acids reduce platelet aggregation and lower blood pressure, and magnesium and potassium support vasodilation. These are genotype-specific rather than generic because the mechanism is elevated platelet production and impaired T-cell braking — not just general cardiovascular risk.

Interactions

SH2B3 rs3184504 sits in a region of strong linkage disequilibrium on chromosome 12q24 together with ATXN2 and BRAP. Many published studies cannot cleanly separate SH2B3 from ATXN2 associations at this locus; however, functional and knockin evidence supports SH2B3 as the causal gene. The autoimmune associations (PTPN22 rs2476601, CTLA4 rs3087243, HLA-DQA1 rs2187668) operate through partially overlapping T-cell activation pathways. A carrier with both rs3184504 T and rs2476601 A (PTPN22 R620W) has two independent lesions in T-cell braking — one reducing SH2B3 suppression of JAK/STAT, the other disrupting LYP-mediated TCR downregulation. The combined effect on autoimmune risk would be expected to be additive or supra-additive, though direct compound evidence is limited to observational co-occurrence data rather than interaction studies.

Carriers of both rs3184504 TT and rs2476601 AG/AA have two independent T-cell activation-amplifying variants and would benefit from combined autoimmune monitoring — including baseline ANA, RF, anti-CCP, and anti-TPO — given the convergent biology.

PTPN22 R263Q — The Protective Brake on Immune Hyperactivation

The PTPN22 gene encodes lymphoid tyrosine phosphatase (LYP), a master brake on T-cell and B-cell activation. While the well-known R620W variant (rs2476601) increases autoimmune risk by disrupting regulatory interactions, the R263Q variant works through an entirely different mechanism — it reduces the catalytic phosphatase activity11 reduces the catalytic phosphatase activity
R263Q is a loss-of-function missense substitution in the PTPN22 catalytic domain, reducing enzymatic dephosphorylation of key T-cell signaling proteins
of the enzyme itself. Carriers of the T allele (representing the Q263 amino acid) have a somewhat weaker LYP phosphatase, which paradoxically helps maintain better immune tolerance and reduces risk for several autoimmune conditions.

The Mechanism

Position 263 lies within the catalytic domain of LYP, where arginine-263 forms a stabilizing contact with glutamine-34. The R263Q substitution — arginine to glutamine — disrupts this contact, displacing the α2′ helix22 displacing the α2′ helix
Crystal structure analysis showed displacement of the α2′ helix in R263Q, resulting in an altered substrate-binding cleft with open WPD-loop configuration versus the half-closed wild-type conformation
and reconfiguring the substrate-binding pocket. The resulting protein retains its basic fold but has reduced ability to dephosphorylate its physiological substrates — primarily the Src-family kinase Lck, ZAP-70, and the TCR-zeta chain. With a less active phosphatase, T-cell receptor signaling proceeds with somewhat less suppression, supporting more robust responses to antigens including self-antigens at tolerance-induction checkpoints.

Crucially, this mechanism is independent of and structurally distinct from R620W. Where R620W sits in the C-terminal proline-rich P1 motif and disrupts CSK binding, R263Q sits directly in the catalytic domain and impairs enzymatic function. Genetic studies confirm complete independence33 Genetic studies confirm complete independence
The R263Q and R620W polymorphisms show no linkage disequilibrium and their autoimmune associations are statistically independent
. Notably, all non-human species examined encode glutamine at position 263 — suggesting the ancestral Q263 form is the evolutionarily conserved variant, and the common R263 (C allele) represents the derived human-specific change.

The Evidence

The protective effect of Q263 (T allele) has been replicated across multiple disease contexts and populations. The founding study by Orrú et al. 200944 Orrú et al. 2009
A loss-of-function variant of PTPN22 is associated with reduced risk of systemic lupus erythematosus. Hum Mol Genet, 2009
demonstrated SLE protection across four cohorts (Spain, Italy, Argentina, US) totaling 2,093 patients and 2,348 controls, with a pooled OR of 0.63 (95% CI 0.47–0.84, p=0.0017). The Q263 allele frequency was approximately 1.7–3.4% across cohorts, consistent with a rare but functionally important variant.

For rheumatoid arthritis, Rodríguez-Rodríguez et al. 201155 Rodríguez-Rodríguez et al. 2011
The PTPN22 R263Q polymorphism is a risk factor for rheumatoid arthritis in Caucasian case-control samples. Arthritis Rheum, 2011
studied 10,971 participants across six European countries and showed the T allele associated with reduced RA risk (OR 0.80, 95% CI 0.67–0.96, p=0.016). Additive analysis confirmed that carrying fewer T alleles at R263Q compounded RA risk: individuals homozygous for the common C allele at R263Q and the risk allele at R620W faced the highest combined risk.

For inflammatory bowel disease, Diaz-Gallo et al. 201166 Diaz-Gallo et al. 2011
Differential association of two PTPN22 coding variants with Crohn's disease and ulcerative colitis. Inflamm Bowel Dis, 2011
showed a striking disease-specific pattern in 6,691 participants: R263Q protected against ulcerative colitis (meta-analysis OR 0.69, 95% CI 0.51–0.93) but showed no association with Crohn's disease. This dissociation — combined with R620W's opposite pattern (associated with Crohn's but not UC) — underscores that the two PTPN22 coding variants affect distinct immunological pathways despite sharing the same gene.

Notably, R263Q does not appear to protect against type 1 diabetes, Graves' disease, systemic sclerosis, or uveitis — conditions where R620W has well-documented associations. This specificity may reflect differences in the immune cell types and signaling thresholds relevant to each disease.

Practical Implications

Carrying one or two copies of the T allele (CT or TT genotype) is associated with modestly lower baseline risk for SLE, RA, and ulcerative colitis. This is a meaningful genetic buffer, particularly for individuals with family history of these conditions. The protective effect is most clearly documented in European-ancestry populations; the variant is essentially absent in East Asian and African populations, so the evidence base largely applies to Europeans and Hispanics.

Reduced genetic risk does not eliminate risk entirely — autoimmune diseases are polygenic and heavily influenced by environmental triggers, HLA alleles, hormonal factors, and microbiome composition. The T allele provides a partial protective signal, not immunity.

Interactions

R263Q and R620W (rs2476601) are the two characterized functional coding variants in PTPN22, and their effects are additive. The combination of CC at R263Q (no Q263 protection) with the risk genotype at R620W yields the highest autoimmune risk. Rodríguez-Rodríguez et al. 201177 Rodríguez-Rodríguez et al. 2011
Additive PTPN22 risk allele analysis: 2 risk alleles OR 1.28, 3 risk alleles OR 2.01, 4 risk alleles OR 3.55
quantified this: carrying four risk alleles (CC at R263Q + TT at R620W) raised RA risk 3.5-fold compared to baseline. Conversely, carrying the T allele at R263Q may partially offset R620W risk.

The mechanistic independence of the two variants — one in the catalytic domain, one in the regulatory P1 motif — means they affect different aspects of LYP function. Neither variant is in linkage disequilibrium with the other, so each must be assessed individually.

MYH7 Arg663His — The Sarcomere Mutation That Triggers Hypertrophic Cardiomyopathy and Atrial Fibrillation

Beta-myosin heavy chain (MYH7) is the molecular engine of the heart. Together with actin filaments, it generates the forceful contraction that drives blood out of the left ventricle with every beat. The Arg663His variant — a substitution of arginine for histidine at position 663 of the protein — sits directly in the myosin motor domain11 myosin motor domain
The ATPase-powered globular head region where myosin physically binds actin and cycles through the power stroke
, the precise region where force is generated. ClinVar classifies this variant as pathogenic, reviewed by an expert panel, and it has been documented in more than 30 unrelated individuals with confirmed hypertrophic cardiomyopathy (HCM).

The Mechanism

At codon 663, arginine — a large, positively charged amino acid — normally interacts with the actin filament at the myosin-actin interface during the cross-bridge cycle. Replacing it with histidine alters the local charge environment and changes how tightly myosin binds actin. The consequence, demonstrated in patient-specific iPSC-derived cardiomyocytes22 patient-specific iPSC-derived cardiomyocytes
A ten-member family cohort carrying Arg663His; iPSC-CMs fully recapitulated HCM pathology in vitro, Lan et al. Cell Stem Cell 2013
, is dysregulated calcium cycling: resting intracellular Ca²⁺ is elevated, Ca²⁺ transients are abnormal, and cardiomyocytes develop contractile arrhythmias at the single-cell level. Crucially, pharmacological restoration of Ca²⁺ homeostasis prevented development of cellular hypertrophy and electrical irregularities in these models — identifying calcium dysregulation as a causal, not merely correlating, mechanism.

The structural consequence at the organ level is asymmetric left ventricular hypertrophy33 asymmetric left ventricular hypertrophy
Predominantly in the proximal interventricular septum, though distribution varies across carriers in the same family
. Unlike MYBPC3-related HCM, which tends to manifest later in life, MYH7 variants cause measurable LVH at younger ages, and the wall thickening seen with Arg663His was documented to increase approximately 40% over a 7-year follow-up — a slow but progressive course.

The Evidence

The foundational clinical study was published by Gruver, Fatkin, Seidman, and colleagues in 199944 Gruver, Fatkin, Seidman, and colleagues in 1999
24-individual kindred; 47% AF prevalence among adults with ventricular hypertrophy; p<0.001 vs. ungenotyped HCM populations
. This family cohort established two distinctive features of Arg663His HCM: near-normal overall survival, and an extraordinarily high atrial fibrillation (AF) rate. Nearly half of adult carriers with ventricular hypertrophy developed AF — a rate far exceeding what is seen in typical HCM.

This AF predisposition was independently confirmed in a multinational HCM registry of 1,040 genotype-positive patients55 multinational HCM registry of 1,040 genotype-positive patients
Mean follow-up 7.2 years; adjusted for age, sex, proband status, LA size, wall thickness, and peak gradient; Lee et al. Circ Heart Failure 2018
. Among all MYH7 pathogenic variant carriers, the hazard ratio for new-onset AF was 1.7 (95% CI 1.1–2.6) compared to MYBPC3 carriers — the highest of any genotype group — after adjusting for established clinical AF risk factors. MYH7 variant carriers also showed earlier disease onset, greater left atrial dilation, and more extensive myocardial fibrosis, all of which structurally predispose to arrhythmia initiation.

Codon 663 is considered a mutational hotspot66 mutational hotspot
A different variant at the same codon, Arg663Cys (rs193922376), is independently pathogenic; 15% of MYH7-positive HCM patients in one cohort had a variant at this codon, ClinVar VCV000042874
in MYH7. Haplotype analysis has shown that Arg663His is not a founder mutation — it has arisen independently multiple times in ethnically diverse populations, which further supports the functional importance of arginine at position 663.

Practical Implications

For a pathogenic MYH7 variant carrier, three clinical priorities stand out. First, cardiac imaging: the absence of detectable LVH at one point in time does not rule out future development — MYH7 penetrance is age-dependent, and repeat echocardiography every 1-2 years77 repeat echocardiography every 1-2 years
ACC/AHA 2020 HCM guidelines recommend periodic clinical evaluation including ECG and echocardiogram for known sarcomere variant carriers
is standard practice. Second, AF vigilance: the Arg663His variant in particular carries an exceptionally high AF rate, and palpitations, breathlessness, or fatigue in a carrier should prompt cardiac rhythm evaluation. Third, family cascade screening: because inheritance is autosomal dominant, each first-degree relative has a 50% chance of carrying the same variant and warrants genetic testing followed by cardiologic evaluation if positive.

Interactions

The most clinically relevant gene-gene interaction involves co-occurrence of pathogenic variants in both MYH7 and MYBPC3 — a "double heterozygous" state reported in a small proportion of HCM families. Double heterozygotes (carrying pathogenic variants in both genes simultaneously) tend to present with more severe hypertrophy, earlier age of onset, and a higher rate of adverse events including sudden cardiac death compared to single-gene carriers. If a carrier of this MYH7 variant also carries a pathogenic MYBPC3 variant (e.g., rs193922376 or other confirmed pathogenic alleles), their risk profile should be escalated accordingly and discussed with a specialized HCM center.

MEFV E148Q — The Most Debated Variant in Autoinflammatory Disease

The MEFV gene encodes pyrin11 pyrin
a scaffolding protein that assembles into the pyrin inflammasome and regulates IL-1β and IL-18 release
, the two cytokines that drive the spiking fevers and serositis of familial Mediterranean fever (FMF). Most pathogenic MEFV mutations cluster in exon 10 — M694V, M680I, V726A — and produce a structurally destabilized pyrin that cannot properly gate inflammasome activation. E148Q is different: it sits in exon 2, changes a glutamic acid to the chemically similar glutamine (a conservative substitution), and is carried by as many as 25–28% of East and South Asian individuals — a population frequency that makes classic Mendelian disease causation implausible.

The Mechanism

Pyrin's B30.2 domain, encoded by exon 10, is the primary site that senses microbial toxins and interacts with the regulatory kinases RhoA-DIRAS3 and PKN1/2. The exon 10 FMF mutations bypass this gating mechanism, leading to unrestrained caspase-1 cleavage of IL-1β precursor. E148Q affects exon 2 instead, where pyrin's PYRIN domain22 PYRIN domain
the N-terminal signaling module that recruits ASC and initiates downstream caspase activation
resides. The substitution p.Glu148Gln is predicted to be tolerated by all computational missense-effect algorithms (SIFT, PolyPhen-2), consistent with the conservative nature of the glutamate-to-glutamine change. In vitro functional studies have yielded mixed results: one assay showed decreased suppression of IL-8 secretion, while a colchicine-response assay found no difference from controls. The net effect on pyrin function appears small compared to exon 10 mutations, which explains both the high population frequency and the debate over pathogenicity.

The Evidence

The debate has been running for over two decades. On the benign side: Tchernitchko et al. (2003)33 Tchernitchko et al. (2003) genotyped 233 FMF patients and 213 disease-free Sephardic Jewish relatives and found identical E148Q allele frequencies (3.62% vs 3.75%, p=0.93), leading them to conclude it "is not implicated in the development of FMF." The gnomAD v4 exome dataset contains 2,140 homozygous individuals — a scale inconsistent with a fully penetrant autosomal recessive disorder.

On the pathogenic side: Topaloglu et al. (2005)44 Topaloglu et al. (2005) examined 26 homozygous E148Q patients and found 77% had recurrent abdominal pain, 66% episodic fever, and 50% arthralgia — features indistinguishable from classic FMF. Most required colchicine therapy. A follow-up study (Topaloglu et al. 201855 Topaloglu et al. 2018) confirmed that E148Q homozygotes have milder disease than exon 10 carriers (50% mild, 46.7% moderate, 3.3% severe) with a high colchicine response rate (73.3% complete response). A 2024 pediatric cohort in Druze patients — where E148Q reaches 56% prevalence — (Awaad et al.66 Awaad et al.) found objective colchicine responsiveness (falling CRP levels) in all treated carriers, providing functional support for pathogenic relevance in that population.

The current consensus leans toward E148Q being a low-penetrance modifier rather than a classical disease allele: symptomatic expression depends on genetic background, compound heterozygosity with other MEFV alleles, and ethnicity. ClinVar (VCV000002542) reflects this tension: 31 active submissions include 1 pathogenic, 16 uncertain significance, 8 likely benign, and 6 benign. Major clinical labs (Labcorp, Invitae, Mayo Clinic) classify it as likely benign or benign.

An intriguing counterpoint: Lidar et al. (2013)77 Lidar et al. (2013) found E148Q significantly overrepresented among Ashkenazi Jewish nonagenarians (19.8% vs 2.6% in the general Ashkenazi population, p<0.0001), raising the hypothesis that chronic low-grade inflammasome priming from E148Q may enhance resistance to infections — an evolutionary trade-off analogous to sickle cell and malaria resistance.

Practical Actions

For heterozygous carriers, no treatment is needed unless there are recurring unexplained fever episodes or inflammatory attacks meeting clinical FMF criteria. For homozygous individuals, the picture is more actionable: colchicine (0.5–1 mg/day) is the standard prophylactic therapy for symptomatic FMF regardless of genotype, and most E148Q/E148Q patients with symptoms respond completely. Monitoring serum amyloid A (SAA) and CRP between attacks is recommended in FMF guidelines because sustained subclinical inflammation is the key driver of amyloidosis — the most serious long-term complication of FMF, and the one complication where E148Q has occasionally been documented in the literature.

For asymptomatic homozygotes (at least 4 of 26 in one study had no symptoms), periodic urinalysis for proteinuria (every 4–6 months) is recommended to screen for early AA amyloidosis, as per European League Against Rheumatism (EULAR) FMF guidelines.

Interactions

E148Q shows its most clinically significant behavior as a compound heterozygote with exon 10 mutations. The complex allele V726A+E148Q (both mutations on the same chromosome) produces disease more severe than V726A alone, suggesting E148Q does modify pyrin function even if it cannot cause FMF in isolation. Compound heterozygotes carrying E148Q with M694V (related SNP rs61752717) or M680I (rs61752720) are significantly less severely affected than M694V homozygotes, but still clinically symptomatic and typically require colchicine. For any E148Q carrier with FMF-like symptoms, complete MEFV panel testing including all common exon 10 variants is essential before concluding the genotype explains the phenotype.

rs3760775

FUT6

Moderate Risk Factor

FUT6 and Vitamin B12 — A Hidden Genetic Lever Common in South Asians

Vitamin B12 is essential for DNA synthesis, myelin formation, and one-carbon metabolism, yet a striking proportion of Indians — estimated at 47–70% of adults — have clinically low circulating B12 levels even when their diets are not severely restricted. Part of this disparity has a genetic explanation: a regulatory variant near the FUT6 gene that subtly suppresses fucosyltransferase11 fucosyltransferase
An enzyme that attaches fucose sugar molecules to glycoproteins and glycolipids, influencing cell-surface signalling and microbial interactions
activity and, through that, lowers circulating B12.

The Mechanism

FUT6 encodes alpha-(1,3)-fucosyltransferase 6, a Golgi enzyme that synthesises sialyl-Lewis X22 sialyl-Lewis X
A carbohydrate structure on cell surfaces that serves as a ligand for E-selectin and controls cell-cell adhesion, immune trafficking, and gut microbial colonisation
glycan structures. These fucosylated glycans line the intestinal epithelium and serve as a biological interface between host cells and gut microbiota.

rs3760775 sits in a regulatory region near the FUT6 promoter/enhancer cluster on chromosome 19p13.3. The G allele (the common allele in Europeans, the risk allele here) is associated with lower FUT6 expression, while the T allele preserves higher expression. The companion variant rs78060698 at the same locus has been shown through luciferase reporter assays and electrophoretic mobility shift assays to exhibit allele-specific promoter and enhancer activity, affecting binding of HNF4α33 HNF4α
Hepatocyte Nuclear Factor 4α — a transcription factor that regulates expression of multiple fucosyltransferases and is highly expressed in the gut and liver
, a key regulator of fucosyltransferase expression in the gut.

The functional pathway linking FUT6 to B12 most likely runs through two routes: first, altered intestinal fucosylation changes the gut microbiome composition, affecting microbial B12 synthesis and uptake; second, fucosylation of haptocorrin44 haptocorrin
Also called transcobalamin I — a glycoprotein that binds ~80% of circulating B12 but is metabolically inert; only the liver can extract B12 from it
, the dominant B12-binding protein in serum, may affect its hepatic clearance rate, altering how much B12 remains in circulation.

The Evidence

The primary evidence comes from a 2017 GWAS by Nongmaithem et al.55 Nongmaithem et al.
Nongmaithem SS et al. GWAS identifies population-specific new regulatory variants in FUT6 associated with plasma B12 concentrations in Indians. Hum Mol Genet, 2017
conducted in 1,001 healthy participants from the Pune Maternal Nutrition Study with replication in 3,418 Indians from independent cohorts (total n = 4,419). The T allele at rs3760775 was associated with higher B12 levels at genome-wide significance (meta-analysis beta = 0.25, P = 1.2×10⁻²³ on log-transformed B12). The effect allele frequency was 0.27 in Indians versus only 0.06 in Europeans (CEU), indicating the variant and its population impact are substantially higher in South Asians.

Conditional analysis showed that rs3760775 captures the primary signal at this locus, with rs3760776 and rs78060698 tagging the same or overlapping signals in the FUT6 region. The association was replicated consistently across sex, age, pregnancy status, and ethnicity subgroups.

A large PLOS Genetics meta-analysis66 PLOS Genetics meta-analysis
Grarup N et al. Genetic Architecture of Vitamin B12 and Folate Levels Uncovered Applying Deeply Sequenced Large Datasets. PLoS Genet, 2013
confirmed the FUT6/FUT3 cluster as one of 11 B12-associated loci across European and Icelandic populations. All identified loci together explain about 6.3% of B12 variance, consistent with polygenic architecture.

The mechanism was further illuminated by studies of the related FUT2 gene: Rogne et al. 201777 Rogne et al. 2017
Rogne S et al. FUT2 secretor variant p.Trp154Ter influences serum vitamin B12 concentration via holo-haptocorrin. Hum Mol Genet, 2017
demonstrated that FUT2 non-secretors have 16–22% higher total B12 but unchanged active B12 (holo-transcobalamin), because fucosylation affects haptocorrin glycosylation and hepatic clearance, not intestinal uptake per se. A similar mechanism likely operates at FUT6 — meaning the effect on total circulating B12 is real, but the clinically relevant fraction (active B12 delivered to tissues) may differ from what total serum B12 tests show.

Practical Actions

If you carry two G alleles at rs3760775 (GG), your FUT6 expression is likely somewhat lower, reducing the fucosylation-dependent mechanisms that support B12 metabolism. Given that Indians already have high background rates of B12 deficiency partly explained by vegetarian diets, this genetic variant can compound dietary inadequacy. The most direct action is to ensure adequate B12 intake and to use the most bioavailable supplemental forms when needed.

For monitoring: serum total B12 may not fully reflect cellular B12 status when fucosylation-dependent haptocorrin clearance is altered. Holo-transcobalamin (active B12) or methylmalonic acid (MMA) testing provides a more accurate picture of functional B12 sufficiency.

Interactions

FUT6 rs3760775 acts additively with FUT2 rs601338 and rs602662 — both affect fucosylation in overlapping pathways. Carrying risk alleles at multiple FUT loci compounds the effect on B12 regulation. Additionally, B12 metabolism intersects with the folate-methylation cycle: low B12 elevates homocysteine, which in turn amplifies risks from MTHFR C677T (rs1801133) heterozygosity. The FUT3-FUT5-FUT6 gene cluster on 19p13.3 shows distinct linkage disequilibrium patterns across ethnicities, which is why this variant's impact is most pronounced in South Asian populations.

NF-κB1 and Chronic Inflammation — A Locus at the Crossroads of Immune Disease

Nuclear factor kappa B (NF-κB) is the master transcription factor of the human inflammatory response. Every time your immune cells detect a pathogen, sense tissue damage, or respond to cytokine signals, NF-κB activates within minutes — switching on hundreds of target genes for cytokines, chemokines, adhesion molecules, and antimicrobial peptides. The NFKB1 gene encodes p105 and p50, the two critical NF-κB subunits that both drive and regulate this response. rs3774937 is a common intronic variant in NFKB1 that has emerged from large-scale genomic studies as a genuine susceptibility locus for chronic inflammatory disease — reaching genome-wide significance across ulcerative colitis and a pleiotropic analysis of five distinct immune-mediated conditions.

The Mechanism

rs3774937 lies within an intron of NFKB1 at chromosome 4 position 102,513,096 (GRCh38) and does not alter the p105 or p50 protein sequence. Intronic variants can nonetheless influence gene function through several mechanisms: altering splicing enhancer or silencer sequences11 splicing enhancer or silencer sequences
Intronic regulatory elements that guide which exons are included in the final mRNA transcript
, modifying transcription factor binding sites embedded within introns, or marking linkage disequilibrium22 linkage disequilibrium
Non-random co-inheritance of alleles at nearby loci; rs3774937 may tag a causal functional variant at the NFKB1 locus without being directly causal itself
with a functional variant that has not yet been fully characterized. The nearby rs28362491 variant — a four-nucleotide ATTG insertion-deletion in the NFKB1 5' regulatory region with documented effects on promoter activity — is the most likely functional candidate at this chromosomal location; rs3774937 has been included in NFKB1 haplotype analyses alongside rs28362491, rs230521, rs230510, and rs4648068, suggesting they tag overlapping regulatory variation.

The biological consequence of reduced NF-κB1 p50 activity is a shift in the immune inflammatory set point. The p50 homodimer acts as a transcriptional repressor of pro-inflammatory genes — it competes with the activating p50/p65 heterodimer for DNA binding at NF-κB response elements. Reduced p50 protein means less of this negative feedback, resulting in higher baseline cytokine production and exaggerated inflammatory responses to immune triggers. This mechanism connects reduced NFKB1 transcriptional activity to the observed associations with conditions driven by chronic or dysregulated inflammation.

The Evidence

The strongest evidence comes from genome-wide association studies of inflammatory bowel disease. Liu et al. (2015)33 Liu et al. (2015)
Association analyses identify 38 susceptibility loci for inflammatory bowel disease and highlight shared genetic risk across populations. Nature Genetics.
— a trans-ancestry meta-analysis — identified rs3774937 at the NFKB1 locus with genome-wide significance for ulcerative colitis (OR approximately 1.10, p=5×10⁻¹⁴, risk allele frequency 0.33 in the discovery cohort). The modest odds ratio paired with extreme statistical significance across tens of thousands of samples confirms this is a robust genetic association.

Ellinghaus et al. (2016)44 Ellinghaus et al. (2016)
Analysis of five chronic inflammatory diseases identifies 27 new associations and highlights disease-specific patterns at shared loci. Nature Genetics.
conducted a pleiotropy analysis in over 86,000 European-ancestry individuals across ankylosing spondylitis, Crohn's disease, psoriasis, primary sclerosing cholangitis, and ulcerative colitis. The NFKB1 locus (tagged by rs3774937) reached p=2×10⁻¹⁸ in the cross-disease analysis — among the most significant signals in the entire study. This pleiotropic signal indicates that the same regulatory variation at NFKB1 contributes to multiple distinct immune-mediated conditions, consistent with NF-κB1's role as a central node in immune homeostasis rather than a disease-specific factor.

Beyond IBD and autoimmune phenotypes, Cabrera et al. (2014)55 Cabrera et al. (2014)
Intronic variants in the NFKB1 gene may influence hearing forecast in patients with unilateral sensorineural hearing loss in Menière's disease. PLoS One.
reported that rs3774937 C allele carriers with unilateral Menière's disease reach hearing stage 3 (greater than 40 dB loss) significantly faster than non-carriers (p=0.009, log-rank corrected). The study enrolled 716 Menière's disease cases and 1,628 controls — the association was specific to unilateral disease, consistent with an inflammatory component in the inner ear lesion driven by NF-κB1 pathway variation.

In the transplant setting, Kuba et al. (2020)66 Kuba et al. (2020)
NFKB1 gene single-nucleotide polymorphisms: implications for graft-versus-host disease in allogeneic hematopoietic stem cell transplantation. Annals of Hematology.
reported strong association of rs3774937 genotype with acute graft-versus-host disease in 109 transplant recipients — evidence that the variant influences post-transplant immune responses, where NF-κB signaling governs allogeneic T cell activation and cytokine production.

Practical Actions

The per-allele OR of approximately 1.10 for ulcerative colitis and similarly modest effect sizes for other conditions mean that rs3774937 does not determine disease fate — it tilts the immune landscape modestly in the direction of dysregulated inflammatory signaling. For CC homozygotes (~9% of Europeans), two copies compound this tilt. The actionable consequence is proactive monitoring of inflammatory burden and support for the NF-κB regulatory mechanisms this locus influences.

Monitoring circulating biomarkers of systemic inflammation — particularly high-sensitivity CRP (hsCRP) and calprotectin if gastrointestinal symptoms arise — provides direct insight into whether this genetic signal is translating into detectable inflammatory activity in the individual. Omega-3 fatty acids at therapeutic doses reduce NF-κB activation and downstream cytokine production through resolvin and protectin pathways that directly intersect with the NF-κB signaling axis.

Interactions

rs3774937 sits within the same NFKB1 genomic region as several other variants studied in this encyclopedia: rs28362491 (the -94ATTG functional indel in the NFKB1 5' regulatory region, associated with cardiovascular and autoimmune risk), rs230523 (intronic NFKB1 variant associated with infection susceptibility), and rs4648127 (a rare protective intronic variant associated with reduced lung cancer risk). These variants are in partial linkage disequilibrium and may partly tag the same causal signal — users who carry risk alleles at multiple NFKB1 intronic positions should treat their combined signal with caution to avoid double-counting the same underlying regulatory effect.

TLR4 (rs4986790) and TLR1 (rs5743708) variants acting upstream of NF-κB activation may compound with NFKB1 variants to produce additive reductions in pathogen-driven immune gene expression when both receptor-level and transcription-factor-level variation co-occur.

IKBKE — The Dual-Pathway Kinase at the Heart of Psoriasis Susceptibility

Psoriasis is not a single disease driven by a single pathway — it is the collision of two inflammatory cascades: the NF-kB signaling axis11 NF-kB signaling axis
NF-kB (nuclear factor kappa-light-chain-enhancer of activated B cells) is a master transcription factor that controls genes for TNF-α, IL-1β, IL-6, and many other pro-inflammatory cytokines; it is the primary driver of the keratinocyte hyperproliferation and immune cell recruitment seen in psoriatic plaques
and the type I interferon axis22 type I interferon axis
Type I interferons (IFN-α and IFN-β) are signaling proteins produced by plasmacytoid dendritic cells and keratinocytes that amplify innate immune responses; they are particularly prominent in early psoriasis lesions and plaque initiation
. The gene IKBKE encodes IKK-epsilon (IKKε)33 IKK-epsilon (IKKε)
Inhibitor of NF-kB kinase subunit epsilon — a serine/threonine kinase related to canonical IKKα and IKKβ but with distinct substrates and a unique position at the crossroads of innate antiviral and inflammatory signaling
, a kinase that sits at the intersection of both pathways and governs the amplitude of both responses simultaneously. The rs41298997 intronic variant in IKBKE was identified in the largest psoriasis GWAS to date as a novel, genome-wide-significant susceptibility locus, with the T allele conferring a modest but statistically robust increase in psoriasis risk.

The Mechanism

IKKε operates at a hub where three signaling pathways converge. First, it phosphorylates IRF3 and IRF744 IRF3 and IRF7
Interferon regulatory factors 3 and 7 — transcription factors that, once phosphorylated by IKKε or its homolog TBK1, undergo conformational change, dimerize, and translocate to the nucleus to activate IFN-β and IFN-α gene transcription; this is the primary amplifier of the antiviral and early innate immune response
, driving type I interferon production. Second, IKKε phosphorylates inhibitors of NF-kB (IkBs), freeing the NF-kB complex to translocate to the nucleus and switch on pro-inflammatory cytokine genes including TNF-α, IL-1β, and IL-6. Third, IKKε phosphorylates STAT1 at serine-70855 STAT1 at serine-708
This phosphorylation selectively promotes assembly of the ISGF3 complex (STAT1:STAT2:IRF9) over GAF (STAT1 homodimer) formation; ISGF3 drives interferon-stimulated gene expression while GAF drives a different transcriptional program; IKKε thus acts as a molecular switch that shapes which arm of the interferon response dominates
, tuning the balance between type I and type II interferon-driven transcriptional programs.

The rs41298997 variant lies within an intron of IKBKE and is presumed to influence gene expression or splicing rather than protein structure, although the exact regulatory mechanism has not been characterized at the molecular level. Its location at 1q32.1 — a chromosomal region enriched for immune regulatory elements — and the genome-wide significance of its association with psoriasis places it in the same functional class as other regulatory psoriasis loci that modulate pathway-level signaling intensity rather than abolishing it.

The Evidence

The primary evidence comes from the Tsoi et al. 2017 meta-analysis66 Tsoi et al. 2017 meta-analysis
Large scale meta-analysis characterizes genetic architecture for common psoriasis associated variants. Nature Communications, 2017. Combined effective sample size >39,000 individuals across eight Caucasian cohorts.
, which identified 16 novel genome-wide-significant psoriasis loci. The rs41298997-T allele reached genome-wide significance (OR=1.13, p=2.37×10⁻⁸) with a risk allele frequency of approximately 19% in cases versus 18% in controls — a small but consistent enrichment across tens of thousands of individuals. The study performed pathway analysis showing that genes at novel loci, including IKBKE, clustered in the "Regulation of I-kB kinase/NF-kB cascade" pathway (11 associated loci) and in T-cell regulatory elements, confirming the biological relevance of the association.

The molecular role of IKKε was established by two mechanistic studies. Fitzgerald et al. 200377 Fitzgerald et al. 2003
Fitzgerald KA et al., IKKε and TBK1 are essential components of the IRF3 signaling pathway. Nature Immunology 2003.
demonstrated that IKKε and TBK1 together are required for IRF3-mediated IFN-β induction in response to viral signals, placing IKKε as an essential gatekeeper of innate antiviral immunity. Ng et al. 201188 Ng et al. 2011
Ng SL et al., IκB kinase ε (IKKε) regulates the balance between type I and type II interferon responses. PNAS 2011.
then showed that IKKε does not merely activate interferons but actively shapes which interferon-driven transcriptome is expressed by determining whether ISGF3 or GAF forms — a distinction with profound implications for whether cells respond in an antiviral versus inflammatory mode.

The OR of 1.13 per T allele is modest by Mendelian disease standards but is consistent with the polygenic architecture of psoriasis, where dozens of loci each contribute small effects that combine to produce the 60–80% heritability of the disease. The T allele frequency is notably higher in East Asian populations (~61%) compared to Europeans (~24%) and Africans (~18%), which may contribute to differences in psoriasis prevalence and phenotype across ancestries.

Practical Implications

For CT heterozygotes and TT homozygotes, the most actionable implication is awareness: psoriasis tends to develop in genetically predisposed individuals following environmental triggers — streptococcal throat infections, mechanical skin trauma (Koebner phenomenon), certain medications (lithium, beta-blockers, antimalarials, rapid steroid withdrawal), and psychological stress. Recognizing early plaque formation and seeking dermatological evaluation promptly improves outcomes significantly.

The IKBKE locus is of particular interest for therapy: IKKε inhibition has been investigated as a therapeutic strategy in inflammatory disease, and IKKε's dual role in NF-kB and interferon signaling means that dysregulation at this node may influence both the initiating events (interferon-driven plasmacytoid dendritic cell activation) and the sustaining events (NF-kB-driven keratinocyte hyperproliferation and cytokine amplification) in psoriasis pathogenesis. While no IKKε-targeted therapy is currently approved for psoriasis, existing biologics that block downstream cytokines (anti-IL-17, anti-IL-23, anti-TNF) all target effectors downstream of the NF-kB pathway that IKBKE feeds into.

Interactions

IKBKE functionally interacts with STAT4 (rs7574865) in psoriasis susceptibility: STAT4 is activated by IL-12 and IL-23 downstream of the NF-kB pathway that IKKε activates. Individuals carrying risk alleles at both loci may have amplified Th1/Th17 polarization, as IKKε drives the upstream cytokine milieu (via IFN and NF-kB signaling) while STAT4 determines how sensitively T cells respond to those cytokines. This interaction has not been formally tested in psoriasis but parallels the well-documented IRF5–STAT4 interaction in lupus susceptibility.

The IKBKE-encoded kinase also functions within a complex containing TBK1 and scaffold proteins including TANK, NAP1, and SINTBAD. Variants affecting other members of this kinase complex (TBK1 variants) may compound with rs41298997 to amplify innate immune signaling, though no interaction data specific to psoriasis yet exist.

rs476828

MC4R MC4R-region BMI variant

Strong Risk Factor

The MC4R Haplotype — A Third Window into Appetite Regulation

The melanocortin-4 receptor (MC4R) gene sits at the center of your brain's appetite control system. Located in the hypothalamus11 hypothalamus
the brain region governing hunger, satiety, and energy balance
, MC4R receives signals from leptin and melanocortin hormones to generate "stop eating" commands. The stretch of chromosome 18 roughly 188 kilobases downstream of MC4R contains a regulatory haplotype block — a cluster of variants inherited together — that controls how much MC4R your neurons produce. rs476828 sits within this block.

This SNP tags the same obesity-associated haplotype as the better-known rs17782313, the second strongest common obesity genetic signal in humans22 second strongest common obesity genetic signal in humans
after FTO rs9939609
. In Europeans, rs476828 and rs17782313 are in perfect linkage disequilibrium33 perfect linkage disequilibrium
r²=1 in the CEU HapMap sample, meaning the two variants are essentially interchangeable as markers for the same underlying haplotype
. In African-ancestry populations the LD is weaker (r²=0.526), meaning rs476828 carries some independent information there. Together with rs12970134, these three variants form the risk haplotype C–C–A, which carries OR=1.796 for obesity.

The Mechanism

rs476828 is an intergenic variant in a regulatory element44 regulatory element
a stretch of DNA that controls gene transcription without encoding protein itself
that modulates MC4R expression in hypothalamic neurons. The full signaling chain works as follows: fat cells secrete leptin55 leptin
a hormone that reports energy stores to the brain
, which activates POMC neurons66 POMC neurons
proopiomelanocortin neurons in the arcuate nucleus that generate appetite-suppressing signals
, which release alpha-melanocyte stimulating hormone (α-MSH), which binds MC4R. When MC4R fires, it suppresses appetite and increases energy expenditure. When regulatory variants reduce MC4R expression, this entire cascade is weakened — fewer receptors means quieter "stop eating" signals and reduced metabolic drive.

The C allele of rs476828 tags a haplotype associated with reduced MC4R promoter activity. Epigenetic studies77 Epigenetic studies
MeQTL analyses examining DNA methylation quantitative trait loci
show that this regulatory block is associated with increased MC4R promoter methylation, further suppressing expression. The net result is a hypothalamic satiety circuit that requires stronger signals to fire — meaning larger portions before fullness registers, and a higher baseline appetite drive.

The Evidence

The original landmark GWAS88 landmark GWAS
Loos et al., Nature Genetics 2008, n=16,876 discovery + 60,352 replication
identified this MC4R-region signal as the second strongest common BMI locus in humans. Per C allele, adults showed a 0.05 Z-score BMI increase (p=2.8×10⁻¹⁵); children aged 7–11 showed a larger 0.13 Z-score effect (p=1.5×10⁻⁸), and severe childhood obesity odds reached OR=1.30 (p=8.0×10⁻¹¹). rs476828 was directly studied in a childhood obesity cohort99 childhood obesity cohort
Grant et al., 728 obese European American children and 3,960 controls
and yielded OR=1.145 (p=0.042), consistent with its role as a perfect proxy for rs17782313 in European-ancestry individuals.

A 2020 haplotype study1010 2020 haplotype study
Wei et al., Mol Med, 1,836 Chinese participants
examined rs476828 directly and found the C allele significantly elevated in obesity cases versus controls (17.1% vs 10.9%, p<0.001), with a dominant-model OR of 1.585 (95% CI=1.176–2.136). The full C–C–A haplotype across rs17782313, rs476828, and rs12970134 yielded OR=1.796 (95% CI=1.447–2.229), illustrating that these three markers collectively tag a single high-risk regulatory configuration.

A prospective cohort study in 5,724 women1111 prospective cohort study in 5,724 women
Qi et al., Hum Mol Genet 2008
found that this MC4R locus (tagged by rs17782313) was associated with higher total energy and dietary fat intake, greater 10-year BMI gain, and 14% increased type 2 diabetes risk per C allele after BMI adjustment — indicating effects on metabolic regulation beyond adiposity alone.

A sex-specific layer was uncovered by neuroimaging in 284 adults1212 neuroimaging in 284 adults
Horstmann et al., PLoS One 2013
: female homozygous CC carriers showed increased gray matter volume in the right amygdala, hippocampus, and orbitofrontal cortex — regions encoding emotional memory and food reward — as well as significantly elevated disinhibition and emotional eating scores. These effects were absent in men, suggesting the MC4R regulatory haplotype has sex-specific impacts on the neural circuits governing food motivation.

Practical Implications

Carrying the C allele at rs476828 means your hypothalamic MC4R system operates with reduced receptor density. Satiety signals that would normally be sufficient to stop eating are muted — your brain requires more food intake before the "full" signal reaches threshold. This creates a persistent biological headwind for weight management, but one that responds well to strategies that compensate for weakened internal signals with external structure.

The fat and energy intake data from the women's cohort are particularly actionable: the genetic effect operates partly through appetite drive for energy-dense foods. Interventions that reduce the palatability and caloric density of available food (meal preparation at home, portion pre-commitment, structured eating environments) directly address the mechanism.

The sex-specific neuroimaging data suggest that women with CC genotype may benefit especially from interventions targeting emotional and reward-driven eating, since this is where the documented brain structural differences concentrate.

Interactions

rs17782313 and rs12970134: rs476828 sits in the same regulatory block as these two MC4R-region variants. In Europeans, rs476828 and rs17782313 are in perfect LD (r²=1) — they are interchangeable markers for the same haplotype. The three-SNP haplotype C–C–A (rs17782313–rs476828–rs12970134) carries the highest obesity risk (OR=1.796), while individual SNP effects are smaller. Carrying risk alleles at all three does not simply stack independent risks; they largely reflect the same underlying regulatory signal.

FTO rs9939609: The combined effect1313 combined effect
documented across multiple populations
of MC4R and FTO risk alleles exceeds either alone — MC4R acts through appetite suppression while FTO acts through thermogenesis, so the two pathways are partially independent and their effects add. Combined risk genotypes in one pediatric study conferred 2.45-fold increased obesity odds. Addressing both pathways simultaneously — appetite structure for MC4R, thermogenic activity for FTO — provides complementary benefit.

TMEM150B 19q13.42 — The Autophagy Modulator and Your Ovarian Clock

On chromosome 19 at band q13.42, within a dense cluster of reproductive-aging loci, lies a variant that has drawn attention for its unexpectedly large effect on early menopause risk. The rs4806660 SNP falls within the gene TMEM150B — also known as DRAM3 (Damage-Regulated Autophagy Modulator 3) — a transmembrane protein expressed in oocytes whose precise role in reproductive aging is still being worked out.

A note on gene attribution: early GWAS publications initially assigned this locus to the nearby BRSK1 gene (also at 19q13.42) based on physical proximity. Subsequent fine-mapping established that rs4806660 and the lead SNP rs1172822 are in near-complete linkage disequilibrium (R²=0.965) and both fall within the TMEM150B gene body. Some older references still label this locus as BRSK1 or "ZNF346" (a chromosome 5 gene unrelated to this locus); the current consensus assigns it to TMEM150B.

The Mechanism

TMEM150B encodes a five-transmembrane-domain protein that belongs to the DRAM family11 DRAM family
Damage-Regulated Autophagy Modulators — a class of lysosomal and endosomal membrane proteins that regulate autophagic flux under conditions of cellular stress
. The protein localises to lysosomes, endosomes, and the plasma membrane. Overexpression studies show that TMEM150B promotes autophagosome accumulation and enhances autophagic flux under baseline conditions; it also promotes cell survival during glucose deprivation through an autophagy-independent mechanism.

Autophagy22 Autophagy
the cellular self-digestion pathway that recycles damaged organelles and proteins to maintain oocyte quality
is increasingly recognised as a critical maintenance system in the dormant oocyte pool. Human oocytes are arrested for decades; their long-term viability depends on continuous quality-control mechanisms to handle accumulated protein aggregates and dysfunctional mitochondria. Defective autophagy in granulosa cells and oocytes has been linked to accelerated follicular atresia and POI in animal models.

rs4806660 is an intronic variant, meaning it does not alter the protein sequence. The most likely mechanism is cis-regulatory33 cis-regulatory
affecting the expression level of a nearby gene rather than its protein sequence
: the C allele may alter TMEM150B expression in ovarian cells in ways that subtly compromise the autophagy-mediated quality control sustaining the follicle pool. However, the precise regulatory effect of this specific intronic change on TMEM150B expression in human ovarian tissue has not yet been directly characterised.

It is important to note that when Tmem150b was completely knocked out in mice, female fertility was normal and hormone levels were unchanged44 female fertility was normal and hormone levels were unchanged
Liu et al. 2020, Scientific Reports
. This may indicate functional redundancy within the DRAM family, species differences in the specific gene's requirement, or that the effect of the common human intronic variant operates through a mechanism distinct from complete gene ablation.

The Evidence

The 19q13.42 locus was first identified by Stolk et al. 200955 Stolk et al. 2009
Loci at chromosomes 13, 19 and 20 influence age at natural menopause. Nature Genetics, 2009
in a two-stage GWAS of 2,979 European women (Rotterdam Study and TwinsUK), with the lead SNP rs1172822 reaching p = 6.3 × 10⁻¹¹. Subsequent fine-mapping using imputed data identified rs4806660 as an additional SNP in the same LD block with even stronger statistical support.

The most informative study for clinical interpretation is the Breakthrough Generations Study66 Breakthrough Generations Study
Murray et al. 2011, Human Molecular Genetics, ~2,000 women with early menopause before age 46 versus controls
. In this case-control design, each C allele was associated with an odds ratio of 1.45 (95% CI 1.32–1.59, p = 8.88 × 10⁻¹⁶) for experiencing menopause before age 46. Notably, the observed odds ratio was substantially larger than the 1.20 expected from quantitative trait estimates — suggesting a non-linear effect where the C allele may have a disproportionate impact on the earliest end of the menopause-age distribution, affecting those most biologically susceptible to early follicular depletion.

Beyond menopause timing, a Brazilian IVF cohort study by Setti et al. 201277 Setti et al. 2012
A chromosome 19 locus positively influences the number of retrieved oocytes during stimulated cycles in Brazilian women. Journal of Assisted Reproduction and Genetics, 2012
found that women carrying the T allele (the protective/common allele) had significantly more antral follicles (+2.54 per cycle, P=0.041) and more retrieved oocytes (+1.41 per cycle, P=0.041) during controlled ovarian stimulation — providing a direct, clinically measurable link between this variant and the functional ovarian reserve available for IVF.

A Mashhad cohort study 202188 Mashhad cohort study 2021
Genetic Determinants of Premature Menopause in a Mashhad Population Cohort. IJFS, 2021
found that the C allele in Iranian women was associated with premature menopause risk (OR 3.09, 95% CI 1.17–8.16 in a recessive model), though this did not survive Bonferroni correction. The direction of effect is consistent with European data, suggesting the association may generalise across Middle Eastern populations, but effect sizes remain uncertain outside the original European GWAS context.

Population differences are striking: the C allele is common in Europeans (~36%), African-ancestry (~37%), and South Asian (~41%) populations, but is found at substantially lower frequency in East Asian populations (~8%). This means the variant's contribution to population-level menopause timing variation is much smaller in East Asian women.

Practical Actions

For women with the common TT genotype, this variant suggests a somewhat lower genetic load for early follicular depletion at this locus. AMH testing remains the most direct measure of current reserve and is appropriate for fertility planning regardless of genotype.

For C allele carriers, the primary implication is a modestly elevated probability of earlier reproductive aging — a probabilistic shift, not a certainty. The variant may be particularly relevant when ovarian response to gonadotropin stimulation is lower than expected, since Setti et al. suggest this locus influences the gonadotropin-responsive follicle pool.

Interactions

TMEM150B rs4806660 + MCM8 rs16991615 (19q13.42 + 20p12.3 menopause loci): MCM8 rs16991615 (E341K, DNA repair helicase) is among the most robustly replicated menopause-timing loci with an effect of approximately 1 year per A allele. Women carrying the common risk GG genotype at MCM8 (absent protective allele) who also carry one or two C alleles at rs4806660 accumulate two independent additive hits on follicular reserve through distinct biological pathways — DNA repair and autophagy regulation. A compound action for this combination emphasising early AMH baseline testing and proactive reproductive timeline discussion is warranted.

TMEM150B rs4806660 + PRRC2A rs1046089 (immune-linked menopause locus): rs1046089 in the HLA-region PRRC2A gene acts through immune modulation of follicular atresia. Combined C allele burden at rs4806660 alongside the A allele at rs1046089 converges on reproductive aging through two further independent mechanisms (autophagy dysregulation and immune-mediated follicle depletion). Women with risk alleles at both loci may benefit from earlier fertility assessment than either variant alone would indicate.