ORIGINAL RESEARCH article

Front. Genet., 22 January 2024

Sec. Human and Medical Genomics

Volume 14 - 2023 | https://doi.org/10.3389/fgene.2023.1301439

Prenatal diagnosis of fetuses with ultrasound anomalies by whole-exome sequencing in Luoyang city, China

  • YW

    Yanan Wang 1† *

  • FY

    Fan Yin 2

  • YC

    Yuqiong Chai 1

  • JJ

    Jiapei Jin 1

  • PZ

    Pai Zhang 1

  • QT

    Qianqian Tan 2

  • ZC

    Zhigang Chen 2

  • 1. Department of Genetics and Prenatal Diagnosis, Luoyang Maternal and Child Health Hospital, Luoyang, China

  • 2. Puluo (Wuhan) Medical Biotechnology Co., LTD., Wuhan, China

Article metrics

View details

13

Citations

4,4k

Views

1,4k

Downloads

Abstract

Background: There is a great obstacle in prenatal diagnosis of fetal anomalies due to their considerable genetic and clinical heterogeneity. Whole-exome sequencing (WES) has been confirmed as a successful option for genetic diagnosis in pediatrics, but its clinical utility for prenatal diagnosis remains to be limited.

Methods: A total of 60 fetuses with abnormal ultrasound findings underwent karyotyping or chromosomal microarray analysis (CMA), and those with negative results were further subjected to WES. The identified variants were classified as pathogenic or likely pathogenic (P/LP) and the variant of uncertain significance (VUS). Pregnancy outcomes were obtained through a telephone follow-up.

Results: Twelve (20%, 12/60) fetuses were diagnosed to have chromosomal abnormalities using karyotyping or CMA. Of the remaining 48 cases that underwent WES, P/LP variants were identified in 14 cases (29.2%), giving an additional diagnostic yield of 23.3% (14/60). The most frequently affected organ referred for prenatal WES was the head or neck system (40%), followed by the skeletal system (39.1%). In terms of pathogenic genes, FGFR3 was the most common diagnostic gene in this cohort. For the first time, we discovered five P/LP variants involved in SEC24D, FIG4, CTNNA3, EPG5, and PKD2. In addition, we identified three VUSes that had been reported previously. Outcomes of pregnancy were available for 54 cases, of which 24 cases were terminated.

Conclusion: The results confirmed that WES is a powerful tool in prenatal diagnosis, especially for fetuses with ultrasonographic anomalies that cannot be diagnosed using conventional prenatal methods. Additionally, newly identified variants will expand the phenotypic spectrum of monogenic disorders and greatly enrich the prenatal diagnostic database.

1 Introduction

Congenital anomalies, ranging from relatively minor to severe multi-system anomalies, affect approximately 3 of 100 live births and are responsible for 21% of perinatal deaths (Ely Driscoll, 2020). Identification of the cause of congenital anomalies is important for genetic counseling and clinical management. Ultrasonography, which provides information on the presence or absence of bones, fetal size, shape, and position, is routinely used in prenatal care (Dhamankar et al., 2020; Macedo et al., 2022). However, ultrasound has difficulties in clearly distinguishing between different types of anomalies, resulting in limited prenatal phenotype information.

Typically, fetuses with ultrasonographic anomalies can be associated with all types of genetic variation, such as chromosome aneuploidies, copy number variations, and single-base mutation. Generally, G-band karyotyping that detects chromosome aneuploidies and unbalanced rearrangements (>5–10 Mb) is recommended as the first-line test for prenatal diagnosis with a diagnostic yield of 32% in fetuses with a structural abnormality. In addition, chromosomal microarray analysis (CMA) that identifies smaller microdeletions and duplications is also offered to improve diagnosis by up to 6% over G-banded karyotyping (Wapner et al., 2012). However, only a proportion of fetal anomalies have clear molecular pathogenesis, and the cause of a large number of fetuses with abnormal ultrasound remains elusive, resulting in challenges in prenatal counseling.

Next-generation sequencing (NGS) technology has been widely used for identification of causative genes in various genetic disorders, including intellectual disability (Vrijenhoek et al., 2018), inherited peripheral neuropathies (Hartley et al., 2018), and epilepsy (Symondsand McTague, 2020). Most recently, a clinical guideline of the American College of Medical Genetics and Genomics (ACMG) highly recommended exome and genome sequencing as a first- or second-tier test for patients with developmental delay or intellectual disability (Manickam et al., 2021). Whole-exome sequencing (WES) plays an important role in postnatal diagnosis, while only recently has WES obtained importance in prenatal diagnosis. Carss et al. (2014), who first introduced the WES strategy into prenatal diagnosis in a cohort of 30 non-aneuploid fetuses, identified 35 de novo single-nucleotide variants, small indels, deletions, or duplications. Causative variants were identified in three out of 30 cases with a diagnostic yield of 10%. Several prior studies have examined the use of prenatal WES, where either karyotype analysis or CMA yielded negative results. These studies reported diagnostic yields ranging from 9.1% to 45.9% (Leung et al., 2018; Sparks et al., 2020; Vora et al., 2020; Huang et al., 2023). These findings highlight the significant contribution of WES toward improving the prenatal diagnostic rates.

Despite the numerous advances in WES technology, conclusions drawn from the majority of pilot studies are biased by small cohorts (Cao et al., 2022; Yang et al., 2022) or are limited to highly selected cases or other certain conditions (Liu et al., 2019; Lei et al., 2022). Further pilot studies involving different disorders, conditions, sample types, and analytic methods are necessary to expand our understanding in the application of WES on prenatal diagnosis and to establish the diagnostic yield of WES for all kinds of prenatal disorders.

Here, we present 60 prenatal cases with various ultrasonographic anomalies. All pregnant women underwent invasive testing, either via karyotyping or CMA or both. Cases with negative results of karyotyping or CMA tests were further subjected to WES for causative variant identification, with the aim of discovering more pathogenic mutations and clarifying the utility of WES in prenatal diagnosis.

2 Materials and methods

2.1 Study cohort

