- Open Access
RAB33B and PCNT variants in two Pakistani families with skeletal dysplasia and short stature
BMC Musculoskeletal Disorders volume 22, Article number: 630 (2021)
Skeletal dysplasia is a heterogeneous group of disorders resulting from different genetic variants in humans. The current study was designed to identify the genetic causes of skeletal dysplasia and short stature in two consanguineous families from Pakistan, both comprised of multiple affected individuals. Patients in one family had proportionate short stature with reduced head circumference while affected individuals in the other family had disproportionate short stature.
Clinical data were obtained and radiological examinations of the index patients were completed. Whole genome sequencing for probands from both families were performed followed by Sanger sequencing to confirm segregation of identified variants in the respective families. In-silico pathogenicity score prediction for identified variant and amino acid conservation analysis was completed.
Whole Genome Sequencing identified a known biallelic variant c.6176_6189delGTCAGCTGCCGAAG; p.(Gln2060ArgfsTer48) in PCNT gene and a novel biallelic variant c.174delC; p.(Asp60ThrfsTer7) in RAB33B gene respectively in affected members of the two families. Clinical imaging revealed platyspondyly and varus deformity in the legs of the affected members in the first family. Radiographs indicated severe platyspondyly, genu valgus deformity of legs and pectus carinatum for the patients in the second family.
In this study we report the phenotypes and genetic variants in two unrelated families with two distinct forms of skeletal dysplasia. This study strengthens the previous findings that patients harboring PCNT variants are phenotypically homogeneous and also extends the genotypic spectrum of RAB33B variants.
Skeletal dysplasias are a heterogeneous group of disorders, each manifesting with variable degrees of skeletal anomalies. In many cases, stunted growth and metabolic abnormalities are also observed. Microcephalic osteodysplastic primordial dwarfism type II and Smith–McCort dysplasia are two distinct autosomal recessive skeletal dysplasias with severe short stature.
Microcephalic osteodysplastic primordial dwarfism type II (MOPDII) (OMIM 210,720) is caused by pathogenic variants of PCNT (OMIM 605,925), which are also associated with other human disorders including Seckel syndrome . MOPDII and Seckel syndrome are distinguishable even though they have some common features such as disproportionate short stature and specific skeletal deformities. Pericentrin (PCNT), encoded by PCNT gene, is a highly coiled coil protein, with a conserved C terminal region responsible for interaction with other protein components .
Variants in RAB33B have been identified in patients with Smith–McCort dysplasia (SMC) negative for DYM mutations [3, 4]. RAB33B belongs to a group in the Rab family which includes small GTP binding proteins. These proteins have a crucial role in cellular trafficking such as exocytosis and endocytosis [3, 4]. Smith–McCort dysplasia (SMC) (OMIM 607,326) and Dyggve–Melchior–Clausen syndrome (DMC) (OMIM 223,800) are two allelic syndromes caused by variants of DYM . Both Smith–McCort dysplasia (SMC) and Dyggve–Melchior–Clausen syndrome (DMC) share same distinctive skeletal features such as platyspondyly, abnormalities in epiphyses and metaphyses, and irregular contours of iliac crest along with spinal cord compression. SMC and DMC are distinguishable on the basis of intellectual abilities with cognition being normal in SMC .
In the present study we investigated two families with multiple affected individuals exhibiting skeletal dysplasia. We performed whole genome sequencing and identified the pathogenic alleles segregating with the disorders.
Two families ZFD-01 and ZFD-02 (Fig. 1A-B) were identified from Shahkot region of Pakistan after study approval by the Institutional Review Board of School of Biological Sciences, University of the Punjab, Lahore. Written informed consents were obtained from the participants and parents of minors. Clinical and family histories were recorded. Radiographs for hand, spine and legs were obtained. DNA was extracted by a standard protocol involving sucrose lysis, salting out and isopropanol precipitation from blood samples of the participants . All procedures were performed in accordance with relevant guidelines.
Whole genome sequencing
Whole genome sequencing was performed on samples of affected individuals IV:1 of ZFD-01 and IV:3 of ZFD-02 respectively to identify the cause of the disorders. DNA libraries were constructed using the IlluminaTruSeq PCR-free method. Pair-end reads (2 × 150 bp) were obtained by sequencing on the HiSeqX instrument (Illumina) at the SciLifeLab facility, Stockholm, Sweden. The average coverage of reads was 30X and data were processed according to the Clinical Genetics laboratory, Centre of Molecular Medicine, Karolinska Institutet in-house pipeline. Burrow-Wheeler Aligner (BWA) was used for mapping the reads to the human genome (assembly b37); Genome Analysis Toolkit (GATK) was used for duplicate marking, variant calling and joint genotyping. Annotation of variants was performed using Variant Effect Predictor (VEP).Variant prioritization was based on: 1) autosomal recessive pattern of inheritance (based on family history and the known consanguinity in the family), 2) minor allele frequency (MAF) less than 0.01 in the public databases gnomAD as well as SweGen (Karolinska Institutet in house database), and 3) impact severity based on the variant being frameshift, introducing or deleting stop codons, affecting splice sites and missense variants with high pathogenicity scores.
The shortlisted variants were checked for segregation in each family by Sanger sequencing in all available DNA samples. Primers were designed using Primer3.0 (https://bioinfo.ut.ee/primer3-0.4.0/).
The two affected individuals in family ZFD-01(IV:1 and IV:3) had proportionate short stature and typical features of MOPDII (Fig. 1C). The heights of the affected individuals, aged 27 and 20 years, were below -6.8 SD (Table 1). The head circumference was below the normal mean for males (58.4 cm), but cognition was not impaired. Their joints and ability to walk were unaffected. Radiographs of individual IV:3 indicated severe platyspondyly in vertebral bodies, along with a slight varus deformity in legs, abnormal metaphyses in distal femurs and slight abnormality in carpal structure (Fig. 1D-F). The phenotypically normal individuals had no skeletal deformities and exhibited normal growth and heights (Table 1).
The affected individuals of family ZFD-02 had disproportionately short statures (Fig. 1G). Individuals IV:2, IV:3, IV:8 and IV:11 were of ages 33, 29, 22 and 10 years respectively (Table 1). The symptoms became evident after the age of 2 years of age. They had restricted joint mobility and valgus deformity in knee joints and barrel shaped chest (Fig. 1H-J). As subjects aged, the disorder grew progressively with movement becoming more difficult. The oldest affected individual had a severely affected and painful gait. Severe platyspondyly, genu valgus deformity of legs, cupping of distal ribs, pectus carinatum, and irregular iliac crests were observed in radiographs of individual IV:11 (Fig. 1K-N).
Whole genome sequencing of the affected individual IV:1 of family ZFD-01 revealed eight variants (Table S1) after analysis and filtration of data. Out of these, only a deletion of 14 bp, c.6176_6189del GTCAGCTGCCGAAG; p.(Gln2060ArgfsTer48) in exon 30 of PCNT (NM_001315529.2) (Fig. 1O) segregated with the disorder. This finding was consistent with a diagnosis of MOPDII for the patients.
In family ZFD-02, whole genome sequencing of the affected individual IV:3 revealed four variants after filtration by applying different criteria (Table S2). A frameshift variant in RAB33B (NM_031296.1), c.174delC; p.(Asp60ThrfsTer7) segregated fully with the phenotype (Fig. 1P) and confirmed the diagnosis of SMC2 (OMIM 615,222). This variant was absent in all public databases. The RAB33B (NM_031296.1), c.174delC; p.(Asp60ThrfsTer7) variant has been deposited in LOVD with ID 0,000,670,684.
The PCNT gene is comprised of 47 coding exons which encodes a 3,336 amino acids long coiled coil protein flanked at terminals by non-helical regions. The mutation identified in family ZFD-01 in PCNT gene is presumed to produce a non-functional truncated protein consisting of 2108 amino acids, if the mRNA escapes nonsense mediated decay. Similarly, the mutation identified in family ZFD-02 in RAB33B gene will severely truncate the protein to 67 amino acids only, if the mRNA is translated. However, both mutant transcripts of PCNT and RAB33B are likely to undergo nonsense-mediated decay and no proteins will be produced.
In this study we used whole genome sequencing of the index subjects in two families with rare forms of skeletal dysplasia to identify the disease-causing gene defects. In family ZFD-01 we identified a homozygous frameshift variant in PCNT, confirming the diagnosis of MOPDII. The affected individuals had severe but proportionate short stature (heights SD < -6.8). Their head circumference was also reduced as observed in previously described patients with MOPDII. Around 50% of affected individuals with MOPDII develop cerebral neurovascular abnormalities . Dentition problems are also common. However, in our patients, cognition was not impaired and both affected individuals had normal dentition.
PCNT is located on chromosome 21 and encodes pericentrin, a centrosomal protein responsible for nucleation of mitotic spindle and binding to calmodulin . It networks with its other components including the γ-tubulin ring complex and provides an active configuration for the proteins, essential for the Microtubular organizing center activities. Disruptions in astral microtubule and spindle orientation can lead to the malfunctioning of Pericentriolar Matrix, which plays a crucial role in organization of spindles and can lead to mitotic arrest or cell death [8, 9].
Various pathogenic variants of PCNT have been identified in MOPDII and Seckel syndrome  (www.hgmd.cf.ac.uk accessed June 2021). The 14 nucleotide deletion, c.6176_6189 del15 in exon 30, revealed by analysis of the whole genome sequencing data in family ZFD-01, has been previously reported in a Pakistani female patient with similar phenotypic features. She had many skeletal manifestations as described for MOPDII as well as polycystic ovaries .
In the second family ZFD-02 a novel homozygous variant c.174delC, p.(Asp60ThrfsTer7) was identified in RAB33B. The affected individuals had typical features of Smith–McCort dysplasia (SMC). Skeletal features on radiographs, including severe platyspondyly, short and broad metacarpals, abnormal carpal bones and lower limb valgus deformity were in concordance with the previously reported cases of SMC phenotypes [4, 5]. These also included barrel chests, limited joint movement, lower limb deformities and normal intellect.
RAB33B (NM_031296.2) is located on chromosome 4q31.1, comprises 2 exons and encodes RAB33B which is localized to Golgi complex . RAB33B is a type of Rab GTPase. Rab GTPase membrane bound proteins are involved in membrane transport. It is involved in retrograde Golgi transport (from Golgi to ER) of proteins . Moreover, it also functions in the formation of autophagosomes by binding with the Atg12-5/16L complex, contributing to macroautophagy .
To date, only seven variants of RAB33B have been identified in patients with a similar phenotype, all from different ethnic backgrounds [3,4,5] (www.hgmd.cf.ac.uk accessed June 2021). A missense mutation in RAB33B can cause marked deficiency of RABB33B inside Golgi and affects the vesicular trafficking of Golgi complex . This novel variant c.174delC p.(Asp60ThrfsTer7) identified in family ZFD-02, will most likely mark the mRNA for nonsense mediated decay, or if the mutant transcript escapes nonsense mediated decay, premature truncation of the protein within the GTPase domain will occur (Fig. 1Q).
The findings of this present study suggest phenotypic homogeneity of patients harboring PCNT or RAB33B variants. The c.174delC; p.(Asp60ThrfsTer7) variant identified in this research adds to the allelic spectrum of RAB33B variants of this rare genetic disorder. Our study further shows that for genetic diagnosis of rare recessive skeletal dysplasias, whole genome sequencing of the index case provides a suitable diagnostic approach.
Availability of data and materials
Data generated in this study is available upon reasonable request to the corresponding authors. We have deposited the sequence variant data in LOVD (ID 0,000,670,684). Due to ethical concerns, our IRB has not approved the deposition of whole genome sequencing or exome sequencing data in a public database.
Willems M, Genevieve D, Borck G, Baumann C, Baujat G, Bieth E, et al. Molecular analysis of pericentrin gene (PCNT) in a series of 24 Seckel/microcephalic osteodysplastic primordial dwarfism type II (MOPD II) families. J Med Genet. 2010;47(12):797–802.
Li F-F, Wang X-D, Zhu M-W, Lou Z-H, Zhang Q, Zhu C-Y, et al. Identification of two novel critical mutations in PCNT gene resulting in microcephalic osteodysplastic primordial dwarfism type II associated with multiple intracranial aneurysms. Metab Brain Dis. 2015;6(30):1387–94.
Salian S, Cho TJ, Phadke SR, Gowrishankar K, Bhavani GS, Shukla A, et al. Additional three patients with Smith-McCort dysplasia due to novel RAB33B mutations. Am J Med Genet A. 2017;173(3):588–95.
Alshammari MJ, Al-Otaibi L, Alkuraya FS. Mutation in RAB33B, which encodes a regulator of retrograde Golgi transport, defines a second Dyggve–Melchior–Clausen locus. J Med Genet. 2012;49(7):455–61.
Dupuis N, Lebon S, Kumar M, Drunat S, Graul-Neumann LM, Gressens P, et al. A Novel RAB33B Mutation in S mith–M c C ort Dysplasia. Hum Mutat. 2013;34(2):283–6.
Grimberg J, Nawoschik S, Belluscio L, McKee R, Turck A, Eisenberg A. A simple and efficient non-organic procedure for the isolation of genomic DNA from blood. Nucleic Acids Res. 1989;17(20):8390.
Lorenzo-Betancor O, Blackburn PR, Edwards E, Vázquez-do-Campo R, Klee EW, Labbé C, et al. PCNT point mutations and familial intracranial aneurysms. Neurology. 2018;91(23):e2170–81.
Wang Y, Dantas TJ, Lalor P, Dockery P, Morrison CG. Promoter hijack reveals pericentrin functions in mitosis and the DNA damage response. Cell Cycle. 2013;12(4):635–46.
Delaval B, Doxsey SJ. Pericentrin in cellular function and disease. The Journal of cell biology. 2010;188(2):181–90.
Starr T, Sun Y, Wilkins N, Storrie B. Rab33b and Rab6 are functionally overlapping regulators of Golgi homeostasis and trafficking. Traffic. 2010;11(5):626–36.
We are thankful to the families for participating in this study. The authors express their gratitude to Clair Wootan and Dr. Edward R Wilcox for critical input to the manuscript. The authors acknowledge support from the National Genomics Infrastructure in Stockholm funded by Science for Life Laboratory, the Knut and Alice Wallenberg Foundation and the Swedish Research Council, and SNIC/Uppsala Multidisciplinary Center for Advanced Computational Science for assistance with massively parallel sequencing and access to the UPPMAX computational infrastructure.
This research was funded by the IRSIP Higher Education Commission (HEC) Pakistan (NA) and HEC grant to School of Biological Sciences, Koshish foundation USA and University of the Punjab (SN) & Swedish Research Council, Academy of Finland, Folkhälsan Research Foundation, Novo Nordisk Foundation, and Sigrid Jusélius Foundation (OM). Open access funding provided by Karolinska Institute.
Ethics approval and consent to participate
The study was approved by the Institutional Review Board of School of Biological Sciences, University of the Punjab, Lahore. Written informed consents were obtained from the participants and parents for their minor children. All procedures were performed in accordance with relevant guidelines.
Consent for publication
Written informed consent was obtained from patients and parents of minors for publication.
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Ain, N.u., Fatima, Z., Naz, S. et al. RAB33B and PCNT variants in two Pakistani families with skeletal dysplasia and short stature. BMC Musculoskelet Disord 22, 630 (2021). https://doi.org/10.1186/s12891-021-04503-2
- Skeletal dysplasia
- Short stature
- Whole genome sequencing
- Smith-McCort dysplasia