Early postnatal docosahexaenoic acid levels and improved preterm brain development.

BACKGROUND: Preterm birth has a dramatic impact on polyunsaturated fatty acid exposures for the developing brain. This study examined the association between postnatal fatty acid levels and measures of brain injury and development, as well as outcomes.
METHODS: A cohort of 60 preterm newborns (24-32 wk gestational age) was assessed using early and near-term magnetic resonance imaging (MRI) studies. Red blood cell fatty acid composition was analyzed coordinated with each scan. Outcome at a mean of 33 mo corrected age was assessed using the Bayley Scales of Infant Development, 3rd edition.
RESULTS: Adjusting for confounders, a 1% increase in postnatal docosahexaenoic acid (DHA) levels at early MRI was associated with 4.3-fold decreased odds of intraventricular hemorrhage, but was not associated with white matter injury or cerebellar haemorrhage. Higher DHA and lower linoleic acid (LA) levels at early MRI were associated with lower diffusivity in white matter tracts and corresponding improved developmental scores in follow-up.
CONCLUSION: Higher DHA and lower LA levels in the first few weeks of life are associated with decreased intraventricular haemorrhage, improved microstructural brain development, and improved outcomes in preterm born children. Early and possibly antenatal interventions in high-risk pregnancies need to be studied for potential benefits in preterm developmental outcomes.

INTRODUCTION

Preterm birth is associated with increased risk for adverse neurodevelopmental impairment in survivors. Brain imaging has provided many clues as to the mechanisms of impaired outcomes, and brain injury patterns such as intraventricular haemorrhage (IVH),(1,2) white matter injury (WMI),(3,4) and cerebellar haemorrhage(5,6) have been shown to be predictive of motor and cognitive deficits. Advanced neuroimaging techniques such as diffusion tensor imaging provide tools to study brain microstructural development.(7) These methods have enabled our understanding of risk factors for impaired brain development, such as infection.(4)

Preterm birth, with premature transition to parenteral nutrition, abruptly changes the nutritional exposures of the developing brain. Placental transfer of fatty acids in the third trimester is low, about 10% of total energy, but comprises a much higher proportion of long-chain polyunsaturated fatty acid (PUFA) than in maternal plasma.(8) Docosahexaenoic acid (DHA), a key essential PUFA for the nervous system, falls to one-third whole blood in utero levels within 2–3 weeks after preterm birth.(9) Even late preterm birth is associated with significant disruption of normal essential PUFA levels, including lower ω-3 fatty acid levels.(10) Essential PUFAs cannot be synthesized de novo, but are important components of cell membranes. Notably, DHA may play a role in decreased brain injury(11,12) and improved development.(13–15) In addition, the ω-6 PUFA arachidonic acid (ARA) is a precursor to hormones that regulate inflammation, including prostaglandins and leukotrienes, with additional potential effects.

Although meta-analyses have not shown overall developmental improvements with postnatal supplementation of DHA and ARA,(16) specific early benefits of DHA with or without ARA have been demonstrated in visual function(17) and some cognitive tests.(18–21) One limitation of these prior studies is the variability in dosage and timing of interventions. Few studies have looked at specific neuroimaging markers for the neurodevelopmental effects of PUFAs to determine the appropriate timing and targeting of intervention to effect improvements. This study aims to assess the associations between early postnatal red blood cell PUFA levels and signs of brain injury (including IVH, WMI, and cerebellar haemorrhage) and development using serial MRI studies. It is hypothesized that higher levels of ω-3 fatty acids, specifically DHA, would be associated with decreased risk for brain injury and improved brain development and neurodevelopmental outcomes in infants born very preterm.

RESULTS

A total of 60 preterm newborns <32 weeks GA were enrolled from the two study centers (Table 1). All newborns completed early MRI scan (mean 31.54±2.26 weeks GA) and associated blood sampling, with 44/60 (73%) completing the near-term age MRI scan (mean 37.71±2.98 weeks GA) and associated blood sampling. Infants who were discharged home prior to term and unable to return for outpatient testing were missed. Early blood sample showed a median DHA level of 4.4% (range 2–6.4%), eicosapentaenoic acid (EPA) level of 0.3% (range 0.1–4%), ARA level of 16.4% (range 10.8–20.4%), and linoleic acid (LA) level of 9.2% (range 5.9–14.1%). Near-term age samples showed a median DHA level of 4.5% (range 0.7–7.7%), EPA level of 0.4% (range 0.1–0.8%), ARA level of 16.6% (3.7–19.2%), and LA level of 10.4% (range 6.2–14.5%).