This study was granted by the Institutional Review Board of the Ethics Committee at the Luoyang Maternal and Child Health Hospital (approval number LYFY-YCCZ-2023006). Pregnant women were recruited at our hospital between November 2019 and December 2022, from which 60 pregnancies with ultrasonographic anomalies that underwent invasive diagnosis were selected for further investigation. All pregnant women and their partners signed written informed consent for invasive procedures, testing, and participation. Phenotypic information and disease classification were provided by the clinicians, according to the ultrasonographic results based on the Human Phenotype Ontology (HPO) and the Online Mendelian Inheritance in Man (OMIM) databases. If a case had multiple-organ abnormalities, then each abnormality was counted separately. The study profile is shown in Figure 1.

FIGURE 1

2.2 Chromosome karyotyping and chromosomal microarray analysis

Amniotic fluid samples (10 mL) were collected from pregnant women and subjected to conventional G-band karyotyping, according to standard operation procedures as previously described (Shi et al., 2019). For CMA analysis, genomic DNA was extracted from amniotic fluid using a QIAamp DNA Mini Kit (QIAGEN Inc., Valencia, CA, United States), followed by amplification, labeling, and hybridization, and then, a CMA test was performed using a 750K microarray chip (Affymetrix Inc., United States), according to the manufacturer’s protocol. The classification of pathogenic copy number variations (CNVs) was determined according to international guidelines of the ACMG.

2.3 Whole-exome sequencing

Fetuses with normal karyotype analysis or CMA results were further subjected to WES for causative variant identification. As previously described (Lippa et al., 2021), WES procedures including library construction, exome capture, and high-throughput sequencing were performed. In brief, 150 ng of genomic DNA was fragmented (250–300 bp), and then, library preparation, including end repair, A-tailing, adapter ligation, and PCR amplification, was conducted using the IDT xGen Exome Research Panel v1.0 (Integrated DNA Technologies, San Diego, United States), according to the manufacturer’s protocols. After getting captured by VCRome version 2.1, all exome libraries were pooled and sequenced using the NovaSeq6000 PE150 platform (Illumina, San Diego, United States), with a target sequencing coverage of 100× and with over 95% of the exonic positions with the depth of coverage >20×.

The QC of the paired-end reads was assessed using FastQC, and raw sequence data were post-progressed on-site using open-source software. After filtering using SAMtools, the sequencing reads were aligned to the human genome reference sequence (hg19/GRCh37) using the Burrows–Wheeler Aligner (BMA, version 0.59). GATK HaplotypeCaller v3.6 was used to detect single-nucleotide variants (SNVs) and indels (Van der Auwera et al., 2013), and then, the variants were annotated using ANNOVAR (Monday, 8 June 2020). Small chromosomal CNVs and single-gene CNVs were analyzed using CNVkit and annotated using AnnotSV. In the present study, the clinical significance of identified variants was classified into the following three categories, namely, pathogenic (P), likely pathogenic (LP), and the variant of uncertain significance (VUS), according to the ACMG guidelines (Green et al., 2013), based on the allele frequency, family segregation, compatibility with the phenotype, in silico prediction, relevant disease databases, and the literature. The allele frequency ≤0.07% in the gnomAD database was considered a variant of clinical significance. Moreover, the Human Gene Mutation Database (HGMD) was used to assess whether the variant identified in the present study had been previously reported.

2.4 Sanger sequencing

All variants considered to be causative for the observed fetal phenotype were validated by Sanger sequencing using both fetal amniotic fluid DNA and parental peripheral blood samples. Further clinical tests were performed, and additional family history information was collected to facilitate interpretation of the Sanger sequencing results.

2.5 The pregnancy outcome

Pregnancy outcomes were collected through follow-up telephone calls until August 2023, including the condition and phenotype of live-born infants. The phenotype of the parents was also recorded.

3 Results

3.1 Detection of fetal chromosomal abnormalities

A total of 60 pregnant women of Han ethnicity with ultrasonographic anomalies were referred to the Luoyang Maternal and Child Health Hospital. Supplementary Table S1 displays the demographics of the pregnant women, revealing that they had a median maternal age of 29.5 ± 4.3 years (range 21–39) and a median gestational age of 23.4 ± 4.9 weeks (range 12–33). Invasive procedures were conducted for all subjects, of which 23 cases received CMA and the other 37 cases received both G-band karyotyping and CMA. Totally, 12 out of 60 (20%) cases were found to have chromosomal abnormalities. One case was a fetus with Down’s syndrome, and another had mosaic aneuploidy. The remaining 10 fetuses with abnormal copy number variants ranging from 150 kb to 4.3 Mb were mainly involved in Duchenne muscular dystrophy, intellectual disability, X-linked 93, glass syndrome, hereditary hearing loss, and renal hypoplasia. All the copy number variants were classified to be pathogenic or likely pathogenic or the variant of uncertain significance, according to ACMG (Table 1).

TABLE 1

