PPARA (Peroxisome Proliferator-Activated Receptor Alpha) is a nuclear receptor that acts as a master regulator of fatty acid oxidation, ketogenesis, and energy homeostasis. It controls how efficiently your muscles burn fat for fuel during prolonged exercise. The intron 7 G/C variant (rs4253778) alters PPARA expression by changing transcription factor binding sites¹¹ The intronic SNP changes binding motifs for the interferon regulatory factor (IRF) family of transcription factors, affecting how much PPARA protein is produced, which cascades into differences in muscle fiber composition, fuel utilization during exercise, and how the heart adapts to physical training.

The Mechanism

The G allele preserves normal PPARA expression, promoting efficient fatty acid oxidation in skeletal muscle and the heart. This drives a higher proportion of type I slow-twitch muscle fibers²² type I slow-twitch muscle fibers
Slow-twitch fibers are oxygen-efficient and fatigue-resistant, ideal for endurance activities like distance running and cycling, which are optimized for sustained aerobic activity. The C allele reduces PPARA function, shifting muscle metabolism away from fat oxidation toward glucose utilization. This favors type II fast-twitch fibers³³ type II fast-twitch fibers
Fast-twitch fibers generate rapid, powerful contractions but fatigue quickly, suited to sprinting and power sports and greater muscle hypertrophy, including of the heart itself.

Critically, this variant also affects cardiac remodeling. Reduced PPARA activity in C allele carriers means the heart relies more on glucose for energy, which promotes greater left ventricular growth in response to exercise or elevated blood pressure.

The Evidence

The landmark Jamshidi et al. study⁴⁴ landmark Jamshidi et al. study
Jamshidi Y et al. Peroxisome proliferator-activated receptor alpha gene regulates left ventricular growth in response to exercise and hypertension. Circulation, 2002 followed 144 British Army recruits through 10 weeks of physical training and found that the effect of the C allele on left ventricular mass was additive: GC heterozygotes gained 11.8 g of left ventricular mass (versus 6.7 g for GG), while CC homozygotes gained 19.4 g -- a nearly 3-fold greater increase. In a separate hypertension cohort (n=1,148), the C allele was also associated with greater left ventricular hypertrophy.

Ahmetov et al.⁵⁵ Ahmetov et al.
Ahmetov II et al. PPARalpha gene variation and physical performance in Russian athletes. Eur J Appl Physiol, 2006 studied 786 Russian athletes and 1,242 controls, finding the GG genotype significantly overrepresented in endurance athletes (80.3% vs 70.0% in controls, P=0.0001). Muscle biopsies confirmed GG homozygotes had a higher percentage of slow-twitch fibers (55.5% vs 38.5%, P=0.003). An increasing linear trend of C allele frequency was observed with increasing anaerobic component of performance (P=0.029).

A meta-analysis of five studies⁶⁶ meta-analysis of five studies
Lopez-Leon S et al. Sports genetics: the PPARA gene and athletes' high ability in endurance sports. A systematic review and meta-analysis. Biol Sport, 2016 pooling 760 endurance athletes and 1,792 controls confirmed the association: the G allele had an OR of 1.65 (95% CI 1.39--1.96) for endurance ability, with no heterogeneity (I^2=0%) or publication bias.

Conversely, Petr et al.⁷⁷ Petr et al.
Petr M et al. PPARA intron polymorphism associated with power performance in 30-s anaerobic Wingate Test. PLoS ONE, 2014 showed that C allele carriers among Czech ice hockey players achieved significantly higher anaerobic peak power (14.6 vs 13.9 W/kg, P=0.036), supporting the C allele's role in power-oriented performance.

A training response study⁸⁸ training response study
Leońska-Duniec A et al. The polymorphisms of the peroxisome-proliferator activated receptors' alfa gene modify the aerobic training induced changes of cholesterol and glucose. J Clin Med, 2019 in 168 women found CC homozygotes had unfavorable metabolic responses to 12 weeks of aerobic training: LDL cholesterol increased (79 to 95 mg/dL) and glucose rose (70.5 to 78.2 mg/dL), while GG and GC carriers showed beneficial decreases.

Practical Implications

Your PPARA intron 7 genotype helps explain which type of physical activity suits your body best. GG carriers are genetically predisposed to excel in endurance sports and benefit from aerobic exercise through efficient fat burning. CG carriers have intermediate characteristics, maintaining some endurance capacity while gaining more from strength-oriented training. CC carriers are oriented toward power and strength, but should pay attention to their cardiovascular and metabolic response to exercise -- particularly monitoring LDL cholesterol and ensuring they include adequate aerobic conditioning.

The cardiac hypertrophy finding is important for any C allele carrier who trains intensely: it represents an exaggerated but physiological adaptation, not a disease state, but it may warrant echocardiographic monitoring for athletes in high-volume training programs.

Interactions

PPARA rs4253778 interacts with the L162V variant (rs1800206) in the same gene. The L162V variant alters the DNA-binding domain of the PPARA protein and affects lipid metabolism independently. Carrying unfavorable alleles at both positions may compound the impact on cholesterol response to exercise. The C allele at rs4253778 combined with the Val162 allele at rs1800206 has been associated with more pronounced adverse lipid changes during training.

PPARA also interacts functionally with ACTN3 (rs1815739), which independently influences muscle fiber composition. The ACTN3 XX genotype (alpha-actinin-3 deficiency) combined with the PPARA GG genotype would strongly favor endurance, while ACTN3 RR with PPARA CC would favor power.

Toll-like receptor 9 (TLR9)¹¹ Toll-like receptor 9 (TLR9)
An endosomal pattern recognition receptor that detects unmethylated CpG dinucleotide motifs in bacterial and viral DNA, triggering innate immune activation through MyD88 and downstream NF-κB and interferon regulatory factor pathways is one of the body's most fundamental alarm systems for microbial invasion. The rs5743836 variant — the third member of the classical TLR9 promoter haplotype trio alongside [rs187084 (-1486T/C) | Located 249 base pairs further upstream of the TLR9 transcription start site; also shows estrogen-responsive regulation and is independently associated with HCV clearance, SLE susceptibility, and osteoarthritis risk] and [rs352140 (exon 2 synonymous) | The most-studied TLR9 variant, altering mRNA stability and TLR9 expression levels without changing the protein sequence] — is positioned 1,237 base pairs upstream of the TLR9 coding sequence in the core promoter region.

In genome files (WGS and consumer chip arrays), this variant is reported as A or G on the plus strand of chromosome 3. Published papers describe it as -1237T/C in coding (minus-strand) notation, because TLR9 lies on the minus strand. The correspondence is: the paper's "T" allele is A on the plus strand (the reference, found in ~79% of people globally), and the paper's "C" allele is G on the plus strand (the alternate at ~21% globally). What makes rs5743836 mechanistically distinct from its companion promoter variant rs187084 is not just its position — it is the nature of the transcription factor binding site it creates.

The Mechanism

