rs804279

GATA4 GATA4/NEIL2 PCOS Susceptibility

Moderate Risk Factor

GATA4/NEIL2 — A Chromosomal Neighbourhood Tied to PCOS Androgen Excess

On the short arm of chromosome 8, at position 8p23.1, a compact genomic neighbourhood hosts two functionally distinct genes: GATA411 GATA4
GATA-binding protein 4, a zinc-finger transcription factor essential for cardiac and gonadal development; expressed in ovarian granulosa and theca cells where it regulates steroidogenic gene expression
and NEIL222 NEIL2
nei-like DNA glycosylase 2, a base excision repair enzyme that removes oxidised purines from DNA
. Between and around them sits rs804279, an intergenic variant that large-scale genome-wide association studies have linked to polycystic ovary syndrome (PCOS) susceptibility in women of European ancestry.

The Mechanism

The rs804279 variant lies approximately 6,400 bp upstream of GATA4 and 3,300 bp from NEIL2 (within the NEIL2 RefSeqGene interval, NG_053043.1:g.1718). As an intergenic regulatory variant, it does not change a protein sequence; instead, it likely modulates the expression or chromatin accessibility of one or both flanking genes in hormone-responsive tissues. GATA4 is a transcription factor expressed in ovarian theca cells, where it drives expression of steroidogenic enzymes including CYP17A1, the key enzyme in androgen biosynthesis. The locus also harbours FDFT133 FDFT1
farnesyl-diphosphate farnesyltransferase 1, the first committed enzyme in cholesterol biosynthesis that provides the substrate for all steroid hormones
, which sits approximately 29 kb downstream — giving this chromosomal neighbourhood three converging mechanisms for steroid hormone biology.

The T allele is the globally common allele (~72% frequency) but is the risk allele for PCOS at this locus. This contrasts with the typical GWAS pattern where the minor allele carries risk, and implies that the protective A allele (reference, ~28% globally) represents a recent or population-specific buffering variant — or that regulatory activity at this locus is required for normal PCOS-associated androgen elevation, with the A allele conferring partial protection.

The Evidence

The key PCOS association at rs804279 was established in a large European meta-analysis44 large European meta-analysis
Day FR et al. Large-scale genome-wide meta-analysis of polycystic ovary syndrome suggests shared genetic architecture for different diagnosis criteria. PLoS Genet 14:e1007813, 2018
combining 10,074 PCOS cases and 103,164 controls. Among 14 genome-wide significant PCOS loci, rs804279 was notable for one specific reason: it was the only locus showing significant heterogeneity across diagnostic criteria (heterogeneity p=2.6×10⁻⁵). The effect was largest in women diagnosed by NIH criteria (hyperandrogenism plus oligo-anovulation) and smallest in self-reported cases — consistent with a variant that specifically tags the androgen-excess component of PCOS rather than the full syndrome as broadly defined.

Computational functional analysis has categorised rs804279 among variants with predicted functional relevance55 variants with predicted functional relevance
Prabhu BN et al. Conceptualization of functional SNPs of PCOS genes: an in silico approach. J Endocrinol Invest 44:1783–1793, 2021
in the GATA4/NEIL2 region, supporting its role as a regulatory rather than neutral variant.

The evidence level is moderate: the locus has genome-wide significance in the Day 2018 meta-analysis, but functional characterisation of what rs804279 specifically does to GATA4 or NEIL2 expression remains incomplete. The diagnostic-criteria heterogeneity means the association is most reliable when PCOS is defined by the NIH standard (requiring both hyperandrogenism and cycle abnormality) rather than by self-report or broader Rotterdam criteria alone.

Practical Actions

Because this locus appears most strongly associated with the hyperandrogenic PCOS subtype, monitoring should focus on androgen biomarkers: free testosterone, DHEAS, and androstenedione are the most directly relevant, with AMH and LH:FSH ratio providing additional characterisation of the reproductive phenotype.

Inositol supplementation (myo-inositol) has the strongest evidence for improving ovarian responsiveness in PCOS with FSH resistance, acting downstream of the gonadotropin signalling cascade that GATA4 expression influences. Women with the TT genotype who have clinical signs (irregular cycles, hirsutism, acne, difficulty conceiving) have clear grounds for proactive evaluation.

Interactions

The GATA4/NEIL2 locus acts independently of the DENND1A (rs7852296) and LHCGR (rs13405728) PCOS susceptibility loci established in the same European meta-analyses. Women carrying risk alleles at multiple PCOS loci accumulate additive PCOS susceptibility. The GATA4 locus appears to operate primarily through the hyperandrogenic axis (linked to androstenedione and testosterone excess) while DENND1A operates through FSH receptor trafficking — distinct but complementary mechanisms both converging on follicular arrest and androgen excess. A compound interaction between homozygous TT at rs804279 and the DENND1A risk haplotype (rs7852296-AA or homozygous risk) would represent convergent failure of both steroidogenic regulation and FSH receptor recycling; this combination warrants investigation as a compound action (no published combined effect size yet).

HLA-DRA rs9268839 — The Class II Sentinel for Rheumatoid Arthritis

The region just upstream of HLA-DRA on chromosome 6p21.32 is home to one of the most powerful common genetic risk signals for rheumatoid arthritis (RA) in the human genome. rs9268839 is an intergenic variant located approximately 16 kilobases upstream of HLA-DRA11 HLA-DRA
HLA class II histocompatibility antigen, DR alpha chain — the gene encoding the alpha subunit of the HLA-DR heterodimer on antigen-presenting cells
, within a dense cluster of HLA class II genes including HLA-DRB5, HLA-DRB9, and TSBP1. The G allele carries an odds ratio of approximately 2.47 per copy for RA in European populations — among the largest effect sizes of any common non-coding variant discovered in autoimmune disease genetics.

The Mechanism

HLA-DRA encodes the invariant alpha chain of the HLA-DR molecule, a class II major histocompatibility complex (MHC) protein expressed on the surface of antigen-presenting cells — dendritic cells, macrophages, and B cells. HLA-DR is a heterodimer: the alpha chain (encoded by HLA-DRA) pairs with a beta chain (encoded by HLA-DRB1, DRB3, DRB4, or DRB5 depending on the haplotype). The combined molecule presents peptide fragments from pathogens and self-proteins to CD4+ helper T cells, initiating adaptive immune responses. The critical insight is that HLA-DRA is largely non-polymorphic22 HLA-DRA is largely non-polymorphic
Unlike HLA-DRB1, which has hundreds of functionally distinct alleles, HLA-DRA varies very little between individuals — its alpha chain is nearly identical across all people
. rs9268839 is therefore not a coding variant in HLA-DRA itself, but a tag SNP marking an extended haplotype block that includes the HLA-DRB5 locus and influences local regulatory architecture through linkage disequilibrium with functional variants in flanking HLA-DRB genes.

