A Novel Splice Variant in the COL1A1 Gene Leads to Exon 46 Skipping and Osteogenesis Imperfecta

Abstract

Background Osteogenesis imperfecta (OI) is a clinical and genetic disorder characterised by bone fragility, growth deficiency and skeletal deformity. Ninety per cent of OI cases are attributable to autosomal dominant variants in the COL1A1 and COL1A2 genes. Methods Candidate variants were identified and verified through trio whole‐exome sequencing (trio‐WES), copy number variation sequencing (CNV‐seq) and Sanger sequencing. Minigene splicing assays were performed in HeLa and HEK293T cells with pcDNA3.1 and pcMINI‐C vectors to investigate the function of the candidate variants. A systematic review of COL1A1 splicing variants and the corresponding genotype–phenotype spectrum was performed. Results Trio‐WES revealed a novel heterozygous variant in the C‐terminal region of the COL1A1 gene: NM_000088.4:c.3423+5G>A. Sanger sequencing confirmed the variant in both the proband (II‐2) and her foetus (III‐1) who were clinically suspected of having OI. The c.3423+5G>A variant causes complete skipping of Exon 46, as demonstrated by a minigene splicing assay. We retrieved 419 COL1A1 splicing variants from PubMed, excluded 15 without phenotypic data and 2 linked to Ehlers–Danlos syndrome and stratified the remaining 402 variants into three types on the basis of splice site location: (1) Variants at canonical splicing sites (77.8%, 313/402) mostly cause mild phenotypes, whereas a minority may be severe. (2) Intron variants in other locations, such as splice region variants (17.9%, 72/402), usually cause mild clinical phenotypes, and deep intronic splice variants (0.4%, 2/402) that may result in severe phenotypes. (3) Other variants (3.7%, 15/402), such as exon variants or fragment loss, are extremely rare. We also preliminarily discuss the mechanisms underlying phenotypic variability and the characteristics of C‐terminal variants. Conclusions This intron variant in COL1A1 was classified as likely pathogenic and was confirmed to disrupt COL1A1 expression. The summary analysis results also revealed a correlation among splicing variants, C‐terminal region variants and disease, suggesting that variant location provides a useful framework for prognosis prediction.