- Open Access
Identification of RUNX2 variants associated with cleidocranial dysplasia
Hereditas volume 156, Article number: 31 (2019)
Cleidocranial dysplasia (CCD) is a rare autosomal dominant disorder mainly characterized by hypoplastic or absent clavicles, delayed closure of the fontanelles, multiple dental abnormalities, and short stature. Runt-related transcription factor 2 (RUNX2) gene variants can cause CCD, but are not identified in all CCD patients.
In this study, we detected genetic variants in seven unrelated children with CCD by targeted high-throughput DNA sequencing or Sanger sequencing.
All patients carried a RUNX2 variant, totally including three novel pathogenic variants (c.722_725delTGTT, p.Leu241Serfs*8; c.231_232delTG, Ala78Glyfs*82; c.909C > G, p.Tyr303*), three reported pathogenic variants (c.577C > T, p.Arg193*; c.574G > A, p.Gly192Arg; c.673 C > T, p.Arg225Trp), one likely pathogenic variant (c.668G > T, p.Gly223Val). The analysis of the variant source showed that all variants were de novo except the two variants (c.909C > G, p.Tyr303*; c.668G > T, p.Gly223Val) inherited from the patient’s father and mother with CCD respectively. Further bioinformatics analysis indicated that these variants could influence the structure of RUNX2 protein by changing the number of H-bonds or amino acids. The experimental result showed that the Gly223Val mutation made RUNX2 protein unable to quantitatively accumulate in the nucleus.
The present study expands the pathogenic variant spectrum of RUNX2 gene, which will contribute to the diagnosis of CCD and better genetic counseling in the future.
Cleidocranial dysplasia (CCD; OMIM #119600) is a rare autosomal dominant disorder mainly characterised by hypoplastic or absent clavicles, delayed closure of fontanelles, multiple dental abnormalities, and short stature [1,2,3]. Variants in runt-related transcription factor 2 (RUNX2) gene (OMIM *600211) can result in haploinsufficiency of the protein and have been related to CCD [1, 2]. The RUNX2 gene is located on chromosome 6p21.1 and encodes a transcription factor with a highly conserved Runt domain [4, 5]. The Runt domain is responsible for binding to a specific DNA motif (TGT/CGGT sequence) in the promoter region of its target genes and heterodimerization with CBFB (core-binding factor subunit beta) [6,7,8]. The former participates in regulating the transcription of multiple genes. The latter increases the DNA-binding affinity as well as protects and stabilizes RUNX2 against proteolytic degradation. The N-terminal side of the Runt domain links a Q/A region consisting of 23 consecutive glutamine residues followed by 17 alanine residues, which acts as a second transactivation domain . The C-terminal side of the Runt domain links a PST (proline/serine/threonine)-rich region, which contains the phosphorylation sites and represents the third transactivation domain [9, 10]. The last five amino acids (VWRPY) of RUNX2 protein compose a conserved motif in all runt proteins, and functions as a transcriptional repression domain [9, 11].
RUNX2 is essential for osteoblastic differentiation and skeletal morphogenesis. In mouse models, the homozygous mutation of RUNX2 gene blocked both intramembranous and endochondral ossification and resulted in a complete lack of bone formation . The heterozygous mutation (RUNX2+/−) caused a similar phenotype to that of human CCD . To date, 184 publicly available mutations in RUNX2 gene have been deposited in the Human Gene Mutation Database (HGMD, www.hgmd.cf.ac.uk). Most of these mutations were missense and clustered in Runt domain. Additionally, nonsense mutations, insertions or deletions are also observed in the RUNX2 gene, which are predominant within the Q/A domain or the PST domain. Although many mutations in the RUNX2 gene have been identified in familial and sporadic cases, novel mutation is still reported recently, suggesting that mutational screening on RUNX2 gene is far from saturation [14,15,16,17,18,19].
In the present study, we conducted genetic evaluation for a cohort of seven Chinese children with CCD by targeted high-throughput DNA sequencing or Sanger sequencing, and found seven different variants in RUNX2 gene, including six pathogenic variants and one likely pathogenic variant. These results will contribute to the diagnosis of CCD and better genetic counseling in the future.
Material and methods
Genomic DNA extraction and genetic testing
A total of seven unrelated children with CCD ranging in age from 1 month to 12 years were enrolled for genetic evaluation (Table 1). Genomic DNA of probands and their family members was extracted from peripheral blood leukocytes using Lab-Aid Nucleic Acid Isolation Kit (Zeesan, China), according to the manufacturer’s instructions.
Among these CCD patients, five patients were firstly detected by targeted high-throughput DNA sequencing, two patients directly by Sanger sequencing (Table 1). For targeted high-throughput DNA sequencing, the preparation of sequencing library was completed using Agilent Inherited Disease panel, Agilent Focused exome panel or xGen Exome research panel v1.0 (Integrated DNA Technologies, Coralville, Iowa). Sequencing was performed on the Illumina HiSeq 2500 or 4000 (Illumina, San Diego, CA), according to the manufacturer’s instructions. Burrows-Wheeler Aligner (BWA, version 0.7.10) was used to mapping reads to the human reference genome (GRCh37/hg19). Base calling, QC analysis and coverage analysis were performed with Picard tools-1.124 and GATK software. Variants were annotated using SnpEff version 4.2. Subsequently, the following variants were filtered out: (i) variants with > 1% frequency in the population variant databases including 1000 Genomes Project, Exome Variant Server (EVS) and Exome Aggregation Consortium (ExAC) or > 5% frequency in our inhouse database (based on 150 exome datasets), (ii) intergenic and 3′/5′ untranslated region variants, none splice-related intronic and synonymous variants.
For Sanger sequencing, all exons of the RUNX2 gene in these probands were amplified by PCR reaction. DNA sequence variants were identified by Mutation Surveyor V4.0.5 software with reference sequences (NG_008020.1).
MutationTaster (http://www.mutationtaster.org), SIFT (http://sift.jcvi.org), and PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2/) were used to assess pathogenic potential of the variants [20,21,22]. Combined with clinical manifestation and modes of inheritance, candidate variants were validated by Sanger sequencing for all family members, and classified according to standards and guidelines of the American College of Medical Genetics and Genomics (ACMG). For the putative pathogenic or likely pathogenic variants, SWISS-MODEL (https://swissmodel.expasy.org) and Swiss-PdbViewer 4.1 software (http://spdbv.vital-it.ch/) were used to analyze the effect of these variants on protein structure [23, 24].
Subcellular localization of the RUNX2 mutant protein
The cDNA of wide-type RUNX2 gene was synthesized by Sangon Biotech (Shanghai) Co., Ltd., and amplified by PCR. The forward primer was 5′-GACACAGATCTCGAGATGGCATCAAACAGCCTCTTCAGC-3′ and the reverse primer was 5′-GTGTCGTCGACTGATATGGTCGCCAAACAGATTCA-3′. The PCR fragment was subcloned into pEGFP-N1 vector with the XhoI and SaII restriction sites. The RUNX2 668G > T (Gly223Val) mutation was introduced into pEGFP-N1 vector with wide-type RUNX2 cDNA by site-directed mutagenesis. The mutant primers were 5′-GCCTTCTGGGTTCCCGAGGTACATCTACTGTAACTTT AAT-3′, and 5′-ATTAAAGTTACAGTAGATGTACCTCGGGAACCCAGAAGGC-3′. All recombinant vectors were fully sequenced to exclude any additional mutations. The empty vector acting as a negative control (NC) and pEGFP-N1 vectors bearing wild type (WT) and mutant (Mut) were transfected into U2OS cells by using lipofectamine 2000 (Invitrogen). The cells were visualized and photographed (magnification 10X and 40X) with a fluorescent microscope (Olympus IX73, Japan).
Clinical features of CCD children
All children underwent a clinical evaluation and were diagnosed as CCD by an experienced pediatrician. The clinical features of these patients including two female and five male patients were summarized in Table 2. Besides the clavicle and skull dysplasia, short stature, scoliosis, enamel hypoplasia, delayed eruption of deciduous teeth, low nasal bridge, delayed mineralization of pubic bone, broad femoral head with short femoral neck, hypoplastic iliac wing, syringomyelia and special faces were also observed in CCD children. Furthermore, hypertelorism was observed in all CCD children, except Family_A_II1. Supernumerary teeth, retention cysts and long second metacarpal were observed in all CCD children, except Family_A_II1 and Family_B_II1.
All patients carried a RUNX2 variant, totally including four novel variants and three reported variants (Figs. 1, 2 and Table 3). Among the seven variants, there were two pathogenic missense variants (c.574G > A, p.Gly192Arg; c.673 C > T, p.Arg225Trp), one likely pathogenic missense variant (c.668G > T, p.Gly223Val), two pathogenic frameshift variants (c.722_725delTGTT, p.Leu241Serfs*8; c.231_232delTG, Ala78Glyfs*82), and two pathogenic stop-gain variants (c.577C > T, p.Arg193*; c.909C > G, p.Tyr303*). The analysis of the variant source showed that all variants were de novo except the two variants (c.909C > G, p.Tyr303*; c.668G > T, p.Gly223Val). The former variant was inherited from the patient’s father with CCD, who carried a de novo heterozygous RUNX2 variant (c.909C > G, p.Tyr303*). The latter variant was inherited from the patient’s mother with CCD, who carried a maternal inherited and heterozygous RUNX2 variant (c.668G > T, p.Gly223Val).
The effect of the RUNX2 variants on protein structure
Among these variants, there were three variants changing the number of H-bonds in RUNX2 protein, including two variants increasing H-bonds (c.574G > A, p.Gly192Arg; c.668G > T, p.Gly223Val) and one variant decreasing H-bonds (c.673 C > T, p.Arg225Trp). In addition, there were four variants (c.722_725delTGTT, p.Leu241Serfs*8; c.231_232delTG, Ala78Glyfs*82; c.577C > T, p.Arg193*; c.909C > G, p.Tyr303*) decreasing the number of amino acids in RUNX2 protein.
Subcellular localization of the RUNX2 mutant protein
To further explore the function of the missense mutation (c.668G > T, p.Gly223Val) not reported, the wide-type and mutant RUNX2 proteins binding green fluorescent protein (GFP) were constructed and transiently transfected into human osteosarcoma U2OS. The result showed that the Gly223Val mutation could affect the subcellular distribution of RUNX2 protein and made RUNX2 protein unable to quantitatively accumulate in the nucleus (Fig. 3).
CCD is a skeletal dysplasia that represents a continuum of clinical findings ranging from classical CCD (dental abnormalities, hypoplastic or aplastic clavicles, and delayed closure of the cranial sutures) to mild CCD to isolated dental anomalies without other skeletal features. To date, no formal clinical diagnostic criteria for CCD have been established. Due to CCD inherited in an autosomal dominant manner, each child of an individual with CCD has a 50% chance of inheriting the pathogenic variant. If the pathogenic variant in the family is known, prenatal diagnosis for pregnancies at increased risk will be possible. Many kinds of molecular testing approaches, including single-gene testing, karyotype analysis and a multigene panel, can be currently used to detect the variants leading to CCD. For single-gene testing, sequence analysis of RUNX2 gene is performed first and followed by gene-targeted deletion/duplication analysis if no pathogenic variant is identified. For karyotype analysis, if RUNX2 testing is not diagnostic and strong suspicion persists in an individual with CCD features who also has multiple congenital anomalies and/or developmental delay, a karyotype analysis may be considered to evaluate complex chromosome rearrangements or translocations that involve RUNX2 locus but do not result in RUNX2 copy number changes [25, 26]. In addition, a multigene panel that includes RUNX2 and other genes of interest may also be considered.
In the present study, we utilized targeted high-throughput DNA sequencing or Sanger sequencing (single-gene testing) techniques to analyze genetic variants in seven CDD children, and found seven different variants in RUNX2 gene, including four novel variants (c.722_725delTGTT, p.Leu241Serfs*8; c.231_232delTG, Ala78Glyfs*82; c.909C > G, p.Tyr303*; c.668G > T, p.Gly223Val) and three reported variants (c.577C > T, p.Arg193*; c.574G > A, p.Gly192Arg; c.673 C > T, p.Arg225Trp) [27,28,29], which were all located in the transactivation region (Fig. 2). The bioinformatics analysis indicated that these variants were disease-causing, damaging and/or probably damaging variants. According to ACMG, six variants (c.574G > A, p.Gly192Arg; c.673 C > T, p.Arg225Trp; c.577C > T, p.Arg193*; c.722_725delTGTT, p.Leu241Serfs*8; c.231_232delTG, Ala78Glyfs*82; c.909C > G, p.Tyr303*) were classified as pathogenic variants, and one variant (c.668G > T, p.Gly223Val) as likely pathogenic variant. In addition, all variants were de novo except the following two variants: c.909C > G, p.Tyr303* and c.668G > T, p.Gly223Val. Thereinto the former variant (c.909C > G, p.Tyr303*) was inherited from the patient’s father, who is also a CCD patient carried a de novo heterozygous RUNX2 variant. The clinical features of the father included short stature and CCD, which were very similar to those of his 3-year-old son. The latter variant (c.668G > T, p.Gly223Val) was inherited from the patient’s mother with CCD, who carried a maternal inherited and heterozygous RUNX2 variant. Both of them also showed similar clinical phenotypes, such as short stature and CCD. By summarizing RUNX2 variants in HGMD and the current study, we found nine variant types, such as missense/nonsense, splicing, small deletions/insertions, gross insertions/duplications. Thereinto missense/nonsense variant was the most common variant type of RUNX2 gene (Table 4). A single amino acid (Gly) substitution at position 332 in RUNX2 protein was found not only in our lab (c.668G > T, p.Gly223Val), but also in Ott’ s study (c.667G > A, p.Gly223Arg) . In addition, protein structure prediction showed that these variants could change the number of H-bonds or amino acids in RUNX2 protein (Fig. 4), suggesting that these variants played an important role in regulating the effective structure and function of RUNX2 protein. The experimental result showed that Gly223Val mutation, located in nuclear localization sequence (NLS) [29, 30], could affect the subcellular distribution of RUNX2 protein. The mutation made RUNX2 protein unable to quantitatively accumulate in the nucleus.
In conclusion, the present study reveals some novel genetic causes of CDD, which not only expands the pathogenic variant spectrum of RUNX2 gene but also will contribute to the diagnosis of CCD and better genetic counseling in the future.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
American College of Medical Genetics and Genomics
Exome Variant Server
Exome Aggregation Consortium
Runt-related transcription factor 2
Ott CE, Leschik G, Trotier F, Brueton L, Brunner HG, Brussel W, et al. Deletions of the RUNX2 gene are present in about 10% of individuals with cleidocranial dysplasia. Hum Mutat. 2010;31(8):E1587–93.
Otto F, Kanegane H, Mundlos S. Mutations in the RUNX2 gene in patients with cleidocranial dysplasia. Hum Mutat. 2002;19(3):209–16.
Otto F, Thornell AP, Crompton T, Denzel A, Gilmour KC, Rosewell IR, et al. Cbfa1, a candidate gene for cleidocranial dysplasia syndrome, is essential for osteoblast differentiation and bone development. Cell. 1997;89(5):765–71.
Levanon D, Negreanu V, Bernstein Y, Bar-Am I, Avivi L, Groner Y. AML1, AML2, and AML3, the human members of the runt domain gene-family: cDNA structure, expression, and chromosomal localization. Genomics. 1994;23(2):425–32.
Yoshida T, Kanegane H, Osato M, Yanagida M, Miyawaki T, Ito Y, et al. Functional analysis of RUNX2 mutations in Japanese patients with cleidocranial dysplasia demonstrates novel genotype-phenotype correlations. Am J Hum Genet. 2002;71(4):724–38.
Ogawa E, Maruyama M, Kagoshima H, Inuzuka M, Lu J, Satake M, et al. PEBP2/PEA2 represents a family of transcription factors homologous to the products of the Drosophila runt gene and the human AML1 gene. Proc Natl Acad Sci U S A. 1993;90(14):6859–63.
Kagoshima H, Shigesada K, Satake M, Ito Y, Miyoshi H, Ohki M, et al. The runt domain identifies a new family of heteromeric transcriptional regulators. Trends Genet. 1993;9(10):338–41.
Ducy P, Zhang R, Geoffroy V, Ridall AL, Karsenty G. Osf2/Cbfa1: a transcriptional activator of osteoblast differentiation. Cell. 1997;89(5):747–54.
Thirunavukkarasu K, Mahajan M, McLarren KW, Stifani S, Karsenty G. Two domains unique to osteoblast-specific transcription factor Osf2/Cbfa1 contribute to its transactivation function and its inability to heterodimerize with Cbfbeta. Mol Cell Biol. 1998;18(7):4197–208.
Pande S, Browne G, Padmanabhan S, Zaidi SK, Lian JB, van Wijnen AJ, et al. Oncogenic cooperation between PI3K/Akt signaling and transcription factor Runx2 promotes the invasive properties of metastatic breast cancer cells. J Cell Physiol. 2013;228(8):1784–92.
Levanon D, Goldstein RE, Bernstein Y, Tang H, Goldenberg D, Stifani S, et al. Transcriptional repression by AML1 and LEF-1 is mediated by the TLE/Groucho corepressors. Proc Natl Acad Sci U S A. 1998;95(20):11590–5.
Komori T, Yagi H, Nomura S, Yamaguchi A, Sasaki K, Deguchi K, et al. Targeted disruption of Cbfa1 results in a complete lack of bone formation owing to maturational arrest of osteoblasts. Cell. 1997;89(5):755–64.
Komori T. Runx2, an inducer of osteoblast and chondrocyte differentiation. Histochem Cell Biol. 2018;149(4):313–23.
Peng YJ, Chen QY, Fu DJ, Liu ZM, Mao TT, Li J, et al. A novel gene mutation of Runx2 in cleidocranial dysplasia. J Huazhong Univ Sci Technolog Med Sci. 2017;37(5):772–6.
Zhang T, Wu J, Zhao X, Hou F, Ma T, Wang H, et al. Whole-exome sequencing identification of a novel splicing mutation of RUNX2 in a Chinese family with cleidocranial dysplasia. Arch Oral Biol. 2019;100:49–56.
Çamtosun E, Akıncı A, Demiral E, Tekedereli İ, Sığırcı A. A cleidocranial dysplasia case with a novel mutation and growth velocity gain with growth hormone treatment. J Clin Res Pediatr Endocrinol. 2018;11(3):301–5. [Epub ahead of print].
Ma D, Wang X, Guo J, Zhang J, Cai T. Identification of a novel mutation of RUNX2 in a family with supernumerary teeth and craniofacial dysplasia by whole-exome sequencing: a case report and literature review. Medicine (Baltimore). 2018;97(32):e11328.
Hashmi JA, Almatrafi A, Latif M, Nasir A, Basit S. An 18 bps in-frame deletion mutation in RUNX2 gene is a population polymorphism rather than a pathogenic variant. Eur J Med Genet. 2019;62(2):124–8.
Zeng L, Wei J, Zhao N, Sun S, Wang Y, Feng H. A novel 18-bp in-frame deletion mutation in RUNX2 causes cleidocranial dysplasia. Arch Oral Biol. 2018;96:243–8.
Schwarz JM, Cooper DN, Schuelke M, Seelow D. MutationTaster2: mutation prediction for the deep-sequencing age. Nat Methods. 2014;11(4):361–2.
Choi Y, Chan AP. PROVEAN web server: a tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics. 2015;31(16):2745–7.
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.
Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, Gumienny R, et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 2018;46(W1):W296–303.
Guex N, Peitsch MC. SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis. 1997;18(15):2714–23.
Purandare SM, Mendoza-Londono R, Yatsenko SA, Napierala D, Scott DA, Sibai T, et al. De novo three-way chromosome translocation 46,XY,t (4;6;21)(p16;p21.1; q21) in a male with cleidocranial dysplasia. Am J Med Genet A. 2008;146A(4):453–8.
Northup JK, Matalon R, Lockhart LH, Hawkins JC, Velagaleti GV. A complex chromosome rearrangement, der (6) ins (6)(p21.1q25.3q27)inv (6)(p25.3q27), in a child with cleidocranial dysplasia. Eur J Med Genet. 2011;54(4):e394–8.
Zhou G, Chen Y, Zhou L, Thirunavukkarasu K, Hecht J, Chitayat D, et al. CBFA1 mutation analysis and functional correlation with phenotypic variability in cleidocranial dysplasia. Hum Mol Genet. 1999;8(12):2311–6.
Puppin C, Pellizzari L, Fabbro D, Fogolari F, Tell G, Tessa A, et al. Functional analysis of a novel RUNX2 missense mutation found in a family with cleidocranial dysplasia. J Hum Genet. 2005;50(12):679–83.
Quack I, Vonderstrass B, Stock M, Aylsworth AS, Becker A, Brueton L, et al. Mutation analysis of core binding factor A1 in patients with cleidocranial dysplasia. Am J Hum Genet. 1999;65(5):1268–78.
Kanno T, Kanno Y, Chen LF, Ogawa E, Kim WY, Ito Y. Intrinsic transcriptional activation-inhibition domains of the polyomavirus enhancer binding protein 2/core binding factor alpha subunit revealed in the presence of the beta subunit. Mol Cell Biol. 1998;18(5):2444–54.
We thank the patients and their families for participating in our study.
This work is supported by grants from the “Youth Research Project of the Shanghai Municipal Health and Family Planning Commission” (No. 20184Y0348, to GXR); the “National Natural Science Foundation of China” (No.81800780 to GXR; No.81670812 and No.81873671, to YYG); the “Jiaotong University Cross Biomedical Engineering” (No.YG2017MS72, to YYG); the “Shanghai Municipal Commission of Health and Family Planning” (No.201740192, to YYG); the “Shanghai Shen Kang Hospital Development Center new frontier technology joint project” (No.SHDC12017109, to YYG).
Ethics approval and consent to participate
This study was carried out in accordance with the Code of Ethics of the World Medical Association (Declaration of Helsinki) and was certified by the Ethics Committee of Xinhua Hospital affiliated with the Shanghai Jiaotong University School of Medicine. Written informed consent was obtained from the patient’s parents.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.