The precise functional variant this SNP tags has not been fully resolved — a common situation in the HLA region, where extreme polymorphism and LD make fine-mapping challenging. What is clear is that the G-allele haplotype is enriched for HLA-DR conformations that present citrullinated self-peptides particularly efficiently to autoreactive CD4+ T cells, driving the production of anti-citrullinated protein antibodies (ACPA)33 anti-citrullinated protein antibodies (ACPA)
The hallmark autoantibody in seropositive RA; detected by the anti-CCP test and predictive of joint erosion
— the defining pathological antibodies in seropositive RA.

The Evidence

Rheumatoid arthritis. The GWAS signal at rs9268839 is among the most replicated findings in autoimmune genetics. In the landmark Okada et al. 2014 meta-analysis44 Okada et al. 2014 meta-analysis
Genetics of rheumatoid arthritis contributes to biology and drug discovery. Nature 2014; 29,880 RA cases and 73,758 controls
spanning European and Asian populations, the HLA class II locus at chromosome 6p21.32 (of which rs9268839 is a primary tag) showed the strongest association in the entire genome. In the Laufer et al. 2019 trans-ethnic fine-mapping study55 Laufer et al. 2019 trans-ethnic fine-mapping study
Genetic influences on susceptibility to rheumatoid arthritis in African-Americans. Hum Mol Genet 2019; 916 AA cases + >100,000 European/East Asian participants
, the HLA-DRB1/DRA region reached p<1×10⁻²⁵⁰ in Europeans — the maximum reportable significance threshold. The rs9268839-G allele carries an OR of approximately 2.47 (95% CI 2.39–2.55) in European cohorts and 1.90 (95% CI 1.81–1.99) in East Asian cohorts, both at p=1×10⁻²⁵⁰.

Sarcoidosis. The G allele at rs9268839 has also been implicated in sarcoidosis susceptibility, a granulomatous inflammatory disease driven by dysregulated antigen presentation in the lung and lymph nodes. In Levin et al. 201366 Levin et al. 2013
Association of ANXA11 genetic variation with sarcoidosis in African Americans and European Americans. Genes Immun 2013; 1,689 cases and 1,252 controls
, rs9268839 near HLA-DRA was identified as a sarcoidosis risk SNP showing significant SNP-SNP interaction with ANXA11 rs1049550, suggesting that the HLA-DRA region risk variant amplifies disease susceptibility when combined with ANXA11 variants that impair calcium-dependent immune regulation.

Practical Actions

Carrying one or two G alleles does not mean RA or sarcoidosis will develop — these are probabilistic risk elevations, not deterministic outcomes. However, G carriers should be alert to early symptoms and pursue earlier investigation when musculoskeletal symptoms appear. The most actionable step for G carriers is anti-CCP antibody testing: ACPA can appear years before clinical RA and predicts severity and progression. G carriers who are also ACPA-positive should seek rheumatology review proactively.

Smoking is the most important environmental co-trigger: the HLA class II region and cigarette smoke interact to produce citrullinated lung antigens that drive ACPA production, multiplying RA risk superadditively in G-allele carriers.

Interactions

rs9268839 operates in the same antigen-presentation pathway as rs660895 (HLA-DRB1 shared epitope tag), rs6910071 (TSBP1/C6orf10 intronic variant), and rs2476601 (PTPN22 R620W). These variants act at distinct steps — peptide groove composition, MHC region haplotype structure, T-cell receptor signalling — and their co-occurrence is associated with multiplicative rather than additive RA risk. rs9268839 and rs660895 are in moderate LD in Europeans and may partly tag the same underlying haplotype, but they are not fully redundant: rs9268839 explains additional variation in RA susceptibility beyond the shared epitope alone.

IL1RL1 (ST2) — A Protective Variant in the IL-33 Alarm Circuit

The IL1RL1 gene11 IL1RL1 gene
Interleukin-1 Receptor-Like 1, encoding the ST2 protein — the primary cell-surface receptor for the alarmin cytokine IL-33, expressed on mast cells, eosinophils, ILC2s, and Th2 lymphocytes
is one of the most replicated loci in allergy and asthma genetics. Most research focuses on variants that increase risk. rs13424006 is notable for the opposite reason: the C allele at this intronic position is protective — carriers have a measurably lower risk of developing late-onset wheeze, the wheezing phenotype most strongly linked to persistent adult asthma.

The IL-33/ST2 axis functions as an epithelial damage alarm. When airway cells are injured by viruses, allergens, or pollutants, they release IL-3322 IL-33
A nuclear alarmin cytokine that, when released from damaged epithelial cells, binds to ST2 on ILC2s, mast cells, and eosinophils, triggering a rapid Th2 cascade
. The intensity of that alarm depends on the ratio between membrane-bound ST2 (which transmits the signal) and soluble sST2 (which mops up free IL-33 before it can reach its target cells). IL1RL1 variants modulate this ratio — some increasing signalling capacity, others buffering it. rs13424006 sits in the intron where this regulatory balance is calibrated.

The Mechanism

rs13424006 maps to position 102,350,776 on chromosome 2 (GRCh38), within the 10th intron of IL1RL1. Like rs10208293 — the other intronic IL1RL1 variant that specifically tags late-onset wheeze in the same study — rs13424006 is presumed to act through an eQTL mechanism, influencing the relative production of membrane-bound IL1RL1 and the soluble decoy sST2 without altering the protein sequence. The C allele tags a haplotype associated with more effective IL-33 buffering: greater production or activity of sST2, which intercepts IL-33 before it activates eosinophils and ILC2s.

rs13424006 and rs10208293 both sit in the 10th intron and both specifically associate with the late-onset wheeze phenotype, but they are genotyped separately and may tag distinct regulatory elements within the same intronic region. Whether their effects are additive, redundant, or mechanistically coupled is not yet established.

The Evidence

The protective effect of the C allele was first defined in a large two-cohort study by Savenije et al. (2014)33 Savenije et al. (2014), which meta-analysed data from 2,007 children in the Dutch PIAMA cohort and 7,247 children in the UK ALSPAC cohort, following wheeze phenotypes from birth to age 8. The analysis distinguished four phenotypes: no wheeze, early-transient wheeze, intermediate-onset, and late-onset wheeze. rs13424006 was one of only two IL1RL1 SNPs that associated specifically with late-onset wheeze (OR approximately 0.74, 95% CI 0.63–0.87), meaning C allele carriers develop this phenotype at substantially lower rates than TT homozygotes.

The late-onset specificity is clinically important. Late-onset wheeze — wheeze that develops after the toddler years, typically in mid-childhood — has a much stronger association with persistent adult asthma than early-transient wheeze, which is largely virus-driven and often resolves. The IL-33/ST2 pathway specifically drives the eosinophilic subtype44 eosinophilic subtype
Eosinophilic asthma is characterised by airway eosinophilia, elevated FeNO, and good response to inhaled corticosteroids and IL-5/IL-13 pathway biologics
of asthma that underlies most late-onset disease.

