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.

rs2073658

USF1 USF1 FCHL Variant

Strong Risk Factor

USF1 — The Transcription Factor at the Heart of Familial Hyperlipidemia

When researchers searched for the genetic root of familial combined hyperlipidemia — the most common inherited lipid disorder, affecting 1–2% of the population and responsible for a disproportionate share of premature coronary disease — the trail led to USF1. This gene encodes upstream stimulatory factor 111 upstream stimulatory factor 1
USF1 is a basic helix-loop-helix leucine zipper transcription factor that binds E-box motifs in the promoters of dozens of metabolic genes
, a master regulator of lipid and glucose metabolism that controls expression of ABCA1, APOA5, APOE, fatty acid synthase, and microsomal triglyceride transfer protein (MTP), among others. The rs2073658 variant sits within an intron of USF1 but has measurable effects on how the gene responds to insulin — effects that ripple outward to triglyceride secretion and cardiovascular risk.

The Mechanism

rs2073658 does not change the USF1 protein directly; instead it sits within a FOXA2 binding site22 FOXA2 binding site
FOXA2 (forkhead box protein A2) is a transcription factor that itself regulates USF1 transcription; the two proteins form a feed-forward regulatory loop
in the USF1 gene. Functional studies by Auer et al. showed that constructs carrying the major (C) allele display higher transcriptional activity than minor (T) allele constructs. When FOXA2 is knocked down, it reduces activity of major allele constructs but not minor allele constructs — indicating that the C allele sustains a feed-forward loop in which FOXA2 activates USF1 transcription and USF1 in turn activates FOXA2, driving expression of MTP and thereby hepatic triglyceride secretion.

The T risk allele disrupts this loop in a different way under metabolic conditions: the Naukkarinen 2009 study33 Naukkarinen 2009 study
Naukkarinen et al. Functional variant disrupts insulin induction of USF1: mechanism for USF1-associated dyslipidemias. Circ Cardiovasc Genet, 2009
profiled fat and muscle biopsies before and after a euglycemic hyperinsulinemic clamp in 47 and 118 individuals respectively. The risk allele of rs2073658 eradicated the normal inductive effect of insulin on USF1 expression in both tissues, leading to perturbed expression of downstream target genes in adipose tissue. The net effect: T allele carriers respond abnormally to insulin at the level of gene regulation, producing a dyslipidemias-prone transcriptional state.

The Evidence

The foundational evidence came from a 2004 Nature Genetics study44 2004 Nature Genetics study
Pajukanta et al. Familial combined hyperlipidemia is associated with upstream transcription factor 1 (USF1). Nature Genetics, 2004
of 60 extended Finnish FCHL families comprising 721 genotyped individuals. Association between USF1 haplotypes and FCHL reached p=0.00002 overall, with a striking p=0.0000009 in males with elevated triglycerides. Carriers of the risk USF1 haplotype showed altered expression of USF1 target genes in fat tissue. This was the first gene definitively associated with FCHL, a disorder previously known only by its phenotype.

Independent replication followed promptly. A study in 314 individuals from 24 Mexican FCHL families55 314 individuals from 24 Mexican FCHL families
Huertas-Vazquez et al. Familial combined hyperlipidemia in Mexicans. Arterioscler Thromb Vasc Biol, 2005
found significant association between rs2073658 and FCHL and triglyceride traits (p=0.0009 for the strongest association), providing independent cross-ethnic evidence. A large Utah pedigree study of 2,195 subjects across 87 families replicated the association with FCHL, LDL cholesterol, and triglycerides (p=0.001–0.05), with the strongest effects in males.

The biological stakes were clarified by a 2016 Science Translational Medicine study66 Science Translational Medicine study
Laurila et al. USF1 deficiency activates brown adipose tissue and improves cardiometabolic health. Sci Transl Med, 2016
showing that individuals carrying variants that reduce USF1 expression have improved insulin sensitivity, a favorable lipid profile (higher HDL, lower TG), and reduced atherosclerosis burden — consistent with T allele carriers having reduced but dysregulated USF1 activity under insulin signaling.

Practical Actions

For T allele carriers, the key levers are reducing hepatic triglyceride production triggers and supporting the insulin-responsive gene regulation pathway that rs2073658 impairs. Limiting dietary saturated and trans fat reduces MTP substrate and hepatic VLDL assembly. Omega-3 fatty acids (EPA/DHA) independently reduce hepatic triglyceride synthesis and VLDL secretion via PPAR-alpha activation, providing a complementary pathway bypass. Monitoring fasting triglycerides, LDL particle number, and apoB is more informative than total cholesterol alone for this variant, since FCHL involves elevated apoB-rich particles across multiple lipoprotein fractions.

TT homozygotes carry two copies of the risk haplotype and have the highest liability for developing the full FCHL phenotype; a lipid specialist assessment is appropriate if triglycerides or LDL remain elevated after dietary optimization.

Interactions

