Efficiency of copy number variation sequencing combined with karyotyping in fetuses with congenital heart disease and the following outcomes

Abstract

Background

Both copy number variant-sequencing (CNV-seq) and karyotype analysis have been used as powerful tools in the genetic aetiology of fetuses with congenital heart diseases (CHD). However, CNV-seq brings clinicians more confusions to interpret the detection results related to CHD with or without extracardiac abnormalities. Hence, we conducted this study to investigate the clinical value of CNV-seq in fetuses with CHD.

Results

A total of 167 patients with fetal CHD including 36 single CHD (sCHD), 41 compound CHD (cCHD) and 90 non-isolated CHD (niCHD) were recruited into the study. 28 cases (16.77%, 28/167) were revealed with chromosomal abnormalities at the level of karyotype. The pathogenic detection rate (DR) of CNV-seq (23.17%, 19/82) was higher than that of karyotyping (15.85%, 13/82) in 82 cases by CNV-seq and karyotyping simultaneously. The DR of pathogenic copy number variations (PCNVs) (31.43%) was higher in niCHD subgroup than that in sCHD and cCHD (9.52% and 23.08%). Conotruncal defect (CTD) was one of the most common heart malformations with the highest DR of PCNVs (50%) in 7 categories of CHD. In terms of all the pregnancy outcomes, 67 (40.12%) cases were terminated and 100 (59.88%) cases were live neonates. Only two among 34 cases with a pathogenic genetic result chose to continue the pregnancy.

Conclusions

CNV-seq combined with karyotyping is a reliable and accurate prenatal technique for identifying pathogenic chromosomal abnormalities associated with fetal CHD with or without extracardiac abnormalities, which can assist clinicians to perform detailed genetic counselling with regard to the etiology and related outcomes of CHD.

Background

Congenital heart disease (CHD) is the most common birth defect, occurring in approximately 0.4–5% of live births [1, 2], and up to 10% of stillbirths [3]. Genetic abnormalities, including aneuploidies (AUP), copy number variants (CNVs) and single nucleotide variants, play a significant role in determining the clinical outcome of CHD [4, 5]. With surgical development, most types of CHD return to normal or near normal after cardiovascular surgery. Once combined with chromosomal abnormalities, the prognosis of fetuses with CHD is poor due to common complications such as severe extracardiac structural abnormalities, mental retardation and developmental delay, etc. Therefore, prenatal genetic diagnosis is strongly recommended for fetuses with CHD. Compared to traditional G-banding karyotyping, CNV-seq can detect additional CNVs of more than 100 kb [6, 7]. Both of the complementary techniques may improve the detection rate (DR) of chromosomal abnormalities in CHD fetuses. However, the drawback of present condition is that CNV-seq detection brings more results for clinicians to interpret the detection results related with the fetal manifestations. Hence, it would of difficulty to get a deeper insight to understand the genetic abnormalities associated with fetuses with CHD through CNV-seq combined with karyotyping. Herein, our group conducted a retrospective study to evaluate the clinical value in terms of the genetic etiology by CNV-seq combined with karyotyping and ultimate pregnancy outcomes of fetuses with CHD. In addition, we also stratified 82 cases by simultaneous CNV-seq and karyotype detection to better understand the DR of chromosomal abnormalities in different types of CHD, and compared the frequency of chromosomal abnormalities in fetuses with single CHD (sCHD), compound CHD (cCHD) and non-isolated CHD (niCHD).

Results

Basic characteristics of study subjects

From January 2016 to October 2022, a total of 167 pregnant women with fetal CHD including 36 sCHD, 41 cCHD and 90 niCHD were recruited into the study. Although there were no significant differences in maternal age (MA), gestational age (GA) and parity history (P > 0.05) among sCHD, cCHD, and niCHD groups, a significant difference of pregnancy outcome (P < 0.01) suggested that the fetuses with niCHD were more likely to be terminated. In the end, 100 cases (59.88%) chose to continue their pregnancy and 67 cases (40.12%) chose to terminate their pregnancy (Table 1).

Table 1 Basic characteristics of study subjects

Karyotype and CNV-seq results in CHD cases

21 AUP and 7 unbalanced chromosomal structural abnormalities were found (Table 2) in 167 fetuses. 21 AUP were consisted of 9 trisomy 21, 7 trisomy 18, 1 trisomy 13, 1 partial trisomy 14 and 9 (47,XY, + 14,der(14) (9;14)(9p23;q22)mat), 2 monosomy X (containing 1 mosaic monosomy X) and 1 47, XXY, demonstrating that trisomy 21 and 18 accounted for the major proportion. The DR of chromosomal abnormalities in the niCHD subgroup (23.33%, 21/90) were obviously higher than that in sCHD group (2.78%, 1/36) (P = 0.013).