The [C allele (plus-strand G) at rs5743836 | Also referred to as -1237C in coding-strand notation throughout the literature] introduces a new regulatory motif at -1237 in the TLR9 promoter: an IL-6-responsive element (IL-6RE)²² IL-6-responsive element (IL-6RE)
A DNA sequence recognized by STAT3, the transcription factor activated downstream of the IL-6 receptor; STAT3 binding to this element drives TLR9 transcription in response to IL-6 signaling. This creates a positive feedback circuit: when TLR9 detects CpG DNA from invading pathogens, it triggers cytokine production including IL-6, which then binds its receptor on immune cells, activates STAT3, and STAT3 binds the newly created IL-6RE to drive further TLR9 transcription. The TT (plus-strand AA) reference genotype lacks this element; the TC (plus-strand AG) and CC (plus-strand GG) genotypes carry it.

[Functional experiments by Carvalho et al. | PLoS One 2011, PMID 22132241] demonstrated this loop directly: TC genotype carriers show higher TLR9 expression when treated with IL-6, and blocking either IL-6 signaling or TLR9 itself reversed the enhanced B-cell proliferation observed in TC individuals upon CpG stimulation. The promoter variant also exhibits [estrogen-responsive regulation | Fischer et al. (Gut 2017, PMID 27196570) showed that both rs5743836 and rs187084 C alleles display estrogen receptor-dependent transcriptional activity, explaining why women carrying these alleles have greater innate immune responses to some pathogens], paralleling the estrogen-dependent effects of the companion rs187084 variant. Together, these features make rs5743836 not merely a quantitative amplifier of TLR9 expression but a qualitative switch that wires TLR9 regulation to two major inflammatory circuits simultaneously: IL-6/STAT3 and estrogen signaling.

The Evidence

The clinical consequences of this molecular wiring are most clearly documented in lymphoma biology, malaria immunity, and sex-specific thrombosis risk.

Non-Hodgkin and Hodgkin lymphoma: The strongest and most replicated association for rs5743836 is with lymphoma. Carvalho et al. (Genes Immun 2012) studied three independent European cohorts totaling over 4,700 subjects: the C allele (plus-strand G) increased NHL risk with OR=1.85 in Portugal (p=7.3×10⁻⁹) and OR=1.84 in Italy (p=6.0×10⁻⁵), though no significant association was seen in the US cohort³³ Carvalho et al. (Genes Immun 2012) studied three independent European cohorts totaling over 4,700 subjects: the C allele (plus-strand G) increased NHL risk with OR=1.85 in Portugal (p=7.3×10⁻⁹) and OR=1.84 in Italy (p=6.0×10⁻⁵), though no significant association was seen in the US cohort
The European-specific replication suggests population-stratified effects, possibly linked to different pathogen exposures and haplotype backgrounds. The mechanistic link is the IL-6/STAT3 pathway: IL-6-driven TLR9 upregulation enhances CpG-induced B-cell proliferation in C allele carriers, and uncontrolled B-cell proliferation is a hallmark of many NHL subtypes. This creates a plausible chain from variant → molecular phenotype → disease.

[A 2022 case-control study in 136 Jordanian Hodgkin lymphoma patients and 238 controls found significantly higher rs5743836 variant allele frequency in cases (p=0.031), with significance across codominant, dominant, and overdominant models | Al-Khatib et al., PLoS One 2022 (PMID 35905120)]. This was synthesized in [a 2025 meta-analysis across multiple cohorts finding OR=1.54 (95% CI 1.03–2.32, p=0.036) in the dominant model for lymphoma overall | Yan et al., BMC Cancer 2025 (PMID 40169945)], confirming the lymphoma association as the most consistent clinical signal for this variant.

Malaria: rs5743836 shows a clear genotype-dose relationship for malaria susceptibility, in the opposite direction from lymphoma. [Omar et al. (Malar J 2012) followed 429 Ghanaian children for one year and found parasitemia levels strongly correlated with rs5743836 genotype: CC (plus-strand GG) carriers had mean parasitemia of 23,532/µL, TC (AG) carriers 14,924/µL, and TT (AA) carriers the lowest at 5,501/µL (p=0.03) | Haplotype analysis showed the TTAG four-SNP haplotype was associated with relative risk of 0.2 for symptomatic malaria (PMID 22594374)]. Independently, [Esposito et al. (Malar J 2012) found CC genotype (plus-strand GG) significantly associated with increased malaria risk in 602 Burundian children (p=0.03) | This was the only TLR9 variant significant for malaria susceptibility in that cohort across 939 subjects (PMID 22691414)].

The paradox — the same allele that drives lymphoma risk also increases malaria susceptibility — resolves when you consider the immunological context. In malaria, the critical host defense is early cytokine-driven clearance of erythrocyte-infected parasites and activation of protective T helper cell responses. Excess IL-6/STAT3-mediated TLR9 upregulation may dysregulate the balance between protective Th1 responses and regulatory immune suppression, paradoxically impairing parasite clearance in the high-IL-6 malaria environment.

Plasmodium vivax vaccine responses: [Carrión-Nessi et al. (PLoS Negl Trop Dis 2025) examined IgG responses against P. vivax circumsporozoite protein (PvCSP) variants in 210 Venezuelan patients. TC (plus-strand AG) heterozygotes produced reduced antibody responses: adjusted OR=0.26 for PvCSP VK247 and OR=0.37 for PvCSP V-like (PMID 40587574)]. High antibody responders had significantly fewer symptoms (p<0.001), implying that genotype-driven antibody differences have clinical consequences for malaria disease burden — and that future P. vivax vaccines may need to account for TLR9 genotype in immunogenicity predictions.

Venous thromboembolism recurrence (sex-specific): [Ahmad et al. (J Thromb Thrombolysis 2017) analyzed 1,050 VTE patients in the Malmö thrombophilia study. In women, rs5743836 was significantly associated with VTE recurrence (HR=3.46, 95% CI 1.06–11.33), rising to HR=5.94 (95% CI 1.25–28.13) for unprovoked VTE. No association was found in men (PMID 28321710)]. The sex specificity aligns with the estrogen-responsive nature of this promoter variant: in women, estrogen-enhanced TLR9 expression may promote thrombo-inflammatory cascades involving platelet TLR9 activation and immune-mediated coagulation dysregulation.

SLE: Despite being a TLR9 promoter variant, rs5743836 does not appear to influence SLE susceptibility. Two large meta-analyses (Wang et al. 2016; Lee & Song 2023; Hu et al. 2017) covering over 5,000 SLE cases found no significant association in any genetic model or ancestry group. The SLE signal in the TLR9 locus appears to be carried primarily by rs187084 in Asian populations, not by rs5743836.

Practical Implications

The G allele's functional impact — creating an IL-6/STAT3 feedback loop in the TLR9 promoter — translates into a heightened and self-amplifying innate immune response when CpG DNA is encountered. The clearest clinical consequences are elevated lymphoma susceptibility (particularly for B-cell lymphomas, where the IL-6-TLR9-B-cell proliferation axis is most directly relevant), increased malaria susceptibility, and female-specific VTE recurrence risk (likely through estrogen-responsive TLR9 amplification of thrombo-inflammatory pathways). There are no established dietary or supplement interventions known to modulate TLR9 promoter activity. The actionable steps for carriers of the G allele center on hematologic vigilance, malaria precautions, and — in women — awareness of the VTE recurrence risk when evaluating anticoagulation decisions.

