# Genetic Analysis in Fetal Skeletal Dysplasias by Trio Whole-Exome Sequencing

Authors: Kai Yang; Ming Shen; Yousheng Yan; Ya Tan; Jing Zhang; Jue Wu; Guangming Yang; Shang Li; Jing Wang; Zhuo Ren; Zhe Dong; Shan Wang; Manli Zhang; Yaping Tian
Journal: BioMed Research International
 (2019)

---

## Abstract

Skeletal dysplasias (SDs) comprise a series of severe congenital disorders that have strong clinical heterogeneity and usually attribute to diverse genetic variations. The pathogenesis of more than half of SDs remains unclear. Additionally, the clinical manifestations of fetal SDs are ambiguous, which poses a big challenge for accurate diagnosis. In this study, eight unrelated families with fetal SD were recruited and subjected to sequential tests including chromosomal karyotyping, chromosomal microarray analysis (CMA), and trio whole-exome sequencing (WES). Sanger sequencing and quantitative fluorescence PCR (QF-PCR) were performed as affirmative experiments. In six families, a total of six pathogenic/likely pathogenic variations were identified in four genes includingSLC26A2, FGFR3, FLNB,andTMEM38B. These variations caused disorders following autosomal dominant or autosomal recessive inheritance patterns, respectively. The results provided reliable evidence for the subsequent genetic counseling and reproductive options to these families. With its advantage in variation calling and interpreting, trio WES is a promising strategy for the investigation of fetal SDs in cases with normal karyotyping and CMA results. It has considerable prospects to be utilized in prenatal diagnosis.


---

## Body

## 1. Introduction

Skeletal dysplasias (SDs), a series of heterogeneous genetic disorders affecting approximately 2.3 to 4.5 of 10,000 births [1–3], are often hereditable and affect the growth, morphometry, and integrity of cartilage and/or bone. SDs are individually rare, but collectively they comprise a large group of disorders ranging from relatively mild anomalies to lethality. According to the 2015 Nosology and Classification of Genetic Skeletal Disorders [4], 436 genetic skeletal disorders were classified into 42 groups, associating with one or more of 364 genes. However, only a small part of these disorders has clear molecular pathogenesis [5]. Moreover, skeletal involvement may also occur in other multisystem syndromes [6]. Therefore, due to the clinical and genetic heterogeneity of SDs, it is very challenging to make a clear diagnosis, particularly in the prenatal diagnosis of fetal SDs.Ultrasonography is still an indispensable first-line screening method. However, it has limitations, mainly in the differential diagnosis of similar SDs. In recent years, the cost of genetic testing techniques has been reduced, while their throughput has dramatically increased, which greatly benefits the precise diagnosis of fetal structural disorders [7–9]. Among all genetic testing methods, the trio WES strategy has the unique advantages in the efficiency of variation calling and the sensitivity of detectingde novo and compound heterozygous variants [9], particularly in cases with obvious structural abnormalities and normal karyotyping and CMA results.In this study, in order to investigate the genetic cause of fetal SDs in eight pregnancies with nonconsanguineous parents, a sequential detection including trio WES was performed to make a clear diagnosis. Thenin silico prediction on the functional impact of the identified novel variants was conducted.


## 2. Materials and Methods

### 2.1. Subjects

Research ethics board approval was obtained from the Human Ethics Committee of Chinese PLA General Hospital (approval no. S2018-066-01), and informed consent forms were signed by all recruited subjects. Between November 2016 and March 2018, we recruited eight families with pregnancies interrupted in their second or third trimesters due to fetal SDs based on clinical and sonographic diagnosis. Detailed information including maternal age, gestational weeks, and obstetric history was documented. Parental peripheral blood and fetal tissue or umbilical cord blood samples were obtained by routine methods during or after the procedure of odinopoeia.


### 2.2. Chromosome Karyotyping and CMA

