rs1144566

RGS16 RGS16 H137R

Moderate Risk Factor

RGS16 H137R — A Coding Variant in the Circadian Clock's Synchronizer

Every cell in your body runs on a roughly 24-hour molecular clock, but those clocks need a conductor to stay synchronized with each other and with the outside world. That conductor lives in the suprachiasmatic nucleus (SCN)11 suprachiasmatic nucleus (SCN)
A cluster of ~20,000 neurons in the hypothalamus directly above the optic chiasm; it receives light signals from the retina and broadcasts a 24-hour timing signal to every organ in the body via neural, hormonal, and temperature rhythms
— and one of its critical molecular regulators is RGS1622 RGS16
Regulator of G-protein Signaling 16; a GTPase-accelerating protein that terminates Gαi/o signaling, thereby controlling when intracellular cAMP can accumulate in SCN neurons
. The rs1144566 variant changes a single amino acid in RGS16 — histidine to arginine at position 137 — in a gene whose protein levels oscillate daily to gate the cAMP pulses that coordinate clock-neuron communication.

The Mechanism

RGS16 operates at a key node in the GPR176–Gz–RGS16 signaling axis33 GPR176–Gz–RGS16 signaling axis
GPR176 is an orphan GPCR enriched in the SCN; it activates the Gz G-protein subunit (a slow-cycling member of the Gαi family), which suppresses cAMP synthesis; RGS16 terminates Gz signaling by accelerating GTP hydrolysis, releasing the cAMP brake
. Each morning, the molecular clock drives a surge in RGS16 expression; this terminates the Gz-mediated cAMP suppression and allows cyclic AMP to accumulate — a biochemical event essential for synchronizing the phase-leading dorsomedial SCN neurons with the light-receiving ventrolateral neurons.

Doi et al. (2011)44 Doi et al. (2011)
Doi M et al. Nature Communications, 2011
showed that mice lacking RGS16 entirely lose their circadian cAMP rhythm in the SCN and develop a lengthened free-running behavioral period. A longer internal period means the clock drifts toward an evening phase — the same directional shift associated with the C allele at rs1144566 in humans.

The H137R change (His→Arg at position 137) falls within the functional RGS domain responsible for G-protein contact. A positively charged arginine replacing a moderately polar histidine at this position is predicted to reduce GTPase acceleration of Gαi/o, impairing the protein's ability to terminate Gz signaling and thereby weakening the daily cAMP gate that sets circadian timing. This is a missense variant on a gene whose null phenotype directly produces the human evening-type circadian shift, making H137R a plausible partial loss-of-function that nudges carriers toward later timing.

The Evidence

The strongest human genetic evidence comes from two large independent GWAS that identified a chronotype signal at the RGS16 locus in high LD with rs1144566.

Hu et al. (2016)55 Hu et al. (2016)
Hu Y et al. GWAS of 89,283 individuals identifies genetic variants associated with self-reporting of being a morning person. Nature Communications, 2016
found rs12736689 (in LD r²=0.89 with rs1144566) to be the single most significant chronotype hit genome-wide at P=7×10⁻¹⁸ in 89,283 participants, and explicitly noted that rs1144566 is the nearby nonsynonymous coding variant in high LD with the lead SNP. This positions H137R as the most compelling functional candidate underlying the association.

Jones et al. (2016)66 Jones et al. (2016)
Jones SE et al. Genome-Wide Association Analyses in 128,266 Individuals Identifies New Morningness and Sleep Duration Loci. PLoS Genetics, 2016
independently confirmed the RGS16 locus in 128,266 UK Biobank participants, with the C-allele at the linked rs516134 conferring an odds ratio of 1.21 (95% CI 1.15–1.27) for morningness at P=3×10⁻¹². The RGS16 locus is the most replicated single locus in human chronotype genetics, confirmed in at least three independent GWAS totaling over 700,000 individuals.

The mechanistic grounding was extended by Nakagawa et al. (2020)77 Nakagawa et al. (2020)
Nakagawa S et al. Time-Restricted G-Protein Signaling Pathways via GPR176, Gz, and RGS16 Set the Pace of the Master Circadian Clock in the SCN. IJMS, 2020
, which showed that RGS16 protein peaks in the early morning in SCN neurons and that this timed expression is indispensable for the cAMP surge that coordinates cellular synchrony. Human genetic variants of RGS16 associating with earlier wake-up times are highlighted as the translational implication of these rodent findings.

The evidence level is rated moderate rather than strong because the functional impact of H137R specifically has not been characterized in cell or animal models — the GWAS signal is well-established, but whether H137R is the causal variant (versus being a proxy for a linked regulatory variant like rs516134) has not been definitively resolved.

Practical Implications

The C allele (H137R / Arg137) is the overwhelmingly common allele globally (~97.5%), making the CC genotype the population default. The H137R form of RGS16 has modestly reduced GTPase activity compared to the ancestral His137 form, but since nearly everyone carries it, this represents the baseline circadian setting rather than a deviation from normal.

The rare T allele (~2.5% globally) encodes the ancestral His137 form with stronger GTPase activity. Carriers of the T allele (CT or TT genotypes, ~5% of people) have a circadian clock that may run with slightly tighter phase advance — the practical expression is a modest morning preference with easier early waking. The T allele is most common in European (~2.5%) and African (~5.5%) ancestry populations and rarest in East Asian populations (~0.1%).

Interactions

rs1144566 sits at the same RGS16 locus as the regulatory variant rs516134 and the linked rs12736689. Both the coding change (rs1144566) and the regulatory signal (rs516134/rs12736689) independently tag the circadian effect at this locus, and they are in high LD (r²=0.89), so most users with the rs1144566 T allele will also carry the morningness-associated allele at rs12736689.

For additive circadian effects, variants in CLOCK (rs1801260) and PER3 (rs5751876) act at different nodes of the same oscillator. Carriers of morningness alleles at both RGS16 (T) and CLOCK may experience stronger morning preference than either variant alone predicts, as the two genes affect different parts of the feedback loop — RGS16 affects intercellular cAMP synchrony while CLOCK affects transcription-factor stability in the core loop.

CYP2J2 rs11572325 — A Cardiac Epoxygenase Variant and MI Risk

Your heart produces its own endogenous cardioprotective molecules. CYP2J2 is a cytochrome P450 enzyme expressed predominantly in cardiomyocytes and vascular endothelial cells, where it converts arachidonic acid into epoxyeicosatrienoic acids (EETs)11 epoxyeicosatrienoic acids (EETs)
Biologically active lipid mediators produced by cytochrome P450 epoxygenases; they dilate coronary arteries, reduce vascular inflammation, inhibit platelet aggregation, and protect cardiac tissue from ischemic injury
. The intronic variant rs11572325 sits within an intron of CYP2J2 on chromosome 1 and forms part of a haplotype block with the promoter variant rs890293 (CYP2J2*7) and the neighboring intronic SNP rs2280275 — all three of which tag the same reduced-CYP2J2-expression signal. Population studies have linked the rs11572325 T allele to both myocardial infarction risk and female-specific essential hypertension.

The Mechanism

