The MGST1 Locus — When Aspirin's Cancer Protection Depends on Your Genotype

Aspirin and NSAIDs are among the best-studied cancer-prevention agents in medicine. Regular aspirin use is associated with a roughly 30% lower risk of colorectal cancer11 roughly 30% lower risk of colorectal cancer
Meta-analyses consistently show OR ~0.69 for regular vs. non-regular aspirin/NSAID users across large population studies
— a finding replicated across dozens of cohort and case-control studies. But this protection is not evenly distributed across people. Your genome contains switches that determine whether aspirin's cancer-fighting effect works for you.

The rs10505806 variant sits in a biologically active intergenic region on chromosome 12p12.3, approximately 726 kilobases downstream of the MGST1 gene. MGST1 encodes microsomal glutathione S-transferase 1, a MAPEG-family enzyme that synthesizes prostaglandin E2 (PGE2)22 microsomal glutathione S-transferase 1, a MAPEG-family enzyme that synthesizes prostaglandin E2 (PGE2)
MAPEG = membrane-associated proteins in eicosanoid and glutathione metabolism; PGE2 is the key prostaglandin that drives colorectal cancer cell proliferation and suppresses anti-tumor immunity
. Aspirin inhibits COX-1 and COX-2 to reduce prostaglandin production — but MGST1 provides an alternative, COX-independent PGE2 synthesis route that aspirin does not directly block.

The Mechanism

The chr12p12.3 locus containing rs10505806 and the nearby rs2965667 variant appears to regulate MGST1 expression through long-range enhancer elements. When aspirin suppresses COX-derived prostaglandins, the MGST1 pathway becomes proportionally more important33 When aspirin suppresses COX-derived prostaglandins, the MGST1 pathway becomes proportionally more important
This is a bypass mechanism — if MGST1 expression is higher, the prostaglandin-driven tumor promotion signal is harder to silence with aspirin alone
. T allele carriers at rs10505806 are in the same genomic neighborhood as T/A allele carriers at rs2965667, the variant with the strongest genome-wide evidence (P = 4.6 × 10⁻⁹) for modifying aspirin-CRC protection.

PGE2 promotes colorectal cancer through multiple mechanisms: it stimulates cancer cell proliferation, suppresses cytotoxic T-cell activity, and promotes angiogenesis. The MGST1 pathway is therefore not a minor backup — it is a key determinant of whether COX inhibition by aspirin translates into actual tumor suppression.

The Evidence

Nan et al. 2015 in JAMA44 Nan et al. 2015 in JAMA conducted a genome-wide interaction study of aspirin/NSAID use with colorectal cancer risk, combining 5 case-control and 5 cohort studies with 8,634 cancer cases and 8,553 controls. The top hit on chromosome 12 (rs2965667 at chr12:17,291,799, 44 kb from rs10505806) reached genome-wide significance (P = 4.6 × 10⁻⁹) for modifying aspirin/NSAID protection. Individuals with the common TT genotype at rs2965667 gained strong protection from aspirin (OR 0.66), while the rare TA/AA carriers (4% of the population) showed a paradoxical increase in risk with aspirin use (OR 1.89, P = .002). rs10505806 is in this same intergenic locus and likely tags the same regulatory variation.

The overall aspirin benefit is substantial: regular use was associated with OR 0.69 (95% CI 0.64–0.74) for colorectal cancer prevention in the full analysis. Drew & Chan's 2021 review55 Drew & Chan's 2021 review and the 2022 JAMA evidence update66 2022 JAMA evidence update both confirm this protective association while noting that individual variation in response is likely genetic.

Practical Actions

For the ~95% of people with the common AA genotype at rs10505806, aspirin's colorectal cancer protection is expected to be robust — consistent with the overall population benefit. For T allele carriers at this locus (AT or TT), the picture is more complex: based on evidence from the neighboring rs2965667, T-carrying genotypes may not benefit from, or could even be harmed by, aspirin used specifically for colorectal cancer chemoprevention.

This does not automatically mean T carriers should avoid aspirin for cardiovascular indications — those benefits operate through different pathways (platelet TXA2 suppression) and are largely preserved. The uncertainty is specifically about the cancer prevention use case. T carriers considering aspirin for colorectal cancer chemoprevention should discuss this genetic profile with their gastroenterologist.

Interactions

The chr12p12.3 region interacts with the IL16-containing locus on chromosome 15q25.2 (rs16973225) — the second genome-wide significant aspirin-CRC interaction identified in Nan et al. 2015. Separately, aspirin's mechanism depends on COX-2 activity, so genetic variants in PTGS2 (COX-2) are relevant co-regulators. The 2024 Science Advances genome-wide interaction study (rs72833769 near TBC1D7/MTOR, and rs350047 near PTGER4) identified additional aspirin-CRC modifiers through the mTOR and PGE2 receptor pathways — functionally linked to MGST1's role in the same prostaglandin axis.

rs10514299

TMEM161B-MEF2C TMEM161B-MEF2C intergenic variant

Strong Risk Factor

The Synapse Sculptor — MEF2C and the Genetics of Depression

Not all depression has the same biological roots. For roughly one in four people of European ancestry who carry at least one copy of the T allele at rs10514299, a portion of that vulnerability may trace back to the MEF2C locus11 MEF2C locus
Myocyte enhancer factor 2C — a transcription factor that acts as a master regulator of synapse number, neuronal survival, and activity-dependent plasticity in the developing and adult brain
on chromosome 5q14.3. This variant sits within a non-coding RNA transcript (TMEM161B-DT) immediately adjacent to MEF2C, where it is thought to influence MEF2C expression levels in brain tissue. The discovery of this locus, confirmed at P = 9.99 × 10−16 in the largest depression GWAS conducted at the time, placed MEF2C at the center of the emerging genetics of mood disorders.

The Mechanism

MEF2C belongs to the MADS-box family of transcription factors and is one of the earliest-expressed MEF2 isoforms in the developing telencephalon. It plays a central role in controlling excitatory synapse number through activity-dependent synapse elimination22 activity-dependent synapse elimination
When neurons fire, calcium signaling activates calcineurin, which dephosphorylates MEF2C, switching it from a repressor to an activator of synapse-pruning genes such as Pcdh10
. In this way, MEF2C acts as a negative regulator of synaptogenesis — constraining the number of excitatory connections to maintain proper excitatory/inhibitory balance.

When MEF2C function is reduced — as the T allele at rs10514299 may cause via altered regulatory element activity — excitatory synapse pruning is impaired, disrupting cortical and hippocampal circuit calibration. Mouse models with conditional Mef2c knockout in excitatory neurons show dramatically reduced network activity, anxiety-like behavior, and cognitive deficits, mirroring phenotypes relevant to depression. MEF2C also regulates neuronal differentiation, axon guidance, and activity-dependent survival — making it a broad-spectrum orchestrator of the neural circuitry that underlies mood regulation.

Pharmacologically, this pathway is relevant beyond genetics: HDAC inhibitors33 HDAC inhibitors
Histone deacetylase inhibitors increase histone acetylation, promoting the transcription of MEF2C target genes. Valproate — used as a mood stabilizer — is a Class I/II HDAC inhibitor; part of its mood-stabilizing effect may operate through MEF2C-dependent transcription
. Valproate (valproic acid) is a known Class I/II HDAC inhibitor, and some of its therapeutic effect in mood disorders may work partly through MEF2C-dependent transcriptional activation.

The Evidence