CaseMain ultrasound findingMicroarray result (ISCN)SizeInheritance or de novoAssociated disorderbClinGen Dosage ID (score)ACMG classification (Green et al., 2013)KaryotypingPregnancy outcomea
1Renal dysplasia and macrocephalyarr[GRCh37]9p22.3 (14803711_14907432)×1103.7 kbDe novoTrigonocephaly 2ISCA-23907 (HI 0)Likely pathogenic (4A)NATOP
PMID:21931569
2Short femur and seroperitoneumarr[GRCh37]18p11.21-q11.2 (15260001_19560000)×14.3 MbUnknownHereditary hearing lossPMID:29100090Pathogenic (2A)NATOP
Renal hypoplasiaPMID:32378186 c
3Hydronephrosis and ventricular septal defectarr[GRCh37]2q33.1 (199379852_201279851)×11.9 MbDe novoGlass syndromeISCA-31904 (HI 3)Pathogenic (2A)46,XXLive birth with groan breathing
4Absent fetal nasal bone and ventricular septal defectNANANADown’s syndromeNAPathogenic (2A)47,XY,+21TOP
5Polyhydramnios, short femur, and abnormality of the gallbladderarr[GRCh37]16p11.2 (28760001_29010000)×1250 KbUnknownChromosome 16p11.2 deletion syndromeISCA-37486 (HI 3)Pathogenic (2A)NATOP
6Increased nuchal translucencyarr[GRCh37]Xq21.1 (80060001_80210000)×0150 KbInherited maternallyIntellectual disability and X-linked 93ISCA-16544 (HI 3)Likely pathogenic (2C-1)46,XYLive birth without any abnormality
7Duodenal ileusarr[GRCh37]16p11.2 (29560001_30210000)×1650 kbUnknownChromosome 16p11.2 deletion syndromeISCA-37400 (HI 3)Pathogenic (2A)46,XXTOP
8Imperforate atrioventricular valve and pulmonary artery atresiaarr[GRCh37]22q11.2 (22250005_22550004)×3300 kbUnknown22q11.2 deletion syndromeNAVUSNATOP
9Partial absence of the cerebellar vermisNANANAMosaic trisomy 2NAPathogenic (2A)47,XY,+2[11]/46,XY[46]Live birth without abnormality
10Short femurarr[GRCh37]Xp21.1 (31860001_32110000)×1250 kbUnknownDuchenne muscular dystrophyISCA-6635 (HI 3)Pathogenic (PVS + PS4+PM2__Supporting)46,XXLive birth without abnormality
11Increased nuchal translucencyarr[GRCh37]16p13.11 (15510001_16360000)×3850 kbInherited paternally16p13.11 recurrent regionISCA-37415 (TS 2)Pathogenic (2A)46,XYLive birth without abnormality
12Short fetal length and bilateral renal echo enhancementarr[GRCh37]17q12 (34800001_36250000)×11.45 MbUnknownChromosome 17q12 deletion syndromeISCA-37432 (HI 3)Pathogenic (2A)NATOP

Characteristics of fetuses with copy number variations detected by chromosomal microarray analysis.

“TOP” denotes the termination of pregnancy.

“NA” denotes not applicable.

a

After termination, the fetal samples were not used for further diagnosis.

b

Associated disorder was determined according to OMIM.

c

LOF, variants of the involved gene GRE1BL, have been reported to be pathogenic (PMID: 29100090, PMID: 32378186, ClinVar ID: 2445432, et al.). In addition, two or more haploinsufficiency predictors suggest it is haploinsufficient.

3.2 Positive diagnostic results identified by WES

After exclusion of 12 fetuses with chromosome abnormalities, the remaining 48 amniotic fluid DNA samples were then subjected to WES. As shown in Table 2, 14 (29.2%) cases were found to be harboring P/LP variants. Accordingly, WES increased the diagnostic yield by further 14 cases (23.3%) above routine genetic testing. In total, 13 P/LP variants were identified in 14 fetuses corresponding to 11 unique genes. Seven (50%, 7/14) fetuses had a de novo mutation, with FGFR3 being the most frequent diagnostic gene carrying the same mutation c.1138G>A (p.G380R) that was identified in four cases. Four (28.6%, 4/14) fetuses inherited the mutation from parents who were heterozygous carriers in an autosomal recessive form, including compound heterozygous SEC24D, EPG5, and FIG4, as well as a homozygous P3H1. The remaining three fetuses (21.4%, 3/14) inherited the mutation in an autosomal dominant way from one of their parents (HNF1B, CRYAA, and CTNNA3). Of the 13 P/LP variants, we have identified five variants that are not currently included in the HGMD. These variants included a de novo missense mutation in the PKD2 gene in fetuses suspected of causing polycystic kidney dysplasia; a maternally transmitted FIG4 canonical splice variant; a maternally transmitted EPG5 missense variant linked to vici syndrome; a maternally transmitted CTNNA3 canonical splice variant associated with arrhythmogenic right ventricular dysplasia; and a maternally transmitted SEC24D canonical splice variant, which may lead to Cole–Carpenter syndrome 2.

TABLE 2