Confirmation came from Chinese Han children (265 asthma cases, 153 controls) in Wu et al. (2021)55 Wu et al. (2021), where CT or CC genotypes were associated with lower asthma susceptibility (adjusted OR 0.584, 95% CI 0.362–0.941, p=0.027) and, in those with asthma, a lower probability of elevated FeNO at baseline (adjusted OR 0.286, 95% CI 0.125–0.652, p=0.003). The FeNO finding directly connects genotype to eosinophilic airway inflammation — the mechanism the IL-33/ST2 axis drives.

Fourteen years of IL1RL1 population genetics, most recently summarised in Savenije et al. (2011)66 Savenije et al. (2011), established that 13 of 15 IL1RL1 SNPs tested were significantly associated with circulating sST2 levels in childhood, confirming the IL1RL1 locus as the primary genetic determinant of the soluble decoy receptor that buffers IL-33 signalling.

Practical Implications

For TT homozygotes, the actionable priority is awareness and monitoring. The TT genotype represents the population baseline — lacking the extra sST2-buffering effect the C allele confers. This does not mean TT individuals are at elevated risk in an absolute sense; it means they do not carry the C allele's protective advantage. If there is any history of childhood wheeze, asthma, or respiratory allergy, checking FeNO is the highest-yield approach to characterise whether eosinophilic airway inflammation is present — and whether eosinophil-targeting treatments (ICS, dupilumab, mepolizumab, benralizumab, itepekimab) would be particularly well-matched.

For CC homozygotes, the protective genotype reduces eosinophilic airway inflammation risk through this particular regulatory element — but this doesn't override other asthma risk factors and shouldn't be interpreted as blanket protection against all respiratory allergy.

Interactions

rs13424006 and rs10208293 both tag the late-onset wheeze phenotype at the IL1RL1 locus, sitting in the same intronic region (10th intron). Their effects may reflect distinct regulatory elements — having the protective C allele at rs13424006 alongside a risk allele at rs10208293 may represent partial compensation rather than full protection.

The IL-33 ligand variant rs992969 determines how much IL-33 is produced upstream. A TT carrier at rs13424006 who also carries the IL-33 production-increasing allele at rs992969 faces elevated ligand meeting reduced buffering — a combined effect on the entire axis. The TSLP variant rs1837253 acts upstream in the same Th2 cascade; multiple risk alleles in this pathway produce substantially higher cumulative eosinophilic inflammation risk than any single SNP predicts.

rs1562444

MTNR1B MTNR1B 3'UTR Melatonin Signaling Variant

Moderate Risk Factor

MTNR1B 3'UTR — A Regulatory Switch for the Melatonin-Glucose Clock

Your body runs on a 24-hour clock, and melatonin is its timekeeper. As darkness falls, melatonin rises — signaling sleep, lowering core body temperature, and, critically, telling your pancreatic beta cells to reduce insulin secretion until morning. The MTNR1B gene encodes the melatonin receptor type 211 melatonin receptor type 2
A G-protein-coupled receptor (MT2) expressed in the brain, retina, and — importantly — pancreatic beta cells, where it directly suppresses glucose-stimulated insulin release via inhibitory cAMP signaling
(MT2) that mediates this signal in pancreatic beta cells.

The rs1562444 variant sits in the 3' untranslated region of MTNR1B — not in the protein-coding sequence, but in the regulatory tail of the gene's messenger RNA. This position places it within a region governing mRNA stability, turnover, and post- transcriptional control, including potential microRNA binding sites. Individuals carrying the G allele — the reference allele at this position, but the global minority allele (~44% worldwide, ~50% in Europeans) — show differences in circulating melatonin levels22 differences in circulating melatonin levels
Wang et al. 2019 (PMID 31815152) reported significant differences in plasma melatonin concentrations between rs1562444 genotypes
compared to AA homozygotes, consistent with altered MTNR1B expression or signaling efficiency.

The Mechanism

Variants in the 3' UTR do not change the receptor protein itself but can alter how much MTNR1B is produced. This region contains AU-rich elements33 AU-rich elements
Sequences in mRNA that govern degradation rate — more instability signals mean lower protein output
that influence mRNA half-life and regulate microRNA binding. A single nucleotide change at position c.*371 (NM_005959.5:c.*371G>A) can shift the binding affinity of endogenous microRNAs or RNA-binding proteins that control MTNR1B transcript abundance.

More MTNR1B protein on beta cells means stronger melatonin-mediated suppression of cAMP and, therefore, reduced glucose-stimulated insulin secretion — particularly in the hours after sunset when melatonin is rising. The rs1562444 G allele appears to influence this expression level in the same direction as the well-established intronic risk variant rs1083096344 rs10830963
The strongest GWAS hit for fasting glucose in the MTNR1B locus, present in 28% of Europeans, with P=3.2×10⁻⁵⁰ in the original discovery; already profiled separately in this encyclopedia
, with which it shares partial linkage disequilibrium55 linkage disequilibrium
A statistical tendency for nearby variants to be inherited together, meaning alleles at one site predict alleles at nearby sites within a population
in the MTNR1B haplotype block.

The Evidence

The primary evidence for rs1562444 comes from several independent contexts. Wang et al. 201966 Wang et al. 2019
Wang P et al. Association of Melatonin Pathway Gene's Single-Nucleotide Polymorphisms with Systemic Lupus Erythematosus in a Chinese Population. J Immunol Res, 2019
genotyped 11 MTNR1B tag SNPs including rs1562444 in 495 SLE patients and 493 controls, reporting that rs1562444 genotype was associated with significant differences in plasma melatonin levels — direct biological evidence that this UTR variant modulates MTNR1B signaling output.

Robeva et al. 202377 Robeva et al. 2023
Robeva R et al. Melatonin Receptor 1B and Corticosteroid Receptor Polymorphisms in Infertile Women with Implantation Failure and Miscarriages. Front Biosci, 2023
found that G-allele- containing genotypes (AG+GG) were significantly enriched in 111 infertile women with recurrent implantation failure compared to 106 controls (19.3% vs. 3.6%, p=0.004), extending the MTNR1B signaling effect into reproductive physiology.

At the wider locus level, the landmark Prokopenko et al. 200988 Prokopenko et al. 2009
Prokopenko I et al. Variants in MTNR1B influence fasting glucose levels. Nat Genet, 2009
and Bouatia-Naji et al. 200999 Bouatia-Naji et al. 2009
Bouatia-Naji N et al. A variant near MTNR1B is associated with increased fasting plasma glucose levels and type 2 diabetes risk. Nat Genet, 2009
GWAS studies established that regulatory variation across the MTNR1B locus drives fasting glucose elevation (beta 0.06–0.07 mmol/L per risk allele) and type 2 diabetes risk (OR 1.09–1.15) in tens of thousands of Europeans.

Evidence for rs1562444's independent metabolic effect is emerging rather than established — its metabolic significance is partly inherited from its position within the MTNR1B haplotype block rather than from direct functional studies of this specific UTR position.

Practical Actions

