Correlation between genotype and the tetrahydrobiopterin-responsive phenotype in Chinese patients with phenylketonuria.

BACKGROUND: A growing body of research has suggested that tetrahydrobiopterin (BH4) responsive phenotype can be predicted by the phenylalanine hydroxylase (PAH) genotype in patients with phenylketonuria (PKU), but data concerning the association between genotype and BH4 responsiveness are scarce in China.
METHODS: A total of 165 PKU patients from China who had undergone a 24-h loading test with BH4 administration were recruited. Genotyping was performed by the next-generation sequencing (NGS) technique. Using the predicted residual PAH activity, we analyzed the association between genotype and BH4-responsiveness.
RESULTS: Among the 165 patients, 40 patients (24.24%) responded to BH4. A total of 74 distinct mutations were observed, including 13 novel mutations. The mutation p.R241C was most frequently associated with response. Two known mutations (p.A322T and p.Q419R) and two novel mutations (p.L98V and IVS3-2A>T) were first reported as responsive to BH4. Residual PAH activity of at least 12.5% was needed for responsive genotypes.
CONCLUSION: Genotype-based predictions of BH4-responsiveness are only for selecting potential responders. Accordingly, it is necessary to test potential responders with a long-term BH4 challenge.

Phenylketonuria (PKU, MIM261600) is an autosomal recessive disorder of amino acid metabolism and is caused by mutations in the phenylalanine hydroxylase (PAH, EC 1.14.16.1) gene. A mutant PAH gene gives rise to the decreased ability of the PAH enzyme to convert phenylalanine (Phe) to tyrosine (Tyr), which subsequently leads to the accumulation of Phe and its secondary metabolites in the body. A variety of clinical problems (e.g., intellectual impairment, autistic behavior, and seizures) manifest in patients with PKU unless treatment starts soon enough after birth (1). PKU falls into the most common inborn error of amino acid metabolism in China, with an incidence of 1:10,000 (2). So far, early identification followed by effective treatment for PKU patients play key roles in preventing the devastating consequences of PAH deficiency.

BH4 is the natural cofactor of PAH, and its pharmaceutical formulation (Sapropterin dihydrochloride) is an efficacious treatment with no severe adverse events for a subset of PKU patients (3). BH4 administration focuses not only on stabilizing the blood Phe level within the therapeutic range but also on improving dietary Phe tolerance by enhancing the catabolism of blood Phe (4,5,6). Not all patients are responsive to BH4 treatment, and most responders have mild biochemical phenotypes (7,8). However, there is no consensus on the gold standard test for ascertaining BH4 responsiveness in PKU patients.

Some evidence has indicated that the BH4-responsive genotype is characterized by substantial residual PAH activity (9,10,11). Moreover, there is a significant difference of PAH activity between BH4 responders and nonresponders (12). These may be related to the molecular mechanism behind BH4 response, including stabilizing the altered conformation of the PAH enzyme and preventing its degradation (13,14). Until now, more than 850 mutations have been reported in the BIOPKU database (University of Zurich, Zurich, Switzerland) (15). Therefore, mutation screening is of value to select potential responders and to exclude nonresponders before BH4 trial.

Due to the allelic heterogeneity in PKU patients from different geographical areas (16), it is worthwhile to detect mutant PAH genotypes for BH4 responsiveness in the Chinese population. Next-generation sequencing (NGS) has emerged as an efficient method for disease diagnosis, with the application of complete genome resequencing, transcriptome sequencing, small RNA sequencing, sequencing of bisulfite-treated DNA and so on (17). Currently, NGS has been performed for the genetic diagnosis of PKU, with increased diagnostic sensitivity and shorter analysis time and lower costs (18,19). This study is using NGS to analyze the PAH gene in a cohort of Chinese PKU patients, with the aim of assessing the association between genotype and BH4 responsiveness.

Results

Patients were classified into three biochemical phenotypes on the basis of the pretreatment blood phenylalanine concentration at the first definite diagnosis: 100 patients (60.61%) with classical PKU (cPKU, Phe levels >20 mg/dl), 49 patients (29.70%) with mild PKU (mPKU, Phe levels 10–20 mg/dl), and 16 patients (9.70%) with mild hyperphenylalaninemia (MHP, Phe levels < 10 mg/dl). Among the 40 patients (24.24%) responsive to BH4, 13 patients were characterized with mild hyperphenylalaninemia, 23 patients with mPKU, and 4 patients with cPKU.

PAH Mutation Spectrum