Interactions

The rs5743836 promoter variant is almost always studied as part of the [TLR9 three-SNP haplotype | The canonical haplotype comprises rs187084 (-1486T/C), rs5743836 (-1237T/C), and rs352140 (exon 2 synonymous) — all three modulate TLR9 expression through different mechanisms and are frequently co-inherited]. The [rs187084 (-1486T/C) | Companion promoter variant 249 base pairs further upstream; shows estrogen-responsive regulation and is significantly associated with SLE in Asians, HCV clearance in women, OA risk, and post-bronchiolitis wheezing] shares the estrogen-responsive regulation with rs5743836, suggesting their combined promoter haplotype may have additive effects on TLR9 expression in women. [Fischer et al. (2017, PMID 27196570) found that both rs187084 and rs5743836 C alleles show estrogen receptor-dependent allele-specific mRNA regulation | The two promoter variants may act synergistically to maximize estrogen-driven TLR9 expression in women]. Individuals carrying C/G alleles at both promoter positions may have the strongest innate immune advantage for viral clearance and the highest liability for TLR9-driven lymphoproliferation.

The pathway context also connects to [TLR4 (rs4986790) | Detects bacterial LPS; a complementary innate immune pathway studied alongside TLR9 in the same Leishmania infantum cohort without significant independent effects] and [TLR2 (rs5743708) | Recognizes bacterial lipoproteins and peptidoglycan; another MyD88-dependent TLR complementary to CpG DNA sensing].

Apolipoprotein L1 (APOL1) is a secreted protein that circulates on HDL particles¹¹ secreted protein that circulates on HDL particles
APOL1 is the only human apolipoprotein with trypanolytic activity — it kills African trypanosomes by forming ion channels in their lysosomal membranes and serves as the innate immune system's weapon against Trypanosoma brucei, the parasite causing African sleeping sickness. The G1 risk haplotype consists of two missense variants in near-perfect linkage disequilibrium: rs73885319 (S342G) and rs60910145 (I384M). Both sit within the SRA-interacting domain²² SRA-interacting domain
Serum Resistance-Associated protein (SRA) is expressed by the human-infective subspecies T.b. rhodesiense; it binds and neutralizes wild-type APOL1 of the APOL1 protein, and together they alter the protein's ability to be neutralized by the parasite — conferring resistance to sleeping sickness but, under a recessive model, dramatically increasing kidney disease risk.

The Mechanism

The I384M substitution (isoleucine to methionine at position 384) occurs in the C-terminal region of the SRA-binding domain. While functional studies show that the S342G component drives the trypanolytic gain-of-function³³ functional studies show that the S342G component drives the trypanolytic gain-of-function
Cooper et al. 2017 demonstrated that S342G alone confers trypanosome resistance in vivo, while I384M alone does not, both variants are inherited together as the G1 haplotype and are required for the full kidney disease risk phenotype. The APOL1 risk variants cause kidney injury through multiple mechanisms: they form active cation channels at the plasma membrane⁴⁴ active cation channels at the plasma membrane
Risk-variant APOL1 inserts into podocyte membranes, creating ion pores that disrupt cellular homeostasis of kidney podocytes, induce mitochondrial dysfunction, and trigger endoplasmic reticulum stress. Crucially, disease requires a recessive model — two risk alleles⁵⁵ recessive model — two risk alleles
G1/G1, G2/G2, or G1/G2 compound heterozygosity; carriers of a single risk allele have minimal kidney risk — meaning one copy is protective against trypanosomes without causing kidney harm, while two copies cross a threshold into cytotoxicity.

The Evidence

The landmark 2010 discovery⁶⁶ landmark 2010 discovery
Genovese et al. Association of trypanolytic ApoL1 variants with kidney disease in African Americans. Science, 2010 identified the G1 and G2 haplotypes as the genetic explanation for the 3- to 5-fold excess kidney disease burden in African Americans. Two-risk-allele carriers face dramatically elevated odds: OR 17 for FSGS, OR 29 for HIV-associated nephropathy (HIVAN), and OR 7–10 for hypertension-attributed end-stage renal disease⁷⁷ OR 17 for FSGS, OR 29 for HIV-associated nephropathy (HIVAN), and OR 7–10 for hypertension-attributed end-stage renal disease
Kopp et al. JASN, 2011. The AASK and CRIC cohort study⁸⁸ AASK and CRIC cohort study
Parsa et al. APOL1 risk variants, race, and progression of chronic kidney disease. NEJM, 2013 confirmed that the high-risk genotype independently accelerates CKD progression (HR 1.88 for composite renal endpoints), regardless of baseline kidney function.

However, penetrance is incomplete — approximately 15–20% of two-risk-allele carriers develop clinical kidney disease, indicating that "second hits" are required⁹⁹ "second hits" are required
Known triggers include HIV infection, interferon signalling (e.g. from COVID-19 or lupus), and hypertension — each upregulates APOL1 expression in podocytes. This means the genotype creates susceptibility, not certainty. A targeted therapy, inaxaplin, reduced proteinuria by 47.6%¹⁰¹⁰ inaxaplin, reduced proteinuria by 47.6%
Phase 2a trial: 13 weeks of inaxaplin in patients with two APOL1 risk alleles and biopsy-proven FSGS in a phase 2a trial of APOL1-associated FSGS, marking the first genotype-directed kidney disease treatment.

Population Context

The G1 haplotype is found almost exclusively in people of African ancestry — approximately 22% allele frequency in African and African American populations, compared to less than 0.01% in Europeans and East Asians. This extreme population stratification reflects positive selection driven by trypanosome resistance¹¹¹¹ positive selection driven by trypanosome resistance
The G1 and G2 variants rose to high frequency in West Africa because heterozygous carriers were protected against T.b. rhodesiense sleeping sickness — a classic example of balancing selection similar to sickle cell trait and malaria. Among African Americans, approximately 13% carry two APOL1 risk alleles (G1/G1, G2/G2, or G1/G2), placing them in the high-risk category.

Practical Implications

For carriers of a single G1 allele (GT genotype), the clinical consequence is minimal for kidney health and may offer innate immunity benefits against trypanosome infection. For two-risk-allele carriers (requiring knowledge of both APOL1 risk variant loci — G1 and G2), the imperative is early and sustained kidney monitoring to detect subclinical disease before irreversible nephron loss. Blood pressure control is critical because hypertension acts as a "second hit" that accelerates APOL1-mediated podocyte injury.

Interactions

Rs60910145 is in near-perfect linkage disequilibrium with rs73885319 — together they define the G1 haplotype. The kidney disease risk requires two APOL1 risk alleles in any combination: G1/G1 homozygosity, G2/G2 homozygosity (rs71785313), or G1/G2 compound heterozygosity. A compound action covering the G1+G2 interaction would be appropriate for users carrying risk alleles at both loci. The recently discovered N264K modifier variant (rs73015316) substantially reduces penetrance when co-inherited with G2-containing haplotypes — an important consideration for future genotyping completeness.

