Development and validation of an expanded targeted sequencing panel for non-invasive prenatal diagnosis of sporadic skeletal dysplasia
BMC Medical Genomics volume 14, Article number: 212 (2021)
Skeletal dysplasia (SD) is one of the most common inherited neonatal disorders worldwide, where the recurrent pathogenic mutations in the FGFR2, FGFR3, COL1A1, COL1A2 and COL2A1 genes are frequently reported in both non-lethal and lethal SD. The traditional prenatal diagnosis of SD using ultrasonography suffers from lower accuracy and performed at latter gestational stage. Therefore, it remains in desperate need of precise and accurate prenatal diagnosis of SD in early pregnancy. With the advancements of next-generation sequencing (NGS) technology and bioinformatics analysis, it is feasible to develop a NGS-based assay to detect genetic defects in association with SD in the early pregnancy.
An ampliseq-based targeted sequencing panel was designed to cover 87 recurrent hotspots reported in 11 common dominant SD and run on both Ion Proton and NextSeq550 instruments. Thirty-six cell-free and 23 genomic DNAs were used for assay developed. Spike-in DNA prepared from standard sample harboring known mutation and normal sample were also employed to validate the established SD workflow. Overall performances of coverage, uniformity, and on-target rate, and the detecting limitations on percentage of fetal fraction and read depth were evaluated.
The established targeted-seq workflow enables a single-tube multiplex PCR for library construction and shows high amplification efficiency and robust reproducibility on both Ion Proton and NextSeq550 platforms. The workflow reaches 100% coverage and both uniformity and on-target rate are > 96%, indicating a high quality assay. Using spike-in DNA with different percentage of known FGFR3 mutation (c.1138 G > A), the targeted-seq workflow demonstrated the ability to detect low-frequency variant of 2.5% accurately. Finally, we obtained 100% sensitivity and 100% specificity in detecting target mutations using established SD panel.
An expanded panel for rapid and cost-effective genetic detection of SD has been developed. The established targeted-seq workflow shows high accuracy to detect both germline and low-frequency variants. In addition, the workflow is flexible to be conducted in the majority of the NGS instruments and ready for routine clinical application. Taken together, we believe the established panel provides a promising diagnostic or therapeutic strategy for prenatal genetic testing of SD in routine clinical practice.
Skeletal dysplasia (SD) is a heterogeneous group of genetic disorders associated with various abnormalities of bones and joints. There are 461 different conditions classified into 42 groups primarily based on the clinical, radiographic, and molecular phenotypes . The clinical manifestations vary in severity from mild phenotypes to severe abnormalities with perinatal mortality due to lung hypoplasia and respiratory complications . Previous studies showed that the incidences of non-lethal and lethal SD are around 1 per 5,000 births [3, 4] and 0.95–1.5 per 10,000 births [5,6,7], respectively. The three most common SD types are thanatophoric dysplasia (around 29%) , osteogenesis imperfecta type 2 (14%), and achondrogenesis (9%) , which account for 40 to 60 percent of all lethal SD [6, 10,11,12]. Importantly, among the perinatal deaths of SD, 23–32% occur during the first week of life [3, 6]. The unexpected loss and highly distressing event result in a hugely psychological burden on bereaved parents. Therefore, it remains in desperate need of precise and accurate prenatal diagnosis of SD in early pregnancy.
Traditionally, the prenatal diagnosis of SD relies on ultrasonography, followed by confirmation using either invasive monogenic testing or post-delivery radiographs and autopsy. Although the ultrasound assessment of SD includes a number of features and criteria, the accuracy of routine ultrasound approach was reported as 40–60%, because the high heterogeneity in genetic defects and large phenotypic variability [13, 14]. The diagnostic confirmation of lethality was still difficult due to the lack of systematic approach . In addition, the traditional genetic testing requires invasive procedures to obtain the fetal specimen such as amniocentesis and chorionic villus sampling, and the Sanger sequencing for each individual gene is costly and time-consuming . Therefore, the earlier prenatal genetic diagnosis is critical issue in maternal–fetal precision medicine.
The discovery of circulating fetal cell-free DNA (cfDNA) has opened a promising opportunity for early detection of fetal genetic defects [16, 17]. Given that the fetal fraction in maternal blood is estimated to be around 4%-30% depending on the gestational age and maternal weight effects [18,19,20], it requires a more sensitive platform to detect the low-frequency variants in association with the fetus using maternal plasma cfDNA. Owing to the advantages of high-throughput, parallel sequencing techniques, the non-invasive prenatal testing (NIPT) has been developed and applied to identify fetal chromosomal aneuploidies using the cfDNA extracted from a blood of pregnant woman [21, 22]. Although the development is still in its infancy, the precise and flexible applications of NIPT in the detection of single-nucleotide variant (SNV) and small deletion/insertion (Indel) have been reported [23,24,25].
To detect low-frequency variants at single-base resolution, targeted sequencing (targeted-seq) that offers ultra-deep coverages on the genomic region of interest has been shown to function accurately and adequately [26, 27]. Furthermore, targeted-seq is cost-effectiveness in comparison with whole genome sequencing or whole-exome sequencing . Two major methods, amplicon- and capture-based that utilize multiplex polymerase-chain reaction (PCR) and hybridization by probes, respectively, have been established to enrich targeted regions of interest for downstream sequencing. Previously, an amplicon-based targeted-seq for prenatally screen of FGFR3 gene [23, 24] and a capture-based targeted-seq for screening 497 genes including several SD-associated genes  have been developed. Due to the restrictions by fewer mutation hotspots or lower coverage depth, it limits the application of these panels in routine clinical practice. Thus it remains in need to establish an expanded targeted sequencing panel for the detection of sporadic SD in early pregnancy.
The advantages of amplicon-based technique are the less requirement of input DNA and higher on-target rate when compared to the capture-based method, thus it is more attractive to be applied in NIPT where the fetal cfDNA is in limited amount and high accuracy is on-demand. Here, we report the development of an expanded amplicon-based targeted-seq panel for the detection of SD at ultra-depth level. This panel examines a total of 87 recurrent hotspots located in 5 most common genes (FGFR2, FGFR3, COL1A1, COL1A2, COL2A1) associated with 11 dominant SD, in which 5 mutations responsible for lethal SD are included. The targeted-seq workflow can be applied to genomic DNA as well as cfDNA, both achieve high sensitivity and specificity of 100%. Most importantly, the workflow is optimized to be carried out on both Illumina and Ion Torrent platforms, which account for more than 90% of instrument marketplace. Collectively, we streamline the targeted-seq workflow of two different systems and confirm its performance in clinical practice. The workflow provides timely detection of SD at early pregnancy and shows the true value of prenatal precision medicine.
A total of 59 DNA specimen, including 36 cell-free (cf) and 23 genomic (g) DNA, were collected from 6 families and 31 pregnant women (Table 1). Among them, family 1 and 2 are with fetus that were diagnosed with suspected dwarfism and osteogenesis imperfecta by ultrasound examination, respectively, while the parents were all asymptomatic. All the others were cases collected during regular pregnancy examination to screen genetic defect at early gestational age and have been followed through the gestation period until baby was delivered. Written and signed informed consent was obtained from all subjects of this study and approved by the Institutional Review Board (IRB) of the National Cheng Kung University.
Genomic DNA collection and processing
Roughly 10 ml of peripheral blood was drawn from all participants and in two cases of Family l and Family 2, umbilical tissue was obtained from the deceased fetus. The buffy coat was separated from peripheral blood through centrifugation at 3000 rpm for 10 min at room temperature. DNA of umbilical cord and buffy coat were isolated by QIAamp DNA Mini Kit (QIAGEN, Hilden, Germany) according to the manufacturer’s protocol. The genomic DNA was stored at -20 °C for long term storage.
The standard DNA “NA11316” harboring homozygous FGFR3 mutation (c.1138 G > A) was purchased from Coriell Institute (Coriell Institute, Camden, New Jersey, USA). To generate the spike-in DNA, NA11316 DNA was mixed with normal maternal gDNA (from Family 2) in different percentage (2.5%, 5%, 10%) to mimic different fetal DNA fractions in the maternal plasma.
Plasma preparation and cfDNA extraction
Maternal plasma was obtained by a two-step centrifugation process. In the first step, the maternal peripheral blood was separated by centrifugation of the collection tube at 2,000 g for 15 min at room temperature and the supernatant was transferred to a 1.5 ml centrifuge tube. The sample was then centrifuged at 14,000g for 10 min at room temperature and the clear plasma was stored at − 80 °C until further processing.
To isolate cfDNA, maternal plasma was centrifuged at 16,000g for 5 min at room temperature, followed by cfDNA extraction using QIAamp Circulating Nucleic Acid Kit (QIAGEN) according to the manufacturer’s protocol. The quality and quantity of cfDNA were determined by LabChip GX Nucleic Acid Analyzer (PerkinElmer, Waltham, Massachusetts, USA) and Qubit Fluorometer (Thermo Fisher Scientific, Waltham, Massachusetts, USA), respectively. The cfDNA was stored at − 20 °C until further process.
Design primer sets for amplicon-based targeted sequencing
To develop an expanded panel to facilitate the detection of SD, we focus on five most relevant genes (FGFR2, FGFR3, COL1A1, COL1A2 and COL2A1) which cover 11 common clinical conditions with dominant inheritance [3, 8, 30,31,32,33]. The clinical conditions include thanatophoric dysplasia, type I-II (TD1-2; OMIM #187600 and # 187601), achondroplasia (ACH; OMIM #100800), osteogenesis imperfecta, type I- IV (OI1-4; OMIM #166200, #166210, #259420, #166220), achondrogenesis, type II (ACG2; OMIM #200610), Apert syndrome (OMIM #101200), Crouzon syndrome (OMIM #123500), and Pfeiffer syndrome (OMIM #101600).
We selected 87 recurrent mutation hotspots that comprises 82 SNVs, 2 multinucleotide variants (MNVs) and 3 indels (Additional file 1: Table S1) based on the information provided from the OMIM database (https://www.omim.org), GeneReviews (https://www.ncbi.nlm.nih.gov/books/NBK1116/), and ClinVar database (https://www.ncbi.nlm.nih.gov/clinvar/). The gene panel covers literature-supported recurrent mutations for 11 common SD clinical conditions at various frequencies: TD1 (~ 3.5%) , TD2 (> 99%) , ACH (> 99%) [34, 36, 37], OI (> 68%) [38,39,40], ACG2 (> 28%) [34, 41], Apert syndrome (~ 99%) [32, 42, 43], Crouzon syndrome (> 39%) [32, 44, 45], Crouzon syndrome w/acanthosis nigricans (100%) , and Pfeiffer syndrome (> 63%) [32, 46, 47]. We applied Ion Ampliseq Designer (version 6.1.3, Thermo Fisher Scientific) to design primers and outputted a set of 51 amplicons in one single-tube multiplexing specification (Additional file 1: Table S1).
Library construction and amplicon sequencing
The amplicon libraries were constructed by Ion AmpliSeq™ Library Kit 2.0 (Thermo Fisher Scientific) according to the manufacturer’s protocol, and sequenced on Ion Proton (Thermo Fisher Scientific) or NextSeq550 instruments (Illumina, San Diego, California, USA). Briefly, approximately 10 ng of gDNA, or 1 ng of cfDNA were amplified using the 51 pairs of pooled primers in a single-tube multiplex PCR setting. The amplified libraries were then subjected to partial digestion, barcode ligation, purification, and run on Ion Proton sequencer (Fig. 1A). To modify the experimental procedure for the NextSeq550 platform, QIAseq Adapter (QIAGEN) was used for the adaptor ligation and followed by one-cycle PCR to join the adaptor and amplicon. Library concentration was determined by Ion Library TaqMan™ Quantitation Kit (Thermo Fisher Scientific) and KAPA Library Quantification Kit Illumina® Platforms (KAPA Biosystems, Wilmington, MA, USA) for Ion Proton and NextSeq550, respectively. As different barcode and index systems were used, the sizes of amplicon library were around 200 bp and 250 bp for Ion Proton and NextSeq550 systems, respectively.
After the sequencing, the low-quality bases (Q score < 20) were trimmed and the filtered reads were aligned to human reference genome (GRCh38) using the Torrent Mapping Alignment Program (TMAP, Thermo Fisher Scientific) and Burrows-Wheeler Aligner (BWA)-v0.7.17  for Ion Proton and NextSeq550 platforms, respectively (Fig. 1B). The output BAM file was used to call variant by using VariantCaller (v22.214.171.124, Thermo Fisher Scientific) under default parameters. To identify low-frequency variants in association with fetal cfDNA, a widely used haplotype-based variant calling tool for the detection of somatic or low-frequency variant , GATK-Mutect2, was applied. Coverage, uniformity, and on-target rate were calculated by using CoverageAnalysis (v126.96.36.199,, Thermo Fisher Scientific) or Samtools-v1.6 . Alternative allele frequency of 2% and 20% were set as the thresholds for calling somatic and germline variants, respectively.
Statistical analyses, including two-way analysis of variance (ANOVA), Pearson correlation coefficient (r), Kruskal–Wallis test, and two-tailed t-test, were conducted by using GraphPad Prism version 7.05 for Windows (GraphPad Software, San Diego, California USA). Data were considered significant when the P values were < 0.05.
Establishment of amplicon-based targeted-seq workflow
To establish amplicon-based targeted-seq workflow (Fig. 1A, B, left), 10 ng of gDNA from 17 individuals were used for library construction and sequencing on Ion Proton platform. After sequencing, the fastq data was trimmed and then mapped to the human genome. We obtained 100% coverage of each amplicon at averaged sequencing depth of 13,160x. The uniformity and on-target rate were 99.23 ± 0.44% and 96.39 ± 3.76%, respectively. Although the mean amplicon read ratio (define by read counts in each amplicon divided to total read counts per sample) of 1.96 ± 0.45% (ratio ranged from 0.63% to 2.96%) is closed to the expected ratio of 1.96%, the amplicon read distributions are significantly different across 51 amplicons due to the characteristic of primers and GC content of each amplicon (P < 0.0001, Fig. 1C). Nevertheless, the pattern of amplicon distribution was no difference among 17 samples (P > 0.9999, Fig. 1C) and the results demonstrated high reproducibility of established workflow.
Next, we aimed to generalize the workflow to be used on NextSeq550 platform (Fig. 1A, B, right). Using slightly modified procedure and under a base coverage of 3,743x, we obtained the uniformity and on-target rate were 99.14 ± 0.92%, and 97.01 ± 0.71%, respectively. The performance is comparable to the data obtained using Ion Proton. Albeit the variation of reads distribution was slightly higher using NextSeq550 (mean ratio of 1.96 ± 1.14%, ranged from 0.49% to 5.39%) than that was using Ion Proton, common features of significant differences of read distribution across 51 amplicons (P < 0.0001, Fig. 1D) and high reproducibility among samples (P > 0.9999, Fig. 1D) were observed. Collectively, these data demonstrated the high amplification efficiency and quality in the established targeted-seq workflows. Furthermore, the flexibility of running the established targeted-seq workflow on both Ion Torrent system (i.e., Ion Proton and Ion PGM) and Illumina system (i.e., HiSeq, MiSeq, and NextSeq), the instrument implemented in the majority of diagnostic laboratories, ensure the flexibility for clinical usage.
Validate the amplicon-based targeted-seq pipeline for the detection of skeletal dysplasia
We then seek to validate the targeted-seq pipeline using DNAs from a family with a fetus diagnosed as dwarfism by ultrasound images at gestational age of 21 weeks (Family 1, Fig. 2A). The targeted-seq workflow was applied to the parents’ and fetal genomic DNAs. Deep sequencing results clearly showed that the abortus carried a heterozygous mutation in FGFR3 (c.1949 A > T, p.Lys650Met), while no mutation was detected from both parents (Fig. 2B). The identified de novo mutation was subsequently confirmed using Sanger sequencing (Fig. 2C). The mutation in FGFR3 (c.1949 A > T) is known to cause thanatophoric dysplasia, a severe short-limb dwarfism syndrome that is usually lethal in the perinatal period (OMIM #187601), which supports the findings of targeted-seq. Taken together, the results revealed that the established targeted-seq workflow effectively identify germline mutation in association with SD.
We next challenged the capability of identifying low-frequency variant using established targeted-seq workflow. We utilized a normal maternal genomic DNA and mixed with a standard DNA sample harboring FGFR3 single base mutation (NM_000142: c.1138G > A, p.G380R) at different percentages of 2.5%, 5%, and 10%. The spike-in DNA was used to mimic the condition in which the fetal cfDNA is mixed with maternal cfDNA at various fractions. Since the read depth has impact on the sensitivity of low-frequency variant detection, amplicon libraries were prepared and sequenced at expected coverages of 10,000x, 15,000x, and 25,000x on Ion Proton platform. We obtained the mean base coverages of 12,580x, 16,059x, and 23,904x, that were not too far away from our expectations. In spite of different read depths, the uniformities of 99.32 ± 0.21%, 98.31 ± 0.32% and 99.61 ± 0.29%, and on-target rates of 98.41 ± 0.20%, 98.43 ± 0.16% and 98.43 ± 0.20% were all similar and represented good quality results. In addition, amplicon read distributions from 12,580x (ranges from 0.71% to 2.96%, mean ratio = 1.96 ± 0.55%), 16,059x (ranges from 0.66% to 2.91%, mean ratio = 1.96 ± 0.52%), and 23,904x samples (ranges from 0.72% to 2.86%, mean ratio = 1.96 ± 0.55%) showed consistency regardless the different depths among samples (P > 0.99, Fig. 2D). The total depth at the spike-in FGFR3 mutation site ranges from 5,645 to 16,835. Most importantly, we identified the spike-in FGFR3 mutation (c.1138G > A) in all samples, even in the lowest 2.5% one with the minimum depth at 12,580x (Table 2). In addition, the FGFR3 mutation (c.1138G > A) was detected in all spike-in samples at the values approximated to the modeled spike-in fractions (P = 0.11, Table 2). These data demonstrated that the workflow enables the detection of low-frequency variant accurately.
We duplicated the targeted-seq procedure with various spike-in DNA samples on NextSeq550. At mean base depth of 6,821x, we achieved uniformity and on-target rate of the samples were 98.67 ± 1.37%, and 98.03 ± 0.77%, respectively. Notably, similar amplicon read distributions among three different spike-in fractions were observed, suggesting cfDNA fraction (as indicated by various spike-in DNA fractions) has no effect on the dynamic of amplicon distribution (P = 0.97, Fig. 2E). As there are additional 55 single nucleotide polymorphisms (SNPs) recorded in the regions covered by the 51 amplicons (Additional file 2: Table S2) in Taiwanese population , we compared the results of variant calling using NextSeq550and Ion Proton systems. There were 5 SNPs identified from 3 spike-in DNA samples sequenced on two platforms. Using calculated minor-allele frequency (MAF) for these SNPs, Fig. 2F showed highly concordant results from 2 platforms (r = 0.79, P < 0.001). Taken together, these data demonstrated the accuracy and high confidence of established targeted-seq workflows for the detection of low-frequency variants, suggesting the possible application of workflow in detecting genetic defects of SD in early pregnancy.
The sensitivity and specificity of established targeted-seq workflow in clinical usage
Given that we have validated the targeted-seq workflow, we sought to evaluate its sensitivity and specificity using cfDNA from early pregnancy stage. Library construction and targeted-seq were performed using 36 cfDNA samples (Table 1). Among them, one sample (Family 2) was suspected to have SD at regular pregnancy examination. The fetus presented short long bones and multiple fractures using ultrasound image at 21 weeks of gestational age, which resembled a case of severe osteogenesis imperfecta (OI)-related symptoms. No mutation was identified in the cfDNA or in gDNAs from fetal umbilical cord tissue. To investigate the underlying genetic defect of the proband, we performed whole-exome sequencing and found novel compound heterozygous mutations in CRTAP gene consisting of SNV and deletion, suggesting the case was a rare autosomal recessive form of OI . The findings supported the true negative result of cfDNA using the established targeted-seq assay. For the rest of 35 cfDNA samples, we found no mutation in the 87 hotspots from the tested SD panel. Concordant to the testing results, follow-up from these women throughout the gestation to the baby delivered showed no indication of fetal abnormality, demonstrating the accurate call for true negative results. Overall, there are 9 true positive detections (3 spike-in samples in 3 different read depths) and 36 true negative detections from our low-frequency variant tests. Although the sample size is still small, our SD targeted-seq assay reaches 100% sensitivity and 100% specificity thus it suggests the established workflow is ready for routine clinical application.
The diagnostics laboratories constantly strive to gain a precise and extensive genetic characterization of patients at increased efficiency, robust and cost-effectiveness. Herein, we developed an amplicon-based targeted-seq panel for the detection of 87 pathogenic mutations, including 82 SNVs, 2 MNVs and 3 indels distributing on 5 most common genes that lead to 11 autosomal dominant SDs. The targeted-seq pipeline can be applied to cell-free and genomic DNA, both of which achieve high sensitivity and specificity of 100%. Given that the turnaround time from blood sample collection to issue test report is ~ 4 days, it is feasible and cost-effective for routine screening for pregnant women at the early pregnancy. As such expanded panel, it can rule out the most common SDs from fetus and reduce the anxiety and stress in pregnancy. For the fetus presenting abnormal skeletal growth, this assay also can be a great help in genetic confirmation using genomic fetal DNA from amniocentesis or chorionic villus sampling.
The accuracy and precision are key determinants of the detection assay used in clinical diagnosis. Several studies have provided different panels to detect monogenic disorders in NIPT recently [23, 24, 29, 53]. For instance, two amplicon-based targeted-seq panels consisting of 18 and 22 hotspots in FGFR3 genes were proposed to screen TD and ACH [23, 24], respectively. In addition, Malcher et al. also reported a capture-based targeted-seq procedure for screening a large panel of 497 genes in NIPT with the mean and median coverage across all sample are 267x and 222x . Since the presence of fetal fraction in the maternal cfDNA can be as small as 4–5%, it needs higher coverage to precisely detect the low-frequency mutation in association with fetal cfDNA. Under the circumstance, lower overall coverage or read depth can easily lead to false negative results. Moreover, Russo et al. developed a amplicon-based targeted-seq for 337 mutation hotspots associated with autosomal recessive and dominate disorders for NIPT, including 42 hotspots in FGFR2 and FGFR3 genes . However, it was a pilot study without reporting statistics of the outcome in terms of accuracy, specificity or sensitivity . In the current study, we established the targeted-seq workflow that shows great performance in deep coverage, high on-target rate and uniformity. Most importantly, the workflow is able to detect variant with as low as 2.5% minor allele frequency. Because the panel is designed to detect autosomal dominant SD bearing a pathogenic mutation in heterozygous state, the minimal requirement of fetal cfDNA should be 5%, which has been set as minimal receiving criteria of fetal fraction in maternal plasma. In this study, all the clinical cfDNA are within the reporting range of fetal fraction (ranging from 5.3% to 30.9%, mean = 16.7 ± 6.7%), suggesting the reliable and accurate interpretation of the genetic testing. According to our assay development, we will run deep targeted-seq with 1 M total reads to reach at least 2000x base depths for all of 87 hotspots. Thus it provides most accurate and robust panel so far in identifying mutations for common dominant SD.
To date, several approaches have been developed for non-invasive prenatal diagnosis of a range of inherited and/or de novo transmission disorders using non-NGS or NGS-based techniques, such as droplet-digital PCR, PCR with restriction enzyme digest (PCR-RED), real-time quantitative PCR, relative dosage haplotype dosage analysis (RHDO) as well as targeted-seq. Digital PCR is useful for precise fetal genotyping by analyzing relative mutation dosages in NIPT and shows capability to detect a monogenic disorder independently of parental origin . However, it requires higher amount of input DNA to reach high sensitivity and has technical difficulty of performing high-throughput in multiplex PCR, which means only for limited number of disorder. In addition, the PCR-RED is implemented in the diagnosis of FGFR3-associated SD in NIPT and shows high accuracy in follow-up of pregnancy outcome . Unfortunately, the PCR-RED relies on subjective analysis using agarose gel electrophoresis and has an inconclusive rate of around 8% , suggesting the impracticable for the clinical practice. The NGS-based RHDO is applied for NIPT for monogenic disorders based on analysis of SNPs in haplotype blocks, and is able to detect complex mutations in gene loci, such as CYP21A2 , and DMD . The limitation of RHDO is the requirement of haplotype construction using parental samples, and therefore the cost is higher and prohibitive to clinical application . Collectively, the above-mentioned limitations from various techniques prohibit them to be applied in NIPT for monogenic disorder detection. Intriguingly, it is reported that advanced paternal age is associated with risk of certain SDs, including TD, OI, Apert syndrome, Crouzon syndrome, and Pfeiffer syndrome [59,60,61,62,63,64], suggesting a hidden risk for transmission of deleterious variants to the offspring when aging. The paternal age-associated mutation hotspots, including FGFR2: c.755C > G, c.755C > T and c.758C > G , and FGFR3: c.1138G > A, c.1138G > C , c.1454A > G, c1948A > G and c.1949A > T [60, 62], are included in the expanded gene panel. Therefore, it further puts emphasis on the routine prenatal screen even when the parents are asymptomatic and the proposed targeted-seq workflow is useful to seize the mutations. Based on our findings, the expanded gene panel for a range of 11 SDs showed true negative results for those pregnant women with normal ultrasonography, suggesting the application is ready for screening for SD in earlier pregnancy. The turnaround time of 4 days and minimum requirement of 1 ng of cfDNA make the established targeted-seq panel an excellent choice to be offered in routine clinical practice.
It is of notice that the mutations result in rare autosomal recessive form of severe OI-related symptoms in family 2  were not detected in the established panel. It is still challenging to detect mutations for autosomal recessive diseases in NIPT as the mutations are contributed from both parents thus it may present at various fractions in the mutation sites. Traditionally it requires invasive approach to obtain the fetal specimen to distinguish the origin of maternal mutation from the fetal cfDNA carries homozygous or compound heterozygous mutations . Recently, Luo et al. reported a capture-based targeted-seq workflow that can screen for chromosomal aneuploidy of chromosomes 13, 18, and 21, microdeletions, and autosomal recessive disorders simultaneously . It shows promising positive prediction rate of 100% for chromosomal aneuploidy and copy number variations, although the value for autosomal recessive disorders is still not satisfied for clinical use yet. Luo et al. adopts pseudo-tetraploid model instead of haplotype-based method (e.g., RHDO) to estimate the fetal genotype at specific locus, and achieved 86.4% accuracy for screening three single-gene recessive disorders . Thus it is worth to develop assay and bioinformatics algorithm to enable genetic diagnosis of autosomal recessive disorders in prenatal period. Nevertheless, at current time, the combination of carrier screening for married couple and targeted-seq panel screening for pregnant women to detect mutations in association with recessive disorders in parents and dominant or do novo mutations in NIPT, respectively, shall provide indispensable value to reduce the deleterious effect of pathogenic mutations on human population.
We established an amplicon-based targeted-seq panel that covers 87 pathogenic hotspot mutations reported in autosomal dominant inherited SD. We demonstrated the workflow is not only with high sensitivity and specificity but also shows remarkable concordance between Ion Torrent and Illumina systems. In addition, the assay requires only 1 ng of cfDNA as input material and takes minimal 4 days to accomplish entire workflow. The overall performance is considered fast and cost-effective for routine clinical practice. We believe the established panel provides a promising diagnostic or therapeutic strategy for prenatal genetic testing of SD in routine clinical practice.
Availability of data and materials
All data generated or analyzed during this study are included in the main text and supplementary information files. The raw data are available upon request.
Achondrogenesis, type II
- ALT AF:
Alternative allele frequency
Genome Analysis Toolkit
Noninvasive prenatal testing
Osteogenesis imperfecta, type I
Osteogenesis imperfecta, type II
Osteogenesis imperfecta, type III
Osteogenesis imperfecta, type IV
Single nucleotide polymorphism
Thanatophoric dysplasia, type I
Thanatophoric dysplasia, type II
Mortier GR, Cohn DH, Cormier-Daire V, Hall C, Krakow D, Mundlos S, Nishimura G, Robertson S, Sangiorgi L, Savarirayan R, et al. Nosology and classification of genetic skeletal disorders: 2019 revision. Am J Med Genet A. 2019;179(12):2393–419.
Mogayzel PJ, Marcus CL. Skeletal dysplasias and their effect on the respiratory system. Paediatr Respir Rev. 2001;2(4):365–71.
Orioli IM, Castilla EE, Barbosa-Neto JG. The birth prevalence rates for the skeletal dysplasias. J Med Genet. 1986;23(4):328–32.
Stoll C, Dott B, Roth MP, Alembik Y. Birth prevalence rates of skeletal dysplasias. Clin Genet. 1989;35(2):88–92.
Rasmussen SA, Bieber FR, Benacerraf BR, Lachman RS, Rimoin DL, Holmes LB. Epidemiology of osteochondrodysplasias: changing trends due to advances in prenatal diagnosis. Am J Med Genet. 1996;61(1):49–58.
Camera G, Mastroiacovo P. Birth prevalence of skeletal dysplasias in the Italian Multicentric Monitoring System for Birth Defects. Prog Clin Biol Res. 1982;104:441–9.
Andersen PE Jr. Prevalence of lethal osteochondrodysplasias in Denmark. Am J Med Genet. 1989;32(4):484–9.
Waller DK, Correa A, Vo TM, Wang Y, Hobbs C, Langlois PH, Pearson K, Romitti PA, Shaw GM, Hecht JT. The population-based prevalence of achondroplasia and thanatophoric dysplasia in selected regions of the US. Am J Med Genet A. 2008;146A(18):2385–9.
Connor JM, Connor RA, Sweet EM, Gibson AA, Patrick WJ, McNay MB, Redford DH. Lethal neonatal chondrodysplasias in the West of Scotland 1970–1983 with a description of a thanatophoric, dysplasialike, autosomal recessive disorder, Glasgow variant. Am J Med Genet. 1985;22(2):243–53.
Krakow D, Alanay Y, Rimoin LP, Lin V, Wilcox WR, Lachman RS, Rimoin DL. Evaluation of prenatal-onset osteochondrodysplasias by ultrasonography: a retrospective and prospective analysis. Am J Med Genet A. 2008;146A(15):1917–24.
Tretter AE, Saunders RC, Meyers CM, Dungan JS, Grumbach K, Sun CC, Campbell AB, Wulfsberg EA. Antenatal diagnosis of lethal skeletal dysplasias. Am J Med Genet. 1998;75(5):518–22.
Kallen B, Knudsen LB, Mutchinick O, Mastroiacovo P, Lancaster P, Castilla E, Robert E. Monitoring dominant germ cell mutations using skeletal dysplasias registered in malformation registries: an international feasibility study. Int J Epidemiol. 1993;22(1):107–15.
Parilla BV, Leeth EA, Kambich MP, Chilis P, MacGregor SN. Antenatal detection of skeletal dysplasias. J Ultrasound Med 2003;22(3):255–58.
Doray B, Favre R, Viville B, Langer B, Dreyfus M, Stoll C. Prenatal sonographic diagnosis of skeletal dysplasias. A report of 47 cases. Ann Genet. 2000;43(3–4):163–9.
Krakow D, Lachman RS, Rimoin DL. Guidelines for the prenatal diagnosis of fetal skeletal dysplasias. Genet Med. 2009;11(2):127–33.
Lo YM, Corbetta N, Chamberlain PF, Rai V, Sargent IL, Redman CW, Wainscoat JS. Presence of fetal DNA in maternal plasma and serum. Lancet. 1997;350(9076):485–7.
Lo YM, Tein MS, Lau TK, Haines CJ, Leung TN, Poon PM, Wainscoat JS, Johnson PJ, Chang AM, Hjelm NM. Quantitative analysis of fetal DNA in maternal plasma and serum: implications for noninvasive prenatal diagnosis. Am J Hum Genet. 1998;62(4):768–75.
Quezada MS, Gil MM, Francisco C, Orosz G, Nicolaides KH. Screening for trisomies 21, 18 and 13 by cell-free DNA analysis of maternal blood at 10–11 weeks’ gestation and the combined test at 11–13 weeks. Ultrasound Obstet Gynecol. 2015;45(1):36–41.
Hou Y, Yang J, Qi Y, Guo F, Peng H, Wang D, Wang Y, Luo X, Li Y, Yin A. Factors affecting cell-free DNA fetal fraction: statistical analysis of 13,661 maternal plasmas for non-invasive prenatal screening. Hum Genomics. 2019;13(1):62.
Wang E, Batey A, Struble C, Musci T, Song K, Oliphant A. Gestational age and maternal weight effects on fetal cell-free DNA in maternal plasma. Prenat Diagn. 2013;33(7):662–6.
Petersen AK, Cheung SW, Smith JL, Bi W, Ward PA, Peacock S, Braxton A, Van Den Veyver IB, Breman AM. Positive predictive value estimates for cell-free noninvasive prenatal screening from data of a large referral genetic diagnostic laboratory. Am J Obstet Gynecol. 2017;217(6):691.e1–691.e6.
Lo KK, Karampetsou E, Boustred C, McKay F, Mason S, Hill M, Plagnol V, Chitty LS. Limited clinical utility of non-invasive prenatal testing for subchromosomal abnormalities. Am J Hum Genet. 2016;98(1):34–44.
Chitty LS, Mason S, Barrett AN, McKay F, Lench N, Daley R, Jenkins LA. Non-invasive prenatal diagnosis of achondroplasia and thanatophoric dysplasia: next-generation sequencing allows for a safer, more accurate, and comprehensive approach. Prenat Diagn. 2015;35(7):656–62.
Terasawa S, Kato A, Nishizawa H, Kato T, Yoshizawa H, Noda Y, Miyazaki J, Ito M, Sekiya T, Fujii T, et al. Multiplex PCR in noninvasive prenatal diagnosis for FGFR3-related disorders. Congenit Anom (Kyoto). 2019;59(1):4–10.
Zhang J, Li J, Saucier JB, Feng Y, Jiang Y, Sinson J, McCombs AK, Schmitt ES, Peacock S, Chen S, et al. Non-invasive prenatal sequencing for multiple Mendelian monogenic disorders using circulating cell-free fetal DNA. Nat Med. 2019;25(3):439–47.
Shin H, Sa JK, Bae JS, Koo H, Jin S, Cho HJ, Choi SW, Kyoung JM, Kim JY, Seo YJ, et al. Clinical Targeted Next-Generation sequencing Panels for Detection of Somatic Variants in Gliomas. Cancer Res Treat. 2020;52(1):41–50.
Kroigard AB, Thomassen M, Laenkholm AV, Kruse TA, Larsen MJ. Evaluation of nine somatic variant callers for detection of somatic mutations in exome and targeted deep sequencing data. PLoS ONE. 2016;11(3):e0151664.
Verhoef TI, Hill M, Drury S, Mason S, Jenkins L, Morris S, Chitty LS. Non-invasive prenatal diagnosis (NIPD) for single gene disorders: cost analysis of NIPD and invasive testing pathways. Prenat Diagn. 2016;36(7):636–42.
Malcher C, Yamamoto GL, Burnham P, Ezquina SAM, Lourenco NCV, Balkassmi S, Antonio DSM, Hsia GSP, Gollop T, Pavanello RC, et al. Development of a comprehensive noninvasive prenatal test. Genet Mol Biol. 2018;41(3):545–54.
Krakow D, Rimoin DL. The skeletal dysplasias. Genet Med. 2010;12(6):327–41.
Krakow D. Skeletal dysplasias. Clin Perinatol 2015;42(2):301–319.
Wenger T, Miller D, Evans K. FGFR Craniosynostosis syndromes overview. In: Adam MP, Ardinger HH, Pagon RA, Wallace SE, Bean LJH, Stephens K, Amemiya A, editors. GeneReviews® [Internet]. Seattle (WA); 1998.
Goncalves L, Jeanty P. Fetal biometry of skeletal dysplasias: a multicentric study. J Ultrasound Med. 1994;13(12):977–85.
Xue Y, Sun A, Mekikian PB, Martin J, Rimoin DL, Lachman RS, Wilcox WR. FGFR3 mutation frequency in 324 cases from the International Skeletal Dysplasia Registry. Mol Genet Genomic Med. 2014;2(6):497–503.
French T, Savarirayan R. Thanatophoric dysplasia. In: Adam MP, Ardinger HH, Pagon RA, Wallace SE, Bean LJH, Stephens K, Amemiya A, editors. GeneReviews® [Internet]. Seattle (WA); 2004.
Legare JM. Achondroplasia. In: Adam MP, Ardinger HH, Pagon RA, Wallace SE, Bean LJH, Stephens K, Amemiya A, editors. GeneReviews® [Internet]. Seattle (WA); 1998.
Ornitz DM, Legeai-Mallet L. Achondroplasia: Development, pathogenesis, and therapy. Dev Dyn. 2017;246(4):291–309.
Steiner RD, Basel D. COL1A1/2 Osteogenesis Imperfecta. In: Adam MP, Ardinger HH, Pagon RA, Wallace SE, Bean LJH, Stephens K, Amemiya A, editors. GeneReviews® [Internet]. Seattle (WA); 2005.
Dalgleish R. The human type I collagen mutation database. Nucleic Acids Res. 1997;25(1):181–7.
Dalgleish R. The Human Collagen Mutation Database 1998. Nucleic Acids Res. 1998;26(1):253–5.
Hattapoğlu S, Durmaz MS. Radiological features of achondrogenesis type 1A: case report and review of the literature. Med J Obstet Gynecol. 2018;6(1):1115.
Wenger TL, Hing AV, Evans KN. Apert Syndrome. In: Adam MP, Ardinger HH, Pagon RA, Wallace SE, Bean LJH, Stephens K, Amemiya A, editors. GeneReviews® [Internet]. Seattle (WA); 2019.
Lajeunie E, Cameron R, El Ghouzzi V, de Parseval N, Journeau P, Gonzales M, Delezoide AL, Bonaventure J, Le Merrer M, Renier D. Clinical variability in patients with Apert’s syndrome. J Neurosurg. 1999;90(3):443–7.
Al-Namnam NM, Hariri F, Thong MK, Rahman ZA. Crouzon syndrome: Genetic and intervention review. J Oral Biol Craniofac Res. 2019;9(1):37–9.
Lin M, Lu Y, Sui Y, Zhao N, Jin Y, Yi D, Jiang M. Extremely severe scoliosis, heterotopic ossification, and osteoarthritis in a three-generation family with Crouzon syndrome carrying a mutant c.799T>C FGFR2. Mol Genet Genomic Med. 2019;7(9):e843.
Lajeunie E, Heuertz S, El Ghouzzi V, Martinovic J, Renier D, Le Merrer M, Bonaventure J. Mutation screening in patients with syndromic craniosynostoses indicates that a limited number of recurrent FGFR2 mutations accounts for severe forms of Pfeiffer syndrome. Eur J Human Genet EJHG. 2006;14(3):289–98.
Kan SH, Elanko N, Johnson D, Cornejo-Roldan L, Cook J, Reich EW, Tomkins S, Verloes A, Twigg SR, Rannan-Eliya S, et al. Genomic screening of fibroblast growth-factor receptor 2 reveals a wide spectrum of mutations in patients with syndromic craniosynostosis. Am J Hum Genet. 2002;70(2):472–86.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60.
Cibulskis K, Lawrence MS, Carter SL, Sivachenko A, Jaffe D, Sougnez C, Gabriel S, Meyerson M, Lander ES, Getz G. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat Biotechnol. 2013;31(3):213–9.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. Genome Project Data Processing S: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Lin JC, Fan CT, Liao CC, Chen YS. Taiwan Biobank: making cross-database convergence possible in the Big Data era. Gigascience. 2018;7(1):1–4.
Tang YA, Wang LY, Chang CM, Lee IW, Tsai WH, Sun HS. Novel Compound Heterozygous Mutations in CRTAP Cause Rare Autosomal Recessive Osteogenesis Imperfecta. Front Genet. 2020;11:897.
Dello Russo C, Cesta A, Longo S, Barone MA, Cima A, Mesoraca A, Sparacino D, Viola A, Giorlandino C. Validation of extensive next-generation sequencing method for monogenic disorder analysis on cell-free fetal DNA: noninvasive prenatal diagnosis. J Mol Diagn. 2019;21(4):572–9.
Perlado S, Bustamante-Aragones A, Donas M, Lorda-Sanchez I, Plaza J, Rodriguez de Alba M. Fetal genotyping in maternal blood by digital PCR: towards NIPD of monogenic disorders independently of parental origin. PLoS ONE. 2016;11(4):e0153258.
Chitty LS, Khalil A, Barrett AN, Pajkrt E, Griffin DR, Cole TJ. Safe, accurate, prenatal diagnosis of thanatophoric dysplasia using ultrasound and free fetal DNA. Prenat Diagn. 2013;33(5):416–23.
New MI, Tong YK, Yuen T, Jiang P, Pina C, Chan KC, Khattab A, Liao GJ, Yau M, Kim SM, et al. Noninvasive prenatal diagnosis of congenital adrenal hyperplasia using cell-free fetal DNA in maternal plasma. J Clin Endocrinol Metab. 2014;99(6):E1022-1030.
Parks M, Court S, Cleary S, Clokie S, Hewitt J, Williams D, Cole T, MacDonald F, Griffiths M, Allen S. Non-invasive prenatal diagnosis of Duchenne and Becker muscular dystrophies by relative haplotype dosage. Prenat Diagn. 2016;36(4):312–20.
Jenkins LA, Deans ZC, Lewis C, Allen S. Delivering an accredited non-invasive prenatal diagnosis service for monogenic disorders and recommendations for best practice. Prenat Diagn. 2018;38(1):44–51.
Carothers AD, McAllion SJ, Paterson CR. Risk of dominant mutation in older fathers: evidence from osteogenesis imperfecta. J Med Genet. 1986;23(3):227–30.
Goriely A, Hansen RM, Taylor IB, Olesen IA, Jacobsen GK, McGowan SJ, Pfeifer SP, McVean GA, Rajpert-De Meyts E, Wilkie AO. Activating mutations in FGFR3 and HRAS reveal a shared genetic origin for congenital disorders and testicular tumors. Nat Genet. 2009;41(11):1247–52.
Goriely A, McVean GA, Rojmyr M, Ingemarsson B, Wilkie AO. Evidence for selective advantage of pathogenic FGFR2 mutations in the male germ line. Science. 2003;301(5633):643–6.
Goriely A, Wilkie AO. Paternal age effect mutations and selfish spermatogonial selection: causes and consequences for human disease. Am J Hum Genet. 2012;90(2):175–200.
Lemyre E, Azouz EM, Teebi AS, Glanc P, Chen MF. Bone dysplasia series. Achondroplasia, hypochondroplasia and thanatophoric dysplasia: review and update. Can Assoc Radio J. 1999;50(3):185–97.
Shinde DN, Elmer DP, Calabrese P, Boulanger J, Arnheim N, Tiemann-Boege I. New evidence for positive selection helps explain the paternal age effect observed in achondroplasia. Hum Mol Genet. 2013;22(20):4117–26.
Dan S, Yuan Y, Wang Y, Chen C, Gao C, Yu S, Liu Y, Song W, Asan, Zhu H, et al. Non-invasive prenatal diagnosis of lethal skeletal dysplasia by targeted capture sequencing of maternal plasma. PLoS ONE. 2016;11(7):e0159355.
Luo Y, Jia B, Yan K, Liu S, Song X, Chen M, Jin F, Du Y, Wang J, Hong Y, et al. Pilot study of a novel multi-functional noninvasive prenatal test on fetus aneuploidy, copy number variation, and single-gene disorder screening. Mol Genet Genomic Med. 2019;7(4):e00597.
The authors would like to thank Center for Bioinformatics and Digital Health at the National Cheng Kung University for providing assistances in NGS data management and analyses.
About this supplement
This article has been published as part of BMC Medical Genomics Volume 14 Supplement 3 2021: Selected articles from the 19th Asia Pacific Bioinformatics Conference (APBC 2021): medical genomics The full contents of the supplement are available at https://bmcmedgenomics.biomedcentral.com/articles/supplements/volume-14-supplement-3.
This study was supported by grants MOST 109-2321-B-006-013 to HSS and MOST 109-2326-B-006-007-MY3 to YAT from the Ministry of Science and Technology, Taiwan. These funding sources had no role in the design, analyses, and interpretation of the data in this study. Publication costs are funded by MOST 109-2326-B-006-007-MY3 to YAT from the Ministry of Science and Technology, Taiwan.
Ethics approval and consent to participate
Written and signed informed consent was obtained from all subjects or their legal guardians of this study and approved by the IRB of National Cheng Kung University Hospital (A-BR-007-003).
Consent for publication
The authors declare that there is no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wang, CY., Tang, YA., Lee, IW. et al. Development and validation of an expanded targeted sequencing panel for non-invasive prenatal diagnosis of sporadic skeletal dysplasia. BMC Med Genomics 14 (Suppl 3), 212 (2021). https://doi.org/10.1186/s12920-021-01063-1