Jakub Maciej
Surmacki
ab,
Laura
Ansel-Bollepalli
ab,
Francesca
Pischiutta
c,
Elisa R.
Zanier
*c,
Ari
Ercole
*d and
Sarah Elizabeth
Bohndiek
*ab
aDepartment of Physics, University of Cambridge, JJ Thomson Avenue, Cambridge, CB3 0HE, UK. E-mail: js2238@cam.ac.uk; laura.bollepalli@cruk.cam.ac.uk; seb53@cam.ac.uk; Tel: +44 (0)1223 769697 Tel: +44 (0)1223 337267
bCancer Research UK Cambridge Institute, University of Cambridge, Robinson Way, Cambridge, CB2 0RE, UK. Tel: +44 (0)1223 769697 Tel: +44 (0)1223 337267
cDepartment of Neuroscience, IRCCS - Istituto de Ricerche Farmacologiche Mario Negri, Via G. La Masa 19, 20156 Milano, Italy. E-mail: francesca.pischiutta@marionegri.it; elisa.zanier@marionegri.it; Tel: +39 02 39 014 721 Tel: +39 02 39 014 204
dDivision of Anaesthesia, University of Cambridge, Addenbrooke's Hospital, Hills Road, Cambridge, CB2 0QQ, UK. E-mail: ae105@cam.ac.uk; Tel: +44 (0)1223 217889
First published on 28th November 2016
Traumatic brain injury (TBI) constitutes a major cause of death and long-term disability. At present, we lack methods to non-invasively track tissue biochemistry and hence select appropriate interventions for patients. We hypothesized that detailed label-free vibrational chemical analysis of focal TBI could provide such information. We assessed the early spatial and temporal changes in tissue biochemistry that are associated with brain injury in mice. Numerous differences were observed in the spectra of the contusion core and pericontusional tissue between 2 and 7 days. For example, a strong signal from haem was seen in the contusion core at 2 days due to haemorrhage, which subsequently resolved. More importantly, elevated cholesterol levels were demonstrated by 7 days, which may be a marker of important cell repair processes. Principal component analysis revealed an early ‘acute’ component dominated by haemorrhage and a delayed component reflecting changes in protein and lipid composition. Notably we demonstrated changes in Raman signature with time even in the contralateral hemisphere when compared to sham control mice. Raman spectroscopy therefore shows promise as a probe that is sensitive to important pathobiological processes in TBI and could be applied in future both in the experimental setting, as well as in the clinic.
Studies to date indicate that TBI is a disease governed by the host response to the injury, which varies dramatically between patients and has a complex spatiotemporal profile. Some commonalities include cell damage and loss of axonal connectivity over the first 24 hours and neuronal phospholipid membrane breakdown, clinically evidenced by early elevation in extracellular glycerol concentrations.1 A variety of molecular processes, such as excitotoxicity/calcium influx, oxidative injury and inflammation, along with cellular changes in resident microglial cells and recruitment of macrophages, may continue to drive injury.2–4 These insights have been provided by traditional experimental techniques, such as immunohistochemistry that must be performed ex vivo and has specificity to predetermined epitopes only. Cerebral microdialysis can provide some chemical information in vivo but has poor spatial resolution and is only sensitive to small molecules and these must be present in the extracellular fluid. Non-invasive tools to track brain tissue biochemistry after injury remains an unmet need, which if met, could shed light on disease pathobiology and furthermore, assist clinically with translation of precise and timely interventions, such as neuroprotective therapies.5
Raman spectroscopy (RS) is a powerful label-free analytical technique that has been widely used to study tissue biochemistry, for example, profiling lipid, protein and nucleic acid composition.6–8 All twelve major and minor human brain lipids can be identified using RS at 785 nm.9 Previous work has revealed that RS can indicate brain tissue damage in response to both radiation10 and penetrating11 injuries. In the latter study, RS was able to delineate injured and healthy brain regions in a mouse model of motor cortex needle injury on the basis of reduced signal intensity at 1660 cm−1, attributed to amide I vibration of proteins and CC stretching of lipids, and at 1002 cm−1, attributed to the ring breathing mode of phenylalanine in proteins.11 However, these data were acquired exclusively in the motor cortex and penetrating injury does not reflect the vast majority of TBI seen clinically. Furthermore, the study was conducted at a single time point immediately after injury, which does not provide crucial information on the evolution of the pathology.2–4
We hypothesized that RS could monitor the progression and resolution of brain injury over space and time in a clinically relevant model of focal TBI. To test this hypothesis, we studied for the first time using RS the spatial and temporal distribution of early changes in brain biochemistry after focal severe TBI modelled by controlled cortical impact in mice. We found that RS provides a sensitive readout of the evolution of tissue injury in our model.
Mice were sacrificed at 2 or 7 days after TBI (n = 3 per time point) under deep anesthesia (ketamine 30 mg/medetomidine 0.3 mg). The brains were then removed, washed once with phosphate buffered saline (PBS) then frozen immediately into liquid nitrogen and stored at −80 °C. Before commencing Raman analysis brains were thawed for 20 min at room temperature.
Whole brain samples were placed on a 25 mm round microscope coverslip in a petri dish and illuminated with 35 mW at 785 nm. Single Raman spectra were collected with 1 s integration time and 10 accumulations. At each position interrogated (as indicated on the figures), at least 10 spectra were collected along a line scan up to 70 μm length. Total numbers of recorded Raman spectra used in the analysis are: contralateral (CL; n(2 days) = 330, n(7 days) = 330, n(control) = 560); pericontusional (PC; n(2 days) = 330, n(7 days) = 180, n(control) = 240); contusion core (CC; n(2 days) = 140, n(7 days) = 240, n(control) = 360) and ipsilateral tissue distant from the contusion core (ID; n(2 days) = 120, n(7 days) = 120, n(control) = 200).
Data processing was performed using Project Plus Four 4.1 (WITec), Origin 2016 (OriginLab) and MATLAB with PLS-Toolbox (Mathworks). All Raman spectra were cosmic ray and baseline corrected then normalized to the area-under-curve. Statistical significance was analysed using two-way ANOVA, followed by Tukey's multiple comparison test. p < 0.05 was considered statistically significant.
To examine the changes recorded by RS in more detail, we computed the differences between the average spectra across all mice at each spatial location and time point (Fig. 3). No significant difference was observed between the CL hemisphere and ID measurements made at ‘uninjured’ tissue away from the contusion core (identified qualitatively from the bright field images) either at 2 or 7 days (Fig. 3A and B). Small differences are apparent when comparing contralateral control tissue between 2 days and 7 days, associated with the peaks at 1140, 1462, 1648, and 1660 cm−1 (Fig. 3C). Striking differences arise between the contralateral controls and the contusion core at 2 days (Fig. 3D). Some of these resolve with time, in particular the spectral shoulder between the peaks at 1562 and 1620 cm−1, which we ascribe to resonance excitation of haemoglobin,17 but the spectral signatures remain very different compared to contralateral also at 7 days (Fig. 3E). Comparing the contusion core between 2 days and 7 days (Fig. 3F), interesting spectral differences can be noted in bands at 426, 701, 801, 1128, 1301, 1440, 1462 and 1660 cm−1, associated with lipids. Peaks at 426, 675, 754, 1003, 1224, 1547, 1562, 1620 cm−1 are likely due to resonance excitation of haemoglobin;17 the extensive area of bleeding present at 2 days after TBI has largely resolved by 7 days resulting in a decreasing signal of haemoglobin (Fig. 1A and B). Pericontusional tissue showed similar changes to the contusion core.
To isolate the origin of the changing lipid signal observed, we acquired average Raman spectra (Fig. 4A) of cholesterol (Fig. 4B), cholesteryl nonanoate (Fig. 4C), cholesteryl oleyl carbonate (Fig. 4D), as well as a mixed lipid solution (intralipid) and protein albumin. The most intense bands from cholesterol are at 426, 546, 605, 701, 741, 1128, 1440, 1462 and 1670 cm−1; the band at 1440 cm−1 can be assigned to CH deformation and 1670 cm−1 to CC stretch vibrations.9 Additional bands at 1735 cm−1 (ester group) as well as 1029 cm−1 and 1301 cm−1 (fatty acid chains) are observed in the spectra of cholesterol derivatives. Characteristic bands of intralipid (phospholipid stabilized soybean oil) are at 718, 845, 879, 970, 1064, 1079, 1264, 1301, 1440, 1656 and 1745 cm−1. Symmetric and antisymmetric stretch vibrations of the choline group N+(CH3)3 in phospholipid9 can be found at 718 and 879 cm−1. Albumin is the most common α-helical protein with the Raman bands at 508, 852, 1003, 1093, 1206, 1339, 1450, 1586, 1604 and 1656 cm−1. Evidence of these Raman bands from cholesterol, cholesterol derivatives, phospholipids and protein are observed in our RS data from the brain samples.
We used the knowledge of these peak identifications to guide our statistical analysis to assess the most significant spatial and temporal changes in the RS data. Changes in lipid and protein metabolism are of particular interest biologically after TBI in relation to cell damage and inflammatory processes. Unfortunately, many lipid and protein Raman signals overlap which makes a clean interpretation difficult. Nevertheless, in an attempt to obtain some mechanistic insight, we performed two-way ANOVA tests using the extracted peak intensity ratios for I(701)/I(718) (cholesterol9/phospholipid;7,9Fig. 5A), I(1301)/I(1250) (mixed fatty acid chains7–10 and amide III protein;8,10Fig. 5B), I(1440)/I(1660) (CH deformation7–10/mixed amide I protein7,8,10 and CC stretching of lipids;7–10Fig. 5C), and I(1301)/I(1620) (mixed amide III protein8,10 and fatty acid chains7–10/haemoglobin;17Fig. 5D). All metrics showed a statistically significant difference between the contralateral and pericontusional/contusional regions at 7 days after TBI; all except I(1440)/I(1660) also showed changes at 2 days after TBI. I(1301)/I(1250) also showed a significant difference between the pericontusional and contusional locations. Interestingly a large effect was observed in the ratio I(701)/I(718) demonstrating a significant relative increase in cholesterol and decrease in phospholipid contributions.
In an attempt to determine a TBI ‘fingerprint’ we subsequently performed principal component analysis on the entire data set. The distribution of scores on principal components 2 and 3 (Fig. 6A) clearly distinguish the contusion and pericontusion regions from the distant contusional region and contralateral hemisphere, as well as the temporal evolution of the data at 2 and 7 days after TBI, based on the illustrated loadings (Fig. 6B). From Fig. 6A, principal component 2 is most effective at distinguishing injured and control tissue whereas component 3 discriminates better at 7 days. From Fig. 6B, PC2 is dominated by intense signals between 1547 and 1620 cm−1, representing resonance Raman shift from haemoglobin and so the acute phase is dominated by the contribution from bleeding at/near the contusion site. An important feature of PC3 the peaks at 701 cm−1 representing an increase in cholesterol signal with respect to the 718 cm−1 phospholipid peak.
One striking spectral feature results from resonant excitation of haemoglobin in the pericontusional and contusion core regions, consistent with haemorrhage resulting from vascular disruption and bleeding at the lesion.18 Haemorrhagic conversion and contusion expansion is a potentially devastating pathology after TBI, which can lead to oxidative damage and cytotoxicity19 so it is particularly important that any non-invasive monitoring technique is sensitive to this process. The haemoglobin signals are no longer apparent by 7 days; the macroscopic appearance of the tissue is consistent with the phagocytic clearance of erythrocytes and endocytosis of haemoglobin and haem by macrophages/microglia.20
A second striking spectral feature appears in the pronounced peaks at 701 cm−1, assigned to cholesterol and 718 cm−1, assigned to phospholipid. Our normalised spectra show a relative increase of cholesterol signal and decrease in phospholipid in the contusion core compared to the contralateral tissue and over time. The relative change at 701/718 cm−1 provides strong evidence of a change in lipid chemistry in the pericontusional tissue and contusion core at 2 days that remains at 7 days. Cholesterol is known to be critical for cellular repair, in neurodegenerative diseases21 and the apolipoproteins involved in cholesterol transport are up-regulated at 4 days after TBI.22 Furthermore, the ε4 allele of apolipoprotein E has been associated with worse outcomes after TBI23 and increased risk of post-TBI dementia.24 Interestingly, local cholesterol exposure may alter the accessibility of amyloid precursor protein (APP) to β-secretase (Bace1), an enzyme responsible for the cleavage of APP to Amyloid-β (Aβ).24,25 The accumulation of Aβ is recognised as a pathognomonic mechanism in Alzheimer's and can also be seen within hours of TBI.26 These considerations suggest that RS could in future provide a probe of these critical mechanisms in TBI pathogenesis.
Another spectral feature characteristic of the TBI appeared at 1440/1660 cm−1. This feature has been used as a marker of the relative concentration of β-sheet protein (amide I peak at 1660 cm−1) – a hallmark of amyloid Aβ.27 While it is hard to disentangle changes in lipid/protein content in regions of the spectrum around 1440 and 1660 cm−1, RS clearly identifies changes in lipid content with respect to protein, consistent with previous findings.11 A final feature of note is the increased heterogeneity of measurements between the contusion and pericontusional tissue reflecting the spatial heterogeneity and temporal evolution of the contusion.3,28 Further studies combining RS and conventional immunohistochemical analysis will allow us to understand the contribution of cell population activation state to the RS brain biochemistry.
There are some limitations to this study. Small differences were observed in the 1440/1462 cm−1 and 1648/1660 cm−1 protein/lipid complex between the spectra recorded in sham mice and contralateral region of injured mice, as well as between the contralateral region at 2 and 7 days. Further work is needed to determine which underlying biological changes could impact brain biochemistry far from the contusion core. Examples include excitotoxicity, inflammation, blood brain barrier breakdown, oedema or intracranial pressure. We are unable, however, from our present data to characterise the cellular population in the volume being probed; correlation with histology and immunohistochemistry over space and time is required to fully understand the pathological underpinnings of our findings. Nonetheless, since RS could be applied in vivo, these findings are promising for tracking the temporal and spatial evolution brain tissue biochemistry after injury.
This journal is © The Royal Society of Chemistry 2017 |