Discovery — Hyde et al. 2016. The first large-scale depression GWAS using 23andMe data44 The first large-scale depression GWAS using 23andMe data
Hyde CL et al. Identification of 15 genetic loci associated with risk of major depression in individuals of European descent. Nature Genetics, 2016
combined data from 75,607 cases and 231,747 controls in a discovery phase, then replicated in 45,773 cases and 106,354 controls. rs10514299 reached a joint p-value of 9.99 × 10−16 — one of the most statistically robust findings in that landmark study. The locus harbored two independent significant signals (rs10514299 and rs454214), suggesting the MEF2C region contains multiple functional elements contributing to depression risk.

Validation — Howard et al. 2019. A meta-analysis of 807,553 individuals55 A meta-analysis of 807,553 individuals
Howard DM et al. Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nature Neuroscience, 2019
(246,363 cases and 561,190 controls) identified 102 independent variants and 269 genes. MEF2C appeared in 10 of the 15 most significant biological gene-sets associated with depression, including GO_EXCITATORY_SYNAPSE, GO_POSTSYNAPSE, and GO_NEURON_SPINE — reinforcing that synaptic regulation at the MEF2C locus is not incidental but central to the biology of MDD.

Functional readout — Muench et al. 2018. A neuroimaging study in alcohol-dependent patients66 A neuroimaging study in alcohol-dependent patients
Muench C et al. The major depressive disorder GWAS-supported variant rs10514299 in TMEM161B-MEF2C predicts putamen activation during reward processing in alcohol dependence. Translational Psychiatry, 2018
demonstrated that T allele carriers (n=45 patients vs. 45 controls) showed significantly greater putamen activation during reward anticipation (p=0.014) and loss anticipation (p=0.024–0.046). The putamen is a key node in the reward circuitry; its hyperactivation in T allele carriers provides a direct functional bridge between the GWAS finding and disrupted reward processing — a hallmark of depressive episodes.

Effect size context. The odds ratio for rs10514299 is approximately 1.05 per T allele in European ancestry populations — modest for any individual but consistent, replicated, and mechanistically coherent. Depression is highly polygenic; this variant represents one of dozens of genome-wide significant contributors.

Practical Implications

Carrying T alleles at this locus does not cause depression, but it does represent a real, biologically grounded increment in susceptibility — particularly through disrupted synaptic calibration and reward circuit regulation. The actionable insight is targeted: strategies that support MEF2C-pathway resilience include interventions with documented neuroplasticity effects. Physical exercise increases MEF2C expression77 increases MEF2C expression
Chen SX et al. demonstrated that voluntary running upregulates Mef2c mRNA in hippocampal neurons in rodent models, suggesting exercise partially restores MEF2C pathway activity
in hippocampal tissue via BDNF-dependent signaling, connecting one of the most evidence-supported depression interventions to this specific locus.

If prescribed a mood stabilizer, valproate (valproic acid) is mechanistically relevant because its HDAC-inhibiting activity upregulates MEF2C target gene transcription — a convergence between pharmacology and the genetic vulnerability encoded at this locus.

Interactions

The TMEM161B-MEF2C locus contains two independent GWAS signals: rs10514299 and the companion variant rs454214 (in MEF2C-AS2, the antisense RNA adjacent to MEF2C). Individuals carrying risk alleles at both loci may experience a cumulative effect from the same pathway. The broader MEF2C biology intersects with dopamine signaling: MEF2C regulates striatal synapse development, and SNPs in dopaminergic genes (rs1800497 DRD2/ANKK1, rs4680 COMT) may compound vulnerability through overlapping reward circuitry mechanisms.

rs10515237

PCSK1 PCSK1 A/G (rs10515237)

Moderate Risk Factor

The Prohormone Scissor: How PCSK1 Variants Shape Insulin and Appetite

Every time your pancreas detects rising blood sugar, it doesn't release insulin directly — it releases proinsulin, a folded precursor that needs to be cut open before it becomes active. The enzyme doing that cutting is prohormone convertase 1/3 (PC1/3)11 prohormone convertase 1/3 (PC1/3)
Encoded by the PCSK1 gene on chromosome 5; cleaves inactive prohormone precursors at specific paired basic amino acid sites to generate active peptide hormones
, the protein encoded by PCSK1. But PC1/3 doesn't work on proinsulin alone — it also cleaves POMC into α-melanocyte-stimulating hormone (α-MSH), a key satiety signal in the brain, and proglucagon into GLP-1, the incretin that amplifies insulin release after meals. A single enzyme sits at the convergence of insulin processing, appetite regulation, and incretin signaling.

The rs10515237 variant sits within an intron of PCSK1 and is in [strong linkage disequilibrium | LD r² ≈ 0.84 in Europeans; variants that are tightly co-inherited — measuring one reliably predicts the other] with rs6234 and rs6235, the non-synonymous variants encoding the Q665E-S690T haplotype in the C-terminal domain of PC1/3. When you carry the G allele at rs10515237, you very likely also carry the functional haplotype that partially reduces PC1/3 activity. This is what makes rs10515237 clinically meaningful despite being intronic itself.

The Mechanism

The Q665E-S690T amino acid pair encoded by the rs6234-rs6235 haplotype alters the C-terminal propeptide region of PC1/3. This region acts as an intramolecular chaperone that guides proper protein folding — changes to it reduce the enzyme's catalytic efficiency without abolishing function entirely. The result is a partial loss-of-function22 partial loss-of-function
Not as severe as the rare biallelic null mutations that cause severe early-onset obesity and malabsorptive diarrhea; common heterozygous variants reduce activity modestly
that propagates through three downstream pathways: slower proinsulin→insulin conversion (raising the proinsulin-to-insulin ratio), reduced POMC cleavage to α-MSH (blunting melanocortin-mediated satiety), and potentially altered proglucagon processing to GLP-1 (modifying incretin responses).

The Evidence

The founding study, Benzinou et al. 2008 (Nature Genetics)33 Benzinou et al. 2008 (Nature Genetics)
13,659 individuals of European ancestry across 8 independent cohorts; the Q665E-S690T haplotype was the strongest association signal at p = 2.31 × 10⁻¹²
, established common PCSK1 variants as genuine obesity risk loci. The rs6234-rs6235 haplotype — which rs10515237 tags — conferred a 22% increased risk of obesity per risk haplotype in the original cohorts. A companion functional experiment confirmed that the N221D mutation (rs6232) significantly impaired PC1/3 catalytic activity in cell-based assays.

The largest meta-analysis to date, Nead et al. 2015 (Human Molecular Genetics)44 Nead et al. 2015 (Human Molecular Genetics)
Up to 331,175 individuals from GWAS consortia; systematic review combining literature and custom array data from diverse ethnic backgrounds
, quantified the effect more precisely: the rs6234/rs6235 haplotype carries OR = 1.07 (95% CI 1.04–1.10, p = 3.00 × 10⁻⁷) for obesity, and rs6232 carries OR = 1.15 (95% CI 1.06–1.24, p = 6.08 × 10⁻⁶). These are modest effect sizes consistent with a polygenic contribution to a complex trait. A HuGE review and meta-analysis by Stijnen et al. 201455 HuGE review and meta-analysis by Stijnen et al. 2014
Comprehensive systematic review of all published PCSK1 association studies; examined associations with BMI, obesity, waist circumference
confirmed the associations and identified that rs6232 effects were stronger in childhood obesity than adult obesity — consistent with the enzyme's role in early growth and metabolic programming.

