rs10984447

DBC1 DBC1 multiple sclerosis susceptibility variant

Moderate Risk Factor

DBC1 9q33 — A Common Variant That Tips the Autoimmune Balance Toward MS

The DBC1 gene11 DBC1 gene
DBC1 (Deleted in Bladder Cancer 1), also known as DBCCR1 and BRINP1 (BMP/Retinoic Acid Inducible Neural Specific 1), sits at chromosome 9q33.1. The gene is expressed at highest levels in the brain and plays roles in DNA damage response, epigenetic regulation via SIRT1 inhibition, and cellular stress signaling
sits at one of the genomic regions repeatedly flagged by multiple sclerosis genome-wide association studies. Its most common variant, rs10984447, illustrates an unusual feature of MS genetics: the reference (major) allele is the susceptibility allele, and rarer protective variant carriers are the genetic minority.

Multiple sclerosis is a chronic neuroinflammatory disease22 chronic neuroinflammatory disease
MS affects approximately 2.8 million people worldwide; autoimmune CD4+ T cells drive demyelination of the central nervous system, leading to relapsing or progressive neurological disability
caused by misdirected immune attacks on central nervous system myelin. Genetic predisposition accounts for roughly 20–30% of MS risk, with the HLA-DRB1*15:01 allele as the dominant single factor (OR ≈ 3.1). The remaining non-HLA genetic risk is distributed across over 100 loci, each contributing modest effect sizes — rs10984447 at 9q33 is among these smaller-effect contributors.

The Mechanism

DBC1 protein is a nuclear scaffold that physically binds and inhibits SIRT133 SIRT1
SIRT1 is an NAD+-dependent deacetylase that regulates inflammation, DNA repair, and immune cell differentiation; DBC1-SIRT1 interaction modulates NF-κB activity and T-cell activation thresholds
, keeping the deacetylase in check under basal conditions. Cellular stress signals (genotoxic damage, oxidative stress) disrupt the DBC1–SIRT1 complex and activate SIRT1-dependent gene programs involved in DNA damage response and apoptosis. Because SIRT1 is a central regulator of T-cell differentiation — suppressing Th17 fate while promoting Treg stability — subtle changes in DBC1 expression or function in immune cells could shift the balance toward pro-inflammatory T-cell programs associated with autoimmune pathology.

rs10984447 is an intronic variant located within a CpG island that is a known hypermethylation target in cancer. Whether this specific variant alters splicing, expression, or epigenetic marking at the DBC1 locus in immune cells has not been directly demonstrated — the MS association is a genome-wide tagging signal, and the precise functional mechanism remains to be established. The variant may mark a regulatory haplotype rather than having a direct intrinsic effect.

The Evidence

The primary MS association comes from the International Multiple Sclerosis Genetics Consortium 2011 GWAS44 International Multiple Sclerosis Genetics Consortium 2011 GWAS
Sawcer et al. analysed 9,772 MS cases and 17,376 controls of European ancestry plus an independent replication cohort, identifying multiple novel loci outside the HLA region with effects in the OR 1.1–1.3 range
. rs10984447 at the DBC1 locus emerged with OR ≈ 1.17 (p = 8 × 10⁻⁶ at genome-wide suggestive threshold, but below the conventional 5 × 10⁻⁸ genome-wide significance threshold).

Replication in different ancestry groups adds some robustness. Al Jumah et al. (2012)55 Al Jumah et al. (2012)
342 subjects: 99 sporadic MS, 22 familial MS, 89 related controls, 132 independent controls in a Saudi Arabian cohort — genetically homogenous population reduces confounding
found a substantially larger effect size in familial MS cases vs independent controls (OR 13.63, 95% CI 1.54–120.83). The wide confidence interval reflects the small familial MS group (n=22), but replication in a non-European population with independent methodology supports the variant's genuine relevance.

The broader 9q33 MS signal sits within the context of Beecham et al. (2013)66 Beecham et al. (2013)
ImmunoChip analysis of immune-related loci in over 80,000 individuals of European ancestry; confirmed 110 non-HLA MS susceptibility variants, establishing the polygenic landscape within which this locus operates
. The modest individual effect (OR 1.17) is consistent with the general architecture of non-HLA MS genetics — many loci, each contributing small increments of risk that compound across the genome.

The evidence level is rated moderate: the 2011 GWAS signals a locus-level association but rs10984447 did not reach genome-wide significance independently, and the causal variant at 9q33 has not been fine-mapped. The Saudi replication supports the signal without resolving the mechanistic question.

Practical Actions

For those carrying one or two A alleles: the individual effect of this variant is small (OR 1.17 per allele from the European cohort). Carrying AA does not mean MS is likely — lifetime MS incidence is approximately 0.3% in the general population, and a 17% per-allele increase brings this to a still-small absolute risk increment. Risk is meaningfully elevated only when this variant combines with high-impact MS risk alleles (notably HLA-DRB1*15:01, IL7R rs6897932, and PTPN22 rs2476601).

The actionable focus is monitoring for early MS symptoms — particularly episodes of optic neuritis, limb weakness, or sensory disturbances — and ensuring vitamin D sufficiency, which has the strongest environmental evidence for reducing MS incidence in genetically susceptible individuals.

Interactions

The DBC1 locus contributes to MS risk additively with other established non-HLA MS susceptibility variants. In particular, IL7R rs6897932 (CC genotype) and PTPN22 rs2476601 (AA genotype) operate through distinct immune cell mechanisms (IL-7 signaling and T-cell activation thresholds, respectively) and compound with 9q33 variants to further elevate polygenic risk. Individuals carrying risk alleles at multiple MS loci have substantially greater susceptibility than the individual SNP odds ratios would predict in isolation.

SCARB1 — Your Gateway for HDL-Carried Vitamin E and Carotenoids

Every fat-soluble micronutrient you absorb --- vitamin E, beta-carotene, lutein, zeaxanthin --- faces a double barrier: the intestinal wall and then the liver. At both sites, a single receptor protein called SR-BI11 SR-BI
Scavenger receptor class B type I, encoded by SCARB1 on chromosome 12q24; an integral membrane protein with a large hydrophobic tunnel that channels lipophilic molecules into cells
(scavenger receptor class B type I, encoded by SCARB1) acts as the primary gateway. SR-BI is expressed in enterocytes lining the small intestine and in hepatocytes, where it performs the selective uptake of HDL22 HDL
High-density lipoprotein particles that carry fat-soluble vitamins and cholesterol from peripheral tissues back to the liver
-associated lipids --- including cholesterol, tocopherols, and carotenoids. The rs11057830 variant lies within an intron of SCARB1. While it does not change the protein's amino acid sequence, it is a tag SNP for a haplotype that influences the receptor's functional expression or activity. People carrying more copies of the G allele show lower circulating alpha-tocopherol than those with the A allele, a relationship confirmed in multiple cohorts totaling over 7,000 individuals.

The Mechanism