CYP2J2 converts arachidonic acid into four EET regioisomers (5,6-EET, 8,9-EET, 11,12-EET, and 14,15-EET). These lipid mediators activate K⁺ channels in vascular smooth muscle to cause vasodilation, suppress nuclear factor-κB signaling to reduce endothelial inflammation, inhibit platelet aggregation, and reduce cardiomyocyte apoptosis after ischemic stress. [| Lai & Chen 2021 (PMID 33716791): EETs from CYP2J2 protect against cardiac hypertrophy, fibrosis, and apoptosis; raising EET levels via CYP2J2 overexpression or sEH inhibition reproducibly reduces myocardial injury in experimental models]. The rs11572325 T allele is in strong [linkage disequilibrium | LD means two alleles are inherited together so frequently that one allele predicts the other across the population] with the CYP2J2*7 promoter variant (rs890293), which disrupts an Sp1 transcription factor binding site and reduces CYP2J2 promoter activity by approximately 50%. T-allele carriers show measurably lower circulating EET metabolite concentrations (14,15-DHET), attenuating the vasodilatory and anti-inflammatory signaling that normally protects the coronary vasculature.

CYP2J2 also preferentially metabolizes omega-3 fatty acids EPA and DHA into 17,18-EEQ and 19,20-EDP22 17,18-EEQ and 19,20-EDP
Omega-3 epoxides with antiarrhythmic and vasodilatory properties; EPA is metabolized by CYP2J2 at approximately 17-fold higher efficiency than arachidonic acid
. This substrate competition means that high dietary omega-3 intake can shift eicosanoid output toward cardioprotective mediators even when EET-generating capacity is reduced.

The Evidence

Myocardial infarction risk: The best-powered study comes from a population-based case-control analysis33 population-based case-control analysis
Marciante et al. Common variation in cytochrome P450 epoxygenase genes and MI risk. Pharmacogenet Genomics, 2008
of 856 incident MI cases and 2,688 controls. The rs11572325 T allele was associated with increased MI risk (OR 1.27, 95% CI 1.08–1.51, P=0.006, Benjamini-Hochberg q=0.090), surviving multiple-testing correction in a panel of 30 tag-SNPs across three CYP epoxygenase genes. Notably, the companion CYP2J2 variant rs10889160 showed an independent association (OR 1.24, 95% CI 1.07–1.43, P=0.004), while neither CYP2C8 nor CYP2C9 variants reached significance — pointing specifically to CYP2J2 as the locus driving cardiovascular EET risk. No association with ischemic stroke was observed.

Hypertension — female-specific effect: A comprehensive Russian population study44 comprehensive Russian population study
Polonikov et al. 2019, n=2,314 (discovery + replication cohorts)
genotyped eight CYP2J2 SNPs including rs11572325. The T allele was associated with essential hypertension in women (OR 1.89, 95% CI 1.22–2.95) in the discovery cohort, with directional replication in an independent Belgorod cohort. Critically, no association was detected in men, and haplotype analysis placed rs11572325 within the highest-risk female haplotype (T–T–G–C–C–C–T–A across eight CYP2J2 SNPs). This sex-dependent effect likely reflects estrogen's modulatory role in arachidonic acid metabolism and CYP2J2 transcriptional sensitivity.

Functional EET reduction: The CYP2J2*7 promoter variant (rs890293), in strong LD with rs11572325, was shown to [reduce promoter activity by ~50% | Liu et al. 2007 (PMID 17126841): CYP2J2*7 T carriers had lower plasma EET metabolites; smoking synergistically increased MI risk 6.7-fold] and produce lower circulating EET concentrations. Smokers carrying the T allele exhibited a 6.7-fold amplification of MI risk relative to non-smoking wild-type carriers (P=0.01 for interaction) — highlighting that environmental exposures compound the genetic EET deficit.

Practical Actions

The most specific actionable intervention is dietary: because CYP2J2 metabolizes EPA and DHA into cardioprotective omega-3 epoxides at far higher efficiency than arachidonic acid, supplementing with concentrated fish oil or algae-based EPA/DHA shifts eicosanoid output toward 17,18-EEQ and 19,20-EDP even when EET-generating capacity is reduced. Cardiovascular monitoring is also warranted, particularly in women, given the replicated hypertension signal. Smoking dramatically amplifies the EET-reduction risk and represents the highest-leverage modifiable factor for carriers.

Interactions

rs11572325 is part of the CYP2J2 haplotype block that includes the promoter variant rs890293 (CYP2J2*7) and the intronic variant rs2280275. Studies of all three variants consistently find similar effect directions — they tag the same reduced-expression signal. The companion MI-associated variant rs10889160 showed an independent association in the same Marciante 2008 study, suggesting the two intronic variants may contribute additively to cardiovascular risk through distinct LD blocks within CYP2J2. The soluble epoxide hydrolase variant rs751141 (EPHX2) modulates EET degradation downstream: reduced EPHX2 activity raises EET levels, potentially offsetting reduced CYP2J2 production. The CYP2C8 epoxygenase variant rs1058932 contributes to the same vascular EET pool via a different gene — together, variants in CYP2J2 and CYP2C8 define an individual's total coronary EET production capacity.

rs1159327

ESR1 ESR1 intron variant

Moderate Risk Factor

ESR1 Intron Variant — Estrogen Receptor Expression and Bone Health

The ESR1 gene encodes estrogen receptor alpha (ERα), the primary nuclear receptor through which estrogen orchestrates bone remodeling, reproductive function, and metabolism. ERα is expressed throughout the skeleton, uterus, ovary, breast, and cardiovascular system — making ESR1 variants some of the most broadly consequential polymorphisms in reproductive-age women. The rs1159327 variant sits deep within an ESR1 intron on chromosome 6 (position 151,726,887 GRCh38) and has been identified at genome-wide significance for bone mineral density11 identified at genome-wide significance for bone mineral density
He et al. Osteoporosis International 2023, n=141,261 UK Biobank participants, p=5×10⁻⁹
.

The Mechanism

As an intronic variant, rs1159327 does not alter the ESR1 protein sequence. Instead, it likely influences gene expression through one of several regulatory mechanisms common to intronic SNPs: disruption of transcription factor binding sites22 transcription factor binding sites
protein complexes that bind DNA and control when and how much a gene is transcribed
, alteration of enhancer activity, or effects on mRNA splicing efficiency. The ESR1 gene is large (~300 kb) with numerous regulatory elements distributed across its introns, and the intron 1 region in particular has been repeatedly implicated in controlling receptor expression levels.

ERα regulates osteoblast (bone-forming cell) proliferation and survival while suppressing osteoclast (bone-resorbing cell) activity. When ERα expression or sensitivity is reduced — as intronic regulatory variants can cause — the balance tips toward net bone resorption. The T allele at rs1159327 is associated with lower average bone mineral density, consistent with a model in which reduced ESR1 expression impairs estrogen's bone-protective signaling.

The Evidence

The primary evidence linking rs1159327 to bone health comes from a 2023 longitudinal GWAS in the UK Biobank33 2023 longitudinal GWAS in the UK Biobank
He et al., Osteoporosis International, n=141,261
, which used trajectory analysis (TrajGWAS) to model heel bone mineral density over time. The ESR1 locus emerged at genome-wide significance (p=5×10⁻⁹) for BMD mean, placing rs1159327 among a handful of ESR1-region variants robustly associated with skeletal health in large populations.