The actionable guidance for MTNR1B G-allele carriers at rs1562444 mirrors the meal-timing interventions validated for the stronger intronic variants: eating earlier in the day, avoiding late dinners, and aligning meals with the low-melatonin window protect beta-cell function when melatonin receptor activity is elevated. Lopez-Minguez and Garaulet et al. 20181010 Lopez-Minguez and Garaulet et al. 2018
Lopez-Minguez J et al. Late dinner impairs glucose tolerance in MTNR1B risk allele carriers: a randomized, cross-over study. Clin Nutr, 2018
demonstrated in a randomized crossover trial that late dinner impairs glucose tolerance specifically in MTNR1B risk allele carriers.

Monitoring fasting glucose and HbA1c at recommended intervals provides the earliest warning if the metabolic effect begins to compound over time, especially for G-allele carriers who also eat late habitually.

Interactions

The MTNR1B locus contains several variants in partial LD that have been studied independently. The intronic rs10830963 (already profiled in Hormones & Sleep) is the strongest metabolic signal at this locus, while rs1562444 and the coding variant rs3781637 (G24E) represent additional dimensions of MTNR1B regulation. Co-carriage of rs1562444 G with the rs10830963 G risk haplotype could compound melatonin receptor overexpression in beta cells, though this specific interaction has not been formally modeled. Individuals carrying risk alleles at both rs1562444 and rs10830963 represent a subset worth considering for compound action guidance (see related SNP rs10830963 for the meal-timing intervention with the strongest published evidence).

rs157582

TOMM40 TOMM40 memory variant

Strong Risk Factor

TOMM40 — The Mitochondrial Gateway to Memory Aging

Every protein your neurons need to maintain their mitochondria must be physically imported across the outer mitochondrial membrane. The protein that guards this gateway is TOM40, encoded by TOMM4011 TOMM40
Translocase of Outer Mitochondrial Membrane 40; one of the core components of the TOM complex, the main protein-import channel for the outer mitochondrial membrane
. When TOMM40 function is compromised in neurons, mitochondrial biogenesis slows, ATP production falls, and the cellular stress that precedes neurodegeneration begins to accumulate. rs157582 sits in an intron of this gene on chromosome 19 — in the same genomic neighborhood as the better-known APOE variants — and is associated in multiple large population studies with accelerated memory decline and faster hippocampal shrinkage in aging.

The Mechanism

rs157582 is a non-coding intronic variant; it does not directly change the TOM40 protein sequence. Its effect is regulatory: the variant likely influences TOMM40 transcription levels, mRNA splicing, or local gene regulation within the densely packed APOE/TOMM40/APOC1 locus22 APOE/TOMM40/APOC1 locus
These three genes on chromosome 19q13.32 are in partial linkage disequilibrium; rs157582 in TOMM40 is physically between APOE and TOMM40 and multiple variants in this region are in LD with APOE ε4, though rs157582 has documented effects independent of APOE ε4 genotype
.

The mechanistic link between TOMM40 variation and neurodegeneration involves two converging pathways. First, Alu retrotransposon insertions in TOMM40 introns33 Alu retrotransposon insertions in TOMM40 introns
Alu elements are primate-specific repetitive DNA sequences; at least one Alzheimer's disease-associated TOMM40 variant originated from an Alu insertion event per Larsen et al. 2017
can disrupt mRNA processing via aberrant splicing and A-to-I RNA editing, reducing functional TOM40 protein levels in neurons. Reduced TOM40 channel activity impairs mitochondrial protein import, lowering the supply of nuclear-encoded subunits needed for oxidative phosphorylation complexes. In neurons — which are uniquely dependent on mitochondria given their high energy demands and inability to rely on glycolysis — even partial TOM40 insufficiency accelerates mitochondrial dysfunction with age.

Second, dysfunctional TOMM40 activates neuroinflammatory cascades. TOMM40 loss-of-function experiments in microglial cells demonstrate NF-κB pathway activation and NLRP3 inflammasome assembly44 NF-κB pathway activation and NLRP3 inflammasome assembly
The NLRP3 inflammasome is a multiprotein complex that cleaves pro-IL-1β into its active form; its chronic activation in microglia is a hallmark of late-stage Alzheimer's disease pathology
, with downstream secretion of IL-1β, IL-6, and TNF-α — pro-inflammatory cytokines that damage synapses and accelerate amyloid-beta and tau accumulation.

The Evidence

The largest genetic epidemiology evidence comes from a GWAS of aging-related verbal memory55 GWAS of aging-related verbal memory
Combined analysis of the Health and Retirement Study (HRS; N=7,486 genotyped) and English Longitudinal Study of Ageing (ELSA; N=6,898), measuring longitudinal immediate and delayed recall performance across multiple waves
. rs157582 reached genome-wide significance for both immediate recall change after age 60 (meta-analysis p=8.3×10⁻¹⁰) and delayed recall level (p=7.0×10⁻⁹). Critically, conditional analyses demonstrated that the signal on immediate recall change was driven by TOMM40, not APOE — providing the strongest evidence to date that rs157582 influences memory aging through a partially independent mechanism.

Brain imaging corroborates this. In 602 non-demented elders from the ADNI cohort66 602 non-demented elders from the ADNI cohort
Alzheimer's Disease Neuroimaging Initiative; participants were non-Hispanic Caucasian adults without dementia at enrollment, followed longitudinally for hippocampal volume changes by MRI
, the T allele of rs157582 was associated with faster hippocampal atrophy rate in a dose-dependent manner (p=1.23×10⁻⁸, genome-wide significant). T allele carriers showed lower Mini-Mental State Examination scores and higher Alzheimer's Disease Assessment Scale cognitive subscale scores, linking the genetic signal directly to measurable brain and cognitive outcomes even before dementia onset.

rs157582 also reaches genome-wide significance in GWAS for Alzheimer's disease itself (OR≈2.73 overall; OR≈2.83 in males, OR≈2.62 in females per Nazarian et al. 2019, PMID 30636644), and for CSF biomarkers of AD pathology including reduced Aβ1-42 (p=1×10⁻²³) and altered tau-to-Aβ ratios — the cerebrospinal fluid fingerprint of amyloid accumulation and tau hyperphosphorylation that precedes clinical dementia by 10–15 years.

Practical Actions

For T allele carriers, the actionable focus is on interventions that specifically support mitochondrial efficiency and limit neuroinflammatory burden. The critical distinction from APOE-focused guidance is timing: because the TOMM40 effect is strongest on rate of change in memory during normal aging (not just dementia risk), interventions aimed at preserving mitochondrial health are relevant decades before any cognitive symptoms emerge.

Coenzyme Q10 (ubiquinol form) directly supports mitochondrial electron transport chain efficiency. Creatine monohydrate supplements the phosphocreatine system that neurons use to buffer ATP supply during periods of high demand. DHA (the omega-3 found predominantly in brain tissue) maintains mitochondrial membrane fluidity, directly affecting the TOM complex's functional environment. These are not generic brain-health supplements — they target the specific mitochondrial import and energy production pathway that TOMM40 variation affects.

