- Open Access
Novel mutation G324C in WNT1 mapped in a large Pakistani family with severe recessively inherited Osteogenesis Imperfecta
Journal of Biomedical Sciencevolume 25, Article number: 82 (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.
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,4,5,6,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.
Material and methods
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
We examined a large consanguineous Pakistani family with nine affected members (Fig. 1a) diagnosed with OI type 3 on the basis of apparent clinical features and radiological findings. All affected individuals manifested similar features (Additional file 1: Table S1). The two index patients are described here in more detail.
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
To gain insight into the structural implications of the mutation p.G324C, we modelled the wild type and mutant proteins using the online I-Tasser server (Fig. 2a, b and c) . Models with the best c-score were selected for further analysis. Interaction sites as predicted by the CASTp results showed a decrease in available pockets in case of mutant compared to the wild type, 40 compared to 49, respectively; clearly indicating the effect of mutation upon the binding capacity of the protein (Fig. 2d). Meta SNP hires with it different predictors (PATHER, PhD-SNP, SIFT and SNAP) to calculate mutation impact on normal protein. Value reported under each prediction is given in (Fig. 2e). All predictors show it as disease causing which validate our result. Obtained calculated values are interpreted as following:
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)
With maximum score of 0.988, Axin 1 has been selected as ligand for further analysis. It is a component of the beta-catenin destruction complex required for regulating catenin beta 1 (CTNNB1) levels through phosphorylation and ubiquitination, and modulating Wnt-signaling. It controls dorso ventral patterning via two opposing effects; down-regulates CTNNB1 to inhibit the Wnt signaling pathway and ventralize embryos, but also dorsalizes embryos by activating a Wnt- independent JNK signalling pathway. Wnt signalling probably facilitates the phosphorylation of CTNNB1 and APC by GSK3B and is likely to function as a tumor suppressor. Ligplot results for the docking interactions of wild and mutant type WNT1 with ligand AXIN 1 are shown in Fig. 3a and b, respectively. Residues involved in hydrogen bonding and hydrophobic interactions are given in Table 1. G324C mutation caused a huge impact on the normal structure and conformation of the protein molecule. This can be well observed through differences in number and position of residues involved in both types of docking interactions. This is due to the alteration of the interaction site after mutation that affected protein ligand interaction. As per CASTp predictions, pocket number 49 and 17 (Fig. 4a) contain maximum residues which are involved in docking of normal protein with ligand while pocket number 12 and 38 (Fig. 4b) contain residues involved in docking interaction of mutant protein. The analysis not only identifies active sites in the WNT1 structure but also validates the docking results.
By performing homozygosity mapping and whole exome sequencing, we identified a novel homozygous mutation in WNT1 as the cause of osteogenic imperfecta in this large consanguineous family. The novel deleterious mutation segregated with the phenotype and in the homozygous individuals resulted in severe clinical findings consistent with deforming type 3OI, thereby providing basis for genetic diagnosis of OI in the family. OI comprises multiple heterogeneous skeletal phenotypes and some forms lead to marked deformities of upper and lower limbs. To date several biallelic missense and nonsense mutations in WNT1 have been identified with multiple overlapping skeletal manifestations (Fig. 5).
The p.G324C mutation is predicted to be damaging to the WNT1 because it occurs at a highly conserved position (Fig. 6). The mutated residue (position 324) is present on a small α-helix with a number of important functional residues in the vicinity. Nearby asparagine is glycosylated while cysteine residues at 315 and 330 form a disulfide linkage. The glycosylation and the disulfide linkage are important for protein stability. Therefore, the mutation at position 324 (Gly to Cys) can be detrimental to protein stability and therefore, functionality. An adjacent cysteine residue (position 325) may well present itself for a potential disulfide linkage with the mutated cysteine residue. On the other hand the mutated cysteine may engage the cystein at 315 instead and compete for a disulfide linkage with the cysteine at 330. In either case the protein stability at this site may be effected which might result in producing a dysfunctional protein. It is also important to note that two other known mutations at position 177 (Gly to Cys) and 298 (Phe to Cys) have been implicated to cause OI [8, 25]. The highly conserved glycine at 177 allow for a compact protein core while the mutated cysteine may disrupt this critical protein site . WNT proteins are characterized by a set of 22 conserved cysteine residues that form the intra-chain disulfide bridges that maintain the tertiary protein structure . This suggests that the mutations involving cysteine residues of WNT1 may have serious impact on the protein stability and function.
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.
Forlino A, Marini JC. Osteogenesis imperfecta. Lancet. 2016;387(10028):1657–71.
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.
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.
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.
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.
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.
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.
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.
Sambrook J, Russell DW. The condensed protocols from molecular cloning: a laboratory manual; 2006.
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.
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.
Li H, Durbin R. Fast and accurate long-read alignment with burrows-wheeler transform. Bioinformatics. 2010;26(5):589–95.
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.
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.
Zhang Y. I-TASSER server for protein 3D structure prediction. BMC Bioinformatics. 2008;9:40.
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.
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.
Binkowski TA, Naghibzadeh S, CASTp LJ. Computed atlas of surface topography of proteins. Nucleic Acids Res. 2003;31(13):3352–5.
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.
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.
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.
Andrusier N, Nussinov R, Wolfson HJ. FireDock: fast interaction refinement in molecular docking. Proteins. 2007;69(1):139–59.
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.
Wallace AC, Laskowski RA, Thornton JM. LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. Protein Eng. 1995;8(2):127–34.
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.
Janda CY, Waghray D, Levin AM, Thomas C, Garcia KC. Structural basis of Wnt recognition by frizzled. Science. 2012;337(6090):59–64.
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.
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.
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.
Table S1. Clinical findings in the nine affected family members with a homozygous WNT1 mutation. The “+” signs indicates the presence and severity of the symptom and “-” sign indicating the absence of the symptom. (DOCX 15 kb)
Table S2. Identification of a single homozygous segment on chr12 overlapping between all cases and none of the controls. (DOCX 15 kb)
Table S3. List of rare homozygous variants in III:5 and III:15 located within the shared homozygous segment on chr12. (DOCX 13 kb)
Figure S1. Sanger sequencing results of affected family. Both parents (II:3 and II:4) (right upper and lower panels) of Patient 1 were heterozygous carriers, arrows indicate the double peaks (heterozygosity). Patient 1 (III:9) (left lower panel) result indicates homozygous change, arrow is indicating the change. While an unaffected individual (III:11) (left upper panel) showed both correct/wild type alleles. (DOCX 106 kb)