Table 2 Chromosomal aberration in 167 fetal CHD detected by karyotyping

In addition, in 82 cases of copy number variation detection combined with karyotype analysis, the pathogenic DR of CNV-seq (23.17%, 19/82) was higher than that of karyotype analysis (15.85%, 13/82) (P = 0.237) (Table 3). The 31.43% (11/35) DR of pathogenic copy number variations (PCNVs) in the niCHD group was higher than 9.52% (2/21) in the sCHD and 23.08% (6/26) in the cCHD, with no significant difference in the DR of chromosomal aberrations (P = 0.156) or in the PCNVs among the three groups (P = 0.142) (Table 3). Moreover, CHD cases were classified to 7 categories according to international statistical classification of diseases and related health problems 10th revision (ICD-10), and the results showed that the DR of PCNVs (50%, 7/14) in conotruncal defect (CTD) is highest, followed by atrioventricular septal defect (AVSD) (Additional file 1: Table S1). Thirty CNVs of 82 cases [19 PCNVs, 6 benign (B) CNVs, 2 like benign (LB) CNVs, and 3 variants of uncertain significance (VOUS)] were revealed by CNV-seq, including several syndromes such as Down's syndrome (n = 4), Edward's syndrome (n = 2), 22q11.2 microdeletion syndrome (n = 1), 1p36 microdeletion syndrome (n = 1), Wolf-Hirschhorn syndrome (n = 1), and recurrent Simpson-Golabi-Behmel syndrome type 1 (SGBS1) (n = 2) (Table 4), plus extra CHD anomalies indicating related syndromes should be considered by clinicians once facing multiple fetal structural abnormalities.

Table 3 Abnormal results in 82 CHD detected by karyotyping and CNV-seq

Table 4 Phenotype and genotype list of CNVs

Among 35 niCHD cases, CHD with cardiovascular structural anomalies and cardiovascular soft index anomalies accounted for 60% (21/35) and 40% (14/35), respectively. Although the DR of PCNVs in the former group (38.1%, 8/21) was higher than that in the latter group (21.43%, 3/14), there were no statistical significance between two groups (P = 0.623) (Additional file 2: Table S2).

Pregnancy outcomes and prognosis classification of all cases

Compared to fetuses without CNVs or with B/LB/VOUS CNVs, more couples chose to induce labor regarding to their fetuses with AUP or P/LP CNVs (94.12% (32/34) vs. 26.32% (35/133), P < 0.001) (Additional file 3: Table S3). 133 fetuses without pathological chromosomal alterations were classified into four grades according to the prognostic grading of fetal CHD (Additional file 4: Table S4) [8], and we found fetuses with CHD of grade III and IV were all terminated. Overall, thirty-five (26.32%, 35/133) were terminated and ninty-eight (73.68%, 98/133) were live neonates in the following day. The fetuses with chromosomal anomalies who were terminated included VOUS [n = 3, one case with a 1.58-Mb 8p21.3p21.2 duplication presented with AVSD, one case with a 0.24-Mb X p21.22 duplication manifested right aortic arch (R-AA), mirror branch of brachiocephalic artery (BA) and hydrocephalus, and one case with a 10.64-Mb 4q13.1 duplication showed single ventricle (SV), persistent truncus arteriosus (PTA) and single umbilical artery], BCNVs [n = 2, one case with a 0.44-Mb 14q11.2 deletion presented with single atrium, mild mitral regurgitation (MR), persistent left superior vena cava (PLSVC), and one case with a 0.36-Mb 11p14.3 deletion showed tetralogy of fallot (TOF), R-AA, left aberrant subclavian artery (L-ASA)]. Other 30 cases chose to terminate the pregnancy due to poor prognosis of severe cardiac defects (with grade 3 and 4) or severe extracardiac defects (with grade 1 and 2).

In the six-month follow-up studies of 100 newborns, including 7 cases with prenatal ventricular/atrial septal defect (VSD/ASD) showed normal cardiovascular structure on ultrasound, 5 cases showed new ASD or more severe VSD than before, 77 newborns were consistent with prenatal cardiovascular structure on ultrasound, 11 newborns without further postpartum cardiac ultrasound examination (two fetuses with pathogenic results were included, and presented normal). Hence, dynamic ultrasound evaluation and follow-up after birth could provide practical guidance to clinicians.

Discussion