Table 1

Patient demographics at UCSF and UBC.

	UCSF (N=33)	UBC (N=27)
Gestational age at birth (weeks)	28.05 ± 1.75	28.49 ± 2.03
5-minute Apgar score	8 (IQR 7–8)	8 (IQR 6–8)
Umbilical artery pH	7.26 ± 0.11	7.44 ± 1.33
Male sex	20 (61%)	17 (63%)
Antenatal glucocorticoid exposure (any)	33 (100%)	17 (63%)
Clinical chorioamnionitis	5 (15%)	5 (19%)
Neonatal sepsis	17 (52%)	18 (67%)
Patent ductus arteriosus	16 (48%)	10 (37%)
Hypotension	21 (64%)	21 (78%)
Intubation (days)	11.64 ± 26.10	10.70 ± 20.25
Parenteral nutrition (days)	15.45 ± 8.17	14.67 ± 10.39
Intralipid supplementation (days)	13.61 ± 7.93	14.15 ± 10.45
Enteral nutrition
-Breast milk (days) / total days from birth to near-term blood sample	78% (IQR 63–92%)	96% (IQR 70–99%)
-Formula (days) / total days from birth to near-term blood sample	12% (IQR 0–57%)	0% (IQR 0–43%)
Postmenstrual age at 1st MRI (weeks)	31.04 ± 1.65	32.14 ± 2.75
Postmenstrual age at 2nd MRI (weeks)	36.13 ± 2.42	40.22 ± 1.85
Intraventricular hemorrhage
-Mild (grade 1–2)	4 (12%)	5 (19%)
-Severe (grade 3–4)	6 (18%)	1 (4%)
White matter injury
-Mild	2 (6%)	5 (19%)
-Moderate/severe	6 (18%)	3 (11%)
Cerebellar hemorrhage	4 (12%)	2 (7%)

Expressed as mean ± standard deviation, mean and intraquartile range (IQR), or frequency and percentage.

PUFAs and brain injury on conventional MRI scanning

Fatty acid levels were compared with findings of IVH, WMI, and cerebellar hemorrhage. Univariate logistic regression analysis showed trends to association with IVH for DHA (OR 0.48, 95% CI 0.22–1.06, P=0.070) and ARA (OR 0.77, 95% CI 0.56–1.05, P=0.098), but no other fatty acids and WMI or cerebellar hemorrhage (P>0.3). After adjusting for confounders, only DHA remained significantly associated with IVH (Table 2). One percent increase in DHA was associated with a 4.3-fold decreased odds of IVH (P=0.024). Of note, brain injury scores did not change from first to second scan, however no association was found between near-term fatty acid levels and brain injury patterns including IVH, WMI, or cerebellar haemorrhage (P>0.2).

Table 2

Logistic regression analysis comparing red blood cell fatty acid levels (percent total fatty acid) and presence of brain injury (of all severity grades) on testing as soon after birth as clinically stable for MRI scanning.

	OR	95% CI	P-value
DHA, 22:6ω3
IVH	0.23	0.06–0.83	0.02
WMI	0.83	0.24–2.86	0.76
Cerebellar hemorrhage	0.69	0.10–5.06	0.72
EPA, 20:5ω3
IVH	2.68	0.38–18.93	0.32
WMI	0.17	0.01–3.66	0.26
Cerebellar hemorrhage	21.72	0.02–20652.18	0.38
ARA, 20:4ω6
IVH	0.82	0.58–1.16	0.26
WMI	0.88	0.59–1.30	0.51
Cerebellar hemorrhage	0.37	0.10–1.27	0.11
LA, 18:2ω6
IVH	1.20	0.84–1.71	0.31
WMI	0.92	0.60–1.42	0.71
Cerebellar hemorrhage	a	a	a