Cognitive baseline testing provides a concrete neurological anchor: by documenting your memory performance in your 40s or 50s (using validated tools like the MoCA, or a comprehensive neuropsychological battery), you create a personalized reference point that makes meaningful future change detectable much earlier than population-average norms allow.

Interactions

rs157582 sits within the APOE/TOMM40/APOC1 locus on chromosome 19, meaning it is in partial linkage disequilibrium with APOE ε4 (rs429358 T allele) and APOE ε2 (rs7412 T allele). However, the GWAS evidence distinguishes the signals: APOE drives delayed recall level, while TOMM40 rs157582 drives the rate of immediate recall change after age 60. Individuals carrying both rs157582 T and APOE ε4 face additive neurodegeneration risk from two partially independent biological pathways — amyloid clearance (APOE) and mitochondrial protein import efficiency (TOMM40). rs2075650, another TOMM40 intron variant, shows similar associations with cognitive aging but some fine-mapping studies suggest its signal is attributable to APOE ε4 LD rather than independent TOMM40 effect. rs10524523, the TOMM40 intron 6 poly-T variable-length polymorphism, has been studied specifically for Alzheimer's age-of-onset prediction and represents a separate functional variant in the same gene.

rs162049

MTRR

Emerging Risk Factor

MTRR rs162049 — The B12 Reactivation Support Variant

MTRR (methionine synthase reductase) is the enzyme that keeps methionine synthase 11 MTR: the enzyme that converts homocysteine to methionine using methylcobalamin (active B12) as a cofactor (MTR) running. During normal operation, MTR oxidizes its methylcobalamin cofactor to an inactive cob(II)alamin form — and it is MTRR's job to reduce it back to active methylcobalamin so MTR can continue. Without functional MTRR, MTR activity declines, homocysteine accumulates, and the methylation cycle slows.

rs162049 is an intronic variant — it sits within an intron of MTRR and does not directly change the protein sequence. Its significance lies in its membership in a functional haplotype: the G-allele risk haplotype is associated with reduced MTRR protein expression, meaning less enzyme is produced rather than less efficient enzyme. The net effect is the same: less B12 reactivation capacity, slower homocysteine remethylation, and altered DNA methylation patterns.

The Mechanism

The G allele at rs162049 co-segregates with a risk haplotype in the MTRR gene. Functional studies with transfected cell lines showed that this haplotype produced significantly lower MTRR protein levels compared to the wild-type haplotype, resulting in elevated homocysteine in culture medium and reduced LINE-1 methylation22 elevated homocysteine in culture medium and reduced LINE-1 methylation
Ohnami S et al. His595Tyr polymorphism in MTRR associated with pancreatic cancer risk. Gastroenterology, 2008
— a marker of global genomic methylation status. This is consistent with MTRR's central role: reduced B12 reactivation → reduced MTR activity → homocysteine accumulation → impaired one-carbon cycle → hypomethylation.

The Evidence

A multicenter Japanese case-control study 33 Ohnami S et al. Gastroenterology 2008 — 317 pancreatic cancer cases vs 1,232 controls identified rs162049 as an independent risk-associated variant for pancreatic cancer (OR 1.33, 95% CI 1.11–1.60; P = 0.0018). The association survived permutation testing under a recessive model (P = 0.024). A separate case-control study 44 Sangrajrang S et al. Breast Cancer Res Treat 2010 — 570 cases / 497 controls in Thai women in Thai women found that the G allele was associated with increased breast cancer risk in postmenopausal women (OR 1.61, 95% CI 1.07–2.44), consistent with impaired methylation-dependent gene regulation. A cross-sectional study 55 Ono H et al. Cancer Science 2012 — 384 Japanese women of 384 Japanese women found no independent effect of rs162049 on global leukocyte DNA methylation, suggesting the variant's functional impact may require haplotype context or additional environmental pressures (low B12, low folate intake).

The evidence overall is emerging: consistent biological plausibility and two independent cancer-risk associations, but no large-cohort homocysteine quantification or randomized intervention data specific to this variant.

Practical Implications

Because the G allele affects MTRR expression rather than enzyme structure, the intervention strategy focuses on reducing downstream demand rather than bypassing the enzyme. Ensuring sufficient methylcobalamin supply gives MTR more substrate to work with, partially compensating for reduced MTRR recycling capacity. Methylfolate (5-MTHF) keeps the methyl-donor pool full upstream. Monitoring plasma homocysteine provides an objective readout of whether the methylation cycle is under strain.

Interactions

rs162049 is most significant in combination with the MTRR A66G missense variant (rs1801394), which reduces enzyme efficiency — stacking reduced expression (rs162049) with reduced efficiency (rs1801394) compounds B12-reactivation impairment. Combined with MTHFR C677T (rs1801133), which limits methylfolate supply upstream, or MTR A2756G (rs1805087), which reduces methionine synthase activity directly, the effect on homocysteine clearance and DNA methylation is substantially amplified. The intronic SNP rs10380 (His595Tyr, MTRR) was co-identified in the same pancreatic cancer haplotype analysis and likely tags the same functional haplotype.

rs1800849

UCP3 -55C>T

Moderate Risk Factor

The Muscle Fat Furnace — UCP3 and Your Metabolic Set Point

Uncoupling protein 3 (UCP3) is a mitochondrial transporter found predominantly in skeletal muscle11 skeletal muscle
the largest metabolically active tissue in the body, accounting for about 40% of body mass and up to 80% of glucose disposal during exercise
, with lower expression in cardiac muscle and adipose tissue. Its primary job is to "uncouple" the proton gradient in the mitochondria from ATP synthesis, dissipating some energy as heat rather than storing it. Beyond thermogenesis, UCP3 plays a central role in fatty acid oxidation — helping the muscle burn fat rather than letting lipid intermediates accumulate and cause insulin resistance.

The -55C>T variant (rs1800849) sits in the core promoter region of UCP3, just 6 base pairs upstream of the TATA box22 TATA box
a DNA sequence that marks where transcription machinery initiates gene reading; variants here directly alter how much protein a gene produces
. Because UCP3 is encoded on the minus strand of chromosome 11, what papers call the "T allele" appears as the "A allele" in 23andMe genotype files — both refer to the same functional variant that increases UCP3 expression.

The Mechanism

This is a regulatory variant: it does not change the UCP3 protein itself, but changes how much of it is produced. Carriers of the T allele (A on plus strand)33 Carriers of the T allele (A on plus strand)
Cassell et al. discovered that skeletal muscle UCP3 mRNA expression was significantly higher in T allele carriers versus CC homozygotes (p < 0.02, n = 18)
produce measurably more UCP3 protein in skeletal muscle.

Higher UCP3 expression has several consequences: greater proton leak across the mitochondrial inner membrane, increased fatty acid oxidation44 fatty acid oxidation
the process by which the body burns fat for fuel, measured by a lower respiratory quotient (RQ)
, and reduced accumulation of toxic lipid intermediates such as diacylglycerol and ceramide. In a landmark mouse study, UCP3 overexpression completely prevented fat-induced insulin resistance55 UCP3 overexpression completely prevented fat-induced insulin resistance
Bézaire et al. showed transgenic UCP3-overexpressing mice fed a high-fat diet maintained normal insulin signaling, whereas wild-type mice developed marked insulin resistance
by keeping diacylglycerol and PKCtheta activity low.