Every cell in your body needs vitamin C, but not all cells are equal in their ability to acquire it. Once dietary ascorbate crosses the gut wall and enters the bloodstream, a second transporter system distributes it into tissues where it is needed most — the brain, adrenal glands, eyes, and metabolically active organs. The gene SLC23A2 encodes SVCT2¹¹ SVCT2
Sodium-dependent Vitamin C Transporter 2 — a high-affinity transporter expressed in metabolically demanding tissues including neurons, adrenal cortex cells, and the aqueous humor of the eye, a high-affinity transporter that pulls ascorbate from the bloodstream into these specialized tissues. The intronic variant rs6133175 sits within SLC23A2 and, despite lying outside the protein-coding sequence, influences circulating plasma vitamin C in a measurable way — with the GG genotype associated with approximately 24% higher plasma ascorbate than the common AA genotype.

The Mechanism

Unlike SLC23A1 (SVCT1), which handles intestinal absorption and renal reabsorption to maintain whole-body vitamin C homeostasis, SVCT2 operates in tissues with high metabolic demand. The brain accumulates vitamin C to concentrations roughly 10-fold higher than plasma — a feat achieved almost entirely by SVCT2 expressed on the blood-brain barrier and neuronal membranes. The adrenal glands similarly use SVCT2 to build the highest vitamin C concentration of any organ in the body, where ascorbate is required for cortisol and adrenaline synthesis.

The rs6133175 variant is an intron variant²² intron variant
Located within a non-coding intervening sequence of the gene; does not directly change the amino acid sequence but may affect gene expression, splicing, or regulatory element activity — its precise molecular mechanism has not been fully characterized. The most likely explanations are altered splicing efficiency³³ splicing efficiency
The process by which intron sequences are removed from pre-mRNA; intronic variants near splice sites can shift the ratio of splice isoforms, changing how much functional protein is made or disruption of a transcriptional regulatory element within the intron. The fact that it influences plasma vitamin C at all — despite not changing the SVCT2 protein sequence — suggests it affects the amount of transporter protein expressed rather than its function.

Because SVCT2 handles redistribution of ascorbate from plasma into tissues, a variant that increases SVCT2 expression or activity could lower plasma levels by pulling more vitamin C into cells, or raise plasma levels by improving renal reabsorption secondary effects. The net effect observed in the EPIC cohort is higher plasma vitamin C in GG homozygotes, though the direction of causality through tissue distribution remains to be mechanistically confirmed.

The Evidence

The primary evidence comes from a nested case-control study in the European EPIC cohort⁴⁴ nested case-control study in the European EPIC cohort
Duell EJ et al. Vitamin C transporter gene (SLC23A1 and SLC23A2) polymorphisms, plasma vitamin C levels, and gastric cancer risk in the EPIC cohort. Genes Nutr, 2013 involving 365 gastric cancer cases and 1,284 matched controls from 10 European countries. Among 311 controls with complete genotyping and plasma vitamin C data, genotype frequencies were AA 37%, AG 45%, and GG 12%. In a recessive model adjusted for age, sex, country, smoking, and season of blood draw, GG homozygotes had plasma vitamin C 24% higher than AA homozygotes (beta = 0.22, 95% CI: 0.029–0.40; P = 0.02). The raw plasma values were AA 39.1, AG 39.4, and GG 45.2 umol/L — a clinically meaningful spread given that adequate status is generally considered to be above 28 umol/L.

Importantly, both rs6133175 (SLC23A2) and rs33972313 (SLC23A1) independently predicted plasma vitamin C levels in multivariable models, suggesting the two genes tag non-overlapping mechanisms: SVCT1 controls gut absorption and renal reclamation, while SVCT2 variants apparently influence a separate step in vitamin C homeostasis.

A Chinese Han population study⁵⁵ Chinese Han population study
Hou H et al. Impact of SLC23A1 and SLC23A2 Polymorphisms on the Risk for Preeclampsia in a Chinese Han Population. J Nutr Sci Vitaminol (Tokyo), 2022 found significant genotypic frequency differences for rs6133175 between preeclampsia cases and controls. Under a recessive model, the A allele (homozygous AA) was associated with protection against preeclampsia (OR = 0.71, 95% CI: 0.55–0.92; P = 0.01), while AG/GG genotypes showed elevated risk. This counterintuitive finding — where the lower-vitamin-C genotype appears protective — may reflect complex tissue-specific redox effects in pregnancy or confounding by population-specific factors.

A case-control study of chronic lymphocytic leukaemia⁶⁶ case-control study of chronic lymphocytic leukaemia
Casabonne D et al. Fruit and vegetable intake and vitamin C transporter gene (SLC23A2) polymorphisms in chronic lymphocytic leukaemia. Eur J Nutr, 2017 found a log-additive association between the G allele and CLL risk (OR = 1.19, 95% CI: 1.00–1.41; P = 0.05), independent of fruit and vegetable intake.

Practical Implications

The key finding for most carriers is straightforward: AA homozygotes — about 50% of the global population — run plasma vitamin C levels roughly 6 umol/L lower than GG homozygotes on the same diet. This is a consistent genetic baseline effect that dietary choices can compensate for, but cannot eliminate. If your dietary vitamin C intake is adequate (above 75–90 mg/day), the genotype effect is unlikely to push you into frank deficiency. But if your diet is limited in vitamin C-rich foods — especially common in winter months or during food restriction — the AA genotype adds a structural disadvantage.

The G allele is notably more common in East Asian populations (~64%) than in Europeans (~37%), and quite rare in African populations (~17%). This means the GG "high-ascorbate" genotype affects about 40% of East Asians compared to roughly 14% of Europeans.

Interactions

This variant operates through a different biological step than rs33972313⁷⁷ rs33972313
SLC23A1 Val264Met — reduces intestinal and renal vitamin C transport capacity in SLC23A1 (SVCT1). The Duell 2013 EPIC study demonstrated that both variants independently predicted plasma vitamin C in the same multivariable model, indicating additive rather than redundant effects. A person carrying the reduced-function SLC23A1 variant (rs33972313 CT/TT) alongside the low-vitamin-C SLC23A2 genotype (rs6133175 AA) faces a dual disadvantage: both absorbing less vitamin C from food and having less favorable tissue distribution.

The closely located variant rs6053005⁸⁸ rs6053005
SLC23A2 intronic variant, ~66 kb downstream of rs6133175, also associated with 24% higher plasma vitamin C in TT homozygotes vs CC homozygotes in the same EPIC cohort in SLC23A2 (approximately 66 kb downstream within the same gene) showed nearly identical effects in the EPIC cohort (TT: +24%, beta = 0.21, 95% CI: 0.058–0.37, P = 0.007). These two SLC23A2 variants likely tag the same haplotype block and may not represent fully independent signals.

Deep within chromosome 2, a modest change in a single DNA letter can shift your lifetime odds of developing an intracranial aneurysm — a balloon-like bulge in a brain artery that, if it ruptures, causes a subarachnoid hemorrhage¹¹ subarachnoid hemorrhage
bleeding into the space surrounding the brain, a life-threatening emergency affecting roughly 500,000 people annually worldwide. The variant rs700651 sits in an intron of the BOLL gene (boule-like RNA binding protein) at chromosome 2q33.1 and was first pinpointed in a landmark 2008 genome-wide association study. It remains one of only a handful of robustly replicated genetic risk factors for this condition.

The Mechanism

BOLL encodes an RNA-binding protein most highly expressed in germ cells, but the 2q33.1 locus appears to influence vascular biology through regulatory mechanisms that extend beyond the BOLL coding sequence itself. Cis-eQTL analysis²² Cis-eQTL analysis
a method linking a DNA variant to changes in nearby gene expression at this locus reveals that the risk-allele genotype is associated with altered expression of SF3B1 (a splicing factor subunit) in thyroid and tibial nerve tissue, and ANKRD44 in testis — suggesting the variant influences RNA splicing regulation in tissues relevant to vascular development and maintenance. Additional genome-wide interaction studies identify strong epistatic interactions between rs700651 and rs1105980 in PTCH1 (the Hedgehog pathway receptor), pointing toward a role for developmental vascular patterning pathways in aneurysm susceptibility.

The mechanism by which intronic variants at 2q33.1 predispose to aneurysm formation is still being characterized. The prevailing model involves altered arterial wall homeostasis — the balance between smooth muscle cell integrity, extracellular matrix remodeling, and inflammatory signaling that determines whether a vessel wall can withstand hemodynamic stress. Risk variants at this locus may subtly shift this balance toward wall weakening, particularly at high-flow branch points in the circle of Willis.

The Evidence

The original discovery came from a multistage GWAS³³ multistage GWAS
genome-wide association study: testing millions of DNA variants simultaneously across thousands of people by Bilguvar et al. (2008, Nature Genetics), conducted in Finnish, Dutch, and Japanese cohorts totalling more than 2,100 intracranial aneurysm cases and 8,000 controls. The 2q33.1 locus was among three genome-wide significant hits, with odds ratios of 1.24–1.36 across the identified loci.

An updated meta-analysis⁴⁴ updated meta-analysis
pooling results from multiple independent studies for greater statistical power by Hong et al. (2019, J Clin Med) specifically refined the rs700651 association across 18,019 individuals spanning European, Japanese, and Korean populations: OR 1.213 (95% CI 1.135–1.296), achieving genome-wide significance. This means each copy of the G allele increases intracranial aneurysm risk by approximately 21%.

A Korean GWAS⁵⁵ Korean GWAS by Hong et al. (2019) independently replicated rs700651 alongside rs6841581 (EDNRA), confirming that the 2q33.1 signal is not population-specific. A separate shared-genetics analysis⁶⁶ shared-genetics analysis examined whether the intracranial aneurysm loci also predispose to abdominal or thoracic aortic aneurysm; limited polygenic overlap was found, suggesting the 2q33.1 variant's primary risk is specific to intracranial vessels rather than systemic arterial fragility.

Practical Actions

Carrying one or two G alleles does not determine fate — most G carriers never develop an aneurysm, and most aneurysms are never detected because they remain small and asymptomatic. However, the elevated risk is real and clinically meaningful, especially in combination with modifiable risk factors: hypertension, smoking, and heavy alcohol use are the strongest environmental amplifiers of aneurysm risk, and all are targetable. For GG homozygotes or AG individuals with additional risk factors (family history of aneurysm, polycystic kidney disease, connective tissue disorders), the evidence supports a conversation with a physician about whether a one-time brain MRA screening scan is appropriate.

Interactions

The most clinically relevant interaction is with rs6841581 in EDNRA (endothelin receptor type A, chromosome 4q31.22), which was identified alongside rs700651 in multiple intracranial aneurysm GWAS studies. EDNRA encodes a receptor for endothelin-1, a potent vasoconstrictor that regulates arterial tone and smooth muscle cell proliferation. The two loci appear to operate through distinct but complementary vascular pathways: the 2q33.1 locus likely influences splicing/regulatory biology, while the EDNRA locus directly affects vasoconstrictor signaling. Neither locus directly interacts with the other at the epistatic level — each contributes independently to aneurysm susceptibility. No compound action is warranted because the individual risk contributions are additive rather than synergistic.

A genome-wide interaction analysis identified a strong epistatic signal between rs700651 and rs1105980 in PTCH1 (the Hedgehog pathway receptor, chromosome 9q22.32): lnOR 1.53, p=6.41×10⁻¹¹. This gene-gene interaction suggests that the BOLL locus risk may be particularly amplified in individuals who also carry certain PTCH1 variants — a finding that merits further investigation.

Your brain's resilience to chronic stress depends partly on its ability to maintain and remodel the microscopic structures through which neurons communicate. GPM6A¹¹ GPM6A
Glycoprotein M6a — a tetraspan proteolipid protein in the neuronal membrane, related to the myelin proteolipid protein (PLP) family encodes a protein that does exactly this: it scaffolds the formation of dendritic spines and filopodia — the tiny protrusions on neurons where synapses form. When GPM6A expression falls, synaptic architecture degrades. When stress silences the gene, the brain physically loses some of its wiring. The rs72704544 variant in the GPM6A gene was identified in a landmark 2024 anxiety GWAS, linking genetic variation at this locus to anxiety disorder risk across five continental ancestry groups.

The Mechanism

GPM6A protein localizes to membrane protrusions on hippocampal neurons, where it drives filopodium formation²² filopodium formation
Filopodia are thin actin-rich protrusions that develop into mature dendritic spines — the structural basis of long-term potentiation and memory encoding. In cell culture, overexpressing GPM6A dramatically increases filopodial density; silencing it with siRNA reduces filopodial structures and synaptophysin clusters — markers of functional synapses — at highly significant levels (p<0.0001). This makes GPM6A a structural determinant of synaptic connectivity in the hippocampus, the brain region most critical for emotional regulation and stress responses.

Under chronic stress, the gene is specifically downregulated in the dentate gyrus and CA3 region³³ dentate gyrus and CA3 region
The dentate gyrus is one of the few brain regions that generates new neurons in adults (adult neurogenesis); CA3 is the primary output of the hippocampus involved in stress-memory encoding of the hippocampus — the circuit most vulnerable to stress-induced atrophy. This downregulation occurs alongside suppression of BDNF, the brain's primary neurotrophic growth factor. miR-124-3p⁴⁴ miR-124-3p
A microRNA highly expressed in neurons that promotes neuronal differentiation by suppressing non-neuronal gene programs, a key regulator of GPM6A expression, is also reduced by chronic stress, providing a molecular mechanism for the gene's silencing. BDNF treatment in vitro rescues miR-124-3p and GPM6A expression together, suggesting the BDNF-miR-124-GPM6A axis is a recoverable pathway.

Rs72704544 is an intronic variant — it does not change the GPM6A protein sequence. Its likely mechanism is regulatory: intronic variants near regulatory elements can alter splicing efficiency, transcription factor binding, or enhancer activity, modifying how much GPM6A protein the hippocampus produces, particularly under stress conditions. The G allele at this locus was identified as the risk-increasing variant for anxiety disorders.

The Evidence

The primary human genetic evidence comes from a 2024 multi-ancestry GWAS⁵⁵ 2024 multi-ancestry GWAS
Friligkou E et al. "Gene discovery and biological insights into anxiety disorders from a large-scale multi-ancestry genome-wide association study." Nature Genetics, 2024. by Friligkou et al. in Nature Genetics. The study enrolled over 1.2 million participants including 97,383 anxiety disorder cases across five continental ancestry groups. Fifty-one genome-wide significant loci were identified, 39 of which were novel; heritability enrichment was concentrated in genes expressed in the limbic system, cerebral cortex, and hippocampus. GPM6A was among the 115 genes associated with anxiety through brain-specific transcriptome analysis. The study also documented genetic overlap with depression, schizophrenia, and bipolar disorder — consistent with GPM6A's broad role in synaptic plasticity across mood disorders.

The human postmortem literature provides direct biological validation. A postmortem hippocampal study of 18 depressed suicides⁶⁶ postmortem hippocampal study of 18 depressed suicides
Fuchsova B et al. "Altered expression of neuroplasticity-related genes in the brain of depressed suicides." Neuroscience, 2015. found statistically significant downregulation of GPM6A in depressed individuals (F=14.55, p=0.0002), alongside CAMK2A and CORO1A — all neuroplasticity-related genes. Normal coexpression patterns among these genes were disrupted in depressed brains, suggesting a systems-level failure of hippocampal plasticity maintenance.

In the stress model literature, three weeks of daily restraint stress in rats⁷⁷ three weeks of daily restraint stress in rats
Cooper B et al. "Expression of the axonal membrane glycoprotein M6a is regulated by chronic stress." PLoS One, 2009. consistently downregulates M6a mRNA in dentate gyrus and CA3 neurons. This finding has been independently replicated and extended: Alfonso et al. 2006⁸⁸ Alfonso et al. 2006
Alfonso J et al. "Regulation of hippocampal gene expression is conserved in two species subjected to different stressors and antidepressant treatments." Biological Psychiatry, 2006. showed stress-induced M6a suppression is conserved across species and stressors, and critically, reversed by the antidepressant tianeptine.

Practical Actions

There is no established pharmacogenomic consequence for rs72704544 — it is a risk modifier, not a drug-response determinant. The practical implications concern neuroplasticity maintenance: factors that support BDNF-driven hippocampal remodeling directly address the pathway this variant affects. Among lifestyle exposures with solid evidence for increasing BDNF and GPM6A expression are aerobic exercise (especially sustained moderate-intensity training), omega-3 fatty acids (EPA/DHA), and avoiding sustained elevation of glucocorticoids. Chronic psychological stress specifically suppresses the BDNF–miR-124–GPM6A axis; interventions that reduce HPA-axis hyperactivation are mechanistically relevant.

For individuals with GG genotype (rare, ~3.4% globally) carrying two risk alleles, both clinician awareness of elevated baseline anxiety risk and proactive monitoring of mood symptoms may be warranted.

Interactions

GPM6A is closely related to GPM6B, its paralogous protein; GPM6B directly interacts with the N-terminal domain of the serotonin transporter (SERT, SLC6A4), decreasing SERT cell-surface expression and serotonin reuptake. While GPM6A's direct interaction with SERT is less established than GPM6B's, both proteins share structural homology and hippocampal expression, suggesting a potential functional convergence in serotonergic modulation.

Rs72704544 is worth considering alongside FKBP5 rs1360780 — the classic stress-axis variant. Both variants affect hippocampal gene expression under chronic stress and both have been associated with anxiety and stress-related phenotypes. The FKBP5 variant impairs glucocorticoid receptor feedback; GPM6A affects the downstream structural consequences of glucocorticoid-mediated neuroplasticity suppression. A compound interaction analysis across these two loci would be biologically motivated.

ETV5 (E-Twenty-Six Version 5) is an obesity-associated transcription factor expressed in key brain regions that regulate energy balance, appetite, and food reward. The rs7647305 variant sits in the regulatory region upstream of ETV5 on chromosome 3, and the C allele has been consistently associated with increased BMI and obesity risk across large GWAS.

The Mechanism

ETV5 is a member of the PEA3 group¹¹ PEA3 group
a subfamily of ETS transcription factors involved in development and neural function of ETS transcription factors. In the brain, it is primarily expressed in the arcuate nucleus²² arcuate nucleus
a hypothalamic region containing hunger-sensing and satiety neurons (AGRP/NPY and POMC/CART), the ventromedial hypothalamus³³ ventromedial hypothalamus
a brain region critical for energy homeostasis and satiety signaling, and the ventral tegmental area⁴⁴ ventral tegmental area
the origin of dopaminergic reward neurons that project to the nucleus accumbens.

The rs7647305 variant maps to the predicted TATA-box⁵⁵ TATA-box
a core promoter element that positions RNA polymerase for transcription initiation of the ETV5 promoter, suggesting it directly affects ETV5 transcription levels.

ETV5-deficient mice have reduced body weight, lower fat mass, and are resistant to diet-induced obesity. The gene's expression in hypothalamic nuclei changes with nutritional state — its transcription in the arcuate nucleus and VTA is altered by diet and food availability⁶⁶ altered by diet and food availability
Gutierrez-Aguilar et al. Nutritional state affects the expression of the obesity-associated genes. Obesity, 2012, linking it directly to feeding behavior.

ETV5 also modulates the HPA axis⁷⁷ HPA axis
hypothalamic-pituitary-adrenal axis, the body's central stress response system that regulates cortisol. ETV5-deficient animals show decreased expression of glucocorticoid receptors, mineralocorticoid receptors, and vasopressin receptors in the hypothalamus, resulting in elevated circulating glucocorticoids. This cortisol dysregulation promotes visceral fat deposition and insulin resistance.

The Evidence

The GIANT consortium⁸⁸ GIANT consortium
Willer et al. Six new loci associated with body mass index highlight a neuronal influence on body weight regulation. Nature Genetics, 2009 meta-analysis of over 32,000 individuals identified the ETV5 locus as one of six new genome-wide significant BMI loci (P < 5 x 10^-8). Simultaneously, Thorleifsson et al.⁹⁹ Thorleifsson et al.
Genome-wide association yields new sequence variants at seven loci that associate with measures of obesity. Nature Genetics, 2009 independently confirmed the association.

In 18,014 Danish adults¹⁰¹⁰ 18,014 Danish adults
Haupt et al. Studies of metabolic phenotypic correlates of 15 obesity associated gene variants. PLoS ONE, 2011, the C allele at rs7647305 was associated with an obesity odds ratio of 1.18 (95% CI 1.08-1.29, P = 1.8 x 10^-4) and a per-allele BMI increase of approximately 0.06 kg/m^2.

The association extends beyond BMI: the variant has been independently linked to childhood hypertension¹¹¹¹ independently linked to childhood hypertension
Wang et al. Two obesity susceptibility loci in LYPLAL1 and ETV5 independently associated with childhood hypertension in Chinese population. Gene, 2017 in a Chinese population (OR 0.654 for the T protective allele under a dominant model).

Practical Actions

ETV5 affects obesity through hypothalamic appetite regulation and reward circuitry rather than through peripheral metabolism. This means strategies targeting appetite signaling and cortisol regulation are more relevant than metabolic interventions for carriers.

