Your heart's pumping power comes from millions of sarcomeres¹¹ sarcomeres
The sarcomere is the basic contractile unit of cardiac muscle. Each sarcomere contains an interdigitating lattice of thick (myosin) and thin (actin) filaments that slide past each other to produce contraction working in concert. The MYL3 gene encodes the essential myosin light chain (ELC), a structural component of the myosin thick filament that stabilizes the lever arm region and regulates the force and speed of cardiac muscle contraction. The Arg94His substitution — a single amino acid change replacing arginine with histidine at position 94 of the ELC — disrupts this regulation in a way that causes the heart muscle to thicken abnormally, a condition called hypertrophic cardiomyopathy (HCM).

HCM is the most common inherited cardiac disease, affecting approximately 1 in 500 people overall. However, MYL3 mutations are a rare cause of HCM — responsible for less than 1% of all sarcomeric HCM cases. The Arg94His variant in particular has been documented in a small number of families worldwide and is classified as pathogenic by eight independent clinical laboratories²² pathogenic by eight independent clinical laboratories
ClinVar VCV000031777: submissions from LabCorp Genetics, Women's Health and Genetics/LabCorp, Ambry Genetics, GeneDx, Mass General Brigham, NIH All of Us, Color Diagnostics, and CHEO Genetics.

The Mechanism

The MYL3 ELC wraps around the myosin heavy chain lever arm³³ lever arm
The lever arm is the portion of myosin that amplifies the small conformational change at the catalytic domain into a large stroke movement that pulls actin filaments. Light chains stabilize this lever arm against mechanical stress. Arginine-94 sits in a functionally critical region of the ELC. Replacing the positively charged arginine with histidine (which carries a smaller, uncharged side chain at physiological pH) is predicted to disrupt local protein folding and the ELC's interaction with the myosin heavy chain, altering the super relaxed state of myosin — a low-energy resting state that is disrupted in HCM-associated ELC mutations. The net result is hypercontractility of the sarcomere, leading over time to pathological left ventricular hypertrophy, diastolic dysfunction, and in some cases left ventricular outflow tract obstruction.

The Evidence

The primary evidence for MYL3 Arg94His comes from a Japanese family study⁴⁴ Japanese family study
Nomura et al., J Cardiol, 2016 (N=7 family members, 5 affected; followed by registry screening of 600 HCM patients) that used whole-exome sequencing combined with bioinformatic filtering to identify the variant. Among carriers in that family, disease penetrance was 88% — very high for an HCM gene. All clinically affected carriers showed asymmetric septal hypertrophy with a maximum left ventricular wall thickness of 18±3mm without outflow obstruction. Two additional carriers were found in a registry of 600 HCM patients, supporting that this is a genuine pathogenic variant rather than a private family mutation.

By contrast, other MYL3 variants (such as the p.Val79Ile variant⁵⁵ p.Val79Ile variant
Andersen et al., Biochem Res Int, 2012 — nine heterozygous carriers in a Danish family) show much lower penetrance (~40%) and late onset. This variability is consistent with the broader GeneReviews HCM data⁶⁶ GeneReviews HCM data
Cirino et al., updated 2025 — MYL3 overall penetrance ~32%, the lowest among sarcomeric HCM genes when all variants are pooled. The Arg94His variant appears to be at the higher end of MYL3 penetrance, based on the Nomura family data, though the evidence is limited to a small number of families. Eight independent clinical laboratories have independently classified it as pathogenic or likely pathogenic, providing strong multi-submitter validation.

In terms of population frequency, Arg94His is extremely rare: gnomAD v4 exomes identified only 7 T alleles among over 730,000 alleles examined globally, exclusively in European and Latino ancestry groups.

Practical Actions

For carriers of MYL3 Arg94His, the central priorities are: (1) confirm the diagnosis with a specialist, (2) establish regular cardiac surveillance, and (3) arrange cascade testing for first-degree relatives.

Unlike many common HCM genes (MYBPC3, MYH7) where genotype-specific prognosis is well-characterized, MYL3 Arg94His has been observed in relatively few families. Standard HCM management — echocardiography every 1–2 years, ambulatory ECG monitoring, exercise testing, and SCD risk stratification — applies. Drug therapy (beta-blockers, calcium channel blockers) addresses symptoms; septal reduction therapy (myectomy or alcohol ablation) is reserved for refractory obstructive cases. The newer selective cardiac myosin inhibitor mavacamten is approved for obstructive HCM and may be relevant if outflow obstruction develops.

Interactions

MYL3 variants have been observed to interact with other sarcomeric HCM genes. Individuals carrying two sarcomeric pathogenic variants (compound or digenic) show substantially worse outcomes — a hazard ratio of 7.5 for adverse events such as cardiac transplantation or ICD placement⁷⁷ hazard ratio of 7.5 for adverse events such as cardiac transplantation or ICD placement
Cirino et al., GeneReviews, 2025 — based on registry data across sarcomeric HCM genes compared to single-variant carriers. This means that if a carrier of MYL3 Arg94His also carries a pathogenic variant in MYBPC3, MYH7, TNNT2, or another sarcomeric gene, clinical severity is likely to be substantially greater.

MYL3 is an essential light chain (encoded by chromosome 3), distinct from the regulatory light chain encoded by MYL2. Variants in ACTC1 (cardiac actin, rs193922680) affect the thin filament partner of myosin and can produce overlapping HCM phenotypes through a different molecular mechanism.

Every time your body burns fat for energy during fasting, sleep, or prolonged exercise, it relies on a chain of enzymes inside the mitochondria. The first and rate-limiting step for very long-chain fatty acids (14–20 carbons) is carried out by very-long-chain acyl-CoA dehydrogenase¹¹ very-long-chain acyl-CoA dehydrogenase
VLCAD — encoded by ACADVL on chromosome 17p13. The enzyme sits on the inner mitochondrial membrane and initiates beta-oxidation of long-chain fatty acids, producing the acetyl-CoA and reduced electron carriers (FADH₂) that feed the Krebs cycle and electron transport chain. This SNP, rs200788251, represents a c.865G>A transition that replaces the glycine at position 289 with the much bulkier, positively charged arginine (p.Gly289Arg). Glycine-289 sits in a structurally conserved region of the enzyme and is under strong evolutionary constraint; the substitution reduces VLCAD activity to approximately 15% of normal in vitro.

The Mechanism

VLCAD²² VLCAD
very-long-chain acyl-CoA dehydrogenase, a homodimeric flavoprotein that removes two hydrogen atoms from the acyl-CoA thioester, generating a trans-2-enoyl-CoA and transferring electrons to electron transfer flavoprotein (ETF) for entry into the respiratory chain depends on precise folding of its active site for catalysis. The p.Gly289Arg substitution introduces a bulky, charged side chain where glycine's compact structure is required for proper protein folding. Functional studies demonstrate that fibroblasts from a compound heterozygous individual carrying the p.Gly289Arg allele alongside a second pathogenic variant showed no detectable enzyme activity and absent protein expression on western blot — indicating the substitution causes protein misfolding and accelerated degradation.

