Massive parallel sequencing as a new diagnostic approach for phenylketonuria and tetrahydrobiopterin-deficiency in Thailand
BMC Medical Genetics volume 18, Article number: 102 (2017)
Hyperphenylalaninemia (HPA) can be classified into phenylketonuria (PKU) which is caused by mutations in the phenylalanine hydroxylase (PAH) gene, and BH4 deficiency caused by alterations in genes involved in tetrahydrobiopterin (BH4) biosynthesis pathway. Dietary restriction of phenylalanine is considered to be the main treatment of PKU to prevent irreversible intellectual disability. However, the same dietary intervention in BH4 deficiency patients is not as effective, as BH4 is also a cofactor in many neurotransmitter syntheses.
We utilized next generation sequencing (NGS) technique to investigate four unrelated Thai patients with hyperphenylalaninemia.
We successfully identified all eight mutant alleles in PKU or BH4-deficiency associated genes including three novel mutations, one in PAH and two in PTS, thus giving a definite diagnosis to these patients. Appropriate management can then be provided.
This study identified three novel mutations in either the PAH or PTS gene and supported the use of NGS as an alternative molecular genetic approach for definite diagnosis of hyperphenylalaninemia, thus leading to proper management of these patients in Thailand.
Phenylketonuria (PKU) is an autosomal recessive metabolic disorder, characterized by progressive intellectual disability, mousy odor, autism, seizures, eczema and motor deficits . The incidence of PKU in Caucasians is approximately 1:10,000 births . Early diagnosis of the disease is beneficial since intellectual disability can be prevented by certain therapeutic intervention . Generally, measuring blood phenylalanine levels provides a reliable and practical diagnosis for PKU. If the blood levels of phenylalanine exceed 120 μmol/l (2 mg/dl), the individual is considered to be hyperphenylalaninemia (HPA) and needs further diagnosis . The majority of Caucasian patients with excessive phenylalanine levels are classified as classic PKU (MIM#261600), which is caused by loss-of-function mutations in the phenylalanine-4-hydroxylase (PAH) gene. About 2% of Caucasian cases are BH4-deficient , caused by mutations in the genes involved in BH4 biosynthesis consisting of 6-pyruvoyl-tetrahydropterin synthase (PTS), GTP cyclohydrolase I (GCH1), quinoiddihydropteridine reductase (QDPR), pterin-4-alpha-carbinolamine dehydratase 1 (PCBD1), and sepiapterin reductase (SPR) [5,6,7]. Interestingly, the incidence of PKU in Thailand is much less with a report of 1:212,535 in newborns  and the proportion of Thai patients with HPA who have BH4 deficiency remains unknown.
Although both PKU and BH4-deficiency patients show high blood phenylalanine levels and progressive intellectual disability, they respond to the treatment differently. Patients with PKU need additional supplement and are required to follow a diet that limits food with high phenylalanine such as dairy products, eggs, meat and fish . This diet therapy is highly recommended to continue throughout their life to prevent behavioral disorders, cognitive, and emotional dysfunction . On the other hand, patients with BH4 deficiency require BH4 supplement, generally 2–20 mg/kg/day [11, 12], as BH4 is also a cofactor in neurotransmitter synthesis, such as dopamine, serotonin, norepinephrine and epinephrine. Thus, deficiency of BH4 leads to not only elevated phenylalanine levels, but also deficiencies of these neurotransmitters . Therefore, it is important for physicians to give an early and precise diagnosis, whether HPA patients are PKU or BH4-deficiency, to provide a proper treatment to the patients.
Neonatal screening for PKU in Thailand was started in 1996 by Department of Medical Sciences, Ministry of Public Health, aiming to early detect and provide treatment to PKU patients through a screening program . The newborns’ blood samples were collected by heel prick on filter paper which was subsequently tested using the Guthrie method and/or the fluorometric method. The positive cases underwent plasma amino acid analysis using high performance liquid chromatography to determine phenylalanine levels . Determination of urinary pterin levels can be used to screen BH4-deficiency. Unfortunately, it is not an integral part of newborn screening in Thailand, which may lead to late diagnosis in some cases .
Conventional genetic diagnosis for PKU and BH4 deficiency relies on Sanger sequencing of the entire coding regions of the responsible genes. A total of 42 exons of PKU, BH4-deficiency associated genes are needed to be analyzed. This approach is time consuming and labor intensive . In contrast to Sanger sequencing, next-generation sequencing (NGS) can provide high throughput information from massive parallel sequencing from a single procedure  which is faster and more efficient. This technology has been widely used for research purpose but increasingly applied for genetic diagnosis, including PKU and BH4-deficiency [14, 16,17,18,19].
In this study, we apply next generation sequencing technique to investigate four unrelated Thai patients with hyperphenylalaninemia. The disease-causing variants were identified in all patients.
Patients I, II, III and IV are children from four unrelated non-consanguineous families. Their neonatal phenylalanine screening levels were 21.4, 53.11, 37.9, and 42.6 mg/dl, respectively (normal <20 mg/dl). After the diagnosis of hyperphenylalaninemia, all patients were treated with low-phenylalanine diet. For those diagnosed with BH4 deficiency, BH4 supplement was given. More clinical details are in Table 1.
Written informed consent for genetic analysis was obtained from the parents of the patients.
Three milliliters of peripheral blood were taken from each patient after informed consent. Genomic DNA was extracted by using Gentra Puregene Blood Kit, Qiagen (Qiagen, Hilden, Germany). The extraction process was done according to the manufacturer’s protocol. Genomic DNA from four patients were sent for exome sequencing using service from Macrogen, Inc. (Seuol, South Korea). The samples were prepared according to Agilent SureSelect Target Enrichment Kit (Agilent Technologies, Santa Clara, CA) preparation guide. The captured libraries were sequenced with Illumina HiSeq 2000 or HiSeq2500 Sequencer. The result sequences were aligned to the human genome reference sequence (UCSC hg19) using Burrows-Wheeler Alignment (bwa-0.7.10, http://bio-bwa.sourceforge.net/) . Picard software (picard-tools-1.118, http://broadinstitute.github.io/picard/) was used for marking and removing duplicated sequences. Genome Analysis Toolkit (GATK3.v4, https://www.broadinstitute.org/gatk/)  was used for data quality assessment, genotyping and variant calling. Finally, SnpEff_v.4.1 (http://snpeff.sourceforge.net/)  was used for variant annotation.
All homozygous and compound heterozygous variants with <1% allele frequency in global population and located on genes associated with hyperphenylalaninemia (PAH, PTS, GCH1, QDPR, PCBD1, SPR and GCHFR) were included. The variants with >1% of the 165 unrelated Thai exomes, our in-house database, were excluded. The novel variants were confirmed by Sanger sequencing. The patients’ genomic DNA was amplified by polymerase chain reaction (PCR), using primers specific to the site of the mutation (Table 2). The PCR products were then sent for Sanger sequencing by Macrogen Inc. (Seoul, South Korea). The sequences were analyzed by Gene Codes Sequencer software (v.5.4.1) (Gene Codes Corporation, MI). The sequencing data were compared to unaffected control sequences. ClustalX 2.1 (http://www.clustal.org/clustal2/)  was used for amino acid conservative analysis. ExAC Browser: Exome Aggregation Consortium Database (http://exac.broadinstitute.org/) was used to confirm the novelty of the variants .
All four patients showed approximately 20,000 DNA variants across the exome (Table 3). Mutation analysis successfully identified all eight mutant alleles in PKU or BH4-deficiency associated genes. Patients I and IV had mutations in the PAH gene, while Patients II and III had mutations in the PTS gene. The detailed characteristics of the mutations are summarized in Table 1. Sanger sequencing of the novel mutations confirmed the presence of the mutations found by exome sequencing (Fig. 1a and b). The amino acid alignment from ClustalX revealed that the two novel missense mutations in the PTS gene identified in Patient III were highly conserved (Fig. 1c).
Newborn screening program for HPA is important since severe intellectual disability from PKU or BH4 deficiency is preventable with proper treatment. Although Thailand has established Neonatal Screening Program for PKU for many years, screening positive patients do not undergo a test for BH4-deficiency due to the unavailability of pterin analysis in Thailand . In addition, BH4 loading test cannot distinguish patients with BH4 deficiency and BH4-responsive PKU . The conventional Sanger sequencing has been the gold standard test for a molecular approach in the genetic diagnosis of inherited disorders. This direct approach is suitable for sequencing hot-spot point mutations or small genes . However, it is rather costly and time consuming for diagnosis of genetic disorders which involve large genes or multiple genes, including genes associated with PKU and BH4-deficiency. The costs associated with NGS are rapidly decreasing. Therefore, we utilize NGS as a diagnostic tool in our study as it has been demonstrated to be effective in several previous reports [14, 16,17,18,19].
All eight mutant alleles were successfully identified. Of these, three have never been previously reported. Patient I has a novel homozygous frameshift mutation (c.686_687insA) in the PAH gene. Patient II has two compound heterozygous missense mutations (c.155A > G and c.200C > T) in the PTS gene. The c.155A > G mutation has been previously reported in Chinese patients  and the c.200C > T in Italian patients . Patient III has two novel compound heterozygous missense mutations (c.274A > C and c.326A > G) in the PTS gene. Both mutations are predicted to be damaging and deleterious from prediction softwares, SIFT and PROVEAN. Patient IV has compound heterozygous, one exonic splice site enhancer (c.611A > G) and one splice acceptor (c.1200-1G > C) mutations in the PAH gene. The c.611A > G mutation has been reported in Asian  and European patients  and the c.1200-1G > C has been recently reported in Chinese patients . Therefore, Patients I and IV were diagnosed as classic PKU and were continued on low phenylalanine diet, while Patients II and III were diagnosed with BH4-deficiency and were given BH4 supplement.
In order to reduce the number of identified variants from NGS data, investigators often need to apply filtering criteria to exclude low quality variants . This parameter could be depth of coverage. Interestingly, one of the variants in this study at position chr11:112,104,166, c.326A > G in the PTS gene has total read depth of just 19 reads, with an alternate read of only 5 reads. This can be considered as very low read depth in other studies, and may have been removed from analysis [32, 33]. We would not be able to find causative mutations in patient III if we discard this variant based on the variant quality alone. In addition, despite no consanguinity history in all four families, Patient I has a frameshift homozygous mutation in the PAH gene. From additional interview, it appeared that both of her parents came from the same district of the same province. It might not be uncommon for patients with non-consanguineous parents to have homozygous mutations.
Patients II and III received BH4 therapy after the definite diagnosis was made. Patient II however had delayed development and recurrent seizure treated with anticonvulsive drugs, phenytoin and levetiracetam. At his last visit when he was three years old, he still had delayed development and occasional abnormal movement (orofacial dyskinesia and chorea). This might be due to the age of Patient II when diagnosed, which was one year and nine months old. On the other hand, Patient III was given BH4 supplement at the age of 6 months after the diagnosis was made. The patient showed normal development at her last visit when she was one year and four months old. This demonstrated the importance of providing early and precise diagnosis to HPA patients.
Only 2% of Caucasian patients with HPA are BH4 deficient. Interestingly, of our four consecutive patients, two are BH4 deficient. The total number of patients are very small. However, since the treatment to prevent intellectual deficit is different between both types, it is warrant to determine the proportion of PAH and BH4 deficiency among Thai patients with HPA. If a significant percentage of Thai patients with HPA have BH4 deficiency, it should be justified to perform genetic testing using NGS as soon as HPA is identified.
This report represents the first study in Thailand to successfully use next-generation sequencing to detect causative mutations in PKU and BH4-deficiency cases, which allow physicians to provide a precise diagnosis and proper effective treatment to the patients. This approach is ready to substitute conventional sequencing for genetic diagnosis of patients with hyperphenylalaninemia.
next generation sequencing
Kawashima H, Kawano M, Masaki A, Sato T. Three cases of untreated classical PKU: a report on cataracts and brain calcification. Am J Med Genet. 1988;29(1):89–93.
Matalon R, Michals K. Phenylketonuria: screening, treatment and maternal PKU. Clin Biochem. 1991;24(4):337–42.
Al Hafid N, Christodoulou J. Phenylketonuria: a review of current and future treatments. Transl Pediatr. 2015;4(4):304–17.
Waisbren SE, Noel K, Fahrbach K, Cella C, Frame D, Dorenbaum A, Levy H. Phenylalanine blood levels and clinical outcomes in phenylketonuria: a systematic literature review and meta-analysis. Mol Genet Metab. 2007;92(1–2):63–70.
Blau N, Hennermann JB, Langenbeck U, Lichter-Konecki U. Diagnosis, classification, and genetics of phenylketonuria and tetrahydrobiopterin (BH4) deficiencies. Mol Genet Metab. 2011;104(Suppl):S2–9.
Thony B, Blau N. Mutations in the BH4-metabolizing genes GTP cyclohydrolase I, 6-pyruvoyl-tetrahydropterin synthase, sepiapterin reductase, carbinolamine-4a-dehydratase, and dihydropteridine reductase. Hum Mutat. 2006;27(9):870–8.
Clot F, Grabli D, Cazeneuve C, Roze E, Castelnau P, Chabrol B, Landrieu P, Nguyen K, Ponsot G, Abada M, et al. Exhaustive analysis of BH4 and dopamine biosynthesis genes in patients with Dopa-responsive dystonia. Brain. 2009;132(Pt 7):1753–63.
Pangkanon S, Ratrisawadi V, Charoensiriwatana W, Techasena W, Boonpuan K, Srisomsap C, Svasti J. Phenylketonuria detected by the neonatal screening program in Thailand. Southeast Asian J Trop Med Public Health. 2003;34(Suppl 3):179–81.
Acosta PB. Nutrition management of patients with inherited metabolic disorders. Sudbury, Mass: Jones and Bartlett Publishers; 2010.
Stemerdink BA, Kalverboer AF, van der Meere JJ, van der Molen MW, Huisman J, de Jong LW, Slijper FM, Verkerk PH, van Spronsen FJ. Behaviour and school achievement in patients with early and continuously treated phenylketonuria. J Inherit Metab Dis. 2000;23(6):548–62.
Shintaku H. Disorders of tetrahydrobiopterin metabolism and their treatment. Curr Drug Metab. 2002;3(2):123–31.
Blau N. Sapropterin dihydrochloride for the treatment of hyperphenylalaninemias. Expert Opin Drug Metab Toxicol. 2013;9(9):1207–18.
Pangkanon S, Charoensiriwatana W, Janejai N, Boonwanich W, Chaisomchit S. Detection of phenylketonuria by the newborn screening program in Thailand. Southeast Asian J Trop Med Public Health. 2009;40(3):525–9.
Trujillano D, Perez B, Gonzalez J, Tornador C, Navarrete R, Escaramis G, Ossowski S, Armengol L, Cornejo V, Desviat LR, et al. Accurate molecular diagnosis of phenylketonuria and tetrahydrobiopterin-deficient hyperphenylalaninemias using high-throughput targeted sequencing. Eur J Hum Genet. 2014;22(4):528–34.
Majewski J, Schwartzentruber J, Lalonde E, Montpetit A, Jabado N. What can exome sequencing do for you? J Med Genet. 2011;48(9):580–9.
Gu Y, Lu K, Yang G, Cen Z, Yu L, Lin L, Hao J, Yang Z, Peng J, Cui S, et al. Mutation spectrum of six genes in Chinese phenylketonuria patients obtained through next-generation sequencing. PLoS One. 2014;9(4):e94100.
Tao J, Li N, Jia H, Liu Z, Li X, Song J, Deng Y, Jin X, Zhu J. Correlation between genotype and the tetrahydrobiopterin-responsive phenotype in Chinese patients with phenylketonuria. Pediatr Res. 2015;78(6):691–9.
Li N, Jia H, Liu Z, Tao J, Chen S, Li X, Deng Y, Jin X, Song J, Zhang L, et al. Molecular characterisation of phenylketonuria in a Chinese mainland population using next-generation sequencing. Sci Rep. 2015;5:15769.
Chen YF, Jia HT, Chen ZH, Song JP, Liang Y, Pei JJ, Wu ZJ, Wang J, Qiu YL, Liu G, et al. Mutational spectrum of phenylketonuria in Jiangsu province. Eur J Pediatr. 2015;174(10):1333–8.
Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25(14):1754–60.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.
Cingolani P, Platts A, Wang le L, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila Melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 2012;6(2):80–92.
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8.
Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, O'Donnell-Luria AH, Ware JS, Hill AJ, Cummings BB, et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature. 2016;536(7616):285–91.
Vatanavicharn N, Kuptanon C, Liammongkolkul S, Liu TT, Hsiao KJ, Ratanarak P, Blau N, Wasant P. Novel mutation affecting the pterin-binding site of PTS gene and review of PTS mutations in Thai patients with 6-pyruvoyltetrahydropterin synthase deficiency. J Inherit Metab Dis. 2009;32(Suppl 1):S279–82.
Katsanis SH, Katsanis N. Molecular genetic testing and the future of clinical genomics. Nat Rev Genet. 2013;14(6):415–26.
Liu TT, Hsiao KJ. Identification of a common 6-pyruvoyl-tetrahydropterin synthase mutation at codon 87 in Chinese phenylketonuria caused by tetrahydrobiopterin synthesis deficiency. Hum Genet. 1996;98(3):313–6.
Oppliger T, Thony B, Kluge C, Matasovic A, Heizmann CW, Ponzone A, Spada M, Blau N. Identification of mutations causing 6-pyruvoyl-tetrahydropterin synthase deficiency in four Italian families. Hum Mutat. 1997;10(1):25–35.
Wang T, Okano Y, Eisensmith RC, Lo WH, Huang SZ, Zeng YT, Yuan LF, Liu SR, Woo SL. Missense mutations prevalent in Orientals with phenylketonuria: molecular characterization and clinical implications. Genomics. 1991;10(2):449–56.
Ellingsen S, Knappskog PM, Eiken HG. Phenylketonuria splice mutation (EXON6nt-96A-->g) masquerading as missense mutation (Y204C). Hum Mutat. 1997;9(1):88–90.
Rehm HL, Bale SJ, Bayrak-Toydemir P, Berg JS, Brown KK, Deignan JL, Friez MJ, Funke BH, Hegde MR, Lyon E, et al. ACMG clinical laboratory standards for next-generation sequencing. Genet Med. 2013;15(9):733–47.
Lim BC, Lee S, Shin JY, Kim JI, Hwang H, Kim KJ, Hwang YS, Seo JS, Chae JH. Genetic diagnosis of Duchenne and Becker muscular dystrophy using next-generation sequencing technology: comprehensive mutational search in a single platform. J Med Genet. 2011;48(11):731–6.
Broman M, Kleinschnitz I, Bach JE, Rost S, Islander G, Muller CR. Next-generation DNA sequencing of a Swedish malignant hyperthermia cohort. Clin Genet. 2015;88(4):381–5.
We would like to thank the patients and family members for participating in this study.
This research is supported by RachadapisekSompote Fund for Postdoctoral Fellowship, Chulalongkorn Academic Advancement into Its second Century Project, Chulalongkorn University, and Thailand Research Fund (BRG5980001).
Availability of data and materials
The data sets supporting the results of this article are available in the NCBI SRA database, which can be accessed via NCBI BioProject repository (accession number: PRJNA401212, https://www.ncbi.nlm.nih.gov/bioproject/401212).
BioSample accession numbers: SAMN07599868, SAMN07599869, SAMN07599870, SAMN07599871.
Ethics approval and consent to participate
Ethical approval was obtained from the institutional review board, Faculty of Medicine, Chulalongkorn University. The written informed consent for genetic analysis was obtained from the parents of the patients.
Consent for publication
The parents gave permission for the publication of the patient’s clinical details and images.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chaiyasap, P., Ittiwut, C., Srichomthong, C. et al. Massive parallel sequencing as a new diagnostic approach for phenylketonuria and tetrahydrobiopterin-deficiency in Thailand. BMC Med Genet 18, 102 (2017). https://doi.org/10.1186/s12881-017-0464-x
- Next generation sequencing
- Tetrahydrobiopterin deficiency
- Newborn screening