This finding sits within a rich broader literature on the ESR1 locus. A meta-analysis of GWAS in premenopausal women44 A meta-analysis of GWAS in premenopausal women
Koller et al. JBMR 2013, discovery n=4,061 + replication n=5,597
identified multiple independent signals in the ESR1 region, and a separate GWAS in 2,468 European men55 a separate GWAS in 2,468 European men
Holliday et al. PLoS One 2011
confirmed the 6q25 locus affects both calcaneal ultrasound BMD and radial volumetric BMD. The largest early meta-analysis — 18,917 individuals across 8 European centers66 18,917 individuals across 8 European centers
Ioannidis et al. JAMA 2004
— found that ESR1 variants affect fracture risk by mechanisms independent of baseline bone density, suggesting the estrogen receptor pathway influences bone quality and geometry as well as mass.

Population data from dbSNP (ALFA) shows the T allele frequency varies from ~17% in South Asian populations to ~45% in African populations, meaning TT homozygosity (~8% globally) is substantially rarer in South/East Asian ancestry groups but more common in individuals of African descent.

Practical Implications

The most direct implications of the T allele relate to bone health across the lifespan. Carriers of one or two T alleles have a tendency toward lower bone mineral density, which translates to heightened importance of bone-protective behaviors — particularly weight-bearing exercise, calcium intake, vitamin D sufficiency, and proactive bone density monitoring. The effect of ESR1 variants on fracture risk appears to operate partly independently of BMD, which means a normal DEXA scan does not fully reassure against elevated fracture risk in T allele carriers.

For women, the ESR1 locus is also relevant to reproductive outcomes. Although the strongest fertility associations in ESR1 involve the nearby intron-1 variants rs2234693 (PvuII) and rs9340799 (XbaI), these three variants are in linkage disequilibrium and likely index overlapping biology. ESR1 polymorphisms have been examined in endometriosis, IVF outcomes, and ovarian reserve studies — the receptor's role in folliculogenesis, endometrial receptivity, and implantation means that even modest changes in expression can influence reproductive efficiency.

Interactions

rs1159327 sits near two of the most studied ESR1 variants: rs2234693 (PvuII, intron 1) and rs9340799 (XbaI, intron 1). These three variants span a regulatory stretch of ESR1 and are in partial linkage disequilibrium — meaning carriers of the rs1159327 T allele are more likely to also carry the rs2234693 T allele and rs9340799 G allele. Haplotype analyses of the ESR1 region show that combined effects on bone, cardiovascular, and hormone-responsive tissues can differ substantially from any single variant's individual contribution. Individuals carrying T alleles at multiple ESR1 loci may have more pronounced reductions in receptor expression.

MCF2L and Osteoarthritis: How a Synovial Gene Shapes Joint Vulnerability

Your joints don't simply wear down mechanically — they're governed by complex molecular signals that determine how the tissues lining them respond to stress, injury, and inflammation. MCF2L, a gene encoding a guanine nucleotide exchange factor11 guanine nucleotide exchange factor
a protein that activates Rho-family GTPases, regulating cytoskeletal dynamics and cell motility
, turns out to be one of those governors. A common intronic variant in MCF2L — rs11842874 — quietly shapes how much of this protein your synovial tissue produces, and with it, your baseline susceptibility to osteoarthritis of the knee, hip, and other large joints.

The Mechanism

rs11842874 sits within intron 4 of the MCF2L gene on chromosome 13 and does not change the protein sequence at all. Instead, it acts as a cis-acting expression quantitative trait locus (eQTL)22 cis-acting expression quantitative trait locus (eQTL)
a variant that modulates how much of a nearby gene's mRNA is produced, without altering the gene's coding sequence
— specifically in synovial membrane tissue, the thin layer lining joint cavities that secretes lubricating fluid and houses key immune and structural cells.

Shepherd et al. (2015)33 Shepherd et al. (2015) demonstrated through quantitative PCR and RNA sequencing of joint tissues that possession of the A allele correlates with increased MCF2L expression in synovial membrane — but not in cartilage. This tissue-specific eQTL pattern was notable: it was the first OA susceptibility locus shown to operate through the synovium rather than cartilage. Luciferase assays confirmed that several SNPs in linkage disequilibrium with rs11842874 display quantitative differences in regulatory activity at the allelic level, pointing to functional regulatory elements within this intronic block.

MCF2L's connection to osteoarthritis biology runs through nerve growth factor (NGF) signaling44 nerve growth factor (NGF) signaling. MCF2L is involved in neurotrophin-mediated regulation of cell motility in the peripheral nervous system, and NGF is a potent mediator of joint pain and inflammation in OA. This mechanistic link gained clinical credibility when humanized anti-NGF antibodies (such as tanezumab) were shown to reduce pain and improve function in knee OA patients — validating the pathway.

The Evidence

The association between rs11842874 and osteoarthritis is one of the most robustly replicated in the field. Day-Williams et al. (2011)55 Day-Williams et al. (2011) conducted a staged GWAS followed by large-scale replication across European cohorts, ultimately encompassing 19,041 OA cases and 24,504 controls. The combined odds ratio was 1.17 (95% CI: 1.11–1.23, p=2.1×10⁻⁸) per A allele — a genome-wide significant signal and the third established locus for OA overall at the time of publication. Given the additive model, AA homozygotes have approximately 1.17² ≈ 1.37-fold elevated risk relative to GG.

The variant also shows evidence of modifying radiographic severity. Valdes et al. (2012)66 Valdes et al. (2012) genotyped three UK cohorts and found rs11842874 was nominally associated with patellofemoral Kellgren-Lawrence grade as a quantitative trait (p=0.027), though the effect on tibiofemoral severity was less consistent. This suggests the MCF2L locus influences both OA susceptibility and the progression of structural joint damage.

The A allele frequency in populations of European descent is approximately 92%, meaning the risk genotype (AA) is by far the most common configuration — around 78% of Europeans carry it. The protective GG genotype, by contrast, is present in only about 1% of Europeans.

Practical Actions

Carrying the A allele doesn't mean OA is inevitable — it means you have population-typical (or slightly elevated) susceptibility that is worth managing proactively, especially if you're physically active or have other joint risk factors. The primary leverage points are reducing cumulative mechanical stress on joint cartilage and supporting synovial tissue health.

High-impact loading patterns — particularly repetitive eccentric loading like downhill running and jumping — generate compressive forces that accelerate cartilage degradation in susceptible joints. For A allele carriers, the MCF2L-driven synovial response to this stress may be dysregulated. Prioritizing low-impact cardio alternatives (cycling, swimming, elliptical) during high-volume training blocks distributes joint load more favorably.

Type II collagen peptides and undenatured type II collagen (UC-II) have specific evidence in supporting cartilage integrity distinct from the type I collagen peptides used for tendon health. The synovial mechanism in MCF2L-associated OA makes this particularly relevant — the synovium is where the eQTL effect is active. Several clinical trials in OA-susceptible individuals have shown benefit from 40 mg UC-II daily (Arthritis & Rheumatism, 2009) and from 10 g hydrolyzed type II collagen.

Interactions