Results are adjusted for potential confounders, including gestational age at birth, chorioamnionitis, antenatal betamethasone exposure, 5-minute Apgar score, patent ductus arteriosus, hypotension, intubation, and neonatal sepsis.

Denotes inadequate data for statistical analysis.

Lower DHA levels were associated with higher IVH severity (Figure 1). Using ordinal logistic regression analysis, 1% increase in DHA was associated with 2.3-fold decreased odds of increasing one degree of IVH severity (95% CI 1.06–5.27, P=0.036), when IVH was categorized as none, mild (grade 1–2), or severe (grade 3–4). Adjusting for confounders, 1% increase DHA was associated with 3.9-fold decreased odds of increasing one degree of IVH severity (95% CI 1.22–12.28, P=0.022).

Figure 1

Percent red blood cell DHA and IVH

Boxplot is divided into infants with no, mild (grade 1–2), or severe (grade 3–4) IVH.

PUFAs and brain microstructure on DTI

Using tract-based spatial statistics, fatty acid levels on early preterm and near-term age blood samples were compared with the diffusion tensor imaging measures of fractional anisotropy (FA), axial diffusivity, and radial diffusivity, adjusting for age at scan. No statistically significant associations were found for any fatty acid levels in the near-term age blood sample. In the early preterm blood sample, only DHA was found to be associated with axial and radial diffusivity (Figure 2). No association was found with fractional anisotropy.

Figure 2

Axial fractional anisotropy maps showing the white matter skeleton in green

Regions of the white matter skeleton showing significant associations between early preterm red blood cell docosahexaenoic acid levels and (A) fractional anisotropy, (B) axial diffusivity, and (C) radial diffusivity are highlighted in shades ranging from yellow (P=0.05) to red (P=0.00). Regions of interest in the PLIC and OR are marked in blue.

To compliment these tract-based spatial statistics findings, regions of interest were drawn in the posterior limbs of the internal capsule (PLIC) and the optic radiations (OR), both functionally relevant pathways, as indicated in blue in Figure 2. Diffusion measures in both time points were compared with early preterm fatty acid levels. Results after adjusting for clinical confounders are described in Table 3. Of note, higher DHA levels were inversely associated with mean diffusivity in the PLIC and OR, and axial diffusivity in the PLIC. Higher LA levels were associated with higher mean diffusivity in the OR.

Table 3

Linear mixed regression analysis comparing early preterm red blood cell fatty acid levels (percent total fatty acid) and diffusion tensor imaging measures on serial MRI scanning in the PLIC and the OR.

	PLIC			OR
	Change	95% CI	P-value	Change	95% CI	P-value
DHA, 22:6ω3
Mean diffusivity (10−3 mm2/s)	−0.0353	−0.0627 – −0.0078	0.01	−0.0535	−0.0956 –−0.0115	0.01
Fractional anisotropy	−0.0021	−0.0218 – 0.0176	0.83	0.0106	−0.0160 – 0.0372	0.44
Axial diffusivity (10−3 mm2/s)	−0.0449	−0.0787 –−0.0111	0.01	−0.0477	−0.1115 – 0.0161	0.14
Radial diffusivity (10−3 mm2/s)	−0.0297	−0.0616 – 0.0023	0.07	−0.0445	−0.0953 – 0.0063	0.09
EPA, 20:5ω3
Mean diffusivity (10−3 mm2/s)	−0.0039	−0.1529 – 0.1451	0.96	−0.0339	−0.2614 – 0.1937	0.77
Fractional anisotropy	−0.0638	−0.1603 – 0.0328	0.20	0.1279	−0.0050 – 0.2608	0.06
Axial diffusivity (10−3 mm2/s)	−0.0881	−0.2679 – 0.0918	0.34	0.1631	−0.1717 – 0.4979	0.34
Radial diffusivity (10−3 mm2/s)	0.0366	−0.1307 – 0.2040	0.67	−0.1285	−0.3937 – 0.1367	0.34
ARA, 20:4ω6
Mean diffusivity (10−3 mm2/s)	−0.0055	−0.0178 – 0.0069	0.39	−0.0025	−0.0214 – 0.0163	0.79
Fractional anisotropy	−0.0024	−0.0107 – 0.0058	0.56	−0.0009	−0.0122 – 0.0104	0.88
Axial diffusivity (10−3 mm2/s)	−0.0077	−0.0227 – 0.0074	0.32	−0.0049	−0.0325 – 0.0227	0.73
Radial diffusivity (10−3 mm2/s)	−0.0043	−0.0183 – 0.0097	0.55	−0.0014	−0.0236 – 0.0208	0.90
LA, 18:2ω6
Mean diffusivity (10−3 mm2/s)	0.0086	−0.0016 – 0.0189	0.10	0.0162	0.0005 – 0.0319	0.04
Fractional anisotropy	−0.0033	−0.0103 – 0.0038	0.36	−0.0037	−0.0137 – 0.0064	0.48
Axial diffusivity (10−3 mm2/s)	0.0060	−0.0070 – 0.0189	0.37	0.0097	−0.0147 – 0.0342	0.44
Radial diffusivity (10−3 mm2/s)	0.0100	−0.0015 – 0.0215	0.09	0.0132	−0.0056 – 0.0321	0.17

Change values represent change in diffusion tensor imaging measures per 1% increase in red blood cell fatty acid level. Results are adjusted for postmenstrual age at time of each scan, IVH severity, WMI severity, sepsis, duration of intubation, and patent ductus arteriosus.

Due to the association between DHA and IVH, subanalysis was performed excluding infants with grade 3–4 IVH (Table 4). Excluding these infants did not meaningfully change the associations found in the PLIC, but did strengthen the associations found in the OR, with additional statistically significant association seen between higher DHA levels and lower radial diffusivity. As in prior analyses, no statistically significant associations were found with EPA or ARA. However after excluding infants with grade 3–4 IVH, LA no longer had a statistically significant association with mean diffusivity in the OR.

Table 4

Linear mixed regression analysis comparing early preterm red blood cell DHA levels (percent total fatty acid) and diffusion tensor imaging measures on serial MRI scanning in the PLIC and the OR, excluding infants with grade 3–4 IVH.

	PLIC			OR
	Change	95% CI	P-value	Change	95% CI	P-value
Mean diffusivity (10−3 mm2/s)	−0.0359	−0.0663 –−0.0055	0.02	−0.0747	−0.1183 –−0.0311	<0.001
Fractional anisotropy	−0.0021	−0.0233 – 0.0190	0.84	0.0149	−0.0152 – 0.0449	0.33
Axial diffusivity (10−3 mm2/s)	−0.0439	−0.0813 –−0.0066	0.02	−0.0609	−0.1320 – 0.0102	0.09
Radial diffusivity (10−3 mm2/s)	−0.0313	−0.0667 – 0.0042	0.08	−0.0642	−0.1194 – 0.0090	0.02

Change values represent change in diffusion tensor imaging measures per 1% increase in red blood cell DHA level. Results are adjusted for postmenstrual age at time of each scan, IVH severity, WMI severity, sepsis, duration of intubation, and patent ductus arteriosus.

PUFAs and postnatal nutrition

Table 5 details the comparisons between fatty acid levels and postnatal nutritional exposures prior to each blood sample. Since intravenous nutrition and lipid supplementation were completed prior to the first scans and blood samples, the number of days of these interventions were compared with the early samples. Since the number of days between birth and near-term testing varied significantly between subjects of varied gestational age at birth, the percent of days exposed to breast milk or formula was compared to near-term samples. There was a trend to lower early DHA levels with increased duration of IV nutrition, however no significant associations were found between early and near-term DHA levels and nutrition modalities. Duration of intravenous nutrition was associated with higher early EPA and LA levels, and lower early ARA levels. Percentage of days of formula feeding was associated with lower term EPA levels.

Table 5

Change in early fatty acid levels compared with intravenous nutrition and late fatty acid levels compared with breast milk versus formula feeding.

	% change	95% CI	P-value
Early preterm red blood cell fatty acid levels
DHA, 22:6ω3
Days intravenous nutrition	−2.19%	−4.62 – 0.25%	0.08
Days lipid supplementation	−1.62%	−4.06 – 0.82%	0.19
EPA, 20:5ω3
Days intravenous nutrition	12.33%	8.72 – 15.96%	<0.001
Days lipid supplementation	12.57%	9.07 – 16.07%	<0.001
ARA, 20:4ω6
Days intravenous nutrition	−2.07%	−3.20 –−0.93%	0.001
Days lipid supplementation	−2.00%	−3.13 –−0.88%	0.001
LA, 18:2ω6
Days intravenous nutrition	1.96%	0.85 – 3.06%	0.001
Days lipid supplementation	1.88%	0.78 – 2.98%	0.001
Near-term age red blood cell fatty acid levels
DHA, 22:6ω3
% days breast milk	−0.01%	−0.06 – 0.03%	0.57
% days formula	−0.02%	−0.05 – 0.003%	0.08
EPA, 20:5ω3
% days breast milk	−0.04%	−0.35 – 0.27%	0.80
% days formula	−0.19%	−0.37 – −0.02%	0.03
ARA, 20:4ω6
% days breast milk	0.003%	−0.02 – 0.02%	0.73
% days formula	0.003%	−0.01 – 0.01%	0.57
LA, 18:2ω6
% days breast milk	−0.01%	−0.04 – 0.03%	0.72
% days formula	−0.001%	−0.02 – 0.02%	0.92

Only nutrition from days prior to each blood sample were counted.

PUFAs and developmental outcome

Developmental assessments were performed on 25 (76%) of the enrolled children at UCSF and 20 (74%) at UBC. Assessment was at a mean of 32.8±3.6 months corrected age. There was no difference in fatty acid levels between those assessed in follow-up and those lost to follow-up (P>0.09). Mean Bayley-3 cognitive scores were 103.6±14.4, language scores 101.0±17.8, and motor scores 98.1±14.3. Although there was no association between early DHA levels and Bayley-3 cognitive scores (P=0.4), a 1% increase in early DHA level was associated with 7.0 point higher language scores (95%CI 1.3–12.7, P=0.017) and 7.5 point higher motor scores (95%CI 3.1–11.9, P=0.001). A 1% increase in early LA level was associated with a 2.4-point decrease in cognitive scores (95%CI 0.16–4.6, P=0.036) and 3.5-point decrease in motor scores (95%CI 1.3–5.6, P=0.003), but no association with language scores (P=0.16). Only grade 3–4 IVH was found to be associated with lower Bayley-3 scores (P<0.02). Adjusting for grade 3–4 IVH, early DHA levels continued to be associated with better motor scores (4.8, 95%CI 0.5–9.1, P=0.03) and early LA levels associated with poorer motor scores (−2.2, 95%CI −4.3 – −0.13, P=0.038). No associations were seen with cognitive or language scores. Early ARA and EPA levels were not associated with outcome (P>0.35).

Using univariate linear regression analysis to assess the contribution of socioeconomic status, early DHA levels were associated with maternal education (0.08% per year, 95%CI 0.0004–0.15, P=0.049), maternal age (0.05% per year, 95%CI 0.01–0.08, P=0.013), and Hispanic ethnicity (−1.10, 95%CI −1.56 – −0.65, P<0.001). Race was not associated with early DHA levels (P>0.17). No demographic variables were associated with early LA levels (P>0.07). Using a backwards selection model to determine the contribution of each variable on outcomes, only DHA remained significantly associated with language scores (7.0, 95%CI 1.3–12.7, P=0.017). However, both DHA (6.7, 95%CI 1.7–11.6, P=0.01) and years of maternal education (1.13, 95%CI 0.03–2.2, P=0.044) were significantly associated with motor scores.

DISCUSSION

This study used serial blood samples, neuroimaging, and developmental assessment to present new evidence for an association between PUFA levels and preterm IVH and neurodevelopment. It has previously been suggested that ω-3 fatty acid exposures, specifically DHA, may have positive effects on neurodevelopmental outcomes after preterm birth, yet the mechanisms of this protective effect have been unknown. Our study shows that higher early preterm red blood cell DHA levels are associated with lower incidence and severity of IVH. Measures of microstructural development suggest opposing associations with the ω-3 fatty acid DHA and the ω-6 fatty acid LA. Outcomes assessed at 30–36 months further support these opposing associations.

Investigating the association between PUFA levels and brain injury, only early DHA levels were found to be associated with decreased incidence and severity of IVH. Although an association between PUFA levels and preterm brain injury has been previously postulated,(12) there is currently no reported literature in humans or animal models for such an association.

Understanding the mechanisms for IVH after preterm birth is important for understanding the timing and potential mechanisms for an association with early DHA levels. Diagnosis of IVH usually occurs in the first day of life, and at least 90% of events are identified within the first 72 hours.(22) Thus, if DHA level was to have an effect on IVH incidence, it would have to be evident by the time the child was born, prior to hemorrhage. Accordingly, this study found that only the early DHA levels (collected on average of 3 weeks after birth) were associated with IVH.

DHA accounts for over 50% of all PUFAs in the central nervous system cell membranes.(23) The specific molecular properties of this highly unsaturated fatty acid result in important changes in phase behaviour, fluidity, permeability, and elasticity of cell membranes.(24) Hemorrhage from IVH is believed to originate from the fragile capillary network of the subependymal germinal matrix, resulting from fluctuations in blood pressure and impaired cerebral autoregulation.(25,26) This highlights the importance of DHA in the central nervous system, and its potential effects on membrane elasticity and resistance to pressure changes that could affect the risk for haemorrhage in the preterm newborn.

Independent of IVH or WMI, early DHA levels were also found to be associated with long-term effects on white matter microstructural development. It has been previously shown in white matter tracts that axial and radial diffusivity decrease, while fractional anisotropy increases, with preterm brain development.(27) Increases in fractional anisotropy may be partially related to increased premyelination and myelination, limiting radial mobility of water. Decreasing axial and radial diffusivity, however, may be due to decreasing brain water content and formation of new barriers to water mobility, including increased macromolecular concentration and outgrowth of neurons, glia, and blood vessels in all directions.(28) Thus, the observed decrease in mean diffusivity with higher DHA and lower LA, without significant change in fractional anisotropy, suggests improved microstructural development independent of myelination. These findings are further supported by studies in macaque monkeys showing impairments in cortical network development after a lifetime of low dietary ω-3 fatty acids.(29)

Investigating the contribution of the neonatal feeding practices, although duration of parenteral nutrition had effects on EPA, LA, and ARA levels, no effects were seen with DHA levels. Lower EPA levels were seen with increased formula-feeding, which is not surprising since Similac Special Care contains no EPA. Data collected in this study in regards to nutritional sources is limited to days exposed to each modality, and thus only estimates the exact exposures for each child. Finding lower EPA levels with formula feeding suggests that our estimated exposures would be able to detect trends in feeding differences. Our results suggest the feeding regimens used during this study have no significant effect on postnatal DHA levels. Indeed, a Cochrane review of lipid supplementation showed no positive or negative association with outcomes.(30)

This study also suggests that higher early DHA and lower early LA levels have long-term associations with improved developmental outcomes at 30–36 months. These associations can be partially, but not fully explained by decreased rates of severe IVH. Markers of low socioeconomic status, including low maternal education, low maternal age, and Hispanic ethnicity were associated with lower early DHA levels, suggesting effects of socioeconomic status may begin antenatally, and not just after hospital discharge. Lack of association with race suggests these differences are not as likely due to genetic variation. These results agree with prior studies showing a relationship between socioeconomic status and DHA levels in adults.(31–33) This highlights the importance of the antenatal maternal diet on perinatal outcomes. Despite adjustment for socioeconomic status however, early DHA and LA levels remained independently associated with outcomes, suggesting that in spite of potential ongoing differences in upbringing between these children, the early perinatal exposures remain an important determining factor.

The inclusion of two study sites enabled not only more rapid subject enrolment, but also an inclusion of a wider range of patient demographics and socioeconomic status. All analyses were adjusted for study site, and significant associations between fatty acid and brain metrics despite adjustment for illness severity suggests that these associations are independent of study site. It is interesting to note, however, that most of the Hispanic subjects in this study come from the UCSF study site, due to the demographics of San Francisco. With lower DHA levels associated with Hispanic ethnicity, it is thus not surprising based on current study findings that UCSF also demonstrated higher rates of IVH, even despite the potentially protective effects of antenatal glucocorticoid treatment at this study site.

For this study, red blood cell PUFA levels were used as a marker of overall fatty acid status. Compared to direct measurements of fatty acid levels within the brain tissue itself, red blood cell fatty acid levels are clearly less accurate. Red blood cell levels were assessed, rather than plasma levels, since plasma levels would be more dependent on day-to-day enteral and parenteral exposures in the neonatal period. However, red blood cell PUFAs are commonly used in animal and human studies as a marker for DHA and neurological development in infants.(17,34) In aged rats supplemented with DHA and EPA, good correlation has been found between red blood cell and brain DHA levels.(35)

The results of this prospective observational study initially may seem to contradict prior studies of postnatal ω-3 fatty acid supplementation, which did not show developmental benefits.(16) However, our current findings may provide clues as to the appropriate timing and selection of targets for intervention. Higher early DHA levels are associated with lower rates and severity of IVH, while higher DHA and lower LA levels are associated with improved microstructural brain development and developmental outcomes at 30–36 months. Since IVH occurs within 72 hours of life, improvements in fatty acid profiles must be effected before this time, likely before birth. In addition, this study found lower early DHA levels to be associated with lower socioeconomic status, suggesting pre-existing and ongoing effects of maternal diet. Thus, further therapeutic trials will need to ensure early intervention and appropriate consideration for families with low socioeconomic status, to ensure appropriate targeting of vulnerable infants. Targeting of intervention and education to expectant mothers with low socioeconomic status may be the most effective for improving fatty acid levels and their implications on functional developmental outcomes.

METHODS

Study participants

A prospective cohort of preterm newborns born between March 2010 and November 2011 admitted to the neonatal intensive care units at either the University of California San Francisco (UCSF) or the University of British Columbia (UBC) was studied. Parents of neonates born <32 weeks GA were approached and consented following a protocol approved by research ethics board at both institutions. Infants were excluded if there was a congenital malformation or syndrome, congenital infection, or the infant was too clinically unstable to transport to the MRI scanner. Clinical history was collected prospectively from the patient charts. Exposure to parenteral nutrition, intravenous lipid supplementation, breast milk nutrition, and formula nutrition were estimated by recording the number of days infants were exposed to each. Infants were first introduced to parenteral nutrition without lipid supplementation, soon followed within a day by introduction of intravenous fat emulsion (Intralipid, Baxter Healthcare Corporation, Deerfield IL), which contains LA (18:2ω6), oleic acid (18:1ω9), and α-linolenic acid (18:3ω3), with smaller amounts of palmitic acid and stearic acid (18:0). Intravenous nutrition was weaned off with tolerance for enteral feeding. Enteral feeding was occasionally supplemented for protein, but not lipids. Breast milk was used whenever available, otherwise formula (Similac Special Care, Abbott Nutrition, Abbott Park IL) was used.

MRI studies

MRI studies were performed on 1.5-T MRI scanners (General Electric Sigma, GE Medical Systems, Milwaukee WI, or Siemens Avanto, Siemens Medical Solutions USA, Malvern PA), as soon after birth as clinically stable and again near near term age. Conventional MRI sequences at both study sites were obtained as previously described, including T2-weighted and volumetric 3-dimensional T1-weighted sequences.(36) Scans were graded by a single neuroradiologist at each site (KJP, AJB) blinded to patient history for brain injury, including IVH, WMI, and cerebellar hemorrhage. IVH was scored using the grading system of Papile et al. (2) WMI was scored using a system shown to be associated with neurodevelopmental outcomes at 12 to 18 months, demonstrated to have good interrater reliability between these two neuroradiologists.(3,37) Injury severity was graded as the highest score from both scans.