Mutations were identified on both alleles for 153 patients and on only a single allele for 12 patients, corresponding to a diagnostic efficiency of 96.36% (318/330). Among the 153 patients, 6 patients were carrying more than two mutations. A total of 75 distinct mutations were discovered (): 48 missense mutations, 11 splicing mutations, 9 nonsense mutations, 6 small deletion mutations, and 1 insertion mutation. These mutations were scattered in all exons and introns, except in introns 1, 2, 9, and 11. Over 65% of the mutations were localized in exons 7, 12, 6, and 11.

Eight mutations (p.R243Q, p.EX6-96A>G, p.R413P, p.R111X, p.Y356X, p.V399V, p.R241C, and IVS4-1G>A), each with a frequency ≥3%, accounted for 59.39% of the alleles. The most prevalent mutation was p.R243Q, with a relative frequency of 18.48%. The frequency for the other seven mutations accounted for 8.18, 7.88, 5.76, 5.15, 4.85, 4.55, and 4.55%.

A total of 13 novel mutations (), each occurring only once, have not been previously reported or listed in the PAHdb or BIOPKUdb. These novel mutations were identified in relevant parents. Meanwhile, these mutations were not present in the control group of 200 healthy individuals of Chinese origin and the >5,000 individuals sequenced and documented in the Exome Variant Server. Details on the pathologic analysis of five missense mutations are provided in .

Correlation Between Genotype and BH4 Responsiveness

In patients with a genotype consisting of two mutations. A total of 147 patients carrying two mutations constituted a spectrum of 105 different genotypic combinations (). The genotypes were divided into homozygous (n = 15) and heterozygous (n = 132). The heterozygous genotypes were further subdivided into three groups: null+null (n = 25), missense+null/functionally heterozygous (n = 74), and missense+missense/compound heterozygous (n = 33).

By analyzing the homozygous genotypes or functionally heterozygous genotypes, eight mutations were associated with responsiveness, including six known mutations (p.R241C, p.A322T, p.F331S, p.S349A, p.R408Q, and p.Q419R) and two novel mutations (p.L98V and IVS3-2A>T). However, four mutations (p.R53H, p.I65T, p.R158Q, and p. R261Q) previously defined as BH4 responsive were negative to BH4.

Three kinds of genotypes (p.[R111*];[EX6-96A>G], p.[R413P]; [Y356*], p.[A434D];[EX6-96A>G]) showed disparate BH4-responsiveness. A good genotype–phenotype correlation was observed in the remaining genotypic types (n = 102), which were unambiguously either responsive (n = 29, 28.43%) or unresponsive (n = 73, 71.57%) to BH4. Among these genotypes with the known predicted residual activity (PRA) in both mutations, the lowest residual PAH activity of the responsive genotype was at least 12.5%, and the highest residual PAH activity of the nonresponsive genotype was up to 46.5%.

As for the data on the BH4-responsiveness, only 6 of the 105 genotypes existed in the BIOPKU database. However, discordant BH4-responsiveness was observed in the two genotypes (p.[R53H];[R243Q] and p.[I65T];[R243Q]).

Patients with a genotype consisting of more than two mutations. Six patients with three or four mutations were observed (). Parents studies demonstrated that one missense mutation (p.R53H) and one splicing mutation (IVS7+2T>A) were always found on a common chromosome. However, these patients having the genotype with at least 26.33% residual PAH activity did not respond to BH4.

Discussion

Sapropterin dihydrochloride (BH4 pharmaceutical formulation) has been approved by the Food and Drug Administration in the United States and by the European Medicines Agency in Europe to treat a subset of PKU patients. Molecular analysis of the PAH gene would assist in predicting which patients could benefit from the BH4 therapeutic approach, which can alleviate the diet burden (20).Thus, the investigation of mutant genotypes for BH4 responsiveness is of great interest for PKU patients in China.

Responsive Mutations

Among the eight responsive mutations, four mutations (p.R241C, p.F331S, p.S349A, and p.R408Q) have previously been defined as BH4-responsive mutations. Similar to our findings on the BH4 response of p.R241C, Shintaku et al. (21) observed regularity in p.R241C regardless of the second mutation in the genotype. The second highest frequency responsive mutation, p.R408Q, was less consistently responsive. One out of seven patients was a nonresponder, having the genotype p.[R408Q];[R413P].

Two known mutations (p.A322T and p.Q419R) and two novel mutations (p.L98V, IVS3-2A>T) were first reported as responsive. Three missense mutations (p.L98V, p.A322T and p.Q419R) in the functionally heterozygous genotypes resulted in mild hyperphenylalaninemia/mPKU and responded to BH4. The splicing mutation IVS3-2A>T associated with cPKU was responsive to BH4, which may be affected by dietary fluctuations during BH4 loading test, resulting in incorrect assignment. However, some splicing mutations (e.g., IVS4−5C>G and IVS10 −3 C>T) were repeatedly observed to be responsive (11,22), so the BH4-responsive mechanism of these splicing mutations may be explained by generating a transcript encoding an enzyme with some activity. Therefore, further studies are needed to clarify the residual activities and BH4-responsive mechanisms of these splicing mutations.