The metabolic consequence of reduced PC1/3 activity was directly measured by Heni et al. 2010 (BMC Medical Genetics)66 Heni et al. 2010 (BMC Medical Genetics)
1,498 German subjects with detailed OGTT and clamp studies; rs6235 minor allele frequency 25.8%
: carriers of the rs6235 risk allele had 8% higher proinsulin area-under-the-curve and elevated proinsulin-to-insulin ratio, confirming that the functional haplotype measurably impairs prohormone processing in vivo.

Practical Implications

Carriers of the G allele at rs10515237 have a modestly elevated proinsulin-to-insulin ratio — their pancreatic beta cells are secreting more precursor and less active hormone per stimulus. Over time, this can translate into reduced beta-cell efficiency and elevated long-term diabetes risk. The satiety pathway impairment compounds this: if POMC-to-α-MSH conversion is reduced, the melanocortin-4 receptor pathway fires less robustly after meals, potentially raising the threshold at which fullness signals terminate eating. This combination — less effective insulin per secretory event plus blunted satiety — is the proposed mechanism linking PCSK1 variants to excess adiposity.

From a dietary standpoint, the primary lever is glycemic load. When meals generate rapid postprandial glucose peaks, beta cells secrete larger proinsulin pulses — and in PCSK1 G-allele carriers, those pulses are converted less efficiently. Spreading glucose load across the day by choosing lower-glycemic foods reduces the secretory burden on an already less-efficient processing system. Monitoring fasting proinsulin (where available) and the standard proinsulin:insulin ratio can detect beta-cell stress before HbA1c rises.

Interactions

The rs10515237/rs6234-rs6235 haplotype interacts additively with rs6232 (the N221D missense variant) within the same gene. Individuals carrying both the rs6232 risk allele and the rs6235 haplotype have two independently acting reductions in PC1/3 activity — compound heterozygosity within PCSK1 is plausible and would confer greater proinsulin accumulation. The melanocortin pathway interaction is also relevant: variants in MC4R (rs17782313) and POMC itself further modify satiety signaling downstream of PC1/3 cleavage, creating a polygenic appetite-regulation score that several GWAS consortia have begun characterizing.

The Endurance Athlete's Genetic Edge

Elite endurance athletes—those competing in marathons, cycling road races, cross-country skiing, and football matches requiring sustained high-intensity output for 90+ minutes—often share a common genetic signature. Among the most significant is rs1052373 in the MYBPC3 gene, which encodes cardiac myosin-binding protein C11 cardiac myosin-binding protein C
a regulatory protein that fine-tunes the force and speed of heart muscle contraction
. A landmark genome-wide association study22 landmark genome-wide association study
Ahmetov et al. 2020. Meta-analysis of 1,206 elite European, Russian, and Japanese athletes
found that individuals with the GG genotype were 2.2 times more likely to become elite endurance athletes compared to those with AA or AG genotypes (P = 1.43 × 10⁻⁸). Among Russian elite athletes, GG homozygotes showed significantly higher VO₂max—the gold-standard measure of aerobic capacity—than AG or AA carriers (P = 0.005).

The Mechanism

Though rs1052373 is a synonymous variant33 synonymous variant
it doesn't change the amino acid sequence of the protein
—substituting one glutamic acid codon (GAG) for another (GAA) at position 1096—it appears to influence gene expression or splicing efficiency in ways that enhance cardiac adaptation to endurance training. MYBPC3 normally acts as a molecular brake on cardiac contraction: when phosphorylated during exercise, it releases myosin heads from their "super-relaxed" state, allowing them to bind actin and generate force. The G allele may subtly modulate this regulatory balance, enabling more efficient oxygen delivery during sustained high-output performance.

MYBPC3 deficiency studies44 MYBPC3 deficiency studies
Mamidi et al. 2022. Shows MYBPC3 loss activates NF-κB pathway, reduces inflammation, and shifts metabolism toward fatty acid oxidation
have revealed that partial loss of cMyBP-C function reduces cardiac inflammation and enhances fatty acid oxidation—a more efficient fuel source during prolonged exercise. This metabolic shift may explain why GG carriers show superior endurance capacity: their hearts can sustain high cardiac output longer without depleting glycogen stores or accumulating lactate.

The Evidence

The initial discovery came from a GWAS55 discovery came from a GWAS
Miyamoto-Mikami et al. 2020. Analyzed 476,728 SNPs in 796 European elite athletes, replicated in Russian and Japanese cohorts
comparing athletes in high-aerobic sports (marathon, cycling, cross-country skiing) versus low/moderate-aerobic sports (sprinting, jumping, throwing). The rs1052373 GG genotype emerged as the strongest genome-wide significant hit. Validation in independent cohorts from Russia (n=410) and Japan (n=466) confirmed the association, with the combined meta-analysis yielding an odds ratio of 2.17 (95% CI: 1.67–2.84).

A 2023 follow-up metabolomics study66 2023 follow-up metabolomics study
Li et al. 2023. Metabolite profiling in 120 elite Chinese athletes
linked the G allele to elevated levels of androstenediol (3β,17β) disulfate (P = 1.82 × 10⁻⁵), a testosterone precursor involved in steroid metabolism. This suggests the variant influences both cardiac contractility and hormonal pathways supporting muscle recovery and adaptation. Four metabolites—quinate, theophylline, decanoylcarnitine, and ursodeoxycholic acid—were also associated with MYBPC3 expression and endurance phenotypes, though their causal roles remain under investigation.

A 2023 comprehensive review77 A 2023 comprehensive review
Pickering & Kiely. Listed rs1052373 G among the seven most promising genetic markers for endurance athlete status
, alongside PPARGC1A rs8192678 (mitochondrial biogenesis) and AMPD1 rs17602729 (purine metabolism). Importantly, the rs1052373 association held across multiple ethnic populations, suggesting a fundamental rather than population-specific effect.

Practical Actions

If you carry one or two G alleles, you possess a genetic advantage for endurance performance. Your heart is likely more efficient at sustaining high cardiac output during prolonged exercise, and your metabolic profile may favor fat oxidation over glycogen depletion. This doesn't guarantee elite status—training, nutrition, psychology, and opportunity all matter—but it does suggest your physiology is well-suited to endurance disciplines.

Training optimization: The GG genotype responds especially well to high-intensity interval training88 high-intensity interval training
4×4 min at 90-95% HRmax with 3 min active recovery improves VO₂max more than steady-state training
. Prioritize intervals of 3–5 minutes at near-maximal aerobic speed, which have been shown to produce the largest VO₂max gains. Your genetic advantage is maximized when you push your heart's upper limits.

Monitoring cardiac adaptation: Use heart rate variability (HRV)99 heart rate variability (HRV)
parasympathetic-mediated recovery indicator
to track training adaptation. Endurance athletes typically show higher resting HRV and faster heart rate recovery. When HRV drops outside your weekly baseline, scale back intensity to avoid overtraining—your cardiac system may be reaching its adaptive ceiling despite your genetic edge.

Fuel strategy: The G allele's association with enhanced fatty acid oxidation suggests you may benefit from periodized carbohydrate intake—training fasted or low-carb to upregulate fat-burning enzymes, then fueling with carbs for competition. However, this remains speculative; no studies have directly tested nutrition interventions based on MYBPC3 genotype.

Interactions