USF1 regulates APOA5, APOE, and ABCA1 — genes with their own common variants in the GeneOps database. The APOA5 rs662799 variant (APOA5*3 haplotype) reduces APOA5 expression and elevates triglycerides; combined carriage of USF1 T and APOA5*3 likely compounds triglyceride burden through distinct but additive mechanisms. The related USF1 SNP rs3737787 (in the 3' UTR) tags the same risk haplotype and was the primary variant studied in several of the replication cohorts; both rs2073658 and rs3737787 track the same USF1 risk haplotype and are in high linkage disequilibrium in European populations.

NRF1 rs2402970 — The Aerobic Baseline Variant

Nuclear respiratory factor 1 (NRF1) is the master transcription factor that executes the mitochondrial biogenesis program — converting the upstream signal from PGC-1alpha into actual transcription of the nuclear genes that build the electron transport chain, import proteins into the mitochondrion, and replicate mitochondrial DNA. NRF1 binds directly to the promoters of TFAM, cytochrome c, and all five respiratory complex subunit genes, making it the essential link between the cell's energy-sensing machinery and the physical manufacture of new mitochondria.

The rs2402970 polymorphism lies within an intron of NRF1 on chromosome 7 at position 129,739,961 (GRCh38 plus strand). The C allele is the major allele globally (~83%) and is associated with higher baseline aerobic efficiency. The T allele is the minor allele (~17% globally; ~12% in Europeans, ~27% in Africans), and is associated with lower ventilatory threshold and poorer running economy at baseline — before any training intervention. This distinguishes rs2402970 from the companion NRF1 variant rs6949152: rs6949152 primarily predicts training response (how much your aerobic capacity improves with endurance training), while rs2402970 predicts baseline aerobic function (where you start from). Together they describe two distinct facets of NRF1 activity in aerobic physiology.

The Mechanism

rs2402970 is an intronic variant with no protein-coding consequence. Its molecular mechanism has not been directly characterized, but intronic variants at positions embedded deep within large introns — as this one is (c.1348+12596C>T per Ensembl annotation) — can influence pre-mRNA splicing efficiency, regulatory element occupancy, or RNA secondary structure in ways that alter mature transcript levels. Consistent with a transcriptional-output effect, the phenotypic pattern in exercise studies is a graded baseline difference across genotypes rather than a binary loss of function.

NRF1's downstream targets explain why a subtle reduction in its transcriptional output manifests as reduced aerobic efficiency specifically: TFAM (the mitochondrial transcription factor A) determines mitochondrial genome copy number; cytochrome c is the electron shuttle between complexes III and IV; and the nuclear-encoded subunits of complexes I–V set the ceiling for oxidative phosphorylation capacity. A T-allele-driven reduction in NRF1 activity would compress the entire downstream cascade, resulting in fewer mitochondria and slightly less efficient oxidative phosphorylation per unit of muscle mass — manifesting as a lower ventilatory threshold and higher metabolic cost at any given running speed.

The Evidence

The primary evidence comes from He et al.11 He et al.
He Z et al. NRF-1 genotypes and endurance exercise capacity in young Chinese men. Br J Sports Med, 2008
, a prospective 18-week endurance training RCT in 102 young Han Chinese male soldiers (mean age 19). Three NRF1 polymorphisms were genotyped: rs2402970, rs6949152, and rs10500120. For rs2402970, a significant genotype effect was seen for ventilatory threshold (VT, p = 0.004) and running economy (RE at 12 km/h, p = 0.027) at baseline — before any training began. These are baseline phenotype differences, not training-response interactions, meaning the genotype predicts the starting aerobic efficiency of individuals rather than how much they improve with exercise. The effect size at p = 0.004 is notably stronger than the rs6949152 signal (p = 0.047 for its training-response interaction), suggesting rs2402970 tags a functional regulatory element with a more direct effect on NRF1 output.

A secondary line of evidence comes from Taherzadeh-Fard et al.22 Taherzadeh-Fard et al.
Taherzadeh-Fard E et al. PGC-1alpha downstream transcription factors NRF-1 and TFAM are genetic modifiers of Huntington disease. Molecular Neurodegeneration, 2011
, which genotyped 15 NRF1 SNPs in more than 400 German Huntington disease patients. NRF1 variants — including rs2402970 — showed nominally significant associations with age of onset of HD motor symptoms. Because HD age of onset is partly determined by how well neurons maintain mitochondrial energy production under the toxic polyglutamine stress of mutant huntingtin, this finding independently supports the hypothesis that NRF1 transcriptional output (influenced by rs2402970) modulates mitochondrial resilience in neuronal tissue — consistent with the aerobic muscle findings but extending to brain energy homeostasis.

A 2024 neuronal study33 2024 neuronal study
Massaro M et al. Nuclear respiratory factor-1 (NRF1) induction drives mitochondrial biogenesis and attenuates amyloid beta-induced mitochondrial dysfunction and neurotoxicity. Neurotherapeutics, 2024
showed that increasing NRF1 expression in neurons under amyloid-beta stress restored mitochondrial mass, improved ATP synthesis, and reduced ROS — reinforcing that even modest variation in NRF1 activity level has functional consequences in post-mitotic cells with high and continuous energy demands.

The Williams et al. 2017 systematic review44 Williams et al. 2017 systematic review
Williams CJ et al. Genes to predict VO2max trainability: a systematic review. BMC Genomics, 2017
identified rs2402970 among candidate variants for aerobic capacity, noting limited independent replication — consistent with the moderate evidence grade assigned here.

Practical Actions

The T allele's association with lower baseline ventilatory threshold and running economy points to interventions that support NRF1 transcriptional output and compensate for reduced baseline mitochondrial density. Unlike rs6949152, where the primary deficit is blunted aerobic adaptation, rs2402970 T-carriers start from a lower aerobic baseline — which affects both endurance performance and the metabolic milieu of skeletal muscle at rest. Lower VT means the muscle shifts to anaerobic metabolism at lower exercise intensities, and poorer running economy means more oxygen is consumed for the same mechanical output.

Mitophagy activators (urolithin A) address mitochondrial quality; NAD+ precursors (NMN or NR) activate the SIRT1/PGC-1alpha pathway that coactivates NRF1; HIIT-style training provides the strongest stimulus for AMPK-driven NRF1 upregulation. For T/T homozygotes, all three approaches in combination are warranted.

Interactions

The closest interaction is with the companion NRF1 variant rs6949152. Both SNPs are intronic in NRF1 and were studied together by He et al. 2008 in the same cohort. rs2402970 predicts baseline aerobic efficiency (VT, running economy), while rs6949152 predicts training-response magnitude (VT gain over 18 weeks). A person carrying T at rs2402970 and G at rs6949152 would start with a lower aerobic baseline and also have a blunted training response — a compound disadvantage in the NRF1 biogenesis axis.

The interaction with PPARGC1A rs8192678 (Gly482Ser) operates one step upstream: the Ser482 allele impairs PGC-1alpha's ability to coactivate NRF1 and MEF2 transcription factors. When PPARGC1A Ser482 reduces the upstream coactivation signal and rs2402970 T independently reduces NRF1 baseline output, the two deficits stack at different points in the same mitochondrial biogenesis cascade.

FOXO3 rs2802292 is a secondary interaction partner through mitochondrial quality control: the FOXO3 G allele enhances mitophagy and stress resilience, partially compensating for reduced NRF1-driven biogenesis. Absence of the protective FOXO3 G allele in a T-carrier at rs2402970 leaves both mitochondrial quantity and quality under-supported.

TGFBR2 Leu308Pro — When the Body's Growth Control Goes Wrong

Every blood vessel in your body is held under tension by the interplay of growth signals and structural proteins. One of the most critical regulators of this balance is TGF-β signaling11 TGF-β signaling
the transforming growth factor-beta pathway controls cell proliferation, extracellular matrix production, and tissue repair in nearly every organ system
. TGFBR2 encodes the type II receptor for TGF-β — the receptor that first captures the TGF-β signal and kicks off a phosphorylation cascade into the cell nucleus. The Leu308Pro variant (rs28934568, ClinVar VCV000012505) substitutes a leucine with a proline in the kinase domain of this receptor, disrupting the intracellular signaling machinery. The result is Loeys-Dietz syndrome type 222 Loeys-Dietz syndrome type 2
LDS2, OMIM 190182 — one of six LDS subtypes, caused by mutations in TGFBR2 and representing approximately 55-60% of all LDS diagnoses
, a multisystem connective tissue disorder in which aortic aneurysm and dissection can occur at unexpectedly small vessel diameters and at younger ages than in comparable conditions like Marfan syndrome.

The Mechanism

TGFBR2 is on chromosome 3 at position 30,672,106 (GRCh38). The T-to-C change at this position converts leucine 308 in the intracellular kinase domain to proline — an amino acid that, due to its cyclic side chain, introduces a rigid kink that disrupts alpha-helical secondary structure. The kinase domain is where TGFBR2 autophosphorylates and phosphorylates its partner receptor TGFBR1, initiating downstream Smad2/Smad3 signaling. Functional studies of related TGFBR2 kinase domain variants confirm that the pathogenic variants reduce Smad2 phosphorylation and TGF-β-induced gene transcription33 the pathogenic variants reduce Smad2 phosphorylation and TGF-β-induced gene transcription
Luo et al. 2020, in vitro assays of a de novo TGFBR2 kinase domain variant
, impairing the growth-factor circuit that normally maintains connective tissue homeostasis.

Paradoxically, affected tissues in LDS show increased TGF-β pathway markers — elevated collagen expression, increased phospho-Smad2 in nuclei — even as the mutant receptor impairs direct signaling. This paradox, first described in the original 2005 discovery paper by Loeys and colleagues44 by Loeys and colleagues
Nature Genetics, ten LDS families with TGFBR1/TGFBR2 mutations
, is thought to reflect compensatory upregulation of alternative TGF-β signaling routes that overshoots the system, driving excessive extracellular matrix remodeling in the aortic wall. This overactive matrix remodeling weakens the structural integrity of the aorta, predisposing it to aneurysmal dilation and catastrophic dissection.

Inheritance is autosomal dominant — one copy of the pathogenic variant is sufficient for disease. Approximately 75% of LDS cases arise from de novo mutations; 25% are inherited from an affected parent.

The Evidence

A systematic review of 3,896 LDS cases by Gouda et al.55 A systematic review of 3,896 LDS cases by Gouda et al.
International Journal of Cardiology, 2022
established that TGFBR1 and TGFBR2-related LDS (types 1 and 2) carry the most severe aortic phenotype among all LDS subtypes. Aortic dissection occurs at smaller diameters than in Marfan syndrome — a critical clinical distinction. The peripartum aortic dissection rate among 222 pregnant LDS patients was 4%, with 1% peripartum mortality.

GeneReviews management guidelines66 GeneReviews management guidelines
Loeys & Dietz 2008, updated 2024; NCBI Bookshelf NBK1133
specify surgical thresholds of approximately 4.0 cm maximal aortic diameter for TGFBR1/TGFBR2-related LDS — lower than the 5.5 cm threshold used for the general population and the 4.5 cm used for SMAD2/SMAD3-related LDS. This conservative threshold reflects the documented tendency of TGFBR2 variant carriers to dissect at smaller sizes.

Velchev and colleagues77 Velchev and colleagues
2021, Advances in Experimental Medicine and Biology
describe the full phenotypic spectrum: arterial tortuosity extends throughout the vascular tree — not just the aortic root — and intracranial, thoracic, and abdominal aneurysms can develop independently. This makes comprehensive arterial imaging essential beyond echocardiography alone.

Practical Actions

Management has four pillars. First, cardiovascular surveillance: annual echocardiography to track aortic root size, with MRA or CT angiography every two years (or annually if growth is detected) to assess the entire arterial tree from head to pelvis. Second, medical therapy: beta-adrenergic blockers or angiotensin receptor blockers (ARBs such as losartan) reduce hemodynamic wall stress and are prescribed from diagnosis in all carriers. Third, activity restriction: contact sports, competitive sports, isometric exercise (heavy weightlifting), decongestants, and triptans (migraine medications) must be avoided. Fourth, family cascade screening: each first-degree relative has a 50% inheritance probability and requires molecular testing or comprehensive cardiovascular evaluation.

Elective surgical repair at aortic diameters approaching 4.0 cm protects against the elevated dissection risk seen at smaller sizes in TGFBR2-related disease. Cardiothoracic surgical planning should begin well before this threshold.

Interactions

TGFBR2 acts in the same TGF-β signaling pathway as TGFBR1 (LDS type 1), SMAD2 (LDS type 4), SMAD3 (LDS type 5), TGFB2 (LDS type 3), and TGFB3 (LDS type 6). While each gene produces a clinically distinct LDS subtype, the downstream pathophysiology — excessive aortic wall remodeling driven by dysregulated TGF-β signaling — is shared. TGFBR2 variants also overlap phenotypically with FBN1 mutations (Marfan syndrome) and COL3A1 mutations (vascular Ehlers-Danlos syndrome), which should be considered in the differential diagnosis during genetic workup when the clinical picture includes significant aortic or arterial disease.

Pregnancy represents a specific interaction: the hemodynamic load of pregnancy combined with the aortic fragility of LDS is a high-risk combination requiring proactive planning with maternal-fetal medicine and cardiology specialists before conception.

PMM2 R141H — The Most Common CDG Mutation and Its Implications for Carriers

Every protein destined for export, membrane display, or passage through the secretory pathway must be decorated with sugar chains — a process called N-linked glycosylation11 N-linked glycosylation
The attachment of oligosaccharide chains to the nitrogen of asparagine residues in proteins destined for the endoplasmic reticulum. "N-linked" refers to the nitrogen atom of asparagine that anchors the glycan.
. This process depends on a supply of activated sugar building blocks, and the enzyme phosphomannomutase 2 (PMM2) is the gatekeeper for one of them: it converts mannose-6-phosphate to mannose-1-phosphate, the precursor of GDP-mannose, which is incorporated into the oligosaccharide core. When PMM2 fails, cells cannot build complete N-glycan chains, and hundreds of glycoproteins are produced in aberrant forms — affecting brain development, liver function, blood coagulation, hormones, and more.

The c.422G>A variant (rs28936415) changes arginine to histidine at position 141 of the PMM2 protein (p.Arg141His, commonly written R141H). It is the most frequent pathogenic PMM2 allele worldwide, particularly in European populations22 particularly in European populations
European carrier frequency approximately 1/79 in the Netherlands and 1/60 in Denmark per Schollen et al. 2000
, and accounts for roughly half of all disease-causing PMM2 alleles identified in CDG patients.

The Mechanism

PMM2 functions as a homodimer — two identical PMM2 subunits join to form the active enzyme. The Arg141 residue sits in a conserved region critical for protein folding and stability. The R141H substitution introduces a smaller, less basic histidine in place of the bulky positively charged arginine, destabilizing the protein structure.

Protein folding studies33 Protein folding studies
Yuste-Checa P et al. The Effects of PMM2-CDG-Causing Mutations on the Folding, Activity, and Stability of the PMM2 Protein. Hum Mutat, 2015
confirmed that R141H is a destabilizing variant that substantially reduces enzyme activity, and that the protein may be partially rescuable with pharmacological chaperones — a finding relevant to emerging treatment strategies.

A critical fact distinguishes R141H from other PMM2 variants: homozygous R141H (two copies of this exact mutation) has never been observed in any living person. Population genetics analysis44 Population genetics analysis
Schollen E et al. Lack of Hardy-Weinberg equilibrium for the most prevalent PMM2 mutation in CDG-Ia. Eur J Hum Genet, 2000
confirmed that given the observed carrier frequency (1/79 in Dutch neonates), homozygous infants should appear in the population if they survived — but they don't. R141H homozygosity leaves too little PMM2 activity for embryogenesis to complete. All living PMM2-CDG patients who carry R141H are compound heterozygotes55 compound heterozygotes
Carry two different pathogenic PMM2 alleles — R141H on one chromosome and a different mutation (e.g. F119L, V44A, I132T) on the other. Neither allele alone causes disease; both are needed.
: they carry R141H on one chromosome and a different, less severe PMM2 mutation on the other that preserves just enough residual enzyme activity to allow survival.

The Evidence

PMM2-CDG (formerly CDG type Ia, or Jaeken syndrome) was the first CDG identified and remains by far the most common, with over 1,000 cases worldwide. The PMM2 gene was discovered in 199766 1997
Matthijs G et al. Mutations in PMM2, a phosphomannomutase gene on chromosome 16p13, in carbohydrate-deficient glycoprotein type I syndrome. Nat Genet, 1997
, with R141H among the original mutations identified.

A landmark French cohort study Schiff et al. 201777 Schiff et al. 2017
Clinical, laboratory and molecular findings and long-term follow-up data in 96 French patients with PMM2-CDG. J Med Genet, 2017
of 96 PMM2-CDG patients documented the full clinical spectrum: cerebellar ataxia in virtually all patients, intellectual disability, neonatal hypotonia, liver involvement, cardiac abnormalities in 20%, and a mortality rate of ~12.5% (12 deaths at mean age 3.8 years, primarily from cardiac complications). R141H was present in approximately 50% of all disease-causing alleles in this cohort. Biochemical markers consistently show elevated liver enzymes, abnormal coagulation parameters (low protein C, antithrombin), elevated TSH, and hypocholesterolemia.

ClinVar (VCV000007706) lists R141H as Pathogenic with criteria provided by over 50 independent submitters including Mayo Clinic Laboratories, ARUP, GeneDx, and Invitae — one of the most extensively classified pathogenic variants in the CDG field.

Practical Actions

For heterozygous carriers: A single R141H allele does not cause PMM2-CDG. Carriers have one fully functional PMM2 gene and maintain normal glycosylation. The primary significance is reproductive: if both partners carry a pathogenic PMM2 allele (any combination), each pregnancy has a 25% chance of inheriting compound heterozygous PMM2 variants and developing CDG.

For compound heterozygous individuals (disease): PMM2-CDG requires multidisciplinary management by metabolic physicians, neurologists, cardiologists, and hepatologists. Biochemical monitoring of transferrin glycoform analysis, liver enzymes, coagulation factors, and thyroid function is essential. While no curative therapy exists, emerging approaches include dietary mannose supplementation and repurposed drugs. Epalrestat, an aldose reductase inhibitor approved for diabetic neuropathy in Japan, showed 30–400% increases in PMM2 enzyme activity88 showed 30–400% increases in PMM2 enzyme activity
Iyer S et al. Repurposing the aldose reductase inhibitor and diabetic neuropathy drug epalrestat for PMM2-CDG. Dis Model Mech, 2019
in patient fibroblasts across multiple PMM2 genotypes, and is the only small-molecule activator of PMM2 identified to date.

Interactions

PMM2-CDG is defined by compound heterozygosity: one R141H allele paired with a second, distinct PMM2 pathogenic variant on the other chromosome. The most common compound heterozygous combinations include R141H + F119L and R141H + V44A, with phenotypic severity correlating broadly with residual enzyme activity — milder second alleles yield milder disease. Any individual carrying R141H who is diagnosed with PMM2-CDG should have full sequencing of both PMM2 alleles to characterize the specific compound genotype, as this informs prognosis and eligibility for emerging genotype-specific treatments.

rs35936514

LHPP LHPP depression risk variant

Strong Risk Factor

The Prefrontal Phosphatase — When Stress Rewires the Brain

Most people associate depression with serotonin. LHPP tells a different story — one about a small enzyme in the prefrontal cortex that quietly prevents stress from reshaping your brain. When LHPP function is reduced, chronic stress finds an open door into the neural circuitry that governs mood, motivation, and resilience.

LHPP encodes phospholysine phosphohistidine inorganic pyrophosphate phosphatase11 phospholysine phosphohistidine inorganic pyrophosphate phosphatase
A HAD-family phosphatase that removes phosphate groups from histidine residues on target proteins. Histidine phosphorylation is an understudied but functionally important post-translational modification in mammalian cells, regulating kinase signaling cascades
. It was essentially unknown in psychiatry until the CONVERGE consortium22 CONVERGE consortium
Collaborative Research on the Genetics of Depression in East Asians — a large-scale genetics consortium focused on recurrent major depressive disorder in Han Chinese women
identified it in a landmark 2015 Nature study33 a landmark 2015 Nature study
CONVERGE consortium. Sparse whole-genome sequencing identifies two loci for major depressive disorder. Nature, 2015
, sequencing 5,303 cases and 5,337 controls to find just two genome-wide significant loci — one near SIRT1 and one in LHPP (P = 6.45 × 10⁻¹²). The rs35936514 variant sits in the 3' UTR and flanking intronic region of LHPP on chromosome 10, likely affecting expression levels rather than protein structure.

The Mechanism

LHPP acts as a brake on stress-induced neuronal dysfunction. In the medial prefrontal cortex44 medial prefrontal cortex
A region of the prefrontal cortex critical for executive function, emotional regulation, and the top-down suppression of the stress response. It communicates extensively with the amygdala and hippocampus to modulate fear and mood
(mPFC), LHPP dephosphorylates a network of histidine kinase substrates55 histidine kinase substrates
Proteins — including NME1/2 — that carry phosphate groups on histidine amino acids. When these kinases are over-active, downstream signaling cascades including MAPK and PI3K/AKT are altered, reducing glutamatergic synaptic transmission
including NME1/2. When LHPP levels fall — whether from reduced expression due to the T allele or from chronic stress itself suppressing LHPP — histidine phosphorylation rises unchecked, weakening glutamatergic synaptic transmission in exactly the circuits that should be most resilient.

Lin et al. (2023)66 Lin et al. (2023)
Lin D et al. LHPP, a risk factor for major depressive disorder, regulates stress-induced depression-like behaviors through its histidine phosphatase activity. Molecular Psychiatry, 2023
showed in mice that mPFC LHPP deletion alone produced no spontaneous depression-like behavior — but on exposure to chronic social defeat stress, LHPP-knockout mice showed dramatically augmented behavioral deficits. Crucially, re-introducing wild-type LHPP reversed the deficits, while a phosphatase-dead LHPP mutant did not, confirming the enzymatic activity is required. A parallel mechanism runs through astrocytes: Sha et al. (2023)77 Sha et al. (2023)
Sha L et al. LHPP-mediated inorganic pyrophosphate hydrolysis-driven lysosomal acidification in astrocytes regulates adult neurogenesis. Cell Reports, 2023
showed that LHPP normally drives lysosomal acidification88 lysosomal acidification
Lysosomes require an acidic environment (pH ~5) to function as the cell's recycling centers. In astrocytes, LHPP hydrolyzes inorganic pyrophosphate to fuel this acidification. Without it, lysosomal function is impaired, the C/EBPβ transcription factor accumulates, and downstream chemokine signaling promoting neurogenesis is triggered
in astrocytes; in its absence, the C/EBPβ pathway triggers compensatory adult neurogenesis that confers stress resilience. This counterintuitive finding (less LHPP → more resilience in knockout mice) may reflect that the full-knockout model removes a gene globally, while the rs35936514 T allele reduces rather than abolishes LHPP expression in specific circuits.

In the ventral hippocampus, a complementary mechanism operates: Zhuang et al. (2024)99 Zhuang et al. (2024)
Zhuang L et al. LHPP in glutamatergic neurons of the ventral hippocampus mediates depression-like behavior by dephosphorylating CaMKIIα and ERK. Biological Psychiatry, 2024
found that LHPP levels increase in glutamatergic neurons during stress-induced depression, and that knocking out LHPP in these neurons enhanced spontaneous activity and stress resilience via the CaMKIIα/ERK pathway1010 CaMKIIα/ERK pathway
CaMKII (calcium/calmodulin- dependent protein kinase II) and ERK (extracellular signal-regulated kinase) regulate synaptic plasticity and BDNF expression. LHPP dephosphorylates both, reducing their activity and impairing the BDNF/PSD95-mediated synaptic strengthening that underlies mood resilience
. Esketamine — but not fluoxetine — reversed LHPP-induced depression-like behaviors, pointing toward glutamatergic rather than serotonergic modulation as the relevant therapeutic axis.

The Evidence

The CONVERGE study represents one of the cleanest GWAS findings in depression genetics. By restricting the sample to Han Chinese women with recurrent, severe MDD — reducing phenotypic and genetic heterogeneity — the consortium achieved genome-wide significance with just 10,640 individuals, far fewer than most psychiatric GWAS require. The LHPP signal (P = 6.45 × 10⁻¹²) was validated in an independent replication cohort. rs35936514 sits in intron 2 of LHPP and the flanking 3' UTR, and the T allele is associated with reduced LHPP expression in the prefrontal cortex.

At the neural circuit level, Cui et al. (2020)1111 Cui et al. (2020)
Cui L et al. Association of LHPP genetic variation (rs35936514) with structural and functional connectivity of hippocampal-corticolimbic neural circuitry. Brain Imaging and Behavior, 2020
examined 122 healthy participants and found that T-allele carriers showed increased hippocampal connectivity to the rostral anterior cingulate cortex and higher fractional anisotropy in the fornix — structural and functional signatures consistent with altered top-down emotional regulation. These brain connectivity differences were present without clinical depression, suggesting the risk phenotype is a continuous trait rather than a categorical diagnosis.

Earlier family-based work by Neff et al. (2009)1212 Neff et al. (2009)
Neff CD et al. Evidence for HTR1A and LHPP as interacting genetic risk factors in major depression. Molecular Psychiatry, 2009
found a chromosome 10 linkage peak (HLOD = 4.4) in Utah pedigrees with familial depression and identified disease-segregating LHPP SNPs in the linkage region. The LHPP effects in that study were contingent on HTR1A genotype, hinting at an interaction between serotonin receptor sensitivity and LHPP-mediated stress signaling.

Practical Implications

The LHPP mechanism converges on a clear biological theme: glutamatergic synaptic function under stress. T-allele carriers do not appear to have impaired mood at baseline — the risk emerges under chronic stress. This has direct implications for how carriers should think about their mental health strategy: minimizing the duration and intensity of chronic stressors matters more than managing acute stress reactions.

The esketamine finding from Zhuang et al. is clinically relevant: if depression does develop in T-allele carriers, the evidence points toward glutamate-targeting treatments (esketamine/ketamine) having particular mechanistic relevance versus standard SSRIs, which do not directly address the LHPP-CaMKII/ERK pathway.

Interactions

The Neff et al. (2009) data suggest a potential interaction between LHPP (rs35936514) and the HTR1A serotonin receptor gene. In that family-based study, the LHPP depression association was modulated by HTR1A genotype, suggesting that individuals carrying risk variants in both genes may have compounded vulnerability. The mechanistic intersection — serotonergic regulation of prefrontal excitatory tone — offers a plausible biological substrate.

The BDNF Val66Met variant (rs6265) is also relevant: since LHPP deficiency reduces CaMKII/ERK activity and thereby impairs BDNF/PSD95-mediated synaptic strengthening, T-allele carriers who also carry the BDNF Met allele (reduced activity-dependent BDNF secretion) may experience compounded vulnerability to stress-induced synaptic weakening in prefrontal and hippocampal circuits.

rs3850641

TNFSF4 TNFSF4 (OX40L) Intron 1 Variant

Moderate Risk Factor

OX40 Ligand and the Immune Gate of Atherosclerosis

TNFSF4 encodes OX40 ligand (OX40L, also called CD252), a co-stimulatory molecule of the TNF superfamily that serves as the key amplifier of T-cell activation in inflammatory settings. OX40L is expressed on antigen-presenting cells — particularly macrophages, dendritic cells, and vascular endothelial cells — and engages OX40 on activated CD4+ T cells to extend their survival, promote cytokine secretion, and drive clonal expansion. Within atherosclerotic plaques, OX40L expression is markedly elevated compared to healthy vascular tissue, driven by oxidized LDL and other inflammatory mediators. The rs3850641 variant sits in the first intron of TNFSF4, in strong linkage disequilibrium with a promoter polymorphism (rs45454293) that regulates OX40L gene expression.

The Mechanism

The functional variant linked to rs3850641 is rs4545429311 rs45454293
A promoter SNP in LD with rs3850641 at r²=0.70; the T-allele reduces TNFSF4 transcription by ~50% compared to the C-allele in luciferase reporter assays, attributed to binding of a transcriptional repressor
. The rs3850641 SNP itself has no direct effect on promoter activity, but because both SNPs travel together on the same haplotype block — the TG haplotype combining rs45454293T and rs3850641G — the G allele at rs3850641 tags a state of reduced OX40L expression in vascular tissue.

Paradoxically, lower OX40L expression appears to increase cardiovascular risk. This may reflect an immunoregulatory role: OX40L signaling is also required for the proper differentiation and persistence of regulatory T cells (Tregs). Disruption of OX40L in the plaque microenvironment could impair Treg-mediated suppression of inflammatory effector T cells, allowing uncontrolled T-cell–driven inflammation to accelerate plaque growth and destabilization. Consistent with this, anti-OX40L antibody blockade in LDL receptor-deficient mice22 anti-OX40L antibody blockade in LDL receptor-deficient mice
Guo et al., Arterioscler Thromb Vasc Biol 2007 — anti-OX40L antibody reduced atherosclerotic lesion area by approximately 53% in Ldlr−/− mice by attenuating Th2-mediated responses and preserving anti-oxLDL IgM antibodies
attenuated atherogenesis, pointing to the pathway as a therapeutic target rather than a simple risk amplifier.

Elevated soluble OX40L (sOX40L) in circulation — released primarily from activated platelets — correlates with acute coronary syndrome severity, plaque instability, and elevated matrix metalloproteinase (MMP-9, MMP-3) levels, reinforcing OX40L as a biomarker of active plaque inflammation.

The Evidence

The variant was first implicated by Wang et al. in Nature Genetics 200533 Wang et al. in Nature Genetics 2005
Positional identification of TNFSF4, encoding OX40 ligand, as a gene that influences atherosclerosis susceptibility
, which used positional cloning from the murine atherosclerosis susceptibility locus Ath1 to identify TNFSF4 as a human candidate gene. The rs3850641 G allele was significantly more frequent in MI cases than controls in two independent Swedish cohorts. Subsequent work by Cunningham et al. PLoS ONE 201144 Cunningham et al. PLoS ONE 2011
A Common Polymorphism in the Promoter Region of the TNFSF4 Gene Is Associated with Lower Allele-Specific Expression and Risk of MI
identified rs45454293 as the functional promoter variant in LD with rs3850641 (r²=0.70), with the TG haplotype significantly enriched in female MI patients vs controls in the SCARF cohort (P=0.01), and demonstrated 50% reduced OX40L promoter activity for the T-allele.

A Chinese Han case-control study (Li et al. Oncotarget 2018, PMID 29921578 — 454 cases, 512 controls)55 (Li et al. Oncotarget 2018, PMID 29921578 — 454 cases, 512 controls) found that homozygous GG carriers had significantly elevated MI risk (crude OR 2.00, 95% CI 1.04–3.86, P=0.039; adjusted OR 2.29, 95% CI 1.20–4.69, P=0.023), while the AG heterozygous genotype conferred no significant independent risk under additive or dominant models.

However, larger meta-analyses show overall null results in aggregate populations. A meta-analysis of 9 studies (Fu et al. 2016, PMID 27008001)66 meta-analysis of 9 studies (Fu et al. 2016, PMID 27008001) found allelic OR 1.10 (95% CI 0.96–1.27, P=0.174), and a systematic review of 11 studies covering 3,865 cases and 6,344 controls (Lu et al. 2018, PMID 29424751)77 systematic review of 11 studies covering 3,865 cases and 6,344 controls (Lu et al. 2018, PMID 29424751) found OR 1.02 (95% CI 0.89–1.17) for the G allele — no overall association. These null findings likely reflect population stratification (original signal primarily in European women), heterogeneity across disease definitions (CHD vs. MI specifically), and statistical power limitations for the rare GG genotype (only ~2% of individuals).

The GG-specific risk signal (recessive model OR ~2) remains coherent across the positive studies and is mechanistically plausible; the overall null meta-analysis findings reflect the rarity of GG homozygotes diluting allele-level effects. Evidence is classified as moderate — replicated in multiple independent studies under the recessive model, with a plausible regulatory mechanism, but no established clinical guidelines and significant heterogeneity across populations.

Practical Actions

For GG homozygotes, the elevated MI risk under the recessive model (OR ~2) warrants proactive cardiovascular monitoring. The risk appears particularly pronounced in women and in populations of East Asian and South Asian descent where OX40L pathway activity in the context of this haplotype may interact differently with other risk factors. Tracking inflammatory biomarkers such as hs-CRP and sOX40L alongside standard lipid panels provides the most relevant picture of vascular inflammatory burden for this genotype.

Interactions

rs3850641 is in linkage disequilibrium with the functional promoter SNP rs45454293 (r²=0.70); the TG haplotype combining both minor alleles carries the highest risk signal in women. rs2205960 and other upstream TNFSF4 promoter region SNPs (rs10912580, rs12039904, rs1234317) form a distinct haplotype block associated with systemic lupus erythematosus susceptibility, illustrating that TNFSF4 haplotypes have different disease consequences depending on which part of the regulatory region they tag.

Individuals carrying both rs3850641 GG and cardiovascular-risk variants in classical inflammatory pathways (e.g., IL-6 signaling genes, CRP variants) may have compounded atherosclerotic risk through synergistic immune activation in the vascular wall — an interaction that has not been formally studied but is mechanistically plausible given the T-cell co-stimulatory nature of the OX40L pathway.