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.

Exome sequencing of cases with neural tube defects identifies candidate genes involved in one-carbon/vitamin B12 metabolisms and Sonic Hedgehog pathway.

Neural tube defects (NTD) result from complex mechanisms between genes, nutrition and environment. The identification of genetic predictors by genome exome sequencing and their influence on genome methylation need further consideration. Gene variants related to 1-CM metabolism (1-CM) could influence the methylation of genes involved in neural tube embryogenesis through impaired synthesis of S-adenosyl methionine. We performed exome sequencing of 6116 genes referenced in OMIM and NTD risk and genome-wide methylation in 23 NTD cases. We replicated the most significant associations in 81 other cases. The analysis of exome sequencing identified one gene of 1-CM, LRP2, and one gene of Sonic Hedgehog (SHH), GLI3, in the 23 NTD cases. The analysis restricted to genes of 1-CM and neural tube embryogenesis identified five gene predictors of 1-CM (LRP2, rs137983840; MMAA, rs148142853; TCN2, rs35838082; FPGS, rs41306702; BHMT, rs763726268) and two of SHH (GLI3, rs35364414; MKS1, rs151023718). We replicated the association of TCN2, BHMT and GLI3 with NTD risk in the 81 cases. We found a significant hemimethylation of CFAP46 that may influence SHH activation in one case, who carried risk alleles in BHMT, LRP2, MMAA and GLI3. In conclusion, we identified new candidate genes and rare variants that highlight an interacting influence of genes involved in SHH and 1-CM in the puzzle of genetic components of NTD risk.

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.