However, four mutations (p.R53H, p.I65T, p.R158Q, and p.R261Q) previously defined as BH4 responsive were unresponsive to BH4. Three mutations, including p.I65T (23,24), p.R158Q (12,25), and p.R261Q (12,23), have been repeatedly reported with high rates of discordant BH4 responsiveness. Moreover, Trefz found mutations with a highly inconsistent rate of BH4-responsiveness in the regulatory domain (25), where p.R53H and p.I65T are located. The BH4-responsiveness of these mutations may be more easily affected by interallelic complementation in the genotype or other nongenetic factors (15). Therefore, they cannot be regarded as dominant in the estimation of BH4 responsiveness.

Genotype and Its Correlation With BH4 Responsiveness

Because most patients are compound heterozygotes, the complete genotype is a precondition for accurately predicting BH4 responsiveness. We found that patients with genotypes that lack residual “PAH activity” showed the highest degree of concordance with unresponsiveness, as previously reported (9). Residual PAH activity of at least 12.5% was required for responsive genotype, similar to the previous finding that a full genotype characterized by a residual activity of at least 15.5% was needed for response (11). However, the same genotypes with disparate BH4 responsiveness and genotypes with high residual PAH activities yet unresponsiveness were observed. Thus, the association between genotypes and BH4 responsiveness in our patients is uncertain, supporting the previous findings that genotype is not an absolute predictor of BH4 responsiveness (9,10,25). These discrepancies were mostly due to the difference of the residual PAH activity in vitro and in vivo. The residual PAH activity in vitro may be subject to the regulating function of gene variances in the introns of the PAH gene (26) and interallelic complementation (27). Accordingly, the residual PAH activity in vivo may be also affected by the genetic factors and the other regulating factors (e.g., BH4 pharmacokinetics, protein catabolic rate, and Phe intake during the test) in the interindividual variances (28,29,30).

As a general rule, the genotype consisted of two mutations in PKU patients. The reason for only one causative mutation found in 12 patients may be a large deletion or a mutation in the intronic sequences far from the exons or in the 5′ and 3′ untranslated regions, yet to be detected, at the PAH locus. However, six patients with atypical genotypes with three or even four mutations were observed. Such unusual situations were recently documented by other researchers (9,12,31), but few studies have investigated the BH4 responsiveness of genotypes with more than two mutations. Our findings agreed with other studies, that patients with the genotype carrying one or two defined responsive mutations and displaying a substantial residual activity did not respond to BH4. One patient was carrying two responsive mutations in cis (p.R241C and p.V245A) in combination with one mutation (p.R408W) (12), and another patient was harboring two mutations in cis (p.Y277C and p.V245A) in combination with one mutation (IVS12+1G>A) (31). In addition to interallelic complementation and many regulating elements, part of these anomalies may be explained by a specific mutation that is dominant in the genotype-based classification of BH4 responsiveness.

As stated above, the strength of our study is providing clinicians with more genotypic information on BH4 responsiveness before the BH4 loading test. Our results highlight the significance of screening the complete PAH genotype to increase the accuracy of genotype-based prediction of BH4 responsiveness. However, some limitations should be noted. First, the BH4 responsiveness of genotypes may be misclassified when only based on the results of the 24-h loading test. Some of the patients who benefited from BH4 therapy responded until 48 h after a BH4 trial (32), whereas initial responders were negative to BH4 after a long-term treatment (3). Second, the effect of a mutant genotype, especially for a genotype harboring one novel mutation or having more than two mutations, on protein structure and molecular function requires further study.

In conclusion, genotype-based predictions of BH4 responsiveness are only designed to identify patients who have a high probability of responding to BH4. Potential responders will be tested with a long-term BH4 challenge.

Methods

Patient Enrollment

From August 2012 to August 2013, we recruited 165 PKU patients (91 males and 74 females) from 21 newborn screening centers in China. All patients fulfilled the diagnostic criteria of PKU, with a blood phenylalanine concentration >2 mg/dl. BH4 deficiency was excluded by analysis of urinary pterins and dihydropteridine reductase activity in erythrocytes. Demographic data (basic information of the patient and family history of PKU) and clinical examination results (physical examinations and laboratory tests) were collected by a detailed and structured questionnaire. The patient ages ranged from 27 d to 13 y old (a median age of 3 y). These patients had undergone the same BH4 loading test/Phe+BH4 loading test in clinic. A single dose of 20 mg/kg exogenous BH4 (sapropterin dihydrochloride, BioMarin Pharmaceutical, CA) was used for the 24-h BH4 loading test performed in these patients, whereas 100 mg/kg Phe (Phenylalanine, Chinese and Western three-dimensional Pharmaceutical, Shanghai, China) was administered 3 h before the BH4 loading test in patients with <10 mg/dl blood Phe level. During the test, the patients maintained their usual diets without restricting Phe intake. Plasma Phe levels were taken at T = 0, 2, 4, 8, and 24 h. A patient is defined as a BH4-responder when there is a reduction of at least 30% in the blood phenylalanine concentration after the 24-h BH4 challenge (7).

