Abstract
Familial hypomagnesemia with hypercalciuria and nephrocalcinosis (FHHNC) is an ultra-rare autosomal recessive renal tubular disease with an incidence of <1/1.000.000 individuals, caused by loss-of-function mutations in CLDN16 and CLDN19. Our study includes a unique cohort representing all known FHHNC patients in Spain, with 90% harbouring mutations in CLDN19. Of these, 70% carry the p.G20D mutation in homozygosis. Despite this high genetic homogeneity, our FHHNC cohort display a high phenotypic variability, even among siblings harbouring identical mutations. Patients were stratified at the extremes of the renal phenotype according to their estimated glomerular filtration rate annual decline and subjected to whole exome sequencing (WES) aiming to find candidate phenotype-modifier genes. Initial statistical analysis by SKAT-O identified numerous variants, which were then filtered based on P-value <0.01 and kidney expression. A thorough prioritization strategy was then applied by an exhaustive disease knowledge-driven exploitation of data from public databases (Human Protein Atlas, GWAS catalog, GTEx) to further refine candidate genes. Odds ratios were also calculated to identify potential risk variants. This analysis pipeline suggested several gene variants associated with a higher risk of developing a more aggressive renal phenotype. While these findings hint at the existence of genetic modifiers in FHHNC, further research is needed to confirm their role and potential clinical significance. Clinical decisions should not be based on these preliminary findings, and additional cohorts should be studied to validate and expand upon our results. This exploratory study provides a foundation for future investigations into the genetic factors influencing FHHNC progression and may contribute to our understanding of the disease's variable expressivity potentially enabling the implementation of more tailored therapeutic strategies.