Variability in the severity of iron overload among homozygotes for the HFE C282Y polymorphism is one of the major problems extant in our understanding of hereditary hemochromatosis (HH). We conducted exome sequencing of DNA from C282Y homozygotes with markedly increased iron stores (cases) and C282Y homozygotes with normal or mildly increased iron stores (controls) to identify rare and common causal variants associated with variability of disease expression in HH. Criteria for cases included serum ferritin >1000 µg/L at diagnosis, and (a) mobilized body iron >10 g by quantitative phlebotomy, and/or (b) hepatic iron concentration >236 µmol/g dry weight. Criteria for controls included (a) serum ferritin <300 µg/L, or (b) age ≥50 y with ≤3.0 g iron removed by phlebotomy or age ≥40 y with ≤2.5 g iron removed by phlebotomy to achieve serum ferritin <50 µg/L. Deep sequencing of the full exome was performed in 33 cases and 14 controls. After quality control filtering, the dataset included 82,068 SNPs and 1,403 insertions/deletions (indels).Our initial analysis tested for differences in the distribution of variants between groups for each gene separately using the Sequence Kernel Association Test (SKAT) that includes rare and common variants but downweights the contribution of common variants to the test statistic. Only non-synonymous variants were included in the by-gene tests. Principal components were constructed from the exome variants to adjust for possible confounding by ancestry and to confirm no ancestral outliers.All study participants were male, and all clustered closely together within a larger group of Europeans in a principal components analysis of ancestry. Mean (SD) ages at presentation were 54 (11.0) y and 56 (9.4) y for cases and controls, respectively. Median serum ferritin was 2788 µg/L in those with increased iron stores and 309 μg/L in those with normal or mildly increased iron stores. The median transferrin saturation (94%) was greater in cases than in the comparison group (70%). In a preliminary analysis, we found 9 genes associated with case-control status. To separate effects of alcohol use and/or alcohol addiction variants, an analysis was conducted to compare the 13 controls and 22 cases who reported never using alcohol or only very light use. The two most significant genes identified in this comparison were GNPAT (p=7.4x10-6) and CDHR2 (p=2.8x10-4). A quantile-quantile (QQ) plot is shown in the Figure, comparing the observed distribution of –(log10p-values) from 10,337 genes to the expected uniform distribution if there were no variants modifying severity of expression, and gives evidence of the effect of the GNPAT gene. [Display omitted] Inspection of the two variants contributing to the GNPAT by-gene p-value revealed one missense variant (rs11558492) for which 0/13 controls had a polymorphism, while 16/22 cases had at least one missense variant, and one case was homozygous for this missense variant. The latter case presented at the early age of 26 with a serum ferritin of 1762 µg/L, 4+ hepatocellular iron and hepatic iron concentration of 284.4 µmol/g dry weight. GNPAT (aka DHAPAT) mutations/deletions have been found in peroxisomal disease, a class of diseases in which increased hepatic iron is observed (Biochim Biophys Acta 1801:272-280, 2010). GNPAT rs11558492 is common among people of European descent but might interact with aberrant HFE to increase risk of hepatic iron overload. Three rare variants in CDHR2 accounted for its low p-value, having a cumulative frequency of 4/13 among controls and 0/22 among cases: rs115050587, rs752138, rs143224505 with minor allele frequencies, MAF = 1.4%, 4.7% and 0.06%, respectively. The first two polymorphisms are predicted to be highly damaging by PolyPhen2 and the third probably damaging. Expression levels of CDHR2 recently have been associated with increased hepatocyte iron and elevated serum ferritin in liver allograft patients (J Clin Invest 122:368-382, 2012).These data indicate associations between iron status in HFE C282Y homozygotes and genes with previous links to iron overload that may modify severity of disease expression. Of note, the data suggest that more than one modifier gene may be involved in determining severity of disease in HFE C282Y homozygotes. Our results identify candidate genes for expanded studies that would examine their functional significance for iron absorption and metabolism. Disclosures:No relevant conflicts of interest to declare.