Previous study has demonstrated that the majority of CHD and other extracardiac defects may result from genetic factors, among which AUP accounts for 9–18% of CHD [4, 9, 10]. In our present paper, 28 cases (16.77%) of chromosomal aberrations were detected by karyotyping in 167 fetuses with CHD. The DR of aneuploidies accounted for 12.57% (21/167) in our study, among which trisomy 21 and 18 presented with the largest proportion (76.19%, 16/21), similar results to those in the literature [11,12,13]. The above information indicates that the dose alterations of genetic materials might increases the risk of malformations of fetuses [12]. Hence, in clinical practice, dynamic ultrasound should be supervised to evaluate systematic and comprehensive structural investigation for the early revealed imbalanced chromosome number or structural alterations. Moreover, 3% to 25% of fetuses with CHD have been reported to be associated with PCNVs [9, 14, 15]. As shown in Table 3, the pathogenic DR of CNV-seq was 23.17% (19/82) similar to that observed in the previous study [13], containing severe syndromes, such as SGBS1, 22q11.2 microdeletion. Relative to traditional karyotyping, an additional 7.32% PCNVs (6/82) by CNV-seq in our study was similar to previous research 7.70 ~ 7.95% by chromosome microarray analysis (CMA) [12, 16]. Hence, CNV-seq could provide an efficient and equivalent to CMA and recognize the microdeletion or microduplication of chromosomes (MMS). Besides above advantages, CNV-seq based on next-generation sequencing is emerging as an alternative methodology due to needed smaller sample size of test, faster experimental cycle and lower cost [11, 17, 18].

It is worth noting that the majority of the CHD fetuses with the MMS were characterized by multiple-system structural anomalies (Table 4) [19]. As can been seen from the basic characteristics of our subjects, fetuses with CHD had no association with MA, GA and parity history (Table 1). In line with the previous report, MMS has been confirmed to be unrelated to the age of the pregnant women [20] Therefore, further comprehensive genetic assessment is required for structural abnormalities of fetuses [21, 22]. Therefore, prenatal diagnosis is strongly recommended for fetuses with CHD whether accompanied with other systematic structure malformations. It is vital for clinicians to comprehensively manage fetuses with CHD based on chromosomal abnormalities, phenotypes in ultrasonic and clinicians’ recommendation.

In 7 subgroups of fetuses with CHD, the statistical analysis demonstrated the major frequencies of chromosomal abnormalities occurred in fetuses with CTD (7/14, 50%), followed by AVSD (1/3, 33.33%), septal defects (6/23, 26.09%) (Additional file 1: Table S1). These DR of chromosomal abnormalities in 7 subgroups were not consistent with the previous reported (73.7% in AVSD, 25.7% in CTD, 17.5% in septal defects) [12], as might be related to the sample sizes and proportions of different types of CHD. The different degrees of DR of chromosomal abnormalities were found in all 7 group, suggesting that fetuses with CHD during pregnancy require chromosome diagnosis. Moreover, 26 fetuses with cCHD presented with 23.08% (6/26) PCNVs including related to 3 kinds of syndromes, significantly higher than those with sCHD 9.52% (2/21) (Table 3), indicating pathogenic CNVs associated syndromes should be considered for compound CHD. In addition, the DR of PCNVs was significantly higher in niCHD with structural anomalies (38.1%, 8/21) than that in niCHD with soft index anomalies (21.43%, 3/14) (Additional file 2: Table S2), consistent with previous studies indicating that CHD fetuses with structural abnormalities are more likely to be related to genetic disorders [23, 24]. Hence, more attentions should be paid to fetuses with multiple structural abnormalities besides CHD.

Once a fetus is diagnosed with CHD, dynamic comprehensive ultrasound, echocardiograms and genetic assessments should be performed properly. Only in this way, could the clinician offer more useful counseling and help for the cases. In our study, there were only 2 normal deliveries in 34 cases with pathogenic genetic results, one fetus with 47, XXY, the other with a 1.34-Mb duplication on chromosome 17p12 inherited from his mother, showing that genetic factors have a significant impact on outcomes of fetuses with CHD. CHD fetuses with chromosomal abnormalities may have a significantly increasing risk of mortality, and require more medical care and medication. Of the remaining 133 cases without pathogenetic chromosomes, 35 cases underwent labor inductions, including 2 BCNVs (1 case with single atrium, mild MR, PLSVC; 1 case with TOF, R-AA; L-ASA), 3 VOUS (1 case with AVSD, 1 case with R-AA, Mirror branch of BA and hydrocephalus, 1 case with SV, PTA and single umbilical artery), suggesting that genetic factors as one aspect to merely be considered and CHD need to be comprehensively evaluated in conjunction with clinical phenotypes (Additional file 4: Table S4). After searching relevant literatures, we have not found the same mutation as VUS fetuses in our report. More data need to be confirmed the relationship between phenotype and VOUS CNVs. For pregnant women who continued their pregnancy, the usual obstetric follow-up of maternal health and the specialized follow-up of fetal health (regular morphological ultrasound) should be ensured. Hence, the destiny of fetuses with CHD often depends on the chromosomal aberrations, the severity of cardiac/extracardiac defects and cognitive level of the parents. It is worth noting that two similar boy fetuses with CHD diagnosed as SGBS1 recurred in a family, and CNV-seq verified it from the mother. Based on SGBS1 as an X-linked recessive disorder, preimplantation genetic diagnosis (PGD) was recommended for the next pregnancy [25, 26]. Therefore, it is helpful to clear genetic history and contribute to the reasonable guidance of second pregnancy.

