Skip to main content

Whole exome sequencing identifies two novel variants in PHEX and DMP1 in Malaysian children with hypophosphatemic rickets



Hypophosphatemic rickets (HR) is a genetic disease of phosphate wasting that is characterized by defective bone mineralization. The most common cause of the disease is mutations in the phosphate regulating gene with homologies to endopeptidases on the X chromosome (PHEX) gene. The aims of this study were to identify the gene variants responsible for HR in three cases of Malaysian origin from three independent families and to describe their clinical, biochemical, and radiological features.


Whole exome sequencing (WES) was performed on all patients and their parents, followed by Sanger sequencing validation. Bioinformatics tools were used to provide supporting evidence for pathogenicity of variants. To confirm that a mutation is de novo, paternity test was carried out. High resolution melting curve analysis was performed to assess the allele frequency in normal controls for mutations that were found in the patients.


The patients showed typical characteristics of HR including lower limb deformity, hypophosphatemia, and elevated alkaline phosphatase. WES revealed two variants in the PHEX gene and one variant in the dentin matrix protein 1 (DMP1) gene. Two of the three variants were novel, including c.1946_1954del (p.Gly649_Arg651del) in PHEX and c.54 + 1G > A in DMP1. Our data suggests that the novel p.Gly649_Arg651del variant is likely pathogenic for HR disease.


This study extends the variant spectrum of the PHEX and DMP1 genes. Our findings indicate that WES is an advantageous approach for diagnosis of genetic diseases which are heterogeneous.


