- 1Department of Pediatrics, Division of Neonatology, Johns Hopkins University School of Medicine, Baltimore, MD, United States
- 2Neuroscience Intensive Care Nursery Program, Division of Neonatology, Johns Hopkins University School of Medicine, Baltimore, MD, United States
- 3Department of Internal Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, United States
- 4Department of Pediatrics, Division of Neonatology, Johns Hopkins All Children's Hospital, St Petersburg, FL, United States
- 5Department of Neuropsychology, Kennedy Krieger Institute, Baltimore, MD, United States
- 6Department of Psychiatry and Behavioral Sciences, Johns Hopkins University School of Medicine, Baltimore, MD, United States
- 7Department of Radiology, Division of Pediatric Radiology and Pediatric Neuroradiology, Johns Hopkins Hospital, Baltimore, MD, United States
- 8Department of Obstetrics and Gynecology, Johns Hopkins University School of Medicine, Baltimore, MD, United States
- 9Department of Pediatrics, Division of Pediatric Cardiology, Johns Hopkins University School of Medicine, Baltimore, MD, United States
Objective: To determine the changes due to therapeutic hypothermia (TH) exposure in the strength of association between traditional clinical and biochemical indicators of severity of neonatal hypoxic-ischemic encephalopathy (HIE) and serum biomarkers. We hypothesized that culmination of TH changes the strength of the relationships between traditional indicators of severity of HIE and serum biomarkers.
Methods: This was a single-center observational cohort study of 178 neonates with HIE treated with TH and followed with serum biomarkers: (i) brain-derived neurotrophic factor (BDNF) and vascular endothelial growth factor (VEGF) (neurotrophins); (ii) tau and glial fibrillary acidic protein (GFAP) (neural cell injury); and (iii) interleukin 6 (IL-6), IL-8, and IL-10 (cytokines), during their first week of life. Adjusted mixed-effect models tested associations with HIE indicators in relation to TH exposure.
Results: At admission, lower Apgar scores and base excess (BE) and higher lactate and nucleated red blood cell (NRBC) count correlated with higher Sarnat scores. These indicators of worse HIE severity, including higher Sarnat score, correlated with lower VEGF and higher tau, GFAP, and IL-10 levels at different time points. Within the first 24 h of life, patients with a Sarnat score >2 had lower VEGF levels, whereas only those with score of 3 also had higher GFAP and IL-10 levels. Tau levels increased during TH in patients with Sarnat score of 3, whereas tau and GFAP increased after TH in those with scores of 2. After adjustments, lower VEGF levels during TH and higher tau, GFAP, and IL-10 levels during and after TH were associated with worse Sarnat scores. Tau and GFAP relationship with Sarnat score became stronger after TH.
Conclusion: Therapeutic hypothermia exerts an independent modulatory effect in the relationships between traditional indicators of severity of HIE and serum biomarkers after adjustments. Thus, the timing of biomarker testing in relation to TH exposure must be carefully considered if biomarkers are proposed for patient stratification in novel clinical trials.
Introduction
Hypoxic-ischemic (HI) encephalopathy (HIE) is the most prevalent type of brain injury in full-term neonates and a main cause of neonatal encephalopathy (NE) (1). Impaired fetal perfusion results in biochemical changes, such as worsening acidosis and ultimately HI injury to the brain (1, 2). Hypoxic-ischemic encephalopathy results from a cascade of excitotoxicity, oxidative stress, and inflammation, which may persist for weeks (3–5). The severity and recovery of HI brain injury do not always align with early assessments using traditional clinical and biochemical indicators because of lack of sensitivity and specificity to brain injury. Therapeutic hypothermia (TH), the only available therapy, reduces death or disability in patients with moderate HIE (6–9). However, sex and other factors (5, 10–13) may play a role in the variable post-natal presentation and response to TH. The lack of specificity of the existing criteria guiding the initiation of TH offers an unsatisfactory assessment of severity of brain injury hindering our ability to identify patients who may benefit of adjuvant therapeutic strategies and to monitor response to therapy (1, 14). Thus, peripheral blood biomarkers may improve the assessment provided by traditional indicators alone.
Brain insults trigger activation of signaling cascades, resulting in a temporally dynamic “biochemical footprint” detectable in circulating blood that have the potential to be used as biomarkers. These biomarkers may improve the ability to grade severity of injury and guide therapy in several neurological conditions (15–17). In neonatal HIE, brain-derived neurotrophic factor (BDNF), glial fibrillary acidic protein (GFAP), and tau have been among those studied against modified Sarnat scores, as well as brain magnetic resonance imaging (MRI) and neurodevelopmental outcomes (18, 19). However, it is unclear how longitudinal changes in these and other proposed HIE biomarkers relate to traditional early clinical and biochemical indicators in the era of TH. Understanding these relationships and how TH exposure may modulate them will inform the optimal timing at which proposed biomarkers may indicate the specific type and severity of HI injury to guide treatment and refine predictive models.
Here, we aim to study the changes attributable to TH exposure in the strength of association between traditional clinical and biochemical indicators of severity of neonatal HIE and serum biomarkers. These biomarkers included (i) BDNF and vascular endothelial growth factor (VEGF), neurotrophins essential for brain development (18, 20, 21); (ii) tau and GFAP, neuronal (19, 22, 23), and astrocytic (24) cytoskeletal markers linked to brain cell injury; and (iii) interleukin 6 (IL-6), IL-8, and IL-10, markers of inflammatory response (18, 20, 21); most of them studied in association to neonatal HIE previously in smaller cohorts. We hypothesized that culmination of TH changes the strength of the relationship between traditional indicators of severity of HIE and serum biomarkers. Here, we report that (i) indicators of worse HIE severity, including higher Sarnat score, correlated with lower VEGF and higher tau, GFAP, and IL-10 levels; and (ii) lower VEGF levels during TH, and higher tau, GFAP, and IL-10 levels during and after TH, were associated with worse Sarnat scores with tau and GFAP relationship with Sarnat score becoming stronger after TH.
Methods
The study received exempt status from the Johns Hopkins University (JHU) institutional review board (IRB), until 2017, when informed consent became required to access medical records and use discarded blood collected for clinical purposes. The Neuroscience Intensive Care Nursery program coordinator (C.P.) identified patients with diagnosis of NE and obtained agreement by the treating clinical team to allow a member of the study team to discuss inclusion of the neonate in the study with the parents. Written consent was obtained from the parents of the participants after 2017, as required by IRB. The study Identification of Diagnostic and Prognostic Biomarkers for Perinatal Hypoxic-Ischemic Brain Injury (BIN study, IRB NA_00026068) was in compliance with the Health Insurance Portability and Accountability Act of 1996 (HIPAA).
Patients
Inclusion Criteria
A prospective cohort of 246 neonates (GA ≥35 weeks) diagnosed with NE, who were initiated in whole-body TH within 6 h of life between April 28, 2009, and November 15, 2019, at a single level IV neonatal intensive care unit (NICU).
Exclusion Criteria
Forty-five patients were excluded for incomplete clinical data, off-label use of TH (<35 weeks GA), partial TH course, non-peripartum events (apparent life-threatening event/brief, resolved, unexplained event), death for causes other than HIE (i.e., severe pulmonary hypertension of the newborn), need for extracorporeal membrane oxygenation, and non-HIE causes of NE (Figure 1). Of 201 patients, 178 had 7 days of biomarker data. Diagnosis of HIE was based on the National Institute of Child Health and Human Development (NICHD) Neonatal Research Network criteria (6).
Whole-Body TH Protocol
Initiation of TH occurred within the first 6 h of life in infants ≥1,800 g birth weight and ≥35 weeks GA with (a) cord or first hour of life pH ≤ 7.0 and/or base excess (BE) ≤-16 or (b) cord or first hour pH > 7.0 but ≤7.15 and/or BE between −10 and −16, plus evidence of a perinatal adverse event and need for assisted ventilation or Apgar <5 at 10 min of life. Although moderate to severe encephalopathy is a criterion used to initiate TH, ~20% of patients included in the analysis had a modified Sarnat score (25, 26) of 1 due to concerns that the neurological examination may worsen after 6 h of life (27–29), particularly in those with severe metabolic acidosis (30, 31). In our cohort, 83% of patients with Sarnat score of 1 who received TH had severe metabolic acidosis. The inclusion of these patients in the analysis better reflects the current state of practice.
Clinical Data
Clinical data were obtained from electronic charts (Table 1). Race was assigned based on maternal race. Absent/reversed end-diastolic flow in umbilical artery was determined by Doppler scan of the umbilical artery. Placenta pathology was used to assign histological chorioamnionitis or funisitis. Assignment of sex of the newborn, confirmation of GA, and Apgar scores were determined by the NICU team during initial assessment. Highest modified Sarnat score during the first 6 h of admission to the NICU was determined by the study team (R.C.V., C.P., and F.J.N.) (25, 26). Most deaths occurring in the cohort were secondary to withdrawal of care in conjunction with the family after discussion about clinical prognosis.
Laboratory Data
Results of the blood gas with lactate (mmol/L) within 1 h of life were adjusted to core body temperature. Lowest (worse) pH and BE measures were used for analysis, as that better represents clinical practice. Biochemistry and hematology data are reported in Table 1. Clinical seizures were confirmed during the first 72 h of life using continuous full-montage video electroencephalogram with interpretation by pediatric neurology. Infection was assigned by clinical or histopathological chorioamnionitis or by proven neonatal sepsis confirmed by positive blood, cerebrospinal fluid, or urine culture with need for >72 h of antibiotics as a surrogate of clinically significant infection.
Perinatal Insult Score
The Perinatal Insult (PI) score was created by assigning a hierarchical value to the NICHD criteria used to decide initiation of TH to treat HIE. The PI score ranges from 0 to 9 and is calculated by assigning points for severity of metabolic acidosis, need for emergency delivery, 10-min Apgar score <5, assisted ventilation at 10 min of life, and Sarnat scores as previously described (32, 33). This score is linked with greater regional injury on brain MRI (32). For our study, the PI score was used as a global assessment of severity of perinatal HI insult based on the NICHD criteria for TH initiation and to determine its relationship with serum biomarkers.
Biomarker Measurements
The set of serum biomarkers reported here was chosen from those linked to brain injury in adults and enough supporting evidence in the pediatric population as specified in the Request for Applications (RFA) from the National Institutes of Health, funding the study. Of 201 patients included in the study, 178 had serum clinical remnants collected at admission (prior to TH initiation) and daily from DOL 0 to 6. All serum samples were aliquoted and stored at −80°C. Each aliquot was exposed to only 0–1 thaw/freeze cycle prior to assaying. A custom multiplex enzyme-linked immunosorbent assay (ELISA) was used to measure BDNF, VEGF, IL-6, IL-8, and IL-10 in 5 μl of serum or GFAP in 12 μl of serum in 96-well plate formats [Meso Scale Discovery (MSD), Rockville, MD] as previously described (34). The interassay coefficient of variation for BDNF, VEGF, IL-6, IL-8, IL-10, and GFAP assays were 3.4, 5.9, 5.1, 6.7, 5.6, and 10.9%, respectively. Tau was measured in 5 μl of serum using a commercial ELISA (N451LAA-1, MSD) with an interassay coefficient of variation of 9.0%.
Statistical Analysis
Data were not normally distributed; thus, results are summarized as median and interquartile range (IQR) for continuous variables. Non-parametric statistical methods were applied to analyze continuous variables including Mann–Whitney U or Kruskal–Wallis (KW) with Dunn post-hoc tests and Spearman ρ correlations. χ2 was applied for categorical variables. Longitudinal biomarker trajectories were analyzed using linear trend of medians (35). Mixed-effect model adjusted for sex to assess the relationships between single indicators of HIE severity (independent variable) and serum biomarkers evaluated for repeated measurements and stratified by time relative to exposure to TH was applied. Additional adjusting for sex and neonatal infection was also modeled. p ≤ 0.05 was considered significant. Statistical analysis was performed using GraphPad Prism [version 8.0.0 (131) 2018; GraphPad Software, San Diego, CA] and Stata (version 15, StataCorp LLC, College Station, TX).
Results
Demographic, Clinical, and Biochemical Variables
Of 246 infants diagnosed with NE and treated with TH during the study period, 45 met the exclusion criteria as detailed in Figure 1, resulting in 201 patients diagnosed with HIE. Before exclusions, 16 deaths occurred in our cohort (6.5%; 16/246), with 2 deaths not linked to HIE and 4 meeting other exclusion criteria. As a result, the mortality rate for our cohort after exclusions was 5% (10/201, Table 1). Of the 201 infants with HIE treated with TH included in the analysis, 23 did not have serial serum samples for biomarker analysis (Figure 1). Patients without available biomarker data (n = 23) had a higher incidence of meconium-stained amniotic fluid (p = 0.02), lower Apgar score at 1 min of life (p = 0.03), and higher death rate (p < 0.01) than patients with biomarkers available (Table 1).
Longitudinal Changes in Serum Biomarkers
The median (IQR) number of serum samples per patient for biomarker analysis was 5 (3, 7). This resulted in 848 samples analyzed from 178 patients (Table 2). Serum BDNF levels decreased after DOL 0 (first 24 h of life) and recovered to levels at and above admission by DOL 5 (p = 0.02, Table 2). Vascular endothelial growth factor levels tended to be lower at admission (before initiation of TH) and recover after 24 h (p = 0.19, Table 2). Both tau and GFAP levels tended to peak by DOL 5–6 (p = 0.19 and p = 0.06, respectively; Table 2). Cytokine levels peaked within the first 24 h of life and decreased thereafter (p < 0.001 in all cases, Table 2).
Correlation Between Clinical Indicators of HIE Severity and Serum Biomarkers
At admission to the NICU, lower 5-min Apgar scores and BE and higher lactate and nucleated red blood cells (NRBCs), all indicators of worse HI insult, correlated with higher Sarnat scores (Table 3). Among the early non-brain-specific indicators of HI insult, higher NRBCs (a marker of prolonged partial in utero hypoxia) correlated with lower BDNF and VEGF and higher tau, and cytokines, but not GFAP, whereas higher lactate (a sub-acute HI marker) mostly correlated with lower VEGF and higher tau, GFAP, and cytokines, but not BDNF, during the first week of life (Tables 4A,B). All other early non-brain-specific indicators of HI insult severity (5-min Apgar score, worse pH, and worse BE) mostly correlated with tau (throughout the first week of life) and IL-10 (during the first 72 h of life), explaining the correlation between tau and IL-10 with the PI score, a global assessment of perinatal insult (Tables 4A,B). Higher Sarnat score (worse encephalopathy) correlated with lower VEGF during the first 24 h of life (DOL 0) and higher tau, GFAP, and IL-10 throughout the first week of life (Tables 4A,B). The integration of the proinflammatory (IL-6 and IL-8) and the anti-inflammatory (IL-10) responses in a patient-by-patient basis was suggested by the lost correlation between the (IL-6 * IL-8)/IL-10 index and most clinical and biochemical indicators of systemic HI severity (i.e., BE, NRBCs, and lactate). Conversely, the persistent correlation between Sarnat scores and this index suggested a disproportionate IL-10 response with worse encephalopathy (Table 4B).
Table 4A. Correlation between indicators of HI severity and BDNF, VEGF, tau, and GFAP during the first week of life.
Table 4B. Correlation between indicators of HI severity and cytokines during the first week of life.
Stratification of Biomarker Trajectories by Sarnat Score
Sarnat scores were inversely correlated with VEGF and directly with tau, GFAP, and IL-10 (Tables 4A,B). Patients with Sarnat scores of 2 and 3 had a VEGF nadir soon after admission, which recovered after 24 h of life [p < 0.001 (Sarnat 2), p = 0.008 (Sarnat 3); Figure 2A]. Unlike VEGF, GFAP discriminated between patients with Sarnat scores of 2 and 3. Only patients with Sarnat scores of 3 had higher GFAP levels since admission to the NICU (p = 0.03 vs. Sarnat 1; p = 0.02 vs. Sarnat 2; KW p = 0.04) and throughout the first 72 h of life (p = 0.003, Figure 2C). Tau, the other marker of neural cell injury, only began to increase after the first 24 h of life (DOL 1) in patients with Sarnat score of 3, reaching levels fourfold to sixfold higher than those in patients with lower scores by 48 to 72 h of life (DOL 2) (p < 0.001 vs. either Sarnat 1 or Sarnat 2; Figure 2B). Conversely, patients with score of 2 only began to increase their tau and GFAP levels around the time of TH completion (72 h of life, DOL 2), resulting in approximately three-fold [p = 0.01 (tau)] and approximately six-fold [p = 0.07 (GFAP)] higher levels than those in patients with scores of 1 [p = 0.006 (tau) and 0.004 (GFAP)] by DOL 5 (Figures 2B,C). At admission, patients with Sarnat scores of 3 also tended to have higher IL-10 levels than those with Sarnat scores of 1 (p = 0.05; Figure 2D). Although IL-10 levels decreased after the first 24 h of life (p < 0.001 for all Sarnat groups), they remained relatively elevated in patients with Sarnat scores of 2 and 3, and by DOL 2 (48–72 h of life) prior to completion of TH, IL-10 levels were more than twice the levels in patients with Sarnat of 1 (p = 0.009; Figure 2D).
Figure 2. Temporal trajectories of Sarnat score–stratified serum biomarkers. Temporal trajectories of VEGF (A), tau (B), GFAP (C), and IL-10 (D) levels during the first 7 days of life (DOL 0–6) stratified by Sarnat score (1, white; 2, gray, and 3, black) are shown. Two specimens were collected during first 24 h of life (DOL 0), the first sample obtained after admission and the one collected in the second 12 h of life. Data are represented as box-and-whiskers plot, where the box represents the IQR limited by the 25th and 75th percentiles, with solid line inside the box representing the median. Whiskers are limited by the last data point sitting within 1.5 times the IQR, whereas outliers are represented as circles beyond the boundaries of the whiskers. The biomarker level in a 10-base log scale is shown in the y axis. Time in hours and DOL for all panels are depicted in the bottom of the figure. TH exposure occurred for 72 h upon admission to the NICU (for simplification DOL 0–2). *p < 0.05 (KW analysis of variance with Dunn Bonferroni post-hoc analysis). DOL, day of life; GFAP, glial fibrillary acidic protein; IL, interleukin; TH, therapeutic hypothermia; VEGF, vascular endothelial growth factor.
Adjusted Relationship Between Severity of HI Insult and Serum Biomarkers
After adjustment for sex [Figures 3A(a−f),B; Supplementary Table 1] and infection (Table 5; Supplementary Table 2) and stratification by TH (DOL 0–2) and post-TH (DOL 4–6) periods, most relationships between indicators of HI insult severity and serum biomarkers persisted as described above. Tau and GFAP had the strongest association with 5-min Apgar score [Figure 3A(a)], worse BE [Figure 3A(c)], and lactate [Figure 3A(e)], but lacked association with NRBCs [Figure 3A(d)], during and after TH. Lower VEGF during TH and higher tau, GFAP, and IL-10 during and after TH were most strongly associated with worse Sarnat score [Figure 3B]. Some of these relationships changed during and after TH. For example, lower 5-min Apgar scores were associated with higher GFAP, but this association was 75% stronger during TH than after TH [Figure 3A(a)], despite higher GFAP levels after TH (Table 2). In addition, the association between lower BE and higher IL-10 levels was six-fold stronger during TH than after TH [Figure 3A(c)]. Conversely, the strength of association between lower 5-min Apgar and higher IL-8 [Figure 3A(a)] and higher NRBCs and higher IL-6 [Figure 3A(d)] became stronger after TH. Lastly, the associations between tau and GFAP with Sarnat scores became 50% stronger after TH [Figure 3B]. Adjustment for infection did not change most relationships between indicators of severity of HI insult and cytokines (Supplementary Table 2).
Figure 3. Forest plots for adjusted associations between HI severity indicators and serum biomarkers by exposure to TH. Forest plots represent relationships adjusted for sex between (A) indicators of severity of HI insult and (B) Sarnat scores with biomarkers levels during (DOL 0–2, red) and after (DOL 4–6, blue) TH. Admission and DOL 3 (TH-rewarming) data were not included in the adjusted analysis intended to assess the effect of TH in trajectories. β-Coefficients reaching p-values < 0.05 in mixed model analysis adjusted for sex do not cross line at zero. *p-value comparing during and after TH periods is shown. The strength of associations, as determined by the β coefficients and CI (95%), is assessed by the distance from zero (reference line). (C) Proposed timeline to maximize utility of serum biomarkers against Sarnat score. Time intervals of maximum utility are shown in horizontal boxes in gray for levels expected to decrease and in black for levels expected to increase with higher HI brain injury. Red and blue zones represent the period during and after TH, respectively. Hashed box for IL-10 represent speculated increased levels with higher HI injury.
Table 5. Summary of mixed-model adjustments for indicators of HI severity, time of biomarker, sex, and neonatal infection status.
Discussion
The strength of the relationships between serum biomarkers and clinical indicators of HIE severity is modulated by TH. Here, we show that after adjusting for sex and infection, lower BE and higher lactate (biochemical markers of decreased perfusion), and worse Sarnat scores (functional outcome) all relate to lower VEGF and higher tau, GFAP, and IL-10. Within the first 24 h of life, patients with Sarnat scores of 2–3 have lower VEGF, whereas those with scores of 3 also have high GFAP and IL-10. Thus, measuring VEGF, GFAP, and IL-10 soon after admission to the NICU may assist to stratify more accurately moderate–severe HI brain injury. During the next 72 h (TH period), tau levels begin to increase. However, while in patients with a Sarnat score of 3 the increase in tau begins on DOL 0–1, in those with a score of 2, this increase is delayed until DOL 2–3. By 72 h after completion of TH (DOL 5), patients with Sarnat scores of 2 reach tau, GFAP, and IL-10 levels similar to those of patients with scores of 3. The increase in tau and GFAP levels after TH may reflect the role of TH in delaying activation of injurious pathways, which are the target of multiple adjuvant therapies under study (36–38). In adjusted models, the associations between Sarnat scores with tau and GFAP become stronger after TH. Whether these associations represent persistent or worsening HI brain injury after completion of TH requires further investigation. In light of our data, we propose a screening schedule of HIE patients using VEGF, tau, GFAP, and IL-10 at specific time points related to TH (Figure 3C).
Hypoxic-ischemic encephalopathy severity is traditionally assessed by a combination of clinical and biochemical indicators, which, despite their lack of brain specificity, are used to guide TH initiation. Although Apgar and Sarnat scores are subjective (25, 39), the Sarnat score remains among the most brain-specific indicator available to assess clinically HI injury (40). Other indicators available at admission to the NICU, such as pH, BE, NRBCs, and lactate levels, are systemic markers of decreased fetal perfusion and are not brain-specific. Scores providing a global assessment of severity of perinatal insult, such the PI score (32), lack sensitivity and specificity to HI brain injury; whereas cerebral oximetry and electroencephalographic activity, which are more specific indicators of brain injury, are not broadly used due to the expertise required for their interpretation (41, 42). Thus, circulating biomarkers for HI brain injury may improve the ability of clinical and biochemical indicators alone to assess severity of injury. However, a better understanding of how TH affects their relationships is needed.
Among the molecules measured for this study, tau and GFAP are the most nervous system-specific proteins, which get spilled to the blood stream upon injury of neurons and astrocytes, respectively (34, 43). Because of their different cellular origins, it is not surprising that they have different temporal trajectories. Glial fibrillary acidic protein can discriminate Sarnat 3 from lower scores at admission to the NICU, whereas tau levels peak until the end of TH in neonates with Sarnat 3 and after TH in those with score of 2. Glial fibrillary acidic protein levels did not discriminate between Sarnat 1 and 2 at any time point. This may be due to the sensitivity of the study to detect differences, as only ~20% of the cohort had a Sarnat of 1, or may reflect similar levels of astrocyte injury in these two groups. When adjusted for sex and infection, worse Sarnat score was associated with higher tau and GFAP levels, more so after TH. Whether patients with more severe HI brain insults develop worsening injury with increasing tau and GFAP levels after completion of TH needs further investigation. Our results agree with other studies that have confirmed the direct association between Sarnat scores and circulating tau (19, 22, 23, 34, 44) and GFAP (21, 24, 34). However, our study is the first to address the potential role of TH/rewarming in the temporal trajectory of these biomarkers.
Vascular endothelial growth factor supports the blood–brain barrier (BBB) after neonatal HI brain injury (45), but its production may be compromised with severe insults. The association between high Sarnat score and decreased VEGF levels during TH shown here may suggest that worse HI brain insults may impair VEGF production. Associations between VEGF and early non-brain-specific indicators, such as low pH and BE and high lactate and NRBCs, persist before and after TH. One possible explanation is that injury of other organ systems may carry a significant influence in circulating VEGF levels, particularly after TH. As a result, the early decrease in VEGF may be primarily cause by HI brain injury with some systemic influence, whereas any deficit seen after TH may be mostly linked to persistent systemic compromise. The influence of the recovery of the BBB microstructure in these relationships after TH (46) remains unexplored. In preclinical models, neonatal HI injury leads to delayed regional decrease in BDNF levels (47, 48), the second neurotrophin studied here. In agreement, we report a decrease in serum BDNF levels in infants after HIE, which shows nonetheless no association with Sarnat scores. In a smaller cohort, Massaro et al. reported correlation between low BDNF levels and worse neurodevelopmental outcomes at age 1, but the analysis was unadjusted, limiting interpretation (19). In a larger cohort, higher Sarnat scores were linked to lower BDNF levels (34), an association likely powered by the inclusion of a group with a Sarnat score of 0. As infants with Sarnat scores of 0 are rarely assessed for TH, from a translation perspective, our results suggest that VEGF is a better candidate than BDNF for predictive modeling in the NICU. Our finding that high NRBCs relate to low BDNF may suggest a link with partial prolonged hypoxic events in utero.
Cytokines, such as IL-6 (15, 29, 49), IL-8 (21, 49), and IL-10 (49, 50), were the first peripheral blood and cerebrospinal fluid biomarkers described for HIE (51), but validation against multiple indicators of severity of HI in a longitudinal infection-adjusted model relative to TH has not been done in a cohort this large. There is a strong association between Sarnat scores with IL-10 during and after TH; however, such associations do not exist with IL-6 and IL-8. Circulating IL-10 may be the best initial biomarker of severity of HI insult to the brain among the chosen cytokines. Conversely, the infection-adjusted associations of non-specific indicators of severe HI insult with proinflammatory cytokines, IL-6 and IL-8, may suggest a systemic inflammatory response, which becomes stronger after TH is completed. Thus, previous associations between higher IL-6 and IL-8 with worse neurodevelopmental outcomes in HIE patients (52) may not be directly linked to severity of brain injury but instead to the degree of multiorgan involvement leading to systemic inflammation, which influences final brain outcomes (51). Altogether, worse global HI injury, not necessarily to the brain, may lead to an early inflammatory response with cytokine release, including IL-10, a response that may persist or worsen after TH because of injury to other organs. Higher IL-6, IL-8, and IL-10 on DOL 0 may help determine severity of global HI insult, whereas persistent elevation of IL-10 before completion of TH and thereafter, in combination with tau and GFAP, may guide the assessment of ongoing brain injury.
Our study has several limitations. Although recent reports have not identified associations between timing of TH initiation and dysfunctional cerebral autoregulation (32), injury on brain MRI or neurodevelopmental outcomes (53), delay in the initiation of TH may influence the trajectory of the proposed biomarkers. In addition, systemic hypotension and disturbed autoregulation, hypoxia, and hypercapnia may further increase brain injury after HI and influence the strength of associations between Sarnat scores with tau and GFAP. Clinical practice drift has occurred over the course of this study (2009–2019), but the influence of these changes on the longitudinal trajectory of serum biomarker presented here is difficult to isolate. We also were limited by blood volume available for analysis preventing quantification of other biomarkers previously studied in association with neonatal HIE. Although this is one of the largest biomarker studies in neonatal HIE and ~90% of the whole cohort had serum samples available for analysis, selection bias is still likely, as a subset of the sicker patients lacked longitudinal blood samples to measure biomarkers. While freezing and thawing of samples were limited, the samples used in the study were frozen for several years, which may have led to sample degradation, which could impact biomarker levels. Lastly, the usefulness of the proposed biomarkers in predicting the outcomes of patients treated with TH for HIE needs future studies.
Conclusions
Understanding the influence of TH on the temporal trajectory of serum biomarkers commonly tested to stratify severity of neonatal HIE is an essential step in developing precision in diagnosis of brain injury and building of predictive models. We show the temporal effect of TH in the relationships between serum biomarkers and traditional indicators of severity of HI insult to the brain and other organs. Admission and the period before completion of TH and rewarming may prove to be optimal time points to use biomarkers to stratify patients to modified and novel adjuvant therapeutic strategies. Considering our analysis, we proposed screening of HIE patients at three critical time points: (i) within the first 24 h of life using VEGF, GFAP, and IL-10; (ii) before completion of TH using tau, GFAP, and IL-10; and (iii) by 72 h after completion of TH with tau and GFAP (Figure 3C). These time-sensitive serum biomarkers would improve the stratification of severity of HI brain injury to assign patients to novel neuroprotective or restorative therapies for HIE and to follow therapeutic response.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Ethics Statement
The studies involving human participants were reviewed and approved by Johns Hopkins School of Medicine IRB. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
Author Contributions
RC-V, DV, CP, GG, AT, SB, EG, FN, and AE: conceptualization and design of the study. RC-V, SM, HS, DV, FN, and AE: methodology. RC-V, CP, BD, EG, FN, and AE: supervision and oversight. RC-V, FN, and AE: funding acquisition. RC-V and DV: formal data analysis. RC-V, DV, CP, BD, EG, GG, AT, SB, FN, and AE: resources. RC-V, SM, HS, and AE: drafting of significant portions of manuscript, tables, and figures. All authors: data acquisition, review, and editing of final manuscript.
Funding
This work was supported by National Institutes of Health RO1HD086058 (AE, FN, EG, DV, CP, AT, GG, and SB); RO1HD070996, AG061643, NS109029, and HD074593-07 (FN); KO8NS096115 and 3K08NS096115-03S1 (RC-V); the JHU-SOM Clinician Scientist Award (RC-V); the Thomas Wilson Foundation (RC-V). Under a license agreement between ImmunArray Ltd. and the Johns Hopkins University, the University and AE are entitled to royalties on an invention described in this study and discussed in this publication.
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.
Publisher's Note
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.
Acknowledgments
We thank the families of the extremely sick neonates included in this study for their willingness to participate. We also thank the nursing and ancillary staff at Johns Hopkins Hospital NICU and core laboratory for their support in the collection of the samples. Lastly, we thank Xueting Tao, MHS from the Department of Pediatrics at JHU for her assistance in formal data analysis, and Ms. Jie Zhu, research specialist at ADE lab for the processing of samples.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2021.748150/full#supplementary-material
References
1. Douglas-Escobar M, Weiss MD. Hypoxic-ischemic encephalopathy: a review for the clinician. JAMA Pediatr. (2015) 169:397. doi: 10.1001/jamapediatrics.2014.3269
3. Northington FJ, Zelaya ME, O'Riordan DP, Blomgren K, Flock DL, Hagberg H, et al. Failure to complete apoptosis following neonatal hypoxia-ischemia manifests as “continuum” phenotype of cell death and occurs with multiple manifestations of mitochondrial dysfunction in rodent forebrain. Neuroscience. (2007) 149:822–33. doi: 10.1016/j.neuroscience.2007.06.060
4. Northington FJ, Chavez-Valdez R, Martin LJ. Neuronal cell death in neonatal hypoxia-ischemia. Ann Neurol. (2011) 69:743–58. doi: 10.1002/ana.22419
5. Chavez-Valdez R, Martin LJ, Northington FJ. Programmed necrosis: a prominent mechanism of cell death following neonatal brain injury. Neurol Res Int. (2012) 2012:257563. doi: 10.1155/2012/257563
6. Shankaran S, Laptook AR, Ehrenkranz RA, Tyson JE, McDonald SA, Donovan EF, et al. Whole-body hypothermia for neonates with hypoxic–ischemic encephalopathy. N Engl J Med. (2005) 353:1574–84. doi: 10.1056/NEJMcps050929
7. Shankaran S, Pappas A, McDonald SA, Vohr BR, Hintz SR, Yolton K, et al. Childhood outcomes after hypothermia for neonatal encephalopathy. N Engl J Med. (2012) 366:2085–92. doi: 10.1056/NEJMoa1112066
8. Azzopardi DV, Strohm B, Edwards AD, Dyet L, Halliday HL, Juszczak E, et al. Moderate hypothermia to treat perinatal asphyxial encephalopathy. N Engl J Med. (2009) 361:1349–58. doi: 10.1056/NEJMoa0900854
9. Tagin MA, Woolcott CG, Vincer MJ, Whyte RK, Stinson DA. Hypothermia for neonatal hypoxic ischemic encephalopathy: an updated systematic review and meta-analysis. Arch Pediatr Adolesc Med. (2012) 166:558–66. doi: 10.1001/archpediatrics.2011.1772
10. Kurinczuk JJ, White-Koning M, Badawi N. Epidemiology of neonatal encephalopathy and hypoxic–ischaemic encephalopathy. Early Hum Dev. (2010) 86:329–38. doi: 10.1016/j.earlhumdev.2010.05.010
11. Hayes BC, McGarvey C, Mulvany S, Kennedy J, Geary MP, Matthews TG, et al. A case-control study of hypoxic-ischemic encephalopathy in newborn infants at >36 weeks gestation. Am J Obstet Gynecol. (2013) 209:29.e1–19. doi: 10.1016/j.ajog.2013.03.023
12. Rossi AC, Prefumo F. Antepartum and intrapartum risk factors for neonatal hypoxic-ischemic encephalopathy: a systematic review with meta-analysis. Curr Opin Obstet Gynecol. (2019) 31:410–7. doi: 10.1097/GCO.0000000000000581
13. Scheidegger S, Held U, Grass B, Latal B, Hagmann C, Brotschi B, et al. Association of perinatal risk factors with neurological outcome in neonates with hypoxic ischemic encephalopathy. J Matern Fetal Neonatal Med. (2019) 34:1020–7. doi: 10.1080/14767058.2019.1623196
14. Dalili H, Nili F, Sheikh M, Hardani AK, Shariat M, Nayeri F. Comparison of the Four Proposed apgar scoring systems in the assessment of birth asphyxia and adverse early neurologic outcomes. PLoS ONE. (2015) 10:e0122116. doi: 10.1371/journal.pone.0122116
15. Matsuo R, Ago T, Kamouchi M, Kuroda J, Kuwashiro T, Hata J, et al. Clinical significance of plasma VEGF value in ischemic stroke - research for biomarkers in ischemic stroke (REBIOS) study. BMC Neurol. (2013) 13:32. doi: 10.1186/1471-2377-13-32
16. Rubenstein R, Chang B, Yue JK, Chiu A, Winkler EA, Puccio AM, et al. Comparing plasma phospho tau, total tau, and phospho tau–total tau ratio as acute and chronic traumatic brain injury biomarkers. JAMA Neurol. (2017) 74:1063. doi: 10.1001/jamaneurol.2017.0655
17. Shang J, Yamashita T, Fukui Y, Song D, Li X, Zhai Y, et al. Different associations of plasma biomarkers in Alzheimer's disease, mild cognitive impairment, vascular dementia, and ischemic stroke. J Clin Neurol. (2018) 14:29. doi: 10.3988/jcn.2018.14.1.29
18. Liu F, Yang S, Du Z, Guo Z. Dynamic changes of cerebral-specific proteins in full-term newborns with hypoxic–ischemic encephalopathy. Cell Biochem Biophys. (2013) 66:389–96. doi: 10.1007/s12013-012-9478-3
19. Massaro AN, Wu YW, Bammler TK, Comstock B, Mathur A, McKinstry RC, et al. Plasma biomarkers of brain injury in neonatal hypoxic-ischemic encephalopathy. J Pediatr. (2018) 194, 67.e1–75.e1. doi: 10.1016/j.jpeds.2017.10.060
20. Aly H, Hassanein S, Nada A, Mohamed MH, Atef SH, Atiea W. Vascular endothelial growth factor in neonates with perinatal asphyxia. Brain Dev. (2009) 31:600–4. doi: 10.1016/j.braindev.2008.09.004
21. Chalak LF, Sánchez PJ, Adams-Huet B, Laptook AR, Heyne RJ, Rosenfeld CR. Biomarkers for severity of neonatal hypoxic-ischemic encephalopathy and outcomes in newborns receiving hypothermia therapy. J. Pediatr. (2014) 164:468.e1–74.e1. doi: 10.1016/j.jpeds.2013.10.067
22. Takahashi K, Hasegawa S, Maeba S, Fukunaga S, Motoyama M, Hamano H, et al. Serum tau protein level serves as a predictive factor for neurological prognosis in neonatal asphyxia. Brain Dev. (2014) 36:670–5. doi: 10.1016/j.braindev.2013.10.007
23. Lv H-Y, Wu S-J, Gu X-L, Wang Q-L, Ren P-S, Ma Y, et al. Predictive value of neurodevelopmental outcome and serum tau protein level in neonates with hypoxic ischemic encephalopathy. Clin Lab. (2017) 63:1153–62. doi: 10.7754/Clin.Lab.2017.170103
24. Ennen CS, Huisman TAGM, Savage WJ, Northington FJ, Jennings JM, Everett AD, et al. Glial fibrillary acidic protein as a biomarker for neonatal hypoxic-ischemic encephalopathy treated with whole-body cooling. Am J Obstet Gynecol. (2011) 205:251.e1–7. doi: 10.1016/j.ajog.2011.06.025
25. Sarnat HB, Sarnat MS. Neonatal encephalopathy following fetal distress: a clinical and electroencephalographic study. Arch Neurol. (1976) 33:696–705. doi: 10.1001/archneur.1976.00500100030012
26. Battin MR, Dezoete JA, Gunn TR, Gluckman PD, Gunn AJ. Neurodevelopmental outcome of infants treated with head cooling and mild hypothermia after perinatal asphyxia. Pediatrics. (2001) 107:480–4. doi: 10.1542/peds.107.3.480
27. Olsen SL, DeJonge M, Kline A, Liptsen E, Song D, Anderson B, et al. Optimizing therapeutic hypothermia for neonatal encephalopathy. Pediatrics. (2013) 131:e591–603. doi: 10.1542/peds.2012-0891
28. Gagne-Loranger M, Sheppard M, Ali N, Saint-Martin C, Wintermark P. Newborns referred for therapeutic hypothermia: association between initial degree of encephalopathy and severity of brain injury (what about the newborns with mild encephalopathy on admission?). Am J Perinatol. (2016) 33:195–202. doi: 10.1055/s-0035-1563712
29. Chalak LF, Nguyen K-A, Prempunpong C, Heyne R, Thayyil S, Shankaran S, et al. Prospective research in infants with mild encephalopathy identified in the first six hours of life: neurodevelopmental outcomes at 18–22 months. Pediatr Res. (2018) 84:861–8. doi: 10.1038/s41390-018-0174-x
30. Takenouchi T, Cuaycong M, Ross G, Engel M, Perlman JM. Chain of brain preservation—a concept to facilitate early identification and initiation of hypothermia to infants at high risk for brain injury. Resuscitation. (2010) 81:1637–41. doi: 10.1016/j.resuscitation.2010.08.001
31. DuPont TL, Chalak LF, Morriss MC, Burchfield PJ, Christie L, Sánchez PJ. Short-term outcomes of newborns with perinatal acidemia who are not eligible for systemic hypothermia therapy. J Pediatr. (2013) 162:35–41. doi: 10.1016/j.jpeds.2012.06.042
32. Gilmore MM, Tekes A, Perin J, Parkinson C, Spahic H, Chavez-Valdez R, et al. Later cooling within 6 h and temperatures outside 33–34°C are not associated with dysfunctional autoregulation during hypothermia for neonatal encephalopathy. Pediatr Res. (2020) 1–8. doi: 10.1038/s41390-020-0876-8
33. Lee JK, Poretti A, Perin J, Huisman TAGM, Parkinson C, Chavez-Valdez R, et al. Optimizing cerebral autoregulation may decrease neonatal regional hypoxic-ischemic brain injury. Dev Neurosci. (2017) 39:248–56. doi: 10.1159/000452833
34. Dietrick B, Molloy E, Massaro AN, Strickland T, Zhu J, Slevin M, et al. Plasma and CSF candidate biomarkers of neonatal encephalopathy severity and neurodevelopmental outcomes. J. Pediatr. (2020) 2020:S0022-3476(20)30832-5. doi: 10.1016/j.jpeds.2020.06.078
35. Machado JAF, Santos Silva JMC. Quantiles via moments. J Econ. (2019) 213:145–73. doi: 10.1016/j.jeconom.2019.04.009
36. Burnsed JC, Chavez-Valdez R, Hossain MS, Kesavan K, Martin LJ, Zhang J, et al. Hypoxia-ischemia and therapeutic hypothermia in the neonatal mouse brain – a longitudinal study. PLoS ONE. (2015) 10:e0118889. doi: 10.1371/journal.pone.0118889
37. Davidson JO, Wassink G, van den Heuij LG, Bennet L, Gunn AJ. Therapeutic hypothermia for neonatal hypoxic–ischemic encephalopathy – where to from here? Front Neurol. (2015) 6:198. doi: 10.3389/fneur.2015.00198
38. Chavez-Valdez R, Emerson P, Goffigan-Holmes J, Kirkwood A, Martin LJ, Northington FJ. Delayed injury of hippocampal interneurons after neonatal hypoxia-ischemia and therapeutic hypothermia in a murine model. Hippocampus. (2018) 28:617–30. doi: 10.1002/hipo.22965
39. Bharti B, Bharti S. A review of the Apgar score indicated that contextualization was required within the contemporary perinatal and neonatal care framework in different settings. J Clin Epidemiol. (2005) 58:121–9. doi: 10.1016/j.jclinepi.2004.04.012
40. Douglas-Escobar M, Yang C, Bennett J, Shuster J, Theriaque D, Leibovici A, et al. A pilot study of novel biomarkers in neonates with hypoxic-ischemic encephalopathy. Pediatr Res. (2010) 68:531–6. doi: 10.1203/PDR.0b013e3181f85a03
41. Dereymaeker A, Matic V, Vervisch J, Cherian PJ, Ansari AH, De Wel O, et al. Automated EEG background analysis to identify neonates with hypoxic-ischemic encephalopathy treated with hypothermia at risk for adverse outcome: a pilot study. Pediatr Neonatol. (2019) 60:50–8. doi: 10.1016/j.pedneo.2018.03.010
42. Mitra S, Bale G, Meek J, Tachtsidis I, Robertson NJ. Cerebral near infrared spectroscopy monitoring in term infants with hypoxic ischemic encephalopathy—a systematic review. Front Neurol. (2020) 11:393. doi: 10.3389/fneur.2020.00393
43. Graham EM, Burd I, Everett AD, Northington FJ. Blood biomarkers for evaluation of perinatal encephalopathy. Front Pharmacol. (2016) 7:196. doi: 10.3389/fphar.2016.00196
44. Toorell H, Zetterberg H, Blennow K, Sävman K, Hagberg H. Increase of neuronal injury markers Tau and neurofilament light proteins in umbilical blood after intrapartum asphyxia. J Matern-Fetal Neonatal Med. (2018) 31:2468–72. doi: 10.1080/14767058.2017.1344964
45. Greenberg DA, Jin K. From angiogenesis to neuropathology. Nature. (2005) 438:954–9. doi: 10.1038/nature04481
46. Saunders NR, Liddelow SA, Dziegielewska KM. Barrier mechanisms in the developing brain. Front Pharmacol. (2012) 3:46. doi: 10.3389/fphar.2012.00046
47. Han BH, D'Costa A, Back SA, Parsadanian M, Patel S, Shah AR, et al. BDNF blocks caspase-3 activation in neonatal hypoxia-ischemia. Neurobiol Dis. (2000) 7:38–53. doi: 10.1006/nbdi.1999.0275
48. Chavez-Valdez R, Martin LJ, Razdan S, Gauda EB, Northington FJ. Sexual dimorphism in BDNF signaling after neonatal hypoxia-ischemia and treatment with necrostatin-1. Neuroscience. (2014) 260:106–19. doi: 10.1016/j.neuroscience.2013.12.023
49. Orrock JE, Panchapakesan K, Vezina G, Chang T, Harris K, Wang Y, et al. Association of brain injury and neonatal cytokine response during therapeutic hypothermia in newborns with hypoxic-ischemic encephalopathy. Pediatr Res. (2016) 79:742–7. doi: 10.1038/pr.2015.280
50. Wu H, Liu G, Yang X, Liu Q, Li Z. Effect of mild hypothermia on the expression of IL-10 and IL-18 in neonates with hypoxic ischemic encephalopathy. Exp Ther Med. (2019) 18:2194–8. doi: 10.3892/etm.2019.7768
51. Foster-Barber A, Ferriero DM. Neonatal encephalopathy in the term infant: neuroimaging and inflammatory cytokines. Ment Retard Dev Disabil Res Rev. (2002) 8:20–4. doi: 10.1002/mrdd.10009
52. Shang Y, Mu L, Guo X, Li Y, Wang L, Yang W, et al. Clinical significance of interleukin-6, tumor necrosis factor-α and high-sensitivity C-reactive protein in neonates with hypoxic-ischemic encephalopathy. Exp Ther Med. (2014) 8:1259–62. doi: 10.3892/etm.2014.1869
53. Guillot M, Philippe M, Miller E, Davila J, Barrowman NJ, Harrison M-A, et al. Influence of timing of initiation of therapeutic hypothermia on brain MRI and neurodevelopment at 18 months in infants with HIE: a retrospective cohort study. BMJ Paediatr Open. (2019) 3:e000442. doi: 10.1136/bmjpo-2019-000442
Keywords: neonatal encephalopathy, cytokines, neurotrophins, GFAP, tau, Sarnat score
Citation: Chavez-Valdez R, Miller S, Spahic H, Vaidya D, Parkinson C, Dietrick B, Brooks S, Gerner GJ, Tekes A, Graham EM, Northington FJ and Everett AD (2021) Therapeutic Hypothermia Modulates the Relationships Between Indicators of Severity of Neonatal Hypoxic Ischemic Encephalopathy and Serum Biomarkers. Front. Neurol. 12:748150. doi: 10.3389/fneur.2021.748150
Received: 27 July 2021; Accepted: 04 October 2021;
Published: 02 November 2021.
Edited by:
Jean-Charles Sanchez, Université de Genève, SwitzerlandReviewed by:
Manish Bhomia, Uniformed Services University of the Health Sciences, United StatesJennifer Munoz Pareja, University of Miami Hospital, United States
Copyright © 2021 Chavez-Valdez, Miller, Spahic, Vaidya, Parkinson, Dietrick, Brooks, Gerner, Tekes, Graham, Northington and Everett. 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: Raul Chavez-Valdez, rchavez2@jhmi.edu