Although CNV-seq is a very valuable detection tool for prenatal diagnosis, it also has limitations. It cannot effectively measure chromosomal inversion, balanced translocation or some special triploids etc., however, karyotype analysis can efficiently complement the above results. Therefore, the combined application of the two technologies can effectively increase the DR of chromosome abnormalities. VOUS is an additional finding of CNV-Seq to detect fetal CHD, and the application of CNV-Seq technology in prenatal diagnosis of fetal CHD brings the biggest clinical difficulty, which is the interpretation of VOUS [27].

Conclusion

The occurrence of CHD is related to chromosomal AUP and CNVs. CNV-seq is an effective adjunct to traditional chromosomal karyotyping. Genetic analysis combined with dynamic ultrasound screening and multidisciplinary counselling can effectively provide valuable information to the clinicians and patients.

Methods

Study subjects

This retrospective study was performed in the prenatal diagnosis center of Xuzhou Central Hospital of Jiangsu Province from January 2016 to October 2022. 167 fetuses with CHD by echocardiogram were enrolled in our study. Written informed consents for data collection and manuscript publication were provided by all the couples. The demographic characteristics were recorded through a comprehensive questionnaire including maternal age (MA), gestational week (GA), pregnancy outcomes, gravidity, parity, histories of abnormal pregnancy, as well as hereditary diseases in Table 1. The average MA was 28.75 ± 4.27 years, the mean GA was 24.05 ± 2.93 weeks; nulliparas accounted for 45.51% and multiparas accounted for 54.49%. Of 167 fetuses with CHD, isolated CHD (n = 77) included sCHD (n = 36) and cCHD (n = 41), and niCHD (n = 90) had additional extracardiac defects according to the classification of fetal CHD. All pregnant women accepted invasive prenatal diagnosis (127 amniocenteses and 40 cordocentesis) for karyotyping, among whom 82 underwent CNV-seq simultaneously. And all subjects denied hereditary diseases. The proposal for our study has been approved by the ethics committee of Xuzhou Central Hospital (XZXY-LK-20230314-034).

Karyotype analysis

The obtained amniotic fluid and umbilical cord blood were cultured and stained with G-banding according to the general operating procedures. 20 well-dispersed, medium-length metaphase phases from each specimen were analyzed according to ISCN (2016,2020) criteria [28, 29], when mosaicism or abnormal karyotypes were found, this was increased to 100 metaphase phases.

CNV-seq

Genomic DNA (gDNA) was extracted from amniotic fluid or fetal cord blood. Then 50 ng of DNA was fragmented and DNA libraries were constructed by end repair, ligated with sequencing adaptors, polymerase chain reaction (PCR) amplifification, and DNA libraries were subjected to massively parallel sequencing using NextSeq 500 platform (Illumina, San Diego, USA), to generate approximately 5 million raw sequencing reads with genomic DNA sequences of 36 base pair in length. Using the hg19 genomic sequence as reference, a total of 2.8–3.2 million reads were uniquely and precisely mapped using the BurrowseWheeler algorithm [30]. Mapped reads were allocated progressively to 20- kilobase (kb) bin sizes from the p to q arms of the 24 chromosomes. And the criteria of sequencing copy number (SCN) results refered to the attached Additional file 5: Table S5. Several public databases such as Database of Genomic Variants (DGV), Online Mendelian Inheritance in Man (OMIM), DECIPHER, University of California, etc., were utilized to interpret the results as gains and losses of copy number. CNVs were interpreted and divided into five categories: pathogenic(P), likely pathogenic(LP), VOUS, likely benign(LB), and benign(B), according to the guidelines outlined by the American College of Medical Genetics (ACMG) [7].

Statistical analysis

The data were analyzed by SPSS software (version 26.0, IBM, Armonk, NY, USA). Continuous variables were expressed as mean ± standard deviation (SD). Enumeration data were expressed as frequency, and chi-square test or Fisher's exact test was used to compare the rates between groups (P < 0.05) was regarded as statistically significant).