Hypophosphatemic rickets (HR) is a genetic disorder in which the proximal renal tubule cannot reabsorb sufficient phosphate (Pi), leading to growth retardation, rickets, and osteomalacia [1]. The most frequent form of hereditary HR is X-linked hypophosphatemic rickets (XLHR; MIM#307800), in which 80% of familial HR and 70% of sporadic HR are diagnosed as XLHR [2]. XLHR is caused by loss-of-function (LOF) mutations in the phosphate regulating gene with homologies to endopeptidases on the X chromosome (PHEX).

The XLHR inheritance pattern is X-linked dominant as both males and females present with the disease, and its prevalence is 1 in 20,000 live births [3]. The number of PHEX mutations reported in the Human Gene Mutation Database (HGMD) [4] is more than 400 mutations, most of which are missense and nonsense mutations, however other types including small deletions/insertions, abnormal splicing, and gross insertions/deletions have also been reported [1, 5].

Other less common forms of hereditary HR include autosomal dominant HR (ADHR), autosomal recessive HR 1 (ARHR1), autosomal recessive HR 2 (ARHR2), hereditary HR with hypercalciuria (HHRH), and X-linked recessive HR (XRHR) which are caused by gene mutations in fibroblast growth factor 23 (FGF23), dentin matrix acidic phosphoprotein 1 (DMP1), ectonucleotide pyrophosphatase/phosphodiesterase 1 (ENPP1), solute carrier family 34, member 3 (SLC34A3), and chloride channel 5 (CLCN5), respectively.

FGF23 plays an important role in regulating phosphate and vitamin D homeostasis. PHEX, DMP1 and ENPP1 negatively regulate the expression of FGF23, which is produced in osteocytes and secreted into the circulation [6, 7]. The deficiencies of these genes lead to excessive secretion of FGF23. Excess FGF23 causes hypophosphatemia in two ways. The expression of phosphate co-transporters in the kidney is suppressed by FGF23, leading to decreased phosphate reabsorption. FGF23 also reduces the production of 1,25-dihydroxyvitamin D (1,25(OH)2D3) and thus decreases the absorption of phosphate in the intestine. Other biochemical findings include elevated serum alkaline phosphatase (ALP) and normal to slightly high parathyroid hormone (PTH) levels [8, 9].

The clinical features are similar in various forms of hypophosphatemia, although this does not imply that the symptoms are identical. Children mostly present with rickets, short stature, dental abscesses, and lower limb bone deformities. Symptoms of the disease in adults include osteomalacia, bone pain, and enthesopathy [5]. Treatment for XLHR, ADHR or ARHR, includes oral phosphate and vitamin D supplementations while, patients with HHRH are only treated with phosphate supplementation [10]. This conventional treatment improves rickets in children. However, a significant number of patients may not show successful results [11] and complications such as nephrocalcinosis and hyperparathyroidism have been reported in some cases [12]. In addition to the conventional treatments, burosumab, which was approved by the US Food and Drug Administration in 2018, is used to treat XLHR. Burosumab is a human monoclonal antibody targeting FGF23, thus, significantly improves rickets, growth, and biochemical abnormalities [8, 12].

The rarity of HR may lead to delayed diagnosis. In addition, the recent developed treatment for this disease cannot be prescribed to the patients without defining and accurately diagnosing the disease. Therefore, correct and early diagnosis of this disease is important for optimal treatment. In the Malaysian population, there is very little information about the clinical, biochemical and, most importantly, genetic features of HR [13, 14]. Previously, only one genetic study of HR disease had been performed in Malaysia in which Sanger sequencing for PHEX, FGF23 and DMP1 genes in four Malaysian HR pediatric patients revealed two PHEX variants [14].

In the present study, we report clinical, radiological, laboratory, and genetic data of three pediatric patients with HR in Malaysia, for which whole exome sequencing (WES) revealed three candidate pathogenic variants. Further analysis was employed to support the pathogenicity of the variants.


Subjects and study approval

Pediatric patients with HR were identified from the pediatric endocrine clinic at Hospital Serdang and Hospital Kuala Lumpur, Malaysia in 2019. The inclusion criteria for the study were patients presenting with clinical and biochemical features of HR including lower limb bowing deformity, short stature, low serum phosphate, elevated alkaline phosphatase, normal serum calcium, and normal 25-hydroxyvitamin D. Informed consent for genetic investigation was obtained from the patients’ parents. Peripheral whole blood samples from three patient-parent trios were collected in lithium heparin-containing tubes and genomic DNA (gDNA) was extracted using the QIAamp DNA Blood Midi Kit (QIAGEN, Hilden, Germany). Approval for ethics was obtained from the Medical Research and Ethics Committee (MREC), Ministry of Health, Malaysia (reference number: NMRR-17-612-33,842 (IIR)).

Genetic analysis

The DNA samples of three affected children and their parents from unrelated families were sent to Macrogen (Seoul, South Korea) for WES. Sequencing libraries were generated using SureSelect V5-post kit (Agilent Tech. Ltd., USA) with 1 μg of gDNA. The size of the PCR-enriched fragments was verified by testing the template size distribution on the 2100 Bioanalyzer (Agilent Tech. Ltd., USA) using a DNA 1000 chip. Quantification of DNA libraries was done using a qPCR technique based on the Illumina qPCR Quantification Protocol Guide. The captured DNA was sequenced as paired end reads (150 bp) on an Illumina Novaseq6000 (Illumina, CA, USA) sequencer. The paired-end sequences were aligned to the human reference genome assembly GRCh37 (hg19) using the Burrows-Wheeler Aligner (BWA) software. Duplicate reads were removed using Picard software. Sequence variant calling was performed using Genome Analysis Tool kit (GATK) version GATK v4.0.5.1 and annotated with effect prediction using the SnpEff tool. Screening of study subjects initially involved examining the WES data for variants in the hereditary HR-related genes including PHEX (NM_000444.5), FGF23 (NM_020638.2), DMP1 (NM_004407.3), ENPP1 (NM_006208.2), SLC34A3 (NM_001177316.1), and CLCN5 (NM_001127898.3). The candidate variants identified through WES were validated by Sanger sequencing using the DNA sequencing service provided by Apical (Apical Scientific Sdn. Bhd., Malaysia). The details of PCR protocol are provided in Additional file 1 Supplementary Table 1. The segregation of variants with disease phenotype in each patient-parent trio was evaluated. Since the parents in the three trios were phenotypically normal, two strategies were used to filter the variants in each trio; de novo strategy to find the de novo variants in affected child, and double-hit strategy to identify homozygous or compound heterozygous variants in affected child, recessively inherited from unaffected parents [15]. Paternity testing was performed for the patient when the parents were not carriers of the variant and there was no family history of the disease, in order to confirm the de novo occurrence of a variant. To do so, Galaxy browser [16] was used to filter the WES data for 50 SNPs which the panel has been developed and validated by Yousefi et al. [17] for the unique identification of an individual. Then, 50 normal female controls of Malay ethnicity were screened for the novel variants using HRM analysis using a LightCycler 480 instrument (Roche, Switzerland) as outlined in Additional file 1 Supplementary Table 2.

Bioinformatics analysis of novel sequence variants

Novel variants were assessed for their allele frequency in the population databases including 1000 Genome Project [18], Genome Aggregation Database (gnomAD) [19], and Exome Sequencing Project (ESP) [20]. Moreover, the bioinformatic tools were used for novel variants to predict their impact and likelihood of causing disease. The effect of in-frame deletion was evaluated using MutationTaster [21] and MutPred-Indel [22]. To predict the impact of splicing variant, Human Splicing Finder (HSF) 3.1 [23], MaxEntScan [24], SpliceAI [25], and NNSPLICE 0.9 [26] were used. Conservation analysis of nucleotides was performed using three programs including PhyloP [27], GERP++ [28], and PhastCons [29] through the UCSC Genome Browser.


Patients’ characteristics

There were three patients of Malay ethnicity identified and recruited for the study. Detailed clinical, biochemical, radiological, and genetic findings are summarized in Table 1. Patient 1 had delayed presentation and diagnosis of HR at 6 years. However, history taking revealed that patient 1 had bowing of legs since 1–2 years of age but the parents did not seek medical attention then. The age of onset of skeletal abnormalities for all patients in this study ranged from 1 to 2 years. Radiographs of all patients showed evidence of rickets with fraying and splaying of the metaphyses of long bones of the limbs (Fig. 1). None of the patients had family members with HR or skeletal deformity. Patient 1 and patient 2 had non-consanguineous parents, while the parents of patient 3 were consanguineous and shared the same maternal grandfather.

Fig. 1
figure 1

Radiograph of lower limbs of patient 1. Radiograph shows features of rickets with fraying and splaying of the metaphyses of the distal femur and distal tibia (red arrows)

Table 1 Clinical, biochemical, radiological, and genetic findings of patients with HR

Variants identified through whole exome sequencing and sanger sequencing confirmation

Three parent-child trios from three unrelated families underwent WES to screen HR genes for mutations. After selecting the de novo variants and variants segregating as recessive traits, three candidate variants were identified which were absent from population databases including 1000 Genome Project, gnomAD, and ESP. A heterozygous variant in PHEX, c.871C > T (p.Arg291Ter), was found in patient 1 which has been reported in HGMD and ClinVar [30] and was absent in both of the patients parents. The other two variants which are novel and unreported, include the PHEX variant c.1946_1954del (p.Gly649_Arg651del) in patient 2 and the DMP1 variant c.54 + 1G > A in patient 3. Sanger sequencing confirmed the presence of the three variants in each patient (Fig. 2). The likelihood of pathogenicity of the novel variants was further assessed.

Fig. 2
figure 2

Sequence confirmation of variants and pedigrees of three families with HR. (A) Heterozygous c.871C > T variant in exon 8 of the PHEX gene in patient 1 and wildtype genotype in her parents. (B) Heterozygous c.1946_1954del variant in exon 19 of the PHEX gene in patient 2 and wildtype genotype in her parents. (C) Homozygous c.54 + 1G > A variant in intron 2 of the DMP1 gene in patient 3 and heterozygous genotype in her parents. Circles and squares in pedigrees represent females and males, respectively. Solid circles with arrow indicate index cases. Double lines between symbols indicate consanguineous marriage. Genotype of each individual is shown under the symbols

Pathogenicity analysis of novel variants

The pathogenicity of the variants was assigned according to guidelines of the American College of Medical Genetics and Genomics (ACMG) [31]. WES in this study detected one previously reported variant and two novel variants in Malay patients with HR disease. Patient 1 was heterozygous for the previously reported nonsense PHEX variant (c.871C > T; p.Arg291Ter) which is located in exon 8 of the gene and the extracellular domain of the PHEX protein (Fig. 3a). The variant was not inherited from her healthy parents (Fig. 2a) and was therefore considered a de novo variant with the confirmation of paternity. The p.Arg291Ter variant is reported in ClinVar as pathogenic with multiple submissions and the collective evidence supports that p.Arg291Ter in the PHEX gene is a pathogenic variant for HR disease (Additional file 1 Supplementary Table 3).

Fig. 3
figure 3

Location of PHEX and DMP1 variants identified in the present study. (A) The bar with numbers represents the 22 exons of the PHEX cDNA. Diagram of the PHEX protein shows the three domains including the N-terminal cytoplasmic domain, the transmembrane domain, and the extracellular domain with conserved cysteine residues (C) and the zinc-binding motifs (Zn). (B) Schematic representation of the DMP1 gene. The numbers represent the 6 exons of the DMP1 gene. The two novel PHEX and DMP1 variants are in bold

A novel de novo variant causing an in-frame deletion in exon 19 of the PHEX gene (c.1946_1954del; p.Gly649_Arg651del) was identified in patient 2 who was heterozygous for this variant. The variant was absent in both parents (Fig. 2b) and the paternity testing confirmed de novo occurrence of p.Gly649_Arg651del in this patient. Annotation of the deletion was absent in disease and population databases including Clinvar, HGMD, 1000 genome, gnomAD, and ESP. Furthermore, the variant was absent in 100 normal control chromosomes of Malay ethnicity (Additional file 1 Supplementary Fig. 1). The bioinformatic tools used for the prediction of in-frame deletion, did not predict a damaging effect on the PHEX protein. Nevertheless, the deleted positions are conserved across species (Additional file 1 Supplementary Table 3). In particular, the glycine residue 649 is highly conserved in the members of a family of zinc-metallopeptidases including the PHEX protein and is located in the zinc-binding motif of PHEX which is important for the catalytic activity [32] (Fig. 3a). Furthermore, missense variants at positions 649 and 651 of the PHEX protein have been reported in HR patients [1, 32, 33], which support that these regions of the protein are important for PHEX function. The evidence suggests p.Gly649_Arg651del can be considered a likely pathogenic variant for HR disease (Additional file 1 Supplementary Table 3).

A second novel variant identified was a canonical splice donor site, c.54 + 1G > A, in intron 2 of the DMP1 gene in patient 3 (Fig. 3b). This splice site variant was of particular interest because the segregation in trio 3 showed that this homozygous variant in the patient was transmitted from the consanguineous parents who were first cousins and carriers of c.54 + 1G > A (Fig. 2c). The disease and population databases have not reported this variant. Fifty ethnically matched (Malay) normal controls were screened for c.54 + 1G > A. No controls were homozygous for the variant and the heterozygous genotype was observed in 12 (24%) of the controls. (Additional file 1 Supplementary Fig. 1). DMP1 is a recessive gene and LOF mutations of the gene cause ARHR disease. Therefore, only the homozygous genotype in patient 3 causes the disease phenotype and the heterozygous carriers are not affected. The c.54 + 1G > A variant occurs at a nucleotide position that is highly conserved and multiple splicing prediction programs predict this variant has a deleterious effect due to the loss of wildtype donor site (Additional file 1 Supplementary Table 3). According to ACMG, canonical splice-site variants result in a null effect, however, functional analysis is required to confirm the impact of variant. Therefore, the c.54 + 1G > A variant in the DMP1 gene is interpreted as a variant with uncertain significance in relation to HR disease since functional study has not been performed and the criteria are not sufficient for assigning the pathogenicity (Additional file 1 Supplementary Table 3) [31].


HR is a rare genetic disease with genetic heterogeneity [34]. There has been scarce genetic confirmation of HR in Malaysia. The present study was established to determine if molecular diagnosis of three Malaysian HR patients could be achieved after screening the common genes for HR. Although clinical, radiographic, and laboratory data can help diagnose HR, genetic testing is more reliable for accurate diagnosis. All patients in this study showed typical clinical, radiographic and biochemical features of HR, with no family history of HR.

In the present study, two variants in PHEX and one variant in DMP1 were identified using WES in three unrelated patients with HR from Malaysia. The two PHEX variants were de novo in the patients who were both females. This reflects previous findings of a higher rate of the PHEX de novo mutations in females than in males [35]. This is likely due to X-linked mutations occurring more frequently on the X chromosome in paternal germ cells [34]. The presence of de novo variants is one of the criteria to consider when establishing the pathogenicity of a variant [31, 36]. The confirmation of paternity testing in the two families with de novo PHEX mutations and the index patients being affected, provide strong criteria towards classifying the pathogenicity of these variants.

The de novo nonsense mutation in PHEX (c.871C > T, p.Arg291Ter) identified in this study has frequently been reported in HR patients from several studies of other populations including European, Chinese, and Japanese which may suggest that the nucleotides of the Arg291 residue is a mutation hot spot in PHEX and the present study confirms that p.Arg291Ter variant is a cause of HR disease in the Malay population. The patients with the p.Arg291Ter mutation in the previous studies had a similar phenotype to patient 1 including hypophosphatemia, bone deformities and radiological signs of rickets which are typical phenotype of HR [1, 32, 37,38,39,40,41,42,43,44,45,46,47,48].

The novel in-frame deletion identified in this study is located in exon 19 of PHEX (p.Gly649_Arg651del). PHEX is a member of the M13 family of membrane-bound zinc-metallopeptidases. There is a sequence motif (642ENXADXGG649) in exon19 of PHEX that is highly conserved in the members of the M13 family. A previous study in neutral endopeptidase (NEP), another member of the M13 family, has shown that this conserved motif plays an important role in zinc binding and catalysis activity [49]. The glycine 649 residue is located in the 642ENXADXGG649 motif, therefore, the deletion of this residue may interfere with zinc binding and catalytic activity of PHEX. A study by Kinoshita et al. [33] also revealed a missense mutation at this position (p.Gly649ASP) in three affected patients with HR from one family. In addition, Francis et al. [32] and Zheng et al. [1] reported a missense mutation at Arg 651 of PHEX (Arg651Pro) in two unrelated patients with HR. Although the prediction tools did not predict this deletion as disease-causing, the previous reports may indicate the importance of glycine and arginine residues at positions 649 and 651 of PHEX, respectively, and that these positions may be mutation hot spots.

The novel splicing variant in the DMP1 gene (c.54 + 1G > A) was also detected with WES. Canonical eukaryotic splicing relies on four conserved nucleotides, i.e. the donor sequence GT at the 5′ end of the intron, and the acceptor sequence AG at the 3′ end [50]. The variant c.54 + 1G > A is located within the canonical site at the 5′ end of intron 2 of the DMP1 gene. Mutations at the canonical splice sequences usually lead to single exon skipping [51]. The prediction of aberrant splicing and the involvement of highly conserved splice site sequences, provided evidence that the variant is likely causing HR in the patient. Of note, the same nucleotide position in DMP1 but with different base change (c.54 + 1G > C) has been reported in a Chinese family with HR disease. Two patients of the family are homozygous and the healthy parents are carriers of the c.54 + 1G > C variant. Ni et al. showed that c.54 + 1G > C in the DMP1 gene leads to the skipping of the whole exon 2 which then disrupts the reading frame of DMP1 [52]. Furthermore, the c.54 + 1G > C variant has been reported in gnomAD in the heterozygous state with a minor allele frequency (MAF) < 0.1% (0.00039%).

Although this research presents novel findings from a unique ethnicity which is rarely studied, it also has several limitations. Small sample size is one of the limitations which is due to the rareness of HR disorder and small area of data collection. Therefore, statistical analysis regarding phenotype and genotype of patients could not be performed. Moreover, genetic analysis was carried out only for parent-child trios and multi-generational genetic investigation was not available. This is particularly important for the family in which the parents of an affected child were carriers of the variant. Another limitation is the absence of functional analysis to fully understand the impact of novel variants.


We have identified candidate disease-causing variants in three Malaysian HR cases of Malay ethnicity. Two of the variants were novel and our data support that the PHEX novel variant is likely pathogenic. However, the evidence was not sufficient for the novel DMP1 variant for assigning the pathogenicity and has been classified as a variant of uncertain significance. This study extends the variant spectrum in PHEX and DMP1 and has shown the advantages of WES as an approach for the diagnosis of genetic diseases that are genetically heterogeneous. However, using conventional sequencing approaches to analyze the common genes prior to WES may identify the causative mutations in a proportion of patients and reduce the cost of genetic analysis. The findings of this study can be useful in genetic counseling and management of patients. Functional analysis of the identified variants is needed to verify the pathogenicity of these variants and will help to understand the underlying mechanisms of the disease and development of targeted therapies.

Availability of data and materials

The datasets used and analyzed during the current study are available from the corresponding authors on reasonable request.



Hypophosphatemic rickets


phosphate regulating gene with homologies to endopeptidases on the X chromosome


dentin matrix acidic phosphoprotein 1


Whole exome sequencing


X-linked hypophosphatemic rickets




Human Gene Mutation Database


autosomal dominant HR


autosomal recessive HR


hereditary HR with hypercalciuria


X-linked recessive HR


fibroblast growth factor 23


ectonucleotide pyrophosphatase/phosphodiesterase 1


solute carrier family 34, member 3


chloride channel 5


1,25-dihydroxyvitamin D


alkaline phosphatase


parathyroid hormone


Medical Research and Ethics Committee


Genome Aggregation Database


Exome Sequencing Project


Human Splicing Finder.


  1. Zheng B, Wang C, Chen Q, Che R, Sha Y, Zhao F, et al. Functional characterization of PHEX gene variants in children with X-linked Hypophosphatemic rickets shows no evidence of genotype-phenotype correlation. J Bone Miner Res. 2020. 35, (9) 1718,

  2. Lin Y, Xu J, Li X, Sheng H, Su L, Wu M, et al. Novel variants and uncommon cases among southern Chinese children with X-linked hypophosphatemia. J Endocrinol Investig. 2020;43(11):1577.

    Article  CAS  Google Scholar 

  3. Bitzan M, Goodyer PR. Hypophosphatemic rickets. Pediatr Clin. 2019;66(1):179–207.

    Article  Google Scholar 

  4. Human Gene Mutation Database. Public database Accessed 11 October 2021.

  5. Tavana N, Thilakavathy K, Kennerson M, Ting TH. Genetic basis of hereditary hypophosphataemic rickets and phenotype presentation in children and adults. Endokrynol Pol. 2021;72(4):366–94.

    Article  CAS  Google Scholar 

  6. Liu S, Guo R, Simpson LG, Xiao Z-S, Burnham CE, Quarles LD. Regulation of fibroblastic growth factor 23 expression but not degradation by PHEX. J Biol Chem. 2003;278(39):37419–26.

    Article  CAS  Google Scholar 

  7. Huang X, Jiang Y, Xia W. FGF23 and phosphate wasting disorders. Bone Res. 2013;1(2):120–32.

    Article  CAS  Google Scholar 

  8. Zhang C, Zhao Z, Sun Y, Xu L, JiaJue R, Cui L, et al. Clinical and genetic analysis in a large Chinese cohort of patients with X-linked hypophosphatemia. Bone. 2019;121:212–20.

    Article  CAS  Google Scholar 

  9. Shimada T, Hasegawa H, Yamazaki Y, Muto T, Hino R, Takeuchi Y, et al. FGF-23 is a potent regulator of vitamin D metabolism and phosphate homeostasis. J Bone Miner Res. 2004;19(3):429–35.

    Article  CAS  Google Scholar 

  10. Acar S, BinEssa HA, Demir K, Al-Rijjal RA, Zou M, Çatli G, et al. Clinical and genetic characteristics of 15 families with hereditary hypophosphatemia: novel mutations in PHEX and SLC34A3. PLoS One. 2018;13(3):e0193388.

    Article  CAS  Google Scholar 

  11. Haffner D, Emma F, Eastwood DM, Duplan MB, Bacchetta J, Schnabel D, et al. Clinical practice recommendations for the diagnosis and management of X-linked hypophosphataemia. Nat Rev Nephrol. 2019;15(7):435–55.

    Article  Google Scholar 

  12. Imel EA, Glorieux FH, Whyte MP, Munns CF, Ward LM, Nilsson O, et al. Burosumab versus conventional therapy in children with X-linked hypophosphataemia: a randomised, active-controlled, open-label, phase 3 trial. Lancet. 2019;393(10189):2416–27.

    Article  CAS  Google Scholar 

  13. Yong S, Aik S. X-linked hypophosphatemic rickets--a report of 2 cases and review of literature. Med J Malaysia. 2000;55:101–4.

    Google Scholar 

  14. Razali NN, Tzer HT, King HL, Kennerson ML, Thilakavathy K. Predicting the impact of PHEX, FGF23 and DMP1 gene variants found in Malaysian Malay patients with Hypophosphataemic rickets through in silico analysis of protein function and mRNA secondary structure. J Biochem Microbiol Biotechnol. 2019;7(2):14–23.

    Article  Google Scholar 

  15. Gilissen C, Hoischen A, Brunner HG, Veltman JA. Disease gene identification strategies for exome sequencing. Eur J Hum Genet. 2012;20(5):490–7.

    Article  CAS  Google Scholar 

  16. Galaxy browser. .

  17. Yousefi S, Abbassi-Daloii T, Kraaijenbrink T, Vermaat M, Mei H, , et al. A SNP panel for identification of DNA and RNA specimens. BMC Genomics 2018;19. (1):90

    Article  CAS  Google Scholar 

  18. 1000 Genome Project. .

  19. Genome Aggregation Database. .

  20. Exome Sequencing Project. .

  21. Schwarz JM, Cooper DN, Schuelke M, Seelow D. MutationTaster2: mutation prediction for the deep-sequencing age. Nat Methods. 2014;11(4):361–2.

    Article  CAS  Google Scholar 

  22. MutPred-Indel. .

  23. Human Splicing Finder. .

  24. MaxEntScan. .

  25. SpliceAI. .

  26. Splice Site Prediction by Neural Network. .

  27. Pollard KS, Hubisz MJ, Rosenbloom KR, Siepel A. Detection of nonneutral substitution rates on mammalian phylogenies. Genome Res. 2010;20(1):110–21.

    Article  CAS  Google Scholar 

  28. Davydov EV, Goode DL, Sirota M, Cooper GM, Sidow A, Batzoglou S. Identifying a high fraction of the human genome to be under selective constraint using GERP++. PLoS Comput Biol. 2010;6(12):e1001025.

    Article  CAS  Google Scholar 

  29. Siepel A, Bejerano G, Pedersen JS, Hinrichs AS, Hou M, Rosenbloom K, et al. Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes. Genome Res. 2005;15(8):1034–50.

    Article  CAS  Google Scholar 

  30. ClinVar. .

  31. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17(5):405–23.

    Article  Google Scholar 

  32. Francis F, Strom TM, Hennig S, Böddrich A, Lorenz B, Brandau O, et al. Genomic organization of the human PEX gene mutated in X-linked dominant hypophosphatemic rickets. Genome Res. 1997;7(6):573–85.

    Article  CAS  Google Scholar 

  33. Kinoshita Y, Saito T, Shimizu Y, Hori M, Taguchi M, Igarashi T, et al. Mutational analysis of patients with FGF23-related hypophosphatemic rickets. Eur J Endocrinol. 2012;167(2):165.

    Article  CAS  Google Scholar 

  34. Yuan L, Wu S, Xu H, Xiao J, Yang Z, Xia H, et al. Identification of a novel PHEX mutation in a Chinese family with X-linked hypophosphatemic rickets using exome sequencing. Biol Chem. 2015;396(1):27–33.

    Article  CAS  Google Scholar 

  35. Durmaz E, Zou M, Al-Rijjal RA, Baitei EY, Hammami S, Bircan İ, et al. Novel and de novo PHEX mutations in patients with hypophosphatemic rickets. Bone. 2013;52(1):286–91.

    Article  CAS  Google Scholar 

  36. Zhang J, Yao Y, He H, Shen J. Clinical interpretation of sequence variants. Curr Protoc Hum Genet. 2020;106(1):e98.

    Google Scholar 

  37. Dixon PH, Christie PT, Wooding C, Trump D, Grieff M, Holm I, et al. Mutational analysis of PHEX gene in X-linked hypophosphatemia. J Clin Endocrinol Metab. 1998;83(10):3615–23.

    Article  CAS  Google Scholar 

  38. Sato K, Tajima T, Nakae J, Adachi M, Asakura Y, Tachibana K, et al. Three novel PHEX gene mutations in Japanese patients with X-linked hypophosphatemic rickets. Pediatr Res. 2000;48(4):536–40.

    Article  CAS  Google Scholar 

  39. Obara-Moszynska M, Rojek A, Kolesinska Z, Jurkiewicz D, Chrzanowska K, Niedziela M. X-linked hypophosphatemic rickets in children: clinical phenotype, therapeutic strategies and molecular background. Endokrynol Pol. 2020;72(2):108.

    Article  Google Scholar 

  40. Holm IA, Huang X, Kunkel LM. Mutational analysis of the PEX gene in patients with X-linked hypophosphatemic rickets. Am J Hum Genet. 1997;60(4):790.

    CAS  Google Scholar 

  41. Gaucher C, Walrant-Debray O, Nguyen T-M, Esterle L, Garabédian M, Jehan F. PHEX analysis in 118 pedigrees reveals new genetic clues in hypophosphatemic rickets. Hum Genet. 2009;125(4):401.

    Article  Google Scholar 

  42. Hernández-Frías O, Gil-Peña H, Pérez-Roldán JM, González-Sanchez S, Ariceta G, Chocrón S, et al. Risk of cardiovascular involvement in pediatric patients with X-linked hypophosphatemia. Pediatr Nephrol. 2019;34(6):1077–86.

    Article  Google Scholar 

  43. Zhang S, Zhang Q, Cheng L, Huang X, Peng Y, Liang Z, et al. Analysis of PHEX gene mutations in three pedigrees affected with hypophosphatemic rickets. Zhonghua Yi Xue Yi Chuan Xue Za Zhi (Chinese Journal of Medical Genetics). 2018;35(5):644–7.

    Article  Google Scholar 

  44. Chesher D, Oddy M, Darbar U, Sayal P, Casey A, Ryan A, et al. Outcome of adult patients with X-linked hypophosphatemia caused by PHEX gene mutations. J Inherit Metab Dis. 2018;41(5):865–76.

    Article  CAS  Google Scholar 

  45. Li S-S, Gu J-M, Yu W-J, He J-W, Fu W-Z, Zhang Z-L. Seven novel and six de novo PHEX gene mutations in patients with hypophosphatemic rickets. Int J Mol Med. 2016;38(6):1703–14.

    Article  CAS  Google Scholar 

  46. Sheng X, Chen X, Lei B, Chen R, Wang H, Zhang F, et al. Whole exome sequencing confirms the clinical diagnosis of Marfan syndrome combined with X-linked hypophosphatemia. J Transl Med. 2015;13:179.

    Article  CAS  Google Scholar 

  47. Freudlsperger C, Hoffmann J, Castrillon-Oberndorfer G, Engel M. Bilateral coronal and sagittal synostosis in X-linked hypophosphatemic rickets: a case report. J Craniomaxillofac Surg. 2013;41(8):842–4.

    Article  Google Scholar 

  48. Rowe PS, Oudet CL, Francis F, Sinding C, Pannetier S, Econs MJ, et al. Distribution of mutations in the PEX gene in families with X-linked hypophosphataemic rickets (HYP). Hum Mol Genet. 1997;6(4):539–49.

    Article  CAS  Google Scholar 

  49. Le Moual H, Devault A, Roques B, Crine P, Boileau G. Identification of glutamic acid 646 as a zinc-coordinating residue in endopeptidase-24.11. J Biol Chem. 1991;266(24):15670–4.

    Article  Google Scholar 

  50. Furlong R. Refining the splice region. Nat Rev Genet. 2018;19(8):470–1.

    Article  CAS  Google Scholar 

  51. Anna A, Monika G. Splicing mutations in human genetic disorders: examples, detection, and confirmation. J Appl Genet. 2018;59(3):253–68.

    Article  CAS  Google Scholar 

  52. Ni X, Li X, Zhang Q, Liu C, Gong Y, Wang O, et al. Clinical characteristics and bone features of autosomal recessive Hypophosphatemic rickets type 1 in three Chinese families: report of five Chinese cases and review of the literature. Calcif Tissue Int. 2020;107(6):636–48.

    Article  CAS  Google Scholar 

Download references


We would like to thank Dr. Lim Poi Giok, Department of Paediatrics, Hospital Kuala Lumpur; patients and parents for their support for the study. We would also like to thank the Ministry of Higher Education Malaysia for funding this study under the grant FRGS/1/2016/SKK02/UPM/02/2.


This work was supported by the Ministry of Higher Education, Malaysia [grant number FRGS/1/2016/SKK02/UPM/02/2].

Author information

Authors and Affiliations



NT, THT, KT, and MLK contributed to the conception and design of the study. NT performed the experiments and drafted the manuscript. NT and KL analyzed the data. MLK, THT, and KT revised the intellectual and clinical content of the manuscript. THT acquired funding for the study. All authors contributed to manuscript revision, read, and approved the submitted version.

Corresponding authors

Correspondence to Tzer Hwu Ting or Karuppiah Thilakavathy.

Ethics declarations

Ethics approval and consent to participate

All methods were carried out in accordance with the ethical standards as laid down in the 1964 Declaration of Helsinki and its later amendments or comparable ethical standards. The study was approved by the Medical Research and Ethics Committee (MREC), Ministry of Health, Malaysia (reference number: NMRR-17-612-33842 (IIR)). Informed consent was obtained from the parents of patients.

Consent for publication

Informed consent for the study and publication was obtained from the parents of patients.

Competing interests

The authors declare that they have no competing interests

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Table 1.

Thermal cycling for PCR amplification of PHEX (Exons 8 and 19) and DMP1 (intron 2). A peqSTAR block thermal cycler (VWR, Radnor, PA, USA) was used. Supplementary Table 2. High resolution melting (HRM) analysis. Thermal cycling followed by Melting curve analysis was performed on a LightCycler 480 (Roche, Switzerland). Supplementary Fig. 1. Temperature shifted difference plots of patients and control group. Supplementary Table 3. Pathogenic prediction of the PHEX and DMP1 variants. Different prediction software tools were used and the pathogenicity of the variants was assigned based on ACMG.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tavana, N., Ting, T.H., Lai, K. et al. Whole exome sequencing identifies two novel variants in PHEX and DMP1 in Malaysian children with hypophosphatemic rickets. Ital J Pediatr 48, 193 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Hypophosphatemic rickets
  • Phosphate
  • PHEX
  • DMP1
  • Mutation
  • Whole exome sequencing
  • Novel