Because ACADVL operates as a homodimer, a single non-functional allele (the heterozygous carrier state) still produces enough functional enzyme dimers — with approximately 50% enzyme activity — to support normal fatty acid oxidation under ordinary dietary conditions. Disease emerges only when both alleles are compromised.

The Evidence

The variant is classified likely pathogenic by the ClinGen ACADVL Variant Curation Expert Panel³³ likely pathogenic by the ClinGen ACADVL Variant Curation Expert Panel
ClinVar VCV000370981, 4-star expert panel review, last evaluated June 2023, based on ACMG/AMP classification criteria. All twelve diagnostic laboratory submitters in ClinVar classify the variant as pathogenic or likely pathogenic. The glycine at codon 289 is conserved across vertebrates, computational tools (SIFT, PolyPhen-2) predict a deleterious effect, and in vitro functional studies confirm the consequence.

A US newborn screening study of 693 individuals⁴⁴ US newborn screening study of 693 individuals
Pena et al. Recurrent ACADVL molecular findings in individuals with a positive newborn screen for very long chain acyl-coA dehydrogenase (VLCAD) deficiency in the United States. Mol Genet Metab, 2016 identified p.Gly289Arg on six alleles from affected individuals, establishing it among the recurrent pathogenic alleles seen in clinical practice in the United States.

VLCAD deficiency occurs in approximately 1 in 30,000–100,000 births. Expanded newborn screening with acylcarnitine profiling — measuring C14:1 acylcarnitine⁵⁵ C14:1 acylcarnitine
tetradecadienoylcarnitine, the characteristic accumulation product of impaired VLCAD-mediated beta-oxidation; a C14:1 level above 1 µmol/L on dried blood spot strongly suggests VLCAD deficiency — now identifies most affected individuals presymptomatically.

Three phenotypes are documented in homozygous or compound heterozygous individuals: - Severe neonatal form: cardiomyopathy (hypertrophic or dilated), pericardial effusion, arrhythmia, and metabolic crisis in the first weeks of life — historically associated with significant mortality. - Infantile hepatic form: hypoketotic hypoglycemia, hepatomegaly, and liver dysfunction, typically presenting during intercurrent illness. - Late-onset myopathic form: the most common presentation in the current NBS era — episodic rhabdomyolysis triggered by prolonged exercise, fasting, or illness; muscle cramps, pain, dark urine (myoglobinuria), and markedly elevated creatine kinase (CK).

Practical Actions

Heterozygous carriers (AG genotype) are clinically asymptomatic. The primary practical significance of carrier status is reproductive: if both parents carry ACADVL pathogenic variants, each child has a 25% risk of VLCAD deficiency. Partner carrier testing before or during pregnancy is the key action.

For individuals with biallelic pathogenic variants (AA genotype, or compound heterozygous for two different ACADVL pathogenic alleles), the consensus-based nutrition management guidelines⁶⁶ consensus-based nutrition management guidelines
Vockley et al. 2021, Mol Genet Metab. PMID 33093005 recommend: avoidance of fasting (using age-appropriate maximum fasting intervals), restriction of long-chain fatty acid intake, MCT supplementation as an alternative fat source bypassing the enzymatic block, and carnitine supplementation to support acylcarnitine clearance. During exercise, consuming easily metabolized carbohydrates before and after activity significantly reduces rhabdomyolysis risk.

Interactions

This variant causes disease in the autosomal recessive setting: a second loss-of-function ACADVL variant on the other chromosome must be present. Most affected patients identified through newborn screening are compound heterozygous (two different pathogenic variants) rather than homozygous for p.Gly289Arg. Genotyping at this single SNP position will detect the AG carrier state but cannot determine whether a second pathogenic variant exists elsewhere in ACADVL — full gene sequencing is required for a complete clinical picture in symptomatic individuals or those with positive NBS results.

Alcohol dehydrogenase 1B (ADH1B) catalyzes the first step of alcohol metabolism, converting ethanol into acetaldehyde in the liver. Most people carry the common Arg370 form (ADH1B*1). The ADH1B*2 allele (His48Arg, rs1229984) is common in East Asian populations and encodes an enzyme roughly 100-fold more active than the common form. But there is a third functional variant — ADH1B*3 (Arg370Cys, rs2066702) — that is found almost exclusively in populations of African descent.

ADH1B*3 encodes a superactive enzyme with substantially higher ethanol oxidation kinetics than the common ADH1B*1 form. Like ADH1B*2, it accelerates the conversion of ethanol to acetaldehyde, producing an aversive physiological response that acts as a natural deterrent against heavy drinking. In African American populations, where this allele reaches frequencies of 15–25%, it is one of the strongest genetic predictors of alcohol consumption patterns ever identified.

This variant is functionally and epidemiologically distinct from rs1229984. The two alleles reside at different positions in the ADH1B protein (residue 48 versus residue 370), arise from independent mutational events, are distributed across different populations, and contribute independently to alcohol use disorder risk. A person of African descent may carry ADH1B*3 with no ADH1B*2 allele, and vice versa for East Asians.

The Mechanism

The rs2066702 A allele, on the plus (forward) strand of chromosome 4, corresponds to the Cys370 substitution in the ADH1B protein — the ADH1B*3 allele. The gene is located on the minus strand, so the A allele in genome files is the complement of the T allele in the coding sequence notation used in older literature.

At the protein level, replacing Arginine with Cysteine at position 370 alters the active-site geometry of the enzyme in a way that increases its catalytic efficiency for ethanol oxidation. The result is accelerated production of [acetaldehyde | A reactive aldehyde intermediate; classified as a Group 1 human carcinogen by the IARC; causes flushing, nausea, and tachycardia] after alcohol ingestion — similar in kind to what ADH1B*2 produces in East Asian populations, though the mechanistic basis (different residue, different kinetic parameters) is independent.

Because the downstream enzyme ALDH2, which clears acetaldehyde, operates at a fixed rate, any increase in acetaldehyde production from faster ADH1B activity creates a transient acetaldehyde surplus. Carriers of ADH1B*3 experience faster and more aversive responses to alcohol, which behaviorally reduces both the amount consumed and the likelihood of developing alcohol use disorder.

The Evidence

Genome-wide significance in African Americans:

The first genome-wide association study for maximum number of alcoholic drinks consumed in a 24-hour period in African Americans¹¹ The first genome-wide association study for maximum number of alcoholic drinks consumed in a 24-hour period in African Americans
Xu K et al. Genomewide Association Study for Maximum Number of Alcoholic Drinks in European Americans and African Americans. Alcohol and Alcoholism, 2015 identified rs2066702 as the peak genome-wide significant signal in African Americans (p = 2.50×10⁻¹⁰). Eight SNPs in the region on chromosome 4 reached significance, all tagging the ADH1B*3 haplotype.

