Skip to main content

Genetic analysis of 55 cases with fetal skeletal dysplasia



Fetal skeletal dysplasia (SD) is a common congenital disability comprising a complex group of skeletal disorders with substantial clinical and genetic heterogeneity. Many of these defects are detected prenatally using ultrasound (US). However, the diagnostic accuracy of the US is limited.


We recruited 55 unrelated fetuses with US-detected skeletal anomalies and performed sequential tests using copy number variation sequencing, targeted skeletal gene panel sequencing, or whole exome sequencing. The detected variants were validated using Sanger sequencing or multiplex ligation-dependent probe amplification. We conducted breakpoint analysis and structural modeling of variants possibly involved in fetal SD.


A definitive diagnosis was achieved in 81.82% of affected fetuses (45/55). We identified chromosomal abnormalities in seven cases and 36 variants, of which 18 were novel pathogenic or likely pathogenic in 11 genes in 38 cases. De novo variants were identified in 27 cases (71.05%, 27/38), and one gonosomal mosaicism variant was found in the mother of one fetus. Our case examples demonstrated the high heterogeneity of fetal SDs and the rare fetal SD-associated challenges.


Careful clinical evaluation of fetuses with SD can guide appropriate molecular testing. Our study extends the SD-associated pathogenic variant spectrum and provides useful genetic counselling guidance and an accurate prenatal diagnosis strategy.


Fetal skeletal dysplasia (SD) is one of the most common congenital disabilities in clinical practice. It has been reported that the prevalence of fetal SD is 2.3 to 4.5 per 10.000 births [1]. SD is a group of heterogeneous disorders associated with skeletal abnormalities, including bone shape, size, and density abnormalities, manifesting as malformed limbs, chest, or skull. According to the 2019 nosology and classification of genetic skeletal disorders, there are 461 types of hereditary bone diseases, divided into 42 classes, involving 437 genes [2]. Moreover, skeletal abnormalities may occur in other multisystem syndromes [3]. Due to the high clinical and genetic heterogeneity of SD, it is challenging to make a precise diagnosis of fetal SD.

Although ultrasound (US) is still an indispensable first-line fetal malformation screening method during pregnancy, it has severe limitations owing to the complex and overlapping phenotypes of fetal SD, especially non-fatal fetal SD. The lack of an accurate prenatal diagnosis makes prenatal counselling and pregnancy management challenging. Therefore, the combination of molecular diagnosis with US screening offers the possibility of identifying the cause of SD. In this study, 55 fetuses with SD were recruited and analyzed using copy number variation sequencing (CNV-seq), targeted skeletal gene panel sequencing, or whole exome sequencing (WES). A qPCR walking strategy and long-range PCR were subsequently performed to identify the breakpoint junction of a novel deletion. Then, in silico prediction of the functional impact of the rare novel variant in GNAI3 was conducted. The present study aimed to enrich the phenotype and genotype of fetal SD, provide a definitive prenatal diagnosis for 55 affected families to improve genetic counselling, reduce the risk of recurrence, and characterize the molecular strategies to diagnose fetuses with SD.


Fetuses and sample collection

This study included 55 fetuses diagnosed with skeletal abnormalities with or without other systemic abnormalities using US at the First Affiliated Hospital of Zhengzhou University between January 2019 and March 2021. The fetuses were not exposed to tobacco, alcohol, radiation, or infectious diseases during pregnancy. We collected fetal abortion tissues (46/55), amniocytes (9/55), and peripheral blood samples from 96 parents. Genomic DNA was obtained using QIAamp DNA Blood Mini Kit 250 (Qiagen, Hilden, Germany) according to standard extraction methods. The present study was approved by the ethics committee (KS-2018-KY-36) of the First Affiliated Hospital of Zhengzhou University, and informed consents were obtained from all the parents of the fetuses.

Quantitative fluorescent (QF)-PCR and CNV-seq