This research was approved by the Ethics Committee of Sichuan University (No.2014018), and informed consent was obtained from the parents of the patients for the molecular study.

DNA Extraction

A blood sample (4 ml) was collected, in ethylene diamine tetraacetic acid, from the patients and their parents by venipuncture. Genomic DNA was extracted using a DNA extraction kit (Qiagen, Hilden, Germany).

Genotype Analysis With NGS

A series of unique primers (BoJin Biotechnology, Shanghai, China) was designed to amplify 13 exons and surrounding introns in patients, covering 200 bp upstream of the exon and 200 bp downstream of the exon of the PAH gene. A pair of primers was used to amplify approximately 150 bp DNA of the human β-globin (HBB) gene (accession number AY260740) as the internal quality control for identifying the false negatives caused by inadequate DNA or failed polymerase chain reaction (PCR).

To improve the throughput of the assay, the PCR primers not only amplify the target DNA but also provide a unique primer index for each of the 96 samples in each plate (multiple index PCR). The result is 96 sets of 10-bp nucleic acid tags at the 5′ end of each PAH and HBB primer (BoJin Biotechnology). The second index for each sample is the 8-bp nucleic acid tags in the library adapter sequence identifying each sample as coming from a specific 96-well plate. This index was attached to amplicons of the samples through the adapter preparation process. Using this “double index system,” hundreds of samples can be mixed together and detected in one sequencing chip at the same time.

PCR was conducted in a GeneAmp PCR system 9700 (Applied Biosystems, Foster City, CA), with a cycling protocol consisting of denaturation at 94 °C for 30 s, 56 °C for 30 s, and 72 °C for 1 min. After 35 cycles, gel electrophoresis was performed to check the quality of the amplified DNA. Only qualified DNA was allowed for the library preparation. PCR amplification products were prepared for DNA pooling, and diverse adapters were added to the amplification products. After the quantitative concentration satisfied the requirement of the library preparation, gene mutations were sequenced by an Illumina Hiseq 2000 sequencing instrument (Illumina, San Diego, CA).

After sequencing the samples, the raw sequence data were analyzed using in-house software. First, all sequence data were traced back to the specimens according to the sequences of primer indices and adapter indices. Second, the amplicon sequences of every sample were aligned with the standard reference for PAH sequences (NC_000012.11) in GenBank; SNPs were found in the target areas and relevant information was noted. The mutation nomenclature is in accordance with the guidelines from the Human Genome Variation Society (33).

Genotypic Validation

All mutations identified in patients were confirmed by analyzing parental DNA. When mutation loci were detected in patients, the same locus of the parental sample was amplified by PCR and analyzed by Sanger automated sequencing.

For newly detected mutations not present in the PAH database (McGill University, Montreal, Canada) (34) and in the BIOPKU database (University of Zurich, Zurich, Switzerland) (15), we randomly selected 200 healthy individuals and referred to the Exome Variant Server (35) to rule out the possibility of polymorphisms. Furthermore, the possible pathologic effect of missense mutations was predicted using PolyPhen-2, SIFT Blink, and SNPs3D (36).

Prediction of Responsiveness to BH4

The predicted residual PAH activity of one mutation in vitro refers to the BIOPKU database. Null mutations contained nonsense, frame-shift, and splicing mutations (only known mutations), and the missense mutations had less than 3% enzyme activity in vitro (36). The assigned BH4 responsiveness of one mutation is based on its effect in the homozygous form or functionally heterozygous form (one mutation in combination with a null mutation) (10).

Statement of Financial Support

This study was supported by the National “Twelfth Five-Year” Plan for Science & Technology Support of China (grant ID: 2012BAI09B04 and 2014BAI06B01). The authors have no other funding or conflicts of interest to disclose.

Table 1.

Spectrum of PAH mutations identified in 165 PKU patients

Table 2.

Pathologic analysis of five novel mutations of the PAH gene

Table 3.

Genotypes and association with BH4 responsiveness in 147 patients with two mutations

Table 4.

Genotypes with additional mutations in cis and BH4 responsiveness in six patients