CaseMain ultrasound findingGeneAssociated disordercAlterationVariant typeHGMD inclusionFunction predictionInheritance/zygosityACMG classification (Green et al., 2013)Pregnancy outcomeb
REVEL/spliceAI effect
13Lateral femoral bowingSEC24DCole–Carpenter syndrome 2NM_001318066.2 c.2241 + 1G>TCanonical splice-NA/YesInherited/compound heterozygousLikely pathogenic (PVS1+PM2_Supporting)Born with large fontanelles
NM_001318066.2 c.941G>A (p.R314H)MissenseCM16177050.685/NALikely pathogenic (PM3_Strong + PM2_Supporting + PP3)
14Spot in the left ventricle and bilateral femoral bowingP3H1Osteogenesis imperfectaNM_001243246.2 c.2041C>T (p.R681a)NonsenseCM086851NA/NAInherited/homozygousLikely pathogenic (PVS_Moderate + PS3_Supporting + PM2_Supporting + PM3)TOP
15Bilateral talipes equinovarusNALCNCongenital contractures of the limbs, face, hypotonia, and developmental delayNM_001350748.2 c.1733A>G (p.Y578C)MissenseCM1768730.976/NADe novo/heterozygousPathogenic (PS2_Very strong + PM2_Supporting, PP3)Born with bilateral talipes equinovarus combined with curved fingers
16Bilateral hyperechogenic kidneysHNF1BaRenal cysts and diabetes syndromeNM_000458.4 c.809 + 1G>ACanonical spliceCS104578NA/YesInherited paternally/heterozygousLikely pathogenic (PVS + PM2_Supporting)Born with bilateral cortical cysts
17Enhanced parenchymal echo in both the kidneysPKD2Polycystic kidney dysplasiaNM_000297.4 c.1034A>G (p.Y345C)Missense-0.855/NADe novo/heterozygousLikely pathogenic (PM1+PM2_Supporting + PM6+PP3)Lost to follow-up
18Lens opacityCRYAAaCongenital cataractNM_000394.4 c.160C>T (p.R54C)MissenseCN0761300.824/NAInherited paternally/heterozygousLikely pathogenic (PS3+PS4__Supporting + PM2_Supporting + PP3)Lost to follow-up
19Short femur and small belly circumferenceFGFR3AchondroplasiaNM_000142.5 c.1138G>A (p.G380R)MissenseCM9407850.696/NADe novo/heterozygousPathogenic ((PS2+PS3_Supporting + PM1+PM5+PM2_Supporting))TOP
20Short femur, short humerus, and polyhydramniosFGFR3AchondroplasiaNM_000142.5 c.1138G>A (p.G380R)MissenseCM9407850.696/NADe novo/heterozygousPathogenic (PS2+PS3_Supporting + PM1+PM5+PM2_Supporting)TOP
21Skeletal dysplasiaFGFR3AchondroplasiaNM_000142.5 c.1138G>A (p.G380R)MissenseCM9407850.696/NADe novo/heterozygousPathogenic (PS2+PS3_Supporting + PM1+PM5+PM2_Supporting)TOP
22Skeletal dysplasiaFGFR3AchondroplasiaNM_000142.5 c.1138G>A (p.G380R)MissenseCM9407850.696/NADe novo/heterozygousPathogenic (PS2+PS3_Supporting + PM1+PM5+PM2_Supporting)TOP
23Double aortic arch and patent ductus arteriosusCTNNAaArrhythmogenic right ventricular dysplasia, 13NM_001127384.3 c.1281 + 1G>ACanonical splice-NA/YesInherited maternally/heterozygousLikely pathogenic (PVS + PM2_Supporting)Live birth without abnormality
24Short femur and femoral bowingCOL1A2Osteogenesis imperfectaNM_000089.4 c.821G>A (p.G274D)MissenseCM0112950.983/NADe novo/heterozygousLikely pathogenic (PM1+PM6+PS4_Moderate + PM2_Supporting + PP3)TOP
25LeukodystrophyEPG5Vici syndromeNM_020964.3 c.7333C>T (p.R2445a)NonsenseCM163403NA/NAInherited/compound heterozygousLikely pathogenic (PVS + PM3_Supporting + PM2_Supporting)TOP
NM_020964.3 c.5714G>A (p.R1905Q)Missense-0.657/NALikely pathogenic (PM1+PM3+PM2_Supporting + PP3)
26Dilated ventricle and fetal pyelectasisFIG4Yunis–Varon syndromec.2376 + 1G>TCanonical splice-NA/YesInherited/compound heterozygousLikely pathogenic (PVS + PM2_Supporting)TOP
NM_014845.6 c.350C>T (p.A117V)MissenseCI22151000.664/NAVUS (PM2_Supporting + PM3+PP3)

Summary of pathogenic and likely pathogenic mutations revealed by whole-exome sequencing.

“-” denotes the variant that was not included in the HGMD.

a

denotes the mutation was inherited from the parent, and the disease displays incomplete penetrance or a highly variable phenotype.

b

After termination, the fetal samples were not used for further diagnosis.

“TOP” denotes termination of pregnancy.

c

Associated disorder was determined according to OMIM.

“NA” denotes not applicable.

“Yes” denotes has effect.

“PM2_Supporting” allele frequency <0.07%.

3.3 Secondary findings and variants of uncertain significance (VUSes)

In addition to P/LP variants, we also identified secondary findings associated with GJB2, DUOX2, and TTN genes in two cases (Supplementary Table S2). In addition, we detected 22 VUSes among 15 cases, three of which have been previously reported at least once (Supplementary Table S3). The remaining 17 cases were WES-negative; neither chromosomal CNVs nor single gene CNVs were identified in these negative cases (Supplementary Table S4).

3.4 Diagnostic rates of WES among ultrasonographic anomalies

We classified the fetal anomaly into seven categories, according to the systems involved, including skeletal, cardiovascular, genitourinary, head or neck, central nervous, gastrointestinal, and respiratory systems. Among the different malformation categories, the skeletal system (55.0%, 33/60) showed the most frequent findings in this cohort by ultrasonic testing, followed by cardiovascular (23.3%), genitourinary (18.3%), and head or neck systems (15.0%) (Figure 2A). The results in Figure 2B showed that diagnostic yield varied between different categories; the highest rate was found in the fetus with head or neck abnormalities (40%, 2/5) and skeletal abnormalities (39.1%, 9/23). In addition, the diagnostic rates in cardiovascular, genitourinary, and central nervous systems were 22.2% (2/9), 33.3% (2/6), and 33.3% (1/3), respectively. However, the sample size was too small that the test failed to detect meaningful differences in diagnostic rates for various systemic diseases.

FIGURE 2

3.5 Pregnancy outcomes

Pregnancy outcomes were available in 54 of 60 (90.0%) cases, and the remaining six cases were lost to follow-up. Twenty-four of 54 (44.4%) pregnancy women preferred termination, including seven fetuses identified with pathogenic copy number variations, eight fetuses with P/LP mutations, one fetus with secondary findings, four fetuses with VUSes, and four fetuses with negative results (Supplementary Table S4). It is noteworthy that the fetal samples were not subjected to further diagnosis post-termination. Of the 30 pregnancy women who chose to continue gestation, adverse pregnancy was observed in eight fetuses. In Table 1, the fetus in case 3 that was diagnosed to have glass syndrome was born with groan breathing. Three fetuses carrying the P/LP variant were born with abnormalities (Table 2), especially case 13, who was suspected to have Cole–Carpenter syndrome 2, was born with large fontanelles. Case 16 showed bilateral cortical cysts at birth, and case 15 presented with bilateral talipes equinovarus combined with curved fingers. The results were consistent with both the ultrasonographic findings and phenotype of the NALCN mutation. Interestingly, bilateral talipes equinovarus after birth was also observed in case 42 who was ultimately classified as having a VUS of the SOX9 gene mutation (Supplementary Table S3). Additionally, two fetuses carrying VUS were born with abnormalities including a live birth with polymicrogyria in case 36 and moderate postnatal growth restriction in case 37. Another fetus with negative results was born with slight postnatal growth restriction in case 56. Of the 17 WES negative cases, 10 (58.8%) fetuses were live births with no abnormalities and one (5.9%) with an abnormality (Supplementary Table S4).