A subsequent GWAS of maximum habitual alcohol intake in 17,029 African American veterans from the VA Million Veteran Program²² A subsequent GWAS of maximum habitual alcohol intake in 17,029 African American veterans from the VA Million Veteran Program
Gelernter J et al. Genome-wide Association Study of Maximum Habitual Alcohol Intake in >140,000 U.S. European and African American Veterans. Biological Psychiatry, 2019 confirmed rs2066702 as the lead locus for African Americans (p = 2.3×10⁻¹²), a far stronger signal than any other variant identified in that ancestry group.

Clinical alcohol use metrics:

In a large electronic health record study of 57,677 African American veterans³³ In a large electronic health record study of 57,677 African American veterans
Justice AC et al. Validating Harmful Alcohol Use as a Phenotype for Genetic Discovery Using Phosphatidylethanol and a Polymorphism in ADH1B. Alcohol and Alcoholism, 2018, 34.3% of participants carried at least one A allele (minor allele frequency 19.0%). Carrying the protective A allele was associated with approximately half the odds of high AUDIT-C scores (aOR = 0.54) and an aOR of 0.51 when combining AUDIT-C with ICD-coded alcohol use disorder — demonstrating that the protective effect is clinically meaningful and detectable across multiple measurement approaches.

AUD diagnostic criteria:

Among African American participants in a detailed AUD phenotyping study⁴⁴ Among African American participants in a detailed AUD phenotyping study
Hart AB et al. Which alcohol use disorder criteria contribute to the association of ADH1B with alcohol dependence? Addiction Biology, 2016, individuals homozygous for the major (non-protective) allele GG endorsed significantly more DSM-IV and DSM-5 AUD criteria (p = 1.9×10⁻⁹). The criterion most strongly linked to this variant was tolerance — the need for increasing amounts of alcohol to achieve the desired effect — suggesting that ADH1B*3 carriers develop tolerance more slowly, consistent with the variant's faster acetaldehyde production making high-dose alcohol less pleasant.

Prenatal exposure and developmental protection:

A longitudinal study in African American families⁵⁵ A longitudinal study in African American families
Dodge NC et al. Protective effects of the alcohol dehydrogenase-ADH1B*3 allele on attention and behavior problems in adolescents exposed to alcohol during pregnancy. Alcoholism: Clinical and Experimental Research, 2014 found that maternal ADH1B*3 carrier status shielded adolescent offspring from behavioral and attention consequences of prenatal alcohol exposure, with allele frequencies of 17.6% in mothers and 21.0% in adolescents — consistent with the expected 15–20% frequency in African American populations. The proposed mechanism is that faster maternal alcohol metabolism reduces peak blood alcohol concentration reaching the fetus.

Practical Actions

For GG carriers (most common in non-African populations): Without the ADH1B*3 protective allele, you lack one of the biological deterrents against heavy drinking. If you are of African descent, your GG genotype means you metabolize alcohol at the common rate — higher-tolerance drinking is biologically more accessible to you, which removes a natural brake on consumption. This is relevant when considering whether your drinking patterns are influenced by behavioral versus genetic factors.

For AG and AA carriers (primarily found in people of African descent): Your ADH1B*3 allele accelerates alcohol-to-acetaldehyde conversion, producing faster aversive responses to alcohol. This is one of the strongest inherited protective factors against alcohol use disorder in populations of African ancestry. The biological deterrent is real — but, as with ADH1B*2 in East Asian populations, it can be overridden socially. Environmental factors (peer norms, social offers) can diminish the genetic protection when drinking pressures are high.

Interactions

ADH1B*3 (rs2066702) and ADH1B*2 (rs1229984, His48Arg) reside in the same gene but at different protein positions. They arise independently and distribute across different populations. In African American populations where both alleles may occasionally co-occur, combined diplotype analysis is informative. However, the primary co-variant to consider alongside rs2066702 in African-ancestry individuals is ALDH2 (rs671), which controls the downstream clearance of the acetaldehyde that ADH1B*3 produces more rapidly.

ADH1C variants (rs1693482, rs698) are pathway partners that also affect the rate of alcohol oxidation. Combined diplotype analyses have documented that ADH1B and ADH1C allele combinations affect liver disease risk and alcohol metabolism outcomes beyond either variant alone.

The gastric epithelium is under constant assault — acid, pepsin, ingested pathogens, and the chronic coloniser Helicobacter pylori¹¹ Helicobacter pylori
A gram-negative bacterium that infects the stomach lining of roughly half the global population; the primary causal agent of peptic ulcers and a major risk factor for gastric cancer. The first line of defence is a thick mucus layer anchored by mucin glycoproteins, the most important of which is MUC1²² MUC1
Mucin 1, a high-molecular-weight, heavily O-glycosylated transmembrane protein expressed on the apical surface of gastric epithelial cells; it provides both steric hindrance against pathogens and acts as a releasable decoy. rs2070803 is a regulatory variant near the MUC1 gene on chromosome 1q22 that influences how much of this protective protein the stomach produces. Carrying one or two copies of the A allele is associated with lower MUC1 surface expression — and, in some of the largest genetic studies ever conducted on gastric cancer, a meaningfully elevated risk of developing diffuse-type gastric cancer³³ diffuse-type gastric cancer
One of two major histological subtypes of gastric cancer (alongside intestinal-type); diffuse-type spreads through the stomach wall without forming a distinct mass and carries a worse prognosis.

The Mechanism

rs2070803 maps approximately 585 base pairs upstream of the MUC1 transcription start site on the reverse strand, placing it in the gene's regulatory region. It is in strong linkage disequilibrium⁴⁴ linkage disequilibrium
Non-random co-inheritance of nearby variants; alleles at linked sites are more often inherited together than expected by chance with the functionally characterised variant rs4072037, which directly alters MUC1 promoter activity and controls the ratio of major MUC1 splice variants produced in the gastric epithelium.

Carriers of the A allele at rs2070803 produce less MUC1 protein at the gastric mucosal surface. This reduction compromises the epithelium in two ways. First, MUC1's extracellular domain normally acts as a steric barrier against bacterial adhesion — its dense glycan coat is physically too large to allow most bacteria to reach the underlying cell membrane. Second, MUC1 functions as a releasable decoy: shed extracellular domain fragments bind pathogen adhesins and carry them away from the epithelium. With less surface MUC1, both mechanisms are weakened.

