- Open Access
A hypomorphic BMPR1B mutation causes du Pan acromesomelic dysplasia
Orphanet Journal of Rare Diseases volume 10, Article number: 84 (2015)
Grebe dysplasia, Hunter-Thompson dysplasia, and du Pan dysplasia constitute a spectrum of skeletal dysplasias inherited as an autosomal recessive trait characterized by short stature, severe acromesomelic shortening of the limbs, and normal axial skeleton. The majority of patients with these disorders have biallelic loss-of-function mutations of GDF5. In single instances, Grebe dysplasia and a Grebe dysplasia-like phenotype with genital anomalies have been shown to be caused by mutations in BMPR1B, encoding a GDF5 receptor.
We clinically and radiologically characterised an acromesomelic chondrodysplasia in an adult woman born to consanguineous parents. We sequenced GDF5 and BMPR1B on DNA of the proposita. We performed 3D structural analysis and luciferase reporter assays to functionally investigate the identified BMPR1B mutation.
We extend the genotype-phenotype correlation in the acromesomelic chondrodysplasias by showing that the milder du Pan dysplasia can be caused by a hypomorphic BMPR1B mutation. We show that the homozygous c.91C>T, p.(Arg31Cys) mutation causing du Pan dysplasia leads to a significant loss of BMPR1B function, but to a lesser extent than the previously reported p.Cys53Arg mutation that results in the more severe Grebe dysplasia.
The phenotypic severity gradient of the clinically and radiologically related acromesomelic chondrodysplasia spectrum of skeletal disorders may be due to the extent of functional impairment of the ligand-receptor pair GDF5-BMPR1B.
The acromesomelic dysplasias (ACD) constitute a rare subgroup of osteochondrodysplasias that are characterised by short stature and shortened limbs with anomalies of the hands and feet . According to the 2010 revision of the nosology and classification of genetic skeletal disorders, five ACD subgroups are recognised, namely the severe Grebe dysplasia (including Grebe (OMIM #200700) and Hunter-Thompson (OMIM #201250) types), the milder du Pan dysplasia (OMIM #228900), a Grebe dysplasia-like phenotype with genital anomalies (OMIM #609441), as well as the ACD types Maroteaux (OMIM #602875) and Osebold-Remondini (OMIM 112910) (the latter two types will not be considered further here) . While Grebe dysplasia is characterised by severe distal limb anomalies with rudimentary fingers and toes, there is mild short stature and milder limb involvement in du Pan dysplasia, including fibular hypoplasia and complex brachydactyly.
Both Grebe dysplasia and du Pan dysplasia are autosomal recessive disorders and can be caused by biallelic loss-of-function mutations of GDF5 (previously known as CDMP1), encoding the growth and differentiation factor 5 [3–6]. GDF5 belongs to the bone morphogenetic protein (BMP) family and binds to the BMP receptors BMPR1A and BMPR1B, with a preference for BMPR1B . Consistent with this ligand-receptor interaction, ACD cannot only be caused by GDF5 mutations, but in rare cases also by loss-of-function BMPR1B mutations. Indeed, Demirhan and colleagues identified a homozygous truncating mutation of BMPR1B underlying a Grebe dysplasia-like phenotype with genital anomalies (absent ovaries, hypoplastic uterus, hypergonadotropic hypogonadism and primary amenorrhea) ; and we previously reported a homozygous missense mutation (p.Cys53Arg) and a homozygous nonsense mutation (p.Trp219*) of BMPR1B in Grebe dysplasia . Functional analysis of the p.Cys53Arg missense alteration indicated that the mutant receptor had reduced localisation at the cell membrane, reduced activation by GDF5 binding and no effect on cell differentiation upon overexpression in an in vitro chondrogenesis assay, consistent with loss of function . Here we have identified and characterised a hypomorphic BMPR1B missense alteration that causes du Pan dysplasia.
Clinical investigation and molecular analyses
The clinical evaluation was performed at the Department of Medical Genetics, Hospital Erasme, Brussels, Belgium, and the Institut de Pathologie et de Génétique, Gosselies, Belgium, and research was performed following a protocol approved by the University of Ulm Ethics Committee. After obtaining written informed consent for participation in the study, including consent to report individual patient data, we sequenced the coding exons and splice sites of GDF5 (NM_000557.2) and BMPR1B (NM_001256794.1) on genomic DNA after PCR amplification. Primer sequences are available on request. PCR products were sequenced on an ABI 3730 DNA Analyzer with BigDye chemistry v3.1 (Applied Biosystems). Sequence traces were assembled, aligned, and analyzed with the Seqman software (DNASTAR Lasergene). The BMPR1B intragenic short tandem repeat (STR) markers D4S2278 and BMPR1B-STR2 (forward primer 5’-TTAAAAGTAAAGACCATATAAAGG-3’; reverse primer 5’-GGACTTAAAACTCTGCAACAAC-3’) were genotyped after PCR amplification on an ABI 3730 DNA Analyzer. Single nucleotide polymorphism (SNP) array genotyping was performed using a Genome-Wide Human SNP Array 6.0 (Affymetrix), according to the manufacturer´s instructions. Linkage analysis and haplotype reconstruction were performed using the program ALLEGRO . Haplotypes were presented graphically with HaploPainter . All data handling was performed using the graphical user interface ALOHOMORA .
The BMPR1B substitutions Arg31Cys (identified here) and Arg31His (present at low frequency in public databases; rs200035802) were modeled on the basis of the crystal structure of the human BMPR1B extracellular domain bound to human GDF5 (PDB entry 3EVS, ). The side chain of Arg31 was exchanged for either a cysteine or histidine using the ProteinDesign tool of the software Quanta2008 (MSI Accelrys) and rotamer searches were performed to identify the most likely side chain conformation. Potential alternative disulfide connectivities were analyzed by performing energy minimization runs in vacuo employing only geometrical energy terms.
Cloning of BMPR1B expression constructs
The expression constructs for HA-tagged Bmpr1b WT and Bmpr1b Cys53Arg were previously described . HA-tagged Bmpr1b WT in pSLAX13 was used as a template for in vitro mutagenesis to introduce the mutation coding for Arg31Cys (Primer forward 5´-ctcggcccaagatcctaTgttgtaaatgccaccac-3´, reverse 5´-gtggtggcatttacaacAtaggatcttgggccgag-3´) or coding for Arg31His (Primer forward 5´-cggcccaagatcctacAttgtaaatgccaccac-3´, reverse 5´-gtggtggcatttacaaTgtaggatcttgggccg-3´). Inserts were subcloned into the expression vector pCS2+ after ClaI restriction.
Luciferase reporter gene assay
NIH/3 T3 cells (ATCC) were seeded in a 96-well plate in growth medium (DMEM 4.5 g/l glucose (Lonza), 10 % FBS superior (Biochrom), 2 mM L-glutamine (Lonza)) and transfected 24 h later using Lipofectamine 2000 (Invitrogen, Life Technologies) following the manufacturer´s instructions. Cells were transfected in growth medium with the control vector pCS2+ or one of the Bmpr1b variants in pCS2+ together with the BRE luciferase reporter construct BRE-pLG3ti  and the Renilla luciferase normalization vector pRL-TK (Promega). After 18 h cells were stimulated with 2 nM human recombinant GDF5 (Biopharm) in serum-reduced medium (DMEM 4.5 g/l glucose (Lonza), 1 % FBS superior (Biochrom), 2 mM L-glutamine (Lonza)). 40 h after transfection cells were lysed in potassium phosphate buffer (9 mM potassium di-hydrogen phosphate, 91 mM di-potassium phosphate, 0.2 % Triton-X-100) and dual luciferase activity was measured as described previously  using the Mithras LB 940 (Berthold Technologies). For statistical analysis GraphPad Prism 5 (GraphPad Software, Inc.) was used.
NIH/3 T3 cells were seeded into a 24-well plate in growth medium on cover glasses (Marienfeld). 24 h later Bmpr1b expression vectors were transfected into the NIH/3 T3 cells using Lipofectamine 2000 (Invitrogen, Life Technologies) following the manufacturer´s instructions. Another 24 h later cells were incubated under serum free conditions (DMEM 4.5 g/l glucose (Lonza), 2 mM L-glutamine (Lonza)) for 1 h, subsequently fixed with 4 % paraformaldehyde in PBS and blocked in PBS containing 10 % FBS superior (Biochrom) over night. Non-permeated cells were incubated with a rabbit anti-HA antibody (H6908, Sigma-Aldrich; diluted 1:100 in PBS containing 10 % FBS superior (Biochrom)) for 1 h. After washing with PBS the cells were incubated with the secondary antibody anti-rabbit-Alexa Fluor 488 (A11008, Molecular Probes Life Technologies) and with DAPI (Invitrogen, Life Technologies) for 1 h. The cover glasses containing stained cells were mounted on microscope slides (‘SuperFrost Plus’, Menzel) using Fluoromount-G (Southern Biotech). Confocal microscopy was performed using Zeiss Axio Imager.M2 equipped with a LSM700 confocal module (Carl Zeiss, 63-fold magnification).
The proposita was a 38 year-old woman who presented to the genetics clinic because of short stature and limb anomalies. She is the 11th child (of 11) born to first cousin parents originating from Morocco (Fig. 1a). One of her sisters and two of her brothers are similarly affected; they live in a remote area of Morocco with very limited access to medical care and no detailed clinical information, skeletal radiographs or DNA were available from them. The parents were unaffected; maternal height was 160 cm and paternal height 180 cm.
The proposita had disproportionate short stature with a height of 148 cm (− 2.7 SD) and mild acromesomelic limb shortening. The fibulae were not palpable. She had short fingers and abnormal finger joints with deviations (Fig. 1b). Fingers II and V were most severely affected. Finger nails were normal. The toes were hypoplastic, in particular toes III to V, with broad and short toe nails (Fig. 1b). Radiological examination for the hands revealed proximal symphalangism and dysplastic middle phalanges (absence of the middle phalanx in the index finger, hypoplasia of the middle phalanges with narrow proximal phalangeal joints in the middle and ring fingers, and absence of the middle phalanx and long proximal phalanx in the little finger), shortening of the first metacarpal and a single long phalanx in the thumb, short proximal phalanx of the index finger and narrow carpal joint spaces (Fig. 1c). Radiographs of the feet showed similar findings, including a single misshapen phalanx in the great toe, hypoplasia of the middle phalanx in the second toe, absence of the middle phalanges in the third, fourth, and fifth toes, and calcaneocuboid fusion (Fig. 1c). The ankle joints were mal-aligned. The fibula was totally missing (Fig. 1c).
The proposita reported normal puberty with menarche at 13 years and regular menstrual cycles. Abdominal ultrasound was normal with the presence of a uterus and ovaries of normal size and shape. There was no hypergonadotropic hypogonadism with normal LH, FSH, estradiol and progesterone.
Because the clinical presentation was compatible with du Pan dysplasia we first sequenced the coding region and exon-intron boundaries of GDF5 and identified no mutation. Upon sequencing of BMPR1B we detected a homozygous c.91C>T (p.Arg31Cys) variant (Fig. 1d). No DNA from other affected and unaffected family members was available to test for co-segregation with du Pan dysplasia. Consistent with BMPR1B being located in an autozygous region, genotyping of two highly polymorphic BMPR1B intragenic STRs showed homozygosity and SNP array genotyping confirmed a large (20.2 Mb) homozygous region at the BMPR1B locus in the proposita, defined by the flanking SNPs rs2131361 and rs17278473 (data not shown). Residue Arg31 is moderately conserved in BMPR1B orthologs. It is conserved e.g., in mouse, rat, cat, dog, and platypus; it is replaced with a histidine in e.g., green monkey, panda and hedgehog, and with a glutamine in e.g., opossum and turtles. None of the >60 species included in the UCSC multiple species alignment has a cysteine at the corresponding position. While the ExAC browser (http://exac.broadinstitute.org/) lists a p.Arg31His substitution (rs200035802) with an allele frequency of 18/121.162, p.Arg31Cys was only identified in 1/121.150 alleles, compatible with p.Arg31Cys being disease causing. Consistently, PolyPhen-2  predicts p.Arg31Cys as “possibly damaging” and p.Arg31His as “benign” (both variants are predicted “not tolerated” by SIFT ). Finally, no variant predicting p.Arg31Cys was identified in exome sequences from 30 Moroccan patients with unrelated disorders (Jaber Lyahyai, personal communication) and no p.Arg31Cys homozygotes or heterozygotes were present in the 5,718 exomes included in the Institut Imagine (Paris) in-house database, which is enriched for exomes of Maghrebian individuals, and contains a very conservatively estimated number of Moroccan exomes of 33.
Structural modeling of the BMPR1B mutation p.Arg31Cys
Arg31 is located at the beginning of the BMPR1B ligand binding domain (LBD), suggesting that the p.Arg31Cys mutation may interfere with GDF5 binding. We modeled p.Arg31Cys using the structure of human BMPR1B bound to human GDF5 . Residue 31 shares no direct contact with the ligand GDF5 and is more than 10 Å remote (Fig. 2a). Thus, a direct impact on the GDF5-receptor binding by loss of direct van der Waals interactions seems unlikely. In our in silico models the exchange of Arg31 by cysteine however places an unpaired thiol group in close proximity to Cys32, Cys47, Cys53, and Cys71 (Fig. 2b) thereby possibly allowing for an alternate disulfide network. The model of such a rearranged disulfide bond network suggested that it will not necessarily destroy the BMP receptor’s three-finger toxin fold, but may introduce conformational changes in the disulfide-bond connected loops in particular of the β1β2-loop (Fig. 2c). As this loop is engaged in contacts with GDF5 and has been implicated in ligand recognition and binding (Fig. 2a)  even such smaller changes may affect ligand-receptor interaction. The effects of p.Arg31His seem subtler as those predicted for p.Arg31Cys (Fig. 2d).
Functional analyses of mutant BMPR1B receptors
We performed a luciferase reporter gene assay using the SMAD dependent BRE reporter to determine the signal transduction capacity of the BMPR1B variant Arg31Cys, which is associated with ACD du Pan type in our patient (Fig. 3a). We compared its activity to the previously described BMPR1B mutation Cys53Arg, which is associated with ACD Grebe type and to Arg31His, another allele listed in genomic databases . The WT receptor induced the luciferase expression approximately 8-fold compared to a control vector and was strongly stimulated by the presence of recombinant GDF5. BMPR1B Arg31Cys showed only moderate activity when no ligand was used for stimulation (less than 2-fold). Arg31His induced the SMAD signaling pathway approximately 5-fold, but remained significantly less active than the WT receptor. SMAD1/5/8 activation via the BMPR1B variants Arg31His and Arg31Cys was significantly increased upon ligand stimulation. However, BMPR1B Arg31Cys still could not accomplish the level of BMPR1B Arg31His or BMPR1B WT, pointing to a partial loss of function. In line with previous findings BMPR1B Cys53Arg exhibited no SMAD activity and could not be stimulated by GDF5 either. Confocal microscopy showed that all investigated BMPR1B variants are translocated to the cell surface (Fig. 3b).
We have identified the homozygous BMPR1B missense mutation p.Arg31Cys in a patient with ACD du Pan type. Previously, GDF5 mutations have been shown to constitute the major cause of ACD Grebe and du Pan types [3–6], and BMPR1B mutations had been identified in single families with Grebe syndrome and a Grebe dysplasia-like phenotype with genital anomalies [8, 9]. However, whether BMPR1B mutations, presumably mutations with milder functional effects, can also cause du Pan dysplasia was unknown. The following evidence supports a diagnosis of du Pan dysplasia in the proposita rather than another ACD or ACD-like skeletal dysplasia: normal appearing radius, ulna and femur including the corresponding joints; mal-shaped but present phalangeal bones of the hands; normal length of the second metacarpal; and milder involvement of the tarsals, metatarsals, and phalangeal bones of the toes. In Grebe and Hunter-Thompson ACD the clinical and radiological manifestations are more severe: the humerus, radius and ulna are short and bowed, the hand bones are more deformed or missing and the elbow and knee joints are affected. Thus, our results show that both Grebe and du Pan dysplasia can be caused by BMPR1B mutations in rare cases. While the number of affected individuals with ACD due to BMPR1B mutations is small, our previous  and present results indicate a possible genotype-phenotype correlation according to which BMPR1B mutations with a strong functional effect would cause Grebe syndrome, while milder mutations would result in the clinically and radiologically milder du Pan dysplasia.
Indeed, when the Arg31Cys mutation or the Arg31His variant are compared with other BMPR1B mutations, e.g., the Cys53Arg exchange , the assumed impact of a missense substitution at position 31 on structure and BMP binding seems rather moderate. The substitution of Cys53 not only disrupts the structure-stabilizing disulfide bond network, but also likely causes local unfolding in the center of the BMPR1B ligand-binding domain due to strong van der Waals overlaps from introducing a bulky, charged arginine side chain into the tightly packed hydrophobic core around Cys53. In contrast, the side chain of Arg31 faces the protein surface being surrounded by bulk solvent. Thus, exchanges against most other amino acid types will likely not alter the structure of the BMPR1B ligand-binding domain. Furthermore the residue at position 31 does not contact GDF5, thus only amino acid substitutions that (indirectly) cause a conformational rearrangement of ligand-binding loops can affect BMP binding.
A luciferase reporter gene assay indeed revealed that both BMPR1B variants, Arg31Cys and Arg31His, can still be stimulated by GDF5 but signal with diminished biological activity compared to BMPR1B WT. This suggests that the exchanges at position Arg31 have for that matter only moderate functional effect and do not completely prevent ligand binding, as it is the case for Cys53Arg.
In conclusion, we provide evidence for a causal role of a hypomorphic BMPR1B missense mutation in du Pan dysplasia, thereby broadening the clinical and molecular overlap between ACD types. This finding also has obvious implications for molecular diagnostic strategies, as BMPR1B sequencing should be considered in patients with ACD harboring no mutations in GDF5.
Krakow D, Rimoin DL. The skeletal dysplasias. Genet Med. 2010;12(6):327–41.
Warman ML, Cormier-Daire V, Hall C, Krakow D, Lachman R, LeMerrer M, et al. Nosology and classification of genetic skeletal disorders: 2010 revision. Am J Med Genet A. 2011;155A(5):943–68.
Thomas JT, Kilpatrick MW, Lin K, Erlacher L, Lembessis P, Costa T, et al. Disruption of human limb morphogenesis by a dominant negative mutation in CDMP1. Nat Genet. 1997;17(1):58–64.
Faiyaz-Ul-Haque M, Ahmad W, Zaidi SH, Haque S, Teebi AS, Ahmad M, et al. Mutation in the cartilage-derived morphogenetic protein-1 (CDMP1) gene in a kindred affected with fibular hypoplasia and complex brachydactyly (DuPan syndrome). Clin Genet. 2002;61(6):454–8.
Faiyaz-Ul-Haque M, Ahmad W, Wahab A, Haque S, Azim AC, Zaidi SH, et al. Frameshift mutation in the cartilage-derived morphogenetic protein 1 (CDMP1) gene and severe acromesomelic chondrodysplasia resembling Grebe-type chondrodysplasia. Am J Med Genet. 2002;111(1):31–7.
Stelzer C, Winterpacht A, Spranger J, Zabel B. Grebe dysplasia and the spectrum of CDMP1 mutations. Pediatr Pathol Mol Med. 2003;22(1):77–85.
Nickel J, Kotzsch A, Sebald W, Mueller TD. A single residue of GDF-5 defines binding specificity to BMP receptor IB. J Mol Biol. 2005;349(5):933–47.
Demirhan O, Turkmen S, Schwabe GC, Soyupak S, Akgul E, Tastemir D, et al. A homozygous BMPR1B mutation causes a new subtype of acromesomelic chondrodysplasia with genital anomalies. J Med Genet. 2005;42(4):314–7.
Graul-Neumann LM, Deichsel A, Wille U, Kakar N, Koll R, Bassir C, et al. Homozygous missense and nonsense mutations in BMPR1B cause acromesomelic chondrodysplasia-type Grebe. Eur J Hum Genet. 2014;22(6):726–33.
Gudbjartsson DF, Jonasson K, Frigge ML, Kong A. Allegro, a new computer program for multipoint linkage analysis. Nat Genet. 2000;25(1):12–3.
Thiele H, Nürnberg P. HaploPainter: a tool for drawing pedigrees with complex haplotypes. Bioinformatics. 2005;21(8):1730–2.
Rüschendorf F, Nürnberg P. ALOHOMORA: a tool for linkage analysis using 10 K SNP array data. Bioinformatics. 2005;21(9):2123–5.
Kotzsch A, Nickel J, Seher A, Sebald W, Muller TD. Crystal structure analysis reveals a spring-loaded latch as molecular mechanism for GDF-5-type I receptor specificity. EMBO J. 2009;28(7):937–47.
Korchynskyi O, ten Dijke P. Identification and functional characterization of distinct critically important bone morphogenetic protein-specific response elements in the Id1 promoter. J Biol Chem. 2002;277(7):4883–91.
Hampf M, Gossen M. A protocol for combined Photinus and Renilla luciferase quantification compatible with protein assays. Anal Biochem. 2006;356(1):94–9.
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7(4):248–9.
Kumar P, Henikoff S, Ng PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009;4(7):1073–81.
We thank the proposita for her participation in the study. We thank Ingelore Bäßmann and Mareen Schmidt-von Kegler for technical assistance; and Jaber Lyahyai as well as the genomics platform and researchers at Institut Imagine Paris for sharing data on mutation frequencies in Moroccan individuals.
The authors declare that they have no competing interests.
KS carried out cloning of BMPR1B constructs, luciferase reporter gene assays and confocal microscopy. JD and DH performed the clinical and radiological characterization. NK, BB, CTG and GB carried out the molecular genetic studies. TM performed the structural analyses. KS, PS and GB drafted the manuscript. JD, PS and GB conceived the study. All authors read and approved the final manuscript.
Katja Stange and Julie Désir contributed equally to this work.
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 https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Stange, K., Désir, J., Kakar, N. et al. A hypomorphic BMPR1B mutation causes du Pan acromesomelic dysplasia. Orphanet J Rare Dis 10, 84 (2015). https://doi.org/10.1186/s13023-015-0299-5
- Acromesomelic dysplasias
- Grebe dysplasia
- du Pan dysplasia