INTRODUCTION {#uog19102-sec-0005}
============

A number of abnormal extracardiac findings have been described in fetuses with congenital heart disease (CHD). They are often small‐for‐gestational age, they have smaller head circumference (HC) and they show signs of abnormal magnetic resonance imaging (MRI), such as smaller brain volume, delayed development and white‐matter injuries, increasing the risk of neurodevelopmental deficiencies[1](#uog19102-bib-0001){ref-type="ref"}, [2](#uog19102-bib-0002){ref-type="ref"}, [3](#uog19102-bib-0003){ref-type="ref"}, [4](#uog19102-bib-0004){ref-type="ref"}, [5](#uog19102-bib-0005){ref-type="ref"}, [6](#uog19102-bib-0006){ref-type="ref"}, [7](#uog19102-bib-0007){ref-type="ref"}, [8](#uog19102-bib-0008){ref-type="ref"}, [9](#uog19102-bib-0009){ref-type="ref"}, [10](#uog19102-bib-0010){ref-type="ref"}, [11](#uog19102-bib-0011){ref-type="ref"}, [12](#uog19102-bib-0012){ref-type="ref"}. Regarding their growth impairment, Rosenthal[12](#uog19102-bib-0012){ref-type="ref"} and Rudolph[13](#uog19102-bib-0013){ref-type="ref"}, [14](#uog19102-bib-0014){ref-type="ref"} proposed two underlying pathophysiological mechanisms: either fetuses with intrinsic growth disturbances are more at risk for developmental errors during cardiogenesis, or CHD might lead to circulatory alterations (such as retrograde flow in the aortic arch), low cerebral blood flow and reduced cerebral oxygen or energy substrate delivery, which are incompatible with optimal growth.

Doppler flow patterns reflecting circulatory alterations in fetuses with CHD have been observed: a lower middle cerebral artery (MCA) pulsatility index (PI), a higher umbilical artery (UA) PI and a lower cerebroplacental ratio (CPR) have been reported in these fetuses compared with healthy ones[1](#uog19102-bib-0001){ref-type="ref"} ^,^ [2](#uog19102-bib-0002){ref-type="ref"} ^,^ [4](#uog19102-bib-0004){ref-type="ref"}, [5](#uog19102-bib-0005){ref-type="ref"}, [6](#uog19102-bib-0006){ref-type="ref"}, [7](#uog19102-bib-0007){ref-type="ref"} ^,^ [15](#uog19102-bib-0015){ref-type="ref"}, [16](#uog19102-bib-0016){ref-type="ref"}, [17](#uog19102-bib-0017){ref-type="ref"}, [18](#uog19102-bib-0018){ref-type="ref"}, [19](#uog19102-bib-0019){ref-type="ref"}, [20](#uog19102-bib-0020){ref-type="ref"}, [21](#uog19102-bib-0021){ref-type="ref"}, [22](#uog19102-bib-0022){ref-type="ref"}, [23](#uog19102-bib-0023){ref-type="ref"}, [24](#uog19102-bib-0024){ref-type="ref"}, [25](#uog19102-bib-0025){ref-type="ref"}. However, the precise relationship between circulatory alterations and fetal growth, particularly of HC and, by implication, brain development, remains unclear. Only a few studies have assessed the relationship between prenatal Doppler flow patterns and fetal growth in CHD fetuses and they were unable to demonstrate a clear association. However, these studies used a single Doppler measurement in a heterogeneous cohort of CHD fetuses[2](#uog19102-bib-0002){ref-type="ref"}, [5](#uog19102-bib-0005){ref-type="ref"}, [7](#uog19102-bib-0007){ref-type="ref"}.

The association between fetal growth and Doppler flow patterns may be influenced by gestational age[1](#uog19102-bib-0001){ref-type="ref"}, [2](#uog19102-bib-0002){ref-type="ref"} and by the nature of the cardiac defect[1](#uog19102-bib-0001){ref-type="ref"}, [5](#uog19102-bib-0005){ref-type="ref"}, [6](#uog19102-bib-0006){ref-type="ref"}, [16](#uog19102-bib-0016){ref-type="ref"}, [19](#uog19102-bib-0019){ref-type="ref"}, [21](#uog19102-bib-0021){ref-type="ref"}. Our aim, therefore, was to define trends in fetal growth and cerebroplacental Doppler flow and to investigate the relationship between fetal head growth and cerebroplacental flow in fetuses with CHD classified according to the expected cerebral arterial oxygen saturation.

METHODS {#uog19102-sec-0006}
=======

Study population {#uog19102-sec-0007}
----------------

This was a retrospective collaborative study in two fetal medicine units in The Netherlands (University Medical Center Groningen and Academic Medical Center Amsterdam). All fetuses with CHD in whom Doppler flow patterns and biometry had been measured serially after 19 weeks\' gestation between January 2010 and November 2016 were included. Fetuses with chromosomal and genetic abnormalities or extracardiac malformations were excluded.

Study design {#uog19102-sec-0008}
------------

As part of routine clinical care, a fetal medicine expert performed fetal biometry and assessed Doppler flow, including HC, abdominal circumference (AC), MCA‐PI and UA‐PI. CPR was calculated as MCA‐PI divided by UA‐PI. All available measurements of HC, AC, MCA‐PI and UA‐PI were retrieved starting from 19 weeks\' gestation, which is the usual referral time after the routine 19--21‐week scan. All fetal measurements were converted into *Z*‐scores to adjust for differences in gestational age based on previously published normative data[26](#uog19102-bib-0026){ref-type="ref"}, [27](#uog19102-bib-0027){ref-type="ref"}, [28](#uog19102-bib-0028){ref-type="ref"}, [29](#uog19102-bib-0029){ref-type="ref"}. In addition, at the first examination following referral, a complete fetal echocardiogram, as well as measurement of the PI of the uterine arteries (UtA‐PI), was performed by a fetal medicine expert and/or an experienced pediatric cardiologist using a standardized protocol to assess fetal cardiac anatomy and function. All cardiac diagnoses were reviewed by a single pediatric cardiologist (S.C.). Postnatal echocardiographic examinations and surgical reports were reviewed for those fetuses which were liveborn and postmortem reports were reviewed in cases of termination of pregnancy or intrauterine fetal demise. Cases with unknown pregnancy outcome were included only if fetal echocardiographic examination was sufficient to confirm the CHD diagnosis. Additional information collected from the maternal and neonatal medical files included gestational age at birth, HC at birth, Apgar score at 5 min and outcome (liveborn, intrauterine fetal demise, termination of pregnancy or neonatal/infant death). We also collected information on maternal complications (e.g. CHD, diabetes, hypertensive disorder, hypothyroidism) and maternal body mass index.

Classification of congenital heart disease {#uog19102-sec-0009}
------------------------------------------

Each fetus was categorized by a pediatric cardiologist into one of three groups according to the expected cerebral arterial oxygen saturation for the particular type of CHD. Much of our knowledge of the fetal circulation is derived from studies performed in fetal lambs[30](#uog19102-bib-0030){ref-type="ref"}. Subsequent human fetal ultrasound studies have shown that blood flow patterns are similar in human and lamb fetuses[13](#uog19102-bib-0013){ref-type="ref"}. Human fetal MRI studies are currently expanding our knowledge of oxygen delivery to the brain and oxygen consumption in normal hearts and in CHD, but precise knowledge of cerebral arterial oxygen saturation in fetal CHD is lacking[11](#uog19102-bib-0011){ref-type="ref"}, [31](#uog19102-bib-0031){ref-type="ref"}. CHD may cause abnormal blood volume and flow patterns, altered chamber and great‐vessel sizes and positions, altered oxygen delivery to organs and increased venous pressure. Blood flow in the ductus venosus and arterial duct is also important. All of these factors were taken into consideration when defining three CHD categories based on the expected cerebral arterial oxygen saturation: expected normal cerebral arterial oxygen saturation, expected mild to moderately reduced cerebral arterial oxygen saturation and expected severely reduced cerebral arterial oxygen saturation. Figure [1](#uog19102-fig-0001){ref-type="fig"} shows the hemodynamics and saturation levels, as reported by Rudolph[13](#uog19102-bib-0013){ref-type="ref"}, [14](#uog19102-bib-0014){ref-type="ref"}, [30](#uog19102-bib-0030){ref-type="ref"}, of three examples of CHD allocated to each of these three groups. Fetuses with different degrees of left ventricular outflow tract obstruction sometimes fell into different categories depending on the intracardiac anatomy and degree of obstruction. In mild aortic stenosis, a normal level of cerebral arterial oxygen saturation is expected. In critical aortic stenosis, retrograde flow over the aortic arch occurs, reducing the cerebral arterial oxygen saturation, hence its allocation to the group with expected mild to moderately reduced cerebral arterial oxygen saturation. Hypoplastic left heart syndrome (HLHS), in which the entire cerebral blood flow may be dependent on retrograde flow in the aortic arch, was allocated to the group with expected severely reduced cerebral arterial oxygen saturation.

![Classification of fetal congenital heart defects according to expected cerebral arterial oxygen saturation, with three examples from each category (adapted from Rudolph)[13](#uog19102-bib-0013){ref-type="ref"}, [14](#uog19102-bib-0014){ref-type="ref"}, [30](#uog19102-bib-0030){ref-type="ref"}. Numbers are saturations (in percent) based on lamb experiments; values may be lower in human fetuses. Arrows indicate flow hemodynamics. In (a--c) are examples of those classified as having expected normal cerebral arterial oxygen saturation: (a) normal heart; (b) coarctation of the aorta; (c) mild to moderate aortic stenosis. In (d--f) are those classified as having expected mild to moderately reduced cerebral arterial oxygen saturation: (d) tetralogy of Fallot; (e) pulmonary atresia with intact septum; (f) critical aortic stenosis. In (g--i) are those classified as having expected severely reduced cerebral arterial oxygen saturation: (g) transposition of the great arteries (TGA) with ventricular septal defect; (h) TGA with intact interventricular septum; (i) hypoplastic left heart syndrome (aortic and mitral atresia). Ao, descending aorta; DA, arterial duct; IVC, inferior vena cava; LA, left atrium; LV, left ventricle; PA, pulmonary artery; PV, pulmonary veins; RA, right atrium; RV, right ventricle; SVC, superior vena cava.](UOG-53-769-g001){#uog19102-fig-0001}

Sensitivity analyses {#uog19102-sec-0010}
--------------------

In the MRI study of Sun *et al*.[11](#uog19102-bib-0011){ref-type="ref"}, the mean aortic oxygen saturation in human fetuses with CHD was lower than that in normal hearts. Furthermore, there were differences in aortic saturation between different types of CHD. Fetuses with tetralogy of Fallot (ToF) had the highest ascending aortic saturation levels of the cyanotic CHDs, followed by transposition of the great arteries (TGA) (including both TGA with intact ventricular septum (IVS) and ventricular septal defect) and single‐ventricle lesions (predominantly HLHS and tricuspid atresia). With this in mind, we repeated the analysis, further subdividing the intermediate (mild to moderately reduced) group into two subgroups, with HLHS reallocated from the severely reduced to the moderately reduced group and only cases of TGA with IVS allocated to the severely reduced group (Table [S1](#uog19102-supitem-0001){ref-type="supplementary-material"}).

As Jansen *et al*.[32](#uog19102-bib-0032){ref-type="ref"} reported that fetuses with ToF had a smaller HC compared with reference values, starting from 20 weeks\' gestation onwards, we also repeated the analyses separately for fetuses with ToF. We performed two different analyses, one with three subgroups (normal, mild to moderately reduced, severely reduced) plus ToF and one with four subgroups (normal, mildly reduced, moderately reduced, severely reduced) plus ToF.

Statistical analysis {#uog19102-sec-0011}
--------------------

All data were entered into an IBM SPSS statistics database (version 23.0, IBM Corp., Armonk, NY, USA). Study population characteristics were described as frequencies (percentage) for categorical variables, mean (± SD) for continuous parameters with an approximately symmetric distribution and median (interquartile range) for continuous data with a skewed distribution. Average trends over time for fetal growth (HC and AC *Z*‐scores) and Doppler flow (MCA‐PI, UA‐PI and CPR *Z*‐scores) were estimated using a linear mixed‐effects model and analyzed with statistical software R (version 3.4.0, Foundation for Statistical Computing, Vienna, Austria). A linear mixed‐effects model takes into account repeat measurements over time and allows the number and timing of the measurements to vary per fetus. To avoid selection bias, all fetuses with at least one measurement were included in the analyses[33](#uog19102-bib-0033){ref-type="ref"}. Average age trends ('fixed effects') were allowed to differ by CHD classification and were modeled by restricted cubic splines. The restricted cubic spline function allowed us to explore the effect of age without making restrictive assumptions about the shape of the time trends. Knots were placed at five fixed quintiles of the predictor\'s distribution, as suggested by Stone and Koo[34](#uog19102-bib-0034){ref-type="ref"}. Fetuses in the study population were considered to be a random sample of the total fetal CHD population. Therefore, we allowed the intercept (i.e. value at birth) and slope to differ per fetus, and assumed these parameters to follow a multivariate normal distribution ('random effects'). In addition, specific time intervals were tested based on the averaged trend. This enabled analysis of the difference in MCA‐PI *Z*‐score trends over time between fetuses with a normal HC and those with an abnormal HC (i.e. HC \< --2.0 SD or \> 2.0 SD) at the last measurement before birth in the total cohort and in the three different CHD categories. It also enabled the analysis of the difference in MCA‐PI trends over time between fetuses with a normal HC/AC ratio and those with an abnormal HC/AC ratio (abnormal HC/AC ratio: \< 5^th^ percentile or \> 95^th^ percentile). Sampling uncertainty was quantified via 95% CI and *P*‐values. *P* \< 0.05 was considered to be statistically significant.

RESULTS {#uog19102-sec-0012}
=======

Patient characteristics {#uog19102-sec-0013}
-----------------------

Included in the study were 181 fetuses with CHD, of which 164 had known outcome. Of these, 144 were liveborn and 122 survived the first 3 months. Maternal characteristics and pregnancy outcome are presented in Table [1](#uog19102-tbl-0001){ref-type="table"}. Six (3%) fetuses suffered an intrauterine death, 13 (7%) pregnancies were terminated and 22 (12%) infants died during the neonatal period or within the first 3 months after birth. Cardiac diagnosis was confirmed by postnatal echocardiographic examination or surgical reports in 156 (86%) cases, by postmortem reports in seven (4%) cases and by fetal echocardiographic examination in 18 (10%) cases. Forty‐four fetuses were allocated to the group with expected normal cerebral arterial oxygen saturation, 84 to the group with expected mild to moderately reduced cerebral oxygen saturation and 53 to the group with expected severely reduced cerebral oxygen saturation (Table [2](#uog19102-tbl-0002){ref-type="table"}). In total, 745 ultrasound examinations were performed in the 181 CHD fetuses. The median number of ultrasound observations per fetus was three, with an interquartile range of two--four and a maximum of 16 observations. Mean (± SD) UtA‐PI *Z*‐scores were --0.44 (± 1.87) for the left and --0.75 (± 2.14) for the right UtA. Fifteen percent of the fetuses had an abnormal UtA‐PI *Z*‐score (\> 2.0).

###### 

Maternal characteristics and pregnancy outcome for study population of 181 cases of fetal congenital heart disease (CHD)

  Characteristic                                                                 Value
  -------------------------------------------------------------------------- -------------
  Maternal                                                                   
  Body mass index at first examination (kg/m^2^)                              24.8 ± 4.6
  Smoker                                                                        17 (9)
  Comorbidity                                                                
  None                                                                         151 (83)
  CHD                                                                            2 (1)
  Diabetes                                                                       6 (3)
  Hypothyroidism                                                                 1 (1)
  Other[\*](#uog19102-note-0002){ref-type="fn"}                                  9 (5)
  Unknown                                                                       12 (7)
  Nulliparous                                                                   61 (34)
  Fetal                                                                      
  Outcome                                                                    
  Survived[†](#uog19102-note-0003){ref-type="fn"}                              122 (68)
  IUFD                                                                           6 (3)
  TOP                                                                           13 (7)
  NND/ID                                                                        22 (12)
  Unknown                                                                       18 (10)
  Gestational age at birth (weeks)[‡](#uog19102-note-0003a){ref-type="fn"}    37.5 ± 4.32
  Birth weight (g)[§](#uog19102-note-0004){ref-type="fn"}                     3055 ± 802
  Head circumference at birth (cm)[¶](#uog19102-note-0005){ref-type="fn"}     33.4 ± 3.16
  5‐min Apgar score[\*\*](#uog19102-note-0006){ref-type="fn"}                  9 ± 1.27

Data are presented as mean ± SD or *n* (%).

Hyperthyroidism, pituitary gland pathology, congenital hepatic fibrosis, migraine.

Beyond 3 months of age.

*n* = 148.

*n* = 142.

*n* = 93.

*n* = 126.

IUFD, intrauterine fetal demise; NND/ID, neonatal/infant death; TOP, termination of pregnancy.

###### 

Classification of congenital heart disease (CHD) according to expected cerebral arterial oxygen saturation in study population of 181 fetuses

  Type of CHD                                                                                                                                                  *n*
  ------------------------------------------------------------------------------------------------------------------------------------------------------------ -----
  Expected normal cerebral arterial oxygen saturation                                                                                                          44
  AVSD (1 with left isomerism, 1 with ASD)                                                                                                                     10
  Isolated VSD/VSDs                                                                                                                                            13
  VSD + ASD (1 with absent ductus venosus)                                                                                                                     2
  Moderate AS (1 with PS)                                                                                                                                      2
  CoA                                                                                                                                                          6
  PS and dysplastic pulmonary valve                                                                                                                            5
  Tricuspid valve abnormalities (antegrade flow over aortic arch)                                                                                              2
  Double discordance with Ebstein\'s anomaly (antegrade flow over aortic arch)                                                                                 1
  Mitral stenosis including small LV (antegrade flow over aortic arch)                                                                                         3
  Expected mild to moderately reduced cerebral arterial oxygen saturation                                                                                      84
  AS (critical, retrograde flow over aortic arch)                                                                                                              1
  CoA with VSD (with or without AS or hypoplastic aortic arch)                                                                                                 10
  Interrupted aortic arch with VSD                                                                                                                             2
  Ebstein\'s anomaly (retrograde flow in arterial duct (with or without PS/PA)) (1 with small VSD)                                                             5
  PS with VSD (*n* = 1)/tetralogy of Fallot (*n* = 14)/DORV with PS (*n* = 4)                                                                                  19
  PA with VSD                                                                                                                                                  4
  Atrioventricular discordance with DORV and CoA                                                                                                               1
  Atrioventricular concordance with mal/transposed great arteries and PS                                                                                       4
  Atrioventricular concordance with mal/transposed great arteries and/or mitral hypoplasia and/or CoA, no PS or Taussig--Bing                                  4
  Transposition of the great arteries with PS and mitral stenosis                                                                                              1
  Common arterial trunk (1 with AVSD)                                                                                                                          7
  PA with IVS                                                                                                                                                  5
  TA (with or without mal/transposition of the great arteries, hypoplastic aortic arch, PS or DORV) (1 with infracardiac TAPVD)                                7
  Hypoplastic right heart (TA or tricuspid stenosis with PA or critical PS with IVS)                                                                           6
  Double inlet left ventricle (with or without mal/transposition of the great arteries with hypoplastic aortic arch or PA) (1 with AVSD and right isomerism)   7
  Complex: right atrial isomerism with unbalanced AVSD, criss‐cross atrioventricular junction, PA with solitary outlet (Ao) from RV and PAPVC                  1
  Expected severely reduced cerebral arterial oxygen saturation                                                                                                53
  HLHS (aortic atresia with or without mitral atresia/stenosis)                                                                                                16
  Transposition of the great arteries with VSD                                                                                                                 10
  DORV with transposition of the great arteries (Taussig--Bing) with or without hypoplastic aortic arch                                                        3
  Transposition of the great arteries with IVS                                                                                                                 24

AS, aortic stenosis; ASD, atrial septal defect; AVSD, atrioventricular septal defect; CoA, coarctation of the aorta; DORV, double outlet right ventricle; HLHS, hypoplastic left heart syndrome; IVS, intact ventricular septum; LV, left ventricle; PA, pulmonary atresia; PAPVC, partial anomalous pulmonary venous connection; PS, pulmonary stenosis; RV, right ventricle; TA, tricuspid atresia; TAPVD, total anomalous pulmonary venous drainage; VSD, ventricular septal defect.

Trends in fetal biometry and Doppler flow {#uog19102-sec-0014}
-----------------------------------------

### *Entire cohort* {#uog19102-sec-0015}

The study included 611 HC, 632 AC, 396 MCA‐PI, 465 UA‐PI and 362 CPR measurements. UtA‐PI was assessed in 86 (48%) cases. All fetuses had at least two HC measurements and 60% had at least two MCA‐PI measurements. There were no differences in baseline characteristics between fetuses with one and fetuses with more than one MCA‐PI measurement. The average trends in *Z*‐scores of fetal biometry (HC and AC) and Doppler flow (MCA‐PI, UA‐PI and CPR) for the entire cohort are shown in Figure [2](#uog19102-fig-0002){ref-type="fig"}. Average trends over time were significant for both HC (*P* \< 0.0001) and AC (*P* \< 0.0001) *Z*‐scores. The HC *Z*‐scores decreased from 20 weeks until 23 weeks (*P* = 0.005), then increased thereafter until 33 weeks (*P* \< 0.0001). After 33 weeks, HC *Z*‐scores decreased again. The AC *Z*‐scores increased progressively during gestation. Doppler flow patterns showed a statistically significant average trend for MCA‐PI (*P* = 0.010) and UA‐PI (*P* \< 0.0001) *Z*‐scores, but not for CPR (*P* = 0.161) *Z*‐scores. The *Z*‐scores of MCA‐PI showed a slight increase between 25 and 30 weeks of gestation, although this trend was not statistically significant (*P* = 0.102), and a decrease after 30 weeks until approximately 35 weeks (*P* = 0.001). Irrespective of the observed trends, the averaged trend lines of all parameters fell within normal ranges (*Z*‐scores \> --2.00 and \< 2.00).

![Averaged trends in *Z*‐scores of fetal growth and Doppler flow patterns in fetuses with congenital heart disease. (a) Head circumference (HC); (b) abdominal circumference (AC); (c) middle cerebral artery pulsatility index (MCA‐PI); (d) umbilical artery pulsatility index (UA‐PI); (e) cerebroplacental ratio (CPR). Gray dots are individual measurements, gray lines are individual trends and dark gray shading shows 95% CIs for fitted model.](UOG-53-769-g002){#uog19102-fig-0002}

### *Subgroups* {#uog19102-sec-0016}

In Figure [3](#uog19102-fig-0003){ref-type="fig"}, the average trends in *Z*‐scores are shown for the three CHD groups categorized according to the expected cerebral arterial oxygen saturation. There were no statistically significant differences in average trends among the CHD groups for HC (*P* = 0.884) or for AC (*P* = 0.879) *Z*‐scores. There were also no differences in the average trends among the groups for UA‐PI (*P* = 0.530) or for CPR (*P* = 0.285) *Z*‐scores. The average trends in MCA‐PI *Z*‐scores were significantly different between the three subgroups (*P* = 0.010): fetuses with an expected mild to moderately reduced or severely reduced cerebral arterial oxygen saturation level showed relatively more fluctuations in MCA‐PI *Z*‐scores (*P* = 0.019 and *P* \< 0.0001, respectively) throughout pregnancy, while fetuses with an expected normal cerebral arterial oxygen saturation level showed no trend in MCA‐PI *Z*‐scores (*P* = 0.831).

![Averaged trends in *Z*‐scores of fetal growth and Doppler flow patterns in three subgroups of fetuses with congenital heart disease, classified according to whether cerebral arterial oxygen saturation would be expected to be normal (![](UOG-53-769-g005.jpg "image")), mild to moderately reduced (![](UOG-53-769-g006.jpg "image")) or severely reduced (![](UOG-53-769-g007.jpg "image")). 95% CIs of fitted model are shaded gray. (a) Head circumference (HC); (b) abdominal circumference (AC); (c) middle cerebral artery pulsatility index (MCA‐PI); (d) umbilical artery pulsatility index (UA‐PI); (e) cerebroplacental ratio (CPR).](UOG-53-769-g003){#uog19102-fig-0003}

Trends in MCA‐PI according to HC before birth {#uog19102-sec-0017}
---------------------------------------------

Figure [4](#uog19102-fig-0004){ref-type="fig"} shows the average trend in MCA‐PI *Z*‐scores in fetuses with normal and those with abnormal HC *Z*‐scores at the last measurement, both for the total cohort and for the three CHD subgroups. Nineteen (10%) fetuses with various types of CHD had an abnormal HC *Z*‐score at the last measurement before birth. There was no statistically significant difference in the average trend of MCA‐PI *Z*‐scores between fetuses with normal HC and those with abnormal HC for the total cohort (*P* = 0.284) or for the CHD subgroups (*P* = 0.363). There was also no difference in the average trend of MCA‐PI *Z*‐scores between fetuses with normal and those with abnormal HC/AC ratio (*P* = 0.505).

![Averaged trends in *Z*‐score of middle cerebral artery pulsatility index (MCA‐PI) for normal head circumference (HC) (solid lines) and for abnormal HC (dashed lines) at last measurement before birth in fetuses with congenital heart disease (CHD): (a) entire cohort; (b) three subgroups according to whether CHD type would be expected to have normal (![](UOG-53-769-g008.jpg "image")), mild to moderately reduced (![](UOG-53-769-g009.jpg "image")) or severely reduced (![](UOG-53-769-g010.jpg "image")) cerebral arterial oxygen saturation. 95% CIs of fitted model are shaded gray for total cohort only; those for CHD subgoups were broad and overlapping.](UOG-53-769-g004){#uog19102-fig-0004}

Sensitivity analyses {#uog19102-sec-0018}
--------------------

The average trends in *Z*‐scores of fetal biometry and Doppler flow patterns from the repeated analysis using four different subgroups (normal, mildly reduced, moderately reduced and severely reduced cerebral arterial oxygen saturation) are shown in Figure S1. The average trends for HC (*P* = 0.656), AC (*P* = 0.707) and CPR (*P* = 0.513) *Z*‐scores did not differ between the four subgroups, whereas trends in MCA‐PI and UA‐PI *Z*‐scores were significantly different between the four subgroups (*P* = 0.028 and *P* = 0.020, respectively). When ToF was analyzed separately from the three subgroups (normal, mild to moderately reduced and severely reduced cerebral arterial oxygen saturation), similar trends in *Z*‐scores of fetal biometry and Doppler flow patterns were observed (Figure S2). The average trends for HC (*P* = 0.950), AC (*P* = 0.701) and CPR (*P* = 0.152) *Z*‐scores did not differ between the subgroups, while the average trends for MCA‐PI (*P* = 0.024) and UA‐PI (*P* \< 0.0001) *Z*‐scores were significantly different between the subgroups. When ToF cases were analyzed separately from the four subgroups (normal, mildly reduced, moderately reduced and severely reduced cerebral arterial oxygen saturation), the results were also similar, with the exception that the MCA‐PI *Z*‐score trend (*P* = 0.073) was no longer significant (Figure S3).

DISCUSSION {#uog19102-sec-0019}
==========

We found that, although HC, AC, MCA‐PI, UA‐PI and CPR remained within normal range in fetuses with CHD, there were significant trends over time in these variables, with the exception of CPR. Trends in MCA‐PI also differed between fetuses grouped according to their expected cerebral arterial oxygen saturation (normal, mild to moderately reduced or severely reduced), while trends in fetal growth did not differ between these subgroups. Fetuses with CHD with expected mild to moderately reduced or severely reduced cerebral arterial oxygen saturation showed more fluctuations in MCA‐PI *Z*‐scores throughout pregnancy than did those with expected normal oxygen saturation levels. Furthermore, we found no significant associations between MCA‐PI and HC *Z*‐scores in fetuses with various types of CHD, subdivided according to the expected cerebral arterial oxygen saturation.

It has been suggested that, in fetuses with CHD, there is preferential blood redistribution to the brain, in order to guarantee an optimal cerebral oxygen supply, as happens in cases of fetal hypoxemia due to placental insufficiency^1,2,4--7,15--25^. Most of these studies used a single Doppler measurement, made during the second or third trimester^4--7,15--25^. We were unable to confirm circulatory redistribution in favor of the brain; in our cohort, Doppler parameters remained within normal ranges. However, in contrast to the situation in chronic hypoxemia, HC *Z*‐scores showed a tendency to decrease at the end of the third trimester, while AC *Z*‐scores continued to increase with advancing gestation. This late head growth impairment suggests that fetuses with CHD may be unable to meet the increased metabolic demands of the developing brain at the end of pregnancy[35](#uog19102-bib-0035){ref-type="ref"}. This hypothesis is further supported by the lack of circulatory compensation, inferred by the constant MCA‐PI *Z*‐scores, despite the lower HC *Z*‐scores.

Previous studies on serial fetal biometry and Doppler measurements in fetuses with CHD show striking discrepancies. While the majority of studies report significant trends over time in both fetal biometry and cerebroplacental Doppler flow, the directions of trends differ significantly between studies[1](#uog19102-bib-0001){ref-type="ref"}, [2](#uog19102-bib-0002){ref-type="ref"}, [32](#uog19102-bib-0032){ref-type="ref"}. Such discrepancies might be caused by differences in study methodology. For instance, while we used the same statistical approach as that of Ruiz *et al*.[1](#uog19102-bib-0001){ref-type="ref"}, they assumed that the effect of gestational age followed a quadratic time trend, which may have produced the tendency for fetal parameters to increase during pregnancy. We made no assumptions regarding the shape of the time trend as we cannot be certain that use of a quadratic trend is justified. Furthermore, discrepancies might be caused by differences in study populations, although the populations are fairly representative of CHD found most commonly during fetal life[1](#uog19102-bib-0001){ref-type="ref"}, [2](#uog19102-bib-0002){ref-type="ref"}, [32](#uog19102-bib-0032){ref-type="ref"}.

Since fetal brain development is dependent on adequate oxygen and nutrient supply, it has been hypothesized that fetuses with CHD in which umbilical venous blood, rich in nutrients and oxygen, is partly shunted away from the brain, have impaired cerebral development[1](#uog19102-bib-0001){ref-type="ref"}, [5](#uog19102-bib-0005){ref-type="ref"}, [6](#uog19102-bib-0006){ref-type="ref"}, [19](#uog19102-bib-0019){ref-type="ref"}, [22](#uog19102-bib-0022){ref-type="ref"}. Several studies have reported lower MCA‐PI and CPR and higher UA‐PI in CHD cases with low expected cerebral oxygen supply[1](#uog19102-bib-0001){ref-type="ref"}, [5](#uog19102-bib-0005){ref-type="ref"}, [6](#uog19102-bib-0006){ref-type="ref"}, [19](#uog19102-bib-0019){ref-type="ref"}, [22](#uog19102-bib-0022){ref-type="ref"}. We were unable to confirm this circulatory redistribution, but found more fluctuations in MCA‐PI *Z*‐scores throughout pregnancy in CHD fetuses with mild to moderately reduced or low expected cerebral oxygen supply compared with fetuses with CHD with normal expected oxygen supply, a potentially important observation. Although we cannot exclude the possibility that this was due to physiological variation, as longitudinal studies in healthy fetuses are lacking, we speculate that, if a manifestation of abnormality, these fluctuations might affect adversely development of the brain. They might suggest impaired hemodynamic autoregulation, which, similar to the situation in preterm infants[36](#uog19102-bib-0036){ref-type="ref"}, [37](#uog19102-bib-0037){ref-type="ref"}, may have a harmful effect on the vulnerable developing brain.

There were no associations between MCA‐PI and HC *Z*‐scores, neither in the entire cohort nor in any of the three CHD subgroups. Previous studies were also unable to demonstrate a clear association between MCA‐PI and HC. Furthermore, Jansen *et al*.[32](#uog19102-bib-0032){ref-type="ref"} found no association between the flow in the ascending aorta and HC. These studies, however, either were based on univariate statistical analyses or used only theoretical hemodynamics. We used multivariate statistical analyses and measured the MCA‐PI. Therefore, we believe that our study is more robust, and we speculate that head growth might be more dependent on other factors, such as maternal or (epi)genetic factors, or nutrients (glucose)[31](#uog19102-bib-0031){ref-type="ref"}, than on blood flow to the brain.

This study has several strengths and limitations. It is the first study to assess the association between trends in Doppler flow and fetal growth in a relatively large cohort of consecutive fetuses with CHD. Furthermore, we included an unselected population (the assessment of fetal growth and cerebroplacental Doppler flow was part of routine clinical care in both institutions) and we excluded all cases with chromosomal abnormalities, including microdeletions. However, this was a retrospective study, operators were not blinded to the presence of CHD, which might have affected Doppler and biometry measurements, and the assessment of Doppler flow was not standardized, resulting in various numbers of measurements across the range of gestational ages. Nonetheless, there were no indications that missing data were not at random. Therefore, our linear mixed‐effects model provides an unbiased estimate of the age trends. In addition, we were unable to confirm cardiac diagnosis using postnatal echocardiograms, surgery reports or postmortem reports in 10% of the cases. Also, comparison with other studies is limited by differences in study populations, reference charts used and study design. We believe that our classification of fetuses into normal, mild to moderately reduced and severely reduced expected cerebral arterial oxygen saturation is based on the best data available[11](#uog19102-bib-0011){ref-type="ref"}, [13](#uog19102-bib-0013){ref-type="ref"}, [14](#uog19102-bib-0014){ref-type="ref"}, [30](#uog19102-bib-0030){ref-type="ref"}, [31](#uog19102-bib-0031){ref-type="ref"}; however, this remains speculative, as it is impossible to measure directly oxygen levels prenatally. We thus also performed a sensitivity analysis with slightly different classification of CHD, but the results remained unchanged. Finally, we were unable to analyze most types of CHD separately due to the relatively small numbers; the exceptions were ToF, previously reported to be associated with smaller HC[32](#uog19102-bib-0032){ref-type="ref"}, and TGA + IVS, for which the lowest level of cerebral arterial oxygen saturation is expected. Multicenter studies providing larger numbers with the same type of CHD are necessary to provide insights into possible mechanisms responsible for suboptimal intrauterine cerebral development in a specific heart lesion.

In conclusion, this study demonstrates that, while there are significant trends in biometry and Doppler flow throughout pregnancy in fetuses with CHD, these measurements are within normal range. Fetal head growth is not associated with the expected cerebral arterial oxygen saturation or fetal Doppler flow, confirming that mechanisms other than circulatory modifications may influence the cerebral development in fetuses with CHD.

Supporting information
======================

###### 

**Table S1** Alternative classification of congenital heart disease according to expected cerebral arterial oxygen saturation, including four subgroups

###### 

Click here for additional data file.

###### 

**Figures S1--S3** Averaged trends in *Z*‐scores of fetal growth and Doppler flow patterns in different subgroups of fetuses with congenital heart disease, classified according to expected level of cerebral arterial oxygen saturation: normal, mildly reduced, moderately reduced or severely reduced (Figure S1); normal, mildly to moderately reduced, severely reduced or tetralogy of Fallot (Figure S2); and normal, mildly reduced, moderately reduced, severely reduced or tetralogy of Fallot (Figure S3).

###### 

Click here for additional data file.

We thank Prof. A. M. Rudolph for his inspiring discussions that helped in finalizing the CHD classification. This study was part of the research program of the Graduate School of Medical Sciences, Research Institutes BCN‐BRAIN and GUIDE. M.J.M. and W.S.K. were supported financially by a University of Groningen Junior Scientific Masterclass grant.
