- Open Access
Novel mutation G324C in WNT1 mapped in a large Pakistani family with severe recessively inherited Osteogenesis Imperfecta
- Mehran Kausar1, 2,
- Saima Siddiqi2,
- Muhammad Yaqoob3,
- Sajid Mansoor4, 10,
- Outi Makitie5, 6,
- Asif Mir7,
- Chiea Chuen Khor8,
- Jia Nee Foo†8, 9Email author and
- Mariam Anees†1Email authorView ORCID ID profile
© The Author(s). 2018
- Received: 10 September 2018
- Accepted: 24 October 2018
- Published: 17 November 2018
The Correction to this article has been published in Journal of Biomedical Science 2019 26:31
Osteogenesis imperfecta (OI) is a clinically and genetically heterogeneous disease with skeletal fragility and variable extra-skeletal manifestations. To date several point mutations in 18 different genes causing different types of OI have been identified. Mutations in WNT1 compromise activity of the osteoblasts leading to disturbed bone mass accrual, fragility fractures and progressive skeletal abnormalities. The present study was conducted to determine the underlying genetic cause of an autosomal recessive skeletal dysplasia in a large consanguineous family from Chinute, Pakistan.
Materials and methods
Blood was collected from 24 individuals of affected family along with clinical data. Homozygosity mapping was performed to confirm consanguinity. SNPs were identified, followed by whole exome and Sanger sequencing. In silico characterization of WNT1 mutation was performed using multiple platforms.
Nine affected family members exhibited severe bone deformities, recurrent fractures, short stature and low bone mineral density. SNP array data revealed homozygous segments > 1 Mb in length accounting for 2.1–12.7% of the genome in affected individuals and their siblings and a single 6,344,821 bp homozygous region in all affected individuals on chromosome 12q12-q13. This region includes two potential OI candidate genes WNT1 and VDR. We did whole-exome sequencing for both genes in two patients and identified a novel damaging missense mutation in exon 4 of WNT1: c.1168G > T (NM_005430) resulting in p.G324C. Sanger sequencing confirmed segregation of mutation with the disease in family.
We report a novel mutation responsible for OI and our investigation expands the spectrum of disease-causing WNT1 mutations and the resulting OI phenotypes.
- WNT signaling
- Whole-exome sequencing
- Osteogenesis imperfecta
Osteogenesis imperfecta (OI) is a genetic disorder with severe bone abnormalities and is the most common form of primary osteoporosis in children . Clinical phenotype may vary from prenatal lethal forms to milder forms with only a few fractures during life time. Severely deforming forms may have lifelong bone fractures along with visible skeletal deformities, short stature, dentinogenesis imperfecta, blue sclerae, and hearing loss . The incidence of OI is approximately 1 in 15,000–20,000 live births and most cases (~ 90%) follow an autosomal dominant pattern of inheritance with mutations in either COL1A1(MIM 120150) or COL1A2 (MIM 120160) . Rest of the 10% cases follow the autosomal recessive pattern of inheritance with mutations in FKBP10, P3H1, LEPRE1, PLOD2, PPIB, SERPINF1, SERPINH1, SP7, BMP1, TMEM38B, CRTAP, CREB3L1, IFITM5 and/or WNT1. LRP5 gene may follow any of the autosomal dominant or recessive patterns; PLS3 gene follows x-linked mode of inheritance [3–7]. Most of these genetic defects induce a qualitative or quantitative change in type I collagen, the main structural protein in the skeleton. Impaired bone-matrix formation and mineralization ultimately devastate the skeletal flexibility and integrity . WNT1 (MIM 164820) has an integral role in osteoblast differentiation and is involved in the regulation of the WNT induced beta-catenin signaling pathway. Biallelic defects in WNT1 result in decreased expression of type I collagen . Mutations encoding other components of the signaling pathway have been reported also in other disorders with abnormal bone mass, suggesting that the pathway has an integral role in the regulation of bone homeostasis .Only a few cases of OI caused by WNT1 mutations have been described thus far and the spectrum of clinical manifestations remains inadequately characterized. Here we present a large consanguineous Pakistani family in which several family members present with severe deforming OI caused by a novel homozygous mutation in WNT1. Wnt1-related osteogenesis imperfecta has symptoms including scoliosis, developmental delay and short stature. WNT1 (Wnt Family Member 1) is an important gene associated with Osteogenesis Imperfecta.
We identified a large consanguineous Pakistani family with multiple family members affected by severe bone deforming OI. Blood samples for genetic analyses were collected from all available family members after informed consent. A pedigree was constructed to draw the relationships and the disease status.
Ethics, consent and permission to publish
This study was approved by the institutional review board of the Institute of Biomedical and Genetic Engineering (IBGE), Islamabad, Pakistan. All adult participants gave written informed consent to participate in this study. For minors, consent was obtained from their legal guardians / parents. Permission was also obtained from the patients to use their data and pictures for publishing.
Clinical examination of patients was performed at the Indoor Patient Department (IPD) of Children’s Hospital and radiography was done at Radiology Department of Children’s Hospital, Lahore, Pakistan.
Blood samples were drawn from 24 individuals from the affected family that included affected (n = 9)and unaffected(n = 15) individuals. Genomic DNA was extracted using standard Sambrook protocol . Quality of DNA was assessed by using Qubit dsDNA HS (High Sensitivity) Assay kit (Invitrogen) and 24 samples were genotyped on the IlluminaOmniExpress 24v1–0-A BeadChip array for a total of 716,503 genetic markers. After removal of single nucleotide polymorphisms (SNPs) that were monomorphic orfailed genotyping in > 1 sample, 392,260 genome-wide SNP markers remained for further analysis. We confirmed the reported familial relationships among genotyped samples using PLINK identity by descent (IBD) analysis (−genome) [10, 11]. We then scanned the data for large homozygous segments (> 1 Mb) which are shared among affected individuals but not for the unaffected individuals. Homozygosity mapping was conducted using PLINK v1.07 using the default settings .
Whole exome sequencing (WES) and sanger sequencing
Targeted enrichment was performed on 1 μg of genomic DNA from two affected individuals (III:5 and III:15) using the Nimblegen SeqCap EZ Exome v3 kit and barcoded for sequencing on a single lane of a multiplexed 2 × 101 bp sequencing run on the Illumina HiSeq 2000platform. Each individual was sequenced to a mean coverage of 45.2–46.7 reads per target base, with 97% of the target exome covered by 10 or more reads. Reads were mapped using BWA v1.7  and variants were called using the GATK v2 Unified Genotyper following the recommended guidelines by GATK ‘Best practices for variant calling v3’ . We used the following primer sequences for the segregation of detected mutations in WNT1 by polymerase chain reaction and Sanger sequencing: WNT1-ForwardGAAATCGCCCAACTTCTGCA and WNT1-Reverse AGTGCTAGCGAGTCTGTTTGG.
In silico predictions and validation of novel WNT1 mutation G324C
I-Tasser [14, 15] was used for three dimensional (3D) model predictions. Protein 3D models reliability was checked using RAMPAGE server . STITCH4 database  was used to predict functional protein partners. Pockets on three-dimensional structures of proteins were identified using CASTp server [18, 19]. Meta SNP  and I-Mutant2.0, were utilized to estimate effect of mutation on stability of protein and to determine whether the mutation has an impact on normal function of protein or not. Docking analysis was carried out using PatchDock server . The refinement of first 10 docked complexes obtained through PatchDock was done using FireDock [22, 23]. Representations (2-Dimensional) and analysis of protein-ligand interaction complexes was done using LIGPLOT . Meta SNP was also utilized to analyze effect of mutation on protein.
OI phenotype examination
Patient 1 (III:9) was a 11 years old male, who presented with a history of multiple fractures after mild trauma since the age of 6 months. His first fracture of right clavicle occurred at the age of 6 months, followed by fracture of right tibial shaft at the age of 8 years, and fracture of right and left femora at 9 years. Pregnancy and birth were uneventful. His parents were distant cousins (IBD proportion < 6%). Presently, the index individual is wheelchair bound because of multiple fractures, deformities and weakness of lower extremities. His weight is 20 kg, height 106 cm and head circumference 52 cm. Physical examination revealed brachycephalic head, flat face, mild blueness of sclerae, right eye squint, short neck, marked thinness of upper and lower extremities, anterior angulation of both femora and right tibia, and flat feet (Fig. 1b). Skeletal radiographs revealed generalized osteopenia, bowing of clavicles, compression of thoracic vertebrae, narrowing of intercostal spaces, and bowing of long bones.
Patient 2(III:15) was a 10 years old male, wheelchair bound, measuring 18 kg in weight, 107 cm in height and 51 cm in head circumference. Pregnancy and delivery were uneventful. Parents were also distant cousins (IBD proportion ~ 6%).His health remained poor since birth. He was able to walk but never became able to run. At the age of 8 years he sustained fractures of right humerus and right tibia and femur after trivial traumas (Fig. 1c). Clinical examination showed triangular face, normal teeth, right eye squint, normal sclerae, and wide protruding chest with increased antero-posterior diameter. Extremities were thin and showed bilateral mild angulation of proximal humeri and marked bowing of right femur. His feet were flat. Radiological examination showed generalized osteopenia, platyspondyli in all vertebrae, narrowing of inter-costal spaces, globular pelvis, bowing of distal parts of right tibia and fibula, and mild bowing of right radius and ulna (Fig. 1d and e).
Homozygous segments of > 1 Mb in length accounted for 2.1–12.7% of the genome in all 9 affected individuals and their full siblings, confirming that these individuals are the offspring of consanguineous unions. We identified a single 6,344,821 base pair segment on chromosome 12q12-q13 (chr12:46,084,699-52,429,520) that was homozygous only in affected individuals (Additional file 2: Table S2). Two candidate genes for bone density were located in the homozygous segment, the vitamin D receptor VDR and the wingless-type integration site family, member 1 WNT1.
Whole exome sequencing and segregation analysis
For identification of causal mutation, we performed whole exome sequencing analysis in two affected cousins III:5 and III:15 (Additional file 3: Table S3). A total of 43,494 variants were found in these two individuals, out of which 29,745 were coding and 13,287 were nonsynonymous, frameshift or splice site variants in well-annotated transcripts. Of these, 1319 were rare, either absent or present in less than 1% of all populations in HapMap, 1000 genomes populations  and the NHLBI exome variant server (EVS) databases (URL: http://evs.gs.washington.edu/EVS/). Of these, 20 were homozygous in both affected cousins, 3 resided within the shared homozygous interval on chromosome 12 and only one was predicted to be damaging and not found in any of the above SNP databases or the Exome Aggregation Consortium (ExAC) database. This was identified as a novel damaging mutation in WNT1 c.1168G > T (RefSeq NM_005430) resulting in p.Gly324Cys. We subsequently Sanger sequenced samples of 24 members of the family and verified that this mutation segregated perfectly with the disease (Additional file 4: Figure S1). This mutation has not been reported previously in any of the publications or online databases including ‘ClinVar’, ‘HGMD(R)’ and ‘ExAC’. No pathogenic mutations were found in CCNT1 or VDR genes; there was only a synonymous SNP in VDR gene with no amino acid change.
In silico characterization of WNT1 mutation G324C
PANTHER: Between 0 and 1. (If > 0.5 mutation is predicted Disease)
PhD-SNP: Between 0 and 1. (If > 0.5 mutation is predicted Disease)
SIFT: Positive Value (If > 0.05 mutation is predicted Neutral)
SNAP: Output normalized between 0 and 1 (If > 0.5 mutation is predicted Disease)
Meta-SNP: Between 0 and 1. (If > 0.5 mutation is predicted Disease)
Summary of residues involved in WNT1 docking interactions (Wild & Mutant)
Hydrogen Bond Interactions
Hydrophobic Bond Interactions
Wild WNT1-Axin 1
Val383, Glu391, Arg395, Val399, Gln400
Arg98, Cys218, Arg344, Cys347, Cys352
Glu384, Pro385, Phe388, Ala389, Leu392, Leu396, Arg401
Arg215, Gln216, Glu217, Arg235, Thr238, Arg240, Asp244, Arg247, Pro320, Thr348, Phe349, His350, Val355
Mutated WNT1-Axin 1
Gln386, Phe388, Glu390, Leu396, Glu397, Gln400
Lys51, Trp95, Arg98, Arg337, Glu368
Val383, Glu391, Leu392, Val399
Asn29, Ser50, Arg91, Lys94, Gly312, Arg335
Mutations in WNT1have been associated with extreme phenotypes of osteoporosis, such as prenatally-onset severe OI, idiopathic juvenile osteoporosis and pregnancy-associated osteoporosis, confirming WNT1’s central role in regulation of bone strength. While biallelic mutations, like in our family, result in severe deforming OI, heterozygous mutations result in childhood-onset or early adulthood-onset osteoporosis with spinal and peripheral fractures, low BMD and low bone turnover, but normal growth and no deformities. Several mutations have been reported. Our patients with homozygous WNT1 mutations had severe growth failure and early-onset skeletal fragility leading to severe deformities and loss of ambulation. This skeletal phenotype is consistent with other cases of OI caused by homozygous WNT1 mutations . No skeletal phenotype was evident in the heterozygous mutation carriers in our family, but no systematic phenotyping with BMD assessment was performed. None of the reported patients received bisphosphonate treatment and it remains unknown whether early-onset therapy would be as efficient as in type I collagen related OI in improving growth and preventing fractures and deformities.
We report a novel mutation G324C in WNT1 gene responsible for severe recessively inherited Osteogenesis Imperfecta. This is being reported for the first time in Pakistan as well as globally. Our report expands the phenotypic and genotypic spectrum of WNT1-related OI. OI is a highly divergent disease with different inheritance patterns and variable severities and further studies on the clinical and functional consequences of different WNT1 mutations may provide important insight into the roles of this gene in bone function and development. Optimal treatment remains to be established in future studies.
We thank all the family members who participated in the study especially Mr. Mohammad Javaid sb who helped and provided all the information required and Mr. Sajid Hussain for helping in collecting the samples. We also thank Wee Yang Meah and Hui-Meng Soo from the Genome Institute of Singapore for their help with the genotyping and sequencing of the samples.
The genetic analysis was supported by the Agency for Science, Technology and Research, Singapore (C.C.K.) and a Singapore National Research Foundation Fellowship (NRF-NRFF2016–03; to J.N.F.). The research is partly supported by Academy of Finland, Sigrid Jusélius Foundation and Folkhälsan Research Foundation, Finland (O.M.). IRSIP Scholarship was awarded to M.K. by Higher Education Commission, Pakistan.
Availability of data and materials
Data and materials related to this work are available upon request.
MK and SS collected the samples and helped in clinical evaluation. SS, JNF, CCK performed the SNP array data analysis, MY and OM did the clinical evaluation of the family. MK, SS, JNF, SM, OM and MA wrote and reviewed the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
This study was approved by the institutional review board of the Institute of Biomedical and Genetic Engineering (IBGE), Islamabad, Pakistan and Nanyang Technological University (IRB-2016-08-011), Singapore.
Consent for publication
All adult participants gave written informed consent to participate in this study. For minors, consent was obtained from their legal guardians/parents. Permission was also obtained from the patients to use their data and pictures for publishing.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Forlino A, Marini JC. Osteogenesis imperfecta. Lancet. 2016;387(10028):1657–71.View ArticlePubMedGoogle Scholar
- Marini JC, Forlino A, Cabral WA, Barnes AM, San Antonio JD, Milgrom S, et al. Consortium for osteogenesis imperfecta mutations in the helical domain of type I collagen: regions rich in lethal mutations align with collagen binding sites for integrins and proteoglycans. Hum Mutat. 2007;28(3):209–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Barnes AM, Chang W, Morello R, Cabral WA, Weis M, Eyre DR, et al. Deficiency of cartilage-associated protein in recessive lethal osteogenesis imperfecta. N Engl J Med. 2006;355(26):2757–64.View ArticlePubMedGoogle Scholar
- Morello R, Bertin TK, Chen Y, Hicks J, Tonachini L, Monticone M, et al. CRTAP is required for prolyl 3-hydroxylation and mutations cause recessive osteogenesis imperfecta. Cell. 2006;127(2):291–304.View ArticlePubMedGoogle Scholar
- Symoens S, Malfait F, D’hondt S, Callewaert B, Dheedene A, Steyaert W, et al. Deficiency for the ER-stress transducer OASIS causes severe recessive osteogenesis imperfecta in humans. Orphanet J Rare Dis. 2013;8(1):154.View ArticlePubMedPubMed CentralGoogle Scholar
- Cho TJ, Lee KE, Lee SK, Song SJ, Kim KJ, Jeon D, et al. A single recurrent mutation in the 5′-UTR of IFITM5 causes osteogenesis imperfecta type V. Am J Hum Genet. 2012;91(2):343–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Semler O, Garbes L, Keupp K, Swan D, Zimmermann K, Becker J, et al. A mutation in the 5'-UTR of IFITM5 creates an in-frame start codon and causes autosomal-dominant osteogenesis imperfecta type V with hyperplastic callus. Am J Hum Genet. 2012;91:349–57.View ArticlePubMedPubMed CentralGoogle Scholar
- Keupp K, Beleggia F, Kayserili H, Barnes AM, Steiner M, Semler O, et al. Mutations in WNT1 cause different forms of bone fragility. Am J Hum Genet. 2013;92(4):565–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Sambrook J, Russell DW. The condensed protocols from molecular cloning: a laboratory manual; 2006.Google Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75.View ArticlePubMedPubMed CentralGoogle Scholar
- Sobel E, Lange K. Descent graphs in pedigree analysis: applications to haplotyping, location scores, and marker-sharing statistics. Am J Hum Genet. 1996;58(6):1323–37.PubMedPubMed CentralGoogle Scholar
- Li H, Durbin R. Fast and accurate long-read alignment with burrows-wheeler transform. Bioinformatics. 2010;26(5):589–95.View ArticlePubMedPubMed CentralGoogle Scholar
- 1000 Genomes Project Consortium, Abecasis GR, Auton A, Brooks LD, MA DP, Durbin RM, et al. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491(7422):56–65.View ArticleGoogle Scholar
- Roy A, Kucukural A, Zhang Y. I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010 Apr;5(4):725–38.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang Y. I-TASSER server for protein 3D structure prediction. BMC Bioinformatics. 2008;9:40.View ArticlePubMedPubMed CentralGoogle Scholar
- Lovell SC, Davis IW, Arendall WB 3rd, de Bakker PI, Word JM, Prisant MG, et al. Structure validation by Calpha geometry: phi,psi and Cbeta deviation. Proteins. 2003;50(3):437–50.View ArticlePubMedGoogle Scholar
- Kuhn M, Szklarczyk D, Franceschini A, von Mering C, Jensen LJ, Bork P. STITCH 3: zooming in on protein-chemical interactions. Nucleic Acids Res. 2012;40(Database issue):D876–80.View ArticlePubMedGoogle Scholar
- Binkowski TA, Naghibzadeh S, CASTp LJ. Computed atlas of surface topography of proteins. Nucleic Acids Res. 2003;31(13):3352–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Dundas J, Ouyang Z, Tseng J, Binkowski A, Turpaz Y, Liang J. CASTp: computed atlas of surface topography of proteins with structural and topographical mapping of functionally annotated residues. Nucleic Acids Res. 2006;34(Web Server):W116–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Capriotti E, Calabrese R, Fariselli P, Martelli PL, Altman RB, Casadio R. WS-SNPs&GO: a web server for predicting the deleterious effect of human protein variants using functional annotation. BMC Genomics. 2013;14(Suppl 3):S6.View ArticlePubMedPubMed CentralGoogle Scholar
- Schneidman-Duhovny D, Inbar Y, Nussinov R, Wolfson HJ. PatchDock and SymmDock: servers for rigid and symmetric docking. Nucleic Acids Res. 2005;33(Web Server):W363–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Andrusier N, Nussinov R, Wolfson HJ. FireDock: fast interaction refinement in molecular docking. Proteins. 2007;69(1):139–59.View ArticlePubMedGoogle Scholar
- Mashiach E, Schneidman-Duhovny D, Andrusier N, Nussinov R, Wolfson HJ. FireDock: a web server for fast interaction refinement in molecular docking. Nucleic Acids Res. 2008;36(Web Server):W229–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Wallace AC, Laskowski RA, Thornton JM. LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. Protein Eng. 1995;8(2):127–34.View ArticlePubMedGoogle Scholar
- Pyott SM, Tran TT, Leistritz DF, Pepin MG, Mendelsohn NJ, Temme RT, et al. WNT1 mutations in families affected by moderately severe and progressive recessive osteogenesis imperfecta. Am J Hum Genet. 2013;92(4):590–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Janda CY, Waghray D, Levin AM, Thomas C, Garcia KC. Structural basis of Wnt recognition by frizzled. Science. 2012;337(6090):59–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Laine CM, Joeng KS, Campeau PM, Kiviranta R, Tarkkonen K, Grover M, et al. WNT1 mutations in early-onset osteoporosis and osteogenesis imperfecta. N Engl J Med. 2013;368(19):1809–16.View ArticlePubMedPubMed CentralGoogle Scholar