H. pylori further compounds genetically reduced expression. The virulence factor CagA⁵⁵ CagA
Cytotoxin-associated gene A protein, injected directly into gastric epithelial cells by H. pylori; once inside, it is tyrosine-phosphorylated and disrupts multiple host signalling pathways is injected into host cells, where it undergoes phosphorylation and promotes binding of MUC1's cytoplasmic tail to β-catenin. The resulting MUC1–β-catenin complex translocates to the nucleus, activating proliferative gene programmes including cyclin-D1. Carriers who already have lower baseline MUC1 expression due to the A allele are therefore more vulnerable to this CagA-driven oncogenic cascade.

The Evidence

The primary evidence comes from a genome-wide association study⁶⁶ genome-wide association study
Saeki N et al. A functional single nucleotide polymorphism in mucin 1, at chromosome 1q22, determines susceptibility to diffuse-type gastric cancer. Gastroenterology, 2011 conducted across three panels of Japanese and Korean individuals (discovery: 606 cases/1,264 controls; validation 1: 304 cases/1,465 controls; validation 2: 452 cases/372 controls). rs2070803 reached genome-wide significance with p = 4.33 × 10⁻¹³ and a meta-analytic odds ratio of 1.71 for diffuse-type gastric cancer. Notably, the association was specific to diffuse-type and not intestinal-type gastric cancer.

A complementary case-control study⁷⁷ case-control study
Xu L et al. Risk of gastric cancer is associated with the MUC1 568 A/G polymorphism. Int J Oncol, 2009 directly measured MUC1 protein in gastric tissue: AA genotype carriers showed significantly lower MUC1 immunostaining (r = −0.179, p = 0.004) and a 1.81-fold increased gastric cancer risk, providing the functional link between genotype and reduced mucosal protection.

A systematic review and meta-analysis⁸⁸ systematic review and meta-analysis
Giraldi L et al. MUC1, MUC5AC, and MUC6 polymorphisms, Helicobacter pylori infection, and gastric cancer. Eur J Cancer Prev, 2018 of 21 studies confirmed the protective effect of the G allele (OR 0.66, 95% CI 0.57–0.78 for dominant model AG/GG vs. AA), with consistent findings across Asian (OR 0.73) and White European (OR 0.48) populations.

The risk is further amplified when rs2070803 co-occurs with the PSCA risk variant rs2294008. Carriers of risk alleles at both loci face a dramatically higher combined gastric cancer risk, as these two genes act in partially overlapping mucosal defence pathways.

Practical Implications

This variant is a genuine gastric cancer susceptibility signal — not a certainty, but a meaningful elevation in lifetime risk that is actionable. The most direct intervention is targeted H. pylori screening, since the genetic risk is substantially mediated through the bacteria's ability to exploit weakened mucosal defences. Current evidence suggests that eradicating H. pylori reduces gastric cancer incidence by roughly 35–45% in infected individuals, and for A allele carriers the absolute benefit may be larger because the residual mucosal vulnerability is genetically amplified.

Dietary and supplement strategies that support mucosal integrity — specifically those with evidence in the gastric context — complement surveillance without replacing it. Antral gastroscopy⁹⁹ gastroscopy
Direct visual inspection of the stomach lining using an endoscope surveillance intervals recommended by gastroenterologists should be followed by anyone who also has chronic gastritis, intestinal metaplasia, or a family history of gastric cancer.

Interactions

rs2070803 is in strong linkage disequilibrium with rs4072037, the functionally characterised MUC1 promoter variant. These two SNPs tag the same haplotype in most studies and their associations are not independent — the risk captured by rs2070803 is largely the same risk captured by rs4072037.

The PSCA variant rs2294008 represents a separate pathway (prostate stem cell antigen, expressed in gastric epithelium) that interacts epistatically with MUC1 variants. Carriers of both MUC1 and PSCA risk alleles face substantially higher combined risk than either alone, as documented in the original Saeki et al. GWAS.

ABCG2 — also called Breast Cancer Resistance Protein (BCRP) — is one of the most versatile efflux pumps in the body. It sits at the apical surface of intestinal cells, renal tubules, and the blood–brain barrier, where it exports a wide range of substances including uric acid¹¹ uric acid
a metabolic waste product of purine breakdown, porphyrins, anticancer drugs, antibiotics, and antiretroviral medications. The rs2231137 variant encodes a valine-to-methionine substitution at position 12 of the protein — a residue in the first transmembrane segment that influences how ABCG2 folds and reaches the cell membrane.

The Mechanism

The V12M change (T allele on the plus strand; the protein letter change is coded on the minus strand as G>A at nucleotide 34 of the coding sequence) sits at the very beginning of the ABCG2 protein's first predicted transmembrane helix. One cell model study found ²² Kamo et al., Intl J Cancer, 2004 that V12M disrupts apical plasma membrane localization — the protein is synthesized normally but fails to traffic to the correct compartment. Cells expressing V12M showed drug efflux less than one-tenth that of wild-type. However, a second study using membrane vesicles from HEK293 cells ³³ Kondo et al., Pharm Res, 2004 found transport activity for estrone sulfate, DHEAS, and methotrexate was comparable to wild-type. The discrepancy likely reflects cell-type and assay differences, making V12M's intrinsic functional impact harder to pin down than the better-characterised Q141K variant (rs2231142), which consistently reduces urate transport by ~53%.

The Evidence

Population genetics paints an interesting picture. The T allele (Met12) occurs at about 4% in Europeans but reaches 27% in East Asians and 23% in Latino/Admixed American populations — a striking frequency difference suggesting different selective pressures across ancestries.

A meta-analysis of 52 studies⁴⁴ meta-analysis of 52 studies
Shi et al., BMC Med Genet, 2020 found that carrying one or two T alleles at rs2231137 was associated with 36–57% lower odds of gout compared to CC homozygotes. This protective direction for gout is somewhat counterintuitive if V12M reduces urate transport, and may reflect incomplete LD with nearby protective haplotypes, population stratification, or context-dependent functional effects.

In drug-handling contexts, the picture is different. A study of 149 HIV-infected patients ⁵⁵ Rodrigues et al., Mol Genet Genomic Med, 2024 found that the T allele and the heterozygous CT genotype were associated with increased risk of antiretroviral-induced hepatotoxicity (OR 1.50 for the allele; haplotype OR up to 2.49 in combination with rs2231142). In a small study of 93 Iranian children, ⁶⁶ Barakat et al., Seizure, 2022 each T allele was associated with approximately 2.4-fold higher odds of drug-resistant epilepsy, suggesting altered antiepileptic drug transport.

An epistatic interaction study in 4,914 Chinese participants ⁷⁷ Tang et al., Hereditas, 2020 found that rs2231137 interacts with the PKD2 variant rs2728121 to influence urate levels and gout risk, with stronger effects in males (Pint = 0.004).

Practical Actions