The MCF2L OA susceptibility locus acts within a broader genetic architecture. Other well-replicated OA GWAS loci include GDF5 (rs143383), which encodes growth differentiation factor 5 and influences joint development and cartilage homeostasis. If you carry risk alleles in both MCF2L and GDF5, the combined load on your joint biology is greater than either alone — though formal compound analyses are limited. The NGF pathway connection also creates potential interaction with pain-processing variants; individuals with both elevated MCF2L expression (A/A at rs11842874) and sensitized pain signaling may have amplified OA symptom burden beyond structural changes alone.

ACADM Arg281Thr — A Pathogenic MCAD Variant Unmasked by Fasting

Every cell that depends on fat for fuel must first strip two-carbon acetyl groups from fatty acids inside mitochondria, a process called beta-oxidation11 beta-oxidation
The sequential enzymatic removal of acetyl-CoA units from the fatty acid chain, producing NADH and FADH2 for the electron transport chain and feeding acetyl-CoA into the citric acid cycle
. The medium-chain step of this pathway — handling fatty acids with 6–12 carbons — is catalysed exclusively by medium-chain acyl-CoA dehydrogenase (MCAD), the protein encoded by ACADM. When MCAD is absent or non-functional, medium-chain fatty acids accumulate, toxic acylcarnitine species build up, and the liver cannot produce enough ketone bodies to spare the brain during fasting. The result — if unrecognised — is hypoketotic hypoglycemia progressing rapidly to seizures, coma, and death.

The c.842G>C variant (rs121434282) exchanges arginine for threonine at position 281 of the MCAD protein (p.Arg281Thr). ClinVar classifies it as Pathogenic (VCV000003596, two-star review status, multiple submitters, no conflicts). It is one of more than 80 pathogenic ACADM variants but one of only three — alongside the common K329E (c.985A>G, rs77931234) and a second minority allele (rs121434281) — that together account for the great majority of disease alleles identified in MCAD deficiency patients.

The Mechanism

Arginine-281 sits within the MCAD active site22 MCAD active site
The catalytic core where FAD (flavin adenine dinucleotide) is positioned to accept electrons from the fatty acid substrate during the dehydrogenation reaction
. Arg281 forms critical salt bridges and hydrogen bonds that stabilise the binding of FAD33 FAD
Flavin adenine dinucleotide — the essential cofactor that shuttles electrons from the fatty acid to the mitochondrial electron transport chain during each oxidation cycle
within the enzyme's active site. Replacing arginine (basic, positively charged) with threonine (neutral, hydroxyl-bearing) disrupts the electrostatic network holding FAD in place. The result is an enzyme with severely impaired or absent catalytic activity: medium-chain substrates cannot be oxidised, and the block in fat oxidation becomes clinically apparent whenever carbohydrate reserves are depleted — during overnight fasting, febrile illness, vomiting, or any state that forces the body to switch to fat as its primary fuel.

Because MCAD deficiency is autosomal recessive, a single Arg281Thr allele paired with one normal allele (heterozygous carrier state) leaves one functional MCAD gene copy, which provides sufficient enzyme activity under normal conditions. Two pathogenic alleles — whether homozygous Arg281Thr or compound heterozygous44 compound heterozygous
Carrying two different pathogenic alleles of the same gene, one on each chromosome — the most common configuration in MCAD deficiency since K329E is the predominant allele and most affected individuals pair it with a second, rarer variant
with another MCAD variant — are required to produce clinical disease.

The Evidence

MCAD deficiency is the most common inherited disorder of mitochondrial fatty acid oxidation. Before universal newborn screening, it carried substantial mortality: a large English screening study55 a large English screening study
Oerton J et al. Newborn screening for medium chain acyl-CoA dehydrogenase deficiency in England: prevalence, predictive value and test validity based on 1.5 million screened babies. J Med Screen, 2011
estimated that up to a quarter of previously undiagnosed children die during their first acute metabolic episode, with a further 16% surviving with severe neurological disability. Prevalence in England is approximately 1 in 10,000 births.

The Arg281Thr variant was first reported by Albers et al.66 Albers et al.
Albers S et al. Compound heterozygosity in four asymptomatic siblings with medium-chain acyl-CoA dehydrogenase deficiency. J Inherit Metab Dis, 2001
in four siblings identified through expanded newborn screening by tandem mass spectrometry. All four were compound heterozygotes — carrying one Arg281Thr allele and one K329E allele — and were asymptomatic at identification despite laboratory confirmation of MCAD deficiency. The authors noted a potentially milder clinical phenotype for this genotype combination and emphasised the importance of cascade testing siblings of any newly screened infant. Despite the milder phenotype observed in this family, ClinVar classifies this variant as pathogenic and clinical management guidelines treat all biallelic ACADM genotypes with the same fasting-avoidance and emergency protocols.

A comprehensive 2023 review by Mason et al.77 Mason et al.
Mason E et al. Medium-chain Acyl-COA dehydrogenase deficiency: Pathogenesis, diagnosis, and treatment. Endocrinol Diabetes Metab, 2023
confirms that MCAD deficiency prognosis is excellent once diagnosed and managed appropriately, but that acute metabolic decompensation must be treated as a medical emergency requiring immediate glucose administration.

Practical Actions

The cornerstone of MCAD deficiency management is preventing prolonged fasting. Age-based fasting limits are: no more than 3–4 hours for infants, 4–6 hours for toddlers, and up to 8 hours for older children and adults. Every affected individual should carry an emergency letter and medical alert identification explaining their diagnosis and the need for immediate glucose administration during acute illness. L-carnitine supplementation (50–100 mg/kg/day) is commonly prescribed to prevent secondary carnitine depletion from acylcarnitine sequestration, though clinical trials have not definitively demonstrated benefit on outcome endpoints.

Medium-chain triglyceride oils (MCT oil, coconut oil in large amounts) should be avoided because they directly load the blocked pathway with medium-chain fatty acid substrates.

Interactions

The Arg281Thr allele is most commonly observed in trans with the K329E allele (rs77931234), the dominant European MCAD disease allele. This compound heterozygous combination (Arg281Thr/K329E) appeared to confer a milder biochemical phenotype in the originally reported family, though clinical management should follow the same guidelines as classical MCAD deficiency.

Any combination of two pathogenic ACADM alleles — whether Arg281Thr/Arg281Thr, Arg281Thr/K329E, or Arg281Thr paired with rs121434281 or another pathogenic allele — is sufficient to cause MCAD deficiency. The interaction between two biallelic pathogenic variants is the defining genetic event; single-allele carrier status does not cause clinical disease.

rs121918393

APOE Christchurch (R136S)

Strong Protective

The Christchurch Mutation — Nature's Blueprint for Alzheimer's Resistance

Apolipoprotein E (ApoE) is the most powerful common genetic determinant of Alzheimer's disease risk in the human genome. The APOE4 allele increases lifetime risk 3–4 fold in heterozygotes and 8–12 fold in homozygotes; APOE2 is modestly protective. But in 2019, a third piece of the puzzle emerged: an ultra-rare variant called APOE3 Christchurch (R136S), carried in homozygous form by a single Colombian woman, appeared to nearly abolish Alzheimer's disease pathology despite her carrying the most penetrant familial AD mutation known — PSEN1 E280A, which causes dementia in virtually all carriers by their late 40s.

Arboleda-Velasquez et al. 201911 Arboleda-Velasquez et al. 2019
Resistance to autosomal dominant Alzheimer's disease in an APOE3 Christchurch homozygote: a case report. Nature Medicine
reported that this woman did not develop mild cognitive impairment until her 70s — nearly 30 years after the expected onset of ~44 years in her 1,200-member kindred. Brain imaging revealed unusually high amyloid burden but remarkably limited tau pathology and preserved metabolism in AD-vulnerable regions. The implication: the Christchurch mutation somehow broke the link between amyloid accumulation and downstream tau-driven neurodegeneration.

The Mechanism

The R136S substitution — arginine to serine at position 136 of the mature protein (154 in the pre-protein, hence the dual nomenclature R136S / p.Arg154Ser) — falls squarely within the heparan sulfate proteoglycan (HSPG) binding domain22 heparan sulfate proteoglycan (HSPG) binding domain
HSPGs are abundant cell-surface and extracellular matrix proteins that act as co-receptors for ApoE
. The Christchurch substitution dramatically weakens HSPG affinity.

This matters because HSPG binding drives several of ApoE's most damaging effects in the AD brain. Wang et al. 2023 (Nature Neuroscience)33 Wang et al. 2023 (Nature Neuroscience) showed that the homozygous R136S mutation rescued APOE4-driven phosphorylated tau accumulation, neuroinflammation, and neurodegeneration in both mouse tauopathy models and human iPSC-derived neurons. Single-nucleus RNA sequencing demonstrated that R136S increased disease-protective and reduced disease-associated cell populations in a gene-dose-dependent manner — homozygotes showed full rescue, heterozygotes partial.

Two additional mechanisms have since been characterized. Zhang et al. 2024 (Neuron)44 Zhang et al. 2024 (Neuron) found that the R136S mutant protein directly binds tau with higher affinity than wild-type ApoE3, blocking tau uptake into neurons and microglia and reducing tau fragmentation by the asparagine endopeptidase AEP — a key step in tau propagation. Separately, Wang et al. 2024 (Immunity)55 Wang et al. 2024 (Immunity) demonstrated that R136S inhibits the cGAS-STING interferon pathway, suppressing the chronic neuroinflammatory signaling that amplifies tau pathology.

On the lipid side, the same HSPG/LDL-receptor-related protein binding domain affected by R136S is required for efficient clearance of triglyceride-rich remnant lipoproteins. This means heterozygotes and homozygotes for the Christchurch allele can develop type III hyperlipoproteinemia66 type III hyperlipoproteinemia
a.k.a. dysbetalipoproteinemia; impaired clearance of IDL and VLDL remnants causing elevated total cholesterol and triglycerides
, particularly when combined with other dyslipidemia-predisposing factors.

The Evidence

The founding 2019 case report was remarkable but N=1. Subsequent research has substantially broadened the evidence base.

Saez-Calveras et al. 2023 (Mol Neurodegeneration)77 Saez-Calveras et al. 2023 (Mol Neurodegeneration) analyzed 455,306 UK Biobank participants and identified 37 heterozygous APOEch carriers (36 European, 1 admixed American; median age 68.6 years). None had developed AD or MCI by the data freeze. Carriers showed lower apolipoprotein B levels and a significantly reduced polygenic risk score for AD (p=0.02). The allele frequency in UK Biobank was approximately 0.004%.

Lopera et al. 2024 (NEJM)88 Lopera et al. 2024 (NEJM) quantified the heterozygous effect in the Colombian PSEN1 E280A kindred: among 27 heterozygous APOEch carriers, median MCI onset was 52 years versus 47 years in matched non-carriers — a 5-year delay, statistically significant (p<0.001). Two heterozygous carriers who underwent brain imaging showed relatively preserved FDG-PET metabolic activity in AD-vulnerable regions. Four autopsy specimens from APOEch carriers showed fewer vascular amyloid deposits. The protective effect is real in heterozygotes, but substantially more modest than in the original homozygous case.

The overall allele frequency is approximately 1–4 per 100,000 in European populations and even rarer in other ancestries. The probability of inheriting two copies (homozygous, as in the Colombian case) is therefore roughly 1 in 600 million to 1 in 6 billion — making the original case effectively unique in human clinical history.

Practical Actions

For the overwhelming majority of people tested, this SNP will return the wild-type CC genotype — the Christchurch allele is too rare to appear in most population samples. The primary clinical relevance is for individuals who carry the AC genotype (heterozygous):

The 5-year delay in PSEN1 E280A carriers is meaningful and suggests that even a single copy of the Christchurch allele provides real neurological protection, likely through partial reduction in HSPG binding and partial tau propagation block. Heterozygous carriers should monitor lipid panels carefully, as the same HSPG/LDL-receptor binding impairment that protects the brain can impair remnant lipoprotein clearance. This is not a reason to treat aggressively, but it is a reason to track.

The Christchurch mutation is now the leading molecular template for next-generation Alzheimer's drug development. Mimetic antibodies that block the HSPG-binding domain of wild-type ApoE are in preclinical development, and the variant itself is being explored for AAV gene delivery as a therapeutic.

Interactions

The Christchurch mutation's protection is most dramatic on an APOE3 backbone (the original Colombian case was APOE3/APOE3-ch). Subsequent work has shown R136S is also protective on an APOE4 background — Wang et al. 202399 Wang et al. 2023 demonstrated R136S on APOE4 rescued most tau and neuroinflammation endpoints in model systems, though with somewhat attenuated effect compared to APOE3 background.

The two canonical APOE isoform SNPs — rs4293581010 rs429358 (the ε4-defining Cys112Arg) and rs74121111 rs7412 (the ε2-defining Arg158Cys) — both fall in the LDL receptor binding region of ApoE and influence AD risk through overlapping but distinct mechanisms. The Christchurch mutation falls in the HSPG-binding subdomain adjacent to the LDL receptor binding site, and its protection is at least partially independent of APOE isoform, as shown by the iPSC data generating APOE2ch, APOE3ch, and APOE4ch constructs.

For individuals with a family history of early-onset Alzheimer's disease, both rs759326281212 rs75932628 (TREM2 R47H) and rs38511791313 rs3851179 (PICALM) influence microglial clearance of amyloid and tau — the same downstream pathways Christchurch modulates — and would be relevant to examine in conjunction.

Protein S Arg561Trp — An Anticoagulant Protein That Cannot Leave the Cell

Protein S is one of the body's key natural anticoagulants. It works as a [cofactor for activated protein C | Activated protein C (APC) inactivates coagulation factors Va and VIIIa; protein S amplifies this process roughly 10-fold (https://pubmed.ncbi.nlm.nih.gov/19809585/)11 https://pubmed.ncbi.nlm.nih.gov/19809585/)], which inactivates the clotting factors that would otherwise sustain a thrombus. When protein S is reduced or absent, activated protein C loses much of its braking power over the coagulation cascade — the balance shifts toward clot formation and thrombosis. The Arg561Trp variant (p.Arg561Trp; c.1681C>T) sits in the C-terminal sex hormone-binding globulin (SHBG)-like domain of protein S, a region essential for the protein's structural stability and its ability to be secreted from liver cells into the bloodstream.

The Mechanism

The arginine-to-tryptophan substitution at position 561 disrupts the folding of the SHBG-like domain22 SHBG-like domain
The C-terminal region of protein S is structurally homologous to sex hormone-binding globulin; it is composed of two laminin G-like (LG) domains that coordinate calcium binding and mediate key protein-protein interactions
. Missense mutations in this domain characteristically cause impaired intracellular secretion33 impaired intracellular secretion
The misfolded mutant protein is retained in the endoplasmic reticulum and degraded rather than exported; this is the predominant mechanism for quantitative Type I protein S deficiency
— the misfolded mutant protein never reaches the circulation. Heterozygous carriers produce approximately 50% of normal protein S levels from their one functional allele, which is sufficient for most daily demands but leaves little reserve when the coagulation system is activated by surgery, immobility, pregnancy, or oral contraceptives.

The Arg561Trp variant was identified by Li et al. (2019)44 Li et al. (2019)
Li L et al. Clinical Manifestation and Mutation Spectrum of 53 Unrelated Pedigrees with Protein S Deficiency in China. Thrombosis and Haemostasis, 2019
as one of three recurrent hotspot mutations in the Chinese protein S deficiency population. In that cohort of 53 pedigrees, over half of probands (52.8%) experienced recurrent or multi-site thrombotic events, with deep venous thrombosis and pulmonary embolism comprising 82.7% of cases.

The Evidence

Hereditary protein S deficiency is an established thrombophilia. Heterozygous PROS1 loss-of-function increases venous thrombosis risk approximately 10- to 20-fold55 approximately 10- to 20-fold
Heeb MJ. Role of the PROS1 gene in thrombosis: lessons and controversies. Thrombosis Research, 2009
, compared to the general population. Among all individuals carrying a protein S deficiency-causing PROS1 variant, approximately 50% develop at least one VTE event in their lifetime, while the other 50% remain asymptomatic — reflecting the incomplete penetrance characteristic of this condition. Annual incidence of first venous thrombosis in affected individuals is approximately 1.9%, rising to 6-10% risk of recurrence after an index event.

A large population-scale study found that likely-deleterious missense variants in PROS1 (including Arg561Trp class changes) were associated with OR 1.98 for VTE66 OR 1.98 for VTE
compared to OR 14.01 for complete loss-of-function variants (nonsense, frameshift, essential splice site)
, while complete loss-of-function variants carried a dramatically higher OR of 14.01.

Across East Asian populations, PROS1 mutations are substantially more prevalent in VTE patients than in Caucasian cohorts. Japanese studies found PROS1 mutation frequencies 5-10 times higher in thrombosis patients77 5-10 times higher in thrombosis patients
Kinoshita et al., Clinical Biochemistry, 2005
versus Caucasian populations, and Chinese cohort data consistently place protein S deficiency among the leading inherited thrombophilias in the region.

Free protein S antigen — not total protein S or the functional assay — is the recommended first-line laboratory test for detecting heterozygous deficiency. The International Society on Thrombosis and Haemostasis recommends confirming two or more abnormal results at least four weeks apart before diagnosing hereditary deficiency.

Practical Actions

Heterozygous carriers should have free protein S quantitation measured under appropriate conditions (not during acute illness, pregnancy, or while taking warfarin, as these artificially suppress levels). A confirmed low result, combined with the genetic finding, supports thrombophilia specialist referral for personalized VTE risk stratification.

For carriers who have had a VTE event, direct oral anticoagulants (DOACs such as rivaroxaban or apixaban) are now preferred over warfarin for long-term anticoagulation. Vitamin K antagonists like warfarin suppress all vitamin K-dependent proteins including protein S itself, complicating monitoring and occasionally triggering warfarin-induced skin necrosis in protein S-deficient individuals. DOACs avoid this interaction.

High-risk situations requiring VTE prophylaxis discussion with a physician include: major surgery, prolonged immobilization, air travel exceeding eight hours, pregnancy and the postpartum period, and initiation of estrogen-containing contraceptives or hormone replacement therapy.

Homozygous carriers (AA) face severely reduced or absent protein S and are at risk for neonatal purpura fulminans, a life-threatening thrombotic emergency requiring immediate specialist management. This genotype is vanishingly rare given the allele frequency of approximately 0.002-0.003% in the global population.

Interactions

Protein S deficiency interacts with Factor V Leiden (rs6025) and the prothrombin G20210A variant (rs1799963) in a compound fashion. Individuals carrying both protein S deficiency and Factor V Leiden have VTE risks that approximate or exceed those of Factor V Leiden homozygotes. Similarly, compound heterozygosity for two different PROS1 mutations — for example Arg561Trp and a second PROS1 null allele — causes severe Type I protein S deficiency with protein S activity below 30% and a markedly elevated risk of early-onset, recurrent, and atypical-site thrombosis.

rs12478601

THADA THADA PCOS/Insulin Resistance

Strong Risk Factor

THADA — The Energy Regulator Linking PCOS to Insulin Resistance

At chromosome 2p21, a gene called THADA11 THADA
thyroid adenoma associated; encodes an armadillo repeat-containing protein that acts as a SERCA uncoupling factor
harbours one of the most robustly replicated genetic loci for polycystic ovary syndrome (PCOS). The rs12478601 C allele was first identified as a PCOS risk variant in 2011 and has since been confirmed across multiple ethnic groups, connecting the syndrome's hormonal features directly to a fundamental mechanism of cellular energy homeostasis.

The Mechanism

THADA functions as a SERCA uncoupling protein22 SERCA uncoupling protein
SERCA = sarco/endoplasmic reticulum Ca²⁺-ATPase, a pump that moves calcium from the cytoplasm into the ER lumen
. When THADA is present and functional, it binds SERCA and reduces the efficiency of calcium pumping — effectively "wasting" some of the energy as heat rather than storing calcium. When THADA activity is reduced (as variants at this locus may cause), SERCA runs unchecked, over-filling ER calcium stores. In pancreatic beta cells, excess ER calcium disrupts the finely tuned calcium signalling required for insulin secretion33 excess ER calcium disrupts the finely tuned calcium signalling required for insulin secretion
THADA knockdown raises ER Ca²⁺ stores and reduces beta-cell responsiveness to GLP-1 and arginine stimulation, suggesting lower effective beta-cell mass
. Impaired insulin secretion in the face of normal or elevated glucose — a hallmark of early insulin resistance — creates the metabolic backdrop on which PCOS develops.

The Evidence

The THADA locus at 2p21 was identified as a genome-wide significant PCOS susceptibility region by Chen et al. 201144 Chen et al. 2011
Genome-wide association study identifies susceptibility loci for PCOS on chromosome 2p16.3, 2p21 and 9q33.3. Nat Genet. 2011
in a Han Chinese cohort of >4,000 PCOS cases and 6,600 controls, with the lead SNP rs13429458 reaching OR 0.67, p=1.73×10⁻²³. rs12478601 is a companion tag SNP in the same THADA LD block, subsequently studied in multiple replication cohorts. The THADA association has been replicated in European-ancestry cohorts by Goodarzi et al. 201255 Goodarzi et al. 2012
Replication of DENND1A and THADA variants with PCOS in European cohorts. J Med Genet. 2012
and confirmed as part of the 14 susceptibility loci reported by the large-scale European meta-analysis by Day et al. 201866 Day et al. 2018
Large-scale GWAS meta-analysis of PCOS suggests shared genetic architecture for different diagnosis criteria. PLoS Genet. 2018
(10,074 PCOS cases, 103,164 controls).

The THADA-SERCA connection was established by Moraru et al. 201777 Moraru et al. 2017
THADA regulates the organismal balance between energy storage and heat production. Dev Cell. 2017
. THADA loss in Drosophila produced obesity, reduced thermogenesis, and elevated ER calcium stores — phenotypes rescued by simultaneously reducing SERCA activity. Human cell experiments confirmed the SERCA uncoupling function. At the clinical level, Tian et al. 202088 Tian et al. 2020
PCOS-GWAS susceptibility variants in THADA, INSR, TOX3, and DENND1A are associated with metabolic syndrome or insulin resistance in women with PCOS. Front Endocrinol. 2020
showed that in 2,082 Han Chinese PCOS women, the CC genotype at rs12478601 was associated with a measurably different metabolic syndrome rate, while Cui et al. 201399 Cui et al. 2013
Genotype-phenotype correlations of PCOS susceptibility SNPs identified by GWAS in a large cohort of Han Chinese women. Hum Reprod. 2013
found CC homozygotes had elevated LDL cholesterol, linking the risk genotype to downstream cardiovascular metabolic consequences.

Practical Actions

Women carrying one or two copies of the C allele face a meaningfully elevated PCOS risk through the insulin secretion pathway. The primary actionable implication is monitoring for features of insulin resistance (hyperinsulinaemia, elevated fasting insulin, impaired glucose tolerance) before full PCOS develops, since the THADA mechanism acts specifically through beta-cell dysfunction rather than primary androgen excess. Dietary strategies that reduce the beta-cell glucose load — particularly limiting refined carbohydrates that trigger large insulin pulses — address the specific mechanism this variant affects. Inositol supplementation (myo-inositol or D-chiro-inositol) has a growing evidence base as an insulin sensitiser specifically in PCOS, distinct from general metabolic interventions.

Interactions

rs12478601 in THADA operates in parallel with two other replicated PCOS loci in the database: rs2479106 and rs7852296 (both in DENND1A), which act through a different pathway — elevated androgen biosynthesis in theca cells. THADA variants tag the metabolic/insulin-secretion arm of PCOS pathogenesis, while DENND1A variants tag the androgenic arm. Women carrying risk alleles at both loci may face the most severe PCOS phenotype, combining impaired insulin secretion with elevated androgen production. This interaction is worth noting for clinicians interpreting combined genomic results.

PTPN22 rs12730735 — The Haplotype Refiner Behind the Autoimmune Switch

PTPN22 encodes lymphoid tyrosine phosphatase (LYP), the master negative regulator of T-cell and B-cell activation. While the R620W missense variant (rs2476601) is the most studied variant in this gene, the PTPN22 locus harbors additional variation that modifies autoimmune risk in ways R620W alone cannot capture. rs12730735 is an intronic variant that serves as a haplotype tag SNP11 haplotype tag SNP
Tag SNPs mark a specific chromosomal segment containing several variants in high linkage disequilibrium, allowing them to be used as proxies for the entire block's effects
, tracking a distinct chromosomal segment of the PTPN22 locus across diverse populations.

The Mechanism

rs12730735 lies within an intron of PTPN22, producing no amino acid change. Its biological effect operates through linkage disequilibrium: the C allele marks a specific chromosomal background within the PTPN22 locus that may carry regulatory variants affecting gene expression, splicing, or transcription factor binding in immune cells. PTPN22 expression is tightly regulated in lymphocytes22 PTPN22 expression is tightly regulated in lymphocytes
Expression studies show PTPN22 is predominantly expressed in lymphoid tissue; intronic variants can affect mRNA processing and expression levels
, and intronic variants in regulatory regions can alter these expression dynamics without changing the protein sequence. In East Asian populations, where R620W is essentially absent (allele frequency ~1%), rs12730735 provides an independent window into PTPN22-linked autoimmune susceptibility — a population that the R620W-centric research framework largely misses.

The Evidence

The most direct evidence for rs12730735 comes from a three-SNP PTPN22 analysis in rheumatoid arthritis33 a three-SNP PTPN22 analysis in rheumatoid arthritis
Bourgey et al. examined rs2476601, rs12730735, and rs11102685 together to model PTPN22 effect on RA; combining all three expanded the risk range from GRR 1–2.7 to 1–4.7 compared to using R620W alone
. When rs12730735 and rs11102685 were added to the R620W-based risk model, the genotypic risk ratio range expanded substantially — from a maximum of approximately 2.7-fold to 4.7-fold — uncovering high-risk subgroups among individuals who appeared identical on R620W status alone. Individuals in the highest genotypic class faced roughly four times the RA risk of the reference group (GRR > 3, confirmed by identity-by-descent excess in sib pairs).

In a Korean population study where R620W is non-polymorphic, rs12730735 minor allele (C) and a five-SNP haplotype (GGCTT) showed significant association with autoimmune thyroid disease44 rs12730735 minor allele (C) and a five-SNP haplotype (GGCTT) showed significant association with autoimmune thyroid disease
Lee et al. found p < 0.01 for the rs12730735 minor allele in AITD, specifically Hashimoto's thyroiditis; the R620W variant was absent from this population
. This study included 389 T1D patients, 212 AITD patients, and 225 controls. Notably, no association was found with type 1 diabetes in this cohort, suggesting rs12730735's primary signal relates to thyroid autoimmunity in populations where R620W does not confound the analysis. This independence makes it especially informative for the ~60% of the world's population with Asian ancestry.

For RA in European populations, a well-powered study of 4,460 RA cases and 4,481 controls55 well-powered study of 4,460 RA cases and 4,481 controls
Wan Taib et al. found no independent PTPN22 effect on RA beyond R620W in 4,460 European cases; rs12730735's RA association in this population is explained by its LD with R620W
found no evidence of RA risk beyond R620W. This suggests rs12730735's RA association in Europeans is largely mediated through co-inheritance with R620W haplotypes rather than an independent causal effect — but the Korean data indicates the variant's chromosomal segment has autoimmune relevance beyond R620W.

More recently, a Polish study of T1D patients66 Polish study of T1D patients
Cichocka et al. found CT genotype of rs12730735 associated with improved glycemic control in 277 T1D patients
found that the CT genotype of rs12730735 was associated with improved likelihood of achieving glycemic control — an intriguing finding that may reflect modulation of immune-mediated beta-cell destruction rates in T1D.

Practical Actions

For carriers of the C allele, the actionable picture is nuanced by ancestry. In European populations, rs12730735 generally co-occurs with R620W haplotypes when present on risk-associated backgrounds, meaning the main autoimmune monitoring guidance is similar to R620W carriers. In Asian populations, rs12730735 provides an independent signal for Hashimoto's thyroiditis risk — the most common autoimmune thyroid condition. Carriers should be attentive to thyroid symptoms and consider periodic thyroid function monitoring. Individuals with East Asian ancestry who carry the C allele and have a family history of thyroid autoimmunity are the clearest beneficiaries of proactive thyroid screening.

