CD58 rs1414273 — The miR-548ac Switch That Defines the MS-Risk Haplotype
Every CD58 intronic variant associated with multiple sclerosis risk — rs2300747, rs12044852, rs1016140, rs1335532 — sits within the same tightly linked haplotype block in the first intron of CD58 on chromosome 1. Among them, rs1414273 holds a mechanistically unique position11 rs1414273 holds a mechanistically unique position
rs1414273 lies directly within the hairpin sequence of hsa-miR-548ac, the microRNA co-encoded from the same primary transcript as CD58 mRNA: it is the functional anchor through which the risk haplotype disrupts the balance between CD58 and an immunoregulatory microRNA.
CD58, also known as LFA-3 (Lymphocyte Function-Associated Antigen 3)22 LFA-3 (Lymphocyte Function-Associated Antigen 3)
LFA-3 is a glycoprotein expressed on antigen-presenting cells and non-immune tissues; it binds CD2 on T cells to stabilise the immune synapse and transmit co-stimulatory signals that licence regulatory T cell induction, governs whether a T-cell encounter with an antigen-presenting cell produces immune tolerance or inflammatory activation. This signalling axis sits at the centre of multiple sclerosis genetics: the CD58 locus is one of the most robustly replicated non-HLA MS susceptibility regions.
The Mechanism
The miR-548ac gene is embedded within CD58 intron 1. Like all microRNAs, it begins as part of a longer primary transcript (pri-miRNA) that is cropped by the nuclear Drosha-DGCR8 endonuclease complex33 Drosha-DGCR8 endonuclease complex
Drosha cleaves the pri-miRNA at the base of the hairpin stem-loop to release the 60–70 nt precursor (pre-miRNA), which is then exported to the cytoplasm and processed by Dicer into the mature ~22 nt miRNA. Because CD58 mRNA and miR-548ac are transcribed from the same genomic locus, Drosha cleavage efficiency has a reciprocal effect: more efficient cleavage favours miR-548ac production at the expense of full-length CD58 mRNA reaching the cytoplasm for translation.
rs1414273 sits at the precise base of the miR-548ac stem-loop. The reference allele is C on the plus strand — corresponding to G on the coding (minus) strand — which pairs with a uridine in the opposite strand of the hairpin to form a G-U wobble base pair. The alternate T allele (coding-strand A) forms an A-U Watson-Crick base pair, which is less efficiently recognised by Drosha. The consequence:
- C allele (risk): G-U wobble → enhanced Drosha recognition → more miR-548ac, less CD58 mRNA
- T allele (protective): A-U Watson-Crick → reduced Drosha efficiency → less miR-548ac, more CD58 mRNA
In cell-culture experiments, the G-containing construct (risk allele) produced 1.5-fold more mature miR-548ac at 24 hours and 3.4-fold more at 48 hours compared to the A-containing construct, after normalisation to precursor RNA. In population-level eQTL analyses spanning HapMap, Geuvadis, and MS patient cohorts, risk allele carriers showed significantly lower CD58 transcript levels and significantly higher miR-548ac levels44 risk allele carriers showed significantly lower CD58 transcript levels and significantly higher miR-548ac levels
The paradoxical inverse relationship between CD58 mRNA and miR-548ac co-encoded from the same locus is explained by Drosha-mediated uncoupling: cleavage of the stem-loop interrupts full-length mRNA synthesis while liberating the microRNA hairpin.
miR-548ac targets validated in this mechanistic framework include SDC4 (syndecan-4)55 SDC4 (syndecan-4)
SDC4 regulates T-cell motility and heparan sulfate proteoglycan signalling, SEL1L66 SEL1L
SEL1L participates in endoplasmic reticulum-associated protein degradation (ERAD) and proteostasis under inflammatory stress, and TNFAIP3 (A20)77 TNFAIP3 (A20)
TNFAIP3/A20 is a master ubiquitin-editing enzyme that terminates NF-κB signalling; miR-548ac-mediated suppression of A20 would sustain inflammatory NF-κB activity. A broader computational screen identified 333 predicted miR-548ac targets enriched in cytokine signalling, MAPK pathways, and protein folding — consistent with a broadly immunomodulatory role for this microRNA.
The Evidence
The functional significance of rs1414273 was first proposed by Galarza-Munoz et al. 201588 Galarza-Munoz et al. 2015
1000 Genomes sequencing data used to identify rs1414273 as the only SNP at the base of the miR-548ac stem-loop in strong LD with the MS-associated haplotype, who recognised that its position within the Drosha recognition element made it a strong candidate for the causal variant within the CD58 MS locus.
The full mechanistic and eQTL evidence was published by Hecker et al. 2019 PLoS Genetics99 Hecker et al. 2019 PLoS Genetics
Expression QTL analysis across >1,000 subjects from HapMap and Geuvadis cohorts combined with in vitro Drosha cleavage experiments and real-time PCR of 32 MS patients. This study demonstrated both the population-level expression effects and the molecular mechanism in a single unified study, making rs1414273 the best-characterised functional variant at the CD58 MS locus.
The variant is in near-complete LD (r² ≈ 1, D' ≈ 1) with rs1335532, the GWAS lead SNP for the CD58 MS locus in European cohorts. Correlated alleles on the forward strand are: rs1414273 C = rs1335532 A (both risk alleles); rs1414273 T = rs1335532 C (both protective alleles). This tight LD means rs1414273 and rs1335532 are functionally interchangeable for GWAS-level MS association signals, but rs1414273 uniquely maps the molecular effect to Drosha recognition.
Population genetics provide an instructive layer. The C (risk) allele is the major allele in Europeans (~87%) — matching the high-frequency risk pattern seen throughout the CD58 locus (rs12044852, rs2300747). In contrast, East Asian and African populations carry the T (protective) allele at ~59% and ~51% respectively. This mirrors the population-level MS burden: the disease is considerably rarer in East Asian populations where the protective T allele predominates.
A pilot replication study in a Malaysian/Kuwaiti cohort found CD58 rs1414273 to be significantly associated with MS in an Arab population (p = 0.00007, OR 2.2, 95% CI 1.5–3.2 in exome analysis), though a separate genotyping-only Kuwaiti replication did not reach significance, likely reflecting population-specific LD structure between rs1414273 and the GWAS signals.
Practical Actions
Since the C allele is the major allele in Europeans (~87%), CC homozygotes represent the largest fraction of the population — and carry the highest miR-548ac burden and lowest CD58 expression at this locus. There is no supplement that directly inhibits miR-548ac or compensates for reduced CD58, but vitamin D is the best-characterised modifiable factor that supports FoxP3+ regulatory T cell function1010 FoxP3+ regulatory T cell function
FoxP3 is the master transcription factor for Tregs; vitamin D drives FoxP3 expression via VDR binding at the FoxP3 promoter, providing a parallel Treg-supporting input to the CD58 co-stimulatory signal through mechanisms independent of CD58.
Interactions
rs1414273 is in near-complete LD with rs1335532 and is strongly correlated with rs12044852 (r² = 0.929 between rs12044852 and rs2300747, with rs1414273 in the same haplotype). The four CD58 intronic variants — rs2300747, rs12044852, rs1016140, rs1414273 — collectively define the CD58 MS-risk haplotype. The rs1016140 G allele has an independent signal for NMO risk and autoimmune thyroid disease via a mechanistically distinct direction, illustrating that the intronic haplotype block contains variants with both shared and divergent functional effects.
miR-548ac's validated target TNFAIP3 (A20) connects this variant to the NF-κB inflammatory pathway, creating a potential interaction with cytokine-pathway SNPs (e.g. TNFRSF1A rs1800693 in the same immune-gut category).
VDR BsmI — How Your Cells Respond to Vitamin D
The vitamin D receptor (VDR) is a nuclear receptor 11 A nuclear receptor is a protein that binds hormones or vitamins inside the cell and directly regulates gene expression that mediates the biological effects of vitamin D throughout your body. When active vitamin D (calcitriol) 22 Calcitriol (1,25-dihydroxyvitamin D) is the hormonally active form of vitamin D binds to VDR, it triggers gene expression changes that affect calcium absorption, immune function, cell growth, and hundreds of other processes. VDR is expressed in nearly every tissue in the body, which is why vitamin D affects so many aspects of health.
The Mechanism
The BsmI variant (rs1544410) is located in an intronic region of the VDR gene. While it does not directly change the protein sequence, it is in linkage disequilibrium 33 Linkage disequilibrium: nearby genetic variants that are inherited together more often than expected by chance with functional variants that affect VDR mRNA stability and expression levels. The T allele is associated with reduced VDR expression, meaning your cells produce fewer vitamin D receptors and are therefore less responsive to circulating vitamin D. The variant frequency varies dramatically by ancestry — 40% in Europeans but only 6% in East Asians.
The Evidence
A meta-analysis of 26 studies44 meta-analysis of 26 studies
Tao S et al. VDR BsmI polymorphism and osteoporosis risk, 2012 and a larger 42-study meta-analysis55 larger 42-study meta-analysis
Zhao L et al. VDR BsmI and osteoporosis in postmenopausal women, 2020
found that VDR BsmI variants are associated with osteoporosis susceptibility
in Caucasians (OR 0.70 for bb vs BB), bone mineral density, and calcium
absorption efficiency. The associations are strongest in populations with lower
baseline vitamin D levels. Additional research has linked VDR variants to
immune function, autoimmune disease risk, and cancer susceptibility, though
these associations are more complex and context-dependent.
The Vitamin D Optimization Challenge
VDR variants create a situation where standard blood levels of vitamin D may not produce standard biological effects. If your cells have fewer vitamin D receptors, you may need higher circulating vitamin D levels to achieve the same cellular response as someone with normal VDR expression. This is why some people with "adequate" blood levels still seem to benefit from higher vitamin D intake.
Practical Implications
If you carry the T allele, maintaining vitamin D levels in the optimal range (30-50 ng/mL) is important, and you may benefit from aiming toward the higher end of that range. Regular testing (1-2 times per year) helps you calibrate your supplementation. Vitamin D3 is preferred over D2, and taking it with a fat-containing meal improves absorption.
Interactions
VDR interacts with CYP2R1 (rs10741657) — if both vitamin D activation and receptor sensitivity are impaired, the combined "double hit" significantly impacts vitamin D status.
The Adiponectin Paradox — When More Isn't Always Better
Adiponectin is your body's master metabolic regulator11 regulator
a hormone secreted by fat tissue that enhances insulin sensitivity, reduces inflammation, and protects against metabolic disease, and the ADIPOQ gene controls how much of it you produce. The rs17300539 variant sits in the gene's promoter region — the control switch that determines transcription activity22 transcription activity
how actively the gene is read and translated into protein. What makes this variant fascinating is its paradoxical effects: the A allele cranks up adiponectin production, yet doesn't always deliver the metabolic protection you'd expect.
Normally, higher adiponectin is protective — it improves insulin sensitivity, lowers inflammation, reduces cardiovascular risk, and guards against type 2 diabetes. People with obesity and metabolic syndrome typically have low adiponectin levels33 low adiponectin levels
adiponectin secretion is impaired in obesity, creating a vicious cycle of worsening insulin resistance, which contributes to their disease. Yet your genotype at rs17300539 introduces a twist: some people produce more adiponectin but still face elevated metabolic risk.
The Mechanism
The rs17300539 SNP is a G-to-A substitution at position -11391 in the ADIPOQ promoter region. In vitro studies44 In vitro studies
laboratory experiments using cell cultures demonstrate that the A allele significantly increases transcriptional activity compared to the G allele, driving higher adiponectin production. The variant likely alters transcription factor binding55 transcription factor binding
proteins that attach to DNA and regulate gene expression at this promoter site, though the exact factors involved haven't been fully mapped.
Adiponectin circulates in your blood in three forms: low molecular weight (LMW) trimers, medium molecular weight (MMW) hexamers, and high molecular weight (HMW) multimers66 low molecular weight (LMW) trimers, medium molecular weight (MMW) hexamers, and high molecular weight (HMW) multimers. The HMW form is the most biologically active — it's the one that enhances insulin sensitivity77 enhances insulin sensitivity
stimulates AMPK activation in muscle and liver, increasing glucose uptake and fatty acid oxidation and delivers cardiovascular protection. Some evidence suggests that rs17300539 may influence the ratio of HMW to total adiponectin88 ratio of HMW to total adiponectin, which could explain why total adiponectin levels don't always predict metabolic outcomes in carriers.
The Evidence
The Framingham Offspring Study99 Framingham Offspring Study
a landmark cardiovascular epidemiology study following multiple generations genotyped 2,543 participants and found that the A allele at rs17300539 showed the strongest association with higher adiponectin levels (P = 2.6 × 10⁻⁸). Each A allele added roughly 1.6 μg/mL to circulating adiponectin. This finding has been replicated across multiple populations — European, Asian, and Latino cohorts all show the same pattern.
But here's the paradox: a 2009 study in obese children1010 a 2009 study in obese children
1,210 Greek children aged 9-13, both obese and non-obese found that A-allele carriers (GA+AA) had higher adiponectin levels but also higher BMI (B = 0.97, P = 0.015) and a 35% increased odds of obesity (OR = 1.35, 95% CI 1.06-1.85). Before adjusting for obesity status, they showed higher fasting insulin and higher HOMA-IR (a measure of insulin resistance). The researchers concluded that "the rs17300539-A variant, though consistently associated with higher adiponectin levels, does not exert any appreciable protective metabolic effect in children."
In adults, the story differs by baseline metabolic health. A 2023 study in 329 obese Caucasian adults1111 A 2023 study in 329 obese Caucasian adults
Spanish cohort with mean BMI 47.8 kg/m² found that GG homozygotes had significantly higher rates of metabolic syndrome (86% vs. 73.9%, P < 0.05), hypertriglyceridemia, hyperglycemia, and insulin resistance (HOMA-IR 7.49 vs. 4.62) compared to A-allele carriers. GG carriers also had lower adiponectin levels (4.27 vs. 6.36 μg/mL). Logistic regression confirmed that the GG genotype independently increased metabolic syndrome risk (OR = 2.52, 95% CI 1.04-6.10) even after adjusting for age, sex, weight, and dietary intake.
The variant also shows strong association with polycystic ovary syndrome (PCOS)1212 strong association with polycystic ovary syndrome (PCOS) in Chinese populations — a family-based transmission disequilibrium test in 197 PCOS families confirmed overtransmission of the risk allele. PCOS is fundamentally a condition of insulin resistance and hyperandrogenism, often accompanied by low adiponectin.
A meta-analysis of 35 studies1313 A meta-analysis of 35 studies
nearly 29,000 participants across multiple ethnicities linked rs17300539 to coronary artery disease (CAD) risk, though effect sizes varied by population and the direction wasn't always consistent — likely reflecting the complex interplay between adiponectin levels, HMW ratio, and other metabolic factors.
Practical Actions
The clearest clinical implication emerges from the bariatric surgery literature1414 bariatric surgery literature: A-allele carriers show better lipid profile improvements after surgery. In 60 extremely obese individuals followed for 32 months post-surgery, those with the A-C haplotype (combining rs17300539-A with rs266729-C) had greater reductions in LDL cholesterol. This suggests that in the context of major metabolic intervention — whether bariatric surgery or intensive lifestyle modification — the A allele's adiponectin-boosting effect finally translates into benefit.
The gene-diet interaction studies1515 gene-diet interaction studies are particularly relevant. In the RISCK study, rs17300539 genotype interacted significantly with dietary fat composition to determine adiponectin levels. Another study in the GOLDN cohort1616 Another study in the GOLDN cohort found that the association between the -11391A allele and lower BMI was modified by monounsaturated fatty acid (MUFA) intake — A-allele carriers who consumed higher MUFA had the lowest BMI and obesity risk.
Fish oil supplementation may be particularly relevant: omega-3 fatty acids activate PPARγ1717 omega-3 fatty acids activate PPARγ, which upregulates adiponectin expression, and one study found that ADIPOQ genotype modified the response to fish oil supplementation in older individuals.
Interactions
The ADIPOQ gene sits at the intersection of several metabolic pathways. Adiponectin signals through two receptors — AdipoR11818 AdipoR1
predominantly expressed in skeletal muscle, activates AMPK pathways and AdipoR21919 AdipoR2
predominantly in liver, activates PPARα signaling. The downstream effects include increased fatty acid oxidation, reduced hepatic glucose production, and improved insulin sensitivity.
Three other common ADIPOQ SNPs show linkage disequilibrium with rs17300539: rs266729 (-11377C>G, also in the promoter, r² = 0.80 with rs17300539), rs2241766 (+45T>G in exon 2, also called Gly15Gly), and rs1501299 (+276G>T in intron 2). These variants may compound or modify effects, particularly regarding the HMW adiponectin ratio. Haplotype analysis sometimes reveals stronger associations than single SNPs alone.
There's emerging evidence for interaction with TCF7L2 variants2020 TCF7L2 variants, the strongest type 2 diabetes risk gene. TCF7L2 regulates adipocyte development and function, and deletion of TCF7L2 in adipocytes impairs glucose tolerance and alters lipid metabolism. The combination of ADIPOQ and TCF7L2 risk variants may identify individuals who benefit most from dietary fat modification.
Finally, the obesity paradox deserves emphasis: if you're lean and metabolically healthy, higher adiponectin from the A allele is likely beneficial. But if you're already obese or insulin-resistant, the A allele may signal a compensatory response — your body is pumping out more adiponectin to counteract metabolic dysfunction, but it's not enough to overcome the underlying problem. In that scenario, the GG genotype's association with lower adiponectin may simply reflect better baseline metabolic health.
FABP2 — Fat Absorption Efficiency
FABP2 (Fatty Acid Binding Protein 2) is expressed in intestinal cells and is responsible for intracellular transport of dietary fatty acids11 Inside enterocytes (intestinal absorptive cells), FABP2 shuttles fatty acids from the cell membrane to the endoplasmic reticulum for processing.
The Mechanism
The Ala54Thr variant (rs1799883) is a missense mutation in exon 2 of FABP2,
where an adenine replaces guanine at the DNA level, substituting alanine with
threonine at position 54 of the protein (p.Ala54Thr). Baier et al.22 Baier et al.
Baier et al. An amino acid substitution in the human intestinal fatty acid binding protein is associated with increased fatty acid binding, increased fat oxidation, and insulin resistance. J Biol Chem, 1995 demonstrated that the
threonine-containing protein has a 2-fold greater affinity for long-chain
fatty acids than the alanine-containing protein, leading to more efficient
fat absorption from the intestine.
The Evidence
The original discovery by Baier et al.33 original discovery by Baier et al.
Baier et al. J Biol Chem, 1995 in Pima Indians showed that
Thr54 carriers had higher fasting insulin, lower insulin-stimulated glucose
uptake, and higher fat oxidation rates. The threonine variant increases the
protein's affinity for long-chain fatty acids by approximately 2-fold.
Carriers of the Thr allele absorb more calories from fat44 Studies estimate Thr carriers may absorb roughly 20-30% more long-chain fatty acids per meal than Ala/Ala individuals, which can contribute to weight gain when fat intake is high.
A meta-analysis by Zhao et al.55 meta-analysis by Zhao et al.
Zhao et al. Association between FABP2 Ala54Thr polymorphisms and T2DM risk: a HuGE review and meta-analysis. Lipids Health Dis, 2014 found significant
associations with type 2 diabetes in Asian populations (OR 1.19, 95% CI
1.05-1.36) but not in Caucasians. The evidence for obesity association
is mixed, with some meta-analyses finding no significant effect on BMI.
Practical Implications
The Thr allele is common across all populations (24-33% frequency), with slightly higher frequency in South Asian and East Asian groups. The practical significance is moderate — this variant matters most when combined with high dietary fat intake, where increased absorption efficiency can contribute to excess calorie intake and insulin resistance.
Interactions
FABP2 Ala54Thr interacts with total dietary fat intake — the variant's metabolic effects are more pronounced on high-fat diets. If you also carry TCF7L2 risk alleles (rs7903146), moderating fat intake becomes doubly important.
TPMT*2 — The Original Thiopurine Deficiency Allele, Independent of the *3 Cluster
TPMT (thiopurine S-methyltransferase) is the enzyme that inactivates thiopurine drugs11 inactivates thiopurine drugs
Azathioprine, 6-mercaptopurine, and thioguanine — used for inflammatory bowel disease, autoimmune conditions, organ transplantation, and childhood leukemia maintenance by methylating them to inert metabolites. When TPMT activity is reduced or absent, these drugs are diverted into a pathway that generates highly toxic thioguanine nucleotides, which incorporate into DNA and cause life-threatening bone marrow suppression at standard doses. TPMT*2 (rs1800462) was the first TPMT deficiency allele ever identified — described by Krynetski and colleagues in 1995 — and it remains one of the four clinically relevant no-function star alleles tested before thiopurine prescription. Unlike the *3 cluster (*3A, *3B, *3C), which dominates TPMT deficiency in Europeans, *2 is a completely independent functional allele on its own chromosomal background and is the rarer of the two European no-function variant classes.
The Mechanism
TPMT*2 is a single-nucleotide substitution at position 238 of the TPMT coding sequence (c.238G>C) that changes alanine 80 to proline22 alanine 80 to proline
p.Ala80Pro — proline introduces a rigid kink into the protein backbone, disrupting local secondary structure near the SAM-binding pocket in the methyltransferase domain. TPMT is on the minus strand of chromosome 6, so the coding-strand c.238G>C mutation corresponds to a plus-strand C→G substitution at chr6:18143724 (GRCh38). Functional studies by Tai and colleagues33 Tai and colleagues
Tai HL et al. Enhanced proteasomal degradation of mutant TPMT (TPMT*3A and TPMT*2) in humans. PNAS 1997 showed that the Ala80Pro substitution destabilizes the folded protein, marks it for ubiquitin-dependent degradation44 ubiquitin-dependent degradation
The mutant enzyme has a half-life roughly 15-fold shorter than wild-type TPMT, with near-complete turnover via the 26S proteasome, and leaves cells with essentially no functional TPMT from the 2 allele. In vitro, TPMT*2 retains only about 1% of wild-type catalytic activity — roughly a **100-fold reduction* — making it functionally equivalent to the *3A and *3C no-function alleles for clinical purposes.
The Evidence
TPMT is the most thoroughly characterised pharmacogene in clinical practice. The CPIC thiopurine dosing guideline55 CPIC thiopurine dosing guideline
Clinical Pharmacogenetics Implementation Consortium — Level A evidence, the highest tier for clinical implementation has been published since 2011 and updated multiple times since. Every CPIC version classifies TPMT*2 as a no-function allele with activity score 0, functionally identical to *3A, *3B, and *3C for dosing purposes. The original 2011 guideline66 original 2011 guideline
Relling MV et al. Clinical Pharmacogenetics Implementation Consortium guidelines for thiopurine methyltransferase genotype and thiopurine dosing. Clin Pharmacol Ther 2011 named *2, *3A, *3B, and *3C as the four variant alleles that should be tested in routine clinical practice — a panel that remains standard. The 2018 update77 2018 update
Relling MV et al. CPIC Guideline for Thiopurine Dosing Based on TPMT and NUDT15 Genotypes: 2018 Update. Clin Pharmacol Ther 2019 extended the framework to NUDT15 (the dominant thiopurine safety gene in East Asians) and refined the dose-reduction recommendations. The original molecular identification of TPMT*288 molecular identification of TPMT*2
Krynetski EY et al. A single point mutation leading to loss of catalytic activity in human thiopurine S-methyltransferase. Proc Natl Acad Sci USA 1995 came from cloning and sequencing TPMT cDNA from a TPMT-deficient patient who had experienced severe myelosuppression on 6-mercaptopurine. Population studies99 Population studies
Otterness DM et al. Human thiopurine methyltransferase pharmacogenetics: gene sequence polymorphisms. Clin Pharmacol Ther 1997 place the TPMT*2 allele frequency at roughly 0.2-0.5% in Europeans, making it 5-10 times less common than the *3A haplotype but still responsible for a measurable fraction of TPMT-deficient individuals of European ancestry. TPMT*2 is essentially absent in East Asian and South Asian populations, where NUDT15 variants (rs116855232 in particular) are the dominant cause of thiopurine toxicity. The FDA includes TPMT status in its pharmacogenomic biomarker table1010 pharmacogenomic biomarker table for azathioprine, mercaptopurine, and thioguanine, with package inserts explicitly recommending pre-treatment genotyping.
Practical Implications
If you carry one or two copies of TPMT*2, you are at risk for severe thiopurine-induced bone marrow suppression at standard doses — the same clinical risk carried by people with TPMT*3A/3B/*3C variants. European populations carry the *2 variant at roughly 0.2-0.5% allele frequency, meaning about 1 in 200-400 people of European ancestry is heterozygous (one functional TPMT allele from the *2), and homozygous *2/*2 is extremely rare. The CPIC 2025 guideline recommends starting at 30-80% of the standard dose for intermediate metabolizers (one no-function allele of any type) and 10% of the standard dose or an alternative medication for poor metabolizers (two no-function alleles — whether *2/*2, *2/*3A, *2/*3C, or any other combination). Because TPMT deficiency classification depends on **all variant alleles in combination*, anyone found to carry *2 must also be genotyped for *3B (rs1800460), *3C (rs1142345), and NUDT15 (rs116855232) to determine the complete metabolizer phenotype. Missing a second variant on the other chromosome is the classic way that a compound heterozygous poor metabolizer gets labeled as an intermediate metabolizer and receives a dangerous dose.
Interactions
The most clinically important interaction for TPMT*2 is with other TPMT no-function alleles in trans (on the opposite chromosome). TPMT*2 is on a completely different chromosomal background from the 3 cluster — it is NOT part of the *3A haplotype — which means someone carrying *2 plus *3A in trans, or *2 plus *3B alone, or *2 plus *3C alone, is a **compound heterozygous poor metabolizer* with essentially zero TPMT activity and the same dosing requirements as a *3A/*3A homozygote (10% of standard dose or thiopurine avoidance entirely). Clinical labs distinguish *2 carriers by direct genotyping at rs18004621111 direct genotyping at rs1800462
Most commercial TPMT panels test *2, *3A, *3B, and *3C simultaneously — but some older panels or raw-data imports may miss *2 entirely, which is why this rsid must be part of any TPMT pharmacogenomic pipeline. If your genome file shows rs1800462(G) combined with rs1142345(C) or rs1800460(T), your prescriber should treat you as a poor metabolizer until haplotype phasing confirms otherwise.
A second critical interaction is with NUDT15 (rs116855232), the other thiopurine safety gene. NUDT15 loss-of-function variants act on a downstream step of thiopurine metabolism, and individuals carrying variants in both TPMT and NUDT15 require larger dose reductions than predicted by either gene alone. Both genes should always be checked together before thiopurine prescription. Finally, co-administration of allopurinol1212 allopurinol
Xanthine oxidase inhibitor used for gout; blocks an alternate thiopurine inactivation pathway or febuxostat with thiopurines creates a double-blockade that is particularly dangerous in TPMT variant carriers — the FDA label instructs reducing azathioprine to 25% of the standard dose when combined with allopurinol, and that reduction must be applied on top of any TPMT-based dose reduction.
The Muscle Growth Brake — How K153R Affects Your Training Response
Myostatin is one of the most powerful negative regulators of skeletal muscle growth
in the human body. Secreted by muscle cells, it acts as a biological brake, preventing
muscles from growing too large. The K153R polymorphism (Lys153Arg) sits within the
mature active peptide11 mature active peptide
the bioactive portion of myostatin after proteolytic
processing of the myostatin protein,
where it can influence both how the protein is processed and how effectively it
binds to its receptor, ActRIIB22 ActRIIB
activin type II receptor B, the primary receptor
through which myostatin signals.
The rare R (arginine) variant appears to reduce myostatin's inhibitory effect, effectively loosening the brake on muscle growth. This has made it a variant of intense interest in sports genetics research.
The Mechanism
Myostatin is synthesized as a latent precursor protein that undergoes proteolytic
processing to become the mature, bioactive peptide. The K153R substitution occurs
at amino acid position 153 in this mature region. The replacement of lysine (K) with
arginine (R) — both positively charged amino acids but with different side chain
properties — may affect either the proteolytic processing of myostatin or its
binding affinity to the ActRIIB receptor33 may affect either the proteolytic processing of myostatin or its
binding affinity to the ActRIIB receptor
Lys and Arg have similar charge but
different side chain structures that could alter protein-protein
interactions.
When myostatin binds to ActRIIB on muscle cells, it activates intracellular
signaling cascades that inhibit both myoblast (muscle precursor cell) proliferation
and differentiation, ultimately limiting muscle mass accumulation. Any variant that
reduces this signaling — whether through altered processing, reduced receptor
binding, or decreased protein stability — would be expected to permit greater
muscle growth in response to mechanical loading (resistance training) or
muscle-building stimuli44 muscle-building stimuli
anabolic signals like IGF-1, testosterone, and
mechanical tension from exercise.
The Evidence
The most comprehensive examination of this variant comes from a 2022 meta-analysis55 2022 meta-analysis
Kruszewski & Aksenov. Association of Myostatin Gene Polymorphisms with Strength
and Muscle Mass in Athletes. Genes, 2022
that analyzed 71 research articles on MSTN polymorphisms. The meta-analysis included
4 studies with 773 athletes and 357 controls across 5 ethnic groups. The key finding:
strength-oriented athletes had a significantly higher frequency of the R variant
compared to controls (OR = 2.02, p = 0.05). Among athletes, those carrying the R
variant showed greater muscle strength and mass gains from power-oriented training
compared to KK carriers.
However, the picture is more nuanced than "R = better muscles." A 2011 study of
281 young non-athletic men66 2011 study of
281 young non-athletic men
Santiago et al. The K153R Polymorphism in the Myostatin
Gene and Muscle Power Phenotypes in Young, Non-Athletic Men. PLoS One,
2011 found that R allele carriers
actually performed worse on vertical jump tests, showing decreased peak power
production during explosive movements. This suggests the R variant may offer
advantages specifically in the context of chronic resistance training adaptation,
but not necessarily in baseline explosive power in untrained individuals.
The strongest evidence for functional impact comes from a training intervention
study in Han Chinese men77 training intervention
study in Han Chinese men
Wang et al. The A55T and K153R polymorphisms of MSTN
gene are associated with strength training-induced muscle hypertrophy. J Sports Sci,
2014. Among 94 previously untrained men
who completed an 8-week strength training program, those with the KR genotype showed
significantly greater muscle thickness gains: +0.30 cm in biceps and +0.42 cm in
quadriceps compared to KK genotype carriers (p < 0.01 for both). This represents
approximately 40-50% greater hypertrophy response from identical training stimulus.
The evidence level is moderate rather than strong due to the rarity of the R allele (3-4% in Europeans, making RR homozygotes extraordinarily rare at <1%) and some contradictory findings across studies. The meta-analysis showed moderate heterogeneity (I² = 33%), suggesting population-specific effects or gene-environment interactions.
Practical Actions
For KR or RR carriers, the variant suggests enhanced potential for muscle hypertrophy in response to progressive resistance training. This doesn't mean you'll automatically build more muscle — training, nutrition, and recovery remain primary determinants — but it suggests your ceiling for muscle growth may be higher than average when these factors are optimized.
Optimal training for leveraging this genetic advantage involves progressive overload with compound movements (squats, deadlifts, bench press, overhead press, rows), training each muscle group 2-3 times per week with sufficient volume (15-25 sets per muscle per week), and ensuring adequate protein intake (1.6-2.2 g/kg body weight daily). Recovery between sessions is crucial — the variant affects adaptation capacity, but adaptation still requires rest.
For AA (KK) carriers, this is the normal, wild-type genotype present in ~94% of Europeans. Myostatin regulation is functioning as designed. You have standard muscle growth potential, which is still substantial when training and nutrition are optimized. The absence of the R variant doesn't limit you to below-average muscle growth — it simply means you lack a rare genetic advantage.
Interestingly, a 2020 Mexican study88 2020 Mexican study
Castro-Rodríguez et al. The Myostatin
rs1805086 variant is associated with obesity in Mexican adults. Gene,
2020 found the R allele associated with
obesity independently of metabolic risk factors, suggesting that reduced myostatin
activity may have trade-offs beyond the muscle compartment. This reinforces that no
variant is universally "good" — context matters.
Interactions
The K153R variant shows documented interaction with another MSTN polymorphism, A55T (rs1805065). The same Han Chinese training study found that individuals carrying variant alleles of both polymorphisms showed the greatest training-induced muscle hypertrophy. The A55T variant is located in exon 1 of myostatin, also in the mature peptide region. These two variants may have additive or synergistic effects on reducing myostatin's inhibitory function.
There is theoretical but less well-documented interaction with ACTN3 R577X
(rs1815739), the "sprinter gene" that determines fast-twitch muscle fiber
composition. Since myostatin preferentially affects fast-twitch (type II) muscle
fibers, and ACTN3 determines the presence of alpha-actinin-3 protein exclusively in
type II fibers, carriers of both the MSTN R allele and the ACTN3 RR genotype might
show enhanced power and strength potential. However, studies specifically examining
this interaction have shown null results for combined effects on
longevity99 null results for combined effects on
longevity
Hirose et al. Muscle-Related Polymorphisms (MSTN rs1805086 and ACTN3
rs1815739) Are Not Associated with Exceptional Longevity. PLoS One,
2016, suggesting the interaction
may be context-dependent.
Nutrition-gene interaction is worth considering: adequate protein intake becomes even more critical for KR/RR carriers to realize the hypertrophic potential. The enhanced capacity for muscle protein synthesis means substrate availability (dietary protein) becomes rate-limiting faster than in TT carriers.
CRY1Δ11 — The Night-Owl Gene That Runs Your Clock 30 Minutes Slow
Inside every cell of your body ticks a molecular clock, cycling with almost
perfect 24-hour precision. The CRY1 protein is one of its master regulators
— a transcriptional brake that keeps the CLOCK:BMAL1 activator complex from
running perpetually. The CRY1Δ11 variant (c.1657+3A>C on the coding strand,
NC_000012.12:g.106992962T>G11 NC_000012.12:g.106992962T>G
genomic HGVS notation, GRCh38 plus strand)
disrupts a splice site at the boundary of intron 11, causing the entire exon
11 to be skipped during mRNA processing. The result is a CRY1 protein missing
24 amino acids from its C-terminal tail — and it is more powerful, not less.
This gain-of-function makes the mutant CRY1 cling more tightly to CLOCK and
BMAL1, prolonging transcriptional inhibition and stretching the molecular
clock period from 24 hours to roughly 24.5 hours. Half an hour per cycle
compounds: carriers' sleep timing drifts 2–2.5 hours later than their
circadian phase actually warrants22 carriers' sleep timing drifts 2–2.5 hours later than their
circadian phase actually warrants
Patke et al. Cell 2017,
producing the signature symptom of Delayed Sleep Phase Disorder.
The Mechanism
The 5' splice site of intron 11 contains the sequence the spliceosome uses
to remove intron 11 and join exon 11 to exon 12. The c.1657+3A>C transversion
weakens this recognition sequence enough that the spliceosome skips exon 11
entirely. The resulting mRNA encodes a CRY1 protein with an in-frame 24-amino
acid deletion in the C-terminal tail domain33 tail domain
The C-terminal tail modulates
CRY1's interaction surface with CLOCK and BMAL1.
Counterintuitively, this deletion enhances rather than impairs CRY1's function.
The mutant protein localizes to the nucleus more readily than wild-type CRY1,
occupies CLOCK:BMAL1 binding sites on chromatin for longer, and suppresses
transcription of circadian target genes — including Per1, Per2, and Dbp —
more persistently. Chromatin immunoprecipitation studies showed the mutant CRY1
displaced CLOCK and BMAL1 from their target gene promoters44 Chromatin immunoprecipitation studies showed the mutant CRY1
displaced CLOCK and BMAL1 from their target gene promoters
Consistent with
a dominant gain-of-function mechanism.
The net effect is a circadian period that runs slow, anchoring the person's
internal clock later and later relative to the external light-dark cycle.
The Evidence
The founding study by Patke et al. in Cell (2017)55 Patke et al. in Cell (2017)
PMID 28388406
identified CRY1Δ11 in affected members of seven unrelated Turkish families
with familial DSPD. In a validation cohort of 70 subjects from six families
(8 homozygous carriers, 31 heterozygous carriers, 31 non-carriers), the
variant segregated with DSPD with a Fisher's exact P < 0.0001 and an odds
ratio of 1,928 — an effect size rarely seen in common-disease genetics. In
temporal isolation experiments, the proband's free-running circadian period
measured 24.52 hours. In vitro, mouse embryonic fibroblasts expressing the
mutant CRY1 showed a period lengthened by approximately 30 minutes relative
to cells expressing wild-type CRY1. The variant frequency in the gnomAD
database is approximately 0.4% globally, rising to 0.65% in non-Finnish
Europeans and ~3% in Ashkenazi Jewish populations, consistent with
Patke et al.'s estimate that roughly 1 in 75 people of certain ancestries
carry this allele66 Patke et al.'s estimate that roughly 1 in 75 people of certain ancestries
carry this allele
Making it one of the most common single-gene causes of
a sleep disorder ever identified.
A follow-up observational study by Smieszek et al. (2021)77 observational study by Smieszek et al. (2021)
Sci Rep, PMID 34635699
enrolled 67 participants (33 CRY1Δ11 carriers, 34 wild-type controls) from
Turkish families and confirmed that carriers had significantly later wake
times, sleep midpoints, and longer sleep-onset latency. Remarkably, the
circadian delay extended to metabolic outputs: bowel movement timing was
approximately 91 minutes later in carriers (p = 0.002), demonstrating that
the lengthened period affects the entire peripheral clock network, not just
the central sleep-wake system.
A 2020 study by Onat et al. in JCI88 Onat et al. in JCI
PMID 32538895
extended the phenotype. Among 96 individuals from 12 Turkish families with
combined ADHD and DSPD, CRY1Δ11 was present in 13% of affected individuals
versus 0% of controls (OR 281, P = 1.99 × 10⁻²¹). A phenome-wide association
study in 9,438 unrelated European adults found the variant associated with
major depressive disorder, insomnia, anxiety, and nicotine dependence. Of 48
CRY1Δ11 carriers with available psychiatric records, 46 (96%) displayed ADHD
symptoms, and 64% had a history of recurrent depression compared with 10% of
non-carriers.
Practical Actions
The circadian delay caused by CRY1Δ11 is mechanistically fixed — the protein is more active than normal — so the goal of treatment is to externally counteract the lengthened period rather than pharmacologically correct it. Two evidence-supported tools phase-advance the circadian clock: morning bright light and evening melatonin.
Morning bright light exposure (10,000 lux for 30 minutes immediately after
waking) suppresses residual melatonin and signals the suprachiasmatic nucleus
to advance the clock. Evening low-dose melatonin (0.5–3 mg taken 5–7 hours
before desired sleep onset) directly phase-advances the melatonin rhythm.
Combined bright light plus melatonin produces larger phase shifts than either
alone99 Combined bright light plus melatonin produces larger phase shifts than either
alone
Wilhelmsen-Langeland et al. J Biol Rhythms 2013, PMID 24132057.
Because CRY1Δ11 carriers have a genuinely longer intrinsic period, they may
need to sustain these interventions indefinitely rather than using them as
a one-time correction.
Light avoidance in the evening is equally important. Evening light — especially blue-wavelength light from screens — delays the circadian clock by suppressing melatonin release. For carriers whose clock already runs late, evening light exposure amplifies the misalignment. Blue-light filtering glasses (amber lenses) from roughly 2 hours before desired bedtime reduce this phase-delaying input.
Homozygous carriers (GG) show no more severe clinical phenotype than heterozygous carriers in the published family data, consistent with the dominant gain-of-function mechanism reaching its ceiling with a single copy.
Interactions
CRY1Δ11 operates at the core of the CLOCK:BMAL1 → PER/CRY negative feedback loop. Other clock gene variants that affect this same loop can modulate the overall period length in concert with CRY1Δ11. The CRY2 gene encodes a paralogous cryptochrome; variants in CRY2 associated with earlier chronotype could theoretically counteract some period lengthening, though no published compound analysis exists for the CRY1Δ11 and CRY2 combination in humans.
The rs2287161 variant in CRY1 (a common intronic SNP) has been associated with major depressive disorder and depression risk in multiple populations, and may modulate baseline CRY1 expression levels independently of the CRY1Δ11 splice defect. Carriers of CRY1Δ11 who also carry rs2287161 risk alleles may have a higher burden of mood symptoms than CRY1Δ11 alone predicts.
The PERIOD genes PER1, PER2, and PER3 interact directly with CRY1 protein in the feedback loop. Common variants in PER3 (particularly the VNTR polymorphism rs57875989) affect sleep architecture and circadian preferences independently; their interaction with CRY1Δ11 has not been systematically studied but represents a plausible compounding pathway.
GDF2 Arg68Leu — When BMP9 Cannot Mature
Your blood vessels are not static pipes. They are living structures continuously
reshaped by molecular signals — and one of the most important of those signals
in the vascular system is BMP911 BMP9
BMP9 (bone morphogenetic protein 9) is a secreted
ligand of the TGF-beta superfamily. Despite the name, it is primarily active in
vascular biology, not bone — it is the physiological ligand for the endothelial
receptor ALK1 and its co-receptor endoglin.
The GDF2 gene encodes BMP9, and the p.Arg68Leu variant — a single amino acid
substitution in the prodomain of the protein — disrupts BMP9's ability to mature
into its active form, impairs signaling through the ALK1 pathway, and has been
linked to a rare vascular malformation syndrome now classified as hereditary
hemorrhagic telangiectasia type 5 (HHT5).
HHT is an autosomal dominant disorder affecting an estimated 1 in 5,000 people globally. It is characterized by abnormal blood vessel formations — telangiectases (tiny dilated vessels visible on skin and mucous membranes) and arteriovenous malformations (AVMs) in internal organs including the lungs, liver, and brain. Recurrent nosebleeds (epistaxis) are almost universally the first symptom. GDF2 variants account for a small fraction of HHT cases, and the clinical phenotype of HHT5 overlaps with but differs from the more common HHT1 (ENG) and HHT2 (ACVRL1) forms in ways that have clinical management implications.
The Mechanism
BMP9 is secreted as a large precursor — a prodomain22 prodomain
The prodomain is a regulatory
segment that must be cleaved off to release the active mature ligand. For BMP9,
the prodomain and the mature domain remain non-covalently associated after cleavage,
forming what is called a 'procomplex'
attached to a mature signaling domain. Cleavage by furin-family proteases releases
the active mature BMP9 dimer, which then binds to ALK1 and its co-receptors
(endoglin, BMPR2) on endothelial cell surfaces.
The Arg68 residue sits within the prodomain. In vitro expression studies of the
p.Arg68Leu mutant showed that the precursor protein is expressed normally but
fails to produce mature BMP9 dimer efficiently33 fails to produce mature BMP9 dimer efficiently
Processing to mature BMP9 was
far less efficient for p.Arg68Leu than wild-type protein; in ALK1 signaling assays
the variant showed approximately 79% of wild-type activity at 10 pg/ml and 79% at
30 pg/ml equivalent dilutions — a modest but consistent and reproducible reduction.
The biological consequence is reduced BMP9/ALK1/SMAD1/5/8 signaling in endothelial
cells, which normally suppresses pathological angiogenesis and maintains AV
differentiation. Loss of this brake promotes disorganized vascular sprouting,
telangiectasia formation, and AV fistula development.
The Evidence
The p.Arg68Leu variant was identified by Wooderchak-Donahue et al. (2013)44 Wooderchak-Donahue et al. (2013)
Wooderchak-Donahue WL et al. Am J Hum Genet 2013 93:530-7 — sequenced 191 patients
with suspected HHT who were negative for ENG, ACVRL1, and SMAD4; identified three
pathogenic GDF2 missense variants including Arg68Leu
in a proband who met clinical HHT criteria (epistaxis and cutaneous telangiectases).
The variant co-segregated with disease in family members: the proband's father and
sister both carried p.Arg68Leu and reported epistaxis. Functional assays in C2C12
and ATDC5 cells transfected with ALK1 confirmed reduced BMP9 signaling. This
variant was absent from 5,400 control exomes, 1000 Genomes, and dbSNP at time
of publication, consistent with a rare pathogenic allele under negative selection.
The variant has since been documented in subsequent case literature. A 2025 pediatric HHT cohort study at CHOP identified the p.Arg68Leu variant in one patient (patient 6 in that series), who presented with epistaxis and mucocutaneous telangiectases; the patient's father, carrying the same variant, also had epistaxis — confirming the family segregation pattern originally described.
ClinVar (VCV000088651) classifies the G>T transversion (p.Arg68Leu) as "Likely Pathogenic" for HHT5, based on functional evidence of altered GDF2 protein processing and the clinical presentation of the carrier and affected family members. A second variant at the same codon — p.Arg68His (G>A, ClinVar VCV000646500) — is classified as "Uncertain Significance," highlighting that not every amino acid change at position 68 has the same functional impact.
Separate work has reinforced the dose-sensitivity of BMP9 signaling. Chomette
et al. (2023)55 Chomette
et al. (2023)
Chomette L et al. Am J Med Genet A 2023 191:2157-2167 — described
a child with homozygous GDF2 missense at the cleavage site who had both PAH and
HHT features; heterozygous parents were entirely asymptomatic, illustrating how
homozygous loss of BMP9 is far more severe than haploinsufficiency
showed that BMP9 processing mutations can cause pediatric PAH. Heterozygous parents
of the affected child remained asymptomatic, consistent with variable expressivity
in GDF2-HHT5 generally.
Practical Actions
HHT5 is actionable: vascular screening detects AVMs before they become symptomatic emergencies. Pulmonary AVMs can cause paradoxical emboli (stroke) and hemoptysis. Brain AVMs risk hemorrhagic stroke. Hepatic AVMs cause high-output heart failure in advanced disease. Screening protocols adapted from HHT1/HHT2 guidelines are applied to GDF2 carriers.
The phenotype of GDF2-HHT5 may be milder and less penetrant than HHT1/HHT2 — some carriers have only epistaxis with no solid organ AVMs. Nonetheless, clinical evaluation and baseline screening are warranted in all confirmed or suspected carriers, as visceral AVMs are present in a meaningful proportion of published cases and can be clinically silent until they cause an acute event.
Genetic counseling is essential: HHT5 is autosomal dominant, meaning each first-degree relative of a carrier has a 50% chance of inheriting the variant. The variable expressivity means a parent with only mild epistaxis can have a child with pulmonary AVMs.
Interactions
GDF2 (BMP9) signals through the same endothelial receptor complex as ENG (HHT1) and ACVRL1/ALK1 (HHT2). Variants in ENG and ACVRL1 that reduce receptor availability or signaling would be expected to compound with GDF2 loss-of-function, though digenic HHT from GDF2 plus ENG/ACVRL1 has not been formally documented in published case series. Similarly, SMAD4 loss-of-function (HHT-juvenile polyposis overlap syndrome) disrupts the downstream effector of the same pathway.
BMPR2 is the type II receptor for BMP9, and pathogenic BMPR2 variants are the most common cause of hereditary pulmonary arterial hypertension (hPAH). GDF2 variants can cause both HHT5 and PAH through reduced BMP9/ALK1 signaling, suggesting that concurrent BMPR2 variants might amplify PAH risk in GDF2 carriers, though this interaction has not been systematically studied.
FOXO3's East Asian Longevity Signal — An Intronic eQTL With Brain Expression Effects
FOXO3 is the most replicated human longevity gene, with protective variants confirmed across European, Asian, and African populations. While the well-characterized rs2802292 G-allele (the HSF1-binding enhancer variant) has been studied primarily in Western and Japanese-American cohorts, rs2253310 captures a distinct signal that may be especially informative for East Asian populations: here, the longevity-protective C allele is the minority allele (~27% in Japanese), while the less-favorable G allele predominates (~73%). This reversed frequency pattern makes rs2253310 a particularly useful genetic marker for East Asian longevity research.
A 2022 longitudinal study11 A 2022 longitudinal study
Ji JS, Liu L, Yan LL, Zeng Y. Comparing Effects of FOXO3 and Residing
in Urban Areas on Longevity: A Gene-Environment Interaction Study. J Gerontol A Biol Sci Med Sci.
2022 followed 3,085 Chinese older adults and found CC
homozygotes had a 19% lower mortality hazard compared to GG homozygotes (HR 0.808, 95% CI
0.667–0.978), a magnitude comparable to the survival benefit of urban versus rural residence in
the same cohort.
The Mechanism
rs2253310 sits in intron 2 of FOXO3 at chromosomal position 108,567,389 (GRCh38, chromosome 6). The gene is on the plus strand, so no strand-complementing is needed — genome files report the same alleles used in publications. The variant does not alter FOXO3 protein sequence; instead, it acts as an [expression quantitative trait locus (eQTL) | a variant that affects how much of a gene product is made, without changing the protein structure], influencing FOXO3 transcription.
A meta-analysis of four longevity cohorts22 A meta-analysis of four longevity cohorts
Bae H, Gurinovich A, Malovini A, et al. Effects of
FOXO3 Polymorphisms on Survival to Extreme Longevity in Four Centenarian Studies. J Gerontol A
Biol Sci Med Sci. 2018 found that among all 17 tested
FOXO3 variants, rs2253310 and rs6911407 showed the most significant effects on FOXO3 expression
in brain tissue — a key finding given FOXO3's role in neuronal stress resistance, autophagy, and
protection against age-related neurodegeneration.
The rs2253310 G allele's adverse biological direction was confirmed by a 2025 Russian study (n=1,365) showing GG homozygotes face nearly double the risk of chronic obstructive pulmonary disease (OR 1.99, p = 5.93×10⁻⁷) , consistent with lower FOXO3-mediated antioxidant defense in lung tissue of G-allele carriers.
The Evidence
Three independent Chinese longitudinal cohort studies using the Chinese Longitudinal Healthy Longevity Survey demonstrate consistent protective effects of the rs2253310 C allele:
Liu et al. 202133 Liu et al. 2021
Liu L, Zhu A, Shu C, Zeng Y, Ji JS. Gene-Environment Interaction of FOXO and
Residential Greenness on Mortality Among Older Adults. Rejuvenation Res. 2021
studied 3,179 adults aged 65 and older, finding CC homozygotes had an HR of 0.803 (95% CI
0.666–0.968) for all-cause mortality. Notably, the protective effect of residential greenness
was amplified in C-allele carriers, suggesting gene-environment synergy in FOXO3 activation.
Ji et al. 2022 cognitive study44 Ji et al. 2022 cognitive study
Ji JS, Liu L, Zeng Y, Yan LL. Effect of FOXO3 and Air Pollution
on Cognitive Function. J Gerontol A Biol Sci Med Sci. 2022
tracked cognitive function over 14 years and found C-allele homozygotes had higher baseline
MMSE scores and lower odds of cognitive impairment over time. The cognitive protection was
strongest in women, older adults, and those in lower air-pollution environments, suggesting that
the C allele's FOXO3 expression boost matters most when cellular stress is high.
The blood pressure findings from Morris et al. 2016 add a cardiovascular dimension: among 843 Japanese Americans, women carrying two C alleles had 6 mmHg lower systolic and 3 mmHg lower diastolic blood pressure than GG homozygotes, with essential hypertension prevalence of 3.3% vs 9.5% (P = 0.03–0.04). This gender specificity echoes the pattern seen with other FOXO3 longevity variants, where protective effects are often stronger in women.
Practical Implications
The practical takeaways from rs2253310 mirror those of other FOXO3 longevity variants — all of them point toward the same lifestyle levers for activating FOXO3 expression: fasting, exercise, and stress management. The unique contribution of rs2253310 is its cognitive aging signal: the C allele's strongest effects appear in brain tissue eQTLs, suggesting that brain health across the lifespan may be particularly relevant to this variant.
For GG homozygotes — especially prevalent in East Asian populations — the cognitive aging data suggest that interventions targeting neurological resilience deserve particular emphasis: regular aerobic exercise (which strongly activates FOXO3 in neural tissue), sleep optimization (FOXO3 expression cycles with circadian rhythms), and minimizing chronic oxidative stressors such as air pollution exposure and smoking.
Interactions
rs2253310 lies in the same FOXO3 intron 2 region as rs2802292 and is in high linkage disequilibrium with the broader FOXO3 longevity haplotype. In East Asian populations, the LD pattern differs from European cohorts — rs2253310's C allele is uncommon where rs2802292's G allele may be more common, meaning the two variants may be partially independent signals in Asian ancestries, which could explain rs2253310's distinct detection in Chinese cohort studies rather than primarily in European studies.
The gene-environment interactions documented for rs2253310 — with residential greenness, urban living, and air pollution — are not replicated for most other FOXO3 variants, suggesting this particular intronic position may regulate FOXO3 expression in response to environmental oxidative stress in a way that differs mechanistically from the HSF1-dependent rs2802292 enhancer.
ABCA1 rs2575876 — A Recessive Dimmer Switch on Cholesterol Efflux
ABCA111 ATP-Binding Cassette Transporter A1 — a large cell-membrane protein that
pumps cholesterol and phospholipids out of cells onto lipid-poor apolipoprotein A-I,
the initiating step for HDL particle assembly is the rate-limiting controller of
reverse cholesterol transport22 reverse cholesterol transport
the pathway by which excess cholesterol is ferried
from peripheral tissues back to the liver for recycling or excretion. When ABCA1
works well, cells shed cholesterol efficiently, nascent HDL forms rapidly, and
cardiovascular risk declines. rs2575876 is a common intronic variant in the ABCA1
gene on chromosome 9 that modulates this process in a dose-dependent way — but its
most clinically significant effect emerges primarily when both copies carry the A allele.
The Mechanism
rs2575876 sits within an intron of ABCA1 at GRCh38 position 104,903,458. The ABCA1 gene transcribes from the minus strand, so the alleles reported by genome sequencing files (plus strand) are G (reference) and A (alternate). This is a non-coding variant with no direct amino acid change. Its location within the intron places it in a region rich in regulatory elements that have been shown to fine-tune ABCA1 transcription in hepatocytes and other tissues.
The neighboring variant rs257587533 rs2575875
a closely related intronic SNP in LD with
rs2575876 and rs4149268 in this ABCA1 locus, studied by Howard et al.44 Howard et al.
Howard et al. Allele-specific enhancers mediate associations between LCAT and ABCA1
polymorphisms and HDL metabolism. PLOS One, 2019,
sits in a transcriptional enhancer that creates an allele-specific STAT355 Signal
Transducer and Activator of Transcription 3 — a transcription factor activated by
cytokines that loops from intronic enhancers to the ABCA1 promoter, driving
hepatic expression binding site. These regulatory variants in tight LD form a
haplotype block in which the cumulative effect of homozygosity matters more than any
single base change. AA homozygotes at rs2575876 sit at the weaker end of this
regulatory spectrum.
The Evidence
The strongest direct evidence comes from a cohort study by Lin et al.66 Lin et al.
Lin et al. Association between high-density lipoprotein and functional outcome of
ischemic stroke patients in a Taiwanese population. Lipids Health Dis, 2024
of 1,310 first-ever acute ischemic stroke patients. Among two ABCA1 SNPs tested,
rs2575876 and rs1883025 were the only two significantly associated with HDL-C levels
in the total population and in sex-stratified subgroups. Under a recessive model —
AA versus GG+GA combined — rs2575876 AA homozygotes showed measurably different
HDL-C levels and a significantly elevated risk of poor functional outcomes at 1 and
3 months post-stroke, particularly when HDL-C was simultaneously abnormal (too low
or too high). This recessive pattern implies that a single A allele is largely
compensated by the G allele, but two A alleles expose the full effect of reduced
ABCA1 regulatory activity.
Genome-wide association data add breadth: multiple large GWAS including the Global Lipids Genetics Consortium (>1.65 million individuals) have identified the ABCA1 locus on chromosome 9 as genome-wide significant for triglycerides and LDL-C, with rs2575876 named as one of the contributing tag SNPs for both traits. Effect sizes are modest — typical of common intronic regulatory variants — but consistent across diverse ancestries.
A broader review by Frikke-Schmidt77 Frikke-Schmidt
Frikke-Schmidt R. Genetic variation in the ABCA1 gene, HDL cholesterol, and risk of
ischemic heart disease. Atherosclerosis, 2010
underscores an important nuance: while ABCA1 variants lower HDL-C, genetically low
HDL per se does not straightforwardly predict ischemic heart disease risk. The
ABCA1-IHD association appears to be partially independent of HDL levels, suggesting
ABCA1 has cardioprotective roles beyond cholesterol efflux capacity alone.
Practical Actions
Because the AA genotype is present in only about 5–6% of the population globally, most people carry at least one G allele and will not face the full recessive effect. For AA homozygotes, the priority is monitoring HDL-C and addressing any confirmed deficit through strategies known to improve ABCA1-driven cholesterol efflux: reducing trans fat intake (trans fats directly suppress ABCA1 expression), substituting monounsaturated and omega-3 fats, and sustaining aerobic exercise. If HDL remains below target despite lifestyle efforts, niacin-based options can raise HDL substantially.
The stroke interaction finding signals that when HDL-C goes outside the normal range in AA carriers, recovery from vascular events is compromised. Maintaining HDL-C within normal bounds is therefore especially important for this genotype.
Interactions
The two ABCA1 SNPs most co-studied with rs2575876 are rs1883025 (a nearby variant on the same haplotype block) and rs4149268 (an upstream intronic variant in the same regulatory region). All three tag overlapping ABCA1 regulatory signals; carrying risk alleles across multiple ABCA1 loci may have an additive dampening effect on ABCA1 expression and HDL biogenesis. The R219K missense variant rs2230806 acts through a different mechanism — reducing efflux protein function rather than expression — and may compound the effect of intronic risk alleles. Any ABCA1 variant in combination with APOE ε4 represents a convergent challenge to brain cholesterol homeostasis, given ABCA1's role in neuronal cholesterol efflux.