4 Discussion

In this study, WES was used to examine fetuses with abnormal ultrasound findings but could not be diagnosed by routine testing. Our findings demonstrated that WES increased the diagnostic yield in fetuses with abnormal sonographic findings by 23.3%, following negative results obtained through karyotype and CMA testing. These results are consistent with previously published series on published prenatal WES, which have reported diagnostic yields ranging from 15% to 25% (Vora et al., 2017; Vora et al., 2020; Lei et al., 2021). Our results further confirmed the potential of WES as a promising tool for extended prenatal diagnosis. The decreasing turnaround time and increasing accessibility, as well as the falling cost of next-generation sequencing will make it more feasible to use the WES strategy clinically. This will remarkably promote parental counseling and pregnancy management. Nevertheless, conventional testing cannot be replaced by WES due to limitations of WES in detecting large insertions/deletions, chromosomal rearrangements, and mutations in regulatory regions (Alotibi et al., 2023). Karyotype analysis or CMA in combination with WES would be a valuable strategy in prenatal testing.

Among seven categories in this cohort, skeletal anomalies and head or neck anomalies were the top two positive predictors of monogenic disorder. Here, we provided a diagnostic rate by WES of 39.1% of fetuses with abnormalities in the skeletal system, which was comparable to previously published prenatal WES studies with 30.4% (Fu et al., 2022) and 32.7% (Kucińska-Chahwan et al., 2022). However, a significantly higher rate of 75% has been reported elsewhere (Liu et al., 2019; Yang et al., 2019). Skeletal anomalies are a heterogeneous group of disorders affecting bones, cartilage, tendons, and joints (Köhler et al., 2017); a highly heterogeneous phenotype may contribute to the differences in detection rates between these studies. In addition, differences in inclusion criteria, sample size, or the definition of a positive result might also be involved. In any case, a higher diagnostic yield would be obtained from these malformations by WES with optimal cost-effectiveness. Therefore, WES is strongly recommended for fetuses with abnormal ultrasound findings of skeletal anomalies and head or neck anomalies.

With regard to pathogenic genes identified by WES, FGFR3 was one of the most frequently diagnosed genes in the current study that is consistent with previous findings (Han et al., 2020; Zhang et al., 2021). However, different variants of FGFR3 were responsible for different disorders with various phenotypes and prognoses that were unable to distinguish from fetal sonographic indicators (Yang et al., 2019). For instance, FGFR3: c.1138G>A (p.G380R) mutation accounted for 90% of the condition in achondroplasia patients (Bellus et al., 1995), while the FGFR3: c.1620C>A (p.N540K) mutation occurred in 70% of individuals with hypochondroplasia (Bellus et al., 2000). Therefore, WES is of importance for diagnosis and genetic counseling to fetuses with variations in identical genes.

A homozygous nonsense variation c.2041C>T (p.R681*) in P3H1 was identified in case 14 with bilateral femoral bowing. The parent denied consanguinity, and the fetus was eventually terminated because this mutation is implicated in osteogenesis imperfecta. The same mutation has been reported in a British family, in which the fetus exhibited abnormal ultrasound findings of short, thick bowed femurs, evidence of fractures, short humerus, and abnormal shape (Chandler et al., 2018). In a consanguineous family of Arabic descent, this mutation was identified in a fetus with severe limbs and chest deformities (Baldridge et al., 2008). All these lines of evidence suggested that the homozygous non-sense variation c.2041C>T (p.R681*) in P3H1 may be a common variant among different races and usually result in severe fetal malformation.

The PKD2 gene, encoding polycystin-2, is associated with polycystic kidney disease, and the clinical symptoms usually do not appear until adulthood, but the disease starts in utero (Janssens et al., 2021). In one study, the author reported on a family carrying a mutation in the PKD2 gene perinatal death due to polycystic kidney disease occurred in the mother’s second and third pregnancies (Bergmann et al., 2008). In case 17, a mutation in the PKD2 gene was identified in a fetus with enhanced parenchymal echo in both the kidneys. This case was lost to follow-up, while the genetic result suggested a high risk of polycystic kidney dysplasia in adulthood; therefore, advanced interfere treatment is a matter of great urgency.

In case 23, the LP variant of CTNNA3 was identified in the fetus with patent ductus arteriosus, but the transmitting mother was not an affected individual. This is largely due to the incomplete penetrance of this mutation in different affected individuals (van Hengel et al., 2013). A maternally transmitted canonical splice variant c.809 + 1G>A in HNF1B was identified in case 16 with bilateral renal cortical cysts after birth, while the same mutation inherited from the father was also reported in a previous study, in which the fetus was born with a unilateral ureteropelvic junction and renal failure, apart from bilateral renal cortical cysts (Heidet et al., 2010). These findings indicated that the severity of the renal disease resulting from the HNF1B mutation was highly variable. In case 18, the fetus inherited a heterozygous missense variation c.160C>T (p.R54C) in the CRYAA gene from the father, and the parents refused to disclose any information about the newborn. Meanwhile, the father reported no notable ocular disease. Khan et al. (2007) reported that three carriers of c.160C>T (p.R54C) in the CRYAA gene were asymptomatic but had similar bilateral discrete punctuate lenticular opacities evident by a careful slit lamp examination (Khan et al., 2007). Thus, in our study, fetal case 18 was highly suspected of congenital cataract, and we recommended further ophthalmic examination to the transmitted father, but it was not adopted. A sustained follow-up is necessary. To the best of our knowledge, five pathogenic or likely pathogenic variants implicated in SEC24D, FIG4, CTNNA3, EPG5, and PKD2 (Table 2) had never been reported prenatally. These newly identified variants demonstrate how prenatal WES expanded the phenotypic spectrum of monogenic disorders in prenatal diagnosis. In a word, regarding the 14 cases with pathogenic or likely pathogenic genes (Table 2), the pregnancy outcomes further supported the main ultrasound findings only in two cases (cases 15 and 16). However, the concordance of the remaining cases was indeterminate, either because of pregnancy termination or loss to follow-up.

