Among the 16 cases in the present study, 12 cases received a definitive molecular diagnosis, including a microdeletion and 8 novel variants.
In case 1, a deletion of 184 kb fragment in the Xp21.1 segment of chromosome X was found in the infant and inherited from his mother. He exhibited abnormal posturing of the lower extremities which were not reported previously. To our knowledge, only one fetus was reported and exhibited fetal growth restriction and oligohydramnios [12]. In 2010, Murugan et al found an in-frame deletion of exon 49 to 53 detected by MLPA and mPCR [13]. According to frame shift hypothesis [14], it usually caused BMD. Unfortunately, our patient was dead and we could not get more information. Therefore, our patient suffered several new symptoms in intrauterine period and expanded the phenotype spectrum of DMD/BMD. In case 2, two compound heterozygous causative mutations in DYNC2H1 gene were detected which was associated with short-rib thoracic dysplasia 3 [15]. Short rib-polydactyly syndrome 3(SRPS 3) was an autosomal recessive disease overlapping with Jeune asphyxiating thoracic dystrophy belonging to the ciliopathy but it was more severe and characterized by early prenatal expression, lethality and variable malformations [16].
Moreover, we found mutations in ALPL in case 3 and 4. Pathogenic variants in ALPL cause hypophosphatasia characterized by defective mineralization of bone and/or teeth in the presence of low activity of serum and bone alkaline phosphatase (ALP) [17]. c.984_986del in case 3 was previously reported in 2012 and 2015 [18, 19], which result in the deletion of Leucine at the β-sheet and decreased activity of its coding protein the tissue-nonspecific isoenzymes of alkaline phosphatase (TNSALP) [20]. Thus, c.984_986del may reduce the enzymatic activity, too. c.1460C > T (p.A487V) and c.1466G > C (p.C489S) were reported early, of which the latter one exhibited a diminished ALP activity, less located on the cell surface and failed to become the mature form [21, 22]. This suggested that p. 487–489 may play an important role in enzymatic activity. c.2T > C caused the absence of the start codon and generated a transcript starting at Met56 as c.3G > A did which does not exhibit enzymatic activity, has no significant effect on the wild type ALPL protein and cannot be attached to the cell membrane [23].
Patient 5 harbored c.742C > T in FGFR3 gene which has been detected in different racial types [24–27]. And patient 16 had a hot mutation c.1138G > A in FGFR3. These fetuses’ ultrasound scan all revealed a narrow chest with shortening of the long bones. Besides, c.8553del and c.12532 + 1G > T in HSPG2 was identified in case 10. HSPG2 is an essential gene and its mutations could lead to Schwartz-Jampel syndrome, type 1(SJS) and severe neonatal lethal Dyssegmental dysplasia, Silverman-Handmaker type (DDSH) [28, 29]. c.8553del produced a truncated protein that terminated at 2878 sites lacking part of domain Ⅳ and the whole domainⅤand c.12532 + 1G > T at the C-terminal region may yield abnormal transcript as predicted by Human Splicing Finder. Unfortunately, we didn’t get the parents’ DNA so that the compound heterozygous condition couldn’t be confirmed. Trio sequencing of patient 7 and the parents detected a heterozygous missense variation c.4813C > T in FLNB which was inherited from the normal father with 62 of 120 (51.7%) reads. Though several cases with family history presenting an autosomal dominant trait has been reported [30, 31], the father in case 7 didn’t have malformations associated with FLNB-Related Disorders like short stature, club feet and facial dysmorphisms [32].
Furthermore, variants related to type I collagen genes and caused osteogenesis imperfect (OI) were identified in 5 cases. Type I collagen is a heterotrimer, containing two α1(I) and one α2(I) chains assembling by procollagen chains with N-terminal and C-terminal globular propeptides flanking the helical domain [33]. The helical domains contain Gly-Xaa-Yaa triplets where glycine was substituted most frequently causing OI [34]. In our study, we detected 4 missense mutations and one splicing mutation variant: c.1921G > A and c.3389G > A in COL1A1 and c.1010G > A, c.2189G > T and c.1764 + 3_1764 + 6del in COL1A2 respectively (Fig. 2b). Four missense mutations were glycine substitutions of Gly-X-Y which would delay helical folding, prolonging access time for modifying enzymes. Former researches has described two infants with c.2188G > T and c.2188G > C in COL1A2 changing the same amino acid site with c.2189G > T and separately exhibiting the same deformities like blue sclera, wormian bones, shortening and bowing of the upper and lower limbs [35, 36]. The similar situation was occurred in c.1010G > A and c.1921G > A. And c.3389G > A in COL1A1 has been listed and defined as likely pathogenic variant in ClinVar. The splicing mutation c.1764 + 3_1764 + 6del in COL1A2 in case 8 may yield abnormal transcript as predicted by Human Splicing Finder. Notably, the variants in case 8 and 12 were maternally inherited. Both of the mothers experienced induced abortion in second trimester twice due to the same skeletal dysplasia malformation, suggesting the mosaicism. The mutated allele c.3389G > A was present in 17 of 66 (25.8%) reads in WES and the mutation-related signal was smaller than that of the wild-type allele in sanger sequencing which suggested that it is a mosaic mutation in the mother of case 12. After that, we tried to recall the asymptomatic mother, however, she didn’t received radiographical examination and we only knew that her sclera and height were normal and didn’t suffer bone fracture before. Similarly, the mutated allele c.1764 + 3_1764 + 6del was present in 22 of 68 (32.3%) reads in the mother of case 8 who presented extremely mild symptoms like short stature compared with the fetus. Nevertheless, the ratio of the fetus was 37.8% (17/45) similar to the mother. We supposed that other factors such as underlying genetic modifiers may affect the phenotypes.