ATG16L1 rs2241879 — Autophagy, Paneth Cells, and Crohn's Disease Risk

ATG16L1 (Autophagy Related 16 Like 1)11 ATG16L1 (Autophagy Related 16 Like 1)
A core scaffold protein required for autophagosome formation — the membrane-bound compartment that engulfs and digests intracellular bacteria, damaged organelles, and protein aggregates
is one of the strongest and most replicated genetic risk loci for Crohn's disease22 Crohn's disease
A form of inflammatory bowel disease (IBD) that causes transmural intestinal inflammation, most commonly affecting the ileum and colon, driven by a dysregulated immune response to gut bacteria
. rs2241879 is an intronic variant within ATG16L1 on chromosome 2q37.1 that sits in tight linkage disequilibrium33 linkage disequilibrium
Two variants are in LD when they are inherited together far more often than expected by chance — they effectively act as proxies for each other
with the coding variant rs2241880 (T300A). It tags the same disease-associated haplotype as T300A and reflects identical IBD risk in the European populations where both were originally characterised.

The Mechanism

ATG16L1 is indispensable for the formation of autophagosomes — the double-membrane vesicles that capture intracellular cargo and deliver it to lysosomes for degradation. In the gut, this xenophagy44 xenophagy
Selective autophagy specifically targeting intracellular pathogens rather than the cell's own components
pathway is the primary mechanism by which intestinal epithelial cells clear bacteria such as Salmonella and Yersinia that breach the mucosal barrier. ATG16L1 cooperates with NOD255 NOD2
Nucleotide-binding oligomerisation domain 2, a bacterial pattern-recognition receptor that detects muramyl dipeptide (MDP) from bacterial cell walls and is itself a major Crohn's disease risk gene
to initiate bacterial autophagy at the site of pathogen entry.

The critical functional consequence of the T300A haplotype (tagged by rs2241879) was revealed by Murthy et al. 201466 Murthy et al. 2014
A landmark Nature paper demonstrating that the threonine-to-alanine substitution at position 300 creates a hypersensitive caspase-3 cleavage motif within ATG16L1
. Under conditions of metabolic stress, apoptotic signalling, or pathogen-induced cell death — all common events in an inflamed intestinal epithelium — caspase-3 becomes activated and cleaves ATG16L1. The T300A variant dramatically accelerates this cleavage, causing the protein to be destroyed faster than it can be replenished. The downstream consequence is an abrupt collapse of autophagy capacity precisely when the gut most needs to clear bacteria.

Paneth cells77 Paneth cells
Highly specialised secretory cells at the base of small intestinal crypts that release lysozyme, defensins, and other antimicrobial peptides into the intestinal lumen to maintain sterility of the stem cell niche
are disproportionately affected. ATG16L1 is selectively critical for Paneth cell biology: loss-of-function in mice produces abnormal cytoplasmic granule morphology, reduced lysozyme secretion, and dysregulated cytokine production (including elevated leptin and IL-1β). Human Crohn's patients homozygous for the ATG16L1 risk allele display the same Paneth cell abnormalities, and this cellular phenotype is associated with dysbiosis of the ileal microbiome — the site where Crohn's most commonly begins.

The Evidence

The original German association study by Glas et al. 200888 Glas et al. 2008
768 CD patients, 507 UC patients, 1,615 healthy controls, 9 ATG16L1 variants genotyped; rs2241879 showed the joint-strongest CD association alongside rs2241880
found rs2241879 to carry the same effect size as the T300A coding variant (OR 0.74, 95% CI 0.65–0.84, p=3.6×10⁻⁶) — consistent with the two variants being in tight LD on the same haplotype. The A allele (non-risk) was protective across all nine ATG16L1 variants tested.

Replication came rapidly. Prescott et al. 200799 Prescott et al. 2007
Independent UK cohorts of 1,236 CD cases and 1,235 controls
confirmed the T300A association with a 1.65-fold overall Crohn's risk and a 2.2-fold risk for ileal disease — the anatomical location where Paneth cells are most abundant. The ileal specificity is mechanistically coherent: Paneth cells exist only in the small intestine and are uniquely dependent on ATG16L1 for granule biology. The landmark Zhang et al. 2009 meta-analysis1010 Zhang et al. 2009 meta-analysis
24 studies, 13,022 CD cases, 17,532 controls
confirmed OR 1.87 (95% CI 1.69–2.05) for GG vs AA genotypes in Caucasian populations, with no significant effect in Asian populations, consistent with the different allele frequency distributions across ancestries. The functional mechanism was completed by Murthy et al. 20141111 Murthy et al. 2014
Nature; knock-in mice with T300A showed defective Yersinia clearance and elevated inflammatory cytokines, rescued by caspase-3 deletion
.

Evidence level is strong: multiple large European cohorts, a definitive meta-analysis, a clear molecular mechanism, and a causal animal model. The absence of association in Asian populations reflects the very different allele frequency spectrum (G allele ~72% in East Asians vs ~49% in Europeans), not a lack of biological plausibility.

Practical Actions

The rs2241879 G allele raises Crohn's disease risk through impaired bacterial clearance and Paneth cell dysfunction. Recommended actions focus on gut microbiome support, early recognition of IBD symptoms, and specific nutritional approaches that can compensate for impaired xenophagy.

For individuals carrying one or two G alleles, the most important modifiable factor is the gut microbiome: a diverse microbiome with high Lactobacillus and Bifidobacterium content reduces the challenge load that ATG16L1-dependent xenophagy must handle. Conversely, broad-spectrum antibiotic use and ultra-processed food diets that deplete microbial diversity increase the risk that uncleared bacteria trigger mucosal inflammation. Fermented foods that deliver live microorganisms directly to the ileum are particularly relevant given the ileal tropism of Paneth cell defects in this genotype.

Specific nutritional research on ATG16L1 variants suggests that reducing dietary patterns that promote intestinal permeability — particularly high intake of emulsifiers (carboxymethylcellulose, polysorbate 80) and refined sugars — may reduce luminal bacterial translocation and the demand for autophagy-mediated clearance. Vitamin D has a documented role in Paneth cell antimicrobial peptide production and ATG16L1 autophagy regulation, making adequate vitamin D status (serum 25(OH)D ≥50 nmol/L) especially relevant for carriers of this variant.

Interactions

rs2241879 (ATG16L1) and rs4958847 (IRGM) are partners in the same xenophagy pathway. IRGM encodes an immunity-related GTPase that acts upstream of ATG16L1, facilitating the selective recruitment of autophagy machinery to bacteria. Individuals carrying risk alleles at both loci face a compounded impairment: IRGM variants reduce autophagy initiation, and the ATG16L1 T300A haplotype accelerates degradation of the resulting ATG16L1 protein under stress. Multiple GWAS of IBD and CD have confirmed that risk allele burden across IRGM, ATG16L1, and NOD2 loci produces substantially greater CD risk than any single locus alone. The interaction is described in prose for compound action processing by the supervisor.

SNX19 — The Endolysosomal Positioning Factor Linked to Coronary Risk

Deep inside every cell, a network of membrane-bound organelles sorts, degrades, and recycles proteins and lipids11 sorts, degrades, and recycles proteins and lipids
Endolysosomes are the cell's recycling centers; their positioning determines how efficiently cargo is processed
. SNX19 — Sorting Nexin 19 — is a molecular tether that anchors endolysosomes to the endoplasmic reticulum (ER), keeping them clustered near the cell nucleus rather than scattered throughout the cytoplasm. A common missense variant in SNX19, rs2298566, alters one of the protein's functional domains and has been associated with elevated coronary heart disease risk in prospective cohort data.

The Mechanism

SNX19 bridges two organelle systems: its two N-terminal transmembrane domains embed the protein in the ER membrane, while its PX domain22 PX domain
Phox Homology domain — a phosphoinositide-binding module found across the sorting nexin family
binds phosphatidylinositol 3-phosphate (PI3P) on endolysosomal membranes. Regulatory PXA and PXC flanking domains act as a molecular governor, preventing excessive tethering under normal conditions. The variant rs2298566 is a missense change that substitutes leucine at position 878 with arginine (p.Leu878Arg) in the major isoform — a significant amino acid property change within the regulatory C-terminal region.

When endolysosomal positioning is disrupted — as demonstrated by Saric et al., Nature Communications 202133 Saric et al., Nature Communications 2021
NIH National Institute of Child Health and Human Development
using SNX19-knockout cells — endolysosomes disperse throughout the cytoplasm and exhibit increased motility. Perinuclear clustering of endolysosomes is important for efficient lysosomal degradation of oxidized LDL, inflammatory signaling complexes, and autophagy substrates. Disruption of this compartment organization is implicated in the accumulation of intracellular lipid and inflammatory cargo relevant to vascular cell biology.

The precise molecular bridge between altered SNX19 tethering and coronary artery disease risk is not yet established at the mechanistic level. SNX19's role in endolysosomal positioning is consistent with the broader evidence that lysosomal dysfunction contributes to macrophage foam cell formation, cholesterol efflux impairment, and vascular inflammation — central processes in atherosclerosis.

The Evidence

The primary evidence for rs2298566's cardiovascular association comes from the Atherosclerosis Risk in Communities (ARIC) study. Bare et al., Genetics in Medicine 200744 Bare et al., Genetics in Medicine 2007
Cohort study with Cox proportional hazards modeling, 9,129 white participants, median 13-year follow-up
assembled five variants — rs20455 (KIF6), rs3900940 (MYH15), rs7439293 (PALLD), rs2298566 (SNX19), and rs1010 (VAMP8) — each of which had prior CHD associations. Among participants with a high composite genetic risk score (the top 4% of carriers), the hazard ratio for incident CHD was 1.57 (95% CI 1.21–2.04; p = 0.001) after adjustment for traditional cardiovascular risk factors; bootstrap validation suggested HR 1.43 in external populations. The paper states rs2298566 had been "confirmed in two studies" prior to ARIC inclusion.

It is important to note the limitations. The published effect size is for the composite five-SNP score — individual effect sizes for rs2298566 alone are not reported in the abstract, and its specific contribution to the composite has not been disaggregated in the literature retrieved. The KIF6 variant (rs20455), which was the most prominent member of this panel, later failed to replicate across large independent cohorts. Whether rs2298566 independently replicates at genome-wide significance has not been established in a GWAS meta-analysis. The evidence level is accordingly rated moderate rather than strong.

Practical Actions

For A-allele carriers, the most actionable steps target the lipid-processing and inflammatory pathways that SNX19's endolysosomal role touches. Monitoring the early biomarkers of cardiovascular risk — lipids and vascular inflammation — allows timely intervention at the stage before clinical coronary disease develops.

There are no specific drug-gene interactions or nutrient interventions documented for rs2298566 in the current literature. Statin response differential has been proposed for genes in this five-SNP panel but has not been established for SNX19 specifically.

Interactions

rs2298566 was studied as one of five CHD-associated variants alongside rs20455 (KIF6), rs3900940 (MYH15), rs7439293 (PALLD), and rs1010 (VAMP8). The composite risk score framework suggests additive rather than synergistic effects across these loci, as each gene operates through largely distinct mechanisms (cytoskeletal dynamics, actin scaffolding, vesicle trafficking). The cardiovascular relevance of SNX19 may also intersect with the 11q25 psychiatric risk locus; SNX19 transcript diversity at this locus is regulated by chromatin accessibility changes relevant to the same cis-regulatory region.

rs2371365

PCLO PCLO Presynaptic Scaffold

Moderate Risk Factor

Piccolo: The Synapse's Master Organizer and Mood Risk

Every time a neuron signals its neighbor, a precisely choreographed event unfolds at the presynaptic terminal: calcium floods in, synaptic vesicles fuse with the membrane, and neurotransmitters flood the synapse. Orchestrating this entire process is a giant scaffolding protein called Piccolo11 Piccolo
encoded by the PCLO gene, located at chromosome 7q11.23
. Piccolo anchors calcium channels, tethers vesicles, and coordinates the actin cytoskeleton that drives vesicle recycling. Without it functioning optimally, the timing and magnitude of neurotransmitter release — including serotonin, dopamine, and norepinephrine — can go subtly awry.

rs2371365 is an intronic variant in PCLO that sits approximately 53 kb from rs2522833, the functional missense variant (Ser4814Ala) that was the top association signal in the first genome-wide study of major depressive disorder to implicate this gene. As an intronic SNP in moderate linkage disequilibrium with the coding variant region, rs2371365 is a proxy marker for the broader PCLO risk haplotype rather than a direct functional mutation itself.

The Mechanism

The PCLO Ser4814Ala substitution (rs2522833), with which rs2371365 is in partial LD22 partial LD
linkage disequilibrium; both variants tag the same 167 kb PCLO risk haplotype identified by Sullivan et al.
, resides near a C2 calcium-binding domain of the Piccolo protein. C2 domains are calcium sensors that trigger phospholipid binding and vesicle fusion events. Functional studies in cultured neurons carrying the Ser4814Ala substitution showed 30% increased excitatory synaptic transmission and elevated Piccolo protein levels at synapses33 30% increased excitatory synaptic transmission and elevated Piccolo protein levels at synapses
Giniatullina et al. 2015; these changes suggest compensatory upregulation in response to altered C2 calcium sensing
.

At the cellular level, Piccolo also regulates presynaptic F-actin assembly by scaffolding actin regulatory proteins including Daam1. When Piccolo is disrupted, boutons show enhanced activity-dependent vesicle exocytosis and reduced F-actin polymerization44 boutons show enhanced activity-dependent vesicle exocytosis and reduced F-actin polymerization
Waites et al. 2011, J Neurosci; loss of Piccolo function paradoxically increases short-term exocytosis while impairing sustainable vesicle recycling
, which may dysregulate monoamine release dynamics across serotonergic, dopaminergic, and noradrenergic synapses in the brainstem and limbic system.

The Evidence

The initial genome-wide association study by Sullivan et al. 2009 in Molecular Psychiatry55 Sullivan et al. 2009 in Molecular Psychiatry
1,738 MDD cases vs. 1,802 controls; identified 11 genome-wide signals in a 167 kb PCLO region; top hit rs2715148 p=7.7×10⁻⁷, rs2522833 p=1.2×10⁻⁶
was the first to implicate PCLO in major depressive disorder. A population-based replication by Hek et al. 201066 Hek et al. 2010
579 depression cases, 912 controls; confirmed rs2522833 association, p=0.0025; meta-analysis across three population-based studies reached p=1.93×10⁻⁹
brought the evidence to genome-wide significance.

The most striking findings concern the brain's emotional processing circuitry. An fMRI study by Woudstra et al. 201277 fMRI study by Woudstra et al. 2012
22 MDD patients and 29 healthy controls; PCLO risk allele carriers showed significantly increased left amygdala activity during angry and sad face processing; effects on fearful faces were specific to MDD
demonstrated that PCLO risk allele carriers — regardless of depression status — show heightened amygdala reactivity. A follow-up study examining emotional memory88 emotional memory
Woudstra et al. 2013, PLoS One; N=89 MDD and 29 controls; risk carriers showed reduced striatal encoding of negative words and blunted amygdalar response to novel positive stimuli
found that risk carriers also show blunted reward signaling, consistent with the anhedonia that characterizes major depression.

Personality correlates were documented by Minelli et al. 201299 Minelli et al. 2012
522 MDD patients, 375 controls; CC homozygotes significantly overrepresented in depressed group, p<0.01; C-allele carriers in controls showed elevated Harm Avoidance and reduced Novelty Seeking on the Tridimensional Personality Questionnaire
, suggesting the risk variant shapes temperamental traits that predispose toward depression before a clinical episode occurs.

At the neuroendocrine level, Schuhmacher et al. 20111010 Schuhmacher et al. 2011
205 depressed inpatients followed through 4-week antidepressant treatment; C-allele carriers showed greater initial HPA axis reactivity and stronger hormonal responsiveness during treatment
linked PCLO genotype to stress hormone regulation, a core feature of melancholic depression.

Practical Actions

Carriers of one or two C alleles at rs2371365 appear to have a subtly altered presynaptic environment, particularly in monoaminergic circuits governing emotional processing. This does not predict depression deterministically — the variant contributes a modest shift in baseline emotional reactivity and stress system tone. Practical steps focus on protecting monoamine system resilience, stabilizing presynaptic calcium signaling, and supporting HPA axis regulation.

The heightened amygdala reactivity documented in risk carriers makes cognitive-reappraisal and mindfulness-based interventions particularly relevant as they act specifically on amygdala-prefrontal circuitry. Monitoring via validated mood self-tracking (PHQ-9 or GAD-7) gives early warning before subclinical symptoms progress. If antidepressants are ever needed, the HPA axis modulation data suggest PCLO genotype may eventually inform treatment selection.

Interactions

The strongest documented interaction involves rs2522833 in PCLO itself: rs2371365 is a proxy marker for the same risk haplotype, so both variants in the same person are likely tagging the same underlying biological signal rather than independent effects. If rs2522833 data is available, that missense variant has more direct functional evidence.

PCLO risk has also been examined in the context of serotonin transporter (5-HTTLPR, rs25531) and BDNF Val66Met (rs6265) variation; epistatic analyses suggest additive or synergistic effects on limbic reactivity when multiple monoamine pathway variants co-occur, though published compound analyses remain limited.

PTPN22 -1123G>C — A Second Voice in the Same Gene

The PTPN22 gene is best known for its coding variant rs2476601 (R620W), the strongest non-HLA risk factor for autoimmunity in Europeans. But the gene harbors a second clinically relevant variant in its promoter region: -1123G>C (rs2488457), located 1,123 base pairs upstream of the transcription start site. Unlike R620W, which changes the protein's structure, this promoter variant sits in regulatory territory — influencing not what PTPN22 does, but [how much of it is made | Promoter variants affect transcription factor binding and thus the rate of mRNA production from the gene]. Its significance depends heavily on ancestry: in Asian populations, where R620W is virtually absent, -1123G>C stands alone as a primary PTPN22 risk signal11 stands alone as a primary PTPN22 risk signal
In Chinese and Japanese cohorts, rs2476601 is non-polymorphic while rs2488457 shows independent association with RA, T1D, LADA, JIA, and UC
. In Europeans, it rides along with the R620W haplotype and cannot be disentangled from it22 cannot be disentangled from it
Norwegian cohort study found -1123G>C association statistically indistinguishable from 1858C>T in RA
in most cohorts.

The Mechanism

The -1123G>C change alters a nucleotide in the 5' regulatory region of PTPN22, a zone that binds transcription factors controlling how actively the gene is transcribed. Expression studies have found a tendency for association between -1123G>C and PTPN22 mRNA levels33 tendency for association between -1123G>C and PTPN22 mRNA levels
Chinese RA study found expression analysis indicated association between -1123G>C and PTPN22 gene expression levels
, though the direction and magnitude have not been pinned down with the same precision as functional coding variants. A Chinese ulcerative colitis study found dramatically elevated PTPN22 mRNA in inflamed colonic tissue44 dramatically elevated PTPN22 mRNA in inflamed colonic tissue
PTPN22 mRNA significantly higher in inflamed vs non-inflamed colon; correlated with CRP r=0.578, p<0.001
, with expression levels correlating with disease severity and systemic inflammation markers. Whether the promoter variant drives this upregulation, or is a passenger tagging another functional element, remains an open question.

A graft-versus-host disease study provided an unexpected functional clue: bone marrow transplant recipients carrying the CC genotype had lower incidence of acute GvHD55 CC genotype had lower incidence of acute GvHD
C/C recipients HR 0.50 for grade II-IV acute GvHD compared to G/G carriers, but higher relapse rate HR 1.78
but higher relapse rates — the classic immune-suppression trade-off. This suggests the C allele (the common protective variant) may dampen immune responsiveness, while the G allele sustains a more active immune phenotype that fights both GvHD targets and autoantigens.

The Evidence

The variant's association spectrum spans multiple autoimmune conditions, but with marked population stratification66 marked population stratification
Variant is polymorphic and shows independent autoimmune associations in Asian populations; in Europeans, tight LD with R620W makes independent contribution undetectable
:

In Asian populations where R620W is absent, -1123G>C shows independent associations. A Japanese/Korean study found OR=1.41-1.42 for acute-onset type 1 diabetes77 found OR=1.41-1.42 for acute-onset type 1 diabetes
Promoter SNP associated with acute-onset but not slow-onset T1D; OR=1.42 in Japanese, combined OR=1.41
, with the authors arguing the promoter SNP is "a more likely causative variant" than R620W in these populations. A Chinese RA study found OR=1.52 for RA in Han Chinese88 OR=1.52 for RA in Han Chinese
494 cases and 496 controls; R620W was non-polymorphic in this cohort
. A Chinese LADA study reported OR=1.99 for latent autoimmune diabetes99 OR=1.99 for latent autoimmune diabetes
Significant association in Chinese Hans; R620W showed no association in same cohort
. In juvenile idiopathic arthritis, a Han Chinese study found OR=2.15 for GC/CC combined vs GG1010 OR=2.15 for GC/CC combined vs GG
137 JIA cases vs 150 controls; C allele OR=1.95
.

A meta-analysis of uveitis1111 meta-analysis of uveitis
8 studies; OR=1.18 overall, OR=1.21 in Asian populations; weaker signal in Europeans
found modest but significant association. For primary immune thrombocytopenia1212 primary immune thrombocytopenia
Meta-analysis of 10 studies, 932 cases; G allele carriers OR=1.23; GG homozygotes 1.51× more susceptible than CC carriers
, a meta-analysis found G allele carriers at modestly elevated risk.

In European populations, a Czech/Azeri study concluded no independent contribution1313 no independent contribution
Only R620W haplotype drove disease risk; -1123G>C alone conferred no additional risk
from -1123G>C in T1D or JIA. The Norwegian RA study could not distinguish1414 could not distinguish
Tight LD between the two variants makes conditional analysis impossible in European samples
between the two variants' contributions. In Europeans, this variant's clinical relevance is primarily as a tag for the R620W haplotype.

Practical Implications

The G allele (minor allele: ~22% in Europeans, ~40% in East Asians) is the risk variant associated with altered PTPN22 expression and modestly elevated autoimmune susceptibility. Risk associations are modest (OR 1.2–2.0 depending on condition and ancestry) and context-dependent. For individuals of East Asian ancestry, this variant has independent predictive value where R620W is not relevant. For individuals of European ancestry, carrying the G allele is primarily informative as a haplotype marker — its combined effect with R620W (rs2476601) on the same haplotype amplifies the overall PTPN22 risk signal, and joint genotyping provides fuller picture than either variant alone.

Monitoring guidance follows the same logic as for other PTPN22 variants: awareness of early autoimmune symptoms (joint inflammation, thyroid changes, skin changes, fatigue) and prompt evaluation if they develop. The variant does not point to a specific actionable intervention but raises the baseline vigilance warranted for immune dysregulation.

Interactions

rs2488457 and rs2476601 (R620W) are the two best-characterized PTPN22 variants, and their relationship depends on ancestry1515 their relationship depends on ancestry
In Europeans, both are on the same high-risk haplotype and are in tight LD; in Asian populations, rs2476601 is absent and rs2488457 acts independently
. Czech and Azeri haplotype analysis showed only haplotypes carrying the R620W minor allele conferred disease risk — adding -1123G>C to a non-R620W haplotype added nothing. This makes the two variants co-travelers in Europeans but independent signals in Asians.

A third PTPN22 variant, rs33996649 (+788G>A), is also non-polymorphic in most Asian populations and tracks with R620W in Europeans, adding to the picture of a PTPN22 "risk haplotype" in Europeans that encodes multiple regulatory and coding changes simultaneously.

Compound interactions between rs2488457 and HLA loci have not been formally characterized but are expected to follow the same pattern as R620W: convergent independent risk from PTPN22 and HLA without strong epistasis.

MEFV M694I — A Founder FMF Mutation with Moderate Severity and High Colchicine Response

Familial Mediterranean fever11 Familial Mediterranean fever
FMF; an autoinflammatory disease caused by dysregulation of the pyrin inflammasome, characterized by recurrent self-limiting attacks of fever, peritonitis, pleuritis, and arthritis
is the most common hereditary periodic fever in the world. The MEFV gene encodes pyrin, a protein expressed in myeloid cells that normally suppresses inflammasome activation in the absence of microbial toxins. Of the hundreds of documented MEFV variants, five founder mutations — M694V, V726A, M680I, M694I, and E148Q — account for roughly 74% of FMF chromosomes in classic cases. M694I (p.Met694Ile) sits at codon 694 in exon 10, just two nucleotides from the more severe M694V mutation, yet the isoleucine substitution confers a distinctly milder phenotype. While M694V dominates the European and Turkish mutation landscape, M694I is the signature exon 10 mutation in East Asian FMF patients and a founder allele in Lebanese kindreds.

The Mechanism

Exon 10 of MEFV encodes part of pyrin's B30.2/SPRY domain, the region that senses microbial toxins and interacts with the regulatory kinases PKN1/2 that keep the inflammasome dormant. Pathogenic exon 10 missense variants impair this gate, allowing unrestrained caspase-122 caspase-1
the protease that cleaves pro-IL-1β and pro-IL-18 into their active inflammatory forms
activation and cytokine release. M694I replaces methionine with isoleucine — a relatively conservative change compared to M694V's valine substitution — but still destabilizes the B30.2 regulatory interface sufficiently to dysregulate pyrin gating.

A 2025 CRISPR/Cas9 knock-in mouse study (Koga et al.33 Koga et al.) showed that M694I knock-in homozygous mice develop a selective enhancement of Th17 cell differentiation alongside elevated serum G-CSF, IFN-γ, IL-1α, IL-5, IL-6, and TNF-α (all P < 0.05), with significantly reduced survival (P < 0.001). This identifies a Th17-driven inflammatory axis as a specific downstream effect of M694I — distinct from the predominantly IL-1β/monocyte-driven mechanism of M694V — and points toward potential IL-17-pathway therapeutic targets for refractory cases.

The Evidence

The clinical picture of M694I emerges most clearly from Arab and East Asian cohort studies, where it appears at clinically significant frequencies. Majeed et al. (2002)44 Majeed et al. (2002) genotyped 278 Arab FMF patients and found that the M694I/M694I genotype carried a mean severity score of 6 ± 1 — compared to 14 ± 2 for M694V/M694V and 10 ± 3 for V726A/V726A (P = 0.003). In this Arab cohort, M694I/M694I represented 14% of identifiable biallelic genotypes. The severity differential is strikingly wide: M694I homozygotes experienced disease roughly half as severe as exon-10-mutation compound heterozygotes.

In Japan, where M694V is essentially absent, M694I is the primary exon 10 mutation. Tsuchiya-Suzuki et al. (2009)55 Tsuchiya-Suzuki et al. (2009) found M694I in 42.5% of 80 Japanese FMF patients, predominantly as E148Q/M694I compound heterozygotes (25%) or heterozygous M694I carriers (17.5%). The national registry study by Migita et al. (2012)66 Migita et al. (2012) confirmed E148Q/M694I (19.8%) and M694I alone (12.7%) as the dominant mutation patterns in 126 Japanese patients, with colchicine effective in 91.8% at a low mean dose of 0.89 mg/day. In a 116-patient Japanese cohort, Kishida et al. (2014)77 Kishida et al. (2014) found M694I carriers had a more severe clinical course than E148Q carriers but a very favorable colchicine response, suggesting the mutation is clinically meaningful but pharmacologically manageable.

Risk of AA amyloidosis exists but is lower than for M694V. Nakamura et al. (2014)88 Nakamura et al. (2014) reported a 51-year-old Japanese male with M694I/M694I who developed biopsy-confirmed renal AA amyloidosis; colchicine resolved both the inflammatory attacks and kidney dysfunction. This case establishes that M694I homozygosity can — given sufficient inflammatory burden and decades without suppressive therapy — progress to the same amyloid complication that defines the natural history of severe FMF.

Practical Actions

For heterozygous carriers, clinical FMF from a single M694I allele is rare; disease typically requires biallelic mutations or compound heterozygosity with a second pathogenic MEFV allele. However, heterozygous carriers should receive genetic counseling because their offspring have a 50% chance of inheriting the allele, and if both partners carry MEFV mutations, the risk of homozygous or compound heterozygous children is significant.

For homozygous or compound heterozygous carriers: colchicine at 0.5–1 mg/day is the standard and highly effective first-line treatment. Japanese registry data show greater than 90% colchicine response rates in M694I-positive patients at lower doses than required for M694V. Monitoring serum amyloid A (SAA) and CRP between attacks provides the earliest signal of sustained subclinical inflammation, the primary driver of long-term amyloidosis risk. Annual urine protein screening is recommended by EULAR FMF guidelines for any biallelic MEFV carrier to detect early renal AA amyloidosis before kidney function declines.

Interactions

The most clinically important interaction for M694I is compound heterozygosity with the exon 2 variant E148Q (rs3743930). This combination — E148Q/M694I — is the dominant mutation pattern among Japanese FMF patients. In a pediatric cohort study (Miyashita et al. 202299 Miyashita et al. 2022), compound E148Q/M694I children had typical FMF attacks with dramatically elevated IL-18 (2,806 ± 2,107 pg/mL during afebrile phases), while single-mutation carriers of either E148Q or M694I alone showed no fever or serositis — pointing to a true synergistic inflammasome interaction. Lebanese founder-effect families carrying M694I as part of a complex haplotype (Medlej-Hashim et al. 2011, PMID 20937419) show that M694I behaves as a fully penetrant disease allele when combined with another pathogenic MEFV mutation. M694I can also compound with M694V (rs61752717), M680I, and V726A (rs28940579), though these combinations are less commonly reported than E148Q/M694I.

ATP7B R778L — The East Asian Wilson Disease Variant

Copper is an essential mineral — a cofactor for enzymes involved in energy production, antioxidant defence, and neurotransmitter synthesis — yet it is toxic when it accumulates. The liver is the body's copper thermostat: it absorbs dietary copper from the portal circulation, uses what is needed, loads the rest onto the plasma protein ceruloplasmin11 ceruloplasmin
The main copper-carrying protein in blood; synthesised in the liver as apo-ceruloplasmin and activated when ATP7B loads copper onto it before secretion. Low serum ceruloplasmin is a hallmark diagnostic marker of Wilson disease. Approximately 95% of plasma copper is bound to ceruloplasmin.
, and exports any surplus into bile for elimination. ATP7B, encoded on chromosome 13, is the transmembrane P-type ATPase responsible for both of these functions. When both copies of ATP7B are non-functional, copper cannot leave the liver — and Wilson disease22 Wilson disease
A rare autosomal recessive disorder of copper metabolism first described by neurologist Samuel Alexander Kinnier Wilson in 1912. Estimated prevalence 1 in 30,000; carrier frequency approximately 1 in 90.
follows.

The R778L variant (p.Arg778Leu, c.2333G>T) replaces arginine with leucine at amino acid position 778 in the ATP-binding domain of ATP7B. Identified as ClinVar Pathogenic variant VCV000003852, it is the most common Wilson disease mutation across East Asia, accounting for 28.96% of all pathogenic ATP7B alleles33 accounting for 28.96% of all pathogenic ATP7B alleles
Zhang et al. 2022, Translational Neurodegeneration, cohort of 1,302 Chinese Wilson disease patients
in Chinese patients — and for approximately 17% of alleles in Hong Kong Chinese, where haplotype analysis reveals a single ancestral founder at least 5,500 years old44 where haplotype analysis reveals a single ancestral founder at least 5,500 years old
Mak et al. 2008, J Hum Genet
.

The Mechanism

The ATP7B gene lies on the minus strand of chromosome 13. At GRCh38 position 13:51,958,333, the plus-strand reference nucleotide is C (corresponding to coding-strand G = arginine at position 778). The pathogenic alternate allele is A on the plus strand (coding-strand T = leucine substitution). Both WGS and consumer genotyping chips report alleles on the plus strand, so the wild-type homozygote appears as CC and the risk homozygote as AA.

Arginine 778 sits within the ATP-binding domain, a region essential for the phosphorylation cycle that drives copper translocation across the membrane. The arginine-to-leucine substitution does not merely reduce enzymatic efficiency — it fundamentally mislocalises the protein55 mislocalises the protein
Zhu et al. 2015, Mol Cell Neurosci, showed R778L uniquely disrupts BOTH subcellular localization AND vesicular trafficking of ATP7B, while P992L only disrupts trafficking; the mutant protein fails to reach the trans-Golgi network where copper loading onto ceruloplasmin occurs
. The ATP7B protein cannot reach its functional destination. As a consequence, copper accumulates in hepatocytes, ceruloplasmin secretion falls, and biliary copper excretion is abolished.

In R778L homozygous iPSC-derived hepatocytes, ATP7B protein expression and ceruloplasmin secretion are both reduced compared to controls66 ATP7B protein expression and ceruloplasmin secretion are both reduced compared to controls
Song et al. 2022, Hum Mol Genet; partial genetic correction (heterozygous CRISPR repair) restores ceruloplasmin expression, suggesting even one functional allele is sufficient to maintain near-normal copper handling
. This haploinsufficiency threshold is why heterozygous carriers are clinically unaffected.

The Evidence

Clinical disease follows copper accumulation in a predictable sequence. Hepatic copper overload causes liver inflammation and fibrosis, sometimes presenting as acute hepatitis or cirrhosis in children and young adults. Copper eventually spills into the systemic circulation and deposits in the brain (causing movement disorders, cognitive changes, and psychiatric symptoms), the cornea (Kayser-Fleischer rings — the pathognomonic brown-gold deposits at the corneal periphery), and the kidneys (Fanconi syndrome). In a meta-analysis of 3,007 Chinese Wilson disease patients across 23 studies77 meta-analysis of 3,007 Chinese Wilson disease patients across 23 studies
Xue et al. 2023, Pediatr Neurol
, R778L carriers presented at a significantly younger age (SMD -0.18, p=0.0004) and had lower serum ceruloplasmin concentrations (SMD -0.21, p=0.03) than non-R778L carriers, consistent with the mutation's severe trafficking defect.

A knock-in mouse model carrying R778L confirmed the liver-to-brain disease pathway: by 3-5 months, mice show motor and cognitive dysfunction driven by hepatic copper overload that triggers systemic neuroinflammation (microglial activation, elevated cytokines), not by direct copper accumulation in brain tissue88 by 3-5 months, mice show motor and cognitive dysfunction driven by hepatic copper overload that triggers systemic neuroinflammation (microglial activation, elevated cytokines), not by direct copper accumulation in brain tissue
Dong et al. 2024, J Neuroinflammation; pharmacological inhibition of hepatic NF-κB normalized cognitive and motor performance, implicating the liver-brain inflammatory axis as the primary driver of neurological symptoms
.

Wilson disease is eminently treatable when diagnosed before end-organ damage sets in. EASL-ERN 2025 guidelines99 EASL-ERN 2025 guidelines
European Association for the Study of the Liver / European Reference Network clinical practice guidelines, J Hepatol 2025
recommend copper chelation (D-penicillamine or trientine) as first-line therapy for hepatic disease, with zinc salts for maintenance or neurological presentations. Liver transplantation is indicated for acute hepatic failure and is curative for the hepatic component of the disease.

Practical Actions

For heterozygous carriers (AC genotype): one functional ATP7B copy is sufficient for normal copper metabolism. No dietary copper restriction or medical treatment is needed. The clinical priority is family planning — if a reproductive partner also carries an ATP7B pathogenic variant, there is a 25% chance per pregnancy that a child will inherit two non-functional copies and develop Wilson disease. Carrier testing of partners and genetic counselling are strongly recommended before conception.

For homozygotes and compound heterozygotes (AA genotype or compound AA with another ATP7B pathogenic variant): Wilson disease should be assumed to be present or imminent. Specialist hepatology evaluation, confirmatory biochemical testing (serum ceruloplasmin, 24-hour urinary copper, liver copper), slit-lamp examination, and neurological assessment are urgent. Lifelong copper-lowering therapy must begin as early as possible — early treatment before symptomatic organ damage prevents all major complications.

Interactions

ATP7B R778L interacts with other pathogenic ATP7B variants in compound heterozygosity — the most common mode of Wilson disease in East Asian patients. Common compound genotypes include R778L/P992L and R778L/A874V. In the Zhang et al. 2022 cohort, R778L/A874V compound heterozygotes had later onset than R778L/R778L homozygotes, suggesting A874V is a partial-function allele. Carriers of R778L who are compound heterozygous with a severe truncating variant (frameshift, splice site) tend to present earliest. The related ATP7B variant rs201038679 (P992L) is the second most common Wilson disease mutation in East Asian populations; rs137853280 (c.1708-1G>C, a splice-acceptor variant) is a common pathogenic allele in European populations. Any of these, when compound-heterozygous with R778L, results in Wilson disease with severity depending on residual function of the trans-allele protein.

rs2922126

GHSR GHSR A/G (rs2922126)

Moderate Risk Factor

The Hunger Hormone's Receptor: When GHSR Expression Tips the Balance

Every time your stomach empties, it sends a chemical distress signal — a 28-amino-acid peptide called ghrelin — up the bloodstream and into the brain. Ghrelin11 Ghrelin
The "hunger hormone," secreted primarily by the stomach fundus before meals and suppressed after eating. Rises during caloric restriction and after weight loss.
docks onto its receptor, GHSR-1a (growth hormone secretagogue receptor type 1a), triggering appetite and promoting fat storage. But ghrelin does far more than tell you it's time to eat. GHSR-1a is also expressed in the brain's mesolimbic reward circuit22 mesolimbic reward circuit
Dopaminergic pathways connecting the ventral tegmental area to the nucleus accumbens, mediating motivation and pleasure responses to food, drugs, and other rewards.
, where it amplifies the desirability of food and drives motivated eating beyond simple caloric need. rs2922126 sits roughly 2 kilobases upstream of GHSR and influences the regulatory region that controls how much receptor is produced.

The Mechanism

rs2922126 is an upstream regulatory variant33 upstream regulatory variant
Located ~2 kb before the GHSR transcription start site, within the promoter/proximal regulatory region that controls GHSR gene expression.
— not a coding change, so the ghrelin receptor protein sequence is identical in all genotypes. The A allele is thought to alter transcription factor binding at the GHSR promoter, potentially shifting baseline receptor expression levels. Higher GHSR expression would amplify ghrelin signaling: stronger pre-meal hunger signals, increased food-reward sensitivity, and enhanced fat storage signaling. Nearby GHSR promoter variants with better-characterized function — particularly rs490683, which disrupts a nuclear factor 1 (NF-1) binding site — show that altered transcription factor access at this promoter44 altered transcription factor access at this promoter
Mager et al. 2008 demonstrated that the rs490683-GG genotype with an intact NF-1 site shows increased GHSR promoter activity in reporter assays.
has measurable metabolic consequences. rs2922126 likely operates through a similar mechanism, though its specific transcriptional effect has not been directly characterized.

The Evidence

The primary evidence for rs2922126 comes from a 2008 case-control study by Li et al.55 Li et al.
Ghrelin receptor gene polymorphisms are associated with female metabolic syndrome in Chinese population. Chinese Medical Journal, 2008.
enrolling 698 metabolic syndrome patients and 762 controls in China (total N=1,460). In women, the A/A genotype was associated with metabolic syndrome (OR 1.41, 95% CI 1.03–1.94), increased waist circumference (OR 1.75, 95% CI 1.26–2.42), and elevated fasting blood glucose (OR 1.49, 95% CI 1.07–2.06). These associations were not observed in men, pointing to a sex-hormone-mediated interaction between ghrelin signaling and central adiposity.

A 2013 prospective study by Yan et al.66 Yan et al.
Adiposity, inflammation, genetic variants and risk of post-menopausal breast cancer findings from a PRoBE design approach. SpringerPlus, 2013.
(180 cases, 732 controls) found that A/A homozygotes had a protective association against post-menopausal breast cancer (OR 0.4, 95% CI 0.18–0.89, p=0.02). This finding may reflect ghrelin's complex relationship with adipose tissue and estrogen signaling in post-menopausal women, though replication is needed before clinical weight is assigned to this result.

Null results also inform the picture. A longitudinal study of 1,362 children Riedl et al.77 Riedl et al.
GH secretagogue receptor gene polymorphisms are associated with stature throughout childhood. European Journal of Endocrinology, 2012.
found no association between rs2922126 and height or BMI across a 10-year follow-up, and a 2021 case-control study Tabaeian et al.88 Tabaeian et al.
NAFLD case-control study, N=310. Journal of Gastrointestinal and Liver Diseases, 2021.
found no association with non-alcoholic fatty liver disease. A systematic review Ghalandari et al.99 Ghalandari et al.
Review of 24 GHRL/GHSR case-control studies. Int J Endocrinol Metab, 2015.
concluded that the overall evidence across GHRL/GHSR polymorphisms and obesity is inconclusive. The metabolic syndrome signal from Li et al. has not been widely replicated in European or African ancestry populations.

Practical Actions

The A allele's primary actionable signal is elevated abdominal adiposity risk in women, mediated through amplified ghrelin/GHSR appetite signaling. Because GHSR drives hunger at the pre-meal level and amplifies food reward signals, targeted strategies that blunt pre-meal ghrelin spikes — rather than generic caloric restriction — are most aligned with the mechanism. High-protein meals are particularly effective at suppressing ghrelin compared to carbohydrate-matched meals, and distributing protein intake across three meals (rather than concentrating it at dinner) gives the most sustained pre-meal ghrelin suppression. For women with A/A genotype, waist circumference and fasting glucose are the most relevant biomarkers to track given the Li et al. findings.

Interactions

GHSR promoter variants rs490683 and rs9819506, studied by Mager et al. 2008 in the Finnish Diabetes Prevention Study, showed associations with body weight and glucose metabolism respectively. rs490683 has direct functional evidence (NF-1 binding site disruption) and may compound with rs2922126 if both variants shift GHSR expression in the same direction. rs572169, studied by Gueorguiev et al. 2009, was associated with obesity (OR 1.73, p=0.007) in an additive model — this is a distinct GHSR variant whose interaction with rs2922126 has not been studied. All four variants are within or near the GHSR promoter/upstream region and may tag related regulatory haplotypes.

IL12B Intronic Variant — Psoriasis Susceptibility and Ustekinumab Response Prediction

The IL12B gene11 IL12B gene
located at 5q33.3, encodes the 40 kDa p40 subunit shared by two functionally distinct cytokines: IL-12 (p40/p35 heterodimer) and IL-23 (p40/p19 heterodimer)
. IL-12 drives Th1 differentiation and IFN-γ production; IL-23 expands Th17 cells and IL-17 production. Both pathways converge on psoriatic skin inflammation and inflammatory bowel disease. rs3213094 is an intronic variant at chr5:159323761 (GRCh38) within the IL12B gene, lying in strong linkage disequilibrium with the established IL12B regulatory haplotype. While earlier literature referred to this region loosely as "promoter-associated," dbSNP classifies rs3213094 as an intron variant (NM_002187.3:c.89-432G>T). Its significance lies not in direct protein change but in its tight linkage with regulatory elements that govern IL12B expression — and, crucially, in its independent value as a pharmacogenomic predictor of ustekinumab response22 pharmacogenomic predictor of ustekinumab response
Ustekinumab (Stelara) is a monoclonal antibody targeting the p40 subunit encoded by IL12B, blocking both IL-12 and IL-23 simultaneously
.

The Mechanism

The T allele at rs3213094 tags a regulatory state of the IL12B locus associated with lower p40 expression. Genome-wide association studies in Chinese Han populations found the T allele to be protective against psoriasis (OR 0.78, combined p = 2.58×10⁻²⁶)33 protective against psoriasis (OR 0.78, combined p = 2.58×10⁻²⁶)
Zhang et al., Nature Genetics 2009, first large Chinese psoriasis GWAS confirming IL12B as a major susceptibility locus
, meaning the C allele is the risk-conferring allele. The protective T allele is associated with reduced IL12B transcriptional output in the regulatory environment defined by this haplotype block. This is pharmacogenomically consequential: individuals with one or two copies of the T allele have reduced baseline p40 availability, which changes the pharmacodynamic landscape for drugs targeting p40. For ustekinumab, the CT heterozygous genotype was associated with significantly better PASI reduction at 3 months compared to CC homozygotes in the Galluzzo 2016 study, suggesting that partial reduction of the IL12B risk-haplotype burden may allow ustekinumab to achieve more complete IL-12/23 pathway blockade.

The Evidence

The primary pharmacogenomic study is Galluzzo et al., Dermatology 201644 Galluzzo et al., Dermatology 2016
IL12B (p40) gene polymorphisms contribute to ustekinumab response prediction in psoriasis
, a retrospective Italian cohort of 64 patients treated with ustekinumab for up to one year. This study genotyped rs3213094 alongside rs2546890 and found that patients with the CT genotype at rs3213094 had significantly greater mean PASI reduction at 3 months compared to CC homozygotes (p = 0.017). The T allele was the favourable allele. The study also examined interactions with HLA-Cw6, the strongest single predictor of ustekinumab response, and found that IL12B genotype provided additional stratification beyond HLA-Cw6 status.

For psoriasis susceptibility, rs3213094 was confirmed in a two-stage Chinese Han GWAS55 two-stage Chinese Han GWAS
Zhang et al. 2009, 1,139 cases + 1,132 controls discovery; 5,182 cases + 6,516 controls + 539 Uygur cases replication
with a combined p-value of 2.58×10⁻²⁶ and an OR of 0.78 for the T allele — placing it among the most statistically robust SNPs at the IL12B locus. European studies of the IL12B haplotype report a complementary picture: the C allele (risk allele) confers elevated psoriasis odds across European populations, with estimates ranging from OR 1.4 to OR 1.9 depending on cohort and analysis model. A meta-analysis of 11 studies confirmed the IL12B risk haplotype in both psoriasis and psoriatic arthritis.

For inflammatory bowel disease, Brinar et al. 201366 Brinar et al. 2013
Multidimensional prognostic risk assessment identifies association between IL12B variation and surgery in Crohn's disease
found that IL12B variation was independently associated with need for early surgery within 5 years of Crohn's disease diagnosis, highlighting the shared pathogenic role of the IL12B locus across skin and gut inflammation. A meta-analysis of IL12B polymorphisms in IBD in Caucasian populations confirmed the association of the IL12B susceptibility haplotype with CD and UC risk.

Practical Actions

For individuals carrying the C allele (the risk genotype), two clinical contexts are most relevant. First, psoriasis susceptibility: the CC genotype marks elevated genetic risk for psoriasis and psoriatic arthritis through the same IL-12/23 axis that drives the condition. Second, and more immediately actionable, is biologic therapy planning: if ustekinumab is being considered for psoriasis or Crohn's disease, rs3213094 genotype provides independent pharmacogenomic information. The CT heterozygous genotype is associated with better early PASI response to ustekinumab; CC homozygotes appear to have attenuated benefit at 3 months, though response improves over longer treatment durations in many patients.

For individuals with one or two T alleles who develop psoriasis requiring biologics, this result supports ustekinumab as a rational first-line biologic choice — particularly in combination with HLA-Cw6 positivity, which additively improves response prediction. Guselkumab and risankizumab target only the p19 (IL-23-specific) subunit rather than the shared p40 encoded by IL12B, so their pharmacodynamic relationship with rs3213094 genotype differs from ustekinumab's.

Gut inflammation implications apply across genotypes: the IL12B locus's shared susceptibility with Crohn's disease and ulcerative colitis means carriers of the CC risk genotype with unexplained gastrointestinal symptoms warrant evaluation for subclinical IBD, mirroring the psoriasis-IBD comorbidity recognized in clinical practice.

Interactions

rs3213094 is in linkage disequilibrium with rs3212227 (3'-UTR IL12B) and rs6887695 (upstream IL12B), which together define the canonical IL12B psoriasis risk haplotype. Signals at rs3213094 and rs3212227 are not fully independent — they tag overlapping haplotype blocks — but rs3213094 has been studied as an independent pharmacogenomic marker for ustekinumab response (Galluzzo 2016), complementing the risk-susceptibility role of the haplotype SNPs. rs12188300 (also IL12B, near-gene regulatory) and rs3213094 together provide broader coverage of the IL12B susceptibility locus. rs11209026 (IL23R R381Q) is a protective loss-of-function receptor variant for IL-23; carrying both the IL12B CC risk genotype and IL23R protective A allele (rs11209026) creates a partially antagonistic combination where elevated p40 supply is partially offset by reduced receptor signaling downstream.

PRF1 A91V — The Perforin Pore Puncher That Falls Short

Perforin is the immune system's master assassin molecule. When a cytotoxic T lymphocyte (CTL) or natural killer (NK) cell locks onto a virus-infected cell or tumor cell, it releases perforin from secretory granules. Perforin polymerizes on the target cell membrane, punching transmembrane pores11 transmembrane pores
Perforin monomers oligomerize into ring-shaped channels ~10 nm in diameter, allowing granzymes to enter and trigger apoptosis
that let granzymes flood in and trigger programmed cell death. Without effective perforin, the immune system cannot terminate its own killing response — and that failure can spiral into a life-threatening inflammatory storm called hemophagocytic lymphohistiocytosis (HLH).

The A91V variant (c.272C>T; rs35947132) substitutes valine for alanine at position 91 of the perforin precursor protein. It is one of the most common functional variants in the PRF1 gene, occurring in approximately 4–9% of chromosomes in European populations and far less often in East Asian and African populations.

The Mechanism

Position 91 sits in the EGF-like domain of the perforin precursor, a region critical for protein folding and post-translational maturation22 protein folding and post-translational maturation
Perforin is synthesized as a 67-kDa precursor that must undergo calcium-dependent conformational change and N-glycosylation before becoming lytically active
. The A91V substitution disrupts an antigenic epitope recognized by the monoclonal antibody deltaG9, a sign that the local tertiary structure is altered. More importantly, the mutation impairs cleavage of the precursor to the active form33 impairs cleavage of the precursor to the active form
Patient cells expressing A91V show only immature and intermediate perforin species; the mature active form is absent or barely detectable on Western blot
. This processing defect translates into a 10-fold reduction in target-cell lysis when purified A91V protein is tested on cells, and roughly a 50% reduction in intact CTL cytotoxicity — because wild-type perforin on the second allele partially compensates in heterozygotes.

The functional deficit is bimodal. Voskoboinik et al. (Blood, 2007)44 Voskoboinik et al. (Blood, 2007) showed defects at both the presynaptic level (reduced perforin secretion into the immunological synapse) and the postsynaptic level (reduced membrane permeabilization even when secretion occurs). The 2015 study by House et al.55 House et al. confirmed these findings translate to primary human NK cells from healthy A91V heterozygous volunteers — a ≥35% reduction in cytotoxicity relative to wild-type individuals, measured in cells that were never exposed to inflammatory conditions.

The Evidence

Functional evidence is strong and replicated across multiple independent laboratories. A91V reduces perforin lytic activity, reduces expression, and causes measurable NK cell cytotoxicity deficits even in healthy heterozygous carriers.

Clinical evidence focuses on compound heterozygosity. In a landmark case series by Clementi et al. (2002)66 Clementi et al. (2002), siblings with adult-onset HLH were found to carry A91V on one allele and a null mutation (W374X) on the other — classic compound heterozygosity. Zhang et al. (Blood, 2011)77 Zhang et al. (Blood, 2011) found A91V in 48% of 25 adult-onset familial HLH patients, many presenting in middle age or later. Documented A91V combinations with W374X, G149S, R104C, and I125T all resulted in HLH with later and milder presentations than pediatric null-allele homozygotes.

Homozygous A91V is an unusual presentation. Mancebo et al. (Haematologica, 2006)88 Mancebo et al. (Haematologica, 2006) described an adult patient who was homozygous for A91V and developed FHL2 triggered by active tuberculosis infection — the first documented A91V homozygous FHL case. Crucially, the patient's monozygotic twin, with identical PRF1 genotype, remained completely healthy, demonstrating that homozygous A91V does not cause disease on its own: environmental triggers matter.

In FHL patient cohorts, Busiello et al. (2006)99 Busiello et al. (2006) found A91V in 26.2% of FHL patients carrying other PRF1 mutations, versus 3.7% in healthy controls (P=0.0002), firmly establishing its role as a disease-modifying susceptibility factor rather than a primary pathogenic mutation.

Beyond HLH, A91V has been detected at elevated frequency in NK/T-cell lymphomas — particularly nasal-origin cases (25% prevalence)1010 NK/T-cell lymphomas — particularly nasal-origin cases (25% prevalence) — and in young severe COVID-19 patients who developed HLH-like hyperinflammatory syndrome. A small study found A91V in 2 of 22 young critical COVID-19 patients (both died), with markedly elevated ferritin suggesting a virally triggered HLH phenotype.

Practical Implications

For heterozygous carriers with a single A91V allele and no known second PRF1 variant, the primary practical implication is awareness: persistent high fever, unexplained cytopenia, and dramatically elevated ferritin are warning signs that warrant urgent HLH evaluation. Standard genetic panels for HLH should include full PRF1 sequencing to identify compound heterozygosity.

For confirmed compound heterozygotes (A91V plus any other PRF1 pathogenic variant) or homozygous A91V individuals, surveillance with serum ferritin during febrile illnesses and avoidance of immunosuppression without specialist oversight are critical. Early recognition and treatment of HLH are the main determinants of survival.

Interactions

The A91V allele interacts with any second loss-of-function PRF1 allele to produce compound heterozygosity — the most clinically significant interaction. Documented combinations include A91V/W374X, A91V/G149S, A91V/R104C, and A91V/I125T, all causing late-onset HLH with milder phenotype than pediatric biallelic null mutations.

A91V also interacts with genes in the exocytosis pathway. Zhang et al. (2011)1111 Zhang et al. (2011) documented adult HLH patients carrying A91V in PRF1 alongside pathogenic variants in MUNC13-4 (UNC13D) or STXBP2, suggesting that cumulative functional deficits across the cytotoxic killing pathway can produce HLH even without biallelic PRF1 mutations.

rs3750243

EIF4EBP1 EIF4EBP1 mTOR Pathway Variant

Strong Risk Factor

EIF4EBP1 — The mTOR Gatekeeper Controlling How Quickly Your Ovaries Age

Inside every ovary sits a dormant population of primordial follicles — the finite reserve of egg precursors you are born with. How quickly this pool depletes determines when your fertility naturally declines and when menopause arrives. One of the most important molecular gatekeepers of this process is 4E-BP111 4E-BP1
eukaryotic translation initiation factor 4E-binding protein 1; encoded by EIF4EBP1 on chromosome 8p11.23; a direct downstream effector of mTORC1 that, when unphosphorylated, inhibits cap-dependent protein synthesis and suppresses cell proliferation and growth
. When mTOR activity is low, 4E-BP1 remains active and keeps primordial follicles quiescent. When mTOR is highly active, 4E-BP1 is phosphorylated and inactivated — follicles start waking up and the reserve depletes faster.

The rs3750243 variant sits approximately 2 kilobases upstream of EIF4EBP1 in a regulatory region. The C allele at this position is associated with later age at menopause and higher anti-Müllerian hormone (AMH) levels — two direct measures of ovarian reserve — suggesting it acts on EIF4EBP1 expression or activity in a way that keeps the mTOR brake engaged for longer.

The Mechanism

mTOR (mechanistic target of rapamycin) integrates nutrient, energy, and growth-factor signals to decide whether a cell should grow, divide, or remain quiescent. In primordial follicles, the PI3K–Akt–mTORC1 axis22 PI3K–Akt–mTORC1 axis
phosphoinositide 3-kinase activates Akt, which phosphorylates and activates mTORC1; mTORC1 then phosphorylates S6K1 and 4E-BP1, promoting cell growth and proliferation; excessive mTORC1 activation in oocytes accelerates follicle recruitment beyond the sustainable rate
is one of the principal determinants of follicle activation rate.

4E-BP1 sits directly downstream of mTORC1. In its unphosphorylated (active) form, it binds eIF4E and blocks cap-dependent mRNA translation — effectively slowing cell growth and keeping follicles dormant. mTORC1 phosphorylates 4E-BP1 at multiple residues, releasing eIF4E and allowing translation to proceed, which promotes follicle activation. The regulatory variant rs3750243 upstream of EIF4EBP1 is hypothesised to influence how readily this brake can be applied: the C allele may sustain higher baseline 4E-BP1 expression or activity, tipping the balance toward follicle quiescence.

The animal-model evidence for mTOR as an ovarian-reserve regulator is strong. Zhang et al. 201333 Zhang et al. 2013
Rapamycin preserves the follicle pool reserve and prolongs the ovarian lifespan of female rats. Eur J Obstet Gynecol Reprod Biol, PMID 23566837
showed that rapamycin (the prototypic mTOR inhibitor) administered every other day for 10 weeks doubled the primordial follicle count in rat ovaries relative to controls (p<0.001), with corresponding decreases in phospho-mTOR and phospho-p70S6K and increases in the longevity regulators SIRT1 and SIRT6. A separate study by Chen et al. 202244 Chen et al. 2022
Rapamycin maintains the primordial follicle pool and protects against chemotherapy-induced damage, PMID 35718464
confirmed that rapamycin inhibited excessive follicle activation and reduced apoptosis in mice exposed to chemotherapy.

The Evidence

rs3750243 emerged as one of the strongest signals in the landmark Ruth et al. 2021 Nature GWAS55 Ruth et al. 2021 Nature GWAS
Genetic insights into biological mechanisms governing human ovarian ageing; 201,323 European-ancestry women; 290 independent signals at genome-wide significance; PMID 34349265
— a study of 201,323 European-ancestry women powered to detect even modest genetic effects on age at menopause. The EIF4EBP1 locus at 8p11.23 showed an exceptional p-value of 5×10⁻²⁶¹ with a beta of +0.376 years per C allele. This means each copy of the C allele is associated with roughly 4.5 additional months of reproductive lifespan — and women who are CC homozygous carry approximately 9 months' advantage over GG homozygotes on average.

The same genomic region influences the most clinically actionable biomarker of ovarian reserve. Pujol-Gualdo et al. 2024 Human Reproduction66 Pujol-Gualdo et al. 2024 Human Reproduction
GWAS meta-analysis of AMH in 9,668 pre-menopausal women; three novel AMH loci including EIF4EBP1; PMID 38815977
identified EIF4EBP1 as one of three novel genome-wide significant AMH loci, meaning that variants in this region measurably alter circulating anti-Müllerian hormone — the glycoprotein secreted by small antral follicles that is used clinically to estimate the remaining follicle pool. The association was further replicated by Erdogan-Yildirim et al. 2025 Genes77 Erdogan-Yildirim et al. 2025 Genes
AMH GWAS in 1,185 Samoan women; replication of EIF4EBP1 AMH locus; PMID 40725450
across an independent non-European cohort, establishing cross-ancestry relevance.

EIF4EBP1's position as a direct mTOR effector gives this locus a compelling mechanistic interpretation. Unlike many GWAS hits that tag unknown regulatory elements near genes of uncertain relevance, 4E-BP1 has a well-characterised role in the exact pathway (mTOR) that animal models have repeatedly shown governs primordial follicle activation rate.

Practical Actions

For women who carry the GG genotype — the most common configuration — the practical implication is awareness: mTOR pathway activity is a meaningful lever on ovarian ageing, and interventions that support the 4E-BP1 brake (mTOR modulation) may have particular relevance. Caloric restriction, time-restricted eating, and CoQ10 supplementation all have mechanistic connections to mTOR activity and mitochondrial function in the ovary.

AMH testing provides a direct window into how this genetic variant is expressing itself in your specific physiology. A single fasting AMH level at any point in the cycle gives a reasonable estimate of remaining follicle pool; trajectories (measured 12–18 months apart) are more informative than a single value.

Interactions

The EIF4EBP1 rs3750243 locus sits in the same mTOR pathway context as other longevity-related mTOR variants. rs2295080 in MTOR itself (the kinase that phosphorylates 4E-BP1) and rs3803304 in PIK3CA (the upstream PI3K activating Akt and mTOR) both influence pathway activity from upstream nodes. The combined picture of rs3750243 (the brake, 4E-BP1) plus upstream mTOR or PI3K variants may determine the net mTOR activity level in ovarian tissue more precisely than any single variant. No compound effect sizes have been published for this combination.

rs2268361 in FSHR (follicle-stimulating hormone receptor) sits at a mechanistically adjacent node: FSH signalling through FSHR activates the same PI3K cascade that feeds into mTOR. FSHR variants affecting receptor sensitivity interact with the mTOR pathway at the point of signal input, while EIF4EBP1 variants affect signal output. The two loci likely modulate ovarian reserve through overlapping but distinct mechanisms.