All fetal specimens underwent conventional G-banded karyotyping test according to standard operation procedures to detect overall chromosomal anomalies. CMA tests with CytoScan 750K (Affymetrix Inc., USA) arrays were performed according to the manufacturer’s manual workflow on all fetal specimens in order to investigate genomic copy number variants with clinical significance. Data was collected and analyzed by GeneChip® Scanner 3000 with AGCC software. The pathogenicity of detected variations was determined according to guidelines issued by the American College of Medical Genetics and Genomics (ACMG) in 2011 [10].


### 2.3. Whole-Exome Sequencing

Trio WES strategy was taken to identify the causal variants. 1μg genomic DNA from 200μl peripheral blood or 5-10 mg fetal tissue was extracted using a Qiagen DNA Blood Midi/Mini Kit (Qiagen GmbH, Hilden, Germany) according to manufacturer’s protocol. DNA fragments were hybridized and captured by IDT’s xGenExome Research Panel (Integrated DNA Technologies, San Diego, USA) according to manufacturer’s protocol. The libraries were tested for enrichment by qPCR, and the size distribution and concentration were determined using an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). The Novaseq6000 platform (Illumina, San Diego, USA), along with 150 bp pair-end reads, was used for the genomic sequencing of DNA. The sequencing reads were aligned to the human reference genome (hg19/GRCh37) using the Burrows-Wheeler Aligner tool and the PCR duplicates were removed by using Picard v1.57 (http://picard.sourceforge.net/). The Verita Trekker® Variants Detection System by Berry Genomics and the third-party software GATK (https://software.broadinstitute.org/gatk/) were employed for variant calling. Variant annotation and interpretation were conducted through the use of ANNOVAR [11] and the Enliven® Variants Annotation Interpretation System authorized by Berry Genomics. During trio analysis, potential monogenetic inheritance patterns, includingde novo, autosomal recessive, autosomal dominant, X-linked recessive inheritance, mitochondrial, and, where possible, imprinted gene variations, were analyzed.In silicoanalysis using Sorting Intolerant from Tolerant (SIFT) (http://sift.bii.a-star.edu.sg/) and Polymorphism Phenotyping V2 (http://genetics.bwh.harvard.edu/pph2/) was performed in order to calculate the pathogenicity index of all novel missense variants with unknown clinical significance. The variants were classified according to the ACMG guidelines for interpretation of genetic variants [12]. For pathogenic or likely pathogenic variations reported by trio WES, Sanger sequencing or quantitative fluorescence PCR (QF-PCR) was performed as a confirmatory experiment (See Supplementary Material 2 for detailed molecular data including primer sequences, reaction systems and amplification conditions). Homological analysis among species was performed using NCBI blast online software (https://blast.ncbi.nlm.nih.gov/Blast.cgi). Three-dimensional structure prediction was conducted through the use of Modeller V9.21 (https://salilab.org/modeller/).


## 2.1. Subjects

Research ethics board approval was obtained from the Human Ethics Committee of Chinese PLA General Hospital (approval no. S2018-066-01), and informed consent forms were signed by all recruited subjects. Between November 2016 and March 2018, we recruited eight families with pregnancies interrupted in their second or third trimesters due to fetal SDs based on clinical and sonographic diagnosis. Detailed information including maternal age, gestational weeks, and obstetric history was documented. Parental peripheral blood and fetal tissue or umbilical cord blood samples were obtained by routine methods during or after the procedure of odinopoeia.


## 2.2. Chromosome Karyotyping and CMA

All fetal specimens underwent conventional G-banded karyotyping test according to standard operation procedures to detect overall chromosomal anomalies. CMA tests with CytoScan 750K (Affymetrix Inc., USA) arrays were performed according to the manufacturer’s manual workflow on all fetal specimens in order to investigate genomic copy number variants with clinical significance. Data was collected and analyzed by GeneChip® Scanner 3000 with AGCC software. The pathogenicity of detected variations was determined according to guidelines issued by the American College of Medical Genetics and Genomics (ACMG) in 2011 [10].


## 2.3. Whole-Exome Sequencing

Trio WES strategy was taken to identify the causal variants. 1μg genomic DNA from 200μl peripheral blood or 5-10 mg fetal tissue was extracted using a Qiagen DNA Blood Midi/Mini Kit (Qiagen GmbH, Hilden, Germany) according to manufacturer’s protocol. DNA fragments were hybridized and captured by IDT’s xGenExome Research Panel (Integrated DNA Technologies, San Diego, USA) according to manufacturer’s protocol. The libraries were tested for enrichment by qPCR, and the size distribution and concentration were determined using an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). The Novaseq6000 platform (Illumina, San Diego, USA), along with 150 bp pair-end reads, was used for the genomic sequencing of DNA. The sequencing reads were aligned to the human reference genome (hg19/GRCh37) using the Burrows-Wheeler Aligner tool and the PCR duplicates were removed by using Picard v1.57 (http://picard.sourceforge.net/). The Verita Trekker® Variants Detection System by Berry Genomics and the third-party software GATK (https://software.broadinstitute.org/gatk/) were employed for variant calling. Variant annotation and interpretation were conducted through the use of ANNOVAR [11] and the Enliven® Variants Annotation Interpretation System authorized by Berry Genomics. During trio analysis, potential monogenetic inheritance patterns, includingde novo, autosomal recessive, autosomal dominant, X-linked recessive inheritance, mitochondrial, and, where possible, imprinted gene variations, were analyzed.In silicoanalysis using Sorting Intolerant from Tolerant (SIFT) (http://sift.bii.a-star.edu.sg/) and Polymorphism Phenotyping V2 (http://genetics.bwh.harvard.edu/pph2/) was performed in order to calculate the pathogenicity index of all novel missense variants with unknown clinical significance. The variants were classified according to the ACMG guidelines for interpretation of genetic variants [12]. For pathogenic or likely pathogenic variations reported by trio WES, Sanger sequencing or quantitative fluorescence PCR (QF-PCR) was performed as a confirmatory experiment (See Supplementary Material 2 for detailed molecular data including primer sequences, reaction systems and amplification conditions). Homological analysis among species was performed using NCBI blast online software (https://blast.ncbi.nlm.nih.gov/Blast.cgi). Three-dimensional structure prediction was conducted through the use of Modeller V9.21 (https://salilab.org/modeller/).


## 3. Results

### 3.1. Clinical Features

In the eight families we recruited, the average age of gravidae was 31 (ranging from 24 to 38), and the average gestational week of these pregnancies was 20.9 (ranging from 16 to 29). None of the couples was consanguineous, and all couples claimed to have no family history of genetic disorders. Major clinical manifestations and information of these pregnancies were listed in Table1 (See detailed clinical data of all eight families in Supplementary Material 1).Table 1
Clinical information of recruited SD pregnancies.
Subject Maternal age (years) Gestational weeks Clinical information※ Fetal sample origin Family 1 31 16 G5P2A1: 2 daughters with normal phenotype; 2 singleton pregnancies with limb shortening Umbilical cord tissue  Family 2 31 22 G1P0: limb shortening Umbilical cord blood  Family 3 29 16 G3P0A2: 2 miscarriages; 1 SD pregnancy diagnosed at 13 gestational weeks after IVF-ET Umbilical cord tissue  Family 4 31 22 G2P0: 1 induced abortion; 1 pregnancy with thick NT(0.30cm at 12 weeks), micrognathia, ulnar and osteogenic dysplasia, abnormal hand shape Umbilical cord blood  Family 5 38 16 G2P1: 1 daughter with normal phenotype; 1 pregnancy with limb shortening Umbilical cord tissue  Family 6 28 23 G1P0: osteogenic dysplasia at 22 gestational weeks Umbilical cord tissue  Family 7 24 23 G3P0A1: 1 miscarriage; 2 pregnancies with thick NT and femur shortening Umbilical cord tissue  Family 8 36 29 G3P1; 1 daughter with normal phenotype; 1 pregnancy with osteogenic dysplasia and angled left femur at 23 gestational weeks Umbilical cord tissue ※G: gravida; P: para; A: abortus; IVF-ET: In Vitro Fertilization-Embryo Transfer.


### 3.2. Genetic Analysis

Results of karyotyping and CMA for all fetal specimens from eight pregnancies were normal. Variations with clinical significance detected by trio WES were listed in Table2, and results of corresponding Sanger sequencing and QF-PCR were demonstrated in Figure 1.Table 2
Information of detected variations.
Subject Gene Variation Disorder Inheritance pattern※ Variation effect (origin) Global MAF※ (dbSNP) SIFT score PolyPhen2 score Variation attribute (ACMG evidence levels) Family 1 SLC26A2 c.292T>C(p.Trp98Arg) Achondrogenesis IB or diastrophic dysplasia AR Missense(mother) C=0.0000/0(TWINSUK); C=0.0000/1(GnomAD_exomes); C=0.0000/1(TOPMED); C=0.0000/1(ExAC); C=0.0003/1(ALSPAC) 0 1 Likely pathogenic (PM2 + PM3 + PP2 + PP3 + PP4) c.1018_1020del(p.Val340del) In-frame deletion(father) / / / Pathogenic (PS1 + PM2 + PM4 + PP4 + PP5) Family 2 FGFR3 c.742C>T(p.Arg248Cys) Thanatophoric dysplasia, type I AD Missense(de novo) / / / Pathogenic (PS1 + PS2 + PM1 + PM2) Family 3 FLNB c.601G>A(p.Ala201Thr) Atelosteogenesis, type I or III AD Missense(de novo) / 0 0.996 Likely pathogenic (PS2 + PM2 + PM5 + PP3) Family 5 FGFR3 c.1138G>A(p.Gly380Arg) Achondroplasia AD Missense variant(de novo) A=0.0000/1(TOPMED) / / Pathogenic (PS1 + PS2 + PM1 + PM2) Family 6 FLNB c.685T>C(p.Ser229Pro) Larsen syndrome AD Missense(de novo) / / / Pathogenic (PS1 + PS2 + PM2) Family 8 TMEM38B c.344C>A (p.S115X) Osteogenesis imperfecta, type XIV AR Stop gained(Mother) / / / Pathogenic (PVS1 + PM2 + PP4) loss 1 (exon:3-4) Exon loss(Father) / / / Pathogenic (PVS1 + PM2 + PM4 + PP4) ※AD: autosomal dominant; AR: autosomal recessive; MAF: minor allele frequency.Figure 1
Results of Sanger sequencing and QF-PCR: (a) a single-base substitution inSLC26A2(c.292T>C); (b) a three-base deletion inSLC26A2(c.1018_1020del); (c) a single-base substitution inFGFR3(c.742C>T); (d) a single-base substitution inFLNB(c.601G>A); (e) a single-base substitution inFGFR3(c.1138G>A); (f) a single-base substitution inFLNB(c.685T>C); (g) a single-base substitution in exon 3 ofTMEM38B(c.344C>A); (h) one copy loss of exon 3 inTMEM38B; (i) one copy loss of exon 4 inTMEM38B.
(a) (b) (c) (d) (e) (f) (g) (h) (i)In Family 1, a compound heterozygous variation inSLC26A2 comprising c.292T>C (Figure 1(a)) and c.1018_1020del (Figure 1(b)) was identified. The two variants were inherited from the mother (c.292T>C) and father (c.1018_1020del) of the proband fetus. Sanger sequencing revealed that the two normal daughters were one carrier of c.1018_1020del as the father and the other one as wild type.Fourde novo variations were identified includingFGFR3:c.742C>T (Figure 1(c)) in Family 2,FLNB: c.601G>A (Figure 1(d)) in Family 3,FGFR3: c.1138G>A (Figure 1(e)) in Family 5, andFLNB: c.685T>C (Figure 1(f)) in Family 6.In Family 8, the compound heterozygous variations detected inTMEM38B, c.344C>A (Figure 1(g)), and loss 1 (exon: 3-4) (Figures 1(h) and 1(i) and Supplementary Figure 2-1), were inherited from the mother and father, respectively. The normal daughter was a carrier of c.344C>A, like her mother.Four novel variations were identified in this study, namely,FLNB: c.601G>A,SLC26A2: c.292T>C,TMEM38B: c.344C>A, andTMEM38B: loss 1 (exon: 3-4).In silico pathogenicity prediction was conducted on two novel missense variants (FLNB: c.601G>A andSLC26A2: c.292T>C), and the results from SIFT and PolyPhen V2 indicated them as “deleterious/probably damaging” (Table 2).The results of NCBI blast showed thatFLNB-Ala201 andSLC26A2-Trp98 amino acids were highly conserved among species (Figure 2(a), also see Supplementary Figure2-2). Additionally, the three-dimensional structure-prediction result of the mutant protein showed thatFLNB: c.601G>A (p.Ala201Thr) may cause the formation of two extra hydrogen bonds (T201-L181 and K152-S177; see Figure 2(e)).Figure 2
Biophysical analyses of theFLNB: c.601G>A(p.Ala201Thr) variation: (a) the conservation ofFLNB Ala201 between multiple species; (b) part of the three-dimensional structure of the wild typeFLNB, with the CH2 subdomain in blue and the actin-binding site 3 in yellow; (c) part of the three-dimensional structure of the mutantFLNB, with the CH2 subdomain in blue and the actin-binding site 3 in yellow; (d) the detailed red block in (b) showing Ala201; (e) the detailed red block in (c) showing Thr201 and the extra hydrogen bonds among T201 and L181, K152, and S177.
(a) (b) (c) (d) (e)


## 3.1. Clinical Features

In the eight families we recruited, the average age of gravidae was 31 (ranging from 24 to 38), and the average gestational week of these pregnancies was 20.9 (ranging from 16 to 29). None of the couples was consanguineous, and all couples claimed to have no family history of genetic disorders. Major clinical manifestations and information of these pregnancies were listed in Table1 (See detailed clinical data of all eight families in Supplementary Material 1).Table 1
Clinical information of recruited SD pregnancies.
Subject Maternal age (years) Gestational weeks Clinical information※ Fetal sample origin Family 1 31 16 G5P2A1: 2 daughters with normal phenotype; 2 singleton pregnancies with limb shortening Umbilical cord tissue  Family 2 31 22 G1P0: limb shortening Umbilical cord blood  Family 3 29 16 G3P0A2: 2 miscarriages; 1 SD pregnancy diagnosed at 13 gestational weeks after IVF-ET Umbilical cord tissue  Family 4 31 22 G2P0: 1 induced abortion; 1 pregnancy with thick NT(0.30cm at 12 weeks), micrognathia, ulnar and osteogenic dysplasia, abnormal hand shape Umbilical cord blood  Family 5 38 16 G2P1: 1 daughter with normal phenotype; 1 pregnancy with limb shortening Umbilical cord tissue  Family 6 28 23 G1P0: osteogenic dysplasia at 22 gestational weeks Umbilical cord tissue  Family 7 24 23 G3P0A1: 1 miscarriage; 2 pregnancies with thick NT and femur shortening Umbilical cord tissue  Family 8 36 29 G3P1; 1 daughter with normal phenotype; 1 pregnancy with osteogenic dysplasia and angled left femur at 23 gestational weeks Umbilical cord tissue ※G: gravida; P: para; A: abortus; IVF-ET: In Vitro Fertilization-Embryo Transfer.


## 3.2. Genetic Analysis

Results of karyotyping and CMA for all fetal specimens from eight pregnancies were normal. Variations with clinical significance detected by trio WES were listed in Table2, and results of corresponding Sanger sequencing and QF-PCR were demonstrated in Figure 1.Table 2
Information of detected variations.
Subject Gene Variation Disorder Inheritance pattern※ Variation effect (origin) Global MAF※ (dbSNP) SIFT score PolyPhen2 score Variation attribute (ACMG evidence levels) Family 1 SLC26A2 c.292T>C(p.Trp98Arg) Achondrogenesis IB or diastrophic dysplasia AR Missense(mother) C=0.0000/0(TWINSUK); C=0.0000/1(GnomAD_exomes); C=0.0000/1(TOPMED); C=0.0000/1(ExAC); C=0.0003/1(ALSPAC) 0 1 Likely pathogenic (PM2 + PM3 + PP2 + PP3 + PP4) c.1018_1020del(p.Val340del) In-frame deletion(father) / / / Pathogenic (PS1 + PM2 + PM4 + PP4 + PP5) Family 2 FGFR3 c.742C>T(p.Arg248Cys) Thanatophoric dysplasia, type I AD Missense(de novo) / / / Pathogenic (PS1 + PS2 + PM1 + PM2) Family 3 FLNB c.601G>A(p.Ala201Thr) Atelosteogenesis, type I or III AD Missense(de novo) / 0 0.996 Likely pathogenic (PS2 + PM2 + PM5 + PP3) Family 5 FGFR3 c.1138G>A(p.Gly380Arg) Achondroplasia AD Missense variant(de novo) A=0.0000/1(TOPMED) / / Pathogenic (PS1 + PS2 + PM1 + PM2) Family 6 FLNB c.685T>C(p.Ser229Pro) Larsen syndrome AD Missense(de novo) / / / Pathogenic (PS1 + PS2 + PM2) Family 8 TMEM38B c.344C>A (p.S115X) Osteogenesis imperfecta, type XIV AR Stop gained(Mother) / / / Pathogenic (PVS1 + PM2 + PP4) loss 1 (exon:3-4) Exon loss(Father) / / / Pathogenic (PVS1 + PM2 + PM4 + PP4) ※AD: autosomal dominant; AR: autosomal recessive; MAF: minor allele frequency.Figure 1
Results of Sanger sequencing and QF-PCR: (a) a single-base substitution inSLC26A2(c.292T>C); (b) a three-base deletion inSLC26A2(c.1018_1020del); (c) a single-base substitution inFGFR3(c.742C>T); (d) a single-base substitution inFLNB(c.601G>A); (e) a single-base substitution inFGFR3(c.1138G>A); (f) a single-base substitution inFLNB(c.685T>C); (g) a single-base substitution in exon 3 ofTMEM38B(c.344C>A); (h) one copy loss of exon 3 inTMEM38B; (i) one copy loss of exon 4 inTMEM38B.
(a) (b) (c) (d) (e) (f) (g) (h) (i)In Family 1, a compound heterozygous variation inSLC26A2 comprising c.292T>C (Figure 1(a)) and c.1018_1020del (Figure 1(b)) was identified. The two variants were inherited from the mother (c.292T>C) and father (c.1018_1020del) of the proband fetus. Sanger sequencing revealed that the two normal daughters were one carrier of c.1018_1020del as the father and the other one as wild type.Fourde novo variations were identified includingFGFR3:c.742C>T (Figure 1(c)) in Family 2,FLNB: c.601G>A (Figure 1(d)) in Family 3,FGFR3: c.1138G>A (Figure 1(e)) in Family 5, andFLNB: c.685T>C (Figure 1(f)) in Family 6.In Family 8, the compound heterozygous variations detected inTMEM38B, c.344C>A (Figure 1(g)), and loss 1 (exon: 3-4) (Figures 1(h) and 1(i) and Supplementary Figure 2-1), were inherited from the mother and father, respectively. The normal daughter was a carrier of c.344C>A, like her mother.Four novel variations were identified in this study, namely,FLNB: c.601G>A,SLC26A2: c.292T>C,TMEM38B: c.344C>A, andTMEM38B: loss 1 (exon: 3-4).In silico pathogenicity prediction was conducted on two novel missense variants (FLNB: c.601G>A andSLC26A2: c.292T>C), and the results from SIFT and PolyPhen V2 indicated them as “deleterious/probably damaging” (Table 2).The results of NCBI blast showed thatFLNB-Ala201 andSLC26A2-Trp98 amino acids were highly conserved among species (Figure 2(a), also see Supplementary Figure2-2). Additionally, the three-dimensional structure-prediction result of the mutant protein showed thatFLNB: c.601G>A (p.Ala201Thr) may cause the formation of two extra hydrogen bonds (T201-L181 and K152-S177; see Figure 2(e)).Figure 2
Biophysical analyses of theFLNB: c.601G>A(p.Ala201Thr) variation: (a) the conservation ofFLNB Ala201 between multiple species; (b) part of the three-dimensional structure of the wild typeFLNB, with the CH2 subdomain in blue and the actin-binding site 3 in yellow; (c) part of the three-dimensional structure of the mutantFLNB, with the CH2 subdomain in blue and the actin-binding site 3 in yellow; (d) the detailed red block in (b) showing Ala201; (e) the detailed red block in (c) showing Thr201 and the extra hydrogen bonds among T201 and L181, K152, and S177.
(a) (b) (c) (d) (e)


## 4. Discussion

The clinical heterogeneity of skeletal dysplasias is strong. Some fetal phenotypes are relatively vague and some fetal phenotypes may not have obvious manifestations until the third trimester, which leads to the challenge in ultrasonic and differential diagnosis. Meticulously designed strategy of genetic testing may help solve this problem. Several studies have discussed the advantage of trio WES with respect to the efficiency of variation screening [9, 13, 14], which makes possible the application to prenatal diagnosis.TheSLC26A2(MIM ∗606718) protein transports ions, particularly sulfate ions, across cell membranes that help cartilage to produce proteoglycans [15]. The impaired function of theSLC26A2 product would be expected to lead to undersulfation of proteoglycans in the cartilage matrix and thereby cause a spectrum of SDs, including achondrogenesis IB (ACG-IB,MIM #600972), atelosteogenesis II (MIM #256050), De la Chapelle dysplasia (MIM #256050), diastrophic dysplasia (DTD, MIM # 222600) and epiphyseal dysplasia, multiple, i.e., 4 (MIM # 226900). Superti-Furga et al. first established an association betweenSLC26A2 and ACG-IB [16]. Among the compound heterozygous mutations detected in Family 1, the variantSLC26A2: c.1018_1020del (p.Val340del) was known as pathogenic, and its homozygous mutation was first reported to cause ACG-IB by Superti-Furga et al. [17]. Contrastingly, the variantSLC26A2: c.292T>C (p.Trp98Arg) (no. Rs753193118 in the dbSNP database) has not yet been reported as pathogenic. Its frequency in the gnomAD database is 4.06×10−6 (http://gnomad.broadinstitute.org/), and it was predicted as deleterious by the SIFT and PolyPhen software programs. We then interpreted c.292T>C as likely pathogenic according to the ACMG criteria with evidence levels PM2+PM3+PP2+PP3+PP4. The relationship between different types of homozygous or complex heterozygous mutations and the severity of corresponding disease phenotypes was discussed [17]. It is believed that DTD is associated with reducedSLC26A2 expression, while ACG-IB results from the null mutations of it. Thus, the ability to predict the specific disease type of the fetus in Family 1 depends on the revelation of the expression level ofSLC26A2.FGFR3(MIM ∗134934) belongs to the fibroblast growth factor family which plays an important role in cell proliferation and differentiation, angiogenesis, wound healing, and embryo development (https://ghr.nlm.nih.gov/gene/FGFR3). It is believed that theFGFR3 protein regulates bone growth by limiting ossification progress, particularly in long bones [18]. The two variants detected in Family 2 (FGFR3: c.742C>T) and Family 5 (FGFR3: c.1138G>A) have been reported many times as pathogenic [19–22]. However, these two variants lead to different disorders (Table 2) that have different phenotypes and prognoses, which are difficult to distinguish from fetal sonographic indications.FGFR3: c.742C>T is one of the two most common mutations in thanatophoric dysplasia, type I (MIM #187600), whileFGFR3: c.1138G>A contributes to more than 90% of the condition in achondroplasia (MIM #100800) patients [23, 24]. These results indicate that genetic analysis is of great significance for prognosis prediction and clinical consultation to these families with variations in identical genes.Filamins, includingFLNB(MIM ∗603381), are actin-binding proteins that also interact with multiple receptors and intracellular proteins, which in turn regulate cytoskeleton-dependent cell proliferation, differentiation, and migration [25]. Previous studies have shown that heterozygous missense variations inFLNBlead to a spectrum of severe SDs including atelosteogenesis type I (AOI, MIM #108720), atelosteogenesis type III (AOIII, MIM #108721), Boomerang dysplasia (MIM #112310), Larsen syndrome (MIM #150250), and spondylocarpotarsal synostosis syndrome (MIM #272460). In our study, the variant in Family 3,FLNB: c.601G>A (p.Ala201Thr), has not been previously reported, but it shares the same amino acid that is affected by a variant,FLNB: c.602C>T (p.Ala201Val), as detected in a neonate with AOIII [26]. We then interpreted c.601G>A as likely pathogenic according to the ACMG criteria (PS2+PM2+PM5+PP3). Sawyer et al. pointed out that missense mutations in particular regions ofFLNB may follow the mechanism of gain of function and enhance its binding affinity with actin [27]. The structure-prediction result in our study is likely consistent with this concept. The identified variationFLNB: c.685T>C (p.Ser229Pro) in Family 6 was previously reported as a pathogenic variant responsible for Larsen syndrome [28]. These two variations in our study are located within the CH2 subdomain of the actin-binding domain (ABD) inFLNB, and they may cause the dysregulation of actin-filamin interaction, which associates to the skeletal phenotype spectrum of the probands.TMEM38B(MIM ∗611236) encodes trimeric intracellular cation-B (TRIC-B) protein, which expresses differently in various tissues and cells of animals. TRIC-B channels act as counter-ion channels that function in synchronization withCa2+ release from intracellular stores [29]. Pathogenic variations inTMEM38Bwere reported to cause a rare autosomal recessive type of osteogenesis imperfecta (OIXIV, MIM #615066). Patients of OIXIV usually develop moderately severe OI. They have various fracture frequencies, mildly to moderately short stature, and gray-to-blue sclera but no occurrence of dental defects [30]. To our knowledge, six different mutations of TMEM38B have been reported in previous studies [30, 31]. The compound heterozygous variation detected in Family 8 consists of two novel variants: one (c.344C>A) can cause premature termination of protein translation, and the other (loss 1 (exon: 3-4)) may result in truncated protein. Each of those variants has a serious impact on the function of TRIC-B protein and is classified as pathogenic according to ACMG criteria. Nevertheless, the mechanism of variable expressivity in different OIXIV cases remains to be studied. Moreover, the detection of these two variants expanded the mutant spectrum of OIXIV and will be very helpful in the continued investigation ofTMEM38B function.In the remaining two families (4 and 7), no variation with clear clinical significance was detected. Thus, further research is essential. Pathogenic variations may be identified after reanalysis over an extensive period of time, particularly given the emergence of new disease-causing genes, and new mechanisms or pathogenesis may be discovered through in-depth investigation.This study helped SD families to identify the cause and accurately assessed the risks inherent with further pregnancies. Each participating family with positive results had a different inheritance pattern of disease and therefore a different risk of recurrence during pregnancy. The families corresponding to autosomal recessive pattern (Families 1 and 8) have a 25% risk in each pregnancy. The families corresponding to autosomal dominant pattern (Families 2, 3, 5, and 6) commonly have minimal risk of recurrence, but it will be relatively higher if there is germinal mosaicism [32, 33]. This issue should be considered in the clinical consultation and subsequent pregnancy examination.


## 5. Conclusions

Specific laboratory diagnosis is difficult with respect to cases that involve skeletal dysplasias, given the low incidence of SDs as well as their strong clinical and genetic heterogeneity, particularly in the field of prenatal diagnosis. Comprehensive application of multiple genetic techniques can effectively improve the diagnosis rate of SDs. Thus, the trio WES strategy provides a robust methodological supplement in case there is lack of clear imageological evidence and sufficient clinical experience.



---

*Source: 2492590-2019-05-14.xml*