Diffusion tensor imaging sequences were acquired using comparable(4,7) multi-repetition, single-shot echo planar sequences (UCSF 6 gradient directions; TR, 7s; TE, 100ms; slice thickness, 3mm and UBC 12 gradient directions; TR, 4900ms; TE, 104ms; FOV, 160ms; slice thickness, 3mm) to determine FA and eigenvalues (λ1, λ2, λ3) in three dimensions describing axial (λ1), radial diffusivity (mean of λ2 and λ3), and mean diffusivity (mean of λ1, λ2, and λ3).(27) Fractional anisotropy describes the directionality of water diffusion. Axial diffusivity has been used as a measure of organization of axonal tracts,(38) while radial diffusivity has been used as a measure of the quality of myelination.(39) Data were analyzed blinded for patient history.

Blood samples

Two blood samples were each drawn as close in time to each MRI scan as possible, timed with clinically-indicated bloodwork and within a week of MRI. Samples were separated into plasma and red blood cell fractions and stored at −70°C until analysis. Proportions of each fatty acid relative to total fatty acids in the red cell membrane were determined by gas chromatography-flame ionization detector.(17) Analysis was then focused to key PUFAs for their impact on brain development, including the ω-3 fatty acids DHA (22:6ω3) and EPA (20:5ω3) and the ω-6 fatty acids LA (18:2ω6) and ARA (20:4ω6).

Neurodevelopmental outcome

Infants were reassessed between 30–36 months corrected age by a developmental psychologist or physiotherapist blinded to the child’s neonatal course. Standardized assessment was performed using the Bayley Scales of Infant Development, 3rd Edition (Bayley-3), and composite cognitive, language, and motor scores were assigned based on the child’s corrected age. Testing was performed in English or Spanish as appropriate. The Bayley-3 has a mean of each composite score of 100 and a standard deviation of 15. Socioeconomic status variables were collected at follow-up from parental self-report, including total years of maternal education, maternal age at childbirth, race, and ethnicity. Race and ethnicity are defined as per the policies of the National Institutes of Health, where ethnicity includes those who are or are not Hispanic/Latino, and race includes American Indian/Alaskan Native, Asian, Black/African American, Native Hawaiian/Other Pacific Islander, or White.

Statistical analysis

Statistical analysis was performed using R version 3.1.1 (The R Foundation for Statistical Computing 2014, Vienna, Austria). All analyses were adjusted for study site to account for site-specific variations. To study the association between fatty acid levels and measures of brain injury (including IVH, WMI, and cerebellar hemorrhage), logistic regression analysis was used to adjust for known clinical risk factors for brain injury, including GA at birth (also associated with fatty acid levels)(10), chorioamnionitis (clinically suspected or confirmed), antenatal glucocorticoid exposure, 5-minute Apgar score, patent ductus arteriosus (diagnosed by echocardiogram, with or without intervention), hypotension requiring medical intervention, intubation, and neonatal sepsis (presumed or confirmed, including necrotizing enterocolitis). Statistical significance was set at P<0.05.

Tract-based spatial statistics (FMRIB Software Library version 5.0, http://www.fmrib.ox.ac.uk/fsl/, Oxford, UK) were used to assess potential associations between diffusion tensor imaging parameters and fatty acid levels in early and near-term blood samples. Data from all infants were fitted into a common space using nonlinear registration, allowing for voxelwise cross-subject statistical analysis of the white matter tract skeleton, using a previously validated method for this age group.(40) Adjusting for age at time of scan, associations were considered significant at P<0.05.

Regions of interest were then selected from the PLIC and OR to compliment the findings from tract-based spatial statistics, using a method previously described.(27) These regions encompass functionally relevant pathways and can be defined with high fidelity.(4) The left and right hemispheric measurements were averaged for each region, as previous.(4) To study the association between fatty acid levels and diffusion tensor imaging measures, linear mixed regression models were used to account for repeated measures. Adjustments were made for confounding factors selected a priori based on previous studies(4) and limited due to sample size limitations, including IVH severity, white matter injury (WMI) severity, patent ductus arteriosus, sepsis, and duration of intubation. Statistical significance was set at P<0.05.

Association between fatty acid level and outcome was explored using linear regression analysis to determine the long-term significance of fatty acid and imaging findings. Further adjustment for brain injury severity was performed to assess its contribution to the relationship. To understand how socioeconomic status may association between fatty acid level and outcomes, linear regression analysis was performed.