The main clinical implications of this variant are in the context of specific drugs whose clearance depends on ABCG2 — particularly antiretrovirals used in HIV treatment and certain antiepileptic drugs. Carriers of the T allele, especially those also carrying risk alleles at the companion ABCG2 variant rs2231142 (Q141K), may have a compounded reduction in ABCG2 transport capacity that warrants monitoring of drug levels and liver enzymes when starting hepatically-cleared medications.

Interactions

The most clinically important interaction is with rs2231142 (Q141K), the dominant functional variant in ABCG2. Q141K alone reduces urate transport by ~53% and drives most of the ABCG2-gout association. Carrying T alleles at both rs2231137 (V12M) and rs2231142 (Q141K) may compound transporter impairment beyond either variant alone — a haplotype combination studied for drug resistance in cancer cells (PMID 40806557). The 12M/141K haplotype showed 10.7-fold resistance to the anticancer drug mitoxantrone, substantially higher than either variant alone.

There is also a documented epistatic interaction with rs2728121 in the PKD2 gene, which co-regulates urate levels and gout susceptibility through separate renal mechanisms.

Every cell in your body depends on a continuous supply of nucleosides — the molecular building blocks of DNA and RNA — and specialized transport proteins ferry these molecules across cell membranes. One such transporter, CNT2¹¹ CNT2
Concentrative Nucleoside Transporter 2, encoded by SLC28A2 on chromosome 15; the "concentrative" label reflects that it moves nucleosides against their concentration gradient, powered by co-transported sodium ions, sits at a critical junction in the intestinal wall and kidney tubules, governing how efficiently purines and certain nucleoside-based drugs enter and leave the body. The rs2413775 variant sits 146 base pairs upstream of the SLC28A2 coding sequence, inside the gene's proximal promoter, and alters how much CNT2 protein the liver and kidney produce.

The Mechanism

The variant creates a T-to-A change at position −146 relative to the SLC28A2 transcription start site, directly within a binding site recognized by hepatic nuclear factor 1²² hepatic nuclear factor 1
HNF1α and HNF1β are transcription factors that control gene expression in liver, kidney, and intestine; they bind specific DNA sequences and recruit the RNA polymerase machinery (HNF1). The A allele enhances HNF1α and HNF1β binding affinity at this site — demonstrated by electrophoretic mobility shift assays showing superior binding of -146A oligonucleotides at all competitor concentrations tested. The result is a gain-of-function³³ gain-of-function
The A allele increases, not decreases, CNT2 expression; this is the globally predominant allele (~62% globally), so "higher expression" is the population norm effect: Yee et al. (2009)⁴⁴ Yee et al. (2009)
Yee SW et al. Identification and characterization of proximal promoter polymorphisms in the human concentrative nucleoside transporter 2 (SLC28A2). J Pharmacol Exp Ther. 2009;328(3):699-707 showed that constructs carrying the A allele drove 1.8-fold greater luciferase activity than T-allele constructs in HepG2 hepatoblastoma cells (p < 0.001), with the effect confirmed in renal and colorectal cell lines and in mouse liver in vivo. Individuals carrying the T allele — the GRCh38 reference base but the global minor allele — produce less CNT2 protein per cell.

The Evidence

The foundational functional work was published in 2009 by the Pharmacogenomics of Membrane Transporters consortium⁵⁵ Pharmacogenomics of Membrane Transporters consortium
PMT — a UCSF-led consortium dedicated to characterizing how genetic variation in drug transporters affects pharmacokinetics in diverse populations. Yee et al. (2009)⁶⁶ Yee et al. (2009) sequenced the SLC28A2 promoter in 272 individuals across five ethnic groups and identified rs2413775 as the only promoter variant with an allele frequency exceeding 20% in all populations studied. In a parallel analysis, Hesselson et al. (2009)⁷⁷ Hesselson et al. (2009)
Hesselson SE et al. Genetic variation in the proximal promoter of ABC and SLC superfamilies: liver and kidney specific expression and promoter activity predict variation. PLoS ONE. 2009;4(9):e6942 confirmed that the T allele was associated with lower promoter activity in hepatic tissue across a broader survey of SLC family promoter variants.

The clinical pharmacogenomic relevance of SLC28A2 expression level has been demonstrated in studies of ribavirin, a purine nucleoside analog used to treat hepatitis C. CNT2 mediates intestinal absorption and renal tubular reabsorption of ribavirin; higher CNT2 expression translates to greater intracellular drug accumulation. A multicenter study by Ampuero et al. (2015)⁸⁸ Ampuero et al. (2015)
Ampuero J et al. Role of ITPA and SLC28A2 genes in the prediction of anaemia associated with protease inhibitor plus ribavirin and peginterferon in hepatitis C treatment. J Clin Virol. 2015;68:8-13 found that the related SLC28A2 variant rs11854484 TT genotype independently predicted clinically significant ribavirin-induced anemia (OR 2.33, 95% CI 1.10–4.95, p = 0.027) — demonstrating that SLC28A2 expression differences have measurable hematological consequences during nucleoside analog therapy.

Beyond drug transport, SLC28A2 plays a role in purine nucleoside homeostasis. Zhou et al. (2019)⁹⁹ Zhou et al. (2019)
Zhou Z et al. Common variants in the SLC28A2 gene are associated with serum uric acid level and hyperuricemia and gout in Han Chinese. Hereditas. 2019;156:7 identified multiple SLC28A2 variants significantly associated with serum uric acid levels and gout risk in 4,015 Han Chinese participants, confirming CNT2's role in purine clearance and urate homeostasis.

Clinical translation of rs2413775 specifically awaits prospective pharmacokinetic/pharmacodynamic studies — the existing evidence is mechanistically compelling but direct outcome data for this exact promoter variant remain limited, warranting an emerging evidence classification.

Practical Actions

For individuals receiving nucleoside analog drugs that are CNT2 substrates — primarily ribavirin and didanosine — the T allele may signal reduced transporter expression and potentially lower intracellular drug accumulation. This could affect both therapeutic efficacy (less drug reaching target cells) and tolerability (potentially less drug accumulation in red blood cells, which drives ribavirin- induced hemolytic anemia). The clinical significance for any individual patient depends on the specific drug, dose, and co-administered agents.

For general health, CNT2 also transports dietary purine nucleosides across the intestinal epithelium. Reduced CNT2 expression in T-allele carriers may modestly affect purine absorption from high-purine foods (red meat, organ meat, shellfish, beer), though the magnitude of this effect on serum uric acid has not been characterized for rs2413775 specifically.

Interactions

CNT2 works alongside the equilibrative nucleoside transporters ENT1 (SLC29A1) and ENT2 (SLC29A2) and the concentrative transporter CNT3 (SLC28A3) to govern nucleoside flux. Studies of ribavirin pharmacogenomics consistently include SLC29A1 variants alongside SLC28A2, as the two transporter families serve partly overlapping substrates in different tissue compartments.