Interactions

The most important interaction is with rs2476601 (R620W). In Europeans, these two variants tend to co-segregate on the same haplotype; the three-SNP combination rs2476601 + rs12730735 + rs11102685 creates more granular risk stratification than R620W alone. When an individual carries the R620W risk allele alongside rs12730735-C, the combined genotypic risk ratio may reach 4-fold or higher — a substantially elevated risk not captured by R620W alone. This interaction is documented in the Bourgey 2007 analysis.

A proposed compound action should be considered for individuals who carry both rs2476601 risk allele (A/G or A/A) and rs12730735 C allele — especially for early and intensive autoimmune monitoring given the synergistic risk elevation.

SLC23A2 rs1279683 — Vitamin C Delivery, Glaucoma Risk, and the SVCT2 Transporter

Vitamin C cannot synthesize itself — it must be obtained from food, absorbed by the gut, and then actively transported into individual cells. That last step is controlled by SVCT211 SVCT2
Sodium-dependent Vitamin C Transporter 2, encoded by SLC23A2 — the primary transporter moving ascorbate from the bloodstream into metabolically demanding tissues such as the brain, adrenal glands, and eyes
. A common intronic variant at rs1279683 in SLC23A2 has been associated with meaningfully lower plasma vitamin C concentrations and a 67% increased risk of primary open-angle glaucoma in carriers of the G allele. Because SVCT2 is the dominant transporter for the retina, brain, and adrenal cortex — tissues that concentrate ascorbate to levels 10–100 times higher than plasma — this variant sits at the intersection of antioxidant defence, intraocular pressure regulation, and cognitive resilience.

The Mechanism

rs1279683 falls within an intron of SLC23A2 on chromosome 20 (GRCh38 position 5,002,446). SLC23A2 is transcribed from the minus strand. The variant's functional consequence is not characterized by in vitro assay, but intronic variants can alter pre-mRNA splicing efficiency22 pre-mRNA splicing efficiency
Intronic variants near splice regulatory elements or branch points can disrupt the recognition signals needed to join exons correctly, potentially reducing transporter protein output
, local enhancer/silencer activity, or mRNA stability. The measurable reduction in plasma ascorbate in G allele carriers implies that this intronic position does influence SVCT2 expression or function in a detectable way.

SVCT2 is most highly expressed in the adrenal gland (RPKM 40.7) and brain (RPKM 20.8). It is the sole active transporter moving ascorbate across the blood-brain barrier and into the retinal ganglion cells. In the aqueous humor of the eye, ascorbate is present at concentrations up to 20-fold higher than plasma — a gradient maintained by SVCT2 activity in ciliary epithelial cells. That gradient is protective against oxidative damage to the trabecular meshwork, which regulates intraocular pressure. Reduced SVCT2 efficiency compresses this protective gradient and may contribute to glaucoma pathogenesis.

The Evidence

The strongest direct evidence comes from a Spanish case-control study33 Spanish case-control study
Zanon-Moreno V et al. Association between a SLC23A2 gene variation, plasma vitamin C levels, and risk of glaucoma in a Mediterranean population. Mol Vis, 2011
of 150 POAG cases and 150 controls. GG homozygotes had plasma vitamin C of 9.0 ± 1.4 µg/ml in cases compared to 10.5 ± 1.6 µg/ml in non-GG genotypes — a statistically significant 14% reduction. Among controls, GG subjects showed 10.9 vs 12.1 µg/ml. The GG genotype conferred an age- and sex-adjusted odds ratio for POAG of 1.67 (95% CI 1.03–2.71, p = 0.010).

The glaucoma association was replicated across a larger evidence base. A meta-analysis of 108 case-control studies44 meta-analysis of 108 case-control studies
Chen M et al. Association of Gene Polymorphisms With Primary Open Angle Glaucoma: A Systematic Review and Meta-Analysis. Invest Ophthalmol Vis Sci, 2019
covering 35,389 POAG patients and 51,742 controls identified rs1279683 among 20 SNPs in 12 genes with consistent genome-wide significance.

A Japanese cohort study of 393 subjects examined vitamin C transporter variants and cognitive decline55 vitamin C transporter variants and cognitive decline
Hayashi K et al. Effects of functional variants of vitamin C transporter genes on apolipoprotein E E4-associated risk of cognitive decline: The Nakajima study. PLoS One, 2021
. Minor allele carriers (i.e., individuals with at least one A allele) did not show elevated cognitive risk. However, among the subset also carrying APOE4, carrying the minor allele (A) was protective: APOE4 carriers who were GG homozygotes had substantially higher cognitive decline risk (OR 2.02, 95% CI 1.05–3.87) compared to APOE4 carriers with at least one A allele. This finding suggests that adequate vitamin C transport is particularly important for neurological resilience when APOE4-mediated inflammation is present.

Large-scale blood cell GWAS data add a complementary dimension: a trans-ethnic analysis of 746,667 individuals66 a trans-ethnic analysis of 746,667 individuals
Chen MH et al. Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations. Cell, 2020
identified rs1279683 as a genome-wide significant signal for mean corpuscular volume (MCV, p = 3×10⁻²²) and mean corpuscular hemoglobin (MCH, p = 3×10⁻¹⁹). Vitamin C supports iron absorption and erythropoiesis — impaired SVCT2 function may reduce intracellular ascorbate available to erythroblasts, subtly affecting red cell development.

Practical Actions

The G allele at rs1279683 is the population-common allele, present at ~69% in Europeans, so GG homozygosity affects roughly 1 in 3 people of European descent. The vitamin C deficit is moderate — GG homozygotes in the Zanon-Moreno cohort averaged about 9–11 µg/ml versus 11–12 µg/ml for other genotypes. This is clinically meaningful at the margin: suboptimal plasma vitamin C is defined below ~28 µmol/L (approximately 5 µg/ml), and GG homozygotes cluster at the lower end of the adequate range. Consistently high dietary or supplemental vitamin C intake can offset the reduced transport efficiency by increasing the plasma substrate available for SVCT2.

Glaucoma management for G allele carriers does not replace standard monitoring but adds a modifiable risk factor: regular intraocular pressure measurement and a confirmed adequate vitamin C status. High-dose vitamin C (1,000–2,000 mg/day) has been shown in small trials to reduce intraocular pressure in some individuals; while not yet a clinical standard, it is a low-risk intervention for GG homozygotes with a family history or borderline IOP.

Interactions

rs1279683 acts in parallel with rs6053005, another intronic SLC23A2 variant associated with plasma vitamin C in the EPIC cohort. Both variants influence SVCT2-mediated tissue delivery; their combined effect has not been formally studied but is likely additive given their independent signals in the same gene. rs6133175, an SLC23A2 variant associated with preeclampsia and also independently predictive of plasma vitamin C, is a third independent signal at this locus.

The interaction with APOE4 (rs429358) is the most clinically significant: GG homozygotes who also carry APOE4 face a compounded risk — APOE4 increases neuroinflammation and oxidative burden in the brain while GG reduces the SVCT2 capacity to supply ascorbate for antioxidant defence. This combination warrants particular attention to optimizing vitamin C status. Propose for compound action: GG × APOE4 → augmented vitamin C supplementation and cognitive monitoring.