The annotation of VUS in the prenatal setting has been debated for years, and it is essential that VUS results are properly explained to avert potential misuse (Cornthwaite et al., 2022). Here, we identified 22 VUSes (Supplementary Table S3), of which three variants (JAG1:c.1136C>T(p.S379F), NBAS:c.6124A>G (p.M2042V), and DYNC2H1:c.7409C>T (p.A2470V)) had been reported previously at least once. We firmly believe that with the promotion of WES in prenatal diagnosis, more VUSes will be identified and reclassified, allowing a more accurate diagnosis will be possible with integrated genotype–phenotype information.

Although WES has many advantages over the traditional test in prenatal diagnosis, some major challenges have been actively debated as prenatal WES clinically is becoming part of routine testing. The lack of an accurate fetal phenotype achieved by ultrasonography, variable disease descriptions, and appearances of ultrasound findings further complicate the interpretation of the WES results (Best et al., 2018). Other challenges including the cost effectiveness of WES (Aaltio et al., 2022) and even social and ethical issues (Horn and Parker, 2018) are also taken into account. There are several limitations in the current study which need to be further improved and supplemented. First, we performed WES only for the proband instead of parental–fetal trios. Since trio-based WES testing would be helpful for variant filtering in initial data analysis and enable rapid identification of de novo variants. In particular, for WES negative cases, trio-based WES is useful for determining the causative genes. Thus, parental–fetal trio analysis is preferable, although it would increase the cost. Second, this is a retrospective study; some of the pregnancy outcomes were obtained through follow-up phone calls rather than confirmation of diagnosis by a clinician; thus, we are not sure whether the parents have given an objective statement of fact. Finally, the relatively small sample size is still a limitation to our study so that we are unable to identify statistical differences of diagnostic rates between different disorders. Therefore, additional efforts are urgent on large cohort studies to better clarify correlations between the fetal genotype and phenotype and provide refined prenatal counseling.

5 Conclusion

In conclusion, we confirmed the potential of WES to improve prenatal diagnosis in a cohort of 60 cases. Our results demonstrated that WES remarkably increased the prenatal diagnosis rate by 23.3% in fetuses with abnormal ultrasound findings but normal G-band karyotyping and CMA results. In addition, we suggested that WES may be recommended when conventional prenatal methods fail to provide a diagnosis in fetuses with ultrasonographic anomalies, particularly for disorders involving the skeletal system and head or neck systems. More importantly, this study identified five P/LP variants that had not been included in the HGMD. Additionally, three VUSes that had been reported at least once previously were recurred in this study. All of these findings further add to our current knowledge of phenotype–genotype relationships, which will improve the ability of genetic counselors to advice families.

Statements

Data availability statement

The datasets presented in this article are not readily available because of privacy restrictions. Requests to access the datasets should be directed to the authors.

Ethics statement

The studies involving humans were approved by the Institutional Review Board of the Ethics Committee in the Luoyang Maternal and Child Health Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.

Author contributions

YW: supervision, writing–original draft, writing–review and editing, conceptualization, formal analysis, funding acquisition, and visualization. FY: investigation, methodology, validation, writing–original draft, conceptualization, data curation, and formal analysis. YC: data curation, formal analysis, funding acquisition, writing–original draft, and project administration. JJ: investigation, methodology, project administration, visualization, and writing–original draft. PZ: conceptualization, formal analysis, validation, writing–original draft, and resources. QT: investigation, methodology, and writing–original draft. ZC: data curation and writing–original draft.

Funding

The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.

Acknowledgments

The authors are very grateful to all the participants who participated in this investigation.

Conflict of interest

FY, QT, and ZC were employed by Puluo (Wuhan) Medical Biotechnology Co., LTD.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1301439/full#supplementary-material