The endurance phenotype is polygenic—no single SNP determines performance. rs1052373 works synergistically with other endurance-associated variants including PPARGC1A rs81926781010 PPARGC1A rs8192678
enhances mitochondrial biogenesis and oxidative capacity
, ACTN3 rs18157391111 ACTN3 rs1815739
XX genotype (absence of alpha-actinin-3 fast-twitch protein) associated with endurance over power
, and PPARA rs42537781212 PPARA rs4253778
regulates fatty acid oxidation pathways
. If you carry beneficial alleles at multiple loci, the cumulative effect may be substantial.

Conversely, the MYBPC3 gene is also implicated in hypertrophic cardiomyopathy1313 hypertrophic cardiomyopathy
pathogenic truncating mutations cause HCM
, a disease characterized by abnormal thickening of the heart muscle. rs1052373 itself is classified as benign1414 classified as benign
all six ClinVar submissions rate it benign for HCM
, but it raises an interesting paradox: variants in the same gene can cause both pathological hypertrophy (disease) and adaptive hypertrophy (athlete's heart). The distinction lies in whether the heart thickens symmetrically with preserved function (adaptive) or asymmetrically with impaired relaxation (pathological). If you have a family history of HCM, genetic counseling is warranted regardless of your rs1052373 status.

rs1054564

GDF15 GDF15 3'UTR rs1054564

Moderate Risk Factor

GDF15 rs1054564 — The MicroRNA Switch That Sets Your Pregnancy Nausea Threshold

Growth Differentiation Factor 15 (GDF15) is perhaps the most dramatic hormone you've never heard of. Outside pregnancy its circulating levels sit quietly in the low hundreds of picograms per millilitre. During the first trimester, largely driven by the feto-placental unit, they surge 10- to 100-fold11 10- to 100-fold
Plasma GDF15 rises from roughly 200–600 pg/mL outside pregnancy to 10,000–100,000 pg/mL in the first trimester; the highest levels are reached around 7–12 weeks
. This hormone acts on GFRAL receptors22 GFRAL receptors
Glial cell-derived neurotrophic factor family receptor alpha-like — expressed almost exclusively in the area postrema and nucleus tractus solitarius of the hindbrain; the brain's primary nausea and vomiting control centre
in the brainstem, triggering nausea, food aversion, and vomiting. Your rs1054564 genotype determines how much GDF15 you make outside pregnancy — and that baseline sets your sensitivity when the first-trimester surge arrives.

The Mechanism

The rs1054564 variant sits in the 3' untranslated region (3'UTR) of GDF15, a stretch of mRNA that does not encode protein but governs how much of the mRNA is translated. The G allele (carried by ~89% of people globally) maintains the canonical stem-loop secondary structure at this position, allowing hsa-miR-1233-3p33 hsa-miR-1233-3p
A microRNA that binds the GDF15 3'UTR and suppresses translation; miR-1233-3p is expressed in a range of tissues including adipose, liver, and cardiovascular tissue
to bind efficiently and suppress GDF15 translation.

The C allele disrupts the stem-loop structure, reducing miR-1233-3p binding affinity and partially abolishing translational suppression44 partially abolishing translational suppression
Reporter assay in HEK293T cells showed significantly higher luciferase activity for the C-allele construct vs G-allele (P=0.034); western blot confirmed reduced endogenous GDF15 protein when miR-1233-3p mimics were transfected into G-allele-bearing cells
. The result: C-allele carriers have measurably higher circulating GDF15 at baseline, across all physiological states, throughout their lives.

The Evidence

Metabolic and cardiovascular effects: In a cohort of metabolic disease patients, Guardiola et al. 202455 Guardiola et al. 2024
Guardiola M et al., "The GDF15 3' UTR Polymorphism rs1054564 Is Associated with Diabetes and Subclinical Atherosclerosis." Int J Mol Sci, 2024
found that rs1054564 variant carriers (C-allele) had significantly higher serum GDF15 levels, higher prevalence of diabetes, and higher frequency of subclinical carotid atherosclerosis compared to GG homozygotes, independently of standard confounders. The direction is consistent with the molecular mechanism: more GDF15 → more of the hormone's downstream cardiometabolic stress signalling.

Ho et al. 201266 Ho et al. 2012
Ho JE et al., "Clinical and Genetic Correlates of GDF15 in the Community." Clinical Chemistry, 2012
analysed 2,991 participants in the Framingham Offspring Study and identified rs1054564 as one of two independent cis-expression quantitative trait loci (eQTL) at the GDF15 locus — meaning the variant directly influences how much GDF15 mRNA is made or translated in blood cells, with genome-wide-significant associations (P = 2.74 × 10⁻³²) at the combined locus signal.

Pregnancy nausea and hyperemesis gravidarum: The connection runs in the opposite direction for GG homozygotes. The landmark Fejzo et al. 2024 Nature paper77 Fejzo et al. 2024 Nature paper
Fejzo M et al., "GDF15 linked to maternal risk of nausea and vomiting during pregnancy." Nature, 2024; 625(7996):760-767
established that women with lower pre-pregnancy GDF15 levels face higher risk of hyperemesis gravidarum (HG). The proposed mechanism is a desensitisation model: women chronically exposed to higher baseline GDF15 (C-allele carriers, and notably patients with β-thalassemia, who have constitutively elevated GDF15) build tolerance to the hormone's emetic effects, such that the first-trimester surge causes less nausea. Women with lower baseline GDF15 (G-allele homozygotes) have not pre-sensitised their brainstem GFRAL receptors and react more severely to the same fetal GDF15 surge.

Consistent with this model, genetic variants that lower circulating GDF15 were specifically enriched in women with recurrent HG in family studies Fejzo et al. 2018 Nat Commun88 Fejzo et al. 2018 Nat Commun
Fejzo MS et al., "Placenta and appetite genes GDF15 and IGFBP7 are associated with hyperemesis gravidarum." Nature Communications, 2018
, and a Mendelian randomisation analysis in the 2024 Nature paper further supported that lower baseline GDF15 is causally linked to HG, not merely correlated with it.

Practical Actions

The rs1054564 genotype informs two distinct clinical contexts depending on biological sex and reproductive status:

For women planning pregnancy, the GG genotype identifies a group with lower baseline GDF15 and therefore higher susceptibility to severe pregnancy nausea. This is actionable before conception: pre-pregnancy interventions that raise GDF15 levels (metformin has been shown to increase GDF15, and a 2025 study found pre-pregnancy metformin use was associated with 70–82% reduced HG risk) may be relevant for women with GG and prior HG history, though this remains to be confirmed in prospective trials specifically targeting rs1054564.

For C-allele carriers, the higher baseline GDF15 reduces pregnancy nausea susceptibility but introduces a different consideration: chronically elevated GDF15 is a biomarker of metabolic stress, and the variant independently increases diabetes and subclinical atherosclerosis risk. Fasting glucose and carotid ultrasound are the specific screening tools that correspond to this genotype's risk profile.

Interactions

rs45543339 and rs1054221 (GDF15, same gene): These are the two common GDF15 haplotype-tagging variants identified in the Fejzo 2024 GWAS of hyperemesis gravidarum. They likely tag partially overlapping haplotypes with rs1054564. The combined effect of carrying low-expression variants at multiple GDF15 regulatory positions has not been formally analysed but would be expected to compound the reduction in baseline GDF15.