The coding-region variant rs11854484 in SLC28A2 has been more extensively studied in clinical cohorts for ribavirin-induced anemia and treatment outcomes than rs2413775 itself, and these variants are not in strong linkage disequilibrium — they tag independent aspects of CNT2 biology (promoter expression level vs. protein function).

CD226 (also known as DNAM-1, DNAX Accessory Molecule-1)¹¹ CD226 (also known as DNAM-1, DNAX Accessory Molecule-1)
A co-stimulatory receptor expressed on T cells, natural killer cells, and myeloid cells that delivers activating signals when it binds PVR (CD155) or Nectin-2 (CD112) on antigen-presenting cells and stressed tissues is one of the few genes in the human genome where both coding and non-coding genetic variation independently contribute to autoimmune risk. rs2469434 sits deep within an intron of CD226 on chromosome 18q22, roughly 12,000 base pairs from the functionally well-characterised coding variant rs763361 (Gly307Ser). The two variants capture distinct signals at the same locus: rs763361 changes the protein's cytoplasmic signalling domain, while rs2469434 is thought to act through regulatory mechanisms — altering CD226 expression levels in immune cells rather than protein structure.

The Mechanism

rs2469434 maps to the notation c.728-3564 in the canonical CD226 transcript (NM_006566.4), placing it 3,564 nucleotides upstream of a major exon boundary — far from splice sites and beyond the reach of conventional splicing effects. The most plausible functional interpretation is that this variant lies within an intronic enhancer or silencer element²² enhancer or silencer element
Regulatory DNA sequences within introns that bind transcription factors and modulate gene expression levels without altering the protein sequence, modulating CD226 transcription in immune cell types. The GWAS signal for altered neutrophil count and total white blood cell count — two haematopoietic phenotypes with high expression in circulating immune cells — is consistent with a cis-regulatory effect on CD226 expression during immune cell maturation and differentiation.

CD226 itself delivers co-stimulatory signals that lower the activation threshold of both innate NK cells and adaptive T cells. Higher expression of the receptor (or dysregulated expression timing) translates into lower thresholds for immune cell activation — the cellular substrate for autoimmune tissue damage. The C risk allele at rs2469434 is associated with increased neutrophil count and elevated WBC in the largest multi-ethnic blood cell GWAS to date (Chen MH et al. 2020, 746,667 individuals)³³ (Chen MH et al. 2020, 746,667 individuals), consistent with a CD226 regulatory effect on myeloid cell development or survival.

The Evidence

The rheumatoid arthritis association at rs2469434 was established by the landmark Okada et al. 2014 RA GWAS meta-analysis⁴⁴ Okada et al. 2014 RA GWAS meta-analysis
Okada Y et al. Genetics of rheumatoid arthritis contributes to biology and drug discovery. Nature, 2014; >100,000 subjects, 29,880 RA cases, 73,758 controls, which mapped 42 novel RA risk loci including the CD226 18q22 region. The C allele shows an odds ratio of approximately 1.07 (95% CI 1.05–1.10) in European populations and 1.11 (95% CI 1.07–1.15) in East Asian populations, with genome-wide significance (P=9×10⁻¹⁰ and 1×10⁻⁸ respectively). These effect sizes are modest — consistent with the polygenic architecture of RA, where each individual variant contributes a small increment to cumulative risk.

The multiple sclerosis association comes from the International MS Genetics Consortium 2019 study⁵⁵ International MS Genetics Consortium 2019 study
Published in Science, 47,429 MS cases and 68,374 controls; established 200 autosomal susceptibility variants outside the MHC, which found rs2469434-C associated with MS susceptibility at OR 1.055 (P=3×10⁻⁸).

Importantly, the haematopoietic phenotype associations — neutrophil count (P=10⁻¹⁶) and total WBC count (P=5×10⁻¹³) from Chen MH et al. 2020 — provide mechanistic support for the autoimmune risk associations. Both phenotypes are cellular traits with well-known links to immune activation and autoimmune disease activity: elevated neutrophil counts are observed in active RA and MS flares, and circulating neutrophil levels are an intermediate phenotype between CD226 regulatory variants and clinical autoimmune disease expression.

Practical Implications

The per-allele ORs from this intronic variant (1.07–1.11 for RA, 1.055 for MS) are modest but real, and they act independently of the better-characterised rs763361 (Gly307Ser) coding variant. For CC homozygotes carrying two C alleles, the cumulative effect is approximately the square of the per-allele OR — roughly 1.14–1.23 × the risk of TT homozygotes, depending on the disease. In the context of RA and MS, which are both T-cell-mediated diseases with well-defined early biomarkers, monitoring for early signs is the primary actionable response.

Interactions

rs2469434 and rs763361 both map to the CD226 locus and independently associate with overlapping autoimmune diseases. Their independence as signals means a person can carry risk alleles at both variants, producing a compound effect at the CD226 locus beyond what either variant predicts alone. More broadly, the CD226 18q22 locus interacts pathogenically with other T-cell checkpoint genes: PTPN22 rs2476601 (LYP phosphatase), CTLA4 rs3087243 and rs231775 (inhibitory co-receptor), and IL2RA rs2104286 (IL-2 receptor alpha). All four genes converge on the T-cell activation/tolerance decision, and elevated CD226 expression from rs2469434 would amplify the effect of any co-existing checkpoint gene risk alleles.

The CNR2 gene encodes the cannabinoid receptor 2¹¹ cannabinoid receptor 2
CB2, a G-protein coupled receptor expressed predominantly on immune cells including macrophages, B cells, T cells, and natural killer cells. Unlike its brain-dominant counterpart CB1, CB2 is the immune system's endocannabinoid tuning knob — when activated by endogenous ligands such as 2-arachidonoylglycerol (2-AG) and anandamide, it suppresses pro-inflammatory cytokine release, reduces immune cell migration into tissues, and dampens excessive immune activation. rs2501401 is an intronic variant in CNR2 that functions as a cis-eQTL²² cis-eQTL
expression quantitative trait locus — a variant that regulates how much of the nearby gene's protein gets made, with the minor A allele associated with lower CB2 expression in whole blood (G allele NES = +0.128, p = 3.02 × 10⁻⁷ in GTEx v10, meaning each G copy raises expression).

The Mechanism

Although rs2501401 does not change any amino acid in the CB2 protein, its intronic location likely affects regulatory elements — enhancers, splice regulatory sequences, or chromatin accessibility signals — that control transcriptional output of CNR2 in immune cells. The G allele (carried by ~80% of Europeans and Africans) is associated with higher CB2 expression. The A allele (minor, ~20% in Europeans but ~65% in East Asians) correlates with lower CB2 levels in blood, which may reduce the brake on inflammatory signalling. CB2 activation normally suppresses NF-κB and MAPK pathways³³ suppresses NF-κB and MAPK pathways
the molecular switches that drive cytokine production in macrophages and T cells, so lower expression could leave immune cells in a more pro-inflammatory baseline state.