SR-BI's lipid transport mechanism relies on a large hydrophobic tunnel that threads through the membrane. Lipid molecules associate with HDL particles docked at the receptor's extracellular domain, then pass through this tunnel into the cell. The selectivity of this tunnel33 selectivity of this tunnel
Li et al. 2023 showed that blocking the cholesterol tunnel with a C384Y mutation abolishes the preferential uptake of lutein and zeaxanthin, while beta-carotene uptake is less affected
explains why SCARB1 shows different uptake efficiencies for different carotenoids: the more hydrophilic xanthophylls (lutein, zeaxanthin) navigate the tunnel most efficiently, while the purely hydrophobic beta-carotene relies on it less. The key functional finding for vitamin E is dramatic: Reboul et al. 200644 Reboul et al. 2006
Reboul E et al. Scavenger receptor class B type I (SR-BI) is involved in vitamin E transport across the enterocyte. J Biol Chem, 2006
showed that blocking SR-BI with a specific inhibitor (BLT1) reduced alpha-tocopherol transport across cultured enterocytes by up to 80%. In mice overexpressing intestinal SR-BI, gamma-tocopherol bioavailability was 2.7-fold higher than in wild-type animals. This positions SR-BI as the dominant --- not merely contributing --- protein for vitamin E intestinal absorption. RS-BI also handles the final delivery step. In brain capillary endothelial cells, adenovirus-mediated SR-BI overexpression55 adenovirus-mediated SR-BI overexpression
Goti D et al. Scavenger receptor class B, type I is expressed in porcine brain capillary endothelial cells and contributes to selective uptake of HDL-associated vitamin E. J Neurochem, 2001
produced a 4-fold increase in HDL-associated alpha-tocopherol uptake, demonstrating that SR-BI shapes vitamin E delivery not just at the gut but at every tissue that relies on HDL-mediated supply. How rs11057830 specifically alters SR-BI function is not fully characterized --- it may tag a haplotype affecting promoter activity, splicing, or a nearby regulatory element. What is established is the functional consequence: fewer copies of the A allele correlates with lower circulating tocopherol.

The Evidence

The primary evidence comes from a genome-wide association study66 genome-wide association study
Major JM et al. Genome-wide association study identifies common variants associated with circulating vitamin E levels. Hum Mol Genet, 2011
conducted across 5,006 men of European descent in two cohorts, including the Alpha-Tocopherol, Beta-Carotene Cancer Prevention (ATBC) Study. After controlling for age, BMI, cholesterol, and supplementation status, rs11057830 was associated with circulating alpha-tocopherol at genome-wide significance (P = 2.0 × 10⁻⁸). The A allele conferred higher levels (beta = +0.04 on a log scale per allele), with mean alpha-tocopherol increasing from 11.8 mg/L (GG) to 12.2 mg/L (GA) to 12.7 mg/L (AA). The finding was replicated in 992 men from the PLCO trial and 2,775 women from the Nurses' Health Study, with the combined meta-analysis reaching P = 8.2 × 10⁻⁹. Together with CYP4F2 rs2108622 and the APOA5 locus rs964184, rs11057830 explained 1.7% of the variance in circulating alpha-tocopherol. For carotenoids, related SCARB1 variants in strong linkage disequilibrium with rs11057830 (particularly rs11057841, r² ≈ 0.7) have been associated with serum lutein levels. A study of 302 healthy subjects77 study of 302 healthy subjects
McKay GJ et al. Investigation of genetic variation in scavenger receptor class B, member 1 (SCARB1) and association with serum carotenoids. Ophthalmology, 2013
found that rs11057841 predicted a 24% increase in serum lutein per T allele (P < 0.01 after permutation correction), with independent replication in the TwinsUK and CAREDS cohorts. Given the strong LD between rs11057841 and rs11057830, this carotenoid association likely applies to the same functional haplotype.

Practical Implications

The 7--8% difference in circulating alpha-tocopherol between GG and AA genotypes is modest in absolute terms. However, for GG carriers --- roughly 73% of people --- it represents the lower end of normal SR-BI-mediated uptake. Because SR-BI accounts for up to 80% of intestinal vitamin E transport, optimizing dietary delivery is directly relevant. Two factors amplify practical impact. First, vitamin E absorption is highly fat-dependent --- consuming it with a meal containing dietary fat improves micellarization and SR-BI access. Second, the natural RRR-alpha-tocopherol form is more efficiently retained by the liver (through the alpha-tocopherol transfer protein, alpha-TTP) than synthetic dl-alpha-tocopherol; for people with already-reduced SR-BI uptake, maximizing what reaches the liver is important. The carotenoid connection is worth noting. Beta-carotene, lutein, and zeaxanthin all use SR-BI for intestinal absorption. GG carriers absorbing less through this pathway may benefit from optimizing carotenoid bioavailability --- eating cooked rather than raw vegetables, pairing with fat, and favoring foods with high carotenoid density.

Interactions