rs888663 (GDF15 promoter): rs888663 is the primary cis-eQTL signal at the GDF15 locus identified in the Ho 2012 Framingham study — a stronger signal than rs1054564, which is a secondary independent signal at the same locus. The two variants likely have independent but additive effects on GDF15 expression through different regulatory mechanisms (promoter activity vs 3'UTR translation).

rs1059369

GDF15 S48T

Moderate Risk Factor

GDF15 S48T — Reading the Body's Distress Signal

GDF15, once called macrophage inhibitory cytokine-1 or MIC-111 macrophage inhibitory cytokine-1 or MIC-1
GDF15 was independently discovered several times and given different names: MIC-1, NAG-1, PLAB, and PDF. All refer to the same gene on chromosome 19
, is one of the most versatile stress signals in human biology. Produced at low levels by most organs under normal conditions, it surges in response to cellular injury, metabolic overload, inflammation, pregnancy, and cancer. Once released into the bloodstream, it crosses into the brainstem where it binds the GFRAL receptor22 GFRAL receptor
Glial cell line-derived neurotrophic factor family receptor alpha-like — expressed almost exclusively in the area postrema and nucleus tractus solitarii, the brainstem regions that coordinate nausea and appetite suppression
, recruiting the co-receptor RET to trigger appetite suppression, nausea, and reduced food intake. It also serves as a potent biomarker: elevated circulating GDF15 consistently predicts all-cause mortality, cardiovascular events, and metabolic dysfunction. The rs1059369 (S48T) variant sits in the coding region of GDF15, and its genetic context — the haplotype it tags — modulates how much of this stress signal your body produces.

The Mechanism

The rs1059369 variant converts serine to threonine at position 48 of the GDF15 preproprotein (NM_004864.4:c.142T>A, p.Ser48Thr). Position 48 falls within the signal peptide and propeptide region33 signal peptide and propeptide region
GDF15 is synthesized as a 308-amino-acid preproprotein. The signal peptide (~28 aa) directs the protein to the secretory pathway; the propeptide (~167 aa) is cleaved by furin-like proteases before the mature 112-amino-acid dimer is secreted
— the portion of the protein that is cleaved before secretion. The serine-to-threonine substitution is chemically conservative (both hydroxyl-bearing), and the S48T change alone is not predicted to substantially alter the mature secreted protein's structure or receptor binding.

The variant's primary role is as a haplotype tag44 haplotype tag
A SNP that reliably tracks another nearby variant through linkage disequilibrium — the non-random co-inheritance of nearby genetic variants
. A three-SNP study of the GDF15 locus in hypertensive patients found that rs1059369 tags a haplotype anchored by the -3148C>G promoter variant (rs4808793), which increases GDF15 transcriptional activity 1.45-fold at baseline and 1.73-fold under adrenergic stimulation. The A allele at rs1059369 co-segregates with the G allele at rs4808793 in a haplotype associated with higher plasma GDF15 and more favorable cardiac remodeling. The functional driver appears to be the promoter variant; rs1059369 is the coding-region marker that travels with it.

The Evidence

Wang et al. (2009)55 Wang et al. (2009)
Wang X et al. The haplotype of the growth-differentiation factor 15 gene is associated with left ventricular hypertrophy in human essential hypertension. Clin Sci (Lond), 2009
genotyped 1,527 hypertensive patients for three GDF15 variants including rs1059369 (+157A>T). The haplotype carrying the A allele was associated with a 25% lower odds of left ventricular hypertrophy (OR=0.75, 95% CI 0.63-0.89, P=0.0009), and carriers showed significantly higher circulating GDF15 levels (P=0.04). This suggests that the A-allele haplotype upregulates GDF15 expression, and that moderately elevated GDF15 is cardioprotective in the hypertensive context — consistent with GDF15's known role in anti-hypertrophic signaling in the heart.

A large integrated analysis by Lemmelä et al. (2022)66 Lemmelä et al. (2022)
Lemmelä S et al. Integrated analyses of growth differentiation factor-15 concentration and cardiometabolic diseases in humans. eLife, 2022
across 14,099 individuals found that circulating GDF15 predicts all-cause mortality, cardiometabolic disease, respiratory illness, and psychiatric disorders, but Mendelian randomization did not support elevated GDF15 as a causal driver of cardiometabolic disease. Instead, GDF15 rises in response to metabolic stress (including obesity, which causally elevates GDF15, IVW pFDR = 0.0040). The picture is that of a biomarker reading, not a culprit — though the GFRAL-mediated appetite suppression pathway is real and pharmacologically targetable.

In pregnancy, where GDF15 rises dramatically and drives nausea via the GFRAL receptor, Pereira et al. (2026)77 Pereira et al. (2026)
Pereira DA et al. Functional characterization of the GWAS lead SNP rs888663 and effects of GDF15 SNPs on GDF15 levels in gestational hypertension and preeclampsia. Mol Biol Rep, 2026
found that TT genotype carriers with preeclampsia or gestational hypertension had lower GDF15 levels than healthy pregnant TT carriers (P<0.05), linking rs1059369 genotype to gestational GDF15 regulation in pathological pregnancies.

Practical Actions

For TT individuals (the common genotype), no specific GDF15-related intervention is needed. The TT genotype produces typical GDF15 stress signaling. For AT and AA carriers, the A allele haplotype is associated with higher circulating GDF15 at baseline, which in the context of metabolic stress (obesity, hypertension, inflammatory conditions) could influence appetite regulation, nausea sensitivity, and cardiac remodeling. Notably, GDF15 is one of the mechanisms by which metformin reduces appetite and body weight — carriers with higher GDF15 may have a more pronounced response to this mechanism.

Because GDF15 is fundamentally a stress-response biomarker, the most actionable insight is using it as a monitoring tool: elevated circulating GDF15 (above ~600 pg/mL in non-pregnant adults) signals metabolic or inflammatory stress and warrants investigation of its cause.

Interactions

rs1059369 is in linkage disequilibrium with rs4808793 (GDF15 promoter) and rs1058587 (H202D, the major GWAS variant for circulating GDF15 levels). The A allele at rs1059369 co-travels with the protective haplotype identified by Wang et al. The most important related variant is rs1058587 (H202D): this variant's G allele is present in approximately 21% of Europeans and substantially affects assay-measured GDF15 levels due to epitope interference — a nuance relevant to interpreting commercial GDF15 blood tests rather than biology per se. For individuals tracking GDF15 as a longevity biomarker, the rs1058587 H202D status is more informative than rs1059369.

CFH Y402H — The Strongest Genetic Risk Factor for Macular Degeneration

The CFH Y402H variant (also called Tyr402His) is the single most important genetic contributor to age-related macular degeneration11 age-related macular degeneration
AMD is the leading cause of irreversible blindness in people over 50 in developed countries
, a progressive disease that destroys the sharp central vision needed for reading and driving. Complement Factor H is a negative regulator of the alternative complement pathway, acting as a brake on inflammatory responses. The Y402H substitution — replacing tyrosine with histidine at position 402 — sits within a critical binding domain where CFH interacts with C-reactive protein and glycosaminoglycans on cell surfaces, particularly in the retina.

The Mechanism

The histidine variant at position 402 reduces CFH's ability to bind to heparan sulfate and other glycosaminoglycans22 heparan sulfate and other glycosaminoglycans
These molecules coat the surface of retinal pigment epithelium cells and Bruch's membrane, where CFH normally regulates complement activation
in Bruch's membrane and on retinal pigment epithelium cells. This impaired binding means CFH-402H cannot effectively suppress complement activation at these sites, leading to chronic low-grade inflammation in the macula. The 402H variant also binds less effectively to malondialdehyde33 malondialdehyde
MDA is a lipid peroxidation product that accumulates with aging and oxidative stress
, a common lipid peroxidation product that accumulates in drusen — the hallmark yellow deposits beneath the retina in AMD. The result is uncontrolled complement-mediated damage to photoreceptors and retinal pigment epithelium, culminating in geographic atrophy (dry AMD) or choroidal neovascularization (wet AMD).

The Evidence

The CFH Y402H association with AMD represents one of the most robust findings in complex disease genetics. The Rotterdam Study44 The Rotterdam Study
5,681 participants with up to 10 years of follow-up
found that CC homozygotes have an 11-fold increased risk of late AMD compared to TT individuals, with cumulative risks of vision-threatening disease by age 95 reaching 48.3% for CC, 42.6% for TC, and 21.9% for TT. Population-attributable risk was calculated at 54%, meaning more than half of AMD cases in populations of European descent can be traced to this variant. A systematic meta-analysis55 A systematic meta-analysis
Combined data from 8 studies
confirmed a multiplicative model where each C allele increases AMD odds by approximately 2.5-fold, with highly consistent effects across Caucasian populations.

The variant's clinical significance extends to treatment response. A meta-analysis of anti-VEGF therapy66 A meta-analysis of anti-VEGF therapy
Included 1,510 patients with neovascular AMD
for wet AMD found OR 1.68 (95% CI 1.09-2.60) for CC vs TT treatment response; multiple subsequent analyses show CC patients require more frequent injections and achieve more variable outcomes. Patients with CC genotype77 Patients with CC genotype
Analysis from multiple treatment cohorts
required a mean of 10.8 intravitreal injections over 12 months versus 7.2 for TC/TT genotypes. The strength of association varies substantially by ethnicity: highly significant in Europeans (C allele frequency ~36%), weaker in East Asian populations88 East Asian populations
C allele frequency ~7% in Japanese and Korean cohorts
where other variants like CFH I62V (rs800292) play a larger role.

Practical Implications

If you carry one or two C alleles, your AMD risk is meaningfully elevated, but AMD is not inevitable — onset typically occurs after age 60, and environmental factors modulate risk substantially. The most critical modifiable factor is smoking: smokers with CC genotype99 smokers with CC genotype
Rotterdam Study data
have a 34-fold increased risk compared to TT non-smokers, versus 11-fold for CC non-smokers. Quitting smoking at any age reduces risk, and the benefit applies regardless of genotype.

For CC homozygotes and high-risk TC heterozygotes, annual dilated eye exams starting at age 50 are prudent, with more frequent monitoring (every 6 months) if early drusen or pigmentary changes appear. Self-monitoring with an Amsler grid1010 Self-monitoring with an Amsler grid
A simple checkerboard pattern test that can detect early distortions in central vision
at home can catch sudden changes indicating conversion to wet AMD, where urgent treatment can preserve vision. AREDS2 supplementation (vitamin C, vitamin E, zinc, copper, lutein, and zeaxanthin) reduces progression risk by approximately 25% in individuals with intermediate AMD, though evidence for benefit in those without existing disease is weaker.

Dietary patterns may play a role: foods rich in lutein and zeaxanthin — particularly leafy greens and colorful vegetables — as well as fatty fish rich in omega-3s, have been studied in the context of AMD risk and CFH genotype1111 the context of AMD risk and CFH genotype, though evidence for a protective effect specifically in high-risk CFH genotypes remains preliminary. Blue light exposure from screens is often cited as a concern, but evidence is weak; far more important is ultraviolet protection1212 ultraviolet protection
Chronic UV exposure contributes to oxidative damage in the retina
through high-quality sunglasses that block UV rays.

If you develop wet AMD, a meta-analysis suggests CC homozygotes may respond better to anti-VEGF therapy than other genotypes, though individual responses vary and close monitoring remains essential. Emerging complement inhibitors targeting the alternative pathway are in late-stage trials and may offer genotype-specific benefits for CFH variant carriers in the coming years.

Interactions

CFH Y402H interacts multiplicatively with variants in ARMS2 (rs10490924 A69S), the second major AMD risk locus. Individuals with high-risk alleles at both loci1313 Individuals with high-risk alleles at both loci
Both CFH CC and ARMS2 TT genotypes
— CFH CC plus ARMS2 TT (risk allele) — have synergistically elevated AMD risk exceeding the product of individual effects, suggesting convergent pathways in complement activation and extracellular matrix regulation. Interestingly, RMD subtype of AMD1414 RMD subtype of AMD
Reticular macular disease, characterized by yellow interlacing networks in the macula
shows an inverse association with CFH 402H but positive association with ARMS2 69S, hinting at distinct pathogenic mechanisms within the AMD spectrum.

The gene-environment interaction with smoking is particularly striking: the combination of CFH risk alleles and smoking1515 the combination of CFH risk alleles and smoking
Data from multiple cohorts including Rotterdam Study
amplifies risk far beyond additive expectations, likely because cigarette smoke components directly activate the complement cascade and generate oxidative stress that overwhelms the already-impaired regulatory capacity of CFH-402H. Other CFH variants including rs1410996 and rs800292 show complex haplotype effects and may refine risk prediction when considered jointly with Y402H.

Compound implications involving CFH Y402H and ARMS2 rs10490924 should be considered when both variants are present, as the combined risk profile may warrant earlier screening and more aggressive preventive measures than either variant alone would suggest.

rs1064608

MTCH2 Pro290Ala

Moderate Risk Factor

MTCH2 Pro290Ala — When Your Mitochondria Struggle to Choose Their Fuel

The mitochondria in your cells do not burn a fixed fuel. Depending on what you have eaten, your activity level, and the time of day, they switch between glucose and fatty acids — a flexibility that is central to healthy metabolism and, it turns out, to restorative sleep. MTCH2 (mitochondrial carrier homolog 211 mitochondrial carrier homolog 2
an outer mitochondrial membrane protein that regulates fuel switching and fusion dynamics
) sits at this switching point. The Pro290Ala variant (rs1064608 C allele) subtly shifts how efficiently MTCH2 performs this role, with downstream effects on body weight, energy regulation, and sleep quality.

The Mechanism

MTCH2 is embedded in the outer mitochondrial membrane and performs two linked functions. First, it regulates carnitine palmitoyltransferase 1 (CPT1), the gate-keeping enzyme that determines how much fatty acid enters the mitochondrion for oxidation. When MTCH2 is functioning normally, it fine-tunes CPT1's sensitivity to malonyl-CoA — a metabolite that signals "enough fat is being burned, slow down." A functional change in MTCH2 shifts this gate, reducing fatty acid oxidation 22 Wu et al., MTCH2 modulates CPT1 activity to regulate lipid metabolism of adipocytes, Nat Commun 2025.

Second, MTCH2 promotes mitochondrial fusion during nutrient stress. When fuel is scarce, mitochondria normally elongate (hyperfuse) to extract more energy from available substrates. MTCH2 enables this response by coupling lipid synthesis flux to the fusion machinery via lysophosphatidic acid signaling 33 Labbé et al., MTCH2 links mitochondrial fusion to lipogenesis, J Cell Biol 2021. The Pro290Ala substitution changes a conserved proline in a transmembrane helix, likely affecting the protein's conformational flexibility and its interaction with CPT1 and fusion proteins.

The Evidence

The MTCH2 locus at chromosome 11 first emerged as a significant obesity signal in 2009, when Willer et al. conducted a GWAS in >91,000 individuals and identified six new BMI loci — TMEM18, KCTD15, GNPDA2, SH2B1, MTCH2, and NEGR144 when Willer et al. conducted a GWAS in >91,000 individuals and identified six new BMI loci — TMEM18, KCTD15, GNPDA2, SH2B1, MTCH2, and NEGR1
Willer CJ et al., Six new loci associated with body mass index, Nature Genetics 2009
. MTCH2 was notable because the associated variants also implicated neuronal pathways, consistent with central regulation of energy balance.

In 2019, the largest insomnia GWAS to date — 1,331,010 individuals across multiple cohorts55 1,331,010 individuals across multiple cohorts
Jansen PR et al., Genome-wide analysis of insomnia, Nature Genetics 2019
— identified 202 loci enriched in genes expressed in hypothalamic and striatal neurons. The MTCH2 chromosomal region appeared among these loci, connecting the mitochondrial fuel-switching biology to the neuronal circuits that regulate sleep-wake transitions.

The functional mechanism was clarified by Fischer et al. 202366 Fischer et al. 2023
Fischer JA et al., Opposing effects of genetic variation in MTCH2 for obesity versus heart failure, Hum Mol Genet 2023
, who showed that the Pro290 allele (G; the reference) associates with reduced MTCH2 expression via an eQTL (rs10838738), impaired glucose oxidation, and higher lactate production — but paradoxically lower BMI. The Ala290 allele (C), carried by about 35% of Europeans, is the higher-BMI genotype. Importantly, this study also identified a cardiac liability: reduced MTCH2 function (Pro290/lower expression) is overrepresented in cardiomyopathy cases, highlighting that the same gene variant has metabolically opposing effects depending on whether the tissue relies primarily on fatty acids (skeletal muscle, resting heart) or glucose (stressed heart).

The sleep connection is biologically coherent: hypothalamic neurons that control sleep-wake cycling depend on precise mitochondrial ATP production. When fuel switching is impaired, these neurons sustain suboptimal energy delivery during the early-morning hours when the body increases metabolic demand before waking — a window where mitochondrial efficiency is most critical for sleep consolidation.

Practical Actions

For carriers of one or two C alleles (Ala290), the actionable priority is supporting mitochondrial metabolic flexibility. Ubiquinol (the active form of CoQ10) directly supports the electron transport chain and has the strongest evidence for mitochondrial supplementation. Acetyl-L-carnitine supports CPT1-mediated fatty acid entry into mitochondria — a pathway MTCH2 variants specifically affect. Time-restricted eating windows of 8–10 hours help synchronize mitochondrial fuel cycling with circadian rhythms, reducing the metabolic load on MTCH2-dependent switching overnight.

For those with two C alleles, metabolic monitoring (fasting glucose, insulin, HbA1c) is warranted, given the documented link between MTCH2 impairment and glucose oxidation deficits.

Interactions

The MTCH2 locus shows genetic correlation with FTO (rs1421085, rs9939609), which independently affects fat mass and appetite regulation. Carrying risk alleles at both loci compounds adiposity risk beyond what either gene predicts alone — a compound action covering both genotypes would be appropriate if strong evidence of interaction exists. MTCH2 also sits in a pathway with PPARGC1A (PGC-1alpha), the master regulator of mitochondrial biogenesis; variants that reduce both MTCH2 function and PGC-1alpha expression may synergistically impair mitochondrial capacity.

PIWIL1 rs10773771 — A Protective 3'UTR Variant in the Ovarian Cancer Genome Sentinel

PIWIL1 (also known as HIWI) is one of four human PIWI proteins, a conserved family of RNA-guided genome guardians that silence transposable elements in germ cells. These mobile genetic elements11 mobile genetic elements
transposons and endogenous retroviruses that can copy themselves and re-insert elsewhere in the genome, causing mutations and genomic instability
are particularly threatening in the ovary, where follicular cells and their oocytes must preserve genomic integrity across a woman's entire reproductive lifespan. PIWIL1 binds piRNAs — a class of small non-coding RNAs derived from transposon sequences — and uses them as guides to recognise and silence their source elements before they can replicate or cause double-strand breaks.

The rs10773771 C>T variant lies in the 3' untranslated region (3'UTR) of the PIWIL1 transcript. The 3'UTR is a regulatory hub where microRNAs bind to control mRNA stability and translation; a single nucleotide change here can substantially alter how tightly a given miRNA represses the gene. The T allele creates a structural shift in the local mRNA folding and modifies miRNA binding affinity to the 3'UTR, with measurable eQTL effects on PIWIL1 expression in multiple tissues.

