- 1National Engineering Laboratory for Digital and Material Technology of Stomatology, Department of Preventive Dentistry, National Center of Stomatology, National Clinical Research Center for Oral Diseases, Peking University School and Hospital of Stomatology, Beijing, China
- 2Department of Stomatology, Xuanwu Hospital Capital Medical University, Beijing, China
- 3National Engineering Laboratory for Digital and Material Technology of Stomatology, Central Laboratory, Department of Oral and Maxillofacial Surgery, National Center of Stomatology, National Clinical Research Center for Oral Diseases, Peking University School and Hospital of Stomatology, Beijing, China
Cleidocranial dysplasia (CCD) is an autosomal dominant inheritable skeletal disorder characterized by cranial dysplasia, clavicle hypoplasia, and dental abnormalities. Mutations involving Runt-related transcription factor 2 (RUNX2) are currently the only known molecular etiology for CCD but are not identified in all CCD patients. No RUNX2 abnormality can be detected in about 20–30% of patients, and the molecular cause remains unknown. The present study includes a family case with typical features of CCD. RUNX2 mutation was first screened by sequencing analysis, and no mutation was detected. Copy number alterations of the RUNX2 gene were then measured by quantitative PCR and multiplex ligation-dependent probe amplification (MLPA). No copy number variation in RUNX2 could be detected. We performed whole-exome sequencing (WES) to identify the underlying genetic mutations. Unexpectedly, no abnormalities could be detected in genes related to the RUNX2 signaling pathway. Therefore, it was supposed that other new unknown gene variations might contribute to the CCD phenotype. We focused on Immunoglobulin superfamily member 10 (IGSF10), a gene related to bone development. An IGSF10 frameshift mutation (c.6001_6002delCT, p.Leu2001Valfs*24) was detected by WES. Sanger sequencing verified that this mutation was only detected in the patient and her affected mother but not in her unaffected father. Bioinformatics studies demonstrated that this mutation could change the 3D structure of the IGSF10 protein and severely damage its function. In addition, alkaline phosphatase (ALP) activity and the ability to form mineralized nodules were inhibited by IGSF10 knockdown compared with normal controls. The expression of bone sialoprotein (BSP) was significantly reduced by IGSF10 knockdown, but not that of other osteogenic markers. Our results provide new genetic evidence that IGSF10 mutation might contribute to CCD.
Introduction
Cleidocranial dysplasia (CCD; MIM 119600) is an autosomal dominant skeletal disorder that is, characterized by delayed closure of the fontanels, hypoplastic or aplastic clavicles, and dental anomalies (Mundlos, 1999). Mid-face hypoplasia, hand abnormalities, short stature, and other skeletal anomalies are also common (Mundlos, 1999). Highly variable expressivity generates phenotypic heterogeneity among CCD patients, even within the same family (Chitayat et al., 1992), which adds to the difficulty of diagnosing CCD.
The causative gene of CCD has been identified as Runt-related transcription factor 2 (RUNX2, also known as CBFA1) (Mundlos et al., 1997), which is essential for osteoblast differentiation and skeletal development (Komori, 2020). Numerous mutations in RUNX2 have been identified in patients with CCD, the majority of which were missense, nonsense, and frameshift mutations (Zhou et al., 1999; Otto et al., 2002). Chromosomal abnormalities, splicing mutations, and intragenic deletions/duplications were also identified in some patients (Lee et al., 2008; Purandare et al., 2008; Ott et al., 2010; Northup et al., 2011; Zhang et al., 2019). However, RUNX2 mutations are only detected in approximately two-thirds of patients with a CCD phenotype (Mundlos et al., 1997; Baumert et al., 2005), and numerous CCD patients have been identified who had no detectable mutations in RUNX2 by sequencing. There are several possible reasons for this observation. First, genetic heterogeneity may exist, such as mutation in RUNX2 gene regulatory elements or in genes encoding proteins that interact with RUNX2, including core-binding factor subunit beta (CBFB) (Goto et al., 2004; Napierala et al., 2005; Khan et al., 2006). Furthermore, other genetic diseases can show a similar phenotype to that of CCD, such as parietal foramina with cleidocranial dysplasia, which results from mutations in Drosophila muscle segment homeobox gene homologue 2 (MSX2) (Garcia-Minaur et al., 2003; Ott et al., 2012).
Among patients with normal RUNX2 sequencing analysis results, detection of RUNX2 copy number changes identified a microdeletion/duplication frequency of 26% (Ott et al., 2010). However, for the remaining patients the molecular cause is unknown. This indicates that other unreported mechanisms may contribute to CCD.
In this study, we present a familial case with typical features of CCD. However, sequencing analysis did not reveal any mutations in RUNX2, and no microdeletion/duplication in RUNX2 was detected by copy number analysis. Further exploration by whole-exome sequencing (WES) and bioinformatics analyses were performed and a new CCD candidate gene was identified, immunoglobulin superfamily member 10 (IGSF10), a member of the immunoglobulin superfamily. A role for IGSF10 in CCD has not been previously reported, though pathogenic IGSF10 mutations have been detected in self-limited delayed puberty (DP) and congenital hypogonadotropic hypogonadism (CHH) (Howard et al., 2016; Amato et al., 2019). The primary manifestation of self-limited DP and CHH is delayed puberty, absent in the present patients (Amato et al., 2019; Raivio and Miettinen, 2019). After identification of the mutation, we preliminarily explored the role of IGSF10 in osteogenic differentiation.
Materials and Methods
Participants
Two CCD patients in a Chinese family were recruited according to the criteria for the clinical diagnosis of CCD (Mundlos, 1999), together with 50 control subjects. Informed consent was obtained from the guardians of the proband and other participants. All study protocols were approved by the Ethical Committee of Peking University Health Science Center (Approval Number: IRB00001052-07100).
Mutation Analysis
Peripheral blood samples were collected from the participants and genomic DNA was extracted using a TIANamp Blood DNA mini kit (TIANGEN, Beijing, China) following the manufacturer’s instructions. RUNX2 mutation was analyzed in the proband and her parents. The exons and exon-intron boundaries of RUNX2 were amplified by polymerase chain reaction (PCR) using intron–exon specific primers as described previously (Quack et al., 1999). We also used other previously reported primers (Napierala et al., 2005) and designed others using software (Primer Premier 5 software, Premier Biosoft, Palo Alto, CA, United States) to amplify the 2.0 kb region upstream of the transcription start site of RUNX2 transcript variant 1 (NM_001024630.4). To confirm IGSF10 mutation detected in WES as pathogenic, we sequenced this variant in the proband, her parents and 50 unrelated control people. IGSF10-specific primers were also designed with the same software (Primer Premier 5 software, Premier Biosoft). Primer sequences are listed in Supplementary Table S1. Genomic DNA was amplified with Premix Taq™ (Takara, Shiga, and Japan). In brief, PCR reactions were performed with a DNA Engine PTC-200 (Bio-Rad Laboratories, Hercules, CA, United States) using a program described elsewhere (Zhang et al., 2017). The amplification products were bi-directionally sequenced using an ABI 3730 XL automatic sequencer (Applied Biosystems, Foster City, CA, United States). DNA sequences were analyzed using NCBI databases and the BLASTN (BLAST nucleotide) program (http://blast.ncbi.nlm.nih.gov/).
Real-Time PCR
RUNX2 copy number was determined by real-time PCR using a SYBR Green PCR kit (Roche Applied Science, Indianapolis, IN, United States). We detected the copy number of the coding region and the flanking untranslated regions (UTRs) of RUNX2 using the primers and method described previously (Ott et al., 2010). The primers were listed in Supplementary Table S2. qPCR was performed in a total volume of 12 µl containing 10 ng genomic DNA, 6 µl FastStart Universal SYBR Green PCR Master (Roche Applied Science), and 1 µl primers (0.2 µmol each). Samples were run on the ABI 7,500 Real-Time PCR System (Applied Biosystems, Foster, CA, United States) in triplicate in separate reactions to determine the variation in copy number between CCD patients and controls. Two controls were analyzed. The copy number of RUNX2 was normalized against those of the Albumin (ALB) gene and were calculated using the 2−ΔΔCt method.
Real-time PCR was also used to evaluate the expression of osteogenic-associated genes. MC3T3-E1 cells are osteoblast precursor cells derived from mouse, we used this cell line to explore the role of IGSF10 in osteogenesis. After MC3T3-E1 cells were cultured in osteogenic medium (OM) for 7 days, total RNA was extracted using Trizol reagent (Invitrogen, Carlsbad, CA, United States) following the manufacturer’s instructions. Next, 500 ng of total RNA were used in a 10 µl reaction volume to synthesize cDNA using a PrimeScript RT reagent kit (Takara). The resultant cDNAs were amplified using specific sets of primers for Igsf10, Runx2, bone sialoprotein (Bsp), alkaline phosphatase (Alp), osteocalcin (Ocn), and osterix (Osx). The housekeeping gene glyceraldehyde-3-phosphate dehydrogenase (Gapdh) was used to normalize RNA expression levels. The primers used are listed in Supplementary Table S3. Primer specificity was confirmed with amplicon dissociation curves. The RT-qPCR was performed using 1 µl cDNA product in a 20 µl reaction volume with a SYBR Green PCR kit (Roche Applied Science) according to the manufacturer’s protocol, and the results were assessed with 7,500 Software 2.0.1 (Applied Biosystems). Relative mRNA expression was calculated using the 2−ΔΔCt method with the formula:
MLPA and Karyotype Analysis
Multiplex ligation-dependent probe amplification (MLPA) was performed with the SALSA MLPA Probemix P080-C2 kit according to the manufacturer’s instructions (MRC-Holland, Amsterdam, Netherlands) by the Department of Medical Genetics, Peking University Health Science Center (Beijing, China). All the peak areas were normalized, and a ratio between 0.8 and 1.2 was considered normal. A heterozygous duplication was identified with a ratio between 1.3 and 1.7, and a heterozygous deletion between 0.3 and 0.7.
Karyotype analysis was performed by Giemsa-Trypsin (GTG) banding using standard procedures by the Department of Medical Genetics, Peking University Health Science Center (Beijing, China). Twenty G-banded metaphase cells were examined.
WES and Data Analysis
Whole-exome sequencing (WES) was performed using genomic DNA from the proband and her parents. DNA libraries were prepared using the Agilent SureSelect system (Agilent Technologies, Santa Clara, CA, United States) and sequenced on the Illumina HiSeq 2000 platform (Illumina, San Diego, CA, United States) to generate 100 bp paired-end reads. Bioinformatics analyses included production of clean data, alignment to human reference genome hg19 using Burrows–Wheeler Aligner (BWA) (Li and Durbin, 2009), variant calling, and annotation using Samtools, Genome Analysis Tool Kit, and Annovar (Hintzsche et al., 2016). We then used in-house AnnoDB software to annotate the confident variant results. Quality control was present throughout the whole procedure. Next, steps were performed to identify potential causal variants. First, the variants were segregated by trait (variants present in affected individuals but not in unaffected family members). Second, variants were further filtered by minor allele frequency (MAF <0.5%) in the context of heterozygous genotype in the 1,000 Genomes database (1000G_ALL and 1000G_EAS) and the NHLBI exome variant server (ESP_ALL). Third, deleterious variants were picked based on SIFT, Polyphen-2, and variant impact. Then, the functions of those genes retained after filtering were reviewed, and genes related to bone development, especially those in the RUNX2 pathway, were selected for further study.
Prediction of Damaging Effects
Multiple sequence alignment for IGSF10 was performed using ClustalW in MEGA-X software to analyze the conservation of affected amino acids (Kumar et al., 2018). To analyze the effect of the mutation on the molecular structure of IGSF10, the three-dimensional (3D) structures of wild-type and mutant IGSF10 were predicted in silico using I-TASSER (https://zhanglab.ccmb.med.umich.edu/I-TASSER/) (Yang et al., 2015).
Cell Culture
MC3T3-E1, a pre-osteoblast murine cell line, was purchased from the Cell Bank of Chinese Academy of Sciences and cultured in a proliferation medium containing α-MEM (Gibco, Grand Island, NY, United States), 10% fetal bovine serum (Gibco), 100 U/ml penicillin (Gibco) and 100 μg/ml streptomycin (Gibco) in an incubator at 37°C with 5% CO2. For osteogenic induction, MC3T3-E1 cells were cultured in proliferation medium plus 100 nM dexamethasone (Sigma-Aldrich, St Louis, MO, United States), 10 mM β-glycerophosphate (Sigma-Aldrich), and 50 μg/ml l-ascorbic acid (Sigma-Aldrich).
Transfection
The pCMV5-Runx2 plasmid expressing full-length wild-type RUNX2 was provided by the School of Dentistry, University of Michigan (Ann Arbor, MI, United States). Plasmid pcDNA3-CMF608FL (Segev et al., 2004), containing a full-length open reading frame for expression of Flag-tagged rat IGSF10 protein, was provided by Quark Pharmaceuticals Inc. (Newark, CA, United States/Ness Ziona, Israel). The sequences of wild-type Runx2 and Igsf10 were confirmed by DNA sequencing. Empty pCMV5 and pcDNA3.1 vectors served as transfection controls. Three small interfering RNAs (siRNAs) for Igsf10 and Runx2 were purchased from Ribo Bio Technology (Guangzhou, China). The most efficient siRNA for each gene was selected for the subsequent experiments: si-Igsf10, 5′-GCACCTTCCTGATTTCAAA-3′; and si-Runx2, 5′-GCACGCTATTAAATCCAAA-3′. Scrambled siRNA was used as a negative control.
For transfection, MC3T3-E1 cells were seeded at a density of 1×105 cells/well in six-well plates. When cells reached 80% confluency for plasmids and 30–50% confluency for siRNAs, transfection was performed with Lipofectamine 3,000 (Invitrogen, Shanghai, China) in accordance with the manufacturer’s instructions. Cells were collected for subsequent experiments 48 h after transfection.
Construction and Transfection of Lentiviral Igsf10-shRNA Expression Vector
For IGSF10 knockdown, we used a lentiviral RNAi small hairpin RNA (shRNA)-encoding system using the pHBLV-U6-Puro lentiviral RNAi vector (Hanbio, Shanghai, China). The targeting sequence of the Igsf10 shRNA was 5′-GCACCTTCCTGATTTCAAA-3′. To establish stable IGSF10-knockdown cell lines, MC3T3-E1 cells were transduced with the lentiviral RNAi vector at a multiplicity of infection (MOI) of approximately 100 in the presence of 10 μg/ml polybrene (Hanbio), and infected cells were selected with 2 μg/ml puromycin (Sigma-Aldrich). The empty lentivector was used as a negative control.
ALP Staining and ALP Activity Assay
MC3T3-E1 cells were cultured in OM 12-well plates for 7 days. Then, the cells were fixed with 4% paraformaldehyde for 15 min. After washing in double-distilled water (ddH2O) three times, the cells were stained using an ALP histochemical staining kit (Cwbiotech, Beijing, China) in accordance with the manufacturer’s protocol.
For ALP activity testing, MC3T3-E1 cells were collected after osteogenic induction for 7 days. The cells were lysed in RIPA buffer containing protease and phosphatase inhibitors (Huaxing bio, Beijing, China). Lysates were centrifuged at 13,523 × g at 4°C for 30 min and the supernatants were collected for analysis. Protein concentration was determined using a Pierce BCA protein assay kit (Thermo Fisher Scientific, Waltham, MA, United States). ALP activity was measured with an ALP activity kit (Nanjing Jiancheng Biotech, Nanjing, China) in accordance with the manufacturer’s instructions.
Alizarin Red Staining and Quantification
MC3T3-E1cells were fixed in 4% paraformaldehyde after osteogenic induction. Alizarin red staining was performed using 2% alizarin red (Sigma-Aldrich). For quantitative assessment of the degree of mineralization, the stain was dissolved in 10% cetylpyridinium chloride (Sigma-Aldrich) and the absorbance at 562 nm was measured.
Western Blotting
After MC3T3-E1 cells were cultured in OM for 7 days, total protein was extracted and quantified with the method described in the ALP activity assay. Samples with equal amounts of denatured protein were analyzed by sodium dodecyl sulphate polyacrylamide gel electrophoresis. Proteins were transferred to polyvinylidene difluoride membranes (Millipore, Bedford, MA, United States) and blocked in 5% skimmed milk. Then, membranes were incubated with primary antibodies against RUNX2 (CST, Danvers, MA, United States), BSP (CST), OSX (Abcam, Cambridge, United Kingdom), and GAPDH (Huaxing bio) at 4°C overnight. After washing in TBS-T, the membranes were incubated with HRP-conjugated secondary antibody (Huaxing bio) for 1 h at room temperature and visualized using an enhanced chemiluminescence blotting kit (Cwbiotech).
Statistical Analysis
All experiments were repeated at least three times and data are expressed as the mean ± SD. Differences between groups were determined by statistical analysis with GraphPad Prism v8.02 (GraphPad Software, Bethesda, MD, United States). Normal distribution and variance homogeneity tests were carried out. For the expression of IGSF10 in the osteogenesis of MC3T3-E1 cells, one-way analysis of variance (ANOVA) was performed and followed by Bonferroni’s post hoc test. Comparison between two groups was performed using unpaired Student’s t-test. p < 0.05 was considered statistically significant.
Results
Clinical Features of the Cleidocranial Dysplasia Patients
A Chinese family was studied. The proband, a 12-year-old girl, showed typical features of CCD, including patent fontanelles and bilateral hypoplastic clavicles on skull and chest radiography (Figures 1A–C). Furthermore, multiple dental anomalies commonly seen in CCD were also found, such as retained deciduous teeth, delayed eruption of permanent teeth, and crossbite (Figures 1D–G). A narrowing upper dental arch and high vault were also present (Figure 1F). Notably, supernumerary teeth were absent.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/748111/fgene-12-748111-HTML/image_m/fgene-12-748111-g001.jpg)
FIGURE 1. Clinical and radiographic findings of the proband and her affected mother. (A) Anteroposterior and (B) Lateral radiographs of the skull of the proband show patent frontal fontanelle and multiple Wormian bones. (C) Chest radiograph of the proband shows cone-shaped thorax and bilateral aplasia of clavicles. (D–G) Panoramic radiograph and Intraoral photographs of the proband show dental anomalies, including retained deciduous teeth and delayed eruption of permanent teeth. (H) Pedigree of the family. The arrow indicates the proband. (I) Anteroposterior and (J) Lateral radiographs of the skull of the proband’s mother show patent frontal fontanelle and multiple Wormian bones. (K) Chest radiograph of the proband’s mother shows cone-shaped thorax and bilateral aplasia of clavicles. (L) Panoramic radiograph of the proband’s mother shows dental anomalies, such as supernumerary teeth in the apices between the left upper canine and the first premolar. The typical signs of patients are pointed out by white arrows.
The pedigree of the family showed that the proband’s mother was also a CCD patient (Figure 1H). The phenotype of the proband’s mother was similar to the proband (Figure 1I–L), except that supernumerary teeth were present in the proband’s mother (Figure 1L). However, the number of supernumerary teeth was uncertain because the tooth extraction history of the mother was not reliable. No obvious signs of any other health problems were present in the two patients. After a detailed inquiry of the family medical history, we found that the proband’s two uncles and grandparent had no obvious CCD phenotypes. We were unable to confirm the genotype because they refused to carry out genetic testing.
No RUNX2 Mutations Were Detected
Sequencing analysis was first performed for the coding region, exon–intron boundaries, and the 2.0 kb region upstream of the transcription start site of RUNX2 by PCR and direct sequencing to explore the genetic basis of the familial case. However, no abnormality was found (data not shown). Because copy number variation of RUNX2 can cause CCD (Lee et al., 2008), we screened for intragenic deletions and duplications of RUNX2 by real-time PCR and MLPA. However, we detected no causative copy number variant of RUNX2 (Figures 2B,C). Chromosomal abnormalities involving RUNX2 can also result in CCD (Purandare et al., 2008), and therefore karyotype analysis was performed for the proband. A normal karyotype was shown for the proband (Figure 2A); therefore, we speculated that genetic variations in genes other than RUNX2 may lead to the disease.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/748111/fgene-12-748111-HTML/image_m/fgene-12-748111-g002.jpg)
FIGURE 2. Variation analysis of RUNX2 gene for the family CCD patients. (A) GTG-banding karyotype of lymphocytes of the proband reveals a normal karyotype of 46 chromosomes. (B, C) Analysis of RUNX2 copy number using RT-qPCR (B) and MLPA (C) shows a normal RUNX2 copy number.
Detection of an IGSF10 Mutation by WES
To identify the genetic cause of the familial case, we performed WES and bioinformatics analyses for the proband and her parents. WES identified 166,648 variants (Figure 3A) and a total of 67 candidate genes were identified after filtering for a MAF <0.5%, predicting deleterious variants, and segregating with the trait in affected family members (Supplementary Table S4). These 67 candidate genes were assessed by a literature review. Unexpectedly, no variants were detected in genes related to RUNX2 signaling. Therefore, we predicted that variation in an unknown gene may contribute to CCD in this family. We found that, of the 67 candidate genes, only IGSF10 was related to osteogenesis based on current knowledge. IGSF10 was finally selected for further investigation because of its potential role in bone development (Segev et al., 2004).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/748111/fgene-12-748111-HTML/image_m/fgene-12-748111-g003.jpg)
FIGURE 3. Mutation analysis of IGSF10 gene for the CCD patients. (A) Flowchart of whole-exome sequencing filtering outcomes. First, only variants that were found in all affected subjects in the family were retained. From among these segregating variants, only those with minor allele frequency (MAF) less than 0.005 in the 1,000 Genomes database and the NHLBI exome variant server were selected for further analysis. Deleterious variants were then picked out based on the SIFT, Polyphen-2 and variant impact. IGSF10 was scanned out by further literature review. (B) Identification of the heterozygous mutation in IGSF10 gene in the family of the proband and unaffected controls. The mutated sequence (upper panel) and control sequence (lower panel) were shown, respectively. Arrows indicated the mutation site. (C) Structure of the IGSF10 transcript and protein. Arrows indicated the mutation site. (D) Partial amino acid sequences alignment in IGSF10 protein among twelve species. The positions of the mutated amino acid in our study are indicated using a black box. (E) The three-dimensional structure of IGSF10 predicted by I-TASSER in wild-type and mutated IGSF10. WT shows the normal last 10 Ig domains, and MT displays the mutated form corresponding to WT. Yellow means the structure from the 2001th aa to the final aa. Red represents the 2001th aa.
A heterozygous two-base deletion in IGSF10 (c.6001_6002delCT, p.Leu2001Valfs*24) was detected in the proband and her mother (Supplementary Table S4). 50 unaffected controls were analyzed by PCR and Sanger sequencing to verify the mutation as pathogenic. Consistent with the WES data, a heterozygous two-base deletion of IGSF10 was detected only in the proband and her affected mother, but not in normal controls or her unaffected father (Figure 3B). This novel mutation was in exon 8 of IGSF10 (NM_178,,822.5) and resulted in a frameshift from codon 2001 to the resultant premature stop codon 2024. This was predicted to encode an IGSF10 protein truncated at the sixth immunoglobulin-like domain (Ig domain) with the loss of 600 amino acids from the C-terminus (Figure 3C).
Partial amino acid sequence alignment of the IGSF10 protein from 12 species showed that the affected amino acid, L2001, has a high level of evolutionary conservation among species (Figure 3D). Furthermore, the detected mutation was predicted by I-TASSER to alter the three-dimensional structure of the last ten Ig domains of IGSF10 (Figure 3E), leading to a loss of function of the predicted protein.
IGSF10 Induces Osteogenic Differentiation and Regulates Bone Sialoprotein
The role of IGSF10 in osteogenic differentiation was further explored. We evaluated the expression dynamics of Igsf10 mRNA after osteogenic induction in MC3T3-E1 cells and found that the levels of Igsf10 mRNA decreased upon OM treatment (Figure 4A). We then constructed an MC3T3-E1 cell line with stable IGSF10 knockdown. IGSF10 knockdown efficiency was first assessed by real-time PCR. The mRNA level of Igsf10 was decreased by more than 60% in MC3T3-E1 cells compared with control cells (Figure 4B). Regrettably, the level of IGSF10 protein could not be determined because a reliable IGSF10 antibody was not available.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/748111/fgene-12-748111-HTML/image_m/fgene-12-748111-g004.jpg)
FIGURE 4. The effect of IGSF10 on the osteogenic capacity of MC3T3-E1 cells. (A) The expression of Igsf10 in the osteogenesis of MC3T3-E1 cells was detected by real-time PCR. (B) The efficiency of Igsf10 knockdown in MC3T3-E1 cells was confirmed by real-time PCR. (C) ALP staining and Alizarin red staining after induction in osteogenic medium (OM). The bottom row indicates microscopic views. Scale bar, 500 μm. (D) ALP activity was analyzed after induction in OM. (E) Quantification of alizarin red staining by spectrophotometry. (F,G) The level of Runx2 expression after Igsf10 overexpression (F) or knockdown (G) transiently. (H,I) The expression of Igsf10 after Runx2 overexpression (H) or knockdown (I). (J,K) Western blotting and quantification of RUNX2, BSP, and OSX in MC3T3-E1 cells with stable Igsf10 knockdown after cultured in OM. (L–P) Quantitative analysis of the mRNA levels of Runx2, Bsp, Alp, Ocn, and Osx in MC3T3-E1 cells with stable Igsf10 knockdown after cultured in OM. NS: not significant. *p < 0.05.
We incubated stable IGSF10-knockdown cells and control cells (named si-Igsf10 and si-Control, respectively) in OM medium and analyzed their osteogenic ability. ALP activity (Figures 4C,D) and the ability to form mineralized nodules (Figures 4C,E) were obviously impaired in IGSF10-knockdown cells compared with control cells.
RUNX2 is a key protein in osteogenesis; therefore, we asked whether the role of IGSF10 in osteogenesis is mediated by RUNX2. Unexpectedly, the level of Runx2 expression did not change, irrespective of overexpression (Figure 4F) or knockdown of IGSF10 (Figure 4G). Similarly, Igsf10 expression was not affected by RUNX2 overexpression (Figure 4H) or knockdown (Figure 4I).
The effect of IGSF10 knockdown on other osteogenic-associated genes was subsequently explored. The level of Bsp mRNA decreased by 30% after IGSF10 knockdown (Figure 4M). The levels of RUNX2 and BSP proteins were similar to their mRNA levels after knockdown of IGSF10, with the level of BSP protein showing a slight reduction (Figures 4J,K). The mRNA levels of Alp, Ocn, and Osx showed no significant change after IGSF10 knockdown in MC3T3-E1 cells (Figures 4N–P).
Discussion
We studied the genetic basis of typical CCD in a family pedigree. We did not detect any sequence or copy number variation in RUNX2 by direct sequencing and copy-number analyses. However, we identified a new candidate gene for CCD, IGSF10, and explored its role in osteogenesis.
The diagnosis of CCD mainly depends on clinical observation and imaging examination. Detection of a pathogenic gene plays an auxiliary role when the clinical characteristics are not typical. The main clinical features of CCD include persistently open skull sutures with bulging calvaria, hypoplasia or aplasia of the clavicles permitting abnormal facility in apposing the shoulders, wide pubic symphysis, short middle phalanx of the fifth fingers, and multiple dental anomalies, such as retained deciduous teeth, delayed eruption of permanent teeth, and supernumerary teeth (Mundlos, 1999). As an autosomal dominant inherited disease with complete penetrance, the generally accepted causative gene of CCD is RUNX2. However, the expressivity of CCD is variable, even among different family members (Chitayat et al., 1992). Three different subtypes are defined for CCD, including classic, mild, and isolated dental anomalies (Zhou et al., 1999). Classic CCD patients show hypoplastic clavicles and delayed closure of the anterior fontanelle, in addition to characteristic craniofacial features. The isolated dental CCD phenotype only shows dental anomalies without cranial or axial skeletal anomalies. The phenotype of mild CCD is between classic CCD and dental CCD, with nominal or absent clavicle findings.
CCD should be differentiated from other similar skeletal diseases. Pycnodysostosis and mandibuloacral dysplasia are two disorders that should be considered as a differential diagnosis of cleidocranial dysplasia. Pycnodysostosis (MIM 265800) also shows hypoplastic clavicles, delayed closure of fontanelles, and delayed tooth eruption. Features that differentiate pycnodysostosis from CCD are acro-osteolysis, bone sclerosis with tendency to fracture, absence of supernumerary teeth, and cathepsin K (CTSK) gene mutations (Elmore et al., 1966; Bizaoui et al., 2019). Mandibuloacral dysplasia (MIM 248370) is an autosomal recessive disorder characterized by growth retardation, craniofacial anomalies with mandibular hypoplasia, skeletal abnormalities with progressive osteolysis of the distal phalanges and clavicles, and pigmentary skin changes (Cenni et al., 2018). Typical progressive osteolysis of the distal phalanges and pigmentary skin changes make it easy to differentiate from CCD (Garavelli et al., 2009). Furthermore, Yunis–Varon syndrome (MIM 216340), Crane–Heise syndrome (MIM 218090), CDAGS syndrome (MIM 603116), and hypophosphatasia (MIM 241500) also need to be differentiated from CCD (Crane and Heise, 1981; Basel-Vanagaite et al., 2008) (Méhes et al., 1972; Zand et al., 2003; Mendoza-Londono et al., 2005). In the present study, the proband showed typical features of CCD, including patent fontanelles, bilateral hypoplastic clavicles, and multiple dental anomalies. Therefore, this patient was confidently diagnosed with CCD.
Heterozygous mutations in RUNX2 have been identified as the cause of CCD (Mundlos et al., 1997). Heterozygous Runx2 mutant mice also display all the representative features of CCD, including open fontanelles and hypoplastic clavicles, but not the dental anomalies because mice only have one set of teeth with no replacement dentition (Selby et al., 1993; Komori et al., 1997; Otto et al., 1997). To date, more than 233 heterozygous mutations in RUNX2 have been identified in CCD patients, the majority of which are missense, nonsense, and frameshift mutations (Jaruga et al., 2016). Additionally, when sequence variants are not detected by direct sequencing, real-time PCR assays or MLPA are good tools to detect RUNX2 gene deletion or duplication by analyzing copy number variants (Lee et al., 2008; Franceschi et al., 2015; Sun et al., 2016). However, in the present proband, no abnormality was detected by direct sequencing in the coding region, flanking intron sequences, or regulatory region of RUNX2. Furthermore, no RUNX2 microdeletions/duplications were detected by real-time PCR or MLPA. Therefore, genetic heterogeneity or other genes in the RUNX2 pathway were alternative causative candidates for this familial case.
WES and bioinformatics analyses were used to explore the genetic basis of CCD in the present patient. We first focused on variants in genes involved in the RUNX2 pathway; however, our screen detected no pathogenic mutations. Subsequently, other potential candidate genes were analyzed and IGSF10 was finally identified as the most promising candidate. IGSF10 (also known as CMF608) belongs to the immunoglobulin superfamily that contain one or more Ig domains, and can bind with other proteins (Natarajan et al.). IGSF10 is a mechanical strain-induced bone-specific protein that is, involved in maintaining the osteochondroprogenitor cell pool (Segev et al., 2004). We detected a heterozygous two-base deletion in exon 8 of IGSF10 in the proband and her affected mother, resulting in a frameshift from codon 2001 to the resultant premature stop codon 2024. This leads to a truncated IGSF10 protein that has a severely affected three-dimensional structure. This is the first report of a mutation in IGSF10 being associated with CCD, although deleterious variants of IGSF10 have been described in DP and CHH patients (Howard et al., 2016; Amato et al., 2019). The primary manifestation of self-limited DP and CHH is delayed puberty, and no signs of delayed puberty or hypogonadism was shown in the proband following our subsequent follow-up. Interestingly, height and bone mineral density can be compromised in some adults with a history of DP (Zhu and Chan, 2017). Parker et al. also found that women with delayed puberty had a higher risk of osteoporosis (Parker et al., 2014). We therefore suggest that IGSF10 is involved in the pathogenesis of some bone-related diseases. However, the underlying mechanism is still unclear and needs further exploration.
IGSF10 plays an important role in bone development. Our research showed that Igsf10 mRNA expression was down-regulated in OM-treated MC3T3-E1 cells, and the osteogenic ability of MC3T3-E1 cells was decreased by IGSF10 knockdown. Moreover, previous research showed that IGSF10 expression was significantly elevated in vivo in response to bone-formation-promoting stimuli such as estrogen administration, and downregulated in response to sciatic neurotomy. The expression of IGSF10 was also confined to skeletal progenitor cells abundant in the regions of active bone modeling and remodeling (Segev et al., 2004). Taking all these results together, we presume that high IGSF10 expression may be associated with a specific developmental stage of MC3T3-E1 cells, and IGSF10 expression may be a precondition to osteogenic differentiation rather than a downstream target.
As a known causative gene of CCD, the expression of Runx2 after IGSF10 knockdown or overexpression was examined. Unexpectedly, IGSF10 over- or under-expression had no effect on the expression of Runx2, and vice versa. Previous studies showed different results. Overexpression of IGSF10 in ROS17/2.8 cells increased the expression of Runx2 (Segev et al., 2004), while overexpression of RUNX2 in NIH3T3 cells downregulated Igsf10 expression (Stephens and Morrison, 2014). There are two possible reasons for this phenomenon. First, the same treatment could result in different or contrasting results from different research backgrounds. In addition, other unknown regulatory mechanisms could be involved in IGSF10-mediated bone development.
In our further study, the expression of other osteoblast-related genes, including Alp, Osx, and Ocn, was not affected by IGSF10 knockdown, although the expression of Bsp was significantly decreased. BSP is a member of the “Small Integrin-Binding Ligand N-linked Glycoproteins” (SIBLING) extracellular matrix protein family of mineralized tissues (Bouleftour et al., 2016). It has long served as an early marker of osteoblast differentiation, as its expression is associated with the onset of mineralization (Chen et al., 1992; Midura et al., 2004). Taken together, these results indicate that the mechanisms by which IGSF10 regulates osteogenesis may be complex. One of the mechanisms may be that IGSF10 regulates early mineralization through pathways involving BSP.
In summary, the present study is a preliminary exploration of the genetic heterogeneity of CCD, and further studies need to be performed in the future to explore the underlying pathogenic mechanisms of IGSF10 mutation. More CCD patients without RUNX2 mutation need to be collected and their IGSF10 mutation status determined. Also, animal Igsf10-knockout models need to be constructed to explore the role of IGSF10 in bone development and osteogenesis. Recently, variations in deep intronic sequence have been reported as a cause of monogenic disorders as well as hereditary cancer syndromes (Vaz-Drago et al., 2017). Deep intronic mutations (i.e., more than 100 base pairs away from exon-intron boundaries) could lead to pseudo-exon inclusion, disrupt transcription regulatory motifs and inactivate non-coding RNA genes (Dhir and Buratti, 2010; Romano et al., 2013; Szafranski et al., 2013; Jafarifar et al., 2014). However, so far, there is no report in CCD caused by the deep intron mutation of RUNX2. Due to the limitations of our current study, we will next detect the RUNX2 deep intron in this family to further eliminate this possibility.
In conclusion, we report a familial case of CCD for whom no variation in RUNX2 was detected by direct sequencing or copy-number analyses. A new candidate gene for CCD, IGSF10, was identified by WES. We explored the function of IGSF10 in osteogenesis and speculate that IGSF10 may regulate early osteogenic differentiation by targeting BSP. Our results provide new genetic evidence that IGSF10 variation may contribute to CCD. While validating the link between IGSF10 and CCD, the mechanism by which IGSF10 mutation leads to CCD needs further exploration.
Data Availability Statement
The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/sra/PRJNA750609
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethical Committee of Peking University Health Science Center (approval number: IRB00001052-07100). Written informed consent to participate in this study was provided by the participants or the participants’ legal guardian/next of kin.
Author Contributions
DDL contributed to conception, design, data analysis, drafted and critically revised the manuscript; YL contributed to data acquisition and interpretation; XLZ contributed to data acquisition and interpretation; YXW contributed to conception, design, data analysis; CYZ contributed to conception, design, data analysis, drafting and critical revision of the manuscript; SGZ contributed to design, data analysis, and critically revised the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (Grant number 81771053), National Natural Science Foundation of China Youth Fund Project (Grant number 82001029) and Peking University School and Hospital of Stomatology Science Foundation for Young Scientists (Grant number PKUSS20160105).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are grateful to all the participants in this study and to the Department of Medical Genetics, Peking University Health Science Center, for their guidance in diagnosing genetic diseases. We thank Dr. Renny T. Franceschi (School of Dentistry, University of Michigan) for providing Runx2 cDNA plasmid and Quark Pharmaceuticals Inc. to donate the expression plasmid comprising the full length Igsf10 cDNA.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.748111/full#supplementary-material
References
Amato, L. G. L., Montenegro, L. R., Lerario, A. M., Jorge, A. A. L., Guerra Junior, G., Schnoll, C., et al. (2019). New Genetic Findings in a Large Cohort of Congenital Hypogonadotropic Hypogonadism. Eur. J. Endocrinol. 181 (2), 103–119. doi:10.1530/eje-18-0764
Basel-Vanagaite, L., Kornreich, L., Schiller, O., Yacobovich, J., and Merlob, P. (2008). Yunis-Varon Syndrome: Further Delineation of the Phenotype. Am. J. Med. Genet. 146A (4), 532–537. doi:10.1002/ajmg.a.32135
Baumert, U., Golan, I., Redlich, M., Aknin, J.-J., and Muessig, D. (2005). Cleidocranial Dysplasia: Molecular Genetic Analysis and Phenotypic-Based Description of a Middle European Patient Group. Am. J. Med. Genet. 139a (2), 78–85. doi:10.1002/ajmg.a.30927
Bizaoui, V., Michot, C., Baujat, G., Amouroux, C., Baron, S., Capri, Y., et al. (2019). Pycnodysostosis: Natural History and Management Guidelines From 27 French Cases and a Literature Review. Clin. Genet. 96 (4), 309–316. doi:10.1111/cge.13591
Bouleftour, W., Juignet, L., Bouet, G., Granito, R. N., Vanden-Bossche, A., Laroche, N., et al. (2016). The Role of the SIBLING, Bone Sialoprotein in Skeletal Biology - Contribution of Mouse Experimental Genetics. Matrix Biol. 52-54, 60–77. doi:10.1016/j.matbio.2015.12.011
Cenni, V., D’Apice, M. R., Garagnani, P., Columbaro, M., Novelli, G., Franceschi, C., et al. (2018). Mandibuloacral Dysplasia: A Premature Ageing Disease with Aspects of Physiological Ageing. Ageing Res. Rev. 42, 1–13. doi:10.1016/j.arr.2017.12.001
Chen, J., Shapiro, H. S., and Sodek, J. (1992). Developmental Expression of Bone Sialoprotein mRNA in Rat Mineralized Connective Tissues. J. Bone Miner Res. 7 (8), 987–997. doi:10.1002/jbmr.5650070816
Chitayat, D., Hodgkinson, K. A., and Azouz, E. M. (1992). Intrafamilial Variability in Cleidocranial Dysplasia: a Three Generation Family. Am. J. Med. Genet. 42 (3), 298–303. doi:10.1002/ajmg.1320420307
Crane, J. P., and Heise, R. L. (1981). New Syndrome in Three Affected Siblings. Pediatrics. 68 (2), 235–237.
Dhir, A., and Buratti, E. (2010). Alternative Splicing: Role of Pseudoexons in Human Disease and Potential Therapeutic Strategies. Febs j. 277 (4), 841–855. doi:10.1111/j.1742-4658.2009.07520.x
Elmore, S. M., Nance, W. E., McGee, B. J., Montmollin, M. E.-d., and Engel, E. (1966). Pycnodysostosis, With a Familial Chromosome Anomaly. Am. J. Med. 40 (2), 273–282. doi:10.1016/0002-9343(66)90108-2
Franceschi, R., Maines, E., Fedrizzi, M., Piemontese, M. R., De Bonis, P., Agarwal, N., et al. (2015). Familial Cleidocranial Dysplasia Misdiagnosed as Rickets Over Three Generations. Pediatr. Int. 57 (5), 1003–1006. doi:10.1111/ped.12692
Garavelli, L., D'Apice, M. R., Rivieri, F., Bertoli, M., Wischmeijer, A., Gelmini, C., et al. (2009). Mandibuloacral Dysplasia Type A in Childhood. Am. J. Med. Genet. 149a (10), 2258–2264. doi:10.1002/ajmg.a.33005
Garcia-Miñaur, S., Mavrogiannis, L. A., Rannan-Eliya, S. V., Hendry, M. A., Liston, W. A., Porteous, M. E. M., et al. (2003). Parietal Foramina With Cleidocranial Dysplasia Is Caused by Mutation in MSX2. Eur. J. Hum. Genet. 11 (11), 892–895. doi:10.1038/sj.ejhg.5201062
Goto, T., Aramaki, M., Yoshihashi, H., Nishimura, G., Hasegawa, Y., Takahashi, T., et al. (2004). Large Fontanelles Are a Shared Feature of Haploinsufficiency of RUNX2 and its Co-Activator CBFB. Congenit. Anom. 44 (4), 225–229. doi:10.1111/j.1741-4520.2004.00043.x
Hintzsche, J. D., Robinson, W. A., and Tan, A. C. (2016). A Survey of Computational Tools to Analyze and Interpret Whole-Exome Sequencing Data. Int. J. Genomics. 2016, 1–16. doi:10.1155/2016/7983236
Howard, S. R., Guasti, L., Ruiz‐Babot, G., Mancini, A., David, A., Storr, H. L., et al. (2016). IGSF 10 Mutations Dysregulate Gonadotropin‐Releasing Hormone Neuronal Migration Resulting in Delayed Puberty. EMBO Mol. Med. 8 (6), 626–642. doi:10.15252/emmm.201606250
Jafarifar, F., Dietrich, R. C., Hiznay, J. M., and Padgett, R. A. (2014). Biochemical Defects in Minor Spliceosome Function in the Developmental Disorder MOPD I. Rna. 20 (7), 1078–1089. doi:10.1261/rna.045187.114
Jaruga, A., Hordyjewska, E., Kandzierski, G., and Tylzanowski, P. (2016). Cleidocranial Dysplasia and RUNX2-Clinical Phenotype-Genotype Correlation. Clin. Genet. 90 (5), 393–402. doi:10.1111/cge.12812
Khan, A., Hyde, R. K., Dutra, A., Mohide, P., and Liu, P. (2006). Core Binding Factor Beta (CBFB) Haploinsufficiency Due to an Interstitial Deletion at 16q21q22 Resulting in Delayed Cranial Ossification, Cleft Palate, Congenital Heart Anomalies, and Feeding Difficulties but Favorable Outcome. Am. J. Med. Genet. 140A (21), 2349–2354. doi:10.1002/ajmg.a.31479
Komori, T. (2020). Molecular Mechanism of Runx2-Dependent Bone Development. Mol. Cell. 43 (2), 168–175. doi:10.14348/molcells.2019.0244
Komori, T., Yagi, H., Nomura, S., Yamaguchi, A., Sasaki, K., Deguchi, K., et al. (1997). Targeted Disruption of Cbfa1 Results in a Complete Lack of Bone Formation Owing to Maturational Arrest of Osteoblasts. Cell. 89 (5), 755–764. doi:10.1016/s0092-8674(00)80258-5
Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: Molecular Evolutionary Genetics Analysis Across Computing Platforms. Mol. Biol. Evol. 35 (6), 1547–1549. doi:10.1093/molbev/msy096
Lee, M. T. M., Tsai, A. C.-H., Chou, C.-H., Sun, F.-M., Huang, L.-C., Yen, P., et al. (2008). Intragenic Microdeletion of RUNX2 Is a Novel Mechanism for Cleidocranial Dysplasia. Hugo J. 2 (1-2), 45–49. doi:10.1007/s11568-008-9024-y
Li, H., and Durbin, R. (2009). Fast and Accurate Short Read Alignment With Burrows-Wheeler Transform. Bioinformatics. 25 (14), 1754–1760. doi:10.1093/bioinformatics/btp324
Méhes, K., Klujber, L., Lassu, G., and Kajtár, P. (1972). Hypophosphatasia: Screening and Family Investigations in an Endogamous Hungarian Village. Clin. Genet. 3 (1), 60–66. doi:10.1111/j.1399-0004.1972.tb01727.x
Mendoza-Londono, R., Lammer, E., Watson, R., Harper, J., Hatamochi, A., Hatamochi-Hayashi, S., et al. (2005). Characterization of a New Syndrome that Associates Craniosynostosis, Delayed Fontanel Closure, Parietal Foramina, Imperforate Anus, and Skin Eruption: CDAGS. Am. J. Hum. Genet. 77 (1), 161–168. doi:10.1086/431654
Midura, R. J., Wang, A., Lovitch, D., Law, D., Powell, K., and Gorski, J. P. (2004). Bone Acidic Glycoprotein-75 Delineates the Extracellular Sites of Future Bone Sialoprotein Accumulation and Apatite Nucleation in Osteoblastic Cultures. J. Biol. Chem. 279 (24), 25464–25473. doi:10.1074/jbc.M312409200
Mundlos, S. (1999). Cleidocranial Dysplasia: Clinical and Molecular Genetics. J. Med. Genet. 36 (3), 177–182.
Mundlos, S., Otto, F., Mundlos, C., Mulliken, J. B., Aylsworth, A. S., Albright, S., et al. (1997). Mutations Involving the Transcription Factor CBFA1 Cause Cleidocranial Dysplasia. Cell. 89 (5), 773–779. doi:10.1016/s0092-8674(00)80260-3
Napierala, D., Garcia-Rojas, X., Sam, K., Wakui, K., Chen, C., Mendoza-Londono, R., et al. (2005). Mutations and Promoter SNPs in RUNX2, a Transcriptional Regulator of Bone Formation. Mol. Genet. Metab. 86 (1-2), 257–268. doi:10.1016/j.ymgme.2005.07.012
Natarajan, K., Mage, M. G., and Margulies, D. H. (2015). Immunoglobulin Superfamily. eLS. Editors John Wiley & Sons, Ltd, 1–7. doi:10.1002/9780470015902.a0000926.pub2
Northup, J. K., Matalon, R., Lockhart, L. H., Hawkins, J. C., and Velagaleti, G. V. N. (2011). 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. 54 (4), e394–e398. doi:10.1016/j.ejmg.2011.03.011
Ott, C. E., Leschik, G., Trotier, F., Brueton, L., Brunner, H. G., Brussel, W., et al. (2010). Deletions of the RUNX2 Gene Are Present in About 10% of Individuals With Cleidocranial Dysplasia. Hum. Mutat. 31 (8), E1587–E1593. doi:10.1002/humu.21298
Ott, C. E., Hein, H., Lohan, S., Hoogeboom, J., Foulds, N., Grünhagen, J., et al. (2012). Microduplications Upstream ofMSX2are Associated With a Phenocopy of Cleidocranial Dysplasia. J. Med. Genet. 49 (7), 437–441. doi:10.1136/jmedgenet-2012-100825
Otto, F., Kanegane, H., and Mundlos, S. (2002). Mutations in theRUNX2 Gene in Patients With Cleidocranial Dysplasia. Hum. Mutat. 19 (3), 209–216. doi:10.1002/humu.10043
Otto, F., Thornell, A. P., Crompton, T., Denzel, A., Gilmour, K. C., Rosewell, I. R., et al. (1997). Cbfa1, a Candidate Gene for Cleidocranial Dysplasia Syndrome, Is Essential for Osteoblast Differentiation and Bone Development. Cell. 89 (5), 765–771. doi:10.1016/s0092-8674(00)80259-7
Parker, S. E., Troisi, R., Wise, L. A., Palmer, J. R., Titus-Ernstoff, L., Strohsnitter, W. C., et al. (2014). Menarche, Menopause, Years of Menstruation, and the Incidence of Osteoporosis: the Influence of Prenatal Exposure to Diethylstilbestrol. J. Clin. Endocrinol. Metab. 99 (2), 594–601. doi:10.1210/jc.2013-2954
Purandare, S. M., Mendoza-Londono, R., Yatsenko, S. A., Napierala, D., Scott, D. A., Sibai, T., et al. (2008). 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. 146a (4), 453–458. doi:10.1002/ajmg.a.31750
Quack, I., Vonderstrass, B., Stock, M., Aylsworth, A., Becker, A., Brueton, L., et al. (1999). Mutation Analysis of Core Binding Factor A1 in Patients With Cleidocranial Dysplasia. Am. J. Hum. Genet. 65 (5), 1268–1278. doi:10.1086/302622
Raivio, T., and Miettinen, P. J. (2019). Constitutional Delay of Puberty Versus Congenital Hypogonadotropic Hypogonadism: Genetics, Management and Updates. Best Pract. Res. Clin. Endocrinol. Metab. 33 (3), 101316. doi:10.1016/j.beem.2019.101316
Romano, M., Buratti, E., and Baralle, D. (2013). Role of Pseudoexons and Pseudointrons in Human Cancer. Int. J. Cel Biol. 2013, 1–16. doi:10.1155/2013/810572
Segev, O., Samach, A., Faerman, A., Kalinski, H., Beiman, M., Gelfand, A., et al. (2004). CMF608-a Novel Mechanical Strain-Induced Bone-Specific Protein Expressed in Early Osteochondroprogenitor Cells. Bone. 34 (2), 246–260. doi:10.1016/j.bone.2003.10.003
Selby, P. B., Bolch, S. N., Mierzejewski, V. S., McKinley, T. W., and Raymer, G. D. (1993). Synergistic Interactions between Two Skeletal Mutations in Mice: Individual and Combined Effects of the Semidominants Cleidocranial Dysplasia (Ccd) and Short Digits (Dsh). J. Hered. 84 (6), 466–474. doi:10.1093/oxfordjournals.jhered.a111373
Stephens, A. S., and Morrison, N. A. (2014). Novel Target Genes of RUNX2 Transcription Factor and 1,25-Dihydroxyvitamin D3. J. Cel. Biochem. 115 (9), 1594–1608. doi:10.1002/jcb.24823
Sun, Y. C., Zhang, T. T., Li, X. M., Zhou, H., Zhang, Z. L., Zhou, J., et al. (2016). A Novel Intragenic Microdeletion in RUNX2 in a Chinese Family With Cleidocranial Dysplasia. Int. J. Clin. Exp. Pathol. 9 (7), 7209–+.
Szafranski, P., Yang, Y., Nelson, M. U., Bizzarro, M. J., Morotti, R. A., Langston, C., et al. (2013). Novel FOXF1 Deep Intronic Deletion Causes Lethal Lung Developmental Disorder, Alveolar Capillary Dysplasia With Misalignment of Pulmonary Veins. Hum. Mutat. 34 (11), 1467–1471. doi:10.1002/humu.22395
Vaz-Drago, R., Custódio, N., and Carmo-Fonseca, M. (2017). Deep Intronic Mutations and Human Disease. Hum. Genet. 136 (9), 1093–1111. doi:10.1007/s00439-017-1809-4
Yang, J., Yan, R., Roy, A., Xu, D., Poisson, J., and Zhang, Y. (2015). The I-TASSER Suite: Protein Structure and Function Prediction. Nat. Methods. 12 (1), 7–8. doi:10.1038/nmeth.3213
Zand, D. J., Carpentieri, D., Huff, D., Medne, L., Napierala, D., Lee, B., et al. (2003). Crane-Heise Syndrome: a Second Familial Case Report With Elaboration of Phenotype. Am. J. Med. Genet. 118a (3), 223–228. doi:10.1002/ajmg.a.10026
Zhang, T., Wu, J., Zhao, X., Hou, F., Ma, T., Wang, H., et al. (2019). Whole-Exome Sequencing Identification of a Novel Splicing Mutation of RUNX2 in a Chinese Family With Cleidocranial Dysplasia. Arch. Oral Biol. 100, 49–56. doi:10.1016/j.archoralbio.2019.02.005
Zhang, X., Liu, Y., Wang, X., Sun, X., Zhang, C., and Zheng, S. (2017). Analysis of Novel RUNX2 Mutations in Chinese Patients With Cleidocranial Dysplasia. PLoS One. 12 (7), e0181653. doi:10.1371/journal.pone.0181653
Zhou, G., Chen, Y., Zhou, L., Thirunavukkarasu, K., Hecht, J., Chitayat, D., et al. (1999). CBFA1 Mutation Analysis and Functional Correlation With Phenotypic Variability in Cleidocranial Dysplasia. Hum. Mol. Genet. 8 (12), 2311–2316. doi:10.1093/hmg/8.12.2311
Keywords: cleidocranial dysplasia, RUNX2, IGSF10, osteoblast differentiation, mutation analysis
Citation: Liu D, Liu Y, Zhang X, Wang Y, Zhang C and Zheng S (2021) An Exploration of Mutagenesis in a Family with Cleidocranial Dysplasia without RUNX2 Mutation. Front. Genet. 12:748111. doi: 10.3389/fgene.2021.748111
Received: 27 July 2021; Accepted: 29 September 2021;
Published: 19 October 2021.
Edited by:
Ramcés Falfán-Valencia, Instituto Nacional de Enfermedades Respiratorias-México (INER), MexicoReviewed by:
Roopashri R. K., A. J. Institute of Dental Sciences, IndiaNancy Monroy-Jaramillo, National Institute of Neurology and Neurosurgery, Mexico
Copyright © 2021 Liu, Liu, Zhang, Wang, Zhang and Zheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chenying Zhang, zhangchy168@163.com; Shuguo Zheng, kqzsg86@bjmu.edu.cn