References

  • 1

    AaltioJ.HyttinenV.KortelainenM.FrederixG. W. J.LönnqvistT.SuomalainenA.et al (2022). Cost-effectiveness of whole-exome sequencing in progressive neurological disorders of children. Eur. J. Paediatr. Neurol.36, 3036. 10.1016/j.ejpn.2021.11.006

  • 2

    AlotibiR. S.SannanN. S.AlEissaM.AldriweshM. G.Al TuwaijriA.AkielM. A.et al (2023). The diagnostic yield of CGH and WES in neurodevelopmental disorders. Front. Pediatr.11, 1133789. 10.3389/fped.2023.1133789

  • 3

    BaldridgeD.SchwarzeU.MorelloR.LenningtonJ.BertinT. K.PaceJ. M.et al (2008). CRTAP and LEPRE1 mutations in recessive osteogenesis imperfecta. Hum. Mutat.29 (12), 14351442. 10.1002/humu.20799

  • 4

    BellusG. A.HefferonT. W.Ortiz de LunaR. I.HechtJ. T.HortonW. A.MachadoM.et al (1995). Achondroplasia is defined by recurrent G380R mutations of FGFR3. Am. J. Hum. Genet.56 (2), 368373.

  • 5

    BellusG. A.SpectorE. B.SpeiserP. W.WeaverC. A.GarberA. T.BrykeC. R.et al (2000). Distinct missense mutations of the FGFR3 lys650 codon modulate receptor kinase activation and the severity of the skeletal dysplasia phenotype. Am. J. Hum. Genet.67 (6), 14111421. 10.1086/316892

  • 6

    BergmannC.BrüchleN. O.FrankV.RehderH.ZerresK. (2008). Perinatal deaths in a family with autosomal dominant polycystic kidney disease and a PKD2 mutation. N. Engl. J. Med.359 (3), 318319. 10.1056/NEJMc0801868

  • 7

    BestS.WouK.VoraN.Van der VeyverI. B.WapnerR.ChittyL. S. (2018). Promises, pitfalls and practicalities of prenatal whole exome sequencing. Prenat. Diagn38 (1), 1019. 10.1002/pd.5102

  • 8

    CaoJ.ChenA.TianL.YanL.LiH.ZhouB. (2022). Application of whole exome sequencing in fetal cases with skeletal abnormalities. Heliyon8 (7), e09819. 10.1016/j.heliyon.2022.e09819

  • 9

    CarssK. J.HillmanS. C.ParthibanV.McMullanD. J.MaherE. R.KilbyM. D.et al (2014). Exome sequencing improves genetic diagnosis of structural fetal abnormalities revealed by ultrasound. Hum. Mol. Genet.23 (12), 32693277. 10.1093/hmg/ddu038

  • 10

    ChandlerN.BestS.HaywardJ.FaravelliF.MansourS.KivuvaE.et al (2018). Rapid prenatal diagnosis using targeted exome sequencing: a cohort study to assess feasibility and potential impact on prenatal counseling and pregnancy management. Genet. Med.20 (11), 14301437. 10.1038/gim.2018.30

  • 11

    CórdobaM.Rodriguez-QuirogaS. A.VegaP. A.SalinasV.Perez-MaturoJ.AmartinoH.et al (2018). Whole exome sequencing in neurogenetic odysseys: an effective, cost- and time-saving diagnostic approach. PLoS One13 (2), e0191228. 10.1371/journal.pone.0191228

  • 12

    CornthwaiteM.TurnerK.ArmstrongL.BoerkoelC. F.ChangC.LehmanA.et al (2022). Impact of variation in practice in the prenatal reporting of variants of uncertain significance by commercial laboratories: need for greater adherence to published guidelines. Prenat. Diagn42 (12), 15141524. 10.1002/pd.6232

  • 13

    DhamankarR.DiNonnoW.MartinK. A.DemkoZ. P.Gomez-LoboV. (2020). Fetal sex results of noninvasive prenatal testing and differences with ultrasonography. Obstet. Gynecol.135 (5), 11981206. 10.1097/aog.0000000000003791

  • 14

    ElyD. M.DriscollA. K. (2020). Infant mortality in the United States, 2018: data from the period linked birth/infant death file. Natl. Vital Stat. Rep.69 (7), 118.

  • 15

    FuF.LiR.YuQ.WangD.DengQ.LiL.et al (2022). Application of exome sequencing for prenatal diagnosis of fetal structural anomalies: clinical experience and lessons learned from a cohort of 1618 fetuses. Genome Med.14 (1), 123. 10.1186/s13073-022-01130-x

  • 16

    GreenR. C.BergJ. S.GrodyW. W.KaliaS. S.KorfB. R.MartinC. L.et al (2013). ACMG recommendations for reporting of incidental findings in clinical exome and genome sequencing. Genet. Med.15 (7), 565574. 10.1038/gim.2013.73

  • 17

    HanJ.YangY. D.HeY.LiuW. J.ZhenL.PanM.et al (2020). Rapid prenatal diagnosis of skeletal dysplasia using medical trio exome sequencing: benefit for prenatal counseling and pregnancy management. Prenat. Diagn40 (5), 577584. 10.1002/pd.5653

  • 18

    HartleyT.WagnerJ. D.Warman-ChardonJ.TétreaultM.BradyL.BakerS.et al (2018). Whole-exome sequencing is a valuable diagnostic tool for inherited peripheral neuropathies: outcomes from a cohort of 50 families. Clin. Genet.93 (2), 301309. 10.1111/cge.13101

  • 19

    HeidetL.DecramerS.PawtowskiA.MorinièreV.BandinF.KnebelmannB.et al (2010). Spectrum of HNF1B mutations in a large cohort of patients who harbor renal diseases. Clin. J. Am. Soc. Nephrol.5 (6), 10791090. 10.2215/cjn.06810909

  • 20

    HornR.ParkerM. (2018). Opening Pandora's box? ethical issues in prenatal whole genome and exome sequencing. Prenat. Diagn38 (1), 2025. 10.1002/pd.5114

  • 21

    HuangW.ZhuX.SunG.GaoZ.KongX. (2023). Whole-exome sequencing in deceased fetuses with ultrasound anomalies: a retrospective analysis. BMC Med. Genomics16 (1), 25. 10.1186/s12920-022-01427-1

  • 22

    JanssensP.DecuypereJ. P.De RechterS.BreysemL.Van GielD.BillenJ.et al (2021). Enhanced MCP-1 release in early autosomal dominant polycystic kidney disease. Kidney Int. Rep.6 (6), 16871698. 10.1016/j.ekir.2021.03.893

  • 23

    KhanA. O.AldahmeshM. A.MeyerB. (2007). Recessive congenital total cataract with microcornea and heterozygote carrier signs caused by a novel missense CRYAA mutation (R54C). Am. J. Ophthalmol.144 (6), 949952. 10.1016/j.ajo.2007.08.005

  • 24

    KöhlerS.VasilevskyN. A.EngelstadM.FosterE.McMurryJ.AyméS.et al (2017). The human phenotype Ontology in 2017. Nucleic Acids Res.45 (D1), D865D876. 10.1093/nar/gkw1039

  • 25

    Kucińska-ChahwanA.RoszkowskiT.NowakowskaB.GeremekM.PaczkowskaM.BijokJ.et al (2022). Extended genetic testing in fetuses with sonographic skeletal system abnormalities. Ultrasound Obstet. Gynecol.59 (5), 660667. 10.1002/uog.23722

  • 26

    LeiL.ZhouL.XiongJ. J. (2021). Whole-exome sequencing increases the diagnostic rate for prenatal fetal structural anomalies. Eur. J. Med. Genet.64 (9), 104288. 10.1016/j.ejmg.2021.104288

  • 27

    LeiT. Y.SheQ.FuF.ZhenL.LiR.YuQ. X.et al (2022). Prenatal exome sequencing in fetuses with callosal anomalies. Prenat. Diagn42 (6), 744752. 10.1002/pd.6107

  • 28

    LeungG. K. C.MakC. C. Y.FungJ. L. F.WongW. H. S.TsangM. H. Y.YuM. H. C.et al (2018). Identifying the genetic causes for prenatally diagnosed structural congenital anomalies (SCAs) by whole-exome sequencing (WES). BMC Med. Genomics11 (1), 93. 10.1186/s12920-018-0409-z

  • 29

    LippaN. C.BaruaS.AggarwalV.PereiraE.BainJ. M. (2021). A novel de novo KDM5C variant in a female with global developmental delay and ataxia: a case report. BMC Neurol.21 (1), 358. 10.1186/s12883-021-02380-9

  • 30

    LiuY.WangL.YangY. K.LiangY.ZhangT. J.LiangN.et al (2019). Prenatal diagnosis of fetal skeletal dysplasia using targeted next-generation sequencing: an analysis of 30 cases. Diagn Pathol.14 (1), 76. 10.1186/s13000-019-0853-x

  • 31

    MacedoM. P.CarvalhoF.GuimarãesS.RebeloS.RamalhoC. (2022). Early gestational diagnosis of lethal skeletal dysplasias: a 15 Year retrospective cohort reviewing concordance between ultrasonographic, genetic and morphological features. Fetal Pediatr. Pathol.41 (1), 3748. 10.1080/15513815.2020.1761915

  • 32

    ManickamK.McClainM. R.DemmerL. A.BiswasS.KearneyH. M.MalinowskiJ.et al (2021). Exome and genome sequencing for pediatric patients with congenital anomalies or intellectual disability: an evidence-based clinical guideline of the American College of Medical Genetics and Genomics (ACMG). Genet. Med.23 (11), 20292037. 10.1038/s41436-021-01242-6

  • 33

    ShiY.MaJ.XueY.WangJ.YuB.WangT. (2019). The assessment of combined karyotype analysis and chromosomal microarray in pregnant women of advanced maternal age: a multicenter study. Ann. Transl. Med.7 (14), 318. 10.21037/atm.2019.06.63

  • 34

    SparksT. N.LianoglouB. R.AdamiR. R.PluymI. D.HollimanK.DuffyJ.et al (2020). Exome sequencing for prenatal diagnosis in nonimmune hydrops fetalis. N. Engl. J. Med.383 (18), 17461756. 10.1056/NEJMoa2023643

  • 35

    SymondsJ. D.McTagueA. (2020). Corrigendum to "Epilepsy and developmental disorders: next generation sequencing in the clinic" [Eur. J. Paediatr. Neurol. (2019) 15-23]. Eur. J. Paediatr. Neurol.24, 170223. 10.1016/j.ejpn.2020.06.011

  • 36

    Van der AuweraG. A.CarneiroM. O.HartlC.PoplinR.Del AngelG.Levy-MoonshineA.et al (2013). From FastQ data to high confidence variant calls: the Genome Analysis Toolkit best practices pipeline. Curr. Protoc. Bioinforma.43 (1110), 11. 10.1002/0471250953.bi1110s43

  • 37

    van HengelJ.CaloreM.BauceB.DazzoE.MazzottiE.De BortoliM.et al (2013). Mutations in the area composita protein αT-catenin are associated with arrhythmogenic right ventricular cardiomyopathy. Eur. Heart J.34 (3), 201210. 10.1093/eurheartj/ehs373

  • 38

    VoraN. L.GilmoreK.BrandtA.GustafsonC.StrandeN.RamkissoonL.et al (2020). An approach to integrating exome sequencing for fetal structural anomalies into clinical practice. Genet. Med.22 (5), 954961. 10.1038/s41436-020-0750-4

  • 39

    VoraN. L.PowellB.BrandtA.StrandeN.HardistyE.GilmoreK.et al (2017). Prenatal exome sequencing in anomalous fetuses: new opportunities and challenges. Genet. Med.19 (11), 12071216. 10.1038/gim.2017.33

  • 40

    VrijenhoekT.MiddelburgE. M.MonroeG. R.van GassenK. L. I.GeenenJ. W.HövelsA. M.et al (2018). Whole-exome sequencing in intellectual disability; cost before and after a diagnosis. Eur. J. Hum. Genet.26 (11), 15661571. 10.1038/s41431-018-0203-6

  • 41

    WapnerR. J.MartinC. L.LevyB.BallifB. C.EngC. M.ZacharyJ. M.et al (2012). Chromosomal microarray versus karyotyping for prenatal diagnosis. N. Engl. J. Med.367 (23), 21752184. 10.1056/NEJMoa1203382

  • 42

    YangK.ShenM.YanY.TanY.ZhangJ.WuJ.et al (2019). Genetic analysis in fetal skeletal dysplasias by trio whole-exome sequencing. Biomed. Res. Int.2019, 2492590. 10.1155/2019/2492590

  • 43

    YangY.WangM.WangH. (2022). Prenatal trio-based whole exome sequencing in fetuses with abnormalities of the skeletal system. Mol. Genet. Genomics297 (4), 10171026. 10.1007/s00438-022-01899-x

  • 44

    ZhangL.PanL.TengY.LiangD.LiZ.WuL. (2021). Molecular diagnosis for 55 fetuses with skeletal dysplasias by whole-exome sequencing: a retrospective cohort study. Clin. Genet.100 (2), 219226. 10.1111/cge.13976

Summary

Keywords

prenatal diagnosis, ultrasonographic anomalies, WES, diagnostic yield, pregnancy outcomes

Citation

Wang Y, Yin F, Chai Y, Jin J, Zhang P, Tan Q and Chen Z (2024) Prenatal diagnosis of fetuses with ultrasound anomalies by whole-exome sequencing in Luoyang city, China. Front. Genet. 14:1301439. doi: 10.3389/fgene.2023.1301439

Received

25 September 2023

Accepted

20 December 2023

Published

22 January 2024

Volume

14 - 2023

Edited by

José María Frade, Spanish National Research Council (CSIC), Spain

Reviewed by

Sateesh Maddirevula, King Faisal Specialist Hospital and Research Centre, Saudi Arabia

Tonia Carter, Marshfield Clinic Research Institute, United States

Updates

Copyright

*Correspondence: Yanan Wang,

†These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics