The prevalence of genetic diagnoses in fetuses with severe congenital heart defects.

PURPOSE: Congenital heart defects (CHD) are associated with genetic syndromes. Rapid aneuploidy testing and chromosome microarray analysis (CMA) are standard care in fetal CHD. Many genetic syndromes remain undetected with these tests. This cohort study aims to estimate the frequency of causal genetic variants, in particular structural chromosome abnormalities and sequence variants, in fetuses with severe CHD at mid-gestation, to aid prenatal counselling.
METHODS: Fetuses with severe CHD were extracted from the PRECOR registry (2012-2016). We evaluated pre- and postnatal genetic testing results retrospectively to estimate the frequency of genetic diagnoses in general, as well as for specific CHDs.
RESULTS: 919 fetuses with severe CHD were identified. After exclusion of 211 cases with aneuploidy, a genetic diagnosis was found in 15.7% (111/708). These comprised copy number variants in 9.9% (70/708). In 4.5% (41/708) sequence variants were found that would have remained undetected with CMA. Interrupted aortic arch, pulmonary atresia with ventricular septal defect and atrioventricular septal defect were most commonly associated with a genetic diagnosis.
CONCLUSION: In case of normal CMA results, parents should be offered exome sequencing sequentially, if time allows for it, especially if the CHD is accompanied by other structural malformations due to the large variety in genetic syndromes.

Introduction

Congenital heart defects (CHDs) are the most common congenital birth defects with a prevalence of 5–8 per 1.000 liveborns.[1] Approximately a third of these cases involve severe CHDs.[2] Although survival rates have increased significantly over the past decades, the extent of (disease-related) morbidity is considerable and life-expectance remains reduced.[3]

It is well known that aneuploidies (numerical chromosome abnormalities), such as trisomy 13, 18, 21 and Turner syndrome, are associated with CHDs.[4,5] Submicroscopic deletions or duplications, commonly referred to as copy number variants (CNVs), have been reported in approximately 10–15% of children with CHD.[6,7] As these CNVs can have a major impact on development, quality of life and life expectancy of children with CHDs,[7,8] chromosome microarray analysis (CMA) is the standard of care in the fetal setting in most developed countries.[9] Nevertheless, it is not uncommon in cases with CHD, that a genetic syndrome is diagnosed after birth or during childhood, as sequence variants remain undetected with CMA.

Although CHDs have been described as part of many genetic syndromes, there are only a few cohort studies available that report on genetic diagnoses (chromosome abnormalities or sequence variants) in fetuses diagnosed with a CHD. The studies that have been performed in fetal cohorts are outdated, as they mainly focus on aneuploidy.[10-13] This restricts prenatal counseling, because currently available data are mostly based on postnatal studies. The prevalence of genetic diagnoses is expected to be lower in the postnatal cohorts, as cases with termination of pregnancy (TOP), intra uterine fetal demise or early neonatal death are often not included in these cohorts.[5] More importantly postnatal sequencing is mostly requested due to the evolving clinical phenotype, which makes it difficult to estimate a true prenatal prevalence of genetic syndromes at mid-gestation.

Antenatal knowledge on the prevalence of genetic variants, which may either be considered the genetic cause for the CHD or a secondary finding, can improve prenatal counseling in cases with severe CHDs. Accurate figures on the occurrence of these genetic diagnoses will aid professionals to make decisions regarding the use of currently available additional tests, such as exome sequencing (ES) or genetic testing for specific genetic syndromes.

This study aims to provide a conservative estimate of the prevalence of chromosome abnormalities and sequence variants in fetuses with severe CHDs, in particular after an aneuploidy is excluded, by assessing results from both pre- and postnatal genetic testing. The potential diagnostic yield of ES for fetal CHDs, and factors that potentially increase the chance of an underlying genetic diagnosis, will be evaluated as well.

Materials and Methods

In the northwestern region of the Netherlands, the care for children with CHDs is centralized in three tertiary referral centers in Amsterdam and Leiden. All fetuses and infants diagnosed with a severe CHD within this region have been registered in the PRECOR registry since 2002. Severe CHD was defined as the need for surgery or therapeutic intervention in the first year of life. Data collection for this registry has been previously described.[14] The Leiden University Medical Center has a general privacy statement informing patients that their data can be used for (retrospective) scientific research.

From this registry, all cases with a diagnosis of a severe CHD, born between 2012 and 2016, were extracted. We chose this period, as CMA became a routine diagnostic test from 2012 onwards, if parents opted for invasive testing in pregnancy. In order to ensure a reasonable follow-up period after birth, we included cases born before 2017 and followed them until September 2019. From a clinical perspective, we chose not to exclude cases in which genetic testing was not performed pre- or postnatally to avoid a substantial selection bias, and consequently an overestimation of prevalence of genetic diagnosis. All numbers reported in this cohort are therefore conservative estimates of the prevalence of chromosome abnormalities and sequence variants. Cases with an aneuploidy were not included in subsequent analyses, as the prevalence of aneuploidy in fetal CHD cases has already been well-documented in the literature.[12] To estimate the minimum prevalence of chromosome abnormalities and sequence variants in fetuses with CHDs, results from pre- and postnatal genetic testing, postnatal clinical evaluation and postmortem reports were assessed.

Antenatal characteristics

The obstetric databases were evaluated to collect information on maternal and fetal characteristics, including maternal age, obstetric history, obesity, smoking, alcohol use, positive family history of CHDs, self-reported consanguinity, multiple pregnancy, the presence of additional fetal extra-cardiac malformations (ECM) and prenatal genetic testing. Routine prenatal genetic testing, which involves Quantitative Fluorescence-Polymerase Chain Reaction (QF-PCR) analysis and chromosome microarray analysis (CMA), was offered in all cases with a prenatal diagnosis of a CHD. Genetic testing for specific genetic diagnoses, not detectable with routine prenatal genetic testing, was only performed in fetuses with features that raised suspicion of a specific syndrome. The CHD was classified as either isolated or non-isolated, as it can be expected that the presence of additional structural malformations in the fetus would increase the probability of an underlying genetic cause. Isolated was defined as the absence of significant structural ECMs or intra-uterine growth restriction (according to consensus-based definition[15]) prenatally, whereas non-isolated comprised cases diagnosed with one or multiple structural ECMs. Minor antenatal extra-cardiac abnormalities, also referred to as ‘soft markers’, such as an echogenic bowel or single umbilical artery, were considered non-significant ECMs.

Postnatal course

Data on pregnancy and neonatal outcome, such as gestational age at birth, birth weight, presence of dysmorphic features postnatally, mortality and need for medication were retrieved from electronic patient records. Until the end of the study period we assessed whether signs of developmental delay, such as neurocognitive or speech and language disorders, were reported in the patient records. Clinical records were evaluated for results of CMA, genetic testing for specific syndromes (e.g. if CHARGE, Noonan, Kabuki or Alagille syndrome was suspected) and focused or exome-wide analysis of ES, if parents chose to delay (additional) genetic testing until after birth. If additional signs for a genetic syndrome, such as dysmorphic features or developmental delay, were apparent but less specific and CMA results were normal, ES was considered. ES was only performed in postnatal cases, as this was not yet available in a prenatal setting during the study period. The specific ES based CHD gene panel used in our facilities for focused analysis of sequence data, consisted of a panel of 129 genes that are known to be associated with CHD.[16] A clinical geneticist was consulted in both pre- and postnatal cases before genetic testing for a specific syndrome or exome sequencing was performed.

Classification of all genetic variants

Structural chromosome abnormalities and sequence variants were classified using the guidelines for interpretation of CNV and sequence variants, developed by The American College of Medical Genetics and Genomics (ACMG).[17,18] This classification allocates abnormal results into the following five categories, based on their expected clinical relevance: pathogenic, likely pathogenic, uncertain significance, likely benign and benign. Likely pathogenic was only used when there was a certainty of at least 90% that a variant was disease-causing.[17] Pathogenic and likely pathogenic variants, later referred to as ‘genetic diagnoses’, will therefore be reported together in the subsequent analysis.

The Affymetrix Cytoscan HD array or Agilent CGH 180 K oligo array (Amadid 023363) platforms were used for CMA, as described earlier by Jansen et al.[19] ES was performed using the Agilent SureSelectXT Human all Exon v5 or Clinical Research Exome v2 capture library kit (Agilent, Santa Clara, USA) accompanied by paired end Sequencing on an Illumina sequencing platform (Illumina, San Diego, USA), generating reads with at least 80x median coverage. The sequence analysis pipeline and tool for annotation of variants has extensively been reported on before.[20] Trio samples of the fetus and both parents were assessed for sequencing to optimize variant filtering, when available. Reported variants were submitted to the DECIPHER database.

Statistical analysis

A conservative estimate of the prevalence of chromosomal abnormalities and sequence variants in mid-gestation fetuses diagnosed with a CHD was determined by complementing genetic diagnoses made in pregnancy with those detected after birth, based on postnatal clinical evaluation and follow-up assessment. We assessed the proportion of sequence variants, not detectable with CMA, as this is currently recommended for all fetuses with CHDs.[9] Information on the presence of clinical features for a genetic syndrome is limited on prenatal ultrasound. This reduces the applicability of genetic testing for a specific syndrome in utero. We therefore evaluated the potential diagnostic yield of ES, either analyzed using our specific CHD gene panel[16] or exome-wide, for the detection of (likely) pathogenic sequence variants in fetal CHD cases. The probability of structural chromosome abnormalities and sequence variants, after aneuploidy is excluded, was determined separately for each heart defect that encompassed at least 1% of this cohort. We also calculated whether the probability of genetic diagnoses in fetuses with CHDs is affected by maternal age, obstetric history, family history of CHDs, consanguinity, multiple pregnancy and additional structural fetal malformations. The clinical impact of structural chromosome abnormalities and sequence variants was assessed by comparing neonatal outcome of fetuses with and without abnormal genetic testing results.

Numeric variables were studied for significant differences using an independent t-test, whereas a χ2-test was used to test associations between categorical variables. A Fisher’s exact test was used, if the expected number was < 5. A p-value <0.05 was considered statistically significant. All statistical analyses were performed using IBM SPSS Statistics 25.0 (IBM, Armonk, NY, USA). This study was approved by the Leiden University Medical Ethics Committee.

Results

A total of 919 cases, diagnosed with a severe CHD between January 2012 and December 2016, were extracted from the PRECOR registry. Parents chose invasive prenatal testing in 542/919 cases (59.0%), whereas in 185/919 cases (20.1%) genetic testing was performed after birth. Aneuploidy testing was performed in all these cases. Genetic testing was not performed pre- or postnatally in 192/919 cases (20.9%). In these cases parents either declined genetic testing or it was not indicated, because dysmorphic features were absent and the CHD was expected not to be associated with genetic syndromes. After clinical assessment by a geneticist, signs of a genetic syndrome were absent after birth in the majority of these cases (173/192, 90.1%), whereas 9.9% (19/192) did have additional ECMs (Fig. 1).

Fig. 1

Flow chart case selection.

*Adjusted prevalence after exclusion of aneuploidy cases. CHD congenital heart defect, CNV copy number variant, T21 trisomy 21, T18 trisomy 18, T13 trisomy 13.

An aneuploidy was found in 211/919 (23.0%) cases. As we were mainly interested in the prevalence of structural chromosome abnormalities and sequence variants, rather than aneuploidy, only the remaining 708 euploid cases were included in further analyses. Genetic testing in these cases involved CMA in 64.7% (458/708), genetic testing for specific syndromes in 13.3% (94/708) and focused (7.2%, 51/708) or exome-wide analysis of exome sequencing data (8.1%, 57/708). Baseline characteristics of these cases are enclosed as Supplementary Material (Table S1).

Prevalence genetic variants

An estimate of the prevalence of chromosome abnormalities and sequence variants amongst fetuses with a CHD in this cohort was determined by complementing all genetic variants reported (including ‘uncertain significance’) in pregnancy with those detected after birth. An overview of all genetic variants encountered in this cohort, is shown in Table 1.

Table 1

Abnormal results chromosome microarray analysis & exome sequencing.

Structural chromosome anomalies (CMA)		Sequence variants (ES)
Pathogenic	63(8.9%)	Pathogenic	33(4.7%)
22q11.2 syndrome	30	CHARGE syndrome	5
1q21.1 deletion syndrome	3	Kabuki syndrome	3
16p11.2 deletion syndrome	2	Noonan syndrome	3
Beckwith-Wiedemann syndr.	1	Tuberous sclerosis complex	3
Cri du chat syndrome	1	Alagille syndrome	2
Phelan-McDermid syndrome	1	CM-AVM syndrome	1
Ringchromosome 20 syndr.	1	Cornelia de Lange syndrome	1
Waardenburg syndrome	1	Jeune syndrome	1
Williams syndrome	1	Rubinstein-Taybi syndrome	1
6p25 microdeletion syndr.	1	Schaaf-Yang syndrome	1
8p23.1 microdeletion syndr.	1	Smith-Lemli-Opitz syndrome	1
13q deletion syndrome	1
15q11-13 duplication syndr.	1
18p deletion syndrome	1
Others	17	Others	11
Likely pathogenic	7 (1.0%)	Likely pathogenic	8 (1.1%)
Uncertain significance	19 (2.7%)	Uncertain significance	20 (2.8%)
Likely benign	23 (3.2%)	Likely benign	0 (0.0%)
Benign	3 (0.4%)	Benign	0 (0.0%)
Total (without (likely) benign)	89 (12.6%)	Total (without (likely) benign)	61 (8.6%)

Data are given as n (%).

Syndr. syndrome,CM-AVM Capillary malformation-arteriovenous malformation syndrome.

The American College of Medical Genetics and Genomics’ and the Association for Molecular Pathology’s guideline for interpreting copy number variations and sequence variants was used to categorize genetic variants into the five categories displayed above.

Genetic diagnoses were encountered in 111/708 (15.6%) euploid cases with severe CHD. These comprised copy number variants (CNVs) diagnosed with CMA in 70/708 cases (9.9%), of which 63/70 were classified as pathogenic and 7/70 likely pathogenic. In 2.7% of cases (19/708) a CNV of uncertain significance was found. The 22q11.2 deletion syndrome was the most common microdeletion syndrome encountered, with a prevalence of 4.2% (30/708), encompassing 42.9% (30/70) of all (likely) pathogenic CNVs detected.

In 41/708 (5.8%) fetuses with a CHD, a potentially causal sequence variant was detected with additional testing. These variants were classified pathogenic in 33/41 and likely pathogenic in 8/41. A sequence variant of uncertain significance was found in 2.8% of cases (20/708). Frequently encountered genetic syndromes included CHARGE (5/41, 12.2%), Kabuki (3/41, 7.3%) and Noonan syndrome (3/41, 7.3%).

Diagnostic yield additional testing

In 65.9% (27/41) of (likely) pathogenic sequence variants, genetic testing for a specific syndrome (e.g. Kabuki or Noonan syndrome) had been performed based on clinical suspicion. These variants were diagnosed postnatally in the majority of cases (16/27, 59.3%). The remaining 34.1% (14/41) of (likely) pathogenic sequence variants were detected with either focused or exome-wide analysis of exome sequencing data.

To compare the strengths and limitations of a gene panel approach to an exome-wide analysis of ES data, we analyzed the diagnostic yield in relation to all reported sequence variants (including ‘uncertain’ significance). In 24.6% (15/61) of all sequence variants, a (likely) pathogenic variant was encountered that did not involve genes known to be related with CHDs. Significant variants, such as a RASA1 sequence variant causing capillary malformation-arteriovenous malformation syndrome or a variant in the TSC1 gene leading to tuberous sclerosis (not suspected prenatally), will thus remain undetected if ES data are only analyzed using a targeted CHD gene panel. The disadvantage of exome-wide analysis of ES data is the increased identification of variants of ‘uncertain significance’, as identified in 19.7% (12/61) of all sequence variants.

Genetic cause for CHD

Eighty-five percent (94/111) of the encountered structural chromosome abnormalities and sequence variants were considered a definite explanation for the development of the heart defect. This means that in 13% (94/708) of euploid CHD cases, a genetic cause for the heart defect was found, which was 18% (94/516) for those who underwent genetic testing. Secondary findings involved genetic causes for severe hemophilia, Charcot-Marie-Tooth Disease and psychomotor retardation.

Structural chromosome abnormalities and sequence variants that were classified pathogenic or likely pathogenic, were both considered a pathogenic diagnosis in subsequent analyses, and therefore reported together. The remaining cases without genetic testing results were regarded as having a normal test result, whereas those with a variant of uncertain significance were reported separately.

Specific heart defects

Heart defects that were most frequently accompanied by a genetic diagnosis included an interrupted aortic arch (IAoA), pulmonary atresia with a ventricular septal defect (PA-VSD), (un)balanced atrioventricular septal defect (AVSD), truncus arteriosus and Tetralogy of Fallot. Genetic diagnoses were, on the other hand, barely encountered amongst fetuses with a diagnosis of transposition of the great arteries with intact ventricular septum (simple TGA), tricuspid atresia, double inlet left ventricle or total anomalous pulmonary venous connection. The genetic diagnoses encountered in two cases with simple TGA were both secondary findings that did not explain the heart defect. In one case hemophilia type A was diagnosed, and the mother turned out to be a carrier, whereas the other comprised a de novo distal 16q11.2 deletion which was very unlikely to be associated with the heart defect. For each common CHD, and isolated cases separately, the probability for genetic diagnoses is presented in Table 2.

Table 2

Prevalence structural chromosome abnormalities and sequence variants for specific diagnoses.

	Genetic diagnosis							Uncertain significance
	Yes		No		OR	95% CI	p
Interrupted aortic archb	10	71.4%	4	28.6%	13.7	4.22–44.57	<0.001a	0	0.0%
isolated	8	66.7%	4	33.3%				0	0.0%
Pulmonary atresia with VSD	5	41.7%	7	58.3%	3.7	1.16–11.92	0.03a	0	0.0%
isolated	1	16.7%	5	83.3%				0	0.0%
AVSD, unbalanced	3	33.3%	6	66.7%	2.6	0.63–10.38	0.18	0	0.0%
isolated	1	20.0%	4	80.0%				0	0.0%
AVSD, balanced	8	28.6%	19	67.9%	2.2	0.94–5.17	0.11	1	3.6%
isolated	4	21.1%	14	73.7%				1	5.3%
Isolated right aortic arch	2	28.6%	5	71.4%	2.0	0.39–10.60	0.33	0	0.0%
isolated	1	16.7%	5	83.3%				0	0.0%
Truncus arteriosus	4	25.0%	11	68.8%	1.9	0.58–5.95	0.29	1	6.3%
isolated	4	36.4%	7	63.6%				0	0.0%
Tetralogy of Fallot	12	21.1%	43	75.4%	1.5	0.74–2.85	0.28	2	3.5%
isolated	8	17.0%	38	80.9%				1	2.1%
Valvular aortic stenosis	5	20.0%	20	80.0%	1.3	0.47–3.46	0.59	0	0.0%
isolated	2	10.0%	18	90.0%				0	0.0%
Ventricular septal defect	19	17.8%	84	78.5%	1.2	0.68–2.01	0.58	4	3.7%
isolated	6	8.5%	63	88.7%				2	2.8%
Left isomerism	2	16.7%	10	83.3%	1.0	0.22–4.65	1.00	0	0.0%
isolated	1	25.0%	3	75.0%				0	0.0%
DORV-Fallot type	5	16.7%	23	76.7%	1.1	0.41–2.95	0.80	2	6.7%
isolated	2	12.5%	14	87.5%				0	0.0%
DORV-Taussig Bing	3	13.0%	18	78.3%	0.8	0.24–2.88	1.00	2	8.7%
isolated	2	10.5%	16	84.2%				1	5.3%
Hypoplastic aortic arch	1	12.5%	6	75.0%	0.8	0.10–7.02	1.00	1	12.5%
isolated	1	20.0%	4	66.7%				1	16.7%
Hypoplastic left heart syndrome	6	11.5%	44	84.6%	0.7	0.28–1.61	0.36	2	3.8%
isolated	5	11.6%	38	84.4%				2	4.4%
Aortic coarctation	5	7.7%	55	84.6%	0.4	0.17–1.10	0.07	5	7.7%
isolated	3	5.8%	49	86.0%				5	8.8%
Valvular pulmonary stenosis	2	6.7%	22	73.3%	0.4	0.10–1.93	0.40	6	20.0%
isolated	2	8.7%	21	77.8%				4	14.8%
Hypoplastic right heart syndrome	1	6.7%	12	80.0%	0.4	0.05–3.21	0.71	2	13.3%
isolated	0	0.0%	10	83.3%				2	16.7%
TGA (with significant VSD or PS)	1	5.0%	18	90.0%	0.3	0.04–2.06	0.34	1	5.0%
isolated	1	5.9%	16	94.1%				0	0.0%
TGA (simple)	2	4.3%	44	93.6%	0.2	0.05–0.90	0.02a	1	2.1%
isolated	2	4.4%	43	93.5%				1	2.2%
Tricuspid valve atresia	1	4.0%	23	92.0%	0.2	0.03–1.58	0.16	1	4.0%
isolated	1	4.8%	20	90.9%				1	4.5%
Double inlet left ventricle	0	0.0%	6	85.7%	n/a		0.60	1	14.3%
isolated	0	0.0%	6	85.7%				1	14.3%
TAPVC	0	0.0%	10	90.9%	n/a		0.38	1	9.1%
isolated	0	0.0%	8	88.9%				1	11.1%
Miscellaneous	14	15.9%	68	77.3%				6	6.8%
Total	111	15.7%	558	78.8%				39	5.5%

Data are given as n (%).

aP-value < 0.05 was considered statistically significant.

b8/10 with a pathogenic variant was diagnosed with 22q11 deletion syndrome.

VSD ventricular septal defect, AVSD atrioventricular septal defect, DORV double outlet right ventricle, TGA transposition of the great arteries, VUS variant of uncertain significance, PS pulmonary valve stenosis, TAPVC total anomalous pulmonary vein connection.

Risk factors for a genetic diagnosis

The probability of a genetic diagnosis was significantly lower for prenatally isolated (11.6%, 63/541) compared to non-isolated cases (28.7%, 48/167) (p < 0.001). Self-reported consanguinity (35.5% vs 2.9%, p = 0.002) and a positive parental history of a CHD (32.3% vs 3.8%, p = 0.03) were also significantly associated with genetic diagnoses. Other variables of interest, such as maternal age, gravidity, obesity or multiple pregnancy were not significantly correlated with genetic diagnoses (Table 3).

Table 3

Genetic variants: risk factors & outcome.

	Pathogenic variants							Uncertain significance
	Yes (n = 111)		No (n = 558)		OR	95% CI	p
Non-isolated CHD	48	28.7%	107	19.2%	0.3	0.20–0.48	<0.001a	12	30.8%
Isolated CHD	63	11.6%	451	80.8%				27	69.2%
Dysmorphic features	36	28.8%	73	13.1%	6.9	3.95–12.07	<0.001a	16	41.0%
Age mother (yr)	31.0	5.38	30.9	4.94			0.76
Gravida	2.6	2.10	2.2	1.43			0.07
Consanguinity, self-reported	11	35.5%	16	2.9%	3.7	1.68–8.25	0.002a	4	10.3%
Positive history of CHD, 1st degree	10	32.3%	21	3.8%	2.4	1.09–5.22	0.03a	1	2.6%
Multiple pregnancy	4	8.5%	41	7.3%	0.5	0.16–1.34	0.15	2	5.1%
Pregnancy outcome
Termination of pregnancy	41	36.9%	136	24.4%	1.8	1.18–2.80	0.01a	4	10,3%
Intrauterine fetal demiseb.c	6	5.4%	23	4.1%	1.1	0.30–3.74	1.00	1	2.6%
Live births	64	57.7%	399	71.5%	0.5	0.36–0.83	0.004a	34	87.2%
Gestational age at birthd	38.0	2.83	38.7	2.57			0.03a	37.9	2.63
< 37 weeks	13	20.6%	50	9.0%	1.8	0.91–3.53	0.09	10	25.6%
Birth weight (g)d,e	2870.2	761.82	3165.7	710.09			0.003a	2821.7	689.80
< 3rd centile	13	20.6%	36	6.5%	2.5	1.23–4.98	0.01a	7	17.9%
Mortalityd	21	32.8%	50	9.0%	3.4	1.87–6.21	<0.001a	9	23.1%
<1 year	20	95.2%	48	8.6%				7	17.9%
Developmental delay	31	75.6%	54	9.7%	15.4	7.15–33.37	<0.001a	11	28.2%
Medication usef	12	30.0%	30	5.4%	4.1	1.89–8.89	0.001a	4	10.3%
Therapy modalities for delayg	31	79.5%	53	9.5%	19.0	8.27–43.66	<0.001a	11	28.2%
Follow-up (yr)d	3.4	2.79	4.4	2.28			0.02a	4.0	2.73

All data are given in n (%) or mean [SD].

aP-value < 0.05 was considered statistically significant.

bIncluded 9 cases that deceased during labour (3 with and 6 without genetic diagnosis).

cOR and P-value calculated for continuing pregnancies (TOP cases not included).

dLive births.

e1 gram equals 0.04 ounces.

fOther than heart medication.

ge.g. speech therapy, physical therapy, psychological therapy or special education.

TOP termination of pregnancy.

Impact on outcome

The effect of genetic diagnoses on the pregnancy outcome of all included cases is depicted in Table 3. The postnatal outcome was compared between cases with and without a genetic diagnosis. First of all, parents opted for TOP significantly more often when the heart defect was accompanied by a genetic diagnosis compared to those without a genetic diagnosis (36.9% vs 24.4%; p = 0.01). The detection of variants of ‘uncertain significance’ did not lead to an increase in TOPs, as the proportion of parents that terminated pregnancy was even lower in these cases compared to those with normal genetic testing results (10.3% vs 24.4%; p = 0.04). This indicates that the specific diagnosis rather influences the parental decision for TOP than the detection of variants of ‘uncertain significance’. The proportion of cases with a birth weight < 3rd centile was higher amongst CHD cases with a genetic diagnosis (20.6% vs 6.5%; p = 0.01). Postnatal mortality was also increased in cases with a genetic diagnosis (32.8% vs 9.0%, p < 0.001). Signs of developmental delay were significantly more often present in children with (75.6%) compared to those without a genetic diagnosis (9.7%) (p < 0.001). Cases with genetic diagnoses required medication for other reasons than the cardiac defect itself, more often (30.0% vs 5.4%, p = 0.001), as well as (medical) support, such as speech therapy, physical therapy or special education (79.5% vs 9.5%, p < 0.001).

Prevalence including aneuploidy

The minimum prevalence of (likely) pathogenic CNVs in all fetuses with a CHD appears 7.6% (70/919). For sequence variants, not detectable with standard micro-array testing, the prevalence of (likely) pathogenic variants is estimated at 4.5% (41/919) of all fetuses with a CHD (Fig. 1). A genetic cause for the CHD was encountered in 33% (305/919) of all CHD cases, which was 42% (305/708) for those who underwent genetic testing (Supplementary Material, Figure S1).

Discussion

This is the largest cohort study on the total prevalence of genetic diagnoses in fetuses with a severe CHD that included results from genetic testing for specific syndromes and ES as well. In the setting of prenatal counseling concerning a fetal heart defect, a 15.7% probability should be counseled for clinically significant genetic diagnoses, other than aneuploidy. These involved CNVs in 9.9%, whereas 5.8% had a sequence variant not detectable with QF-PCR and CMA. These numbers, however, comprise conservative estimates, as not all patients underwent genetic testing and exome sequencing was not performed in the majority of cases. More importantly, a genetic diagnosis worsened the prognosis significantly, both on surgical outcome and, not unexpectedly, (neuro-)development.[21]

Our results show that a genetic diagnosis has a substantial impact on neonatal outcome in fetuses with a severe CHD. Not only does it significantly increase the risk of mortality, but also morbidity, as these cases showed more often developmental delay (75.6% vs 9.7%) and required more (medical) support and medication, other than for the heart defect itself. The estimated prevalence of pathogenic CNVs in 7.6% of our entire cohort corresponds to three previously described cohorts of fetuses with CHDs, as these report pathogenic CNVs in 8–11%.[22-24] The proportion of CNVs appeared slightly lower in postnatal cohorts that reported CNVs in 5–8% of neonates with CHDs.[25,26] One recent fetal cohort reported a prevalence of 16% pathogenic CNVs.[27] This proportion, however, reflects a selected population, as they only included those referred for invasive genetic testing.

As CMA is routinely offered in pregnancy, there remains a residual probability of at least 6.6% (41/616) to identify (likely) pathogenic sequence variants after birth. Due to the large variety of genetic syndromes, which can occur in less than 1:1000 CHD cases, exome sequencing should be considered in fetuses with CHDs. Especially, as it may change the prognosis considerably, which is why parents opt for invasive genetic testing in the first place. Our results suggest that the diagnostic yield of ES will be no less than 6.6% if ES is offered routinely for all fetal CHDs, and even higher for those with additional structural malformations. This association with ECMs is in line with previous literature,[14,22,23] and therefore an important factor for prenatal counseling and the decision to perform additional genetic testing in some cases. Another recently published fetal cohort found diagnostic sequence variants in as many as 13.6% of fetuses with any cardiac malformation with normal results after aneuploidy testing and CMA.[28] These findings can again be explained by the selection of their cohort; ES was not offered to all women and only those with ES results were eligible for inclusion. The estimated diagnostic yield encountered in our cohort is therefore conservative.