The Evidence

Direct clinical association studies for rs2501401 specifically have not been published. The evidence connecting CNR2 variation to autoimmune risk comes primarily from studies of the functional coding variant Q63R (rs35761398).

A 2025 Iranian cohort study (n=240) found that CNR2 Q63R RR carriers showed >2.5-fold increased rheumatoid arthritis risk⁴⁴ CNR2 Q63R RR carriers showed >2.5-fold increased rheumatoid arthritis risk
Nateghi et al. International Journal of Genomics, 2025 under multiple inheritance models. A 2020 case-control study (n=200)⁵⁵ 2020 case-control study (n=200) in Iran reported OR 2.70 (95% CI 1.47–4.97, p=0.001) for multiple sclerosis under the dominant model for the same Q63R variant. Studies in pediatric immune thrombocytopenic purpura similarly found OR ~2.35 for the minor allele, replicating across Egyptian and Italian cohorts⁶⁶ replicating across Egyptian and Italian cohorts. In chronic HCV infection, the variant was the only independent predictor of immune-mediated disorders⁷⁷ only independent predictor of immune-mediated disorders
cryoglobulinemia, autoimmune thyroiditis, ANA positivity, p=0.005 (49.4% vs 24.1% in controls).

The eQTL relationship for rs2501401 is confirmed in GTEx whole-blood data (p = 3.02 × 10⁻⁷), establishing that this intronic variant meaningfully influences CNR2 transcript levels. Mechanistically, CB2 has been shown to directly suppress T-cell JAK/STAT signalling⁸⁸ directly suppress T-cell JAK/STAT signalling
Xiong et al. Signal Transduction and Targeted Therapy, 2022, and increased CB2 expression is seen in inflamed intestinal tissue of children with IBD, suggesting a compensatory anti-inflammatory role.

The weight of evidence places this variant's evidence level at emerging: the eQTL is statistically robust, but a direct clinical trial or large GWAS using rs2501401 as the index SNP has not been published.

Practical Actions

For people carrying the A allele (lower CB2 expression), the primary strategy is supporting endocannabinoid tone through precursor nutrients and direct CB2 agonism, combined with early monitoring of inflammatory markers. Omega-3 fatty acids (EPA and DHA) are substrates for 2-arachidonoylglycerol (2-AG) synthesis, the most abundant endogenous CB2 ligand, making omega-3 sufficiency especially relevant when receptor expression is reduced. Palmitoylethanolamide (PEA) — an endogenous CB2 agonist — has clinical trial evidence in inflammatory and neuropathic pain conditions and can supplement the reduced endocannabinoid tone directly. For those with AA genotype and a personal or family history of autoimmune disease, periodic monitoring of inflammatory markers (CRP, ESR) provides early warning before symptoms become established.

Interactions

The most relevant interaction is with the functional coding variant rs35761398 (Q63R) in the same CNR2 gene. rs35761398 changes glutamine to arginine at position 63, altering receptor signalling efficacy. Carrying both rs2501401 AA (lower expression) and rs35761398 RR (altered receptor function) could theoretically compound impaired CB2 pathway activity, though no published study has examined this specific combination. rs35761398 is the variant with direct clinical association data and should be interpreted alongside rs2501401 when available.

NEGR1 (Neuronal Growth Regulator 1) encodes a GPI-anchored cell adhesion molecule¹¹ GPI-anchored cell adhesion molecule
Glycosylphosphatidylinositol — a lipid anchor that attaches proteins to the outer surface of cell membranes belonging to the IgLON family. It is strongly expressed in the hypothalamus — the brain region that governs hunger, satiety, and energy balance. The rs2815752 variant sits upstream of NEGR1 and alters its expression, making it one of the first obesity-associated genes identified by genome-wide association studies that pointed to a neuronal, rather than metabolic, mechanism for weight gain.

The Mechanism

NEGR1 promotes neurite outgrowth, synapse formation, and cell-cell adhesion in hypothalamic neurons, particularly in the arcuate nucleus²² arcuate nucleus
The arcuate nucleus contains hunger-promoting (AgRP/NPY) and satiety-promoting (POMC) neurons that integrate signals from hormones like leptin and ghrelin and ventromedial hypothalamus. It accumulates at GABAergic inhibitory synapses and promotes clustering of GAD65³³ GAD65
Glutamic acid decarboxylase 65 — the enzyme that synthesizes GABA, the main inhibitory neurotransmitter at synaptic membranes. When NEGR1 expression is reduced — as occurs with the A risk allele — hypothalamic appetite circuits are disrupted: the GABAergic synaptic architecture weakens, and the brain's ability to regulate food intake diminishes.

A linked 45-kb deletion polymorphism upstream of NEGR1 removes several conserved regulatory elements, likely explaining how the common risk haplotype reduces NEGR1 expression. ⁴⁴ Willer et al. Six new loci associated with body mass index highlight a neuronal influence on body weight regulation. Nat Genet, 2009

The Evidence

The NEGR1 locus was identified in the landmark GIANT consortium GWAS⁵⁵ GIANT consortium GWAS
Willer et al. Meta-analysis of 15 GWAS for BMI (n>32,000) with replication in 14 cohorts (n>59,000). Nat Genet, 2009 as one of six new BMI loci, with the A allele increasing BMI by 0.10 kg/m² per copy (P = 6.0 x 10^-8, combined n = 83,499). The largest BMI meta-analysis to date (339,224 individuals⁶⁶ 339,224 individuals
Locke et al. Genetic studies of body mass index yield new insights for obesity biology. Nature, 2015) confirmed NEGR1 among 97 genome-wide significant loci and highlighted its central nervous system expression pattern.

Mouse studies validate the human genetics. NEGR1-knockout mice show reduced body mass (8-13%)⁷⁷ reduced body mass (8-13%)
Lee et al. Functional inactivation of the GWAS obesity gene NEGR1 in mice causes a body mass phenotype. PLOS ONE, 2012, decreased food intake, and reduced lean mass — confirming NEGR1's role in appetite and body composition. In rats, decreased hypothalamic NEGR1 expression⁸⁸ decreased hypothalamic NEGR1 expression
Boender et al. The obesity-associated gene Negr1 regulates aspects of energy balance in rat hypothalamic areas. Physiol Genomics, 2014 directly increases body weight and food intake, particularly affecting carbohydrate preference.

In a human dietary study of 26,107 participants, the A risk allele was associated with lower fat intake but higher carbohydrate and fiber intake⁹⁹ A risk allele was associated with lower fat intake but higher carbohydrate and fiber intake
Rukh et al. Genetic susceptibility to obesity and diet intakes in the Malmö Diet and Cancer Study. Genes Nutr, 2013, suggesting NEGR1 influences not just how much we eat but what we prefer to eat.

Practical Actions

