- Research article
- Open Access
A novel COMP mutation in a Chinese family with multiple epiphyseal dysplasia
BMC Medical Genetics volume 21, Article number: 115 (2020)
Multiple epiphyseal dysplasia (MED) is a skeletal disorder characterized by delayed and irregular ossification of the epiphyses and early-onset osteoarthritis. At least 66% of the reported autosomal dominant MED (AD-MED) cases are caused by COMP mutations.
We recruited a four-generation Chinese family with early-onset hip osteoarthritis, flatfoot, brachydactyly, and mild short stature. An assessment of the family history, detailed physical examinations, and radiographic evaluations were performed on the proband and other family members, followed by the performance of whole-exome sequencing (WES). The pathogenicity of the candidate mutation was also analyzed.
An AD-MED family with 10 affected members and 17 unaffected members was recruited. The main radiographic findings were symmetrical changes in the dysplastic acetabulum and femoral heads, irregular contours of the epiphyses, a shortened femoral neck, and flatfoot. Lower bone density was also observed in the ankle joints, wrist joints, and knees, as well as irregular vertebral end plates. In the proband, we identified the missense mutation c.1153G > T (p. Asp385Tyr), located in exon 11 of the COMP gene. This mutation was assessed as 'pathogenic' because of its low allele frequency and its high likelihood of co-segregation with disease in the reported family. Sanger sequencing validated the novel heterozygous mutation c.1153G > T (p. Asp385Tyr) in exon 11 of COMP in all affected individuals in the family.
Our results underlined a key role of the Asp385 amino acid in the protein function of COMP and confirmed the pathogenicity of the COMP (c.1153G > T; p. Asp385Tyr) mutation in AD-MED disease. We have therefore expanded the known mutational spectrum of COMP and revealed new phenotypic information for AD-MED.
Multiple epiphyseal dysplasia (MED; MIM# 132400) is a skeletal disorder characterized by delayed and irregular ossification of the epiphyses as well as early-onset osteoarthritis . To date, six genes have been associated with MED, including five genes that cause autosomal dominant MED (AD-MED; COMP, COL9A1, COL9A2, COL9A3 and MATN3) and one gene that causes autosomal recessive MED (rMED; SLC26A2) [2,3,4]. The incidence of AD-MED is estimated to be 1 in 10,000 individuals, and at least 66% of reported AD-MED cases are caused by COMP mutations . These cases are known as EDM1 (or COMP-MED), and are characterized by mild short stature, premature osteoarthritis of load-bearing joints, and abnormalities of the epiphyses of hands, long bones, and hips [6, 7].
The COMP gene encodes cartilage oligomeric matrix protein (COMP) made up of 757 amino acids . COMP is a 552 kDa pentameric adhesive glycoprotein that is mainly found in synovium, tendons, ligaments, and the extracellular matrix of cartilage [9, 10]. The binding of COMP to extracellular matrix proteins is essential for the integrity of the cartilage and extracellular matrix . Since the 1990s, more than 80 novel mutations involved in the pathogenesis of MED have been reported in COMP [9, 11]. The locations of these mutations are predominantly concentrated in the highly conserved type III (T3) calcium-binding repeat domain, and it has been demonstrated that the T34 mutation is significantly associated with MED compared with other T3 repeats. These mutations affect the secretion of extracellular matrix proteins and extracellular matrix integrity, often leading to skeletal abnormalities including pseudoachondroplasia (PSACH) and MED .
Here, we report a four-generation Chinese family with early-onset hip osteoarthritis, flatfoot, brachydactyly, and mild short stature. An assessment of the family history, detailed physical examinations, and radiographic evaluations were performed on the proband and other family members, followed by the performance of whole-exome sequencing (WES). The pathogenicity of the candidate mutation was then analyzed.
The proband (III-10) was initially diagnosed with Legg–Calvé–Perthes disease based on the observed radiographic changes, including uneven density of the bilateral femoral head and bilateral femoral head collapse, which indicates avascular necrosis of the femoral head. After an investigation of the proband’s family history, a four-generation pedigree with 27 family members was recruited from Shanxi Province, China. Clinical evaluation including medical history, physical examination, and radiographic assessment, as well as peripheral venous blood (EDTA-K2 anticoagulant) was collected from all recruited individuals. For patients who cannot visit our hospital for clinical evaluation, we collected their clinical information by telephone interviews. Their consent forms and anteroposterior radiographs of pelvis were delivered by the express. Written informed consent was obtained from each participant; if the participant was younger than 16 years old, written informed consent was obtained from their parents or legal guardians. The Ethics Committee of Peking Union Medical College Hospital (PUMCH) approved this study.
DNA preparation and WES
According to the manufacturer’s protocols, genomic DNA samples were extracted from peripheral blood leukocytes of each family member by peripheral blood DNA extraction kit (QIAamp DNA Blood Mini Kit; Qiagen, Germany). Purified DNA was qualified by Nanodrop2000 (Thermo Fisher Scientific, Waltham, MA, USA) and quantified by Qubit 3.0 using the dsDNA HS Assay Kit (Life Technologies, Carlsbad, CA, USA).
WES was performed on genome DNA of the proband (III-10) and individual II-6, II-9, III-6, III-18, and IV-12. The genomic DNA was broken into 180 to 280 base pair (bp) fragments by ultrasonoscope. Illumina paired-end libraries were prepared from DNA samples. The exome sequenced captured by the SureSelect Human All Exon V6 + UTR r2 core design (91 Mb, Agilent) were sequenced on an Illumina HiSeq 4000 platform (Illumina, San Diego, CA, USA). The raw sequencing data were analyzed through Genome Analysis Toolkit (GATK, Version 3.4.0). The error assessment, variant calling and annotation were performed through in-house developed Peking Union Medical College Hospital Pipeline (PUMP) .
The interpretation of variants was performed according to the American College of Medical genetics and Genomics (ACMG) guideline . The variants were filtered through the following procedures: (1) the allele frequency of variants is required to be less than 1% or absent from databases like 1000 Genomes Project (The 1000 Genomes Project Consortium 2015) and Genome Aggregation Database (gnomAD) (http://gnomad-old.broadinstitute.org/). (2) Variants were filtered out when they were synonymous mutation or located in introns without influence on splicing and biological function. (3) Variants were reviewed by taken into consideration of phenotype evaluation, inheritance model, reported documents as well databases like Human Gene Mutation Database (HGMD) and Online Mendelian Inheritance in Man (OMIM) (available at: https://omim.org/).
A pair of PCR (Polymerase chain reaction) primers (F: 5′-CCATGAAGTTGGGACTCTGT-3′, R: 5′-GGTCATTTCTCTGGCAGTGT-3′) were designed using Primer 3 (http://primer3.ut.ee/) to amplify the Exon 11 of COMP gene. The PCR program was 95 °C for 3 min, followed by 38 cycles at 94 °C for 30 s, 58–60 °C for 30 s, 72 °C for 50 s, and a final extension at 72 °C for 8 min. All the collected samples in this family were used as DNA temple, the products of PCR analyzed in 2% agarose gels and purified by QIA quick PCR purification kit (Qiagen, Germantown, USA). Sequencing was performed on ABI3700 sequence Detection System (Applied Biosystems, Inc., Foster City, CA, USA). The reference sequence of the candidate gene was obtained from UCSC Genome Browser (http://genome.ucsc.edu) and compared with sequencing data through CodonCode Aligner (version 220.127.116.11; CondonCode, Centerville, MA, USA).
3-D protein structures of both wild and variant type COMP protein were predicted using an online modeling server, SWISS-MODEL program (https://swissmodel.expasy.org/), which were then viewed and edited by the molecular visualization system PyMOL (PyMOL Molecular Graphics System, Version 2.3.3, Schrödinger, LLC).
We recruited an AD-MED family with 10 affected members and 17 unaffected members (Fig. 1a). The proband (III-10) was a 38-year-old woman who had developed a waddling gait at the age of 6 years (Fig. 2). At around the age of 22 years, the patient reported pain in the bilateral hip, lumbar vertebrae, and bilateral knees. The height of the proband was normal (152 cm), and there was no evidence suggesting growth retardation in her childhood. The main radiographic findings were symmetrical changes in the dysplastic acetabulum and femoral heads, irregular contours of the epiphyses, a shortened femoral neck, and flatfoot (Fig. 2). Decreased bone density of the ankle joints, wrist joints, and knees was observed, as well as irregular vertebral end plates. Her fingers and elbows were normal. When examining the other affected family members, we found one patient (III-2) with a mild phenotype (only presenting with the symptoms of hip pain), one patient (III-18) with mild sacroiliitis, and two patients with brachydactyly (II-9, II-13) (Fig. 3). All affected individuals in this family had flatfoot, except for one patient (III-2) (Table 1). Furthermore, compared with the general Chinese male population (average height = 170 cm), male patients in this family have generally shorter stature (≤165 cm).
After the raw sequencing data were processed through the analytical pipeline [13, 15], we identified that the individuals who underwent WES in total carried 3185 rare single nucleotide variants (SNVs), including missense, frameshift, splicing, and nonsense SNVs, as well as those with unknown influence, such as the synonymous or non-coding variants. All pathogenic and likely pathogenic variants were manually reviewed according to ACMG guidelines and the OMIM database. As a result, we identified the missense mutation, c.1153G > T (p. Asp385Tyr), located in exon 11 of COMP. This mutation was assessed as 'pathogenic' based on its low allele frequency and its high likelihood of co-segregation with disease in the reported family. Sanger sequencing was performed in 24 family members, including 7 affected members and 17 unaffected members (Fig. 1a). All affected family members carried the heterozygous mutation c.1153G > T (p. Asp385Tyr) in exon 11 of COMP, while the unaffected family members did not (Fig. 1b). These results further indicated that this mutation was co-segregated in our family and provides strong evidence for the pathogenicity of this mutation. The three-dimensional structure of the COMP protein provided further evidence of pathogenicity, with the mutation resulting in the replacement of the long side chain of Asp385 by a phenolic hydroxyl of tyrosine (Fig. 4).
In a previous study by Mabuchi et al., a different mutation at the same position (c.1153G > A, p. Asp385Asn) has been reported to cause MED . Furthermore, Jackson et al.  also identified this recurrent mutation (c.1153G > A, p. Asp385Asn) in two British families and one Dutch family with MED. In addition, Liu et al. recently described an AD-MED family in which multiple members had been diagnosed with ANFH. Using WES, they identified a heterozygous variant in COMP (c.1153G > A) that contributed to the disease phenotype .
Consistent with the previously reported phenotypes, our pedigree cases also showed symptoms of classical AD-MED, such as femoral head necrosis, mild short stature, early-onset osteoarthritis of the knee and hip, and brachydactyly. These symptoms highlight the key role for Asp385 amino acid in the protein function of COMP and confirm the pathogenicity of COMP (c.1153G > T; p. Asp385Tyr) in AD-MED disease.
In the current study, we investigated a four-generation family with early-onset hip osteoarthritis caused by a heterozygous c.1153G > T mutation. Their prominent symptoms included severe osteoarthritis, knee pain, and femoral head necrosis, which is consistent with COMP-associated MED. According to the ACMG guideline and our study, this variation: a) is absent from the records of any frequency database, such as the 1000 Genomes Project and gnomAD, indicating that this mutation is rare (ACMG pathogenicity criteria: PM2); b) is located in a well-established functional domain (PM1); c) is co-segregated with disease in multiple affected family members (PP1); and d) has computational evidence supporting a deleterious effect on the gene (PP3). Based on these lines of evidence, we classified this mutation as pathogenic, with an important role in the genetic etiology of our family.
A different mutation at the same position (c.1153G > A, p. D385N) has been reported to cause MED . In this previous study, the reported sporadic case had mild short stature and early onset osteoarthrosis, and was diagnosed with a severe form of MED—the “Fairbank type”. In another study, the authors reported a pedigree with severe hip osteoarthrosis. Using WES, a c.1153G > A mutation of the COMP gene was identified to be related to the MED phenotype of this family . The phenotypes of the pedigrees reported by Liu et al.  were partially similar to that of our patients. They mainly manifested as necrosis of the femoral head, but their phenotype did not include flatfoot. In the current study, the c.1153G > T mutation was predicted to lead to an amino acid substitution from Asp to Tyr, implying that other substitutions at this position can also lead to the classical MED phenotype.
The precise function and pathogenic mechanisms of mutant COMP in MED have not been fully defined. However, compelling evidence indicates that COMP proteins play an important role in maintaining cartilage and extracellular matrix integrity. The misfolding of mutant COMP affects its normal secretion from the endoplasmic reticulum of chondrocytes, and this intracellular retention is toxic to these cells, resulting in premature chondrocyte death [8, 18]. These events reduce the number of chondrocytes in the growth plate, ultimately reducing linear growth, and the phenotypic outcome is dwarfism. Moreover, the reduction of COMP secretion also affects the assembly of collagen fibers , leading to a decrease in articular cartilage mechanical strength and the occurrence of early-onset osteoarthritis [12, 20]. In addition, the calcium-binding T3 repeat of COMP has been found to provide support for chondrocyte attachment . Changes in the three-dimensional calcium-dependent structure of the mutant COMP may therefore alter chondrocyte attachment, thereby contributing to MED phenotype development. Briggs et al.  previously confirmed that both PSACH and MED-associated mutations are predominantly located within the T3 repeat domain of COMP (90% of mutations). They also reported that missense mutations and in-frame insertions/deletions of single residues in T35–7 usually cause PSACH, while missense mutations in T33–4 are more likely to cause MED. Our novel c.1153G > A mutation is located in the T33–4 repeat of COMP. The aforementioned mechanisms may explain some of the phenotypes in our family, such as mild short stature and early-onset osteoarthritis.
Clubfoot is a rare radiological finding that has been observed mostly in association with rMED. In a previous study , rMED patients with clubfoot were reported to carry homozygous/compound heterozygous mutations in SLC26A2 at birth. In addition, Superti-Furga et al.  described rMED patients with normal stature, clubfoot, and double-layered patellae caused by a DTDST mutation. In our study, all affected individuals in the family had flatfoot except for a single patient (III-2). However, other intra-familial differences in phenotype were observed in our family. For example, one patient (III-18) had mild sacroiliitis, with lower back pain at the age of 18, and also had gait abnormalities, mild short stature, flatfoot, and hip osteoarthritis. In addition, two male patients (II-9, II-13) had brachydactyly, but the other patients in this family did not have this phenotype. Furthermore, the male patients in this family were generally shorter in height, but it is unclear if this is related to the c.1153G > T mutation. In a pedigree study conducted by Sakamoto et al. , intra-familial differences in severity were also observed in their four-generation family: radiological manifestations in the knees were more severe in the proband’s father than in the proband, and the proband’s young sister was of shorter stature than the proband. Liu et al.  also reported intra-familial differences; for example, a twin brother in their reported family had more severe walking limitations than other family members. These intra-familial phenotypic differences may be difficult to explain by genetic factors but may be explained by the effects of environmental factors .
We identified a novel heterozygous pathogenic mutation in COMP from an AD-MED family that exhibited COMP-associated MED, and other phenotypes including flatfoot. Our results expanded both the mutational and phenotypic spectra of COMP and suggested that this mutation of a key amino acid residue is disease-causing.
Availability of data and materials
The datasets generated during the current study are available in the Mendeley repository, https://doi.org/10.17632/gcfxv4yws3.1 (DOI: https://doi.org/10.17632/gcfxv4yws3.1). The datasets analysed during the current study include the variant frequency data from the gnomAD Browser repository; https://storage.googleapis.com/gnomad-public/release/2.1.1/vcf/exomes/gnomad.exomes.r2.1.1.sites.vcf.bgz; the variant frequency data from the 1000 Genomes Project repository: ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/; the gene-disease relationship data from OMIM repository: https://www.omim.org/static/omim/data/mim2gene.txt; the hg19 reference genome from the UCSC Genome Browser repository: http://hgdownload.soe.ucsc.edu/goldenPath/hg19/bigZips/; the Genome Analysis Toolkit (GATK, Version 3.4.0) software: https://github.com/broadinstitute/gatk/releases/tag/18.104.22.168. The variant-disease relationship data from the Human Gene Mutation Database (HGMD) is not publicly available (the public version of this database is available only to registered users from academic institutions / non-profit organisations or commercial users who purchase a license of HGMD Professional), but can be required by contacting the corresponding author (Nan Wu) at email@example.com.
Multiple epiphyseal dysplasia
Autosomal dominant MED
Cartilage oligomeric matrix protein
Peking Union Medical College Hospital
American College of Medical genetics and Genomics
Online Mendelian Inheritance in Man
Human Gene Mutation Database
Polymerase chain reaction
Single nucleotide variants
Autosomal recessive MED
Steven A, Richard M, William S, Michael M. Multiple epiphyseal dysplasia. J Am Acad Orthop Surg. 2015;23:164–72.
Dasa V, Eastwood J, Podgorski M, Park H, Blackstock C, Antoshchenko T, et al. Exome sequencing reveals a novel COL2A1 mutation implicated in multiple epiphyseal dysplasia. Am J Med Genet A. 2019;179:534–41.
Kausar M, Mäkitie R, Toiviainen-Salo S, Ignatius J, Anees M, Mäkitie O. Recessive multiple epiphyseal dysplasia - clinical characteristics caused by rare compound heterozygous SLC26A2 genotypes. Eur J Med Genet. 2019;62:103573.
Seo S, Song HR, Kim H, Yoo W, Shim J, Chung C, et al. Comparison of orthopaedic manifestations of multiple epiphyseal dysplasia caused by MATN3 versus COMP mutations: a case control study. BMC Musculoskelet Dis. 2014;15:84.
Anthony S, Munk R, Skakun W, Masini M. Multiple epiphyseal dysplasia. J Am Acad Orthop Surg. 2015;23:164–72.
Liu FX, Li YX, Zhang XD, Ren CA, Huang SZ, Yu MX. EDM1: a novel point mutation in cartilage oligomeric matrix protein gene in a Chinese family with multiple epiphyseal dysplasia. Chin Med J (Engl). 2013;126:1103–7.
Liu HY, Xiao JF, Huang J, Wang Y, Wu D, Li T, et al. Diagnosis with multiple epiphyseal dysplasia using whole-exome sequencing in a chinese family. Chin Med J (Engl). 2017;130:104–7.
Posey KL, Coustry F, Hecht JT. Cartilage oligomeric matrix protein: COMPopathies and beyond. Matrix Biol. 2018;71–2:161–73.
Briggs MD, Brock J, Ramsden SC, Bell PA. Genotype to phenotype correlations in cartilage oligomeric matrix protein associated chondrodysplasias. Eur J Hum Genet. 2014;22:1278–82.
Yang X, Ruan G, Xu J, Zheng S, Wang K, Ding C. Associations between suprapatellar pouch effusion-synovitis, serum cartilage oligomeric matrix protein, high sensitivity C-reaction protein, knee symptom, and joint structural changes in patients with knee osteoarthritis. Clin Rheumatol. 2020. https://doi.org/10.1007/s10067100190490510067.
Cohn DH, Briggs MD, King LM, Rimoin DL, Wilcox WR, Lachman RS, et al. Mutations in the cartilage oligomeric matrix protein (COMP) gene in pseudoachondroplasia and multiple epiphyseal dysplasia. Ann N Y Acad Sci. 2010;785:188–94.
Acharya C, Yik JH, Kishore A, Van Dinh V, Di Cesare PE, Haudenschild DR. Cartilage oligomeric matrix protein and its binding partners in the cartilage extracellular matrix: interaction, regulation and role in chondrogenesis. Matrix Biol. 2014;37:102–11.
Wang K, Zhao S, Zhang Q, Yuan J, Liu J, Ding X, et al. Whole-exome sequencing reveals known and novel variants in a cohort of intracranial vertebral–basilar artery dissection (IVAD). J Hum Genet. 2018;63:1119–28.
Liu J, Zhou Y, Qi X, Chen J, Chen W, Qiu G, et al. CRISPR/Cas9 in zebrafish: an efficient combination for human genetic diseases modeling. Hum Genet. 2017;136:1–12.
Wang K, Zhao S, Liu B, Zhang Q, Li Y, Liu J, et al. Perturbations of BMP/TGF-β and VEGF/VEGFR signalling pathways in non-syndromic sporadic brain arteriovenous malformations (BAVM). J Med Genet. 2018;55:675–84.
Mabuchi A, Manabe N, Haga N, Kitoh H, Ikeda T, Kawaji H, et al. Novel types of COMP mutations and genotype-phenotype association in pseudoachondroplasia and multiple epiphyseal dysplasia. Hum Genet. 2003;112:84–90.
Jackson GC, Mittaz-Crettol L, Taylor JA, Mortier GR, Spranger J, Zabel B, et al. Pseudoachondroplasia and multiple epiphyseal dysplasia: a 7-year comprehensive analysis of the known disease genes identify novel and recurrent mutations and provides an accurate assessment of their relative contribution. Hum Mutat. 2012;33:144–57.
Briggs MD, Bell PA, Pirog KA. The utility of mouse models to provide information regarding the pathomolecular mechanisms in human genetic skeletal diseases: the emerging role of endoplasmic reticulum stress (review). Int J Mol Med. 2015;35:1483–92.
Schulz J, Nüchel J, Niehoff A, Bloch W, Schönborn K, Hayashi S, et al. COMP-assisted collagen secretion--a novel intracellular function required for fibrosis. J Cell Sci. 2016;129:706–16.
Wang C, Liu G, Zhang W, Wang W, Ma C, Liu S, et al. Cartilage oligomeric matrix protein improves in vivo cartilage regeneration and compression modulus by enhancing matrix assembly and synthesis. Colloids Surf B Biointerfaces. 2017;159:518–26.
Chen FH, Thomas AO, Hecht JT, Goldring MB, Lawler J. Cartilage Oligomeric matrix protein/Thrombospondin 5 supports chondrocyte attachment through interaction with Integrins. J Biol Chem. 2005;280:32655–61.
Makitie O, Geiberger S, Horemuzova E, Hagenas L, Mostrom E, Nordenskjold M, et al. SLC26A2 disease spectrum in Sweden - high frequency of recessive multiple epiphyseal dysplasia (rMED). Clin Genet. 2015;87:273–8.
Superti-Furga A, Neumann L, Riebel T, Eich G, Steinmann B, Spranger J, et al. Recessively inherited multiple epiphyseal dysplasia with normal stature, club foot, and double layered patella caused by a DTDST mutation. J Med Genet. 1999;36:621.
Sakamoto Y, Yamamoto T, Kajino Y, Kabata T, Tsuchiya H, Miyake N, et al. Multiple epiphyseal dysplasia mimicking osteoarthritis due to acetabular dysplasia: a report of a familial case with a COMP mutation. J Orthop Sci. 2017;22:967–71.
We are grateful to the patients, their families, clinical surgeons, and genetic counselors for providing samples and clinical histories. We thank geneseeq inc. forexome sequencing technical support. We thank ekitech ltd. (Beijing) for providingmachine learning solutions.
This work was supported by the National Natural Science Foundation of China (81822030 to N.W., 81772299 to Z.W.), Beijing Natural Science Foundation (7191007 to Z.W.), Tsinghua University-Peking Union Medical College Hospital Initiative Scientific Research Program, the National Key Research and Development Program of China (No. 2018YFC0910500 to Z.W. and N.W.), and the CAMS Initiative Fund for Medical Sciences (2016-I2M-3-003 to N.W. and 2017-I2M-2-001 to Z.W.). All the sequencing, sample collection, and manuscript editing costs for this study were derived from these funds.
Ethics approval and consent to participate
The study was approved by the Ethical Review Board of the Peking Union Medical College Hospital (2012/03/14, JS-908). Written informed consent was provided by each participant. When the individuals younger than the age of 16, the written informed consent is obtained from their parents or legal guardians.
Consent for publication
Written informed consent for publication of clinical details and/or clinical images was obtained from the all the participants. When the individuals younger than the age of 18, the written informed consent is obtained from their parents or legal guardians.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Shao, J., Zhao, S., Yan, Z. et al. A novel COMP mutation in a Chinese family with multiple epiphyseal dysplasia. BMC Med Genet 21, 115 (2020). https://doi.org/10.1186/s12881-020-01040-y
- Cartilage oligomeric matrix protein (COMP)
- Multiple epiphyseal dysplasia
- Femoral head necrosis
- Whole-exome sequencing