Abstract
Despite therapeutic hypothermia, survivors of neonatal encephalopathy have high rates of adverse outcome. Early surrogate outcome measures are needed to speed up the translation of neuroprotection trials. Thalamic lactate (Lac)/N-acetylaspartate (NAA) peak area ratio acquired with proton (1H) magnetic resonance spectroscopy (MRS) accurately predicts 2-year neurodevelopmental outcome. We assessed the relationship between MR biomarkers acquired at 24–48 h following injury with cell death and neuroinflammation in a piglet model following various neuroprotective interventions. Sixty-seven piglets with hypoxia–ischemia, hypoxia alone, or lipopolysaccharide (LPS) sensitization were included, and neuroprotective interventions were therapeutic hypothermia, melatonin, and magnesium. MRS and diffusion-weighted imaging (DWI) were acquired at 24 and 48 h. At 48 h, experiments were terminated, and immunohistochemistry was assessed. There was a correlation between Lac/NAA and overall cell death [terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL)] [mean Lac/NAA basal ganglia and thalamus (BGT) voxel r = 0.722, white matter (WM) voxel r = 0.784, p < 0.01] and microglial activation [ionized calcium-binding adapter molecule 1 (Iba1)] (BGT r = −0.786, WM r = −0.632, p < 0.01). Correlation with marker of caspase-dependent apoptosis [cleaved caspase 3 (CC3)] was lower (BGT r = −0.636, WM r = −0.495, p < 0.01). Relation between DWI and TUNEL was less robust (mean diffusivity BGT r = −0.615, fractional anisotropy BGT r = 0.523). Overall, Lac/NAA correlated best with cell death and microglial activation. These data align with clinical studies demonstrating Lac/NAA superiority as an outcome predictor in neonatal encephalopathy (NE) and support its use in preclinical and clinical neuroprotection studies.
Introduction
Neonatal encephalopathy (NE) secondary to intrapartum hypoxia–ischemia is a significant cause of brain injury in term infants affecting 2–3 per 1,000 live births in the UK (). Therapeutic hypothermia (HT) has reduced mortality and disability in survivors of NE [relative risk (RR) 0.75, 95% CI 0.68–0.83, number needed to treat (NNT) = 7] (). However, despite treatment, there remains a 24–30% mortality rate and 22–44% risk of moderate to severe disability at 18 months following moderate to severe NE (, ). NE has a complex and multifactorial etiology; however, over the last decade, preclinical () and clinical () studies suggest that coexisting infection and inflammation with hypoxia–ischemia (HI) exacerbate brain injury. A strong association exists between fetal infection/inflammation (e.g., chorioamnionitis, funisitis), perinatal brain damage, and neurodisability ().
In single () and multicenter () studies of NE babies who have been cooled, the 1H magnetic resonance spectroscopy (MRS) thalamic lactate (Lac)/N-acetylaspartate (NAA) peak area ratio acquired within 15 days of birth accurately predicts neurodevelopmental outcomes. Refinements in the spectral fitting including threonine (Thr) and N-acetylaspartylglutamate (NAAG) in the fitting function can improve the analysis of the spectrum in the regions close to Lac and NAA, respectively, and better signal to noise at 3 Tesla (3T) have optimized the predictive accuracy of Lac/NAA (). Using a threshold of 0.39, the sensitivity and specificity of BGT Lac/NAA for 2-year motor outcome was 100% and 97%, cognition 90% and 97% and language 81% and 97%, respectively (). In the TOBY Xenon early-phase clinical neuroprotection trial, adverse outcomes were correctly identified in 95.65% of cases by basal ganglia and thalamus (BGT) Lac/NAA, whereas prediction of adverse outcome using fractional anisotropy (FA) was 78.79% (). Using Lac/NAA peak area ratio as a qualified biomarker in the clinical context in a small proof-of-concept neuroprotection trial therefore avoids substantial financial and opportunity costs associated with large randomized controlled trials (RCTs).
Over the last two decades, we have used BGT and white matter (WM) Lac/NAA as one of our primary outcome markers in neuroprotection studies of adjunct therapies with HT in our piglet model (–). The piglet model allows for regional assessment of brain immunohistochemistry at 48 h with analyses including quantification of terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL)-positive cells, assessment of neuroinflammation [ionized calcium-binding adapter molecule 1 (Iba1) ramification index], and quantification of cleaved caspase 3 (CC3), a marker of caspase-dependent apoptosis.
Given the importance of MRI biomarkers in neonatal clinical neuroprotection trials and the translational pathway from preclinical to clinical RCTs, our aim was to assess: (i) the relationship between MR biomarkers [1H MRS metabolite ratios, mean diffusivity (MD), FA], acquired at 24 and 48 h following injury, and brain cell death and neuroinflammation at 48 h in the piglet following various neuroprotective interventions; (ii) brain immunohistochemistry differences related to the Lac/NAA peak area ratio clinical threshold of 0.39 (this ratio accurately predicts 2-year motor, cognitive, and language outcomes in babies with NE) (). In this study, we included retrospective data from different injuries (hypoxia–ischemia, hypoxia, inflammation-sensitization) and neuroprotective interventions (HT alone and with magnesium or melatonin) to assess the relation between MR biomarkers and immunohistochemistry in the piglet model.
Methods
Animal Experiments, Surgical Preparation, and Intensive Care Management
All animal experiments were approved by the UCL Ethics Committee and performed according to UK Home Office Guidelines [Animals (Scientific Procedures) Act, 1986]. The study complies with Animal Research: Reporting of in vivo Experiments (ARRIVE) guidelines.
Piglets were anesthetized and surgically prepared as described previously (, , , , ). In brief, all piglets were sedated with intramuscular midazolam and anesthetized with inhaled 3–4% v/v isoflurane. A tracheostomy was performed, and piglets were intubated (Smiths Medical, Ashford, Kent, UK) and ventilated (SLE 2000 Infant Ventilator, Surrey UK) for the duration of the experiment. Carotid vascular occluders (OC2A, in vivo Metric, Healdsburg, CA, USA) were sited for all studies, except study 3. Umbilical venous and arterial access were obtained (arterial catheter Vygon 2.5Fr, venous catheter−2Fr double lumen), and a peripherally inserted central venous catheter (Vygon 2Fr Nutriline) was sited in the proximal forelimb for infusion of intravenous drugs. Piglets were transferred onto a specialized incubator following surgery where continuous vital signs, multichannel electroencephalography (EEG) (Nicolet EEG, Natus), and cerebral near-infrared spectroscopy (NIRS) were monitored. Sedation was maintained with infusion of fentanyl (4 μg/kg/h) and inhaled isoflurane.
Piglets were cared for in accordance with local neonatal intensive care guidelines throughout the experiment. Following insult, maintenance fluid was restricted to 40 ml/kg/day. Ventilation settings were titrated according to arterial blood gas measurements. Mean arterial blood pressure (MABP) was maintained >35 mmHg using infusions of dopamine, dobutamine, noradrenaline, and adrenaline as required. Electrolytes, urea and creatinine, and blood glucose were monitored. All piglets received benzylpenicillin and gentamicin. 10% calcium gluconate (0.5 ml/kg) and salbutamol (4 μg/kg) were used to treat hyperkalemia. Seizures were treated with intravenous phenobarbitone followed by phenytoin if persistent.
Study Selection and Variations in Study Design
This study was a retrospective, secondary analysis of four preclinical neuroprotection piglet studies (, , , ). The study protocols evolved, reflecting optimization and development of study designs over the years, and are shown in Figure 1. For full details of the study methodology and results, please refer to publications (, , , ).
Figure 1
All studies lasted 48 h; however, studies varied according to brain injury protocols, duration of HT, and neuroprotective agents used. Primary outcome measures for all studies were identical; MRS was acquired at 24 and 48 h after insult, and immunohistochemistry was assessed at 48 h using the same methodology. Acquisition using the clinical 3T scanner (Philips Achieva) was introduced during study 1 to enhance the translational relevance of our preclinical model. Prior to this, 1H MRS was acquired using a 9.4T MRI scanner. Only piglets with MRS data at 3T and immunohistochemistry data were included in this secondary analysis. Piglets scanned at 9.4T or with no 3T MRS data were excluded.
Brain Injury
Piglets in studies 1 (
Piglets in studies 3 (
Neuroprotective Interventions
All piglets in studies 1 and 2 were cooled to 33.5°C using a servo-controlled water mattress (Tecotherm); however, protocols varied between the studies. In study 1, piglets were cooled from 2 h after HI over a duration of 24 h. Piglets in study 2 were cooled from 1 h after HI for a total duration for 12 h. No piglets in study 3 received HT. In study 4, piglets in the HT treatment arm were cooled for 12 h. All piglets that received HT were rewarmed at a controlled rate of 0.5°C/h to the target temperature of 38°C. Normothermia at 38.5°C was maintained by the water mattress.
Various neuroprotective agents were used in these studies. Piglets in study 1 received either an intravenous melatonin infusion at 2 and 26 h after HI at 5 mg/kg over 6 h or vehicle at the same volume and infusion rate. In study 2, piglets received (i) magnesium as a loading bolus of 180 mg/kg followed by continuous infusion 8 mg/kg/h at 1 h after HI; (ii) melatonin at 18 mg/kg over 2 h at 1 h and 25 h after HI; or (iii) vehicle at the same volume and rate. No additional agents were used in studies 3 or 4.
Magnetic Resonance Imaging
Piglets were transferred to the 3T MRI scanner at 24 and 48 h post insult. Imaging was performed with similar protocols as those used in NE babies on the same 3T scanner (
Figure 2

3T Magnetic Resonance Imaging. Piglets were scanned at 24 and 48 h in the Philips Achieva clinical magnet. 1H- Magnetic Resonance Spectroscopy (MRS) (A) was acquired using chemical shift imaging (CSI) with 8 × 8 matrix and 8 × 8 × 10 mm3 voxels, TR/TE was 2,000 ms/288 ms. The spectral width was 2 kHz with 2,048 points. Voxels selected were over the left thalamus [basal ganglia, thalamus (BGT)] and left subcortical white matter (WM) at the level of the centrum semiovale. Spectra were processed with Tarquin included in basis set, and the ratio of lactate (Lac)/N-acetylaspartate (NAA) was calculated from the amplitude of the fitted components (Lac+Thr/NAA+NAAG). For diffusion-weighted imaging (DWI) (B), regional of interest (ROI) for mean diffusivity (MD) and fractional anisotropy (FA) in the deep gray matter (DGM) (red) and WM (green) were automatically identified by atlas label propagation. The internal capsule was selected to represent WM, and the caudate, putamen, globus pallidus, thalamus, and hypothalamus were used for DGM.
Diffusion-weighted imaging (DWI) was acquired using a protocol similar to clinical studies (Figure 2B) (
Histology
Brain histology and immunohistochemistry samples were prepared as previously described (
TUNEL was used to assess cell death. As previously described (
For each piglet, eight regions of the brain were examined. In seven regions, TUNEL-positive nuclei were counted (Figure S1) from three fields in each of the R0 and R1 slices at 40× magnification. The hippocampus was present in the R1 section only. The counts were converted into cell counts per mm2.
To assess microglia activation, slides were prepared as previously described by Martinello et al. (
For CC3 immunohistochemistry, brain sections were pretreated as for Iba1 staining, incubated in rabbit anti-CC3 (1:100) (Cell Signaling 9661L) for 32 min followed by swine anti-rabbit immunoglobulin for 44 min. Sections were mounted on Vectrashield + DAPI as described above. CC3 cells were counted at 20× magnification in three fields per brain region and converted to counts per mm2.
Data and Statistical Analysis
Data analysis was carried out using SPSS Statistics 24 (IBM). The overall whole-brain TUNEL-positive cell counts, CC3-positive cell counts, and Iba1 ramification index were deduced from the sum of the average counts in eight regions of the brain (Figure S1).
The 24 and 48 h MRS data were collected for each of the BGT and WM regions and separately correlated with average whole-brain TUNEL, CC3, and Iba1 counts. In addition, the overall mean MRS FA and MD values were deduced from the 24 and 48 h scans and compared with the three immunohistochemistry markers.
The MRS, DWI, and histology count values were log10 transformed to normalize the distribution. The correlation was assessed using Pearson's rank coefficient, and scatter plots were created with GraphPad Prism v8 to illustrate the trend. P-values were calculated with two-tailed test to indicate statistical significance. As we compared multiple independent tests, the threshold for statistical significance was corrected to preserve a type 1 error rate (where p < 0.05 is significant) using Bonferonni correction. A p < 0.01 denotes statistical significance. Logistic regression modeling in infants with NE identified Lac/NAA of 0.39 as the optimal cutoff value for sensitivity and specificity to predict adverse neurodevelopmental outcomes at 18 months (
Results
Sixty-seven male large white piglets were recruited from four neuroprotection studies including 11 (16.4%) from Robertson et al. (
Figure 1 summarizes the treatment regimens across the four studies. In total, 59 of 67 (88.1%) piglets were subjected to cerebral injury. The remaining eight animals were either naive (n = 3, 4.5%) or LPS inflammation sensitized without cerebral injury (n = 5, 7.5%). Cerebral injury included: HI (39/67, 58.2%), LPS inflammation sensitized hypoxic injury (8/67, 11.9%), and LPS inflammation sensitized HI injury (12/67, 17.9%).
Neuroprotective interventions also varied with 45 of 67 (67.2%) piglets receiving HT for either 12 h (34/67, 50.7%) or 24 h (n = 11/67, 16.4%). Melatonin was given to 11 (16.4%) animals, and magnesium was given to eight (11.9%) piglets.
Lactate/N-Acetylaspartate to Terminal Deoxynucleotidyl Transferase dUTP Nick End Labeling
Mean Lac/NAA and TUNEL-positive cell counts of the whole brain correlated in the BGT (r = 0.722, p < 0.001) and WM voxel (r = 0.784, p < 0.001) (Figures 3A,B). The positive correlation was present at both 24 and 48 h (Table 1).
Figure 3

The relationship between thalamic [basal ganglia, thalamus (BGT)] and white matter (WM) 1H magnetic resonance spectroscopy (MRS) lactate (Lac)/N-acetylaspartate (NAA) and whole-brain cell death [terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL)] (A,B), microglia activation [ionized calcium-binding adapter molecule 1 (Iba1) ramification index] (C,D), and cleaved caspase 3 (CC3) (E,F). All data in the scatterplots were log10 transformed, and the Pearson's correlation coefficients (r-values) are illustrated. All values p < 0.001. Representative micrographs showing TUNEL, Iba1, and CC3 immunohistochemistry stains in piglets with Lac/NAA <0.39 (G–I) and Lac/NAA ≥0.39 (J–L) are shown.
Table 1
| MRS Lac/NAA | DWI Mean Diffusivity | DWI Fractional Anisotropy | ||||
|---|---|---|---|---|---|---|
| r | p-value | r | p-value | r | p-value | |
| 24 h BGT region | 0.709 | <0.001 | −0.488 | <0.001 | 0.424 | <0.001 |
| 48 h BGT region | 0.661 | <0.001 | −0.635 | <0.001 | 0.240 | 0.065 |
| Mean BGT region | 0.722 | <0.001 | −0.615 | <0.001 | 0.523 | <0.001 |
| 24 h WM region | 0.669 | <0.001 | −0.511 | <0.001 | 0.452 | <0.001 |
| 48 h WM region | 0.729 | <0.001 | −0.633 | <0.001 | 0.495 | <0.001 |
| Mean WM region | 0.748 | <0.001 | −0.635 | <0.001 | 0.342 | 0.005 |
Pearson's correlation coefficient comparing magnetic resonance spectroscopy (MRS) lactate (Lac)/N-acetylaspartate (NAA), diffusion-weighted imaging (DWI) mean diffusivity and fractional anisotropy at 24 and 48 h and the mean of the two time points with whole-brain terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL)-positive count.
All values were log10 transformed prior to correlation analysis. Correlation is significant if p < 0.01 (corrected with Bonferroni correction where p < 0.05 is considered significant). BGT, basal ganglia and thalamus; WM, white matter.
Lac/NAA ≥0.39 was associated with significantly higher TUNEL-positive cells in the whole brain across both voxels and both the 24 and 48 h scans (Table 2, Figures 3G,J) (p < 0.001) [mean BGT Lac/NAA ≥0.39, mean TUNEL count = 103 cells/mm2 (SD 2.12); mean BGT Lac/NAA <0.39, mean TUNEL count = 15.2 cells/mm2 (SD 2.56); p < 0.001].
Table 2
| Whole-Brain TUNEL Count | p-value | Whole-Brain Iba1 Ramification Index | p-value | Whole-Brain CC3 | p-value | ||||
|---|---|---|---|---|---|---|---|---|---|
| Lac/NAA <0.39 | Lac/NAA ≥0.39 | Lac/NAA <0.39 | Lac/NAA ≥0.39 | Lac/NAA <0.39 | Lac/NAA ≥0.39 | ||||
| BGT voxel at 24 h | 18.2 (SD 2.96) | 109 (SD 2.03) | <0.001 | 1.96 (SD 1.30) | 0.88 (SD 2.08) | <0.001 | 25.8 (SD 4.94) | 4.89 (SD 4.94) | <0.001 |
| BGT voxel at 48 h | 15.1 (SD 2.84) | 83.4 (SD 2.42) | <0.001 | 2.05 (SD 1.29) | 1.05 (SD 2.05) | <0.001 | 27.9 (SD 3.75) | 6.78 (SD 4.83) | 0.001 |
| Mean BGT voxel | 15.2 (SD 2.56) | 103 (SD 2.12) | <0.001 | 1.94 (SD 1.29) | 0.96 (SD 2.13) | <0.001 | 25.0 (SD 3.47) | 5.82 (SD 5.19) | 0.001 |
| WM voxel at 24 h | 15.9 (SD 3.26) | 84.7 (SD 2.27) | <0.001 | 1.87 (SD 1.61) | 1.11 (SD 1.96) | 0.004 | 24.2 (SD 4.76) | 7.71 (SD 4.39) | 0.008 |
| WM voxel at 48 h | 12.2 (SD 2.61) | 58.3 (SD 3.15) | <0.001 | 2.09 (SD 1.36) | 1.29 (SD 1.96) | 0.003 | 27.4 (SD 2.23) | 10.7 (SD 5.90) | 0.076 |
| Mean WM voxel | 11.3 (SD 2.27) | 71.6 (SD 2.75) | <0.001 | 2.04 (SD 1.29) | 1.17 (SD 1.99) | 0.003 | 33.0 (SD 2.37) | 8.08 (SD 5.41) | 0.001 |
Histology cell counts using a clinical threshold for lactate (Lac)/N-acetylaspartate (NAA) of 0.39.
Values shown are the geometric mean and standard deviation (SD). Using this threshold, we showed significant differences in cell counts for cell death [terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL)] and microglial activation [ionized calcium-binding adapter molecule 1 (Iba1) ramification index] in the piglets. Correlation is significant if p < 0.01 (corrected with Bonferroni correction where p < 0.05 is considered significant). BGT, basal ganglia and thalamus; WM, white matter.
Lactate/N-Acetylaspartate to Ionized Calcium-Binding Adapter Molecule 1
Iba1 ramification index was used to assess microglial activation. Activated microglia become ameboid with fewer processes, represented by a lower ramification index. There was a strong negative correlation between mean Lac/NAA and whole-brain Iba1 ramification index (Figures 3C,D). The negative correlation was strongest with the mean BGT voxel (r = −0.786, p < 0.001) but also present in the mean WM voxel (r = −0.632, p < 0.001).
Using a Lac/NAA threshold of 0.39, we noted significant differences in the Iba1 ramification index between piglets at all time points and voxels (p < 0.001) (Table 2, Figures 3H,K). Lac/NAA ≥0.39 was associated with lower Iba1 ramification (mean BGT Lac/NAA ≥0.39, Iba1 ramification index 0.96 vs. 1.94 with Lac/NAA <0.39) (p < 0.001).
Lactate/N-Acetylaspartate to Cleaved Caspase 3
The correlation between Lac/NAA and CC3 was negative in the BGT voxel (r = −0.636; p < 0.001) but was weaker in the WM voxel (r = −0.495; p < 0.001) (Figures 3E,F).
Lac/NAA ≥0.39 was associated with lower CC3 counts at 24 h (CC3 count 4.89 vs. 25.8 cells/mm2, p < 0.001) and 48 h scans (CC3 count 6.78 vs. 27.9 cells/mm2, p < 0.01) in the BGT voxel. The CC3 count was also significantly lower with Lac/NAA ≥0.39 at 24 h in the WM voxel (CC3 count 7.71 vs. 25.0, p < 0.01) (Table 2, Figures 3I,L).
Magnetic Resonance Spectroscopy Metabolite Peak Ratios to Terminal Deoxynucleotidyl Transferase dUTP Nick End Labeling
Pearson's correlation coefficients were deduced to compare other proton MRS metabolite ratios with whole-brain TUNEL count (Table 3). We observed strong positive correlations between the total whole-brain TUNEL count and mean BGT Lac/Cho (r = 0.765, p < 0.001) and BGT Lac/Cr (r = 0.765, p < 0.001). There was a lesser correlation in the corresponding WM voxels (WM Lac/Cho r = 0.701, p < 0.001; Lac/Cr 0.671, p < 0.001). There was a weak correlation between BGT NAA/Cho (r = −0.530, p < 0.01) and BGT NAA/Cr (r = −0.565, p < 0.001) with TUNEL. There was no correlation between TUNEL and Cho/Cr (BGT voxel, r = 0.019, p = 0.88; WM voxel, r = −0.051, p = 0.68).
Table 3
| r | p-value | ||
|---|---|---|---|
| Mean Lac/NAA | BGT | 0.722 | <0.001 |
| WM | 0.748 | <0.001 | |
| Mean Lac/Cho | BGT | 0.765 | <0.001 |
| WM | 0.701 | <0.001 | |
| Mean Lac/Cr | BGT | 0.766 | <0.001 |
| WM | 0.671 | <0.001 | |
| Mean NAA/Cho | BGT | −0.530 | <0.001 |
| WM | −0.565 | <0.001 | |
| Mean NAA/Cr | BGT | −0.565 | <0.001 |
| WM | −0.624 | <0.001 | |
| Mean Cho/Cr | BGT | 0.019 | 0.877 |
| WM | −0.051 | 0.684 |
Pearson's correlation coefficient comparing magnetic resonance spectroscopy (MRS) metabolite ratios with whole-brain terminal deoxynucleotidyl transferase dUTP nick end labeling.
Correlation is significant if p < 0.01 (corrected with Bonferroni correction where p < 0.05 is considered significant). BGT, basal ganglia and thalamus; Cho, choline; Cr, creatine; Lac, lactate; NAA, N-acetylaspartate; WM, white matter.
When compared with other metabolic ratios Lac/NAA was most consistent in yielding similarly strong correlation co-efficient values in the BGT and WM voxels.
Diffusion-Weighted Imaging to Terminal Deoxynucleotidyl Transferase dUTP Nick End Labeling
The correlation between DWI MD and TUNEL was negative (mean DGM MD to TUNEL r = 0.615, p < 0.001; mean WM MD to TUNEL r = −0.635, p < 0.001) as illustrated in Figures 4A,B; however, the correlation was not as a strong as between Lac/NAA and TUNEL-positive cells (Table 3). The correlation between FA and TUNEL was weak at 24 and 48 h (mean DGM r = 0.523, p < 0.001; mean WM r = 0.342, p < 0.01) (Figures 4C,D; Table 1). Representative T2-weighted images (T2W), MD and FA maps are shown in Figures 4E–J.
Figure 4

Relationship between diffusion-weighted imaging (DWI) and whole-brain cell death [terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL)-positive cells]. Scatterplots showing the correlation between whole-brain TUNEL and average mean diffusivity (MD) (A,B) and fractional anisotropy (FA) (C,D) localized in the deep gray matter (DGM), white matter (WM) regions. Data were log10 transformed, and the Pearson's correlation coefficients (r-values) are illustrated. *p < 0.01, **p < 0.001. Representative T2-weighted images (T2W), MD and FA maps are shown for a piglet with low TUNEL-positive counts (E–G, respectively) and high TUNEL-positive counts (H–J).
Discussion
1H MRS Lac/NAA peak area ratio correlated with overall TUNEL-positive cell death and microglial activation in a piglet model of term perinatal brain injury. Compared to other MRS metabolite ratios and DWI, Lac/NAA demonstrated the best correlation to TUNEL-positive cell death at 24 and 48 h. This concurs with studies in cooled infants with NE where BGT Lac/NAA peak area ratio has a high level of accuracy for outcome prediction compared to other MR methods (
The combined increased lactate and reduced NAA on MRS (translating to a high Lac/NAA peak area ratio) suggest brain mitochondrial impairment and impaired oxidative metabolism during “secondary energy failure.” It is possible that this ongoing injury may be amenable to late therapies and thus is an important measure to direct therapies. Woo et al. (
Lac/NAA peak ratio in our piglet studies represents more precisely Lac+threonine/total NAA. Mitra et al. (
The source and exact mechanism of raised cerebral lactate remain unknown. Both neurons and astrocytes produce lactate in hypoxic conditions in vitro (
There is increasing interest in the role of lactate as a neuroprotective agent. In adults with traumatic brain injury, administration of 13C-labeled lactate via the microdialysis catheter and simultaneous collection of the microdialysates, with 13C NMR analysis, revealed 13C labeling in glutamine, consistent with lactate metabolism in the TCA cycle (
NAA is a metabolite produced by aspartate N-acetyl transferase and found in neurons. NAA is transported from neurons to oligodendrocytes where it is metabolized into aspartate and acetate and used for energy production and myelin synthesis, respectively (
We noted a negative correlation between DWI MD and TUNEL-positive cells in the brain. In babies with NE, MD pseudonormalizes at around 7 days in non-cooled infants or 10 days in infants who received HT (
In our experience, although highest levels of Lac are seen in the first few days after birth, brain Lac persists for months in babies with adverse outcome after NE; this persisting brain Lac is associated with abnormal MRI and brain alkalosis (
Microglia are one of the first inflammatory cells to be activated following HI (
We observed a poor correlation between Lac/NAA and CC3, and this relationship was unexpectedly negative. This finding concurs with our previous neuroprotection studies where we see little relation between CC3 and other markers of injury severity, particularly TUNEL-positive cells (
There are limitations to this study. These data were retrospective and obtained from different studies with differing insults and neuroprotective interventions. However, this is also a strength of the study as the strong correlation of Lac/NAA to TUNEL-positive cells supports the validity of this biomarker across perinatal brain injury which is frequently multifactorial and heterogeneous in nature. In this study, some animals were cooled for 12 h, rather than 24 h (
In conclusion, we describe a strong correlation between MRS Lac/NAA and TUNEL-positive cells and microglial activation across WM and gray matter in male piglets after a range of perinatal insults and neuroprotective interventions. These preclinical data concur with clinical studies that have demonstrated the utility of BGT Lac/NAA as a surrogate marker that best predicts outcome in NE and can be used to expedite early-phase clinical neuroprotection trials in NE.
Statements
Data availability statement
The datasets generated for this study are available on request to the corresponding author.
Ethics statement
The animal study was reviewed and approved by UCL Ethics Committee.
Author contributions
RP organized and analyzed the data and drafted the manuscript with the help of KM. KM, AA-B, CM, IL, and TM undertook the experiments. CM undertook microscopy and cell quantification and organized the histology results. MS and AB scanned the piglets and collected MRS and DWI data. XG assisted with MR physics aspects of the study and reviewed the manuscript. NR obtained funding for the studies, designed the studies, and reviewed the manuscript. All authors reviewed the manuscript and approved the final version as submitted and agreed to be accountable for all aspects of the work.
Funding
These studies were funded by the Medical Research Council MR/M006743/1, Chiesi Pharmaceuticals (research grant), and Action Medical Research for Children (GN2295).
Acknowledgments
We thank Debbie Kraus for statistical support and Mariya Hristova for her assistance with histology preparation and fluorescent microscopy. This work was undertaken at University College London Hospitals/University College London, which received a proportion of funding from the UK Department of Health's National Institute for Health Research Biomedical Research Centres funding scheme.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2020.00883/full#supplementary-material
- 3T
3-Tesla
- ADC
apparent diffusion coefficient
- ANLS
astrocyte neuron lactate shuttle
- AUC
area under the curve
- BGT, basal ganglia
thalamus
- CC3
cleaved caspase 3
- Cho
choline
- Cr
creatine
- DGM
Deep gray matter
- DWI
diffusion-weighted imaging
- FA
fractional anisotropy
- HI
hypoxia–ischemia
- HT
therapeutic hypothermia
- Iba1
ionized calcium-binding adapter molecule 1
- Lac/NAA
lactate/N-acetylaspartate
- Lac
lactate
- LPS
lipopolysaccharide
- MD
mean diffusivity
- MRS
magnetic resonance spectroscopy
- NAA
N-acetylaspartate
- NAAG
N-acetylaspartylglutamate
- NAD
nicotinamide adenine dinucleotide
- NE
neonatal encephalopathy
- RCT
randomized controlled trial
- ROI
region of interest
- T2W
T2 weighted
- Thr
threonine
- TUNEL
terminal deoxynucleotidyl transferase dUTP nick end labeling
- WM
white matter.
Abbreviations
References
1.
GaleCStatnikovYJawadSUthayaSNModiNGroupBIew. Neonatal brain injuries in England: population-based incidence derived from routinely recorded clinical data held in the National Neonatal Research Database. Arch Dis Child Fetal Neonatal Ed. (2018). 103:F301–6. 10.1136/archdischild-2017-313707
2.
JacobsSEBergMHuntRTarnow-MordiWOInderTEDavisPG. Cooling for newborns with hypoxic ischaemic encephalopathy. Cochrane Database Syst Rev. (2013) 2013:CD003311. 10.1002/14651858.CD003311.pub3
3.
AzzopardiDStrohmBMarlowNBrocklehurstPDeierlAEddamaOet al. Effects of hypothermia for perinatal asphyxia on childhood outcomes. N Engl J Med. (2014) 371:140–9. 10.1056/NEJMoa1315788
4.
ShankaranSLaptookAREhrenkranzRATysonJEMcDonaldSADonovanEFet al. Whole-body hypothermia for neonates with hypoxic-ischemic encephalopathy. N Engl J Med. (2005) 353:1574–84. 10.1056/NEJMcps050929
5.
EklindSMallardCLeverinALGillandEBlomgrenKMattsby-BaltzerIet al. Bacterial endotoxin sensitizes the immature brain to hypoxic–ischaemic injury. Eur J Neurosci. (2001) 13:1101–6. 10.1046/j.0953-816x.2001.01474.x
6.
NelsonKBWilloughbyRE. Infection, inflammation and the risk of cerebral palsy. Curr Opin Neurol. (2000) 13:133–9. 10.1097/00019052-200004000-00004
7.
TannCJNakakeetoMWilleyBASewegabaMWebbELOkeIet al. Perinatal risk factors for neonatal encephalopathy: an unmatched case-control study. Arch Dis Child Fetal Neonatal Ed. (2018) 103:F250–6. 10.1136/archdischild-2017-312744
8.
MitraSKendallGSBainbridgeASokolskaMDinanMUria-AvellanalCet al. Proton magnetic resonance spectroscopy lactate/N-acetylaspartate within 2 weeks of birth accurately predicts 2-year motor, cognitive and language outcomes in neonatal encephalopathy after therapeutic hypothermia. Arch Dis Child Fetal Neonatal Ed. (2019) 104:F424–32. 10.1136/archdischild-2018-315478
9.
LallyPJMontaldoPOliveiraVSoeASwamyRBassettPet al. Magnetic resonance spectroscopy assessment of brain injury after moderate hypothermia in neonatal encephalopathy: a prospective multicentre cohort study. Lancet Neurol. (2019) 18:35–45. 10.1016/S1474-4422(18)30325-9
10.
AzzopardiDChewATDeierlAHuertasARobertsonNJTusorNet al. Prospective qualification of early cerebral biomarkers in a randomised trial of treatment with xenon combined with moderate hypothermia after birth asphyxia. EBioMedicine. (2019) 47:484–91. 10.1016/j.ebiom.2019.08.034
11.
RobertsonNJMartinelloKLingamIAvdic-BelltheusAMeehanCAlonso-AlconadaDet al. Melatonin as an adjunct to therapeutic hypothermia in a piglet model of neonatal encephalopathy: a translational study. Neurobiol Dis. (2018) 121:240–51. 10.1016/j.nbd.2018.10.004
12.
RobertsonNJFaulknerSFleissBBainbridgeAAndorkaCPriceDet al. Melatonin augments hypothermic neuroprotection in a perinatal asphyxia model. Brain. (2013) 136:90–105. 10.1093/brain/aws285
13.
LingamIMeehanCAvdic-BelltheusAMartinelloKHristovaMKaynezhadPet al. Short-term effects of early initiation of magnesium infusion combined with cooling after hypoxia-ischemia in term piglets. Pediatr Res. (2019) 86:699–708. 10.1038/s41390-019-0511-8
14.
MartinelloKAMeehanCAvdic-BelltheusALingamIRagabSHristovaMet al. Acute LPS sensitization and continuous infusion exacerbates hypoxic brain injury in a piglet model of neonatal encephalopathy. Sci Rep. (2019) 9:10184. 10.1038/s41598-019-46488-y
15.
BroadKDFierensIFleissBRocha-FerreiraEEzzatiMHassellJet al. Inhaled 45-50% argon augments hypothermic brain protection in a piglet model of perinatal asphyxia. Neurobiol Dis. (2016) 87:29–38. 10.1016/j.nbd.2015.12.001
16.
LorekATakeiYCadyEWyattJPenriceJEdwardsAet al. Delayed (“secondary”) cerebral energy failure after acute hypoxia-ischemia in the newborn piglet: continuous 48-hour studies by phosphorus magnetic resonance spectroscopy. Pediatr Res. (1994) 36:699–706. 10.1203/00006450-199412000-00003
17.
RobertsonNJLingamIMeehanCMartinelloKAAvdic-BelltheusASteinLet al. High-dose melatonin and ethanol excipient combined with therapeutic hypothermia in a newborn piglet asphyxia model. Sci Rep. (2020) 10:3898. 10.1038/s41598-020-60858-x
18.
YushkevichPAPivenJHazlettHCSmithRGHoSGeeJCet al. User-guided 3D active contour segmentation of anatomical structures: significantly improved efficiency and reliability. Neuroimage. (2006) 31:1116–28. 10.1016/j.neuroimage.2006.01.015
19.
ConradMSSuttonBPDilgerRNJohnsonRW. An in vivo three-dimensional magnetic resonance imaging-based averaged brain collection of the neonatal piglet (Sus scrofa). PLoS ONE. (2014) 9:e107650. 10.1371/journal.pone.0107650
20.
ItoDTanakaKSuzukiSDemboTFukuuchiY. Enhanced expression of Iba1, ionized calcium-binding adapter molecule 1, after transient focal cerebral ischemia in rat brain. Stroke. (2001) 32:1208–15. 10.1161/01.STR.32.5.1208
21.
ThayyilSChandrasekaranMTaylorABainbridgeACadyEBChongWKet al. Cerebral magnetic resonance biomarkers in neonatal encephalopathy: a meta-analysis. Pediatrics. (2010) 125:e382–95. 10.1542/peds.2009-1046
22.
AlderliestenTde VriesLSStaatsLvan HaastertICWeekeLBendersMJet al. MRI and spectroscopy in (near) term neonates with perinatal asphyxia and therapeutic hypothermia. Arch Dis Child Fetal Neonatal Ed. (2017) 102:F147–52. 10.1136/archdischild-2016-310514
23.
ShibasakiJAidaNMorisakiNTomiyasuMNishiYToyoshimaK. Changes in brain metabolite concentrations after neonatal hypoxic-ischemic encephalopathy. Radiology. (2018) 288:840–8. 10.1148/radiol.2018172083
24.
WooCWLeeBSKimSTKimKS. Correlation between lactate and neuronal cell damage in the rat brain after focal ischemia: an in vivo 1H magnetic resonance spectroscopic (1H-MRS) study. Acta Radiol. (2010) 51:344–50. 10.3109/02841850903515395
25.
AridasJDYawnoTSutherlandAENitsosIDitchfieldMWongFYet al. Detecting brain injury in neonatal hypoxic ischemic encephalopathy: closing the gap between experimental and clinical research. Exp Neurol. (2014) 261:281–90. 10.1016/j.expneurol.2014.07.009
26.
ChoiCCouplandNJKalraSBhardwajPPMalykhinNAllenPS. Proton spectral editing for discrimination of lactate and threonine 1.31 ppm resonances in human brain in vivo. Magn Reson Med. (2006) 56:660–5. 10.1002/mrm.20988
27.
SonnewaldUWangAYSchousboeAEriksonRSkottnerA. New aspects of lactate metabolism: IGF-I and insulin regulate mitochondrial function in cultured brain cells during normoxia and hypoxia. Dev Neurosci. (1996) 18:443–8. 10.1159/000111439
28.
StovellMGMadaMOHelmyACarpenterTAThelinEPYanJLet al. The effect of succinate on brain NADH/NAD. Sci Rep. (2018) 8:11140. 10.1038/s41598-018-29255-3
29.
CarpenterKLJallohIHutchinsonPJ. Glycolysis and the significance of lactate in traumatic brain injury. Front Neurosci. (2015) 9:112. 10.3389/fnins.2015.00112
30.
PellerinLMagistrettiPJ. Glutamate uptake into astrocytes stimulates aerobic glycolysis: a mechanism coupling neuronal activity to glucose utilization. Proc Natl Acad Sci USA. (1994) 91:10625–9. 10.1073/pnas.91.22.10625
31.
LamaSAuerRNTysonRGallagherCNTomanekBSutherlandGR. Lactate storm marks cerebral metabolism following brain trauma. J Biol Chem. (2014) 289:20200–8. 10.1074/jbc.M114.570978
32.
RobertsonNJCoxIJCowanFMCounsellSJAzzopardiDEdwardsAD. Cerebral intracellular lactic alkalosis persisting months after neonatal encephalopathy measured by magnetic resonance spectroscopy. Pediatr Res. (1999) 46:287–96. 10.1203/00006450-199909000-00007
33.
Uria-AvellanalCRobertsonNJ. Na+/H+ exchangers and intracellular pH in perinatal brain injury. Transl Stroke Res. (2014) 5:79–98. 10.1007/s12975-013-0322-x
34.
WuTWTamraziBHsuKHHoEReitmanAJBorzageMet al. Cerebral lactate concentration in neonatal hypoxic-ischemic encephalopathy: in relation to time, characteristic of injury, and serum lactate concentration. Front Neurol. (2018) 9:293. 10.3389/fneur.2018.00293
35.
CarpenterKLJallohIGallagherCNGricePHoweDJMasonAet al. (13)C-labelled microdialysis studies of cerebral metabolism in TBI patients. Eur J Pharm Sci. (2014) 57:87–97. 10.1016/j.ejps.2013.12.012
36.
RoumesHDumontUSanchezSMazuelLBlancJRaffardGet al. Neuroprotective role of lactate in rat neonatal hypoxia-ischemia. J Cereb Blood Flow Metab. (2020). 10.1177/0271678X20908355. [Epub ahead of print].
37.
RoskoLSmithVNYamazakiRHuangJK. Oligodendrocyte bioenergetics in health and disease. Neuroscientist. (2018) 25:334–43. 10.1177/1073858418793077
38.
SchmitzBWangXBarkerPBPilatusUBronzlikPDadakMet al. Effects of aging on the human brain: a proton and phosphorus MR spectroscopy study at 3T. J Neuroimaging. (2018) 28:416–21. 10.1111/jon.12514
39.
CheongJLCadyEBPenriceJWyattJSCoxIJRobertsonNJ. Proton MR spectroscopy in neonates with perinatal cerebral hypoxic-ischemic injury: metabolite peak-area ratios, relaxation times, and absolute concentrations. Am J Neuroradiol. (2006) 27:1546–54.
40.
SijensPEWischniowskyKTer HorstHJ. The prognostic value of proton magnetic resonance spectroscopy in term newborns treated with therapeutic hypothermia following asphyxia. Magn Reson Imaging. (2017) 42:82–7. 10.1016/j.mri.2017.06.001
41.
BednarekNMathurAInderTWilkinsonJNeilJShimonyJ. Impact of therapeutic hypothermia on MRI diffusion changes in neonatal encephalopathy. Neurology. (2012) 78:1420–7. 10.1212/WNL.0b013e318253d589
42.
IgarashiHSuzukiYHuberVJIdaMNakadaT. N-acetylaspartate decrease in acute stage of ischemic stroke: a perspective from experimental and clinical studies. Magn Reson Med Sci. (2015) 14:13–24. 10.2463/mrms.2014-0039
43.
AncoraGTestaCGrandiSTononCSbravatiFSaviniSet al. Prognostic value of brain proton MR spectroscopy and diffusion tensor imaging in newborns with hypoxic-ischemic encephalopathy treated by brain cooling. Neuroradiology. (2013) 55:1017–25. 10.1007/s00234-013-1202-5
44.
RobertsonNJCowanFMCoxIJEdwardsAD. Brain alkaline intracellular pH after neonatal encephalopathy. Ann Neurol. (2002) 52:732–42. 10.1002/ana.10365
45.
MillarLJShiLHoerder-SuabedissenAMolnárZ. Neonatal hypoxia ischaemia: mechanisms, models, and therapeutic challenges. Front Cell Neurosci. (2017) 11:78. 10.3389/fncel.2017.00078
46.
Rocha-FerreiraEHristovaM. Antimicrobial peptides and complement in neonatal hypoxia-ischemia induced brain damage. Front Immunol. (2015) 6:56. 10.3389/fimmu.2015.00056
47.
FleissBGressensP. Tertiary mechanisms of brain damage: a new hope for treatment of cerebral palsy?Lancet Neurol. (2012) 11:556–66. 10.1016/S1474-4422(12)70058-3
48.
Charriaut-MarlangueCBessonVCBaudO. Sexually dimorphic outcomes after neonatal stroke and hypoxia-ischemia. Int J Mol Sci. (2017) 19:61. 10.3390/ijms19010061
49.
BurguillosMADeierborgTKavanaghEPerssonAHajjiNGarcia-QuintanillaAet al. Caspase signalling controls microglia activation and neurotoxicity. Nature. (2011) 472:319–24. 10.1038/nature09788
50.
VillapolSBonninPFauSBaudORenolleauSCharriaut-MarlangueC. Unilateral blood flow decrease induces bilateral and symmetric responses in the immature brain. Am J Pathol. (2009) 175:2111–20. 10.2353/ajpath.2009.090257
51.
AbrahamMCShahamS. Death without caspases, caspases without death. Trends Cell Biol. (2004) 14:184–93. 10.1016/j.tcb.2004.03.002
52.
ZhouTLinHJiangLYuTZengCLiuJet al. Mild hypothermia protects hippocampal neurons from oxygen-glucose deprivation injury through inhibiting caspase-3 activation. Cryobiology. (2018) 80:55–61. 10.1016/j.cryobiol.2017.12.004
53.
Hellström-WestasLRosénISvenningsenNW. Predictive value of early continuous amplitude integrated EEG recordings on outcome after severe birth asphyxia in full term infants. Arch Dis Child Fetal Neonatal Ed. (1995) 72:F34–8. 10.1136/fn.72.1.F34
Summary
Keywords
neonatal encephalopathy, magnetic resonance spectroscopy, hypoxia–ischemia, piglet, therapeutic hypothermia, neuroprotection
Citation
Pang R, Martinello KA, Meehan C, Avdic-Belltheus A, Lingam I, Sokolska M, Mutshiya T, Bainbridge A, Golay X and Robertson NJ (2020) Proton Magnetic Resonance Spectroscopy Lactate/N-Acetylaspartate Within 48 h Predicts Cell Death Following Varied Neuroprotective Interventions in a Piglet Model of Hypoxia–Ischemia With and Without Inflammation-Sensitization. Front. Neurol. 11:883. doi: 10.3389/fneur.2020.00883
Received
21 January 2020
Accepted
10 July 2020
Published
04 September 2020
Volume
11 - 2020
Edited by
Deirdre M. Murray, University College Cork, Ireland
Reviewed by
Yohan van de Looij, Université de Genève, Switzerland; Diego Iacono, Biomedical Research Institute of New Jersey, United States
Updates

Check for updates
Copyright
© 2020 Pang, Martinello, Meehan, Avdic-Belltheus, Lingam, Sokolska, Mutshiya, Bainbridge, Golay and Robertson.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nicola J. Robertson n.robertson@ucl.ac.uk
This article was submitted to Pediatric Neurology, a section of the journal Frontiers in Neurology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.