Primary hyperoxaluria (PH) is a rare autosomal recessive disease characterized by oxalate accumulation in the kidneys and other organs. Three loci have been identified: AGXT (PH1), GRHPR (PH2), and HOGA1 (PH3). Here, we compared genotype to phenotype in 355 patients in the Rare Kidney Stone Consortium PH registry and calculated prevalence using publicly available whole-exome data. PH1 (68.4% of families) was the most severe PH type, whereas PH3 (11.0% of families) showed the slowest decline in renal function but the earliest symptoms. A group of patients with disease progression similar to that of PH3, but for whom no mutation was detected (11.3% of families), suggested further genetic heterogeneity. We confirmed that the AGXT p.G170R mistargeting allele resulted in a milder PH1 phenotype; however, other potential AGXT mistargeting alleles caused more severe (fully penetrant) disease. We identified the first PH3 patient with ESRD; a homozygote for two linked, novel missense mutations. Population analysis suggested that PH is an order of magnitude more common than determined from clinical cohorts (prevalence, approximately 1:58,000; carrier frequency, approximately 1:70). We estimated PH to be approximately three times less prevalent among African Americans than among European Americans because of a limited number of common European origin alleles. PH3 was predicted to be as prevalent as PH1 and twice as common as PH2, indicating that PH3 (and PH2) cases are underdiagnosed and/or incompletely penetrant. These results highlight a role for molecular analyses in PH diagnostics and prognostics and suggest that wider analysis of the idiopathic stone-forming population may be beneficial.
PH3 is the least severe form, with good preservation of kidney function in most patients. The typical presentation is recurrent urolithiasis and marked hypercalciuria in the first decade, but less active stone formation later.1,2,36,37 HOGA1 encodes the liver-specific mitochondrial enzyme 4-hydroxy-2-oxogluterate aldolase (HOGA), and mutations cause hydroxy-2-oxogluterate aldolase build-up, inhibiting GR/HPR function.11,38 PH3 accounts for approximately 10% of genetically characterized cases, with some carriers found to be idiopathic stone formers, suggesting sensitivity to haploinsufficiency.1,39 There are 19 described mutations with c.700+5G>T accounting for about 50% of all HOGA1 alleles; p.E315del is found predominantly in Ashkenazi Jews.22,36,37,39 No allelic correlations have been established for PH2 or PH3.
Here, we evaluated genotype-phenotype correlations at the genic and allelic level using a large collection of PH patients from the Rare Kidney Stone Consortium Primary Hyperoxaluria registry (RKSC PH registry). In addition, we used publicly available whole-exome sequencing data to provide population-based estimates of PH prevalence. Overall, these studies provide new insight into PH phenotypes and the significance of particular alleles. In addition, they indicate a carrier rate far higher than estimated from clinical populations, suggesting a significant level of underdiagnosis and/or incomplete penetrance/variable expressivity.
We mutation-screened the three known PH loci in 301 PH families (355 patients) using Sanger sequencing. Of these, 68.4% (206 pedigrees, 247 patients) were PH1, 9.3% (28 pedigrees, 35 patients) were PH2, and 11.0% (33 pedigrees, 38 patients) were PH3. Two mutant alleles were detected in each case. The remaining 11.3% (34 pedigrees, 35 patients) had a clinical phenotype consistent with PH according to the RKSC PH registry entry criteria (Concise Methods), yet no mutation was detected (NMD) in the known genes.
In total, we identified 121 different mutations within this cohort (83 AGXT, 21 GRHPR, and 17 HOGA1), of which 36 have not been previously described in the Human Gene Mutation Database (HGMD 2013.3 Professional) or the PH Mutation Database22 (PHMD; 18 AGXT, 13 GRHPR, and 5 HOGA1) (Supplemental Figure 1, Table 1). One of the novel mutations, GRHPR c.[-4G>A, -3C>T], is upstream of the transcription start codon generating a frameshifted novel start site with a stronger Kozak consensus (Supplemental Figure 2A).
The genotypic and allelic breakdown varied considerably between the different genes (Figure 1). For instance, 50.0% of PH2 pedigrees contained two truncating mutations, an allelic combination that accounted for only the minority of PH1 genotypes (14.1%) and was absent in PH3 (Figure 1A). GRHPR was more prone to frameshifting InDels (insertion, duplication, deletion, or insertion+deletion), accounting for 44.6% of all mutant alleles, whereas most AGXT mutations were missense (67.0%) (Figure 1B). These enrichments were driven by common mutations; the AGXT mistargeting change p.G170R was found in 50.5% of PH1 families and GRHPR c.103delG in 35.7% of PH2 families. HOGA1 showed considerably less allelic heterogeneity with the two most common alleles (c.700+5G>T and p.E315del), accounting for 74.2% of the total (Figure 1B).
Genotypic and allelic breakdown varies considerably between the different PH types in the mutation resolved 267 pedigrees. (A) Genotypic breakdown by PH type showing patients with two truncating alleles (nonsense, splice, and frameshifting InDels [insertion, duplication, deletion, or insertion+deletion]), two nontruncating alleles (missense and inframe InDels), or a truncating plus a nontruncating allele. Number in parentheses represents number of families with that genotype. (B) Allelic analysis by PH type of different mutation types. Common alleles of each PH type plus their frequencies are also highlighted. AGXT p.G170R was found homozygously (hom) in 29 families and heterozygously (het) in 75 families, GRHPR c.103delG (hom, 5 families; het, 5 families), HOGA1 c.700+5G>T (hom, 8 families; het, 12 families), HOGA1 p.E315del (hom, 7 families; het, 7 families).
Renal Survival plots showing (A) poorer renal survival for PH1 patients followed by PH2, and (B) better renal survival for PH1 patients with two MiR alleles. (A) Kaplan-Meier renal survival plot of the PH1, PH2, PH3, and NMD cohorts. (B) Kaplan-Meier renal survival plot of PH1 patients categorized as having two MiR alleles (homozygous, or compound heterozygous; AGXT p.G41R, p.F152I, p.G170R, p.I244T, 89 patients), one MiR allele in combination with a non-MiR allele (97 patients), or two alleles of which neither was a MiR allele (no MiR, 81 patients). Tables below Kaplan-Meier plots show survival estimates with number of patients at risk in parentheses.
Renal survival plots showing that only homozygous p.G170R patients but not homozygotes of other MiR genotypes had a renal survival advantage. Kaplan-Meier renal survival plot plus survival estimate table for PH1 patients homozygous for p.G170R (34 patients), homozygous or compound heterozygous for another MiR allele (p.G41R, p.F152I, p.I244T, 19 patients), or having no MiR allele (81 patients). Pairwise comparisons show greater renal survival in patients homozygous for AGXT p.G170R but not the other MiR alleles compared with patients without a MiR allele (Hom p.G170R versus no MiR; P<0.0001 [P<0.001], Hom other MiR versus no MiR P=0.27 [P=0.35], Hom p.G170R versus Hom other MiR P=0.04 [P=0.05]). P values in square brackets are adjusted for the effect of pyridoxine treatment, reducing the overall significance but upholding the survival advantage of AGXT p.G170R homozygotes.
To evaluate genotype-phenotype correlations in our PH2 cohort, patients were grouped according to mutation type (two truncating, two nontruncating, or a truncating plus a nontruncating allele). However, no significant associations were observed with age at symptoms or ESRD, nephrocalcinosis, and urine chemistry (Supplemental Table 3).
Patients with PH3 were grouped by the common mutations (homozygous for c.700+5G>T or p.E315del; heterozygous, compound heterozygous for c.700+5G>T/p.E315del, or lacking these mutations). As for PH2, no significant genotype-phenotype associations in our relatively small PH3 cohort were apparent (Supplemental Table 4). Interestingly, one PH3 patient developed ESRD at 8 years. This patient was homozygous for a pair of closely linked, nonconservative missense variants predicted to be highly pathogenic, rather than the common HOGA1 alleles (Supplemental Figure 2B, Table 1). This patient was symptomatic at age 18 months, with multiple calcium oxalate stones that were surgically removed; CKD was first noted at age 6 years, at least partly attributed to obstructive uropathy.
Intrafamilial variable disease expression and phenotypic heterogeneity among patients with the same allelic combinations have been observed in PH1.5,25 In the current cohort, five PH1 families had differences of >20 years between siblings regarding the timing of ESRD onset, and considerable phenotypic variability was present between patients homozygous for specific AGXT alleles (Supplemental Table 5A). Within the PH2 and PH3 cohorts, homozygotes for common mutations also displayed considerable variability of disease expression (Supplemental Table 5A). Within one PH2 family homozygous for GRHPR c.103delG, one of three affected siblings reached ESRD at age 21 years, whereas the other two manifested with only occasional stone disease, with eGFRs of 40 and 85 ml/min per 1.73 m2 at 56 and 53 years, respectively.
From the National Heart, Lung, and Blood Institute (NHLBI) Exome Sequencing Project (ESP) data, we extracted allelic counts of 22 known PH mutations (12 AGXT, 5 GRHPR, and 5 HOGA1) found within the population of 4300 European Americans (EA) and 2203 African Americans (AA) (Supplemental Figure 1, Supplemental Table 6).40 Further, we included all frameshifting or nonsense changes and scored all rare, nonsynonymous variants not previously identified as PH mutations using in silico tools, and thus identified 20 additional likely pathogenic variants (4 AGXT, 6 GRHPR, and 10 HOGA1) (Supplemental Figure 1, Supplemental Table 7).
With use of the Hardy-Weinberg equation, an overall PH carrier frequency (CF) of 1:71 with a predicted prevalence of 1:58,243 was calculated from the 22 known mutant alleles (Supplemental Table 6, Table 4). Including the predicted pathogenic alleles increased the CF and prevalence to 1:58 and 1:38,630, respectively (Supplemental Tables 7 and 8). PH was predicted to be approximately 60% less prevalent in AAs than in EAs (EA: CF, 1:61, prevalence, 1:40,276; AA: CF, 1:104, prevalence, 1:100,380) (Table 4). This difference was largely driven by the common EA allele HOGA1 c.700+5G>T (CF, 1:165; AA: CF, 1:1102) (Figure 4, Supplemental Table 6). The prevalence of PH1 was similar in EAs and AAs, even though AGXT p.G170R (EA: CF, 1:429) was unique to the EA population, but was balanced by AGXT p.R289H, the most prevalent AA allele, which accounted for approximately 33% of all AA mutant alleles (CF, 1:289; EA: CF, 1:1956) (Figure 4, Supplemental Table 6). Interestingly, the AGXT “minor” allele variants p.P11L/p.I340M were 15 times less prevalent in AAs than EAs (EA: CF, 1:3, prevalence, 1:15; AA: CF, 1:8, prevalence, 1:236) (Supplemental Table 6). As with PH1, PH2 prevalence was also similar in EAs versus AAs as GRHPR c.103delG was relatively common in both populations (EA: CF, 1:375; AA: CF, 1:426) (Figure 4, Supplemental Table 6).
Carrier frequencies of common PH alleles found in the NHLBI ESP account for significant racial differences. Depiction of the most common PH alleles separated by EA and AA CF. AGXT p.G170R (CF in EAs, 0.23%) and p.R289H (CF in AAs, 0.35%) account for 45.1% and 68.8% of the total PH1 CF, respectively. Similarly, GRHPR c.103delG accounts for 75.3% and 63.4% of the total EA or AA PH2 CF, respectively, and HOGA1 c.700+5G>T accounts for 78.4% and 100% of the total EA or AA PH3 CF, respectively.
When separated by PH type, PH1 and PH3 had a similar CF (1:195 versus 1:185) and prevalence (1:151,887 versus 1:135,866) (Table 4), whereas in our patient screen PH1 was six times more common (68.4% of families versus 11.0%). These data suggest a distinct underdiagnosis and/or lack of penetrance of PH3. The frequency of PH3 was predominantly driven by HOGA1 c.700+5G>T, which had an approximately 2.8 times greater CF (1 in 232) than AGXT p.G170R (1 in 649) in the ESP data, yet they accounted for 5.2% and 24.9% of all PH alleles, respectively, within our patient cohort. On the basis of these calculations, PH2 was predicted to be about 2.2 times less prevalent than PH1 or PH3 (CF, 1:279; prevalence, 1:310,055) (Table 4), but in our disease cohort its prevalence was similar to that of PH3 and approximately seven times rarer than PH1.
Here we present genotype-phenotype correlations in one of the largest fully genotyped PH cohorts to date and for the first time calculate PH prevalence based on large population data. Because of the population size, we have significantly increased the number of unique PH mutations (AGXT by 10.1%, GRHPR by 46.4%, and HOGA1 by 26.3%) (Supplemental Figure 1). This analysis highlights the allelic heterogeneity of these disorders, although the population importance of a few common alleles is also reinforced (Figure 1). While elevated urinary glycolate,5 l-glycerate,4 or HOGA38 levels can often provide strong evidence regarding the PH gene involved (Table 2), there is significant overlap between PH types and normal individuals. Thus, we emphasize here the value of molecular analyses as a gold standard method to diagnose PH. As sequencing costs fall and new methods are adopted, it is increasingly practical to implement a protocol that fully screens the coding regions of all genes, not just commonly inherited variants,5 and thoroughly evaluates all detected variants.
Our genic prevalence in mutation resolved families—77.1% PH1, 10.5% PH2, 12.4% PH3—coincides with previous estimates (PH1, 80%; PH2, 10%; PH3, 10%) and highlights PH1 as the most common form, while PH3 is now at least as common as PH2.1–3 Interestingly, 11.3% of our clinically defined cohort (34 pedigrees) had no mutation detected in the known PH genes. Given the phenotype in these families (see below) and the reliability of molecular diagnostics (with no single PH mutations detected in our NMD group), further genetic heterogeneity in PH seems likely. Although our NMD cases met the RKSC PH registry criteria (Concise Methods), first symptoms occurred 4–9 years later than for genetically defined PH (Table 2), urine oxalate excretion was lower than PH1 and PH2 (but similar to PH3), and ESRD was rare, although there was considerable interfamilial variability (Supplemental Table 5B). These data indicate that rigorous clinical/molecular assessment is required to identify a population suitable for further gene discovery efforts in PH.
We found no difference in phenotype between patients with two non-MiR nontruncating versus two truncating alleles, suggesting that all these AGXT alleles are fully inactivating (Supplemental Figure 3, Supplemental Table 2). This contrasts with a recent publication that reported better renal survival of patients with two non-p.G170R nontruncating variants versus two truncating variants (16.9 years versus 9.9 years).47 Although it is possible we did not observe this correlation because of our smaller cohort size (247 patients in the RKSC PH registry versus 410 patients in the OxalEurope Consortium), and the significance of pairwise comparisons was not provided.47 The OxalEurope investigators also analyzed allelic effects between specific missense variants,47 but given the observed intrafamilial/intragenotype phenotypic heterogeneity (Supplemental Table 5A)5,25,47 due to genetic/environmental modifier effects, caution is required with interpretation in these small populations. More detailed allelic analysis of PH1 will be possible as population sizes increase and through pooled analysis of available cohorts.
In our cohort, PH2 was similar to PH1 in terms of age at symptoms and oxalate levels; however, progression to ESRD was significantly slower (Figure 2, Table 2).2,4 Elevated urinary l-glycerate levels separated PH2 cases from PH1 and PH3 cases, but some overlap with the genetically unresolved group was seen (minimum, maximum urine l-glycerate mg/g creatinine: PH1, 0.0, 129.0; PH2, 151.0, 4355.0; PH3, 0.0, 40.0; NMD, 0.0, 152.0).4 Similarly, dramatic intrafamilial/intragenotype phenotypic heterogeneity was observed (Supplemental Table 5A), suggesting significant genetic/environmental modifier effects in PH2. No significant associations of phenotype with mutation type were detected in patients with PH2 (Supplemental Table 3), although our cohort was probably too small to identify subtle differences in disease severity. Interestingly, most GRHPR mutant alleles were frameshifting, largely driven by the common c.103delG variant.
In our cohort, no PH3 family had two truncating alleles (Figure 1), with only one such patient described.36 Phenotypic analysis of patients with varying or no detected HOGA levels will prove interesting and provide insight into whether most HOGA1 alleles are incompletely or fully penetrant, consistent with structural analyses.48 Interestingly, some HOGA1 carriers present with mild hyperoxaluria or idiopathic urinary stone disease39 (unlike other forms of PH), and the first reported patient with PH3 to reach ESRD was homozygous for two novel, linked missense variants (Supplemental Figure 2); suggesting that more penetrant mutations may have a stronger phenotypic effect. However, the ESRD in this single patient with PH3 may have, at least in part, been aggravated by renal injury from recurrent stone removal surgeries. Paradoxically, patients with PH3 present with symptoms at the earliest age but subsequently have relatively mild disease (Figure 2, Table 2), suggesting decreased sensitivity to HOGA levels with aging. The lack of a detected genotype-phenotype correlation between the two most common HOGA1 alleles (c.700+5G>T and p.E315del) may again be due to the small populations and other genetic/environmental modifier effects (Supplemental Figure 5A).
The most important finding from this study is that the overall carrier frequency (approximately 1:70) and the inferred prevalence (approximately 1:58,000) calculated from population data are close to an order of magnitude greater than estimates from clinical studies. Because these new estimates come from counting known mutant alleles and are largely driven by a few common alleles, this study is unlikely to overestimate the carrier frequency. The estimates of PH1 prevalence are about twice the previous estimates, but the difference is much more notable for PH3, which has an observed population carrier frequency higher than PH1 but is 6-fold less common in currently identified clinical populations. PH2 has a population carrier frequency half that of PH1 but is more than 7-fold less common in clinical populations. It is notable that PH is determined to be about 2.5-fold less prevalent in AAs than EAs, a difference mainly driven by the HOGA1 c.700+5G>T mutation (Figure 4). Interestingly, AGXT p.G170R is exclusively found in EAs, and the AGXT MTS variants (p.P11L, p.I340M, minor haplotype) are much rarer in AAs than EAs (Figure 4, Supplemental Table 6).49 This finding suggests that mutations predicted to mistarget AGT might generally remain asymptomatic in AAs because presence of the minor haplotype is required for the MTS to be active.29 Interestingly, PH1 prevalence is nonetheless similar between EAs and AAs because of the most common AA mutation, AGXT p.R289H (reported in the PHMD22). This variant is significantly more common than any PH allele (excluding EA: HOGA1 c.700+5G>T) (Figure 4) and scores poorly,40 calling into question its true pathogenic nature. Removal of this allele would decrease PH1 prevalence by 48% (CF, 1:270; prevalence, 1:291,256) and overall PH prevalence by 18% (CF, 1:79; prevalence, 1:71,333, Supplemental Table 9). Further, PH1 prevalence in AAs would decrease by 87% (CF, 1:549; prevalence, 1:1,203,908) and PH prevalence overall within the AA population by 56% (CF, 1:156; prevalence, 1:225,971). Consequently, PH would be about 82% less prevalent in AAs than EAs.
The difference between the expected and observed prevalence, especially for PH3 and to a lesser extent for PH2, is likely due to underdiagnosis of these diseases, which have overall milder phenotypes and are much less likely than PH1 to result in ESRD. In addition, relatively mild disease in a subset of PH1 patients (especially p.G170R homozygotes) and a wide variability in age at first symptoms (Table 2), may account for the 2-fold underdiagnosis of PH1. Further, misdiagnosis of PH for other renal diseases, as illustrated by homozygosity for an AGXT mutation in a patient with nephronophthisis,50 may also play a role in underestimation of PH prevalence. Therefore, PH genotype analysis of idiopathic stone disease populations and other patients with PH phenocopies but no known mutations are likely to be informative for undiagnosed PH cases. Indeed, within a stone clinic setting the diagnosis of PH should be carefully considered among patients with unexplained hyperoxaluria and/or CKD, according to the published diagnostic algorithm.51 This consists of an initial comprehensive analysis of urine oxalate and chemistries plus plasma oxalate, followed by molecular testing of the PH genes in suspected cases.
The relevant institutional review boards and ethics committees approved the study, and all participants gave informed consent. The RKSC PH registry is a secure web-based registry that enables international contributions.3 Patients within the registry met at least one of the following criteria: (1) Liver biopsy confirmed specific enzyme deficiency of AGT, GR/HPR, or HOGA, (2) molecular diagnostic testing of AGXT, GRHPR, and HOGA1 identified two mutant alleles in one of the genes, (3) urinary oxalate excretion was >0.8 mmol/1.73 m2 per day in the absence of an identifiable secondary cause of hyperoxaluria, and/or (4) ESRD had developed, along with (a) predialysis plasma oxalate>60 μmol/L and renal biopsy confirming extensive oxalate deposition or (b) systemic oxalosis (retinal oxalate deposits; oxalate deposits in bone marrow, skin, or other tissue; nephrocalcinosis; or calcium oxalate nephrolithiasis). As of June 2014, 422 patients from 37 countries had been enrolled in the RKSC PH registry, of whom 355 were screened for mutations in the known PH genes.
Patients enrolled in the RKSC PH registry who consented to molecular testing had blood drawn for DNA isolation. Sanger sequencing was performed of all coding exons±20 bp (AGXT, NM_000030.2; GRHPR, NM_012203.1; HOGA1, NM_138413.3) using M13-tailed primers (Beckman Coulter). Primer sequences can be obtained by request. All Sanger chromatograms were analyzed using Mutation Surveyor, version 4.06 (Softgenetics), and identified variants were categorized as nonsense, missense, frameshifting, or inframe InDels; typical splice (±2 bp from exon boundary); or atypical splice (>±2 bp from exon boundary). A total of 115 families (approximately 38%) were screened for the coding regions of all three genes simultaneously, while the remainder were screened sequentially until two pathogenic alleles were identified (first AGXT, second GRHPR, third HOGA1). None of the NMD cases had a single pathogenic allele in any of the three genes, and none of the patients with PH1, 2, or 3 had a third pathogenic allele in the same or another PH gene.
Amino acid substitutions not previously reported in HGMD or PHMD22 were evaluated using the in silico prediction programs AlignGVGD (http://agvgd.iarc.fr/index.php), PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2/), and SIFT (http://sift.jcvi.org). The novel GRHPR atypical splice site mutation (c.404+5G>A) was evaluated using the Berkeley Drosophila Genome Project (http://www.fruitfly.org) with the DNA sequence of the altered ±2 exons as input. The novel GRHPR Kozak sequence mutation (c.[−4G>A, −3C>T]) was assessed using NetStart 1.0 (http://www.cbs.dtu.dk/services/NetStart/) with the full mRNA+5′ untranslated region as input. Additionally, for amino acid substitutions, multiple sequence alignments using AGT, GR/HPR, and HOGA protein orthologs of bull (Bos taurus), mouse (Mus musculus), rat (Rattus norvegicus), dog (Canis lupus), chicken (Gallus gallus), and zebrafish (Danio rerio) were used to evaluate evolutionary conservation. Amino acid substitutions were scored as pathogenic if two of three in silico tools predicted them as damaging, the change was evolutionarily conserved, and/or the variant segregated appropriately.
Phenotype analyses used the first available 24-hour urine (oxalate, calcium, citrate, glycolate, l-glycerate) and first renal imaging (nephrocalcinosis). More than 80% of these tests were done from 1 year before to 3 years after PH diagnosis. Comparisons of genotypic groups with continuous and nominal phenotypic factors at disease presentation were done using the Kruskal–Wallis rank test and chi-squared tests, respectively. Mixed-effects models (urine chemistry ranks) and generalized estimating equations (nephrocalcinosis) were also used to test genotypic effects while adjusting for possible pedigree effects. Time from birth to development of renal failure was estimated using Kaplan–Meier curves, with comparisons based on the log-rank test and censoring subjects without renal failure at last follow-up. The Cox model with robust variance estimator was used to test associations of phenotype with time to ESRD, with adjustment for potential within-pedigree correlations. The adjustment for the effect of pyridoxine treatment was performed using a time-dependent covariate in the Cox model. All statistical tests were two sided, with an α level of 0.05; only P values adjusted for potential within-pedigree effects as described above were reported.
The data of the NHLBI ESP40 were downloaded from http://evs.gs.washington.edu/EVS/ in VCF file format and analyzed using SNP & Variation Suit, version 7 (GoldenHelix). Variants within AGXT, GRHPR, and HOGA1 that were known as PH mutations (HGMD, PHMD, RKSC PH registry) were annotated (Supplemental Figure 1, Supplemental Table 6). All other variants with an MAF <0.01% were scored for potential pathogenicity using the criteria described above (Supplemental Figure 1, Supplemental Table 7). The only PH mutation reported as homozygous in ESP was GRHPR c.103delG. However, the variant is located in a poly-G stretch, which often results in misalignment; hence, genotype and allele counts were adjusted to represent this homozygous call as heterozygous. Using alternate and total allele counts, prevalence (1/p2) and CFs [1/(1−p2−q2)=1/(2pq)] were calculated according to Hardy-Weinberg equilibrium for each PH type using the sum of all alternate PH1, 2, or 3 alleles (p, known [Table 4] or known+scored as pathogenic [Supplemental Table 8]) and all WT alleles (1−p=q). Calculations for PH overall were done by adding the p2 values of each PH type, with bigenic PH individuals counted as carriers.
We thank all of the patients and families who have participated in the RKSC PH registry as well as the many physicians who collected the detailed clinical records (see below and Supplemental Material for the list of RKSC coordinators and participation/contribution sites, respectively). Further, we thank the study coordinators who collected clinical data and biological samples, as well as Ramila A. Mehta for her statistical advice.
The RKSC is a part of the National Institutes of Health Rare Diseases Clinical Research Network. Funding and/or programmatic support for this project has been provided by U54-DK083908 from the National Institute of Diabetes and Digestive and Kidney Diseases, National Center for Advancing Translational Sciences, Mayo Hyperoxaluria Center, and Mayo Foundation. This work was further supported by the Oxalosis and Hyperoxaluria Foundation and the American Society of Nephrology Ben J. Lipps Research Fellowship (K.H.).
This work was partially presented as two abstracts at the Annual Meeting of the American Society of Nephrology, 2013 (Atlanta, GA), and the 11th International Primary Hyperoxaluria Workshop, 2014 (Chicago, IL).
RKSC Coordinators and Coordination Sites: Craig Langman (Children's Memorial Hospital, Chicago), Lawrence Copelvitch (Children's Hospital of Philadelphia), Yaacov Frishberg (Shaare Zedek Medical Center, Jerusalem, Israel), Vidar Edvardsson (Landspitali University Hospital, Reykjavik, Iceland), David Goldfarb (New York University), Dean Assimos (University of Alabama, Birmingham), and Todd Lowther (Wake Forest University School of Medicine, Winston-Salem).
K.H. and A.G.C. contributed equally to this work.
This article contains supplemental material online at http://jasn.asnjournals.org/lookup/suppl/doi:10.1681/ASN.2014070698/-/DCSupplemental.