The prevalence of sequence variants appeared also higher in fetuses (6.6%) than neonates, as one postnatal cohort study reported genetic syndromes in 5.1% of neonates with normal chromosomes.[26] This risk is probably higher in fetal cohorts compared to cohorts that focus exclusively on postnatal cases, as cases with TOP, intra uterine fetal demise or early neonatal death are often not included in postnatal cohorts. It is therefore important that our data are evaluated from a prenatal perspective to enable prenatal counseling at mid-gestation. This is confirmed by the fact that parents in this cohort opted for TOP more often, if a genetic cause was identified, which is similar to our previous findings in prenatal exome sequencing.[20]

This study shows that if a focused approach is chosen as the method of advanced genetic testing in the setting of normal CMA, 36.6% (15/41) of (likely) pathogenic variants would remain undetected. This method is less preferable in a prenatal setting, as essential clinical symptoms may be impossible to detect in the fetus, which may lead to analysis of the ‘wrong’ gene panel. Exome-wide analysis of ES data can detect changes in the entire exome, but at the expense of the turnaround time or costs. Sequence variants in genes not associated with CHDs particularly, may also increase the risk of additional morbidity instead of being an explanation for the development of the heart defect itself. This might be important for prenatal counseling, as it can affect prognosis and neonatal management significantly. Exome-wide analysis of sequence data is also imperative to identify novel pathogenic genes and consequentially improve currently available gene panels for CHDs. It did lead to the detection of variants of uncertain significance in 19.7% (12/61) of all sequence variants, which may complicate prenatal counseling and parental decision making. We believe, however, that the advantages of an exome-wide analysis of exome sequencing data for fetal CHDs may outweigh its difficulties,[20] particularly in the presence of additional structural malformations.

A limitation of this study is that in 20.9% of our entire cohort, genetic testing was not performed. The majority of these cases, however, comprised isolated cases without dysmorphic features after birth that showed a normal development, as all the children do have follow-up visits in our center due to relatively short travel distances and a very low threshold for genetic testing after birth. We therefore chose not to exclude cases, to avoid substantial selection bias, as this might lead to an overestimation of the prevalence.

Due to the large sample size and completeness of our regional CHD registry, we were able to stratify the probability of genetic diagnoses according to the specific heart defect, using not only results from karyotyping and FISH for 22q11.2, but CMA, genetic testing for specific syndromes and exome sequencing in selected cases as well. Although several recent cohorts have studied aneuploidy or 22q11.2 deletion syndrome in CHD cases,[22,23,26,29] evidence on the prevalence of other structural chromosome abnormalities and sequence variants for specific heart defects is limited. IAoA, PA-VSD and AVSD were most associated with the presence of genetic diagnoses. The particularly strong correlation between IAoA and submicroscopic genetic changes (69% probability in IAoA cases) was also demonstrated by a large study that evaluated the results of karyotyping and fluorescence in situ hybridization (FISH) for 22q11.2 in a population of infants with CHDs.[30] Although we evaluated results from additional diagnostic modalities as well, genetic diagnoses were encountered in a similar proportion (71%) of IAoA cases. This stresses that IAoA is mainly associated with 22q11.2 deletion syndrome. One study on CMA in fetuses with a VSD reported pathogenic CNVs in 12%.[13] This might be an overestimation, as our cohort comprised twice as many VSD cases with pathogenic CNVs in 9.3% (10/107). A cohort of Tetralogy of Fallot infants found genetic diagnoses in 25% of cases, after exclusion of aneuploidy, which comprised a heterogenous set of genetic syndromes.[31] This is consistent with the 21.1% probability of pathogenic CNVs or sequence variants in our cohort.

In conclusion, this cohort study shows that, after an aneuploidy is excluded, structural chromosome abnormalities and sequence variants are identified in a substantial proportion of cases with severe CHDs. In 5.8% of euploid fetuses with a CHD, the genetic diagnosis would not have been found, if only CMA had been performed. ES should therefore be considered for fetal CHDs, especially if accompanied by other structural malformations, because genetic diagnoses can affect neonatal outcomes significantly. Future research, which offers ES to all fetuses with a CHD, is however needed to obtain more reliable estimates.

Supplementary Table S1

Supplementary Figure S1