The Mechanism

Reporter gene and bioinformatics analyses by Liu et al. 201322 Liu et al. 2013
Potentially functional genetic variants in microRNA processing genes and risk of HBV-related hepatocellular carcinoma. Mol Carcinog 2014
demonstrated that rs10773771 changes the binding ability of miRNAs to the 3'UTR of PIWIL1 in cell-based assays. The T allele reduces repression of PIWIL1 by one or more targeting miRNAs, potentially allowing higher PIWIL1 expression. Consistent with this, GTEx eQTL data shows the CC genotype is significantly associated with altered PIWIL1 transcript levels in skin (P = 1.5×10⁻¹⁰) and thyroid tissue (P = 4.9×10⁻⁶), with the T allele correlating with expression patterns distinct from C/C homozygotes. Computational modelling confirms that rs10773771 alters the mRNA secondary structure at the 3'UTR, modifying the accessibility of regulatory motifs.

In the ovary, PIWIL1 is expressed in follicular cells and is thought to protect the genomic stability of both oocytes and the somatic cells surrounding them. Higher or more appropriately regulated PIWIL1 activity, as the T allele may support, could translate into more effective transposon silencing and reduced mutagenic stress in these cells — a plausible mechanism linking the T allele to lower epithelial ovarian cancer (EOC) risk.