The population frequency of this variant shows a striking geographic gradient, with higher T allele frequency in colder northern climates — consistent with selection pressure for thermogenic capacity. Northern Asian populations carry the T allele at ~45% frequency versus ~7% in sub-Saharan African populations.

The Evidence

BMI and obesity: In a UK Caucasian study of 1,009 individuals, the -55T allele was negatively correlated with body mass index66 the -55T allele was negatively correlated with body mass index
Beekman et al. Uncoupling protein 3 genetic variants in human obesity. Int J Obes, 2001
— T carriers had, on average, lower BMI than CC homozygotes, consistent with the higher fat-burning capacity conferred by increased UCP3 expression.

Type 2 diabetes: Results are ethnicity-dependent and directionally complex. A French cohort found the T allele was associated with roughly 50% reduced risk of developing type 2 diabetes77 the T allele was associated with roughly 50% reduced risk of developing type 2 diabetes
Meirhaeghe et al. An uncoupling protein 3 gene polymorphism associated with a lower risk of T2DM in a French cohort. Diabetologia, 2001
(T allele frequency 22% in controls versus 13% in T2D patients, replicated in a second cohort). However, a meta-analysis of 12 studies88 meta-analysis of 12 studies
Yu et al. Associations between UCP polymorphisms and susceptibility to T2DM. Diabetologia, 2013
found that the C allele (GG genotype in 23andMe) was associated with T2DM risk in Asian populations (OR 1.22, 95% CI 1.04–1.44) but not in European populations, and a 2021 meta-analysis found no overall association after ethnic stratification. A large Chinese rural cohort found the AA genotype associated with prediabetes99 AA genotype associated with prediabetes
Li et al. UCP2 and UCP3 variants associated with prediabetes and T2DM. BMC Med Genet, 2018
(aOR 1.68, 95% CI 1.02–2.78), particularly under a recessive model.

Dietary fat response: A clinical intervention study found that T allele carriers showed blunted improvements in insulin resistance, LDL-cholesterol, and glucose after a high-protein/low-carbohydrate diet1010 showed blunted improvements in insulin resistance, LDL-cholesterol, and glucose after a high-protein/low-carbohydrate diet
Molina-Vega et al. Effect of -55CT polymorphism of UCP3 on insulin resistance and cardiovascular risk after a high protein diet. Ann Nutr Metab, 2016
, while GG homozygotes showed robust metabolic improvements on the same diet.

Lipid profile: Paradoxically, despite the protective effects on BMI and diabetes risk, the TT genotype has been associated with higher total cholesterol and LDL-cholesterol in some studies — suggesting that the increased fat-burning may shift circulating lipid dynamics.

Practical Implications

The overall evidence picture is nuanced. The common GG genotype (coding-strand CC) is associated with lower UCP3 expression, potentially less efficient fat oxidation in skeletal muscle, and — particularly in Asian populations — greater susceptibility to insulin resistance and type 2 diabetes. For GG individuals, dietary fat composition is particularly important: diets higher in saturated fat may be less well-tolerated because the reduced UCP3 expression impairs the muscle's ability to safely oxidize incoming fatty acids, leading to greater accumulation of intramyocellular lipid intermediates.

The AG heterozygote has intermediate UCP3 expression and a moderate metabolic profile. The AA homozygote has the highest UCP3 expression and the strongest fat-oxidation capacity, though this does not provide blanket protection against all metabolic risk — and some dietary interventions (high protein, low carb) appear less effective for AA carriers.

Interactions

UCP3 interacts functionally with UCP2 (rs659366, -866G>A), which is expressed in many tissues including pancreatic beta cells and regulates insulin secretion differently from UCP3's skeletal-muscle-dominant effects. Individuals carrying both UCP2 and UCP3 promoter variants may experience compounded effects on energy balance and glucose metabolism. The UCP3 gene cluster on chromosome 11q13 is also near UCP2, and variants in this cluster have been studied as a haplotype unit in diabetes prevention cohorts. PPARGC1A (rs8192678), the master regulator of mitochondrial biogenesis and a co-activator of UCP3 expression, interacts with this variant: reduced PGC-1alpha activity from the rs8192678 Ser variant would further limit UCP3 upregulation in individuals who also carry the GG genotype at rs1800849. A compound action covering rs1800849 GG + rs8192678 TT would be appropriate if sufficient evidence exists for the combined phenotype.

CYP2C8 rs1934953 — Epoxygenase Pathway Variant

The CYP2C8 gene encodes one of the major phase I drug-metabolizing enzymes in the liver, responsible for clearing a clinically important set of medications including the chemotherapy drug paclitaxel, diabetes medications rosiglitazone and pioglitazone, and the antimalarial amodiaquine. Beyond drug metabolism, CYP2C8 plays a second, often underappreciated role: it is the primary hepatic and vascular enzyme that converts arachidonic acid into epoxyeicosatrienoic acids (EETs)11 epoxyeicosatrienoic acids (EETs)
EETs are lipid signaling molecules with vasodilatory and anti-inflammatory properties
, a family of lipid mediators that relax blood vessel walls, protect the heart, and modulate inflammation. rs1934953 sits in an intron of CYP2C8 and appears to influence this epoxygenase function.

The Mechanism

rs1934953 is an intronic variant — it does not change the CYP2C8 protein sequence. Instead, it likely acts as a regulatory variant, influencing CYP2C8 expression levels or splicing efficiency. The C allele (present at approximately 33% frequency in Europeans) has been linked to altered EET production. CYP2C8-derived EETs promote vasodilation via hyperpolarization of vascular smooth muscle and have protective effects in cerebrovascular and cardiovascular contexts. When CYP2C8 epoxygenase activity is reduced, EET levels fall and the balance shifts toward more vasoconstriction and inflammation. The C allele has also been studied in the context of carcinogen metabolism — CYP2C8 processes procarcinogens in the bladder, and altered expression may change how efficiently those compounds are activated or detoxified.

The Evidence

A 2017 Russian cohort study22 2017 Russian cohort study
Polonikov A et al. Contribution of CYP2C gene subfamily involved in epoxygenase pathway to hypertension. Clin Exp Hypertens, 2017
of 816 participants found that rs1934953 showed borderline significant association with essential hypertension risk (P ≤ 0.04), alongside a stronger signal from the nearby CYP2C8 variant rs7909236 (OR 2.99, 95% CI 1.39–6.44). The same group studied CYP2C8 rs1934953 in coronary heart disease33 studied CYP2C8 rs1934953 in coronary heart disease
Polonikov A et al. Polymorphisms of CYP2C8, CYP2C9, CYP2C19 and CHD risk. Gene, 2017
in 1,255 participants but found no significant independent association.

