Original Research ARTICLE
Early Changes in Glutamate Metabolism and Perfusion in Basal Ganglia following Hypoxia-Ischemia in Neonatal Piglets: A Multi-Sequence 3.0T MR Study
- 1Department of Radiology, Shengjing Hospital of China Medical University, Shenyang, China
- 2Department of Imaging Systems Clinical Science, Philips Healthcare, Beijing, China
The excitotoxicity of glutamate metabolism as well as hemodynamic disorders of the brain are both risk factors for neonatal hypoxic–ischemic brain damage (HIBD). In the present study, changes in glutamate metabolism in the basal ganglia were detected by proton magnetic resonance spectroscopy (1H-MRS) at 0–6, 8–12, 24–30, and 48–60 h after the induction of hypoxia-ischemia (HI) in newborn piglets. Meanwhile, correlation analysis was performed by combining the microcirculatory perfusion informations acquired by intravoxel incoherent motion (IVIM) scan to explore their possible interaction mechanism. The results suggested that Glu level in the basal ganglia underwent a “two-phase” change after HI; perfusion fraction f, an IVIM-derived perfusion parameter, was clearly decreased in the early stage after HI, then demonstrated a transient and slight recovery process, and thereafter continued to decrease. The changes in f and Glu level were in a significant negative correlation (r = −0.643, P = 0.001). Our study results revealed that Glu level is closely associated with the microcirculatory perfusion changes in the acute stage of HIBD.
Basal ganglia injury (BGI) is a common type of hypoxic–ischemic brain damage (HIBD) and also the main cause of permanent dysneuria and cerebral palsy in perinatal full-term neonates. The basal ganglia are very sensitive to hypoxia–ischemia (HI) damage and susceptible to selective neuronal injury (Martin et al., 1997; Rocha-Ferreira and Hristova, 2016). Glutamate (Glu), as the most important excitatory neurotransmitter in animals, plays a crucial role in maintaining the function of glutamatergic neurons in the basal ganglia (Alexander and Crutcher, 1990). The excitotoxicity caused by the massive accumulation of Glu after HI is the central link for HIBD and is also the initiator and executer of brain injury (Hagberg et al., 1987; Choi and Rothman, 1990; Coyle and Puttfarcken, 1993; Rego et al., 1996).
Normally the release and reuptake of Glu are in a dynamic balance. The efficient Glu uptake system enables low, extracellular concentrations of Glu in order to avoid excitotoxicity, and only a small amount of Glu is involved in signal transduction as an excitatory neurotransmitter (D'souza and Slater, 1995; Cooper and Jeitner, 2016; Danbolt et al., 2016). However, HI leads to an increase in Glu release and/or damages the Glu uptake system, thus causing a sharp elevation of extracellular Glu; its level and the excessive activation of relevant receptors contributes to nervous excitotoxicity and results in generalized pathological brain lesions. Of these excessively activated receptors, the N-methyl-D-aspartic acid receptor plays an important role in neuronal injury after mediating HI as outlined in our preliminary study (Wang X. Y. et al., 2012).
Perinatal HIBD has a complex pathophysiological mechanism. In addition to the excitotoxicity of Glu, hemodynamic disorders of the brain are also a risk factor for HIBD (Pryds et al., 1990; Howlett et al., 2013; Massaro et al., 2015). Hemodynamic disorders of the brain can contribute to secondary energy dysmetabolism after HI, and thus have a role in pathological mechanisms involved in brain damage secondary to HI. The pathogenesis of brain damage can be revealed by evaluating cerebral perfusion. For example, a pivotal question arises concerning cerebral metabolism in a HI environment in response to changes in cerebral perfusion: do mutual effects and interactions exist between such variables? In order to understand HIBD mechanisms, it is essential to explore the relationship between changes in Glu metabolism and changes in microcirculatory perfusion after HI.
Proton magnetic resonance spectroscopy (1H-MRS), as a non-invasive technique, provides information on changes of metabolites in brain tissues. At present, two major functional magnetic resonance imaging (MRI) techniques are available to measure cerebral perfusion: perfusion-weighted imaging (PWI) and arterial spin labeling (ASL). Considering the specificity of the newborns, invasive PWI cannot be used as it requires the use of exogenous contrast media. With regard to ASL, this can only evaluate a single parameter—cerebral blood flow (CBF)—which does not allow the all-round comprehension of a hemodynamic reserve. This leads us to the question of whether an MRI technique exists that can not only meet the special requirements of non-invasive examinations in newborns, but also support the comprehensive evaluation of cerebral perfusion? Intravoxel incoherent motion (IVIM), a non-invasive MR perfusion imaging technique, can quantify the microcirculation of blood in the capillary network (perfusion) and the diffusion composition of true water molecules in tissues in vivo using a bi-exponential model (Le Bihan et al., 1986, 1988). Measurements using IVIM are based on the premise that the microcirculation of blood, or perfusion, is a kind of non-uniform, irregular random motion (i.e., incoherent motion) (Le Bihan et al., 1988). A linear correlation exists between IVIM-derived perfusion parameters (pseudo-diffusion coefficient D* and perfusion fraction f) and conventional perfusion parameters (Le Bihan and Turner, 1992; Federau et al., 2014a,b). In theory, therefore, IVIM-derived perfusion parameters can be used to evaluate the microcirculatory perfusion of the brain after HI and thus provide more comprehensive information. Currently, the successful application of IVIM in the central nervous system focuses on studies to recognize the ischemic semi-dark band (Federau et al., 2014b; Hu et al., 2015), identify the benignity and malignancy of tumors(Bisdas et al., 2014; Hu et al., 2014; Suh et al., 2014) and evaluate the efficacy of radiochemotherapy (Hauser et al., 2014; Cui et al., 2015; Xiao et al., 2015). In this context, IVIM, perhaps, will become a new MR perfusion technique useful in the study of perfusion of the microcirculation in newborns.
Studies are urgently required on the critical changes that occur at the Glu level during the development of HIBD and its relationship with microcirculatory perfusion. In this study, we measured Glu-related metabolites and microcirculatory perfusion separately by 1H-MRS and IVIM in a HI animal model. We investigated the dynamic changes in Glu and microcirculatory perfusion after HI and preliminarily explored their possible interaction mechanism.
Materials and Methods
Preparation of Experimental Animals and Establishment of HIBD Model
All animal experiments were reported in accordance with Animal Research: Reporting of In vivo Experiments guidelines. This study was approved by the Institutional Animal Care and Use Committee of Shengjing Hospital of China Medical University. Twenty-five newborn male or female Yorkshire piglets (P3–5 d; weight: 1.5–2 kg) were randomly selected from the Laboratory Animal Center of Shengjing Hospital of China Medical University. All animal models were established according to the Regulations for the Animal Care and Use published by the Shengjing Hospital of China Medical University.
The newborn piglets were anesthetized by an intramuscular injection of 0.6 mL/kg Su-Mian-Xin (xylazine hydrochloride) in the buttocks. During anesthesia, the animals' vital signs were closely observed. When the piglets were found to become comatose, as evidenced by muscle relaxation, decreased limb muscular tension and a delayed corneal reflex, the animals were placed in a supine position on a bench prior to each operation. A laryngoscopic tracheal cannulation (φ 2.5 mm) was performed, and each animal was then connected to a TKR-200C small animal ventilator (Jiangxi Teli Anesthesia and Respiratory Equipment Co., Ltd, Jiangxi, China) for mechanical ventilation with 100% oxygen. Ventilator parameters were as follows: respiration ratio inspiration/expiration (I/E) = 1:1.5; and respiration frequency = 30 bpm. Heart rate and peripheral oxygen saturation were monitored using a TuffSat handheld pulse oximeter (GE, Boston, MA, USA). The incision area and adjacent skin were then disinfected and a median incision was made in the neck. Bilateral common carotid arteries and adjacent internal jugular veins and vagus nerves were dissected, and a 5.0 mm silk suture was indwelled. After the condition of animals stabilized for 30 min, the bilateral common carotid arteries were clipped with small artery clamps to interrupt their blood flow, and 6% oxygen–containing mixed gas was delivered mechanically for 40 min. Thereafter, the HI induction procedure was completed, the small artery clamps removed, and the blood flow of the bilateral common carotid arteries was recovered. Oxygen (100%) was mechanically delivered again, and the incision was sutured. After the operation, the animals were transferred to an incubator (37°C) to ensure a body temperature in the normal range during postoperative recovery. This well-established model was used to induce bilateral HI injury, as documented in our previous studies (Wang H. et al., 2012; Zhang et al., 2012).
Grouping of Experimental Animals
All newborn piglets were randomly divided into the control group (sham-operation group, n = 5) and the HI model group (n = 20), and the HI model group was then further divided into 4 subgroups according to different time points after HI: 0–6, 8–12, 24–30, and 48–60 h (n = 5 per group). The animals in the pseudo-operation group underwent preoperative preparation similar to those in the model group but without a HI induction procedure.
MR Scans and Data Post-Processing
At different time points after HI, a conventional MR scan was performed on all animals using a 3.0T MRI system (Achieva 3.0T TX; Philips Healthcare Systems, Best, The Netherlands) with an eight-channel phased array head coil. Images (including axial and sagittal T1-weighted images [T1WI] and axial T2-weighted images [T2WI] of the head) were acquired by fast gradient echo. The relevant scan parameters used were: T1WI Repetition time (TR) 200 ms; Echo time (TE) 2.3 ms; matrix 224 × 162; slice thickness 5 mm; T2WI TR 5000 ms; TE 80 ms; matrix 224 × 162; and slice thickness 5 mm.
Measurements by single-voxel 1H-MRS were completed with a point-resolved spectroscopy (PRESS) sequence using one 90° and two 180° radio frequency pulses (13.2224 ms; bandwidth = 1,231 Hz). The following parameters were used: TR = 2,000 ms; TE = 37 ms; bandwidth = 2,000 Hz; number of signal acquisitions = 64; and a volume of interest (VOI; 10 × 10 × 10 mm) was positioned in the left basal ganglia. A plain scan was performed on newborn piglets to acquire coronal scan images of the basal ganglia as a region of interest (ROI). After the completion of a 1H-MRS scan, the raw data were then entered into LCModel (Linear Combination Model) software (version 6.3-1B) (Provencher, 2001) for quantification, with an unsuppressed water signal used as an internal reference, as well as automatic processing and analysis. The LCModel software supported automatic baseline corrections and the smoothing of spectrum raw data, as well as the absolute concentration measurement of several metabolites, including alanine (Ala), aspartic acid (Asp), creatine (Cr), phosphocreatine (PCr), γ-aminobutyric acid (GABA), glucose (Glc), glutamate (Glu), glutamine (Gln), glycerylphosphorylcholine (GPC), phosphatidylcholine (PCh), glutathione (GSH), inositol (Ins), lactic acid (Lac), N-acetylaspartate (NAA), N-acetylaspartate glutamate (NAAG), scyllo-inositol (Src) and taurine (Tau). Glu is the most abundant excitatory neurotransmitter in the brain and is essential for normal brain function; its complex signals were generated at 2.04–2.35 ppm (ppm: 10−6) and 3.75 ppm. The peak of Gln overlapped with that of Glu at approximately 2.35 ppm. LCModel software can automatically distinguish the overlapping signals of different metabolites with a similar chemical shift in the same frequency area (Wisnowski et al., 2013). The spectral fitting was performed with a range from 0.2 to 4.0 ppm (see Figure 1). In this study, we quantified the concentrations of Glu, Gln and Glx (Glu + Gln complex). Cramér–Rao lower bounds (CRLBs), which were provided by LCModel software, were used in assessing the precision of the quantification of metabolites. Spectrum data were included in the statistical analysis, which met the following criteria: (1) signal–to–noise ratio (SNR) ≥ 5; and (2) CRLBs for the concentrations of metabolites obtained were < 50%, and generally < 25%.
Figure 1. Representative T2-weighted image and corresponding 1H-MRS of a normal piglet. A MRS voxel was placed in the left basal ganglia (black square). LCModel fitting results (red), fit from 0.2 to 4.0 ppm, are presented. Glu, glutamate; Gln, glutamine; NAA, N-acetylaspartate; Lac, lactate; Cho, choline; Cr, creatine.
An IVIM scan was performed with a single-shot spin echo planar imaging (EPI) sequence; TR/TE 4,000 ms/90 ms; matrix 228 × 231; slice thickness 5 mm; and 16 b-values (0, 10, 20, 40, 80, 110, 140, 170, 200, 300, 400, 500, 600, 700, 800, and 900 s/mm2). IVIM analysis was performed using in-house MATLAB software (R2010a, The MathWorks, Inc., Natick, MA, USA) and a bi-exponential model. The relation between the signal attenuation and b is expressed by the following formula (Le Bihan et al., 1988, 1989):
Where Sb and S0 stand for signal intensity at b ≠ 0 and b = 0, respectively. The relevant parameters finally obtained by least squares (Pfeuffer et al., 1999) are as follows: diffusion coefficient D (given in units of × 10−3 mm2/s) denotes the true molecular diffusion; the pseudo-diffusion coefficient D* (given in units of × 10−3 mm2/s) represents the diffusion linked to microcapillary perfusion; and the perfusion fraction f (given as a percentage) indicates the proportion of microcirculatory perfusion-related diffusion in the total diffusion. As there is a variation in the speed of motion of different molecules, D* is markedly greater than D. At b < 200 s/mm2, the perfusion effect is dominant and the signal attenuation detected reflects information from perfusion; at b > 200 s/mm2, the microcirculatory perfusion effect can be neglected, and the signal attenuation detected nearly reflects the true diffusion of water molecules in voxels.
Firstly, IVIM raw data were imported into post-processing software, and then post-processed using a bi-exponential model to give pseudo-color images of D, D*, and f. Secondly, the ROI was manually marked on the left basal ganglia in images at b = 0 s/mm2 (see Figure 2) by combining T2WI scan images (avoiding the adjacent cerebrospinal fluid [CSF], blood vessels and noise areas), and then D, D*, and f values were calculated. If the calculated results were < 0, f, D, and D* were set to 0 to comply with actual physiology. If f was > 0.3 and D* was > 0.05 as calculated, f and D* were also set to 0 because such calculated results may be caused by the effects of SNR or CSF flow, and were therefore not physiological values (Federau et al., 2012). Finally, the relevant parameters of ROI in model and control groups were measured repeatedly three times and results then averaged. Differences in various parameters between model and control groups were analyzed.
Figure 2. Representative IVIM images of newborn piglets before and after HI. (a–d) showed the axial diffusion-weighted imaging (DWI) image (b = 0 s/mm2) of a normal newborn piglet (a) (where the blue ellipse highlights the region of interest (ROI) marked in the left basal ganglia region), and the corresponding D, D*, and f images obtained by IVIM post-processing software (b–d; D = 0.642 × 10−3 mm2/s, D* = 13.207 × 10−3 mm2/s, f = 12.753%). (e–h) showed the axial DWI image (b = 0 s/mm2) of a newborn piglet at 6 h after HI (e), and the corresponding D, D*, and f images (f–h; D = 0.269 × 10−3 mm2/s, D* = 5.412 × 10−3 mm2/s, f = 7.197%). From the axial DWI images, we could see that at 6 h after HI, the signals of cerebral parenchyma were evidently enhanced, the cerebral cortex was obviously swollen and edematous, and the cortical sulci and gyrus became shallower; in the D and f images, weakened signals of cerebral parenchyma were observed after HI correspondingly, but the difference in the image of D* was not significant.
SPSS 20.0 statistical software was used for analyses. All data were expressed as mean ± standard deviation (SD). Multiple comparisons of data with a homogeneity of variance were performed by one-way analysis of variance (ANOVA), and those of data with a heterogeneity of variance by a Kruskal–Wallis H test. Correlations between the Glu level and IVIM-derived perfusion parameters (D* and f) were analyzed by Spearman rank correlation analysis. The correlation level indicated by the correlation coefficient was defined as follows: r > 0.8, very high correlation; 0.6 < r ≤ 0.8, significant correlation; 0.4 < r ≤ 0.6, ordinary correlation; 0.2 ≤ r ≤ 0.4, low correlation; and 0 ≤ r < 0.2, weak or no correlation (Fujima et al., 2014). P < 0.05 suggested that a difference was statistically significant.
In this study, 1H-MRS and IVIM scans were performed on the newborn piglets in the control group and HI model group at different time points of 0–6, 8–12, 24–30, and 48–60 h after HI. The results are described below.
The representative 1H-MRS in the basal ganglia at different time points after HI are shown in Figure 3. Figure 4 shows the mean levels of Glu, Gln, and Glx in the control group and model group. After HI, there was a sharp increase in Glu level in the basal ganglia at 0–6 h, followed by a transient decrease at 8–12 h, to a level that was still higher than that in the control group; thereafter, it increased again, demonstrating a “two-phase” change. Compared with the control group, the Glu level was increased at different time points after HI, and differences were statistically significant (0–6 h group, P < 0.001; 8–12 h group, P = 0.016; 24–30 h group, P < 0.001; 48–60 h group, P < 0.001). There was a statistically significant difference between the HI 0–6 and 8–12 h groups (P = 0.007) and between HI 8–12 and 24–30, or 48–60 h groups (P < 0.001, P = 0.021, respectively), but no statistically significant difference was observed between other groups.
Figure 3. 1H-MRS images of newborn piglets at different times points after HI injury. (A–C) and (D) represented four different time points after HI operation: 0–6, 8–12, 24–30, and 48–60 h. In the spectra, the Glu peak was elevated at 0–6 h (A), then slightly depressed at 8–12 h (B), and thereafter elevated again at 24–30 h (C), and 48–60 h (D). However, the Gln peak did not show significant changes in the above image. Glu, glutamate; Gln, glutamine; NAA, N-acetylaspartate; Lac, lactate; Cho, choline; Cr, creatine.
Figure 4. Quantification results measured by 1H-MRS in control and model piglets at different time-points after HI insult, (A) Glu, (B) Gln, and (C) Glx. Errors bars indicate standard deviation values. *Compared with the control group, P < 0.05. Glu, glutamate; Gln, glutamine; Glx, glutamate + glutamine.
Figure 4B demonstrated that Gln level tended to increase slightly and then decrease after HI, but there were no statistically significant differences between the different groups.
Figure 4C showed that there was a similar change in Glx level as in Glu level, where the difference between the control group and the HI 24–30 h group was statistically significant (P = 0.009).
IVIM Data Processing
Figure 5 shows the changes in D, D*, and f in the control group and model group as detected by IVIM. After HI, D in the basal ganglia was markedly decreased at 0–6 h (compared with the control group, P < 0.001) and then gradually recovered over time, but it was still slightly lower than that in the control group (no statistically significant differences between different HI time point subgroups and the control group, P > 0.05).
Figure 5. D (A), D* (B), and f (C) parameters measured by IVIM scanning in control (black circles) and model piglets at different time-points, 0–6 h (black squares); 8–12 h (black regular triangles); 24–30 h (black inverted triangles); and 48–60 h (black diamonds). The horizontal lines in the scatter plots represent medians.
Figure 5B showed the changes in D* in the different groups. After HI, D* was decreased at 0–6 and 8–12 h, then increased at 24–30 h and thereafter decreased again at 48–60 h. However, there were no statistical differences in D* between the different groups.
As shown in Figure 5C, f was clearly decreased at 0–6 h after HI and then began to recover at 8–12 h, but it was still lower than that in the control group; at 24–30 and 48–60 h, f continued to decrease again. The difference in f was statistically significant between the control group and the HI model group at different time points (P < 0.001, for both). There was a statistically significant difference between 0–6 and 8–12 h groups and between 8–12 and 24–30, or 48–60 h groups (P < 0.001, for both). The specific data for D, D*, and f are shown in Table 1.
Analysis of Correlations between Metabolites of Glu Metabolism and IVIM-Derived Perfusion Parameters, D* and f
After HI, the Glu concentration in the basal ganglia showed a significant negative correlation with f (r = −0.643, P = 0.001; Figure 6A); similarly, there was an ordinary negative correlation between the Glx concentration and f (r = −0.478, P = 0.016; Figure 6C). However, no significant correlation was observed between the Gln concentration and f (Figure 6B). In addition, D* did not correlate with the concentration of Glu, Gln, or Glx (Figures 6D–F).
Figure 6. Scatter plots of IVIM-derived f (%) and D* (×10−3mm2/s) parameters over concentrations (mmol/kg) of Glu, Gln, and Glx. A significant negative correlation between f values and the Glu concentration was observed (A), whereas an ordinary negative correlation between f values and the Glx concentration was observed (C). However, there was no correlation between f values and the Gln concentration (B). Unlike for f values, no correlation between D* values and Glu (D), Gln (E), or Glx (F) concentrations was observed. Spearman's rank correlation coefficient r was calculated. Glu, glutamate; Gln, glutamine; Glx, glutamate + glutamine.
In the present study, we established a HI newborn piglet model and then investigated the level of Glu metabolism and the microcirculatory perfusion changes in the brain by multi-sequence MRI. The results showed that the change in Glu level was closely related to the microcirculatory perfusion level in the acute stage of HIBD. Since the occurrence of HIBD is a result of the interactions and influences of multiple factors that contribute to an abnormal pathological environment in the brain. On one hand, the development of dysmetabolism in cerebral energy causes a massive extracellular accumulation of Glu. The resultant excitotoxicity is one of the key factors inducing neuronal injury and death. On the other hand, hemodynamic disorders of the brain can also lead to damaged brain cells. In the brain, the neurons are located nearby small blood vessels, so the neurotransmitters released by synapses can also regulate the functions of these vessels (Huang et al., 1994) and nerve activities are tightly coupled with the degree of vascular perfusion (Busija et al., 2007; Jackman and Iadecola, 2015). In this study, we performed a preliminary exploration of relevant pathological mechanisms by multi-sequence MRI.
Our study results showed that Glu level underwent a “two-phase” change after HI. The sharp rise in Glu level in the 0–6 h group can be explained by the following mechanism. Energy dysmetabolism in nerve cells occured in the early stage, and the neurons were then depolarized, which caused a substantial release of excitable Glu into the extracellular spaces (Phillis et al., 1994; D'souza and Slater, 1995; Guo et al., 2011; Rocha-Ferreira and Hristova, 2016). Meanwhile, increased NADH/NAD+ ratio may directly promote the de novo synthesis of Glu (Ottersen et al., 1996), and the Glu level in astrocytes was then increased accordingly. In the 8–12 h group, the Glu concentration reached a transiently low peak due to the work of Glu transporters (GluTs), responsible for taking up extracellular Glu into astrocytes. On this basis, excessive Glu is transformed into Gln and then returns back to neurons, completing the “Glu (in neuron)–Gln (in astrocyte)” cycle (Bacci et al., 2002). Excessive Glu then suppresses the excessive excitation of basal ganglia, in a possible self-protective mechanism within human astrocytes (Sofroniew and Vinters, 2010; Ouyang et al., 2014). In the 24–30 h group, the Glu concentration increased again due to two causes: cell rupture caused by reperfusion injury led to an increase in Glu release; and the Glu reuptake mechanism was damaged in the late stages of the disease's course and thus resulted in the massive extracellular accumulation of Glu (Matsumoto et al., 1996). Meanwhile, the severe shortage of ATP inhibited the activity of glutamine synthetase (which is responsible for transforming Glu into non-excitatory Gln in the Glu-Gln cycle) (Dao et al., 1991), thus causing the accumulation of residual Glu in astrocytes (Torp et al., 1993).
Moreover, we did not find statistically significant differences in Gln concentration between the different groups in the present study, and this was perhaps associated with the PRESS sequence (TE = 37 ms) used for the 1H-MRS scan. A PRESS sequence (TE = 37 ms) is not the best choice to study Gln. Henry et al. (2011) believed that two-dimensional (2D) J-resolved MRS could completely disperse the peak overlapping on the 2D plane, which is caused by the approximate coupling constants in a one-dimensional plane, and had good stability and no susceptibility to linewidth fluctuation. Therefore, 2D J-resolved MRS is the optimal spectrum study method available for Gln quantification. In the future, we may be able to further investigate and explore Gln by 2D J-resolved MRS.
However, changes in local microcirculatory perfusion in the brain and their association with Glu level have not yet been fully clarified. In this study, we quantitatively evaluated changes in microcirculatory perfusion after HIBD using IVIM. As is well known, water molecules flow with blood (except for Brownian motion) within the capillary network of the microcirculation, and thus capillary blood flow can be regarded as another form of water diffusion. At the voxel level, the water flowing in randomly oriented capillaries can be regarded as an irregular random motion (Le Bihan et al., 1988) known as “pseudo-diffusion” due to the pseudo random organ distribution of the capillary network, and is related to the structure and blood flow rate of the capillary network. In the present study, perfusion was defined as the incoherent motion of water molecules in capillaries at the voxel level (Hu et al., 2015). The analysis of multiple b-value diffusion-weighted imaging (DWI) with an IVIM model is applicable to the quantification of two motion components, including water molecule diffusion and microcirculatory perfusion. Three parameters (D, D*, and f) are then finally calculated using a Levenberg–Marquardt non-linear least square fitting routine, of which D* and f both provide information on microcirculatory perfusion. The perfusion fraction f represents the volume percentage of perfusion-related diffusion effect in total diffusion effect, and it is positively correlated with the cerebral blood volume of the brain (Federau et al., 2014a,b); the pseudo-diffusion coefficient D*, i.e., the capillary perfusion-related diffusion coefficient, is determined depending on the blood flow rate and geometry of capillaries. Therefore, the changes in f and D* can be regarded as the changes in microcirculatory perfusion.
This study showed that after HI, f decreased at 0–6 h and then transiently recovered at 8–12 h but was still lower than that in the control group, and thereafter, it continued to decrease (see Table 1 and Figure 5C). Several investigators (Pulsinelli et al., 1982; Qiao et al., 2004; Ohshima et al., 2012) have measured changes in local perfusion in the HI brain of animal models using 4-iodo-[14C]-antipyrine, ASL and laser speckle flowmetry (LSF), respectively. They found a transient recovery in cerebral perfusion volume, followed by a decrease in perfusion, which is consistent with our study results. A transient increase in cerebral perfusion volume in the early stages of HIBD may be attributed to the activation of a cerebral vascular self-regulatory mechanism (Grant et al., 2005) (a self-protective mechanism of the body) after the recovery of blood reperfusion in the bilateral common carotid arteries. Blood perfusion decreased in the late stages, after HI. Such delayed hypoperfusion in the clinic is controversial, but most evidence suggests that it is closely related to secondary brain damage after HI (Jensen et al., 2006).
In the present study, the Glu concentration in the basal ganglia of newborn piglets significantly correlated with f in a negative manner (r = −0.643, P = 0.001). At 8–12 h after HI, f increased transiently (Table 1 and Figure 5C), corresponding to a transient low peak of Glu release (see Figure 4A), which indicates that the excessively accumulated Glu was cleared due to the recovery of perfusion in the early stages after HI. (Yamaguchi, 1977) The majority of current studies (Busija et al., 2007; Longo and Goyal, 2013) have shown that Glu and its synthetic analogs can dilate small arteries and veins in the brain and thus increase f. The Glu concentration increased while f decreased again, perhaps because delayed hypoperfusion can mediate the release of this excitatory neurotransmitter (Matsumoto et al., 1996). In turn, the abovementioned vasodilatory effect of Glu may be weakened with the progression of disease, and Glu may even constrict small arteries in the cerebral pia mater. As the Glu concentration increased, cerebral vasospasms may have become aggravated and may have even induced the breakage of postcapillary venules, thus causing severe microcirculatory disturbances (Huang et al., 1994), decreasing f further. Based on this finding, we postulate that the disturbances in microcirculatory perfusion and Glu release may result from, and contribute to, each other, and both may induce neuronal injury following HI.
However, our study results showed no statistically significant difference in D*. A possible reason may be the big limitation of D*, which results from uncertainty and very poor reproducibility (Wu et al., 2015; Nougaret et al., 2016; Yang et al., 2016). We found that D* displayed a large degree of dispersion during post-processing computations, maybe because of a high variation in D* within brain tissues due to the blood–brain barrier; f usually shows a small value in brain tissues. That is to say, D* is more likely to have an error when the proportion of perfusion is lower (King et al., 1992; Bisdas et al., 2013). In addition, the movements of animals and the partial volume effect of CSF during scanning can influence the measuring accuracy of D*. Some investigators have proposed that the accurate measurement of D* was also affected by the cardiac cycle, with the measurement value of D* significantly higher in the systole than in the diastole (Federau et al., 2013; Xu et al., 2016). Relative to D*, f was less affected by these physiological factors, had lower noise and was more uniform.
In addition, this study also demonstrated that the diffusion coefficient D, a parameter reflecting the true diffusion motion of water molecules, was markedly decreased at 0–6 h after HI and then gradually recovered over time, but was still slightly lower than that in the control group (see Figure 5A). This was perhaps due to increased glycolysis in the tissues in the early stage of HI, with more and more lactate accumulated, resulting in intracellular acidosis and cytotoxic edema. This subsequently led to the reduction of extracellular spaces (Tuor et al., 2014), which limited the diffusion of water molecules and decreased D. Subsequently, cellular swelling caused by cytotoxic edema compressed the capillaries and then resulted in the further hypoxia of brain tissues. Such hypoxia acted on vascular endothelial cells to increase vascular permeability, resulting in vascular edema. As a result, water molecules were retained extracellularly, the extracellular space became enlarged (Wang H. et al., 2012), the diffusion of water molecules was enhanced, and D increased. Conventional DWI reflects the microstructural changes in brain tissues with changes of the apparent diffusion coefficient (ADC). However, diffusions in biotic tissues actually measured by ADC include both the true diffusion of water molecules and capillary perfusion effects, which can be distinguished by IVIM (Le Bihan et al., 1986, 1988). Therefore, multiple b-value DWI, based on an IVIM bi-exponential model, can provide a value of D without the effects of perfusion factors, and thus can more accurately measure the diffusion of water molecules. In short, the early pathological changes in brain tissues after HI are indicated by changes in D.
Of course, our study had several limitations. Firstly, a small sample size was used, which may have caused a bias in results. However, we hope our results lay the basis for further, larger studies. Secondly, the coil diameter used in this study was relatively bigger than the head of the newborn piglets, and the SNR of images was not high enough. Therefore, work must focus on seeking a better SNR for future studies. Finally, the IVIM sequence was scanned with 16 b-values (including 9 b-values within a range of 0–200 s/mm2) in the present study, which took a long time, though more accurate and detailed perfusion information was provided; the motion displacement during scanning may have affected the accuracy of results. However, in order to overcome this limitation, data with severe motion was excluded.
In this study, we evaluated the changes in Glu metabolism and microcirculatory perfusion in the HI brain by 1H-MRS combined with IVIM. Glu concentration was increased in the early stage after HI, then transiently recovered and finally increased again, showing a “two-phase” change; perfusion-related parameter f showed a clear decrease, then transiently and slightly recovered, and thereafter continued to decrease. There was a significant negative correlation between the two parameters. Our data highlight the potential of combining changes in Glu concentration and f to explore the close relationship between cerebral dysmetabolism and microcirculatory disturbance after HI.
XW and YD participated in conceiving and designing of the idea. KS contributed to providing the post-processing assistant. YD and KS contributed to analyzing the experiment results. YD also contributed to drafting and editing of the manuscript. XW and YD also participated in critically revising the paper. All authors have read and approved the final manuscript for publication.
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (grant no. 30570541, 30770632, 81271631) and Outstanding Scientific Fund of Shengjing Hospital (item no. 201402).
Conflict of Interest Statement
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.
We would like to thank our colleagues in the Department of Radiology, Shengjing Hospital of China Medical University and MR Engineers whose efforts make this work possible.
Alexander, G. E., and Crutcher, M. D. (1990). Functional architecture of basal ganglia circuits: neural substrates of parallel processing. Trends Neurosci. 13, 266–271. doi: 10.1016/0166-2236(90)90107-L
Bacci, A., Sancini, G., Verderio, C., Armano, S., Pravettoni, E., Fesce, R., et al. (2002). Block of glutamate-glutamine cycle between astrocytes and neurons inhibits epileptiform activity in hippocampus. J. Neurophysiol. 88, 2302–2310. doi: 10.1152/jn.00665.2001
Bisdas, S., Braun, C., Skardelly, M., Schittenhelm, J., Teo, T. H., Thng, C. H., et al. (2014). Correlative assessment of tumor microcirculation using contrast-enhanced perfusion MRI and intravoxel incoherent motion diffusion-weighted MRI: is there a link between them? NMR Biomed. 27, 1184–1191. doi: 10.1002/nbm.3172
Bisdas, S., Koh, T. S., Roder, C., Braun, C., Schittenhelm, J., Ernemann, U., et al. (2013). Intravoxel incoherent motion diffusion-weighted MR imaging of gliomas: feasibility of the method and initial results. Neuroradiology 55, 1189–1196. doi: 10.1007/s00234-013-1229-7
Busija, D. W., Bari, F., Domoki, F., and Louis, T. (2007). Mechanisms involved in the cerebrovascular dilator effects of N-methyl-d-aspartate in cerebral cortex. Brain Res. Rev. 56, 89–100. doi: 10.1016/j.brainresrev.2007.05.011
Cooper, A. J., and Jeitner, T. M. (2016). Central role of glutamate metabolism in the maintenance of nitrogen homeostasis in normal and hyperammonemic brain. Biomolecules 6, pii: E16. doi: 10.3390/biom6020016
Cui, Y., Zhang, C., Li, X., Liu, H., Yin, B., Xu, T., et al. (2015). Intravoxel incoherent motion diffusion-weighted magnetic resonance imaging for monitoring the early response to ZD6474 from nasopharyngeal carcinoma in nude mouse. Sci. Rep. 5:16389. doi: 10.1038/srep16389
D'souza, S. W., and Slater, P. (1995). Excitatory amino acids in neonatal brain: contributions to pathology and therapeutic strategies. Arch. Dis. Child. Fetal Neonatal Ed. 72, F147–F150. doi: 10.1136/fn.72.3.F147
Federau, C., Hagmann, P., Maeder, P., Muller, M., Meuli, R., Stuber, M., et al. (2013). Dependence of brain intravoxel incoherent motion perfusion parameters on the cardiac cycle. PLoS ONE 8:e72856. doi: 10.1371/journal.pone.0072856
Federau, C., Maeder, P., O'Brien, K., Browaeys, P., Meuli, R., and Hagmann, P. (2012). Quantitative measurement of brain perfusion with intravoxel incoherent motion MR imaging. Radiology 265, 874–881. doi: 10.1148/radiol.12120584
Federau, C., O'Brien, K., Meuli, R., Hagmann, P., and Maeder, P. (2014a). Measuring brain perfusion with intravoxel incoherent motion (IVIM): initial clinical experience. J. Magn. Reson. Imaging 39, 624–632. doi: 10.1002/jmri.24195
Federau, C., Sumer, S., Becce, F., Maeder, P., O'Brien, K., Meuli, R., et al. (2014b). Intravoxel incoherent motion perfusion imaging in acute stroke: initial clinical experience. Neuroradiology 56, 629–635. doi: 10.1007/s00234-014-1370-y
Fujima, N., Yoshida, D., Sakashita, T., Homma, A., Tsukahara, A., Tha, K. K., et al. (2014). Intravoxel incoherent motion diffusion-weighted imaging in head and neck squamous cell carcinoma: assessment of perfusion-related parameters compared to dynamic contrast-enhanced MRI. Magn. Reson. Imaging 32, 1206–1213. doi: 10.1016/j.mri.2014.08.009
Grant, D. A., Franzini, C., Wild, J., Eede, K. J., and Walker, A. M. (2005). Autoregulation of the cerebral circulation during sleep in newborn lambs. J. Physiol. 564, 923–930. doi: 10.1113/jphysiol.2005.083352
Hagberg, H., Andersson, P., Kjellmer, I., Thiringer, K., and Thordstein, M. (1987). Extracellular overflow of glutamate, aspartate, GABA and taurine in the cortex and basal ganglia of fetal lambs during hypoxia-ischemia. Neurosci. Lett. 78, 311–317. doi: 10.1016/0304-3940(87)90379-X
Hauser, T., Essig, M., Jensen, A., Laun, F. B., Munter, M., Maier-Hein, K. H., et al. (2014). Prediction of treatment response in head and neck carcinomas using IVIM-DWI: Evaluation of lymph node metastasis. Eur. J. Radiol. 83, 783–787. doi: 10.1016/j.ejrad.2014.02.013
Henry, M. E., Lauriat, T. L., Shanahan, M., Renshaw, P. F., and Jensen, J. E. (2011). Accuracy and stability of measuring GABA, glutamate, and glutamine by proton magnetic resonance spectroscopy: a phantom study at 4 Tesla. J. Magn. Reson. 208, 210–218. doi: 10.1016/j.jmr.2010.11.003
Howlett, J. A., Northington, F. J., Gilmore, M. M., Tekes, A., Huisman, T. A., Parkinson, C., et al. (2013). Cerebrovascular autoregulation and neurologic injury in neonatal hypoxic-ischemic encephalopathy. Pediatr. Res. 74, 525–535. doi: 10.1038/pr.2013.132
Hu, L. B., Hong, N., and Zhu, W. Z. (2015). Quantitative measurement of cerebral perfusion with intravoxel incoherent motion in acute ischemia stroke: initial clinical experience. Chin. Med. J. 128, 2565–2569. doi: 10.4103/0366-6999.166033
Hu, Y. C., Yan, L. F., Wu, L., Du, P., Chen, B. Y., Wang, L., et al. (2014). Intravoxel incoherent motion diffusion-weighted MR imaging of gliomas: efficacy in preoperative grading. Sci. Rep. 4:7208. doi: 10.1038/srep07208
Jensen, E. C., Bennet, L., Hunter, C. J., Power, G. C., and Gunn, A. J. (2006). Post-hypoxic hypoperfusion is associated with suppression of cerebral metabolism and increased tissue oxygenation in near-term fetal sheep. J. Physiol. 572, 131–139. doi: 10.1113/jphysiol.2005.100768
Le Bihan, D., Breton, E., Lallemand, D., Aubin, M. L., Vignaud, J., and Laval-Jeantet, M. (1988). Separation of diffusion and perfusion in intravoxel incoherent motion MR imaging. Radiology 168, 497–505. doi: 10.1148/radiology.168.2.3393671
Le Bihan, D., Breton, E., Lallemand, D., Grenier, P., Cabanis, E., and Laval-Jeantet, M. (1986). MR imaging of intravoxel incoherent motions: application to diffusion and perfusion in neurologic disorders. Radiology 161, 401–407. doi: 10.1148/radiology.161.2.3763909
Le Bihan, D., Turner, R., and Macfall, J. R. (1989). Effects of intravoxel incoherent motions (IVIM) in steady-state free precession (SSFP) imaging: application to molecular diffusion imaging. Magn. Reson. Med. 10, 324–337. doi: 10.1002/mrm.1910100305
Longo, L. D., and Goyal, R. (2013). Cerebral artery signal transduction mechanisms: developmental changes in dynamics and Ca2+ sensitivity. Curr. Vasc. Pharmacol. 11, 655–711. doi: 10.2174/1570161111311050008
Martin, L. J., Brambrink, A. M., Lehmann, C., Portera-Cailliau, C., Koehler, R., Rothstein, J., et al. (1997). Hypoxia-ischemia causes abnormalities in glutamate transporters and death of astroglia and neurons in newborn striatum. Ann. Neurol. 42, 335–348. doi: 10.1002/ana.410420310
Massaro, A. N., Govindan, R. B., Vezina, G., Chang, T., Andescavage, N. N., Wang, Y., et al. (2015). Impaired cerebral autoregulation and brain injury in newborns with hypoxic-ischemic encephalopathy treated with hypothermia. J. Neurophysiol. 114, 818–824. doi: 10.1152/jn.00353.2015
Matsumoto, K., Lo, E. H., Pierce, A. R., Halpern, E. F., and Newcomb, R. (1996). Secondary elevation of extracellular neurotransmitter amino acids in the reperfusion phase following focal cerebral ischemia. J. Cereb. Blood Flow Metab. 16, 114–124. doi: 10.1097/00004647-199601000-00014
Nougaret, S., Vargas, H. A., Lakhman, Y., Sudre, R., Do, R. K., Bibeau, F., et al. (2016). Intravoxel incoherent motion-derived histogram metrics for assessment of response after combined chemotherapy and radiation therapy in rectal cancer: initial experience and comparison between single-section and volumetric analyses. Radiology 280, 446–454. doi: 10.1148/radiol.2016150702
Ohshima, M., Tsuji, M., Taguchi, A., Kasahara, Y., and Ikeda, T. (2012). Cerebral blood flow during reperfusion predicts later brain damage in a mouse and a rat model of neonatal hypoxic-ischemic encephalopathy. Exp. Neurol. 233, 481–489. doi: 10.1016/j.expneurol.2011.11.025
Ottersen, O. P., Laake, J. H., Reichelt, W., Haug, F. M., and Torp, R. (1996). Ischemic disruption of glutamate homeostasis in brain: quantitative immunocytochemical analyses. J. Chem. Neuroanat. 12, 1–14. doi: 10.1016/S0891-0618(96)00178-0
Ouyang, Y. B., Xu, L., Liu, S., and Giffard, R. G. (2014). Role of astrocytes in delayed neuronal death: GLT-1 and its novel regulation by microRNAs. Adv. Neurobiol. 11, 171–188. doi: 10.1007/978-3-319-08894-5_9
Pfeuffer, J., Provencher, S. W., and Gruetter, R. (1999). Water diffusion in rat brain in vivo as detected at very large b values is multicompartmental. MAGMA 8, 98–108. doi: 10.1016/s1352-8661(99)00013-7
Phillis, J. W., Smith-Barbour, M., Perkins, L. M., and O'regan, M. H. (1994). Characterization of glutamate, aspartate, and GABA release from ischemic rat cerebral cortex. Brain Res. Bull. 34, 457–466. doi: 10.1016/0361-9230(94)90019-1
Pulsinelli, W. A., Levy, D. E., and Duffy, T. E. (1982). Regional cerebral blood flow and glucose metabolism following transient forebrain ischemia. Ann. Neurol. 11, 499–502. doi: 10.1002/ana.410110510
Qiao, M., Latta, P., Foniok, T., Buist, R., Meng, S., Tomanek, B., et al. (2004). Cerebral blood flow response to a hypoxic-ischemic insult differs in neonatal and juvenile rats. MAGMA 17, 117–124. doi: 10.1007/s10334-004-0058-4
Rego, A. C., Santos, M. S., and Oliveira, C. R. (1996). Oxidative stress, hypoxia, and ischemia-like conditions increase the release of endogenous amino acids by distinct mechanisms in cultured retinal cells. J. Neurochem. 66, 2506–2516. doi: 10.1046/j.1471-4159.1996.66062506.x
Suh, C. H., Kim, H. S., Lee, S. S., Kim, N., Yoon, H. M., Choi, C. G., et al. (2014). Atypical imaging features of primary central nervous system lymphoma that mimics glioblastoma: utility of intravoxel incoherent motion MR imaging. Radiology 272, 504–513. doi: 10.1148/radiol.14131895
Torp, R., Arvin, B., Le Peillet, E., Chapman, A. G., Ottersen, O. P., and Meldrum, B. S. (1993). Effect of ischaemia and reperfusion on the extra- and intracellular distribution of glutamate, glutamine, aspartate and GABA in the rat hippocampus, with a note on the effect of the sodium channel blocker BW1003C87. Exp. Brain Res. 96, 365–376. doi: 10.1007/BF00234106
Tuor, U. I., Morgunov, M., Sule, M., Qiao, M., Clark, D., Rushforth, D., et al. (2014). Cellular correlates of longitudinal diffusion tensor imaging of axonal degeneration following hypoxic-ischemic cerebral infarction in neonatal rats. Neuroimage Clin. 6, 32–42. doi: 10.1016/j.nicl.2014.08.003
Wang, H., Wang, X., and Guo, Q. (2012). The correlation between DTI parameters and levels of AQP-4 in the early phases of cerebral edema after hypoxic-ischemic/reperfusion injury in piglets. Pediatr. Radiol. 42, 992–999. doi: 10.1007/s00247-012-2373-7
Wang, X. Y., Wang, H. W., Fu, X. H., Zhang, W. Q., Wu, X. Y., Guo, Q. Y., et al. (2012). Expression of N-methyl-d-aspartate receptor 1 and its phosphorylated state in basal ganglia of a neonatal piglet hypoxic-ischemic brain injury model: a controlled study of (1)H MRS. Eur. J. Paediatr. Neurol. 16, 492–500. doi: 10.1016/j.ejpn.2012.01.005
Wisnowski, J. L., Bluml, S., Paquette, L., Zelinski, E., Nelson, M. D. Jr., Painter, M. J., et al. (2013). Altered glutamatergic metabolism associated with punctate white matter lesions in preterm infants. PLoS ONE 8:e56880. doi: 10.1371/journal.pone.0056880
Wu, W. C., Chen, Y. F., Tseng, H. M., Yang, S. C., and My, P. C. (2015). Caveat of measuring perfusion indexes using intravoxel incoherent motion magnetic resonance imaging in the human brain. Eur. Radiol. 25, 2485–2492. doi: 10.1007/s00330-015-3655-x
Xiao, Y., Pan, J., Chen, Y., Chen, Y., He, Z., and Zheng, X. (2015). Intravoxel incoherent motion-magnetic resonance imaging as an early predictor of treatment response to neoadjuvant chemotherapy in locoregionally advanced nasopharyngeal carcinoma. Medicine (Baltimore) 94:e973. doi: 10.1097/MD.0000000000000973
Xu, X. Q., Choi, Y. J., Sung, Y. S., Yoon, R. G., Jang, S. W., Park, J. E., et al. (2016). Intravoxel incoherent motion mr imaging in the head and neck: correlation with dynamic contrast-enhanced MR imaging and diffusion-weighted imaging. Korean J. Radiol. 17, 641–649. doi: 10.3348/kjr.2016.17.5.641
Yang, S. H., Lin, J., Lu, F., Han, Z. H., Fu, C. X., Lv, P., et al. (2016). Evaluation of antiangiogenic and antiproliferative effects of sorafenib by sequential histology and intravoxel incoherent motion diffusion-weighted imaging in an orthotopic hepatocellular carcinoma xenograft model. J. Magn. Reson. Imaging 45, 270–280. doi: 10.1002/jmri.25344.
Zhang, Y. F., Wang, X. Y., Guo, F., Burns, K., Guo, Q. Y., and Wang, X. M. (2012). Simultaneously changes in striatum dopaminergic and glutamatergic parameters following hypoxic-ischemic neuronal injury in newborn piglets. Eur. J. Paediatr. Neurol. 16, 271–278. doi: 10.1016/j.ejpn.2011.05.010
Keywords: hypoxic–ischemic brain damage, 1H-MRS, IVIM, glutamate, perfusion
Citation: Dang Y-x, Shi K-n and Wang X-m (2017) Early Changes in Glutamate Metabolism and Perfusion in Basal Ganglia following Hypoxia-Ischemia in Neonatal Piglets: A Multi-Sequence 3.0T MR Study. Front. Physiol. 8:237. doi: 10.3389/fphys.2017.00237
Received: 30 January 2017; Accepted: 05 April 2017;
Published: 25 April 2017.
Edited by:Tijana Bojić, Institut za Nuklearne Nauke Vinča, University of Belgrade, Serbia
Copyright © 2017 Dang, Shi and Wang. 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) or licensor 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: Xiao-ming Wang, firstname.lastname@example.org