The Evidence

The primary clinical association comes from a three-center case-control study by Liu et al. 202333 three-center case-control study by Liu et al. 2023
Relationship between PIWIL1 gene polymorphisms and epithelial ovarian cancer susceptibility among southern Chinese women. BMC Cancer 2023
examining 288 EOC cases and 361 healthy controls from southern China. Using TaqMan genotyping and multinomial logistic regression, the study found the rs10773771 CT/TT genotypes were associated with decreased EOC susceptibility. The protective effect was most pronounced in subgroups including women aged 53 years and older, those with advanced clinical stage, high pathological grade tumours, and those with specific biomarker expression profiles. Haplotype analysis confirmed that carriers of the GTG haplotype (incorporating the T allele at rs10773771) showed significantly decreased susceptibility relative to the reference GCG haplotype.

Supporting functional evidence from the same research group and others shows that this 3'UTR variant consistently modulates PIWIL1 post-transcriptional regulation across multiple cancer contexts, including hepatocellular carcinoma, stroke, and paediatric leukaemia — suggesting the variant's functional impact on PIWIL1 expression is real, even if the downstream phenotypic consequences vary by tissue and disease context.

Important limitations: The EOC association comes from a single case-control study of 649 participants from southern China. Replication in independent cohorts, particularly in non-East-Asian populations where T allele frequencies differ (East Asian ~63%, European ~60%, African ~31%), has not been published. The evidence level is therefore emerging and the association should be interpreted cautiously until replicated.

Practical Actions