A 2015 study of subarachnoid hemorrhage44 2015 study of subarachnoid hemorrhage
Donnelly MK et al. EET metabolic pathway variants and aneurysmal subarachnoid hemorrhage outcomes. J Cereb Blood Flow Metab, 2015
demonstrated that CYP2C8 variants in the EET pathway significantly affected outcomes, with the CYP2C8*4 allele associated with 44–36% lower CSF EET/DHET levels and 2.2–2.5x higher risk of delayed cerebral ischemia — establishing the clinical relevance of CYP2C8 epoxygenase function in cerebrovascular biology.

The most striking association comes from bladder cancer: a 2022 case-control study55 2022 case-control study
Qu W et al. Impact of CYP2C8 genetic variants on bladder cancer susceptibility. Front Endocrinol, 2022
found the TT genotype strongly protective against bladder cancer (OR 0.26, 95% CI 0.14–0.47, p = 1.20E-05, codominant model), with the T allele showing consistent protection across dominant (OR 0.62) and recessive (OR 0.31) models. A 2023 Mexican population study66 2023 Mexican population study
Ambrocio-Ortiz E et al. CYP2C8 SNPs and COPD from biomass-burning smoke. Curr Issues Mol Biol, 2023
linked CYP2C8 variants including rs1934953 to COPD susceptibility in the setting of biomass-burning smoke exposure.

Practical Implications

The direct pharmacogenomic relevance of rs1934953 for specific drug dosing is not established by CPIC or DPWG guidelines — these focus on coding variants (CYP2C8*2, *3, *4). However, the broader context of CYP2C8 activity through EET production has implications for cardiovascular health and for individuals on CYP2C8-metabolized drugs. Individuals with the CC genotype have the lowest EET-producing capacity among common genotypes and may warrant closer cardiovascular monitoring. The variant's association with COPD from biomass smoke exposure suggests environmental interactions — individuals with the C allele who have heavy biomass/occupational smoke exposure may face heightened respiratory risk.

Interactions

rs1934953 operates in the same epoxygenase pathway as CYP2J2 (rs10509681) and EPHX2 variants, which together govern EET production and degradation. The functional CYP2C8*3 variant rs11572080 (p.Arg139Met) is a missense variant in strong linkage disequilibrium with haplotype blocks in the same region; individuals with both rs1934953 C allele and reduced-function coding variants in CYP2C8 may have compounded reduction in EET output. No CPIC compound-genotype recommendations currently exist for this combination, but the pathway logic is well established.

rs2016105

ELK3 ELK3 TSH regulatory variant

Strong Risk Factor

ELK3 and the Genetic Baseline of Thyroid Stimulating Hormone

Your thyroid is calibrated to a set point — a target circulating level of thyroid-stimulating hormone (TSH)11 thyroid-stimulating hormone (TSH)
TSH is produced by the pituitary gland to signal the thyroid to produce T3 and T4; it rises when thyroid output is low and falls when it is adequate
that your hypothalamic-pituitary-thyroid axis defends continuously. That set point is not the same for everyone: a substantial fraction of variation in TSH levels across individuals is genetically determined, and rs2016105 in the ELK3 gene is one of the contributors. Carriers of the rare A allele have a modestly elevated tendency toward hypothyroidism — not because their thyroid gland is diseased, but because a transcriptional regulator that fine-tunes thyroid signaling operates differently in their cells.

The Mechanism

ELK322 ELK3
ETS transcription factor ELK3; also known as NET, SAP-2, or ERP; a member of the ETS domain family recruited by serum response factor to bind serum response elements in gene promoters
encodes an ETS-domain transcription factor that alternates between repressor and activator modes depending on the cellular signaling context — specifically, it suppresses transcription in the absence of Ras activity and switches to activation when Ras signaling is present. This Ras-dependent toggle positions ELK3 at the intersection of growth factor signaling and gene expression in multiple secretory cell types, including thyroid follicular and parafollicular cells.

The rs2016105 variant sits within an intron of ELK3 on chromosome 12q23.1. Intronic variants at this position can influence gene expression by disrupting or creating regulatory elements — splicing enhancers, intronic enhancers, or RNA secondary structures — without changing the protein sequence itself. The precise molecular mechanism by which this variant modulates ELK3 activity in thyroid-relevant tissues has not been characterized at functional resolution, but the strength and consistency of its GWAS signal across multiple large cohorts establishes that the variant meaningfully alters thyroid function at the population level.

In the context of thyroid endocrinology, ELK3 participates in the Ras-Raf-1-ELK3 signaling cascade33 Ras-Raf-1-ELK3 signaling cascade
Demonstrated in medullary thyroid carcinoma cells by Ma et al. 2022: RREB1 regulates C-cell differentiation and calcitonin secretion via this pathway
that governs cell differentiation and hormone secretion in thyroid C cells. Normal ELK3 activity in these cells helps calibrate the output of calcitonin and, indirectly, the pituitary-thyroid axis set point.

The Evidence

The association between rs2016105 and hypothyroidism was identified in the VA Million Veteran Program GWAS44 VA Million Veteran Program GWAS
Verma A et al. Science 2024 — diversity and scale analysis of 2,068 traits in 635,969 U.S. veterans across four ancestry groups
, one of the most ethnically diverse GWAS cohorts ever assembled. The G allele (carried by ~98% of the population) showed a protective effect of approximately β = −0.25 (p = 3×10-26 in the strongest association), meaning the rare A allele confers approximately 28% increased odds of hypothyroidism per copy (OR ≈ 1.28, derived from the logistic beta coefficient). Four independent association signals at this locus were identified across ancestry-stratified analyses, with p-values consistently between 2×10-16 and 3×10-26.

The 2025 Nature Genetics hypothyroidism mega-meta-analysis55 2025 Nature Genetics hypothyroidism mega-meta-analysis
Rand SA, Ahlberg G et al. GWAS and polygenic risk prediction of hypothyroidism; 113,393 cases, 1,065,268 controls; December 2025
identified 350 loci, including 179 previously unreported, with 29 linked through TSH. This study also analyzed 482,873 individuals for circulating TSH levels directly, establishing that many hypothyroidism loci alter the TSH set point before clinical disease develops. ELK3 is among the loci identified in this trans-ethnic effort, further corroborating the MVP finding.

The A allele is notably absent in East Asian populations (frequency ~0%) and rare in African populations (~0.6%), with the highest frequency in Europeans (~2.7%). This ancestry specificity means the variant is almost exclusively clinically relevant for individuals of European descent.

Practical Implications

For AG heterozygotes — the relevant genotype for 95% of A-allele carriers given the rarity of AA homozygosity — the absolute risk increase for hypothyroidism is modest (~28% relative). Hypothyroidism is common (lifetime prevalence ~5-10% in women, ~2-3% in men in European populations), so this translates to a shift from roughly 7% baseline lifetime risk to approximately 9% for AG carriers — an additional ~2 percentage points in absolute terms.