Interactions

ETV5 rs7647305 contributes to polygenic obesity risk alongside FTO rs9939609, MC4R rs17782313, KCTD15 rs29941, and MTCH2 rs10838738. The ETV5 mechanism is distinct — it operates through central appetite regulation and HPA axis modulation, while FTO affects thermogenesis and MC4R directly modulates satiety neurons. In genetic risk score analyses, individuals carrying risk alleles across multiple loci show cumulative BMI increases of 2-3 kg/m² compared to those in the lowest risk category. The combination of ETV5 (appetite/reward dysregulation) with MC4R (satiety impairment) risk alleles may compound appetite-related effects particularly strongly.

Paraoxonase-1 (PON1) is a calcium-dependent enzyme bound exclusively to HDL particles¹¹ HDL particles
High-density lipoprotein, the "good cholesterol" that transports cholesterol from tissues back to the liver and carries anti-atherogenic enzymes in the bloodstream. Its primary cardiovascular role is preventing LDL from oxidizing — the crucial first step in atherosclerotic plaque formation. rs854555 is an intronic variant within PON1 at chromosome 7q21.3 that does not change the protein sequence itself, but travels in linkage disequilibrium with functional PON1 variants. It is a haplotype tag²² haplotype tag
A variant whose allele can be used to infer the genotype of nearby variants in strong LD, acting as a proxy for a cluster of co-inherited alleles for a low-activity PON1 haplotype, meaning carriers of the A allele tend to have lower overall PON1 enzymatic activity on their HDL particles.

The Mechanism

PON1 activity in plasma is determined by two independent factors: how much enzyme the liver produces (controlled by promoter variants such as rs854571 at −108C>T) and how efficiently the available enzyme works (controlled by coding variants such as rs662 Q192R and rs854560 L55M). rs854555, lying within an intron, does not directly alter transcription or protein function. Instead, its alleles are inherited together with combinations of functional variants — the A allele co-segregating preferentially with low-activity haplotype configurations across diverse populations.

The consequence of reduced PON1 loading on HDL is a decrease in the enzyme's capacity to hydrolyze lipid peroxides accumulating on LDL and HDL particles. When PON1 activity falls, oxidized LDL — the key driver of foam cell formation³³ foam cell formation
Macrophages engulf oxidized LDL to form foam cells, the cellular building blocks of atherosclerotic plaques and arterial plaque — accumulates unchecked. The GWAS Catalog records a genome-wide significant association between rs854555-A and altered response to TNF antagonist therapy (p = 2 × 10⁻⁶), consistent with PON1's documented role in modulating inflammatory signaling via oxidized phospholipid hydrolysis.

The Evidence

Direct evidence for rs854555 as an independent risk variant is limited; its importance lies in haplotype context. A 2017 case-control study in Han Chinese found that the A-A haplotype at rs854555 and rs662⁴⁴ A-A haplotype at rs854555 and rs662
Li et al. Medicine (Baltimore) 2017 (combining the rs854555 A allele with the rs662 Q variant) was significantly associated with increased disease susceptibility (OR 2.74, 95% CI 1.28–5.84), while rs854555 alone showed no independent effect — the hallmark of a haplotype-tagging variant.

The pathway consequence is well-established even where this specific SNP's independent contribution is modest. A meta-analysis of 43 studies comprising 20,629 subjects⁵⁵ 43 studies comprising 20,629 subjects
Zhao et al. Mol Genet Metab 2012 found that reduced PON1 activity is a significant risk factor for coronary heart disease (SMD −0.78, 95% CI −0.98 to −0.57, P<0.001 for paraoxonase activity; SMD −0.50 for arylesterase activity). A second meta-analysis of 20 studies (n=5,417)⁶⁶ 20 studies (n=5,417)
Zuin et al. Dis Markers 2022 confirmed that PON1 arylesterase activity is significantly lower in CAD patients versus controls (SMD −0.587, P<0.0001). Carriers of the rs854555 A allele, by co-inheriting low-activity haplotype configurations, participate in this pathway.

Practical Actions

The most evidence-supported strategy for individuals with low-activity PON1 haplotypes is increasing dietary polyphenol intake. Pomegranate juice consumption for 12 months increased serum PON1 activity by 83% in a controlled study⁷⁷ 83% in a controlled study
Aviram et al. Clin Nutr 2004 of carotid artery stenosis patients, while simultaneously reducing LDL basal oxidative state by 90%. Extra virgin olive oil — through its oleic acid content and minor phenolic compounds — has been shown in multiple trials to upregulate hepatic PON1 mRNA expression and directly stimulate arylesterase activity on HDL particles.

Because rs854555 is a haplotype tag rather than a functional variant, its predictive value is strengthened when interpreted alongside the coding and promoter PON1 variants (rs662, rs854560, rs854571). Individuals with the A allele at rs854555 combined with unfavorable genotypes at those functional sites carry the most reduced total PON1 activity. Direct measurement of serum PON1 arylesterase activity provides a functional readout that integrates all genetic and non-genetic determinants.

Interactions

rs854555 sits within a well-characterized PON1 haplotype structure. Its A allele tends to co-inherit with configurations that reduce functional PON1 enzyme delivery to HDL. The three key functional sites are: rs662 (Q192R)⁸⁸ rs662 (Q192R)
Amino acid substitution that trades LDL antioxidant efficiency for organophosphate hydrolysis speed; the 192R allele is associated with higher CAD risk in multiple meta-analyses, rs854560 (L55M, which reduces PON1 protein stability and serum concentration by >50%), and rs854571 (−108C>T promoter, which controls total transcriptional output of the gene). The combined genotype across all four sites determines an individual's effective PON1 activity far more precisely than any single variant alone.

The IL4R gene¹¹ IL4R gene
Interleukin-4 receptor alpha chain, located at chromosome 16p12.1 encodes the alpha subunit of the IL-4 receptor, which sits at the heart of the body's allergic response circuitry. The rs1805011 variant (also called Glu375Ala in many publications, or Glu400Ala in the canonical full-length transcript) is a missense change in the extracellular domain of this receptor — the region that physically contacts the IL-4 cytokine during binding. Unlike the well-characterized intracellular Q576R variant (rs1801275), which amplifies downstream signaling, Glu375Ala sits upstream in the ligand-interaction zone, with the potential to subtly alter how efficiently IL-4 engages the receptor.

The Mechanism

The IL-4Rα extracellular domain folds into two immunoglobulin-like subdomains that form a binding cradle for [IL-4 | The cytokine interleukin-4, a key driver of Th2-type immune responses and IgE class switching in B cells]. Glutamic acid at position 375/400 sits within this ligand-binding region. Substituting glutamic acid (charged, negatively) for alanine (neutral, nonpolar) changes the local electrostatic environment near the binding interface.