For women carrying two C alleles (CC genotype), who may have relatively lower PIWIL1 3'UTR regulatory efficiency, the primary consideration is awareness of the modest EOC susceptibility signal and whether standard gynaecological surveillance is in place. This does not warrant intensive screening beyond age-appropriate guidelines — the evidence is from one study in a specific population and the absolute risk elevation is uncertain. CT and TT carriers appear to have typical-to-reduced EOC susceptibility from this variant.

Because PIWIL1 is a transposon-silencing genome guardian, oxidative and genotoxic stressors that generate DNA double-strand breaks (tobacco smoke, ionising radiation, certain endocrine-disrupting chemicals) may compound the impact of any reduction in PIWIL1 effectiveness. Minimising genotoxic exposures is prudent for CC carriers.

Interactions

rs10773771 is independent of rs28416520, a regulatory promoter variant in the same PIWIL1 gene associated with earlier age at natural menopause via a recessive mechanism. These two variants affect PIWIL1 through distinct molecular mechanisms — rs28416520 alters promoter methylation and transcription initiation, while rs10773771 modifies 3'UTR miRNA binding and post-transcriptional regulation. Women carrying both CC at rs10773771 and AA at rs28416520 may have additive reduction in PIWIL1 activity through two independent mechanisms, though no published study has examined this combination.

rs10848087 (PIWIL1 G>A) is a companion SNP from the same 2023 EOC study that increases EOC susceptibility; it and rs10773771 may tag partially independent functional elements within the PIWIL1 locus.

The Chromatin Bridge: How One Variant Controls Both Barrier Immunity and Skin Repair

A regulatory variant on chromosome 1p36.11 sits quietly 3,757 base pairs upstream of IFNLR111 IFNLR1
interferon lambda receptor 1, also known as IL28RA — the unique subunit of the receptor for type III (lambda) interferons, critical for antiviral defense at epithelial barriers
. For years this locus was attributed to IFNLR1 based on proximity alone. Then a 2021 chromatin interactome study22 a 2021 chromatin interactome study
Shi et al. used H3K27ac HiChIP to map active enhancer-promoter contacts in primary keratinocytes and T cells, revealing long-range loops from psoriasis GWAS loci to distal gene targets
rewrote the story: in keratinocytes, the psoriasis risk signal at rs10794648 loops not to IFNLR1 but to GRHL3 — a transcription factor that directs epidermal barrier repair after immune attack. This locus may therefore exert its effect through both genes, coordinating two complementary responses when the skin faces infection or inflammatory injury.

The Mechanism

IFNLR1 angle. Type III interferons (IFN-λ1–4) activate the JAK-STAT cascade through a heterodimeric receptor: IFNLR1 pairs with IL-10Rβ to transduce the antiviral signal. Unlike type I interferons (which signal ubiquitously), IFNLR1 expression is concentrated at epithelial barriers33 concentrated at epithelial barriers
skin, gut, lung, and reproductive mucosa — precisely the surfaces where pathogen contact occurs first
. This restricted expression makes IFNLR1 the gatekeeper of mucosal-layer antiviral immunity. In psoriatic skin, IL28RA expression is significantly decreased in lesional vs non-lesional tissue44 significantly decreased in lesional vs non-lesional tissue, and IFNLR1 has been shown to directly inhibit keratinocyte proliferation through cell cycle arrest — implying that when IFNLR1 signaling is reduced, keratinocytes may proliferate pathologically.

GRHL3 angle. Grainyhead-like 3 (GRHL3) is a transcription factor essential for epidermal development in the embryo that is re-activated in adult skin after barrier disruption. In psoriasis lesions, GRHL3 expression is elevated 2.62-fold55 elevated 2.62-fold
Gordon et al. 2014, JCI: GRHL3 binds known epidermal differentiation gene targets and correlates with IL-17 and IL-22 cytokine activity in lesional skin
, representing an active epidermal repair response to immune attack. GRHL3 downstream targets include IVL (involucrin, a cornification scaffold protein), KLF4 (a keratinocyte differentiation regulator), and barrier lipid-processing enzymes. When GRHL3 is knocked out in mice66 When GRHL3 is knocked out in mice
Grhl3-conditional knockout: exacerbated imiquimod-induced psoriasiform lesions, delayed healing, resistance to anti-IL-22 therapy
, disease becomes worse and harder to treat, confirming a protective barrier-repair role. Loss of GRHL3 in differentiated keratinocytes further triggers TARC/CCL17 secretion77 triggers TARC/CCL17 secretion
a Th2 chemokine that drives basal keratinocyte hyperproliferation through a paracrine loop
, linking GRHL3 loss to both Th2 immune skewing and the hyperproliferative phenotype of inflamed skin.

The Evidence

The 1p36 locus was first identified by Strange et al. 2010 (Nature Genetics)88 Strange et al. 2010 (Nature Genetics)
"Genome-wide association study identifies new psoriasis susceptibility loci and an interaction between HLA-C and ERAP1"
, a GWAS of 2,622 psoriasis cases and 5,667 controls that reached genome-wide significance (OR ~1.22, p=6.89×10⁻⁸) with the lead SNP rs4649203 in the same LD block as rs10794648. The GWAS Catalog records rs10794648-C as the risk allele for psoriasis with OR ~1.2 and p-values in the range of 1×10⁻⁶–2×10⁻⁹ across replication studies.

The chromatin looping evidence (PMID 33607115) repositioned the functional target: HiChIP data in HaCaT keratinocytes showed direct physical contact between the 1p36 risk locus and the GRHL3 promoter, while no loop was detected to IFNLR1. GRHL3 was independently associated with psoriasis susceptibility before the chromatin data — it is among only eight genes simultaneously associated with both altered psoriasis gene expression and a GWAS susceptibility signal, underscoring a causal connection.

Practical Actions

Carriers of the CC or CT genotype at rs10794648 carry a modest increase in psoriasis risk. The C allele is the population-majority allele (~75% European frequency), so most people have at least one copy. The absolute OR (~1.2 per allele) translates to a meaningful shift in background lifetime risk only when combined with other psoriasis susceptibility loci and environmental triggers (skin trauma, streptococcal infection, certain medications, stress). The TT genotype — two copies of the T (protective) allele — is present in only ~6% of Europeans and is associated with the lowest locus-attributable psoriasis risk.

For individuals at elevated psoriasis risk (family history, early lesions, confirmed streptococcal triggers), the biology at this locus points to two actionable areas: (1) supporting IFN-λ mucosal antiviral defenses, particularly against skin-tropic viruses; and (2) protecting the physical barrier against the disruption events — skin trauma, harsh detergents, wet work — that activate the GRHL3 repair pathway.

Interactions

The IFNLR1/GRHL3 locus interacts biologically with the broader type III interferon pathway. The IFNL3 locus (rs12979860, rs8099917) on chromosome 19q13 affects IFN-λ production; variants reducing IFN-λ production at the source (IFNL3/4) compound with reduced receptor expression at IFNLR1 to further diminish mucosal antiviral signaling. For psoriasis specifically, this locus interacts with HLA-C (the major MHC psoriasis gene) and ERAP1 (aminopeptidase controlling peptide loading onto HLA-C): the Strange et al. study showed HLA-C × ERAP1 epistasis, and non-MHC loci including 1p36 contribute additively to polygenic psoriasis risk. Clinically, the 1p36 locus may modulate severity more in psoriasis triggered by infection (viral or streptococcal) than in purely trauma-induced disease.