This variant acts at the absorption step of the vitamin E pathway. A functionally distinct step is handled by rs6994076 in TTPA, which governs how much alpha-tocopherol the liver redistributes to tissues via VLDL particles. Individuals who carry both rs11057830 GG (reduced SR-BI absorption) and rs6994076 TT (reduced alpha-TTP expression) face compounded limitations: less vitamin E absorbed through the intestine AND less efficiently retained and distributed by the liver. This combination represents the most actionable double-hit in the vitamin E pathway, warranting both dietary and supplementation attention. The compound effect of GG at rs11057830 + TT at rs6994076 should be modeled as a compound action recommending higher-dose natural tocopherol supplementation with fat-rich meals. For carotenoids, rs12934922 in BCO1 (beta-carotene 15,15'-oxygenase) interacts with SCARB1 at a downstream step: rs11057830 determines how much beta-carotene enters the body; BCO1 rs12934922 determines how efficiently it is converted to retinol. Poor converters (BCO1 TT) who also carry SCARB1 GG absorb less beta-carotene AND convert less of it to vitamin A.

IL-6 Receptor Signaling and Atrial Fibrillation Risk

The interleukin-6 receptor (IL-6R) encoded by IL6R is the gateway through which interleukin-6, one of the body's primary inflammatory messengers, communicates its instructions to target cells. When IL-6 binds its receptor, it triggers STAT3 phosphorylation11 STAT3 phosphorylation
Signal Transducer and Activator of Transcription 3 — a transcription factor that turns on inflammatory gene programs
, driving the liver to produce acute-phase proteins including C-reactive protein (CRP), fibrinogen, and serum amyloid A. Chronically elevated IL-6 signaling through IL-6R accelerates atherosclerosis, promotes atrial fibrosis, and contributes to the pro-inflammatory milieu that sustains cardiac arrhythmias. The rs11265611 variant sits within an intron of IL6R and serves as a tag for a functional haplotype block that modulates IL-6 receptor signaling efficiency.

The Mechanism

rs11265611 is an intronic variant that does not directly change the IL-6 receptor protein sequence. Instead, it marks — via linkage disequilibrium22 linkage disequilibrium
LD: the tendency for certain allele combinations to be inherited together more often than chance predicts
— a haplotype block spanning the IL6R locus that influences receptor expression or signaling capacity. This locus is in strong LD (r²=0.69) with rs4845625, an IL6R intronic variant whose T allele is associated with elevated CRP, LDL-C, and apolipoprotein B, and with rs2228145 (Asp358Ala), a missense variant that alters receptor shedding efficiency and thereby controls the balance between membrane-bound (pro-inflammatory) and soluble (anti-inflammatory) forms of IL-6R.

The G allele at rs11265611 corresponds to the higher-signaling haplotype: carriers tend to have higher baseline IL-6 pathway activity, producing more CRP and fibrinogen and maintaining a more pro-inflammatory systemic state. This persistent low-grade inflammation is mechanistically linked to atrial fibrosis — the structural remodeling of atrial tissue that creates the electrical substrate for atrial fibrillation to initiate and sustain itself.

The Evidence

The primary evidence for rs11265611 comes from the CHARGE Targeted Sequencing Study33 CHARGE Targeted Sequencing Study
Lin et al. Targeted sequencing in candidate genes for atrial fibrillation. Heart Rhythm, 2014
, which sequenced candidate genes in 948 AF cases and 3,330 controls from four major cohorts (ARIC, CHS, Framingham, MGH). The study found rs11265611 associated with a 30% reduction in AF risk per protective allele (OR 0.70, 95% CI 0.58–0.85, p=1.70×10⁻⁶ after Bonferroni correction) — a finding that prior genome-wide studies had not captured because this specific variant was not genotyped or imputed.

The nearby rs4845625 provides additional mechanistic evidence. Its T allele — which the rs11265611 G allele tags — was independently associated with AF recurrence after catheter ablation44 AF recurrence after catheter ablation
Wu et al. A variant of IL6R is associated with the recurrence of atrial fibrillation. Heart Rhythm, 2014
in 278 Chinese Han patients (early recurrence OR 1.84, late recurrence OR 1.92). These findings from European and East Asian populations support a cross-ancestry role for IL6R variation in AF susceptibility.

At the Mendelian randomization level, Cupido et al. (2022)55 Cupido et al. (2022)
Dissecting the IL-6 pathway in cardiometabolic disease: a Mendelian randomization study. Br J Clin Pharmacol, 2022
demonstrated that genetically lower CRP mediated by IL6R variants causally reduces AF risk (OR 0.90, 95% CI 0.86–0.95 per 1 mg/L lower CRP). This establishes that the IL6R-inflammation pathway is not merely associated with AF but is likely in the causal chain — a finding with direct implications for anti-inflammatory cardiovascular prevention strategies.

Practical Actions

For G allele carriers, the primary implication is elevated baseline IL-6 pathway activity contributing to higher atrial fibrillation risk. This genotype reinforces the case for monitoring systemic inflammation via high-sensitivity CRP (hs-CRP) and for targeted anti-inflammatory strategies where evidence supports them. Colchicine has demonstrated reduction in post-ablation AF recurrence in meta-analyses Kommu et al., Cureus 202366 Kommu et al., Cureus 2023 (RR 0.57–0.58 at 3 and 12 months), making it a genotype-informed option to discuss with a cardiologist before ablation procedures.

Interactions

rs11265611 is in LD with rs4845625 (r²=0.69) and rs2228145 (Asp358Ala), and functions as part of a broader IL6R haplotype block. The Asp358Ala variant (rs2228145) directly alters IL-6 receptor ectodomain shedding and has independent mechanistic evidence for changing soluble IL-6R levels, CRP, and cardiovascular risk. Users carrying risk haplotypes across multiple IL6R markers may experience compounded effects on IL-6 signaling. Downstream, IL-6R risk haplotypes interact with other inflammatory loci — particularly CRP gene variants (rs1205, rs1130864) — because CRP production is itself a readout of IL-6 receptor signaling.

ANKRD30A rs1148259 — A Metabolomics Signal in the Sphingolipid Layer

When geneticists search for the DNA variants that shape our metabolic chemistry, they sometimes find signals in unexpected places. rs1148259 sits in the 3′ untranslated region11 3′ untranslated region
The 3′ UTR is the non-coding tail of a messenger RNA that follows the stop codon. It does not change the protein sequence but often carries regulatory sequences that control mRNA stability, transport, and translation efficiency.
of the ANKRD30A gene on chromosome 10 — technically a synonymous position that leaves the encoded protein unchanged, yet it emerged as one of the strongest metabolite associations in the first large-scale metabolomics genome-wide association study, with circulating sphingolipid levels differing by genotype at p = 3.04 × 10⁻⁹.

The Gene

ANKRD30A (ankyrin repeat domain 30A, also catalogued as the breast cancer antigen NY-BR-1) encodes a nuclear protein bearing multiple ankyrin repeat domains22 ankyrin repeat domains
Ankyrin repeats are 33-amino-acid motifs that mediate protein-protein interactions and are found in hundreds of signalling and regulatory proteins.
. The gene shows highest expression in mammary glandular epithelium and testis, with secondary expression in adipose tissue. It acts as a DNA-binding transcription factor33 DNA-binding transcription factor
A protein that attaches to specific DNA sequences and regulates whether nearby genes are switched on or off, and at what level.
and has been identified as a host factor for HIV-1 replication and as an interferon-stimulated gene, pointing to roles in both cellular regulation and antiviral defence. No direct enzymatic role in sphingolipid synthesis or degradation has been described for ANKRD30A, making the metabolomics signal difficult to interpret mechanistically.

The Mechanism

rs1148259 is annotated as a synonymous variant — all possible nucleotide substitutions at this position encode the same amino acid (alanine at codon 1270). In at least one transcript isoform, the variant falls in the 3′ UTR rather than the coding sequence. Variants in 3′ UTRs can alter mRNA stability, affect binding sites for microRNAs44 microRNAs
Small non-coding RNA molecules (~22 nucleotides) that bind to 3′ UTR sequences and suppress gene expression post-transcriptionally, often by destabilising the mRNA or blocking its translation.
or RNA-binding proteins, and thereby change the amount of protein produced without altering its sequence. Whether this is the mechanism behind the sphingolipid association remains uncharacterised.

Sphingomyelins55 Sphingomyelins
A class of sphingolipids found abundantly in cell membranes, especially myelin sheaths around nerve fibres and in lipid rafts. Sphingomyelin consists of a ceramide backbone linked to a phosphocholine head group.
are structural membrane lipids and signalling molecules involved in lipid raft formation, apoptosis signalling, and inflammation. Circulating sphingomyelin concentrations have been associated with cardiovascular risk, insulin resistance, and neurological outcomes in epidemiological studies.

The Evidence

The association was identified by Gieger et al. (2008)66 Gieger et al. (2008)
Gieger C et al. Genetics meets metabolomics: a genome-wide association study of metabolite profiles in human serum. PLoS Genet, 2008
in the first ever genome-wide metabolomics study. The cohort consisted of 284 European men enrolled in the KORA (Cooperative Health Research in the Region of Augsburg) study in Bavaria, Germany. Using a metabolomics platform measuring 363 serum metabolites and their ratios, rs1148259 reached p = 3.04 × 10⁻⁹ for sphingolipid amount — just below the conventional genome-wide threshold — making it one of roughly a dozen suggestive loci in the study.

This discovery-cohort signal has not, to date, been independently replicated in published literature. The primary metabolomics GWAS literature from this period converged on well-characterised loci — FADS1 for polyunsaturated fatty acids, LIPC for glycerophospholipids, SCAD and MCAD for acylcarnitines — as its core findings, and rs1148259 was not among the loci carried forward into replication efforts that followed. Its evidence status therefore remains emerging: a single discovery association in a small, male-only European cohort with no established biological mechanism.

Practical Implications

Given the emerging and unreplicated nature of this finding, rs1148259 does not currently support strong, genotype-specific dietary or supplementation recommendations. The most clinically actionable step for individuals carrying the C allele is awareness that their sphingolipid metabolism may differ from the average, warranting attention to the dietary factors and biomarkers known to influence sphingolipid levels more broadly.

Dietary choline from eggs, liver, and soybeans is the direct precursor to the phosphocholine head group of sphingomyelin and supports membrane phospholipid synthesis. Adequate serine intake matters because serine is the carbon backbone donor for de novo sphingolipid synthesis. Monitoring standard lipid panels — which include some sphingomyelin-rich lipoprotein fractions indirectly — provides a baseline picture of sphingolipid-related cardiovascular risk.

Interactions

The strongest metabolomics locus in the Gieger 2008 study was rs174548 in FADS1, which controls desaturation of polyunsaturated fatty acids and influences glycerophospholipid and sphingomyelin composition via shared membrane substrates. The PLEK locus (rs9309413) showed associations with multiple sphingomyelin species in the same study. Because sphingolipid metabolism is interconnected with the broader phospholipid network, FADS1 and PLEK variants together with rs1148259 may jointly shape an individual's sphingolipid profile, though no compound analysis has been published for this specific three-way combination.

A Third Independent Signal in the Vitamin D Carrier Gene

The GC gene11 GC gene
Also called the group-specific component gene or DBP gene. It encodes vitamin D binding protein (VDBP), a 58-kDa glycoprotein made mainly by the liver that carries 85-90% of circulating 25(OH)D and 85% of 1,25(OH)₂D in the bloodstream
has long been known as the dominant genetic determinant of circulating vitamin D levels in the blood. Most research has focused on two functional missense variants — rs7041 (Asp432Glu) and rs4588 (Thr436Lys) — and one intronic GWAS tag (rs2282679). rs1155563 is a fourth GC region variant that behaves as a third independent signal22 third independent signal
Independent signal means it retains statistical significance even after conditioning on the other GC variants in the same statistical model, indicating it captures genetic variation not already tagged by rs2282679, rs7041, or rs4588
at the locus, confirmed in the Ahn et al. 2010 GWAS33 Ahn et al. 2010 GWAS
Ahn J et al. Genome-wide association study of circulating vitamin D levels. Hum Mol Genet, 2010
of 4,501 Europeans (P = 3.8 × 10⁻²⁵).

The Mechanism

rs1155563 sits within an intron of the GC gene on chromosome 4 (position 71,777,771, GRCh38) and lies slightly further from the two coding missense variants (rs7041 at 71,752,617 and rs4588 at 71,752,606) than the previously described tag SNP rs2282679 (at 71,742,666). Like all intronic SNPs, it does not change the amino acid sequence of VDBP. Instead, it likely acts as a tag SNP44 tag SNP
A genetic variant that, through linkage disequilibrium — the tendency for nearby alleles to be inherited together — reliably marks a haplotype block carrying one or more functional variants that directly alter protein expression or function
for an independently segregating haplotype within the GC locus that influences VDBP expression level, mRNA stability, or splicing efficiency.

The C allele (minor allele, ~28% frequency in Europeans) is the effect allele associated with lower circulating 25(OH)D. Because GC is transcribed on the minus strand, the coding-strand notation used in many publications differs from the plus-strand alleles reported here (plus-strand T is reference; C is the effect allele). The effect is additive: each additional C allele reduces serum 25(OH)D, and CC homozygotes show the greatest reduction.

The Evidence

The Ahn et al. 2010 GWAS55 Ahn et al. 2010 GWAS
Ahn J et al. Genome-wide association study of circulating vitamin D levels. Hum Mol Genet, 2010
of 4,501 European-ancestry participants identified rs1155563 as a genome-wide significant GC locus variant (P = 3.8 × 10⁻²⁵), with the C allele associated with lower circulating 25(OH)D, remaining significant after conditioning on the two established coding variants. The parallel SUNLIGHT consortium GWAS66 SUNLIGHT consortium GWAS
Wang TJ et al. Common genetic determinants of vitamin D insufficiency: a genome-wide association study. Lancet, 2010
of 33,996 Europeans confirmed the GC locus dominance (P = 1.9 × 10⁻¹⁰⁹ for the lead variant rs2282679) and showed that individual variants at this locus can independently tag different VDBP haplotypes.

Cross-ethnic replication came from a study of 3,210 Chinese Hans77 study of 3,210 Chinese Hans
Lu L et al. Associations between common variants in GC and DHCR7/NADSYN1 and vitamin D concentration in Chinese Hans. Hum Genet, 2012
, which confirmed GC variants including rs1155563 as significant determinants of 25(OH)D (beta −0.036 to −0.076 per risk allele across GC SNPs).

A pediatric GWAS in 761 Finnish infants88 pediatric GWAS in 761 Finnish infants
Kämpe A et al. Genetic variation in GC and CYP2R1 affects 25-hydroxyvitamin D concentration and skeletal parameters: A genome-wide association study in 24-month-old Finnish children. PLoS Genet, 2019
found rs1155563 to be the lead SNP at the GC locus in this young cohort (beta = −9.49 nmol/L per C allele), with the effect being substantially larger in infants receiving high-dose supplementation (30 µg/day) than in those receiving standard doses (10 µg/day), demonstrating that this variant modifies the biological response to vitamin D3.

In the Western Australian Raine Study99 Western Australian Raine Study
Anderson D et al. Genome-wide association study of vitamin D levels in children. Genes Immun, 2014
of 1,140 adolescents, rs1155563 reached genome-wide significance for vitamin D levels at age 14 (P = 3.9 × 10⁻⁹), providing independent replication in a different ancestry context.

Practical Actions

Because rs1155563 is a GC locus intronic variant, its primary clinical meaning mirrors that of the other GC tag SNPs: C allele carriers have lower total serum 25(OH)D on standard blood tests. The total/free vitamin D paradox still applies — reduced VDBP binding capacity means a larger fraction of vitamin D circulates free and bioavailable, partly buffering the apparent deficit on total 25(OH)D assays.

The most actionable implication is for supplementation dosing. The Kämpe 2019 pediatric GWAS showed that C allele effects were amplified under high-dose vitamin D3 — CC carriers receiving 30 µg/day still had lower 25(OH)D than non-carriers — suggesting these individuals need to supplement more consistently and may need to retest more frequently to verify dose adequacy.

Cholecalciferol (vitamin D3) taken with a fat-containing meal remains the standard recommendation. Monitoring total 25(OH)D annually, and interpreting borderline results in light of this genotype, gives the most actionable picture.

Interactions

rs1155563 is one of at least three partially independent signals at the GC locus — the others being rs2282679 (the lead GWAS tag) and the two missense variants rs7041 (Asp432Glu) and rs4588 (Thr436Lys). The combined haplotype structure of this locus determines the net VDBP isoform (Gc1s, Gc1f, or Gc2) and its associated transport capacity.

Upstream and downstream pathway partners compound GC effects. CYP2R1 (rs10741657) governs hepatic 25-hydroxylation; DHCR7/NADSYN1 (rs12785878) affects skin synthesis efficiency. Multi-SNP genetic scores combining GC variants with CYP2R1 and DHCR7 variants confer approximately two-fold increased vitamin D deficiency risk in susceptible haplotype carriers.

KIF2B rs116098458 — Spindle Assembly, Oocyte Quality, and Menopause Timing

The fidelity of chromosome segregation during female meiosis depends critically on the construction and function of the meiotic spindle — the microtubule apparatus that physically moves chromosomes to opposite poles of the oocyte. KIF2B (Kinesin Family Member 2B) is a member of the kinesin-13 family of microtubule depolymerases11 KIF2B (Kinesin Family Member 2B) is a member of the kinesin-13 family of microtubule depolymerases
KIF2B, also known as KIF2B/Kif2b, belongs to the kinesin-13 subfamily along with KIF2A and KIF2C (MCAK). Unlike motile kinesins, kinesin-13 members do not walk along microtubules — they bind microtubule ends and induce depolymerization to regulate spindle dynamics
. rs116098458 is an intronic variant in a long non-coding RNA (lncRNA) located antisense to KIF2B at chromosome 17q22, identified in the 2021 Ruth et al. Nature genome-wide association study as associated with variation in age at natural menopause.

The Mechanism

KIF2B is expressed at low levels in dividing cells but plays a distinct, non-redundant role in the early stages of spindle bipolarity. Depletion of KIF2B in human cells causes over 70% of mitoses to produce monopolar or disorganized spindles, with chromosome movement velocity reduced to approximately 20% of normal22 Depletion of KIF2B in human cells causes over 70% of mitoses to produce monopolar or disorganized spindles, with chromosome movement velocity reduced to approximately 20% of normal
Manning et al. 2007 Mol Biol Cell — the three kinesin-13 paralogs have distinct, non-overlapping functions during mitosis
. In oocytes, this matters acutely: meiotic spindle assembly occurs without centrosomes — an acentrosomal process that is particularly dependent on coordinated microtubule depolymerase activity to achieve bipolar organization from a disordered starting state. Failure of spindle bipolarity leads to chromosome segregation errors and oocyte aneuploidy33 chromosome segregation errors and oocyte aneuploidy
Aneuploidy — the wrong number of chromosomes in an egg — is the most common cause of embryo arrest, miscarriage, and implantation failure in IVF; it increases steeply with maternal age as spindle machinery becomes less efficient
.

The rs116098458 variant lies within a lncRNA transcribed antisense to KIF2B. Antisense lncRNAs can regulate their sense-strand gene partners through several mechanisms — transcriptional interference, chromatin remodeling, RNA–RNA base-pairing, or by acting as competing endogenous RNAs. Whether this specific variant alters KIF2B expression in oocytes or granulosa cells has not yet been characterized mechanistically. Its detection in a large-scale GWAS of menopause timing places it at a biologically plausible locus: the chromosome 17q22 region encompasses KIF2B and its antisense regulatory non-coding transcripts, and spindle assembly fidelity is a documented contributor to oocyte quality across the reproductive lifespan.

The Evidence

The primary evidence for rs116098458 comes from the Ruth, Day et al. 2021 Nature GWAS44 Ruth, Day et al. 2021 Nature GWAS
Genetic insights into biological mechanisms governing human ovarian ageing. Nature 596:393–397, 2021
, which analysed age at natural menopause in approximately 200,000 women of European ancestry and identified 290 genome-wide significant loci. The study implicates DNA damage response, meiotic recombination, and spindle assembly genes as the biological drivers of ovarian reserve depletion — the process that ultimately determines when menopause occurs. Loci at spindle assembly genes fit a coherent mechanistic model: accumulating errors in meiotic chromosome segregation over thousands of ovarian cycles deplete the functional follicle pool.

It is important to note that the T allele at rs116098458 is extremely rare: approximately 0.5% globally in gnomAD, essentially absent in individuals of European and East Asian ancestry, and present at approximately 1.3% in individuals of African ancestry. This frequency distribution means the GWAS association was likely driven by African-ancestry participants or discovered in trans-ethnic analyses. Effect size estimates from the primary publication are not publicly resolved for this specific variant at the time of writing.

The biological context of KIF2B in spindle assembly is mechanistically solid. The kinesin-13 proteins KIF2A, KIF2B, and KIF2C have distinct non-redundant roles in bipolar spindle establishment and chromosome movement55 kinesin-13 proteins KIF2A, KIF2B, and KIF2C have distinct non-redundant roles in bipolar spindle establishment and chromosome movement
Manning et al. 2007 Mol Biol Cell — each kinesin-13 paralog localizes to different spindle components and acts at different cell-cycle stages
. KIF2B acts specifically in early spindle assembly, and its absence cannot be compensated by KIF2A or MCAK.

Practical Implications

Because this variant is extremely rare and absent in most ancestry groups, its direct relevance is narrow. Individuals of African ancestry who carry the T allele have a slightly altered KIF2B-locus regulatory configuration whose net effect on oocyte spindle assembly and menopause timing is not yet precisely quantified. The general principle — that efficient meiotic spindle assembly is protective for long-term ovarian function — applies regardless of genotype and is supported by strong mechanistic data. Exposures that impair spindle microtubule dynamics (persistent organic pollutants with microtubule-binding properties, heavy metals such as cadmium, bisphenol A) are plausible environmental modifiers.

Interactions

rs10804920 (TP63 intronic variant): TP63 encodes the master DNA-damage checkpoint transcription factor in primordial follicle oocytes. While TP63 and KIF2B operate through mechanistically distinct pathways (DNA damage response vs. spindle assembly), both ultimately affect the rate of follicle pool depletion. Spindle defects during meiosis can generate DNA damage as a secondary consequence of chromosome missegregation, potentially activating the TAp63 checkpoint downstream. Women carrying risk alleles at both a spindle assembly locus and a DNA-damage checkpoint locus may have compounded vulnerability, though no published study has examined this interaction directly.

rs2307449 (BRSK1) and rs244715 (ZNF346/UIMC1): Other menopause-timing loci in the database that operate through DNA double-strand break repair (BRSK1) and BRCA1 complex function (UIMC1). Combined polygenic burden from spindle assembly and DNA repair loci could warrant proactive ovarian reserve assessment.

rs1165205

SLC17A3 SLC17A3 variant (NPT4)

Moderate Risk Factor

SLC17A3 rs1165205 — NPT4 and the Apical Urate Exit

Every day your kidneys filter and then selectively reclaim or discard uric acid as they fine-tune the serum urate setpoint. Most people know about the reabsorptive side of this equation — transporters like URAT111 URAT1
SLC22A12, encoded by the SLC22A12 gene, is the principal apical urate reabsorber in the kidney proximal tubule; it retrieves uric acid from the tubular lumen back into the bloodstream
that pull uric acid back from the urine into the blood. But secretion — actively pushing uric acid from the tubular cell into the urine — is equally important, and it depends on a distinct set of apical transporters on the urine-facing surface of the proximal tubule. One of these is NPT422 NPT4
sodium-phosphate transporter 4, encoded by SLC17A3 on chromosome 6p22; an apical multispecific organic anion efflux transporter that drives uric acid out of tubular cells into the urine for excretion
.

rs1165205 is an intronic variant in SLC17A3 that sits in the same gene cluster as SLC17A1 (encoding NPT1, the related apical urate efflux transporter). The A allele at this position is associated with higher serum uric acid in multiple populations and with reduced protection against gout. Because it is intronic, the variant likely influences NPT4 expression or mRNA processing rather than changing the protein directly — the secretory efficiency of the proximal tubule is tuned up or down depending on which version of the haplotype you carry.

The Mechanism

NPT4 operates as a voltage-driven organic anion efflux transporter33 voltage-driven organic anion efflux transporter
the membrane potential gradient across the apical membrane provides the driving force; an increase in extracellular potassium (simulating depolarization) enhances transport activity
. It transports uric acid from the interior of the proximal tubule cell into the tubular lumen, working alongside NPT1 (SLC17A1) to provide the secretory counterbalance to URAT1-mediated reabsorption. The net direction of urate movement — secretion or reabsorption — determines whether serum urate rises or falls.

When genetic variation at the rs1165205 locus impairs this secretory arm, the reabsorption-to-secretion ratio tilts toward retention, elevating the serum urate setpoint. This mechanism also explains why loop and thiazide diuretics44 loop and thiazide diuretics
furosemide and bumetanide are competitive inhibitors of NPT4-mediated urate transport; this provides the molecular explanation for the well-known phenomenon of diuretic-induced hyperuricemia
raise uric acid: they directly inhibit NPT4, reproducing the pharmacological equivalent of reduced-function genetic variants. Individuals with the A allele who also use diuretics are therefore doubly disadvantaged on urate secretory capacity.

The Evidence

rs1165205 was identified as a genome-wide significant locus for serum uric acid in the Dehghan et al. landmark GWAS published in The Lancet55 The Lancet
Dehghan A et al. Association of three genetic loci with uric acid concentration and risk of gout: a genome-wide association study. Lancet, 2008
, with p=3.3×10⁻²⁶ in the Framingham cohort (7,699 participants). The protective T allele was associated with OR 0.85 (95% CI 0.77–0.94, p=0.002) for gout risk in European participants. The association was not significant in African-ancestry ARIC participants, consistent with the very high A allele frequency in African populations (~88%) leaving little statistical power to detect effects.

The GWAS Catalog reports a pooled effect of approximately 0.09 unit decrease in serum urate per T allele (p=4×10⁻²⁹) across a total sample exceeding 15,000 participants. A multi-cohort analysis (n=4,492)66 multi-cohort analysis (n=4,492)
Brandstätter A et al. Sex and age interaction with genetic association of atherogenic uric acid concentrations. Atherosclerosis, 2010
confirmed the A allele as a component of the composite genetic risk score for elevated serum urate across European cohorts, with the SLC17A3 locus contributing modestly but consistently to the variance in serum uric acid levels.

Functional studies by Jutabha et al.77 Jutabha et al.
Jutabha P et al. Human sodium phosphate transporter 4 (hNPT4/SLC17A3) as a common renal secretory pathway for drugs and urate. J Biol Chem, 2010
established the biological plausibility: NPT4 is the exit route for uric acid from the proximal tubule apical membrane, and loss-of-function mutations in SLC17A3 found in hyperuricemia patients abolish urate efflux capacity in cell systems. The GWAS signal at rs1165205 is consistent with a regulatory variant that modulates this secretory capacity.

Practical Actions

Elevated urate from reduced renal secretory capacity is particularly responsive to reducing dietary purine load (since the kidneys cannot compensate by increasing secretion) and to avoiding substances that further inhibit NPT4, principally loop and thiazide diuretics. Individuals with the AA genotype on diuretics should specifically flag this to their prescriber, as the pharmacological inhibition of NPT4 compounds genetically reduced secretory capacity.

Vitamin C supplementation at 500–1,000 mg/day has modest evidence for lowering serum urate by approximately 0.5 mg/dL through competitive inhibition of renal urate reabsorption — a compensatory mechanism for reduced secretory capacity. Serum uric acid monitoring establishes the individual setpoint and guides when dietary changes need reinforcement.

Interactions

rs1165205 operates at the NPT4-mediated arm of renal urate secretion. The most clinically significant interaction is with rs2231142 in ABCG2, the breast cancer resistance protein that mediates the other major apical secretory route for urate. ABCG2 Q141K (rs2231142) reduces ABCG2 activity by approximately 50%, and when combined with reduced NPT4 capacity both secretory pathways are simultaneously compromised, producing substantially higher serum urate than either variant alone. This combination is particularly relevant in East Asian populations where both risk alleles are common.

rs1165205 is also in high linkage disequilibrium (r²=0.97) with rs1183201 in the adjacent SLC17A1 gene, which encodes NPT1 — the sibling apical urate secretory transporter. The two signals tag the same haplotype block spanning the SLC17A1–SLC17A3–SLC17A4 gene cluster. Individuals carrying risk alleles at both loci may experience additive reduction in total apical secretory capacity, though these SNPs are sufficiently correlated that both rarely appear as independent effects in the same model.

rs11652075

CARD14 CARD14 Arg820Trp (R820W)

Strong Risk Factor

CARD14 Arg820Trp — The Psoriasis Biologic Response Predictor

CARD14 is a scaffold protein11 scaffold protein
A non-enzymatic protein that organizes signaling complexes at cellular membranes
expressed primarily in keratinocytes — the epidermal cells that form the skin barrier. Under inflammatory conditions, CARD14 recruits BCL10 and MALT1 to form an activation complex that triggers NF-κB22 NF-κB
Nuclear Factor kappa-light-chain-enhancer of activated B cells; a master regulator of inflammation that switches on dozens of proinflammatory genes
signaling, producing TNF-α, IL-17, IL-23, and other cytokines that drive skin inflammation. The rs11652075 variant changes arginine to tryptophan at position 820 (p.Arg820Trp), altering the protein's regulatory domain in a way that affects how readily this NF-κB amplification cascade activates. This SNP carries two clinically distinct pieces of information: it slightly raises the odds of developing psoriasis, and it strongly predicts whether anti-TNF biologic therapy will achieve remission.

The Mechanism

Position 820 of CARD14 sits in the coiled-coil domain33 coiled-coil domain
A structural motif formed by two or more alpha-helices wound around each other; in CARD14, this domain controls protein-protein interaction with BCL10 and governs complex formation efficiency
that controls interaction with BCL10. The arginine-to-tryptophan substitution replaces a positively charged, hydrophilic residue with a bulky aromatic one. Functional assays by Jordan et al. demonstrated that CARD14 variants affecting this region produced NF-κB activation levels >2.5-fold above wild-type44 >2.5-fold above wild-type
Measured by luciferase reporter assay in HEK293 cells transfected with CARD14 expression constructs; two variants required TNF-α stimulation to show full activation
in cell-based assays. The variant also destroys a CpG dinucleotide methylation site55 destroys a CpG dinucleotide methylation site
CpG sites are targets for DNA methylation; when the C>T substitution occurs at a CpG, the cytosine can no longer be methylated, potentially altering epigenetic silencing of the region
, which may affect transcriptional regulation of nearby sequences. The biological net effect is a keratinocyte NF-κB pathway that operates at a mildly lower activation threshold, releasing more TNF-α and IL-17A under the same inflammatory stimuli — explaining both the modest psoriasis susceptibility and the robust response to TNF blockade.

The Evidence

The Jordan et al. 2012 study in American Journal of Human Genetics established the foundational evidence: across seven psoriasis cohorts with more than 6,000 cases and 4,000 controls66 more than 6,000 cases and 4,000 controls
Cohorts included European, North American, and Australasian ancestry groups
, rs11652075 reached genome-wide significance for psoriasis association (p=2.1×10⁻⁶). Notably, adjustment for the major HLA-Cw*0602 psoriasis risk allele strengthened the CARD14 signal, confirming it operates through an independent pathway.

A subsequent meta-analysis by Shi et al. pooled five studies totaling 32,807 cases and 45,458 controls77 five studies totaling 32,807 cases and 45,458 controls
Ancestry breakdown: European and East Asian populations both represented
and confirmed the T allele is protective against psoriasis (pooled OR=0.877, 95%CI 0.834–0.922, P<0.001), with consistent effects in both European (OR=0.883) and Asian (OR=0.872) populations. The effect is modest per allele — the primary clinical utility of this variant lies in pharmacogenomics, not risk stratification.

The pharmacogenomic evidence is more striking. Coto-Segura et al. sequenced the entire CARD14 gene in 116 psoriasis patients treated with TNF inhibitors88 116 psoriasis patients treated with TNF inhibitors
79 responders, 37 non-responders; response defined as PASI 75 reduction at week 24; anti-TNF agents included adalimumab, etanercept, and infliximab
. The CC genotype (no T allele) was significantly enriched among responders (OR=3.71, 95%CI 1.30–10.51, P=0.01). Patients with CC genotype were nearly four times more likely to achieve PASI 75 response by week 24. The mechanistic interpretation is straightforward: if CARD14-mediated NF-κB activity is the dominant driver of a patient's psoriasis, TNF-α blockade more effectively disrupts that pathway; patients with the protective T allele may have psoriasis driven by other mechanisms less responsive to anti-TNF therapy.

Practical Actions

For people without psoriasis, the CC genotype conveys only modest susceptibility — well under 1% absolute lifetime risk increase from this variant alone. For those who do develop psoriasis and are considering biologic therapy, CC status is meaningful: it identifies them as likely responders to adalimumab, etanercept, or infliximab before a single injection is given. For CT and TT carriers who develop psoriasis, anti-TNF therapy remains an option but response is less reliably predicted by this variant; IL-17 or IL-23 inhibitors may be comparably or more effective choices depending on other clinical factors.

The T allele's destruction of a CpG methylation site does not currently have specific management implications — no dietary or supplement intervention has been shown to compensate for epigenetic dysregulation at this locus.

Interactions

CARD14 rs11652075 operates in the PSORS2 psoriasis susceptibility locus and interacts with the major PSORS1 locus (HLA-Cw*0602). The Jordan 2012 data showed that conditioning on HLA-Cw*0602 status actually increased the CARD14 signal, suggesting the two loci contribute to psoriasis through partly non-overlapping mechanisms. For pharmacogenomics purposes, rs61751629 (another CARD14 coding variant) has been examined alongside rs11652075; the combination of CARD14 rare variants also predicted favorable anti-TNF response in the Coto-Segura dataset. These interactions are candidates for compound action assessment pending larger pharmacogenomic study replication.

ACADVL K382Q — A Pathogenic Fatty Acid Oxidation Variant and Its Carrier Implications

Inside every mitochondrion, a molecular relay strips energy from fat molecules two carbons at a time in a process called beta-oxidation11 beta-oxidation
Beta-oxidation is the main pathway by which cells convert dietary and stored fat into ATP, particularly during fasting, prolonged exercise, and periods of high energy demand
. Very long-chain acyl-CoA dehydrogenase (VLCAD), encoded by the ACADVL gene on chromosome 17, catalyzes the critical first step in this relay for fatty acids with chain lengths of 14–20 carbons — the very long-chain fats found abundantly in foods like meat, dairy, and vegetable oils. When VLCAD fails, these long-chain fats accumulate as toxic acylcarnitines and cannot be converted to energy, causing VLCAD deficiency (OMIM #201475)22 VLCAD deficiency (OMIM #201475)
An autosomal recessive inborn error of metabolism affecting 1 in 30,000–100,000 births; listed on all US newborn screening panels since the early 2000s
in its most severe forms.

The c.1144A>C variant (K382Q, p.Lys382Gln) substitutes a positively charged lysine for a neutral glutamine at position 382, within a region critical for FAD cofactor binding. This single change is sufficient to abolish VLCAD function. The variant was first identified by Souri et al. in 33 Souri M et al., Am J Hum Genet, 1996 in a patient with VLCAD deficiency and confirmed pathogenic through expression experiments: CHO cells transfected with K382Q cDNA showed no detectable VLCAD enzyme activity and produced a protein with abnormal dimer assembly — structurally broken rather than merely impaired.

The Mechanism

VLCAD functions as a homodimer anchored to the inner mitochondrial membrane. Lysine-382 sits within the FAD-binding domain; its positive charge is required for proper folding of the subunit interface. The K382Q substitution eliminates this charge, disrupting dimer assembly. The resulting misfolded protein is rapidly degraded — VLCAD activity drops to effectively zero in homozygous or compound heterozygous affected individuals 44 ClinGen ACADVL Expert Panel review (Dec 2022): functional data shows 19% residual activity via enzymatic assay; REVEL pathogenicity score 0.95. Carriers with one functional copy produce enough VLCAD to oxidize long-chain fats normally — heterozygotes are asymptomatic and have normal VLCAD activity in lymphocyte and fibroblast assays.

The Evidence

ClinVar classifies K382Q as likely pathogenic (Variation ID 1628), reviewed by the ClinGen ACADVL Variant Curation Expert Panel (4-star expert review, December 2022), with supporting submissions from Labcorp Genetics, Baylor Genetics, and Myriad Genetics. The OMIM allelic variant entry (609575.0008) documents the original 1996 pathogenic characterization.

In a study of Japanese VLCAD patients, K382Q accounted for 12.7% of mutant alleles in the pre-newborn-screening cohort55 K382Q accounted for 12.7% of mutant alleles in the pre-newborn-screening cohort
Osawa et al., Mol Genet Metab, 2022 — frequency fell to 3.1% in the expanded NBS cohort, reflecting ascertainment bias toward milder variants in a screened population
, establishing it as one of the more clinically significant ACADVL variants.

Functional fibroblast studies66 Functional fibroblast studies
Schiff et al., Mol Genet Metab, 2013
demonstrate that heterozygous carriers show normal VLCAD enzyme activity — confirming that a single functional copy is sufficient for normal fatty acid oxidation. This is the biochemical basis for why VLCAD deficiency follows autosomal recessive inheritance and why carrier individuals require no clinical management for themselves.

The variant is extremely rare in gnomAD v4 (1 in ~149,000 alleles, global), consistent with strong negative selection against loss-of-function VLCAD variants.

Practical Actions

Carriers of K382Q (AC genotype) are healthy and require no dietary restrictions or treatment for themselves. The clinical relevance is reproductive: if both partners carry a pathogenic ACADVL variant (from any combination of the 200+ known pathogenic/likely-pathogenic alleles), each pregnancy has a 25% chance of being affected. VLCAD-affected infants are identified by newborn screening via elevated C14:1 acylcarnitine; early dietary management (MCT-enriched, long-chain fat-restricted diet) prevents the most severe outcomes including cardiomyopathy and hypoglycemic crises.

Carriers may wish to confirm their partner's ACADVL carrier status through comprehensive gene sequencing — a panel approach is more informative than single-variant testing given the allelic heterogeneity of this gene.

Interactions

VLCAD deficiency (the disease) is a compound heterozygous or homozygous condition in most patients. K382Q has been documented in affected individuals in compound heterozygosity with other ACADVL variants such as the common p.V283A (c.848T>C) variant and various splice-site and truncating mutations. The severity of the resulting deficiency depends on the residual activity of the second allele: null+null combinations cause severe neonatal cardiac disease, while hypomorphic+null combinations (like many c.848T>C compound heterozygotes) cause milder myopathic presentations. K382Q is a functional null, so its severity in an affected child is determined primarily by the partner allele.

SLC28A2 Pro22Leu — The Ribavirin Transporter Variant and Anemia Risk

Every dose of ribavirin you swallow must cross the intestinal wall before it can reach the bloodstream, and the protein doing most of that work is CNT211 CNT2
Concentrative Nucleoside Transporter 2, a sodium-coupled symporter expressed on the apical (luminal) membrane of jejunal enterocytes that preferentially transports purine nucleosides — including ribavirin — against their concentration gradient, driven by the intestinal sodium electrochemical gradient
, encoded by SLC28A2 on chromosome 15q21.1. The rs11854484 variant changes a proline to leucine at position 22 of this transporter, altering how efficiently it accumulates ribavirin in the enterocyte and, downstream, how much drug ends up trapped in red blood cells.

The Mechanism

rs11854484 sits at GRCh38 position chr15:45,253,279 (NC_000015.10) within the coding sequence of SLC28A2. The C>T substitution converts Pro22Leu in the transporter's N-terminal cytoplasmic domain — a region that influences membrane trafficking and transporter turnover. Because SLC28A2 is on the plus strand, the plus-strand alleles match the coding strand directly: C is the reference, T is the Pro22Leu variant.

Ribavirin is a structural analog of guanosine that enters cells via CNT2 in the small intestine. Once inside enterocytes and erythrocytes, it is phosphorylated to ribavirin triphosphate, which cannot exit the cell easily. Red blood cells have no de-phosphorylation capacity, so ribavirin accumulates and disrupts membrane integrity — the direct cause of hemolytic anemia. The Pro22Leu variant appears to enhance CNT2 activity or increase transporter surface expression, leading to greater ribavirin uptake per unit of drug ingested.

Beyond ribavirin, CNT2 is the primary intestinal transporter for dietary and salvage-pathway purine nucleosides (adenosine, inosine, guanosine). Variants in this transporter may subtly shift purine nucleoside bioavailability and the balance between de-novo synthesis and the salvage pathway, though the clinical implications outside of drug therapy are not yet well characterised.

The Evidence

The clearest clinical evidence comes from a prospective cohort study of 216 Swiss HCV patients22 prospective cohort study of 216 Swiss HCV patients
Rau et al. J Hepatol 2013
treated with pegylated interferon-α plus ribavirin (with a subset receiving triple therapy including telaprevir or boceprevir). The TT genotype was associated with significantly higher weight-adjusted ribavirin serum levels at week 4 (p=0.02). Most strikingly, clinically significant anemia (haemoglobin drop requiring dose reduction) occurred in 56% of TT carriers versus only 33% of CC/CT carriers (p=0.006). In multivariate analysis, rs11854484 TT was an independent predictor of clinically significant anemia. Patients receiving triple therapy with protease inhibitors showed the same pattern, with TT genotype identifying a subgroup at substantially higher anaemia risk.

A pharmacokinetics study of 174 HCV-1 and HCV-4 Italian patients33 pharmacokinetics study of 174 HCV-1 and HCV-4 Italian patients
D'Avolio et al. Ther Drug Monit 2012
identified rs11854484 TT genotype as one of three independent predictors of sustained virological response (alongside IL28B rs8099917 and CYP27B1 rs4646536), with the number of "favourable" variant alleles correlating inversely with treatment failure — suggesting that higher ribavirin exposure in TT carriers improves antiviral efficacy while simultaneously raising the risk of anaemia.

A secondary analysis of 169 HCV-1 patients treated with standard peg-IFN/ribavirin44 secondary analysis of 169 HCV-1 patients treated with standard peg-IFN/ribavirin
Doehring et al. Pharmacogenet Genomics 2011
examined the full nucleoside transporter gene family (SLC28A2, SLC28A3, SLC29A1, SLC29A2) for ribavirin response; SLC28A2 variants (including rs11854484) modulated both ribavirin levels and haemoglobin outcomes.

Practical Actions

For TT homozygotes the key implication is raised awareness before any ribavirin-based treatment. Modern hepatitis C treatment is largely interferon-free and often ribavirin- free, but ribavirin is still used in some DAA (direct-acting antiviral) regimens for genotype 3 or retreatment cases. Knowing the TT genotype in advance allows clinicians to start at lower ribavirin doses, monitor haemoglobin more frequently (weekly for the first 4 weeks rather than every 2 weeks), and prepare for dose adjustment earlier — which preserves treatment completion rather than forcing discontinuation.

CT heterozygotes have intermediate ribavirin exposure and a modest elevation in anemia risk. Standard monitoring applies, but the genotype can inform shared decision-making with the treating hepatologist.

CC homozygotes have the reference transporter activity and the lowest anemia risk from this locus. Ribavirin dosing and monitoring follow standard protocols.

Interactions

rs11854484 (SLC28A2) operates within a genetic risk matrix for ribavirin-induced anaemia. The most important interaction partner is ITPA rs1127354, which encodes inosine triphosphatase — the enzyme that metabolises ribavirin triphosphate in erythrocytes. ITPA-deficient patients (rs1127354 CC) accumulate less ribavirin phosphate in RBCs and are paradoxically protected from haemolysis; the protective ITPA genotype partially counteracts the elevated ribavirin load in TT carriers (rs11854484). SLC28A3 rs56350726 and rs10868138 encode the related CNT3 transporter and have been associated with sustained virological response in some cohorts. An interaction between SLC28A2 TT and SLC28A3 variants would represent compounded transporter effects on ribavirin bioavailability.