Because the A allele is extremely common (carried by 86% of Europeans and over 99% of East Asians as at least one copy), this is not a rare high-impact variant — it is a common, moderate-effect contributor to obesity susceptibility. The practical value lies in understanding your position on the appetite regulation spectrum.

Carriers of two A alleles should focus on structured meal timing and protein-forward meals to support satiety signaling that their hypothalamic wiring handles less efficiently. Monitoring waist circumference rather than just body weight captures the metabolic risk more accurately.

Interactions

NEGR1 contributes to a polygenic obesity risk profile alongside FTO (rs9939609) and MC4R (rs17782313). While each variant has a modest individual effect, carrying risk alleles at multiple loci has a cumulative impact on BMI. A combined FTO + NEGR1 risk genotype may warrant more structured appetite management strategies than either alone.

NEGR1 is also a cross-trait hit for major depressive disorder, with integrative genomic analysis¹⁰¹⁰ integrative genomic analysis
Li et al. Integrating GWAS and eQTL data identifies NEGR1 as a causal risk gene of major depression. J Affect Disord, 2020 identifying it as a causal MDD risk gene. The shared neurodevelopmental mechanism — disrupted synaptic connectivity in hypothalamic and limbic circuits — may explain why obesity and depression so often co-occur. Carriers of the A allele at rs2815752 who also carry depression-associated variants may experience compounded effects on both appetite regulation and mood.

Deep inside every primordial follicle, an arrested oocyte faces a biological paradox: it must preserve its genetic integrity across potentially decades before fertilisation, while its genome is under constant threat from parasitic DNA elements called transposable elements¹¹ transposable elements
mobile genetic sequences that can copy themselves and insert elsewhere in the genome, causing mutations if unchecked. The piRNA pathway — named for the PIWI-interacting RNAs that guide it — is the germline's primary defence against this threat. PIWIL1 is one of four human PIWI proteins, and it is expressed in meiotic oocytes in primordial follicles, where it helps silence endogenous retroviruses and other transposons that would otherwise compromise egg quality and genome stability.

The rs28416520 variant sits in a CpG-rich regulatory region of the PIWIL1 promoter. The A allele modifies local CpG methylation patterns, altering PIWIL1 transcription in germ cells. Critically, its effect on ovarian ageing follows a recessive inheritance pattern²² recessive inheritance pattern
recessive: two copies of the risk allele are required to produce the full effect; one copy alone has little impact.

The Mechanism

PIWIL1 belongs to the Argonaute superfamily of RNA-guided silencing proteins. In oocytes, it binds piRNAs — a class of small non-coding RNAs (23–29 nucleotides) derived from transposon sequences — and uses these guides to recognise and cleave complementary transposon transcripts. In golden hamsters, which express PIWIL1 throughout oogenesis analogously to humans, disruption of PIWIL1 causes oocyte dysfunction through transposon de-silencing and widespread transcriptomic dysregulation, leading to embryonic arrest³³ disruption of PIWIL1 causes oocyte dysfunction through transposon de-silencing and widespread transcriptomic dysregulation, leading to embryonic arrest
Lim et al. 2021, Nature Cell Biology — PIWIL1-null female hamsters are sterile due to piRNA-pathway failure in oocytes. In humans, PIWIL1 nuclear foci are observed specifically in meiotic oocytes within primordial follicles.

The rs28416520 promoter variant likely reduces PIWIL1 expression in a dosage-dependent manner: one A allele can be compensated by the remaining G allele, but AA homozygosity results in insufficient PIWIL1 activity to maintain full transposon surveillance. Over the long dormant period of primordial follicle maintenance, progressive transposon-mediated DNA damage may accelerate follicle loss and bring forward the timing of menopause.

The Evidence

The primary evidence comes from a landmark genome-wide association study by Ruth et al. 2021⁴⁴ genome-wide association study by Ruth et al. 2021
Genetic insights into biological mechanisms governing human ovarian ageing. Nature 596:393–397 examining approximately 200,000 women of European ancestry. This study identified 290 genetic determinants of ovarian ageing, measured as variation in age at natural menopause (ANM). rs28416520 at the PIWIL1 locus was among the variants that showed significant departure from the additive allelic model, exhibiting instead a recessive pattern — meaning AA homozygotes showed a substantially larger effect on menopause timing than would be predicted if the effects simply added up per allele. The study also demonstrated that experimental manipulation of piRNA-pathway genes highlighted by this GWAS increases fertility and extends reproductive life in mice, directly validating the biological relevance of this pathway to ovarian reserve.

A promoter functional study of rs28416520⁵⁵ promoter functional study of rs28416520
Zhang et al. 2020, assessing CpG-region SNPs in the PIWIL1 promoter in a Chinese gastric cancer cohort confirms that this SNP alters methylation-sensitive regulatory activity in the PIWIL1 promoter, establishing that it has measurable effects on PIWIL1 transcription.

A comprehensive study of inherited piRNA pathway defects⁶⁶ comprehensive study of inherited piRNA pathway defects
Wyrwoll et al. 2024, Nature Communications — 39 infertile men carrying biallelic variants in 14 piRNA pathway genes including PIWIL1 demonstrates that PIWIL1 loss-of-function follows an autosomal recessive pattern in humans and causes germline transposon de-repression with consequent germ cell failure. This human genetic data is entirely consistent with the recessive model seen for rs28416520 in the ovarian ageing GWAS.

Practical Actions

For the roughly 14% of women who carry two A alleles, the implications are most relevant to reproductive planning. The recessive architecture means GA heterozygotes (approximately 46% of women) are not meaningfully affected — only AA homozygotes carry the elevated risk. AA women may experience earlier decline in ovarian reserve than the population average, warranting earlier baseline measurement of anti-Müllerian hormone (AMH) and antral follicle count if family planning is a consideration.

Because the mechanism involves piRNA-pathway surveillance of transposon activity in oocytes, and because oxidative DNA damage is known to accelerate follicle depletion generally, strategies to protect oocyte genome integrity — particularly avoiding genotoxic exposures such as tobacco smoke — are especially relevant for AA homozygotes.

Interactions

PIWIL1 functions within a larger piRNA biogenesis network that includes PIWIL2, PIWIL3, and PIWIL4, plus scaffold proteins (TDRD1, TDRD9, HENMT1, MAEL) and RNA-modifying enzymes. Women carrying variants in multiple piRNA pathway genes may have compounded effects on piRNA biogenesis. The ovarian ageing GWAS (Ruth et al. 2021) also implicated other DNA damage response loci, including rs10183486 (TLK1) and rs16991615 (MCM8), that may act through convergent pathways with PIWIL1 to influence the cumulative rate of follicle depletion. A formal compound analysis of PIWIL1 + TLK1 genotypes has not been published but is biologically plausible given both genes act to protect oocyte DNA integrity.