Whole-exome sequencing identifies de novo mutation in the COL1A1 gene to underlie the severe osteogenesis imperfecta
© Maasalu et al.; licensee BioMed Central. 2015
Received: 10 February 2015
Accepted: 30 April 2015
Published: 10 May 2015
Osteogenesis imperfecta (OI) comprises a clinically and genetically heterogeneous group of connective tissue disorders, characterized by low bone mass, increased bone fragility, and blue-gray eye sclera. OI often results from missense mutations in one of the conserved glycine residues present in the Gly-X-Y sequence repeats of the triple helical region of the collagen type I α chain, which is encoded by the COL1A1 gene. The aim of the present study is to describe the phenotype of OI II patient and a novel mutation, causing current phenotype.
We report an undescribed de novo COL1A1 mutation in a patient affected by severe OI. After performing the whole-exome sequencing in a case parent–child trio, we identified a novel heterozygous c.2317G > T missense mutation in the COL1A1 gene, which leads to p.Gly773Cys transversion in the triple helical domain of the collagen type I α chain. The presence of the missense mutation was confirmed with the Sanger sequencing.
Hereby, we report a novel mutation in the COL1A1 gene causing severe, life threatening OI and indicate the role of de novo mutation in the pathogenesis of rare familial diseases. Our study underlines the importance of exome sequencing in disease gene discovery for families where conventional genetic testing does not give conclusive evidence.
KeywordsOsteogenesis imperfecta Type I collagen OI genotype–phenotype COL1A1 De novo mutation
Osteogenesis Imperfecta (OI), or “brittle bone” disease, is a heritable disorder of collagen type I metabolism with a generalized involvement of connective tissues. Collagen is the most abundant protein in mammals, constituting a quarter of the total protein weight . Collagens are grouped into families based on their structural and functional features. Type I collagen is the major protein in bone, skin, tendon, ligament, sclera and cornea tissues, blood vessels, and hollow organs ENREF 2 . OI is mostly caused by quantitative or qualitative collagen type I defects. The condition is characterized by low bone mass, bone fragility, and often short stature. Extraskeletal manifestations may include blue-gray eye sclera and dental abnormalities. The clinical severity varies widely from nearly asymptomatic forms with a mild predisposition to fractures, normal stature, and normal lifespan to profoundly disabling and even lethal [3, 4].
The pathogenetic approach to OI is changed with the recent identification of non-collagenous genes, mutations in which may cause OI. In general, a clear genotype-phenotype correlation does not exist. General rules for genotype-phenotype correlations have been published only in COL1A1/2-related OI . Approximately 90 % of individuals affected with OI are heterozygous for a causative variant in one of the two genes, COL1A1 or COL1A2, which encode the pro-1(I) and pro-2(I) chains of type I procollagen, respectively .
The proportion of cases caused by a de novo COL1A1 or COL1A2 mutation varies according to the severity of the disease. Approximately 60 % of cases of classic non-deforming OI with blue sclerae or common variable OI with normal sclerae, virtually 100 % of perinatally lethal OI, and close to 100 % of progressively deforming OI are caused by de novo mutations [7, 8].
In 1979, a classification of OI was introduced by David Sillence and the disease was divided into four types with a wide spectrum of clinical features, where OI type II is the most severe and prenatally lethal form of the disorder . Initial classification, based on clinical and histological manifestations, was extended into five distinct types of OI . By now, genetic studies described various OI phenotypes, and genetic OI classification is broadened up to 15 different OI types, according to the affected gene. In addition to collagen genes, OI is caused by mutations in the CRTAP (OI VII), LEPRE 1 (OI VIII), BMP1 (OI XIII), TMEM38B (OI XIV), IFITM5 (OI V), SERPINH1 (OI X), WNT1 (OI XV), SP7 (OI XII), PPIB (OI IX), SERPINF1 (OI VI), FKBP10 (OI XI) genes [11–14].
In addition to big genotype diversity, the phenotype manifestations may vary not only among representatives of the same OI type but also among carriers of the same mutation and even affected members of the same family. Some of the OI forms represent transitional phenotypes between different forms. Therefore, the genotype-phenotype relationships and clinical manifestations of OI are still difficult to explain.
The aim of present study is to describe the phenotype of OI II patient, with severe life threatening bone fragility, and report a novel mutation, causing current phenotype. We strongly believe that new additional information on current transitional OI form will deepen the understanding of the pathogenesis of the brittle bone disease.
Materials and methods
The proband (716) was the first child (first pregnancy, first delivery) of non-consanguineous Estonian parents; age of the mother and father were 23 and 27 years, respectively. The parents were healthy without history of chronic or clinically significant diseases.
The girl was born with totally soft skull, the head was large (diameter 43 cm), and bones in occipital region were not palpable (skinhead). Her head was flat, being collapsed from front to back direction. She has exophthalmia and blue sclerae. The newborn had disproportional growth retardation (limbs > body): short, bowed arms and legs, and both hips were hyperflexed and turned outward.
Total skinhead and only partially developed upper part of parietal and mandibular bones were detected in Babygram X-ray investigation on date of birth. All long bones were extremely osteopenic; she had accordion-type ribs, and several fractures in different healing stages were detected in every bone. In addition, fresh right humeral fracture and left tibial fracture were confirmed.
In genetic counseling, at the age of 3 days, the baby was diagnosed with osteogenesis imperfecta type II based on clinical signs and X-ray data. As family history was negative for the OI, de novo autosomal dominant mutation was suspected.
Whole-exome sequencing was performed on an affected child and both unaffected parents at the NGS core facility of the Estonian Genome Center, University of Tartu. Exome capture was performed using the TruSeq Exome Enrichment kit (Illumina) following the manufacturer’s protocol. The captured libraries were sequenced with Illumina HiSeq2000 with 100-bp paired-end reads. Over 10 Gb of sequence was generated from each individual, resulting in a coverage depth of 84× for both parents and 87× for an affected child. Sequence reads were aligned to the human reference genome (hg19, GRCh37) with the Burrows-Wheeler Aligner (BWA, version 0.6.1) . Single-nucleotide substitutions and small indel variants were called with SAM tools (version 0.1.18), Picard tools (version 1.60), and a Genome Analysis Toolkit (GATK, version 1.5.21) [16, 17]. Genotypes were called at all positions with high-quality sequence bases and filtered to retain SNPs and insertion-deletions with Phred-like quality scores of at least 20. We focused on non-synonymous and canonical splice-site variants being absent from public datasets (including dbSNP135 and the 1000 Genomes Project) and in-house exome and full-genome data. We used the PolyPhen-2, SIFT, and Condel software tools to predict the functional effects of mutations [18–20]. Mutation analysis was performed with Sanger sequencing on an affected child (716), on both parents (710, 711), and an unaffected brother (715).
Here, we report a heterozygous p.Gly773Cys mutation in COL1A1, affecting a highly conserved residue in the triple helical domain of collagen type I α chain, as a novel causative variant for severe OI with the lethal outcome. Therefore, given the large number of different genes responsible for OI forms, the genotype-phenotype relationships and clinical manifestations are difficult to explain and present report gives additional information of pathophysiology of OI.
The severity of clinical signs (blue sclerae, totally soft, large and collapsed head; disproportional growth retardation, short and bowed arms and legs, severe breathing failure) and clinical findings (extremely osteopenic bones, several fractures in different healing stages in all bones, deformities, bone development retention) of our patient resembles those of the patients described with the lethal type of OI.
The p.Gly773Cys is located within a long stretch of helical residues 691–823, in which essentially mostly lethal substitutions are identified and also corresponds to Bodian’s “COL1A1 decision tree” about the lethality of mutation [8, 21, 22]. This region correlates with the major ligand binding region 2 (MLBR 2) that extends from helical positions 682–830 in the type I collagen fibril and includes sites that are crucial for collagen self-assembly . This region is important also for interactions of collagen monomers or fibrils with α1β1/α2β1 integrins, matrix metalloproteinases (MMPs), fibronectin, and cartilage oligomeric matrix protein . The Gly773 residue lies within overlapping MMP interaction domain and cell interaction domain with human type I collagen fibril and within ligand binding sites for secreted protein, acidic, and rich in cysteine (SPARC), fibronectin, and MMP1, 12, and 13 . Therefore, mutations in this region interfere severely with the function of protein.
Substitutions of the same glycine residue in COL1A1 often have independent lethal and non-lethal outcomes. The described p.Gly773Cys mutation was lethal, excluding possibility that this family manifests OI type III. Moreover, neighboring Gly to Cys substitutions are reported to have diverse outcomes. P.Gly770Cys was reported to cause OI type II with a lethal outcome, whereas p.Gly788Cys represents a non-lethal variant causing OI type III/IV .
The present study describes a de novo p.Gly773Cys mutation in COL1A1 related to the severe OI that gives additional information of pathophysiology of OI. This shows that phenotyping together with genotyping is important to identify special patients and relevant for counseling of parents.
Our study underlines the importance of exome sequencing in disease gene discovery for families where conventional genetic testing does not give conclusive evidence.
This study was supported by the Estonian Science Agency project IUT20-46 (TARBS14046I), the European Regional Development Fund and Archimedes Foundation to the Centre of Excellence on Translational Medicine, the Targeted Financing from the Estonian Ministry of Education and Research (grant SF0180142s08), the HypOrth Project funded by the European Union 7th Framework Programme grant agreement n° 602398, the EU FP7 grant BBMRI-LPC (#313010), the Development Fund of the University of Tartu (grant SP1GVARENG), and the grant 3.2.0304.11-0312 funded by the European Regional Development Fund to the Centre of Excellence in Genomics (EXCEGEN).
- Kielty CM, Grant ME. Connective Tissue and Its Heritable Disorders. Hoboken, NJ, USA: John Wiley & Sons, Inc; 2002. p. 159–221.View ArticleGoogle Scholar
- Chan T-F, Poon A, Basu A, Addleman NR, Chen J, Phong A, et al. Natural variation in four human collagen genes across an ethnically diverse population. Genomics. 2008;91:307–14.View ArticlePubMed CentralPubMedGoogle Scholar
- Ben Amor IM, Glorieux FH, Rauch F. Genotype-phenotype correlations in autosomal dominant osteogenesis imperfecta. J Osteoporos. 2011;2011:540178.PubMed CentralPubMedGoogle Scholar
- Amor IMB, Rauch F, Gruenwald K, Weis M, Eyre DR, Roughley P, et al. Severe osteogenesis imperfecta caused by a small in-frame deletion in CRTAP. Am J Med Genet A. 2011;155A:2865–70.View ArticlePubMedGoogle Scholar
- Ben Amor M, Rauch F, Monti E, Antoniazzi F. Osteogenesis imperfecta. Pediatr Endocrinol Rev. 2013;10(2):397–405.PubMedGoogle Scholar
- Van Dijk FS, Cobben JM, Kariminejad A, Maugeri A, Nikkels PGJ, van Rijn RR, et al. Osteogenesis imperfecta: a review with clinical examples. Mol Syndromol. 2011;2:1–20.PubMed CentralPubMedGoogle Scholar
- Byers PH, Steiner RD. Osteogenesis imperfecta. Annu Rev Med. 1992;43:269–82.View ArticlePubMedGoogle Scholar
- Marini JC, Forlino A, Cabral WA, Barnes AM, San Antonio JD, Milgrom S, et al. Consortium for osteogenesis imperfecta mutations in the helical domain of type I collagen: regions rich in lethal mutations align with collagen binding sites for integrins and proteoglycans. Hum Mutat. 2007;28:209–21.View ArticlePubMed CentralPubMedGoogle Scholar
- Sillence DO, Senn A, Danks DM. Genetic heterogeneity in osteogenesis imperfecta. J Med Genet. 1979;16:101–16.View ArticlePubMed CentralPubMedGoogle Scholar
- Rauch F, Glorieux F. Osteogenesis imperfecta. Lancet. 2004;363(9418):1377–85.View ArticlePubMedGoogle Scholar
- Switch gene. Osteogenesis Imperfecta Variant Database. Leiden Open Variation Database [https://oi.gene.le.ac.uk/home.php]
- Laine CM, Wessman M, Toiviainen-Salo S, Kaunisto MA, Mäyränpää MK, Laine T, et al. A novel splice-mutation in PLS3 causes X-Linked early-onset low-turnover Osteoporosis. J Bone Miner Res. 2015;30(3):510–8.View ArticlePubMedGoogle Scholar
- Keupp K, Beleggia F, Kayserili H, Barnes AM, Steiner M, Semler O, et al. Mutations in WNT1 cause different forms of bone fragility. Am J Hum Genet. 2013;92:565–74.View ArticlePubMed CentralPubMedGoogle Scholar
- Fahiminiya S, Majewski J, Mort J, Moffatt P, Glorieux FH, Rauch F. Mutations in WNT1 are a cause of osteogenesis imperfecta. J Med Genet. 2013;50:345–8.View ArticlePubMedGoogle Scholar
- Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.View ArticlePubMed CentralPubMedGoogle Scholar
- McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297–303.View ArticlePubMed CentralPubMedGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Kumar P, Henikoff S, Ng PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009;4:1073–81.View ArticlePubMedGoogle Scholar
- 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:248–9.View ArticlePubMed CentralPubMedGoogle Scholar
- González-Pérez A, López-Bigas N. Improving the assessment of the outcome of nonsynonymous SNVs with a consensus deleteriousness score. Condel Am J Hum Genet. 2011;88:440–9.View ArticleGoogle Scholar
- Bodian DL, Madhan B, Brodsky B, Klein TE. Predicting the clinical lethality of osteogenesis imperfecta from collagen glycine mutations. Biochemistry. 2008;47:5424–32.View ArticlePubMedGoogle Scholar
- Bodian DL, Chan T-F, Poon A, Schwarze U, Yang K, Byers PH, et al. Mutation and polymorphism spectrum in osteogenesis imperfecta type II: implications for genotype-phenotype relationships. Hum Mol Genet. 2009;18:463–71.View ArticlePubMed CentralPubMedGoogle Scholar
- Di Lullo GA, Sweeney SM, Korkko J, Ala-Kokko L, San Antonio JD. Mapping the ligand-binding sites and disease-associated mutations on the most abundant protein in the human, type I collagen. J Biol Chem. 2002;277:4223–31.View ArticlePubMedGoogle Scholar
- Sweeney SM, Orgel JP, Fertala A, McAuliffe JD, Turner KR, Di Lullo GA, et al. Candidate cell and matrix interaction domains on the collagen fibril, the predominant protein of vertebrates. J Biol Chem. 2008;283:21187–97.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.