When IL-4 binds IL-4Rα and the receptor complexes with either the common gamma chain (type I receptor)²² common gamma chain (type I receptor)
Forms on lymphocytes and drives Th2 differentiation or IL-13Rα1 (type II receptor)³³ IL-13Rα1 (type II receptor)
Forms on non-hematopoietic cells and airway epithelium, drives mucus and airway remodeling, the intracellular JAK1/TYK2 kinases are activated, phosphorylating STAT6⁴⁴ STAT6
Signal transducer and activator of transcription 6, which moves to the nucleus to drive IgE class switching, Th2 differentiation, and mucin production. The Glu375Ala substitution is hypothesized to alter binding kinetics or receptor conformation in a way that modifies the threshold or magnitude of this signaling cascade — though the precise structural consequence has not been characterized in functional studies to the same degree as the intracellular variants.

The Evidence

Association evidence across several populations supports rs1805011 as a modest risk modifier for atopic disease.

Bottema et al. (2010), analyzing rhinitis and asthma in three cohorts totaling over 700 trios and case-control pairs⁵⁵ Bottema et al. (2010), analyzing rhinitis and asthma in three cohorts totaling over 700 trios and case-control pairs
Family-based and case-control design; European-ancestry populations found that IL4R Glu375Ala (rs1805011) was associated with asthma. Critically, they also found a gene-gene interaction⁶⁶ gene-gene interaction
When the effect of one genetic variant depends on the genotype at a second locus, producing a combined risk larger than either alone between rs1805011 and IL13 Arg130Gln — the two variants together conferred greater asthma risk than either alone, pointing to convergent Th2 pathway amplification.

A meta-analysis of case-control studies by Zhu et al. (2013)⁷⁷ meta-analysis of case-control studies by Zhu et al. (2013)
Six studies examining IL-4 and IL-4R polymorphisms across Asian and European populations confirmed that individuals homozygous for the A allele at rs1805011 were significantly less likely to develop asthma, with the combined C-carrier genotypes (CC+AC) showing an overall OR of 0.39 versus AA (P=0.04), meaning the C allele roughly doubles asthma susceptibility in a dominant model.

In Polish children, Narożna et al. (2016)⁸⁸ Narożna et al. (2016)
177 asthmatic children versus 194 healthy controls found rs1805011 to be the strongest IL4R association in their dataset, reaching genome-wide- suggestive significance for mild asthma (p=0.00005) and a significant signal for atopic dermatitis (p=0.0056).

A genome-wide pharmacogenomic analysis of asthma exacerbations Anderson et al. (2013)⁹⁹ Anderson et al. (2013)
Data from four salmeterol clinical trials, 199 exacerbators vs 502 controls identified rs1805011 among a cluster of IL4R coding variants with consistent genetic effects across three independent studies (P<0.0006), suggesting this variant influences exacerbation biology in addition to initial susceptibility.

Evidence on atopic dermatitis and eczema is more nuanced. A Japanese women's cohort study (Miyake et al. 2013)¹⁰¹⁰ Japanese women's cohort study (Miyake et al. 2013)
188 eczema cases vs 635 controls found a protective association for the C allele in the context of eczema (OR 0.55, 95% CI 0.31-0.99), contrasting with asthma data. This divergence likely reflects condition-specific interactions with other immune loci, environmental exposures, and population differences.

Practical Actions

For carriers of the C allele, the core concern is modestly amplified IL-4 signaling leading to higher IgE class switching and Th2 immune polarization. Monitoring IgE levels quantifies the downstream consequence of altered receptor function. Screening for common aeroallergens and food allergens provides actionable information about which exposures are most likely to trigger clinical responses.

Quercetin directly targets the IL-4/STAT6 signaling axis. In vitro studies confirm that quercetin at 5 micromolar concentrations suppresses IL-4-induced STAT6 phosphorylation¹¹¹¹ In vitro studies confirm that quercetin at 5 micromolar concentrations suppresses IL-4-induced STAT6 phosphorylation
Achievable with oral doses of 1,000 mg reaching plasma levels of 5-12 micromolar, making it a molecularly targeted intervention for Th2 pathway overactivation.

Vitamin D modulates Th1/Th2 immune balance. Optimal levels support regulatory T cell function and suppress Th2-skewed cytokine production, directly counteracting the immune shift that this variant may facilitate.

Interactions

rs1805011 interacts functionally with IL13 Arg130Gln — together they yield greater asthma risk than either variant alone (Bottema et al. 2010), consistent with both variants affecting the same IL-4/IL-13 type II receptor complex.

Within the IL4R gene itself, rs1805011 forms haplotypes with [I75V (rs1805010) | Ile75Val, extracellular domain variant with independent atopy associations] and [Ser503Pro (rs1805015) | Extracellular domain variant also associated with atopic phenotypes in Japanese cohorts]. The combined haplotype across these three extracellular-domain positions may have additive effects on IL-4 binding efficiency. The intracellular variant Q576R (rs1801275) operates through a distinct gain-of-function mechanism; carrying risk alleles at both extracellular and intracellular positions would be expected to compound the overall shift toward amplified Th2 signaling.

Methionine synthase (MTR), also known as MS, catalyzes the final step that converts homocysteine back to methionine using methylcobalamin (active B12) as a cofactor and methylfolate as the methyl donor. This reaction sits at the crossroads of the methylation cycle and is essential for keeping homocysteine levels in check.

The Mechanism

The A2756G variant (rs1805087) causes an aspartic acid-to-glycine substitution ¹¹ Aspartic acid-to-glycine substitution at position 919 of the protein (p.Asp919Gly) at position 919 of the MTR protein. The G allele produces an enzyme with altered activity that tends to favor the active (reduced) state of B12. Paradoxically, this may seem beneficial, but the altered enzyme kinetics can lead to disrupted methylation cycling under certain conditions, particularly when B12 or folate levels are suboptimal. ClinVar classifies this variant as benign given its population frequency.

The Folate Trap

MTR is at the center of what biochemists call the "methyl-folate trap." ²² When MTR is impaired, methylfolate accumulates unusably — a functional folate deficiency despite normal blood levels When MTR activity is impaired, methylfolate accumulates because it cannot donate its methyl group to homocysteine. This creates a functional folate deficiency even when total folate levels appear adequate. Understanding your MTR status helps explain why some people with "normal" folate levels still show signs of impaired methylation.

Clinical Significance

Studies have linked the G allele to altered homocysteine metabolism, though the effects are typically modest. A meta-analysis³³ meta-analysis
Zhao D et al. MTR A2756G and cancer risk, 2010 examined the variant's association with cancer risk across multiple study types. The variant becomes more clinically relevant when combined with MTRR variants (which affect B12 reactivation) and MTHFR variants (which affect methylfolate production). This triad of enzymes works as a coordinated system ⁴⁴ MTR + MTRR + MTHFR form a triad: folate provides the methyl group, B12 carries it, and MTRR keeps B12 active, and weakness at multiple points compounds the effect.

Practical Implications

If you carry the G allele, ensuring generous B12 intake is important since your MTR enzyme has altered B12 handling. Active B12 forms are preferred. Combined with adequate folate (as methylfolate if you have MTHFR variants), this supports optimal homocysteine conversion and methylation cycling.

Interactions

MTR works directly with MTRR (rs1801394) — MTR performs the reaction and MTRR reactivates it. Both interact with MTHFR (rs1801133) as the provider of the methylfolate substrate.