QF-PCR and CNV-seq were performed as previously described [4, 5]. The identified CNVs were then mapped to the reference genome GRCh37/hg19. According to previously reported guidelines, we interpreted and classified the clinical significance of candidate CNVs using databases of Genomic Variants (, Database of Chromosomal Imbalance and Phenotype in Humans Using Ensembl Resources (, ClinVar (, Online Mendelian Inheritance in Man (OMIM) database ( and ClinGen (

Next generation sequencing (NGS)

For fetuses 11–25 (2019.1–2019.6), a special NGS gene panel analysis (including 811 congenital skeletal anomalies genes, Additional file 1: Table S1) was performed by a commercial company (MyGenostics, Inc., Beijing, China). For fetuses 8–10 and 26–55 (2019.7–2021.3), WES was performed and enriched for exonic sequences using the IDT xGen Exome Research Panel V1.0 (Integrated DNA Technologies, San Diego, USA). Sequencing results were mapped to the human genome (GRCh37/hg19). All variants were annotated using databases including the 1000 genomes project (1000G,, dbSNP (, the genome aggregation database (GnomAD,, ClinVar, the Human Genomic Mutation Database (HGMD, and OMIM. Variants were classified according to the guidelines recommended by the American College of Medical Genetics and Genomics (ACMG).

Verification of gene variants

Pathogenic or likely pathogenic variants were verified using Sanger sequencing. Multiplex ligation-dependent probe amplification (MLPA) was performed to detect the copy number of the EVC2 gene in case 49 and family members using the P456 probe mix (MRC-Holland, Amsterdam, The Netherlands) on an ABI PRISM 3100 Genetic Analyzer according to the manufacturer’s instructions. The obtained data were analyzed using the Coffalyser software.

Breakpoint analysis

To identify the junction of the gross deletion in case 49, a qPCR walking strategy and long-range PCR were performed. Multiple primers were designed to be regularly spaced every ~ 25 kb (step 1), ~ 5 kb (step 2), and ~ 1.5 kb (step 3) within the breakpoint regions (Fig. 4, Additional file 1: Table S2). qPCR analysis was performed using the FastStart Universal SYBR Green Master kit (Roche Applied Science, Germany) in ABI QuantStudio 5 according to the manufacturer’s instructions and the 2Ct method, as previously described [6]. Long-range PCR was performed with the forward primer of D131 and the reverse primer of U4 using Phusion High-Fidelity DNA polymerases (Thermo Scientific, Schwerte, Germany). PCR products were analyzed using electrophoresis in 1% agarose gel, and the subsequent sequences were aligned to the reference genome hg19 using the BLAT tool of the UCSC Genome Browser ( to identify the breakpoint. Finally, the 200 bp reads upstream and downstream of the breakpoints were analyzed to search for repetitive elements using the “RepeatMasker” program (UCSC).

Conservation analysis and molecular model

Protein sequencing data used in the conservation analysis was downloaded from the NCBI HomoloGene database ( The UGENE software was used for the visualization of multiple sequence alignments. A structural model of the novel variant in GNAI3 was constructed using the homology-modeling server SWISS-MODEL ( and visualized using PyMOL software, as previously described [7].


Clinical features of cases

Among the 55 families in the present study, the median age of pregnant women was 28 years, and the median gestational age was 23 weeks. According to the clinical characteristics detected using US, fetuses were divided into three groups: 10 (18.18%) had multiple malformations, 41 (74.55%) exhibited short limb deformities (short long bones less than -2standard deviation) with or without long bone bending, fracture, “telephone receiver-shaped” changes, or thoracic hypoplasia, and 4 (7.27%) isolated SD (hypomineralization, premature closure of cranial suture, small mandible, and syndactyly). Detailed information about the clinical features is provided in Additional files 1 and 2: Table S3 and Fig. S1–S12.

Genetic analysis

A total of 55 fetuses with suspected SD were studied, and 45 cases were diagnosed with hereditary diseases, with a diagnostic rate of 81.82% (45/55) (Fig. 1, Table 1). CNV-seq was performed for 1–10 fetuses with multiple malformations, and chromosomal abnormalities were identified in seven fetuses (70%, 7/10), including cases 1–5 with trisomy-18, case 6 with pathological microdeletions and microduplications (Xp22.33p22.12(2700000–19680000) × 1; 11p15.5p15.4 (180000–7100000) × 3; 20q13.2 (50380000–50640000) × 1), and case 7 with triploidy (69-XXX).

Fig. 1
figure 1

Flowchart showing the methods and results of 55 cases enrolled in our study. QF-PCR: quantitative fluorescent PCR; CNV-seq: copy number variation sequencing; Panel: a special NGS gene panel analysis (including 811 congenital skeletal anomalies genes); WES: whole exome sequencing

Table 1 Summary of main malformation classification and molecular results for 55 cases with skeletal abnormalities

To provide a molecular diagnosis, targeted skeletal gene panel sequencing accompanied by CNV-seq was applied to fetuses 11–25, and WES to fetuses 8–10 and 26–55. Thirty-eight cases were diagnosed with hereditary diseases, yielding a diagnostic rate of 79.17% (38/48). Thirty-six causative variants were detected in 11 bone development-related genes, including FGFR3 (15/38), COL1A1/COL1A2(9/38), DYNC2H1(5/38), COL2A1(3/38), ALPL (1/38), EVC2(1/38), FGFR2(1/38), GNAI3(1/38), NPR2(1/38) and RMRP (1/38) (Fig. 2, Table 1). Of these, 18 were novel variants. De novo mutations (DNMs) were identified in 27 (71.05%, 27/38) cases with autosomal dominant (AD) inheritance, compound heterozygotes were identified in eight (21.05%, 8/38) cases with autosomal recessive inheritance, and in the remaining three (7.9%, 3/38) cases with AD inheritance, the variants were inherited from the mother. Detailed information is provided in Table 1.

Fig. 2
figure 2

Genes in the present study. Eleven genes were diagnosed in 38 fetuses, including FGFR3(15/38), COL1A1/COL1A2(9/38), COL2A1(3/38), DYNC2H1(5/38), ALPL (1/38), EVC2(1/38), FGFR2(1/38), GNAI3(1/38), NPR2(1/38) and RMRP (1/38)

Case examples

Case 42

Case 42 was referred to our hospital because the fetus showed short limbs and fetal hydrops evaluated at 14 weeks and 4 days using US examination (Fig. 3A). Pregnancy was terminated at 17 weeks of gestation. The couple’s second pregnancy was evaluated at 22 weeks for suspected abnormal development of fetal limbs (micromelia; US image unavailable) and terminated. WES analysis confirmed that case 42 was heterozygous for the COL2A1 (NM_001844.4) c.2303G > T variant, resulting in a protein change: p.G768V, consistent with severe type 2 collagenopathy. Sanger sequencing confirmed that the fetus carried the mutation, which was also detected in the mother (the peak of the T allele was much lower than that of the G allele) (Fig. 3E, F). Further detailed clinical assessment of the mother at age 25 revealed mild symptoms compared with those of her fetus, including short stature (height 152 cm, father 175 cm, mother 162 cm), slightly shorter length of the left lower limb than that of the right, and a protruding short fourth toe on her left foot (Fig. 3C). In addition, she showed visual abnormalities, including severe myopia. When she was 1 year old, her left lower limb was shorter than her right lower limb by 2 cm. Her condition seemed to worsen, and she had difficulty walking due to genu varus of the left knee at age 7. Physical examination at age 16 revealed left lower extremity varus deformity, shortening (the length of the right lower limb was 66.8 cm while that of the left lower limb was 58.6 cm), internal rotation deformity of the left leg, and pelvic inclination (Fig. 3B). The patient presented limb claudication and had undergone femoral osteotomy to correct the 8-cm shortening of the left lower extremity. Due to the asymmetrical phenotype and two similar experiences of abortion, we supposed that the mother was affected by mosaicism. For an accurate diagnosis, WES was applied to the mother (II:2), and Sanger sequencing was performed on her parents. The WES results exhibited the mutated allele was present in 11 of 52 (21.2%) reads (Fig. 3G), consistent with the results of Sanger sequencing. COL2A1 c.2303G > T was not detected in her parents and has not yet been recorded in any database. The novel de novo pathogenic variant (COL2A1 c.2303G > T) was highly conserved across different species (Homo sapiens, Pan troglodytes, Macaca mulatta, Canis lupus familiaris, Bos taurus, Mus musculus, Rattus norvegicus, Gallus gallus, Danio rerio, and Xenopus tropicalis) (Fig. 3H).

Fig. 3
figure 3

COL2A1 variant caused severe type 2 collagenopathy in case 42 and mother (II2). A The image of US examination of case 42 revealed short limbs. B, C Clinical photograph of the mother with her pelvis inclined on the left side, shortened left lower limb (before extension osteotomy), and shortened fourth toe in the left foot (arrow). D The pedigree map of case 42. E, F Sanger sequencing of case 42 and the mother (II2). G The variant allele fraction and depth of coverage of the COL2A1 pathogenic variant detected by WES in case 42 and the mother. H The variant (COL2A1 c.2303G > T) was highly conserved across different species

Case 49

Due to positive pregnancies of fetal SD, case 49 came to our attention. The pregnancy was terminated because of the abnormal US at 12 weeks of gestation, presenting shortened and curved long bones (Fig. 4A). Two years ago, a scan at 23 weeks of gestation of the first affected fetus exhibited short long bones (< 5standard deviation), polydactylism, and a “common” atrioventricular  valve ring with the absence of the cardiac crux. Compound heterozygous variants (c.871-2A > G and whole EVC2 gene deletion) in EVC2 were found in case 49 using WES, inherited from the mother and father, respectively, which were validated using Sanger sequencing and MLPA analysis (Fig. 4C–F). Both variants were PVS1 + PM2-supporting + PM3 and classified as “pathogenic” according to the ACMG guidelines. Collectively, the affected fetuses were diagnosed with Ellis-van Creveld syndrome (EVC, MIM#225500). WES-based CNV analysis revealed heterozygosity at telomeric SNP rs2301856 (chr4:5461801) and centromeric SNP rs899691 (chr4:5754544), suggesting that the deleted region was smaller than 290 kb. To determine the precise size and breakpoints of the gross deletion, we estimated the number of copies (one or two) using qPCR analysis. We performed long-range PCR at 1.5-kb consecutive intervals to amplify the upstream and downstream regions of the EVC2 gene of the proband. The strategy and results are presented in Fig. 4G, H. Long-range PCR was performed for the affected fetus using the forward primer of D131 and reverse primer of U4, generating a PCR product of 1.5 kb. The aligned sequences revealed a deletion of 205 kb (Chr4:5,548199–5753996) involving the entire EVC2 gene and a part of the EVC gene. To investigate the potential underlying mechanism, the sequences between the two missing boundaries were aligned, and no microhomologous sequence was found. After searching for repetitive elements in the vicinity of the breakpoint, the telomeric breakpoint was found to be in a long interspersed nuclear element (LINE, L1PA4) belonging to the L1 family (Fig. 4I).

Fig. 4
figure 4

EVC2 variants caused EVC in case 49. A The image of US examination of case 49 revealed short and curved long bones. B The pedigree map of case 49. C MLPA results of EVC2 in case 49: a heterozygous deletion of EVC2-22 ~ EVC-8. DF Sanger sequencing showed that the proband had a novel variant inherited from the mother. G Schematic presentation of the gross deletion. Arrows indicate genes; Solid horizontal lines indicate retained regions; Broken horizontal lines indicate deleted regions; Multiple primers were designed to be regularly spaced every ~ 25 kb (step 1), ~ 5 kb (step 2), and ~ 1.5 kb (step 3) within breakpoint regions. Long-range PCR with the forward primer of the D131 and the reverse primer of the U4 generates a PCR product of ~ 1.5 kb. H The breakpoint sequences of case 49 and Sanger-sequence of the gross deletion. I UCSC Genome Browser information around the telomeric breakpoint regions (chr4: 5,547,999–5,548,399)

Case 50

Case 50 was admitted to our hospital at 13 weeks of gestation and was diagnosed as micrognathia by US (mandibular length: 8 mm, biparietal diameter: 27 mm). Besides, oligohydramnios was found in case 50 (the single deepest pocket of amniotic fluid was 40 mm). WES demonstrated a heterozygous variant, c.119G > T (p.G40V), in the GNAI3 gene. The verification results of the parents suggested that the variant was a de novo mutation (Fig. 5B–D). The novel variant was classified as “likely pathogenic” based on the interpretation of PM1 + PM2-supporting + PM5 + PM6 + PP3 according to the ACMG guidelines. Structural modeling showed that glycine 40 was located in the random coil structure of the protein secondary structure and formed a hydrogen bond with lysine 46, which possessed a bond length of 2.9 Å (Fig. 5E). After changing to valine, the hydrogen bond was slightly lengthened to 3.3 Å (Fig. 5F). Glycine 40 of Gnai3 is highly conserved among species (Homo sapiens, Pan troglodytes, Macaca mulatta, Canis lupus familiaris, Bos taurus, Mus musculus, Rattus norvegicus, Gallus gallus, Danio rerio, and Xenopus tropicalis) (Fig. 5H). In line with the aforementioned results, the etiology of case 50 was attributed to the mutation in GNAI3.

Fig. 5
figure 5

GNAI3 variant caused auriculocondylar syndrome (ACS) in case 50. A The pedigree map of case 50. BD Sanger sequencing showed that the proband had a de novo variant. E, F The 3D molecular structure of Gnαi3. The magnified views of the wild-type Gly40 (E) and mutant Val40 (F) are shown respectively. The H-bonds are shown as green dashed lines, and H-bond distances (Å) are shown in red numbers G Gnαi3 domains are depicted in blue: boxes G1–G5. The variant reported here is in red, and previously described variants are in black. H Protein alignment showing conservation of residues GNAI3 p.Gly40Val across multiple species. This mutation occurred at evolutionarily conserved amino acid positions


In this study, molecular testing of 55 cases yielded a total detection rate of 81.82% (45/55), which was slightly higher than 76.67% (23/30), 64% (35/55), 24% (8/34), and 15% (10/65) reported in previous studies on fetal bone abnormalities in recent years [3, 8,9,10]. The detection rate is affected by many factors, such as the cohort size, study inclusion criteria, and methods. The discrepancy in the detection rates may be explained by sample differences. In our study, most of the samples were fetal abortion tissue samples, which reflected the severity of the phenotypes. The diagnosis involved seven cases with chromosomal abnormalities and 38 with gene variants. Notably, the diagnostic rate of chromosomal abnormalities in fetuses with multiple malformations in our study was as high as 70% (7/10), comparable to 50% (1/2) and 62.5% (5/8) reported by Jin Han et al. [11] and Liu et al. [8] respectively. Nonetheless this observation must to be validated in a larger cohort. In addition, compared with fetuses with only skeletal malformation, cases with multiple malformations had the highest karyotype and chromosomal microarray analysis (CMA) detection rates according to Fu et al. [12] (karyotype 4.9% vs. 26.9%, and CMA 5.9% vs. 12.6%). These findings suggest that CNV analysis is an accurate and cost-effective method for identifying the etiology of SD, especially in affected fetuses with multiple malformations. Variants, of which 18 were novel, were detected in 38 cases in this study, at a rate of 79.17% (38/48). These results enrich the spectrum of genetic variants associated with fetal SD. The rate of DNMs was 71.05% (27/38), while a rate of 60% (21/35) was reported in a related study [9], which re-emphasized the importance of prenatal diagnosis and gene analysis for SD-afflicted fetuses.

The most common nosology groups in our study were the FGFR3 chondrodysplasia group (15/38, 39%), osteogenesis imperfecta and decreased bone density group (9/38, 24%), and ciliopathies with major skeletal involvement (5/38, 13%), consistent with the results of other studies [13, 14]. Furthermore, FGFR3 variants are the most common cause of fetal SD. In our cohort, cases with the FGFR3 c.1138G > A (6/38, 15.79%) and c.742C > T (6/38, 15.79%) variants accounted for one-third of the total cases.

Mosaicism may have been more common than previously estimated. In a 2019 study by Cao et al. [15], disease-related mosaic variants were present in 1.5% of probands, and 0.3% were parental mosaic variants in cases with a definite molecular diagnosis based on 12,000 whole-exome data. In our study, the mother of case 42 had two fetuses affected by severe SD. She carried a novel somatogonadal mosaic variant, COL2A1 c.2303G > T, which was detected using WES (21.2% mosaicism in her peripheral blood DNA) and confirmed using Sanger sequencing. It is worth noting that the abnormal skeletal symptoms in the mother were asymmetrical. Her left lower limb was shorter than her right, and a protruding short fourth toe was also present on her left foot, while her right lower limb and foot were normal. A father with somatogonadal mosaicism of COL2A1 c.1403G > A (22% mosaicism in blood DNA) exhibited average stature (height 173 cm) and had two affected fetuses with severe type 2 collagenopathy carrying the heterogynous mutation identified by Morrison et al. [16]. Yamamoto et al. described a similar phenomenon in mild type 2 collagenopathy with paternal somatogonadal mosaicism [17]. Based on these findings, we suspected that there were many complex mechanisms between the phenotype and mosaicism levels, not just the threshold effect [16]. Moreover, it has been suggested that because of the differences in the time and region of somatic mosaicism, the clinical phenotype presented individuation [18]. The intrafamilial recurrence rate in gonosomal mosaicism was high [19], so it is of great significance to detect gonosomal mosaicism and provide recommendations for performing prenatal diagnosis in subsequent pregnancies.

Based on positive family history and the US findings, the clinical diagnosis of EVC syndrome can also be made for the first affected fetus of family 49 [20]. EVC syndrome is a chondral and ectodermal dysplasia characterized by short stature, thoracic hypoplasia, dysplastic fingers, teeth, nails, and congenital heart disease (occurring in approximately 60% of cases), especially atrioventricular septal defect and single atrium [21]. In addition, patent ductus arteriosus in fetuses with EVC syndrome has previously been reported [20]. A common atrioventricular valve with the absence of the cardiac crux was found at 22 weeks of gestation in the first affected fetus of family 49 using the three-dimensional US. However, the second affected fetal intracardiac structures could not be revealed by the US at 12 weeks of gestation. Our findings on the affected fetuses were consistent with the existing literature. A novel compound heterozygous variant (c.871-2A > G) and a gross 205 kb deletion in EVC2 were identified in case 49. Finally, we sought to define the boundaries of the deletion and decipher the underlying mechanism. Using a qPCR walking strategy to sequence the breakpoint junction, the telomeric breakpoint was considered to be derived from repeating elements that overlapped with LINE (L1PA4). Rearrangement of breakpoints may result from distinct mutational mechanisms: (1) homologous recombination, (2) non-homologous end joining, (3) microhomology-mediated replication-dependent recombination, (4) LINE-1-mediated retrotransposition, and (5) telomere healing [22]. LINE-1 sequences constitute a significant factor in genome arrangements [23]. These regions may be fragile and prone to chromosomal abnormalities.

Previous studies have implicated variants of the GNAI3 gene in auriculocondylar syndrome (ACS, MIM#602483), a rare AD disorder characterized by micrognathia, mandibular condyle hypoplasia, and question mark auricular malformation [24]. A novel de novo variant, c.119G > T (p.G40V) in GNAI3, was identified in case 50 using WES analysis because of “micrognathia” at 13 weeks of gestation detected using US examination. To the best of our knowledge, this is the fourth fetus afflicted with ACS with a GNAI3 mutation. Liu et al. summarized the clinical features of the other three fetuses, suggesting that severe micrognathia and mandibular hypoplasia accompanied by polyhydramnios were prenatal indicators of ACS [24]. Amniotic fluid was normal in case 50 (amniotic fluid volume was 40 mm) in the present study. A possible explanation for this is the incomplete expression of the phenotype in the fetus in the first trimester. In the future, prenatal characteristics must be identified in more affected fetuses to promote ACS awareness. To date, only seven variants of GNAI3 have been reported [24], and all eight variants were missense, including p.Gly40Val reported in the present study (Fig. 5G). Gnai3, encoded by GNAI3, consists of five highly conserved motifs (G1–G5) that are involved in the rat sarcoma superfamily and G α protein binding to guanosine diphosphate (GDP)/guanosine triphosphate (GTP) [25]. Most variants (7/8, 87.5%) were located in the G1 or G4 box, and the p.Gly40Val variant was located in the G1 box (Fig. 5G). Although the p.Gly40Val variant did not disrupt the overall structure of the protein as demonstrated using structural modeling, it may indirectly disrupt GDP/GTP binding [26].

Our study has some limitations. First, our sample size was relatively small, and most tissue samples were derived from fetal abortion. Thus, the results of a larger sample size would be more objective and accurate regarding the diagnosis rate of fetal SDs. Second, different tissues of the mother of case 42 were not collected to explore the relationship between the mosaicism level and asymmetrical phenotype. Third, because case 50 was detected in the first trimester, we could not obtain more useful clinical information, except for micrognathia. In addition, further functional studies should be performed to analyze the effects of these variants, especially those of the novel and VUS variants.


In summary, the detection rate in our study was 81.82% (45/55), and the detection rate of CNV analysis in multiple malformations was 70% (7/10), indicating that CNV analysis is essential in SD fetuses with multiple malformations. In addition, 38 different variants were identified, of which 18 were novel. Moreover, we demonstrated gonosomal mosaicism in a patient with asymmetrical symptoms and confirmed the breakpoints of a novel gross deletion. Our study is informative for genetic counselling and future pregnancies of affected families and extends the clinical features and genetic spectrum of fetal SD.

Availability of data and materials

The original data belonging to this study are available from the corresponding author upon request.



Skeletal dysplasia




Copy number variation sequencing


Whole exome sequencing


Quantitative fluorescent-PCR


Online Mendelian inheritance in man


Next generation sequencing


1000 genomes project

G nomAD:

Genome aggregation database


Human genomic mutation database


American College of Medical Genetics and Genomics


Multiplex ligation-dependent probe amplification


De novo mutations


Autosomal dominant


Ellis-van creveld syndrome


Long interspersed nuclear element


Chromosomal microarray analysis


Auriculocondylar syndrome


Guanosine diphosphate


Guanosine triphosphate


  1. Yang K, Shen M, Yan Y, et al. Genetic analysis in fetal skeletal dysplasias by trio whole-exome sequencing. Biomed Res Int. 2019;2019:2492590.

    PubMed  PubMed Central  Google Scholar 

  2. Mortier GR, Cohn DH, Cormier-Daire V, et al. Nosology and classification of genetic skeletal disorders: 2019 revision. Am J Med Genet A. 2019;179(12):2393–419.

    Article  PubMed  Google Scholar 

  3. Lord J, McMullan DJ, Eberhardt RY, et al. Prenatal exome sequencing analysis in fetal structural anomalies detected by ultrasonography (PAGE): a cohort study. Lancet. 2019;393(10173):747–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Sun Y, Kong X, Zhao Z, et al. Mutation analysis of 419 family and prenatal diagnosis of 339 cases of spinal muscular atrophy in China. BMC Med Genet. 2020;21(1):133.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Li Q, Zhu X, Wang C, et al. Identification of a rare case with nagashima-type palmoplantar keratoderma and 18q deletion syndrome via exome sequencing and low-coverage whole-genome sequencing. Front Genet. 2021;12:707411.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Bai Y, Chen Y, Kong X. Contiguous 22.1-kb deletion embracing AVPR2 and ARHGAP4 genes at novel breakpoints leads to nephrogenic diabetes insipidus in a Chinese pedigree. BMC Nephrol. 2018;19(1):26.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Fang X, Zhu C, Zhu X, et al. Molecular analysis and novel variation identification of Chinese pedigrees with mucopolysaccharidosis using targeted next-generation sequencing. Clin Chim Acta. 2022;524:194–200.

    Article  CAS  PubMed  Google Scholar 

  8. Liu Y, Wang L, Yang YK, et al. Prenatal diagnosis of fetal skeletal dysplasia using targeted next-generation sequencing: an analysis of 30 cases. Diagn Pathol. 2019;14(1):76.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Zhang L, Pan L, Teng Y, et al. Molecular diagnosis for 55 fetuses with skeletal dysplasias by whole-exome sequencing: A retrospective cohort study. Clin Genet. 2021;100(2):219–26.

    Article  CAS  PubMed  Google Scholar 

  10. Petrovski S, Aggarwal V, Giordano JL, et al. Whole-exome sequencing in the evaluation of fetal structural anomalies: a prospective cohort study. Lancet. 2019;393(10173):758–67.

    Article  CAS  PubMed  Google Scholar 

  11. Han J, Yang YD, He Y, et al. Rapid prenatal diagnosis of skeletal dysplasia using medical trio exome sequencing: Benefit for prenatal counseling and pregnancy management. Prenat Diagn. 2020;40(5):577–84.

    Article  CAS  PubMed  Google Scholar 

  12. Fu F, Li R, Li Y, et al. Whole exome sequencing as a diagnostic adjunct to clinical testing in fetuses with structural abnormalities. Ultrasound Obstet Gynecol. 2018;51(4):493–502.

    Article  CAS  PubMed  Google Scholar 

  13. Barkova E, Mohan U, Chitayat D, et al. Fetal skeletal dysplasias in a tertiary care center: radiology, pathology, and molecular analysis of 112 cases. Clin Genet. 2015;87(4):330–7.

    Article  CAS  PubMed  Google Scholar 

  14. Stevenson DA, Carey JC, Byrne JL, et al. Analysis of skeletal dysplasias in the Utah population. Am J Med Genet A. 2012;158A(5):1046–54.

    Article  PubMed  Google Scholar 

  15. Cao Y, Tokita MJ, Chen ES, et al. A clinical survey of mosaic single nucleotide variants in disease-causing genes detected by exome sequencing. Genome Med. 2019;11(1):48.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Morrison PJ, Znaczko A, Jansson M. Paternal somatogonadal COL2A1 mosaicism causing recurrence of severe type 2 collagenopathy. Am J Med Genet A. 2020;182(9):2191–4.

    Article  CAS  PubMed  Google Scholar 

  17. Yamamoto K, Kubota T, Takeyari S, et al. Parental somatogonadal COL2A1 mosaicism contributes to intrafamilial recurrence in a family with type 2 collagenopathy. Am J Med Genet A. 2020;182(3):454–60.

    Article  CAS  PubMed  Google Scholar 

  18. Forzano F, Lituania M, Viassolo A, et al. A familial case of achondrogenesis type II caused by a dominant COL2A1 mutation and “patchy” expression in the mosaic father. Am J Med Genet A. 2007;143A(23):2815–20.

    Article  CAS  PubMed  Google Scholar 

  19. Campbell IM, Stewart JR, James RA, et al. Parent of origin, mosaicism, and recurrence risk: probabilistic modeling explains the broken symmetry of transmission genetics. Am J Hum Genet. 2014;95(4):345–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Chen CP, Su YN, Hsu CY, et al. Ellis-van Creveld syndrome: prenatal diagnosis, molecular analysis and genetic counseling. Taiwan J Obstet Gynecol. 2010;49(4):481–6.

    Article  CAS  PubMed  Google Scholar 

  21. Ozturk O, Bagis H, Bolu S, et al. Ellis-van Creveld syndrome novel pathogenic variant in the EVC2 gene a patient from Turkey. Clin Case Rep. 2021;9(4):1973–6.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Chen JM, Cooper DN, Ferec C, et al. Genomic rearrangements in inherited disease and cancer. Semin Cancer Biol. 2010;20(4):222–33.

    Article  CAS  PubMed  Google Scholar 

  23. Beck CR, Collier P, Macfarlane C, et al. LINE-1 retrotransposition activity in human genomes. Cell. 2010;141(7):1159–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Liu X, Sun W, Wang J, et al. Prenatal diagnosis of auriculocondylar syndrome with a novel missense variant of GNAI3: a case report. BMC Pregnancy Childbirth. 2021;21(1):780.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Gordon CT, Vuillot A, Marlin S, et al. Heterogeneity of mutational mechanisms and modes of inheritance in auriculocondylar syndrome. J Med Genet. 2013;50(3):174–86.

    Article  CAS  PubMed  Google Scholar 

  26. Romanelli TV, Gordon CT, Zechi-Ceide RM, et al. Novel variants in GNAI3 associated with auriculocondylar syndrome strengthen a common dominant negative effect. Eur J Hum Genet. 2015;23(4):481–5.

    Article  Google Scholar 

Download references


We deeply thank all the families for participation in the present.


This study was supported by grants from the National Natural Science Foundation of China (No. 81501851 to BY and No. 81701533 to LN) and the Henan Province Medical Science and Technique Foundation (No. 201702003 to BY).

Author information

Authors and Affiliations



BY designed the study, performed data analysis, and revised the manuscript. SY performed protein structure analysis and wrote the manuscript. HYQ, DHK, LQQ, ZXF, and MJJ performed lab experiments. LN, WL, JZH collected the blood samples and clinical data. KXD collected clinical data and revised the manuscript. All authors contributed to data interpretation, approved the final, and submitted the version of the manuscript.

Corresponding authors

Correspondence to Ying Bai or Xiangdong Kong.

Ethics declarations

Ethics approval and consent to participate

This research was approved by the ethics committee of the First Affiliated Hospital of Zhengzhou University (KS-2018-KY-36). Written informed consents were obtained from all the participants.

Consent for publication

All participants/families agree to the publication of the presented work.

Competing interests

The authors have declared no conflicting interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1. Table S1:

811 genes of a genetic skeletal disease panel; Table S2: Primer sequence in this study; Table S3: Summary of clinical information and molecular results for 55 cases with skeletal abnormalities.

Additional file 2.

Images of US examination of some fetuses in this study.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Bai, Y., Sun, Y., Liu, N. et al. Genetic analysis of 55 cases with fetal skeletal dysplasia. Orphanet J Rare Dis 17, 410 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Skeletal dysplasia
  • Molecular diagnosis
  • Genotype–phenotype analysis
  • Breakpoints
  • Novel variants