The clinical value lies primarily in interpretation: if you carry the A allele and have TSH levels at the upper end of the reference range, this genotype provides biological context supporting earlier treatment consideration. It also argues for periodic TSH monitoring rather than a single-timepoint assessment, since individuals with this variant trend toward hypothyroidism rather than against it.

Interactions

ELK3 sits in the thyroid-function gene network alongside FOXE1 (rs965513), DIO2 (rs225014), and DIO1 (rs11206244). Each of these influences thyroid hormone set point through different mechanisms — transcriptional regulation, receptor sensitivity, and T4-to-T3 conversion respectively. A carrier with multiple thyroid-axis risk alleles across these loci may have a compound shift in TSH baseline that individual SNP effects underestimate. No formal compound analysis has been published for ELK3 combined with these variants, but the convergent biology makes interaction effects plausible.

rs2016520

PPARD +294T>C

Moderate Risk Factor

PPARD +294T>C — The Fat-Burning Regulator at the Heart of Endurance

PPARδ11 PPARδ
Peroxisome Proliferator-Activated Receptor delta — a nuclear receptor transcription factor that binds fatty acids and drives gene expression programs for fat oxidation, mitochondrial biogenesis, and muscle fiber remodeling
is often called the "exercise factor in a bottle" — researchers found that activating it in sedentary mice produced animals with dramatically improved endurance without any training. In humans, PPARδ governs how efficiently skeletal muscle burns fat during prolonged exercise. The +294T>C variant (rs2016520) sits in the 5'UTR regulatory region of the PPARD gene and alters the binding of a transcription factor that controls how much PPARδ protein is made. It is one of the most consistently replicated genetic markers for endurance athletic performance, identified across Russian, Polish, Israeli, and Chinese athlete cohorts.

The Mechanism

The +294 position in PPARD's 5'UTR (also described as -87 relative to the start codon) is a putative Sp-1 binding site22 putative Sp-1 binding site
Sp-1 (Specificity Protein 1) is a ubiquitous transcription factor that activates gene expression by binding GC-rich motifs in promoter and regulatory regions
. The C allele alters this binding motif, increasing Sp-1 affinity and driving higher PPARD transcriptional output. In vitro reporter assays have confirmed that the C allele produces significantly higher PPARD expression than the T allele.

The downstream consequences are substantial: elevated PPARδ promotes a transcriptional program in skeletal muscle that shifts fuel use toward fatty acid oxidation33 shifts fuel use toward fatty acid oxidation
PPARδ directly regulates genes for fatty acid uptake (CD36, FABP), beta-oxidation (CPT1, ACADM, HADH), and uncoupling (UCP2, UCP3), while suppressing glucose-dependent pathways during sustained effort
. It also drives the development of type I (slow-twitch) oxidative muscle fibers, increases mitochondrial density, and improves lactate clearance efficiency. The net effect in trained C-allele carriers is a metabolic phenotype suited to prolonged aerobic effort: higher fat oxidation rates, preserved glycogen, and greater endurance capacity.

The Evidence

The landmark 2009 study by Ahmetov and colleagues44 2009 study by Ahmetov and colleagues
Ahmetov II et al. The combined impact of metabolic gene polymorphisms on elite endurance athlete status and related phenotypes. Hum Genet, 2009
genotyped 1,423 Russian athletes and 1,132 controls for 15 gene polymorphisms, identifying PPARD rs2016520 C as one of ten discrete "endurance alleles." A meta-analysis combining the Caucasian cohorts yielded OR 1.57 (95% CI 1.30–1.91, p < 10⁻⁵) for elite endurance athlete status in C-allele carriers. Notably, the frequency of the C allele increased with competitive level among endurance-sport athletes, suggesting a dose-response relationship between the allele and elite performance.

A haplotype study of 660 Polish elite athletes55 haplotype study of 660 Polish elite athletes
Cieszczyk P et al. Genomic haplotype within the Peroxisome Proliferator-Activated Receptor Delta (PPARD) gene is associated with elite athletic status. Scand J Med Sci Sports, 2015
found that rs2016520 was individually associated with overall elite athletic performance (p = 0.00002) and particularly with strength-endurance sports. Analysis of three PPARD haplotypes revealed that the A/C/C haplotype (rs2267668/rs2016520/rs1053049) was dramatically underrepresented in all elite athletes compared with controls (p < 0.000001), indicating that the T allele at rs2016520 is part of a haplotype protective against elite performance in endurance sports.

An Israeli athlete cohort study66 Israeli athlete cohort study
Eynon N et al. Is there an interaction between PPARD T294C and PPARGC1A Gly482Ser polymorphisms and human endurance performance? Int J Sports Med, 2009
found that while PPARD rs2016520 alone did not reach significance in a cohort of 155 athletes, the compound genotype of PPARD CC + PPARGC1A Gly/Gly (at rs8192678) was dramatically overrepresented in elite endurance athletes versus national-level athletes (OR 8.32, 95% CI 2.2–31.4), underscoring the importance of gene-gene interactions in elite endurance capacity.

At the clinical level, a 12-week training intervention in 168 women77 12-week training intervention in 168 women
Leońska-Duniec A et al. The polymorphisms of the PPARD gene modify post-training body mass and biochemical parameter changes in women. PLOS One, 2018
demonstrated that PPARD C-allele carriers showed significant decreases in total cholesterol and triglycerides following aerobic training — a favorable metabolic response not seen in TT homozygotes — confirming that the allele's effects are exercise-dependent and emerge with training.

Practical Actions

If you carry the C allele (CT or CC), your muscles are primed to respond to endurance training with enhanced fat-burning capacity and favorable lipid changes. Prioritize aerobic training sessions at moderate intensity (60–75% of maximal heart rate) where fat oxidation is maximized, and allow sufficient volume for the training-induced lipid benefits to emerge (studies show effects after 12+ weeks of consistent aerobic work).

If you are TT homozygous, you have the common ancestral genotype. Evidence from one study suggests TT carriers may be better responders to aerobic training in terms of VO2max improvement from a lower baseline — meaning consistent training still produces substantial aerobic gains, even though you may not carry the elite endurance advantage of the C allele.

Dietary fat quality matters for all PPARD genotypes: omega-3 fatty acids (EPA and DHA) are natural PPARδ ligands that activate the receptor, potentially amplifying the fat-oxidation program. Ensuring adequate omega-3 intake is relevant regardless of genotype.

Interactions

PPARD rs2016520 interacts powerfully with PPARGC1A rs8192678 (Gly482Ser): the compound genotype of PPARD CC and PPARGC1A Gly/Gly showed an OR of 8.32 for elite endurance status versus national-level athletes in the Israeli cohort, far exceeding what either variant contributes alone. PPARGC1A encodes PGC-1alpha, the transcriptional coactivator that physically interacts with PPARδ to drive mitochondrial biogenesis in response to exercise. PPARA (rs4253778) is a closely related nuclear receptor in the same fat-oxidation pathway — individuals carrying favorable variants at both PPARA and PPARD may have additive endurance advantages.