Population and evolutionary genetics of the PAH locus to uncover overdominance and adaptive mechanisms in phenylketonuria: Results from a multiethnic study.

BACKGROUND: Phenylketonuria (PKU) is the most common inborn error of amino acid metabolism in Europe. The reasons underlying the high prevalence of heterozygous carriers are not clearly understood. We aimed to look for pathogenic PAH variant enrichment according to geographical areas and patients' ethnicity using a multiethnic nationwide cohort of patients with PKU in France. We subsequently appraised the population differentiation, balancing selection and the molecular evolutionary history of the PAH locus. METHODS: The French nationwide PKU study included patients who have been referred at the national level to the University Hospital of Nancy, and for whom a molecular diagnosis of phenylketonuria was made by Sanger sequencing. We performed enrichment analyses by comparing alternative allele frequencies using Fisher's exact test with Bonferroni adjustment. We estimated the amount of genetic differentiation among populations using Wright's fixation index (Fst). To estimate the molecular evolutionary history of the PAH gene, we performed phylogenetic and evolutionary analyses using whole-genome and exome-sequencing data from healthy individuals and non-PKU patients, respectively. Finally, we used exome-wide association study to decipher potential genetic loci associated with population divergence on PAH. FINDINGS: The study included 696 patients and revealed 132 pathogenic PAH variants. Three geographical areas showed significant enrichment for a pathogenic PAH variant: North of France (p.Arg243Leu), North-West of France (p.Leu348Val), and Mediterranean coast (p.Ala403Val). One PAH variant (p.Glu280Gln) was significantly enriched among North-Africans (OR = 23·23; 95% CI: 9·75-55·38). PAH variants exhibiting a strong genetic differentiation were significantly enriched in the 'Biopterin_H' domain (OR = 6·45; 95% CI: 1·99-20·84), suggesting a balancing selection pressure on the biopterin function of PAH. Phylogenetic and timetree analyses were consistent with population differentiation events on European-, African-, and Asian-ancestry populations. The five PAH variants most strongly associated with a high selection pressure were phylogenetically close and were located within the biopterin domain coding region of PAH or in its vicinity. Among the non-PAH loci potentially associated with population divergence, two reached exome-wide significance: SSPO (SCO-spondin) and DBH (dopamine beta-hydroxylase), involved in neuroprotection and metabolic adaptation, respectively. INTERPRETATION: Our data provide evidence on the combination of evolutionary and adaptive events in populations with distinct ancestries, which may explain the overdominance of some genetic variants on PAH. FUNDING: French National Institute of Health and Medical Research (INSERM) UMR_S 1256.

Genotype-phenotype associations in French patients with phenylketonuria and importance of genotype for full assessment of tetrahydrobiopterin responsiveness.

BACKGROUND: Mutations in Phenylalanine Hydroxylase (PAH) gene cause phenylketonuria. Sapropterin (BH4), the enzyme cofactor, is an important therapeutical strategy in phenylketonuria. However, PAH is a highly polymorphic gene and it is difficult to identify BH4-responsive genotypes. We seek here to improve prediction of BH4-responsiveness through comparison of genotypes, BH4-loading test, predictions of responsiveness according to the literature and types and locations of mutations. METHODS: A total of 364 French patients among which, 9 % had mild hyperphenylalaninemia, 17.7 % mild phenylketonuria and 73.1 % classical phenylketonuria, benefited from a 24-hour BH4-loading test and had the PAH gene sequenced and analyzed by Multiplex Ligation Probe Amplification. RESULTS: Overall, 31.6 % of patients were BH4-responsive. The number of different mutations found was 127, including 26 new mutations. The mutations c.434A > T, c.500A > T, c.529G > C, c.1045 T > G and c.1196 T > C were newly classified as being BH4-responsive. We identified 261 genotypes, among which 46 were newly recognized as being BH4-responsive. Even though patients carry 2 responsive alleles, BH4-responsiveness cannot be predicted with certainty unless they present mild hyperphenylalaninemia. BH4-responsiveness cannot be predicted in patients carrying one responsive mutation only. In general, the milder the phenotype is, the stronger the BH4-response is. Almost exclusively missense mutations, particularly in exons 12, 11 and 8, are associated with BH4-responsiveness and any other type of mutation predicts a negative response. CONCLUSIONS: This study is the first of its kind, in a French population, to identify the phenotype associated with several combinations of PAH mutations. As others, it highlights the necessity of performing simultaneously BH4 loading test and molecular analysis in monitoring phenylketonuria patients.