A Volumetric Method for Quantification of Cerebral Vasospasm in a Murine Model of Subarachnoid Hemorrhage

Summary

The goal of this article is to present a method that allows a 3-dimensional reconstruction of the cerebrovascular tree in mice after micro computed tomography and determination of volumes of entire vessel segments that can be used to quantify cerebral vasospasm in murine models of subarachnoid hemorrhage.

Abstract

Subarachnoid hemorrhage (SAH) is a subtype of hemorrhagic stroke. Cerebral vasospasm that occurs in the aftermath of the bleeding is an important factor determining patient outcome and is therefore frequently taken as a study endpoint. However, in small animal studies on SAH, quantification of cerebral vasospasm is a major challenge. Here, an ex vivo method is presented that allows quantification of volumes of entire vessel segments, which can be used as an objective measure to quantify cerebral vasospasm. In a first step, endovascular casting of the cerebral vasculature is performed using a radiopaque casting agent. Then, cross-sectional imaging data are acquired by micro computed tomography. The final step involves 3-dimensional reconstruction of the virtual vascular tree, followed by an algorithm to calculate center lines and volumes of the selected vessel segments. The method resulted in a highly accurate virtual reconstruction of the cerebrovascular tree shown by a diameter-based comparison of anatomical samples with their virtual reconstructions. Compared with vessel diameters alone, the vessel volumes highlight the differences between vasospastic and non-vasospastic vessels shown in a series of SAH and sham-operated mice.

Introduction

Aneurysmatic subarachnoid hemorrhage (SAH), a subtype of hemorrhagic stroke, is a common disease in neurointensive care units. Besides early brain injury (EBI), which comprises cerebral damage caused by the bleeding event itself, another important factor determining patient outcome is delayed cerebral ischemia (DCI), defined by clinical deterioration through impaired cerebral perfusion or cerebral infarction not associated with interventional or surgical procedures^1,2,3. Important mechanisms contributing to DCI are vasospasms of large cerebral vessels on the one hand; on the other hand, microcirculatory dysfunction with vasospasm of microvessels and microthrombosis, and ischemia related to cortical spreading depressions play a role (reviewed in Madonald 2014¹). Therefore, diagnosis of vasospasm of the large cerebral vessels is crucial in clinical practice and displays an important endpoint in many clinical and experimental studies.

Despite the fact that the features of vasospasm in murine SAH models are not directly transferable to the human patient, murine models of SAH related vasospasm have been of growing significance in the last years. In these models SAH is induced by endovascular filament perforation^4,5,6,7,8, transection of cisternal vessels⁹, or injection of blood into the CSF^10,11,12. In contrast to large animal models of SAH which were traditionally designed to study vasospasm¹³, murine models have the great advantage that numerous transgenic mice strains are available. This makes them an excellent tool for studying molecular mechanisms leading to vasospasm and DCI. However, determination of cerebral vasospasm in mice is challenging. This is because in contrast to large animal models in which vasospasm can be examined using clinical imaging techniques, in vivo imaging to analyze cerebral vasospasm in mice is not yet available. Therefore, vasospasm is commonly determined using either histological sections^10,11 or microscopically after casting of cerebral vessels^7,9,12. However, these techniques have the disadvantage that vessel diameters are examined at defined points only.

Based on a previous study⁷, this manuscript presents a method for objective and reproducible analysis of vasospasm in a murine SAH model. The method is based on perfusion and casting of the cerebral vessels, ex vivo micro-CT scanning, digital reconstruction of the vessel tree, and subsequent evaluation of volumes of entire cerebral vessels.

Protocol

The animal experiments were approved by the responsible animal care committee (Landesuntersuchungsamt Rheinland-Pfalz) and carried out in accordance with the German Animal Welfare Act (TierSchG). All applicable international, national, and institutional guidelines for the care and use of animals were followed.

In this study, male C57BL6 mice (age 10-12 weeks) were used. In brief, subarachnoid hemorrhage was induced by endovascular filament perforation under anesthesia with isoflurane. The left external carotid artery was prepared surgically. Then, a filament was inserted into the external carotid artery and advanced intracranially through the internal carotid artery which was perforated at the carotid T, inducing a subarachnoid hemorrhage. A rise in intracranial pressure was taken as an indicator of successful endovascular perforation. A detailed protocol of an endovascular filament perforation model of SAH in mice has been published by others^8,14.

1. Perfusion and Endovascular Casting

1. In this study, perfusion was performed 72 hours after induction of SAH. Induce anesthesia by intraperitoneally injecting 5 µg/g body weight (bw) midazolam, 30 ng/g bw fentanyl, and 0.5 µg/g bw medetomidin. Continue only after a sufficient anesthesia level has been reached, which is confirmed by the absence of reactions to pain stimuli.
2. Open the thorax, puncture the left ventricle with a 21G cannula, open the right atrium, and clamp the descending aorta as described elsewhere¹⁵.
3. Perform a transcardiac perfusion using the following solutions: (i) Dulbecco's Phosphate Buffered Saline containing MgCl₂ and CaCl₂ at pH 7.4 with 1 g/L glucose, and (ii) 4% paraformaldehyde solution.
   1. Start the perfusion with solution (i) for 2 minutes and continue with solution (ii) for 4 minutes.
   2. Infuse the solutions at a temperature of 37 °C and using a pressure-controlled pump with variable perfusion rate to perfuse with a constant pressure of 70 mmHg, which was found to be the optimal perfusion pressure to analyze vasospasm in mice¹⁶. Avoid a loss of pressure when switching from solution (i) to solution (ii).
4. After perfusion with solutions (i) and (ii), continue perfusion for 20 minutes at room temperature with radiopaque casting agent (see Table of Materials) at a constant rate of 0.2 ml/min.
5. Allow for curing of the radiopaque casting material at 4 °C overnight. Then remove the brain from the skull as previously described¹⁷, transfer the sample to 4% PFA solution, and store the sample at 4 °C until micro CT scanning.

2. Micro Computed Tomography

1. Place the brain in the center of a plastic tube with a blunt anatomical forceps. Choose a tube with a slightly larger diameter than the sample to ensure that the object does not move during image acquisition. Use gauze to close the tube.
2. Attach the plastic tube to the micro-stepping motor on the computer-navigated-control (CNC) positioning system in the X-ray cabin, in which the object is rotated around its horizontal axis.
3. Align the sample in the field-of-view under X-ray radiography. To achieve the maximal magnification, place the object as close as possible to the x-ray source and maximize the distance to the detector as far as possible.
4. Use a step-and-shoot image acquisition protocol with the following scan parameters: set exposure time to 1 s for each projection to optimize the signal-to-noise ratio (SNR), tube voltage 80 kV (current 38 µA), 360° rotation resulting in 1,000 projections.
5. For reconstruction of RAW-Data use a filtered back projection algorithm applying the Shepp-Logan filter with a matrix of 1024 x 1024 x 1024 voxels using reconstruction software (see Table of Materials). For further analysis import the resulting DICOM data into 3D visualization software (see Table of Materials).

3. 3-dimensional Reconstruction of the Intracranial Vascular Tree and Determination of Vessel Volumes

Note: Background information on the functions of the visualization software can be found using the function help.

1. Import Dicom data into visualization software using the function Import.
2. Visualize the vessel tree with the function Volren. Choose the visualization threshold so that the large cerebral arteries are depicted in sharp outlines. It is important to use the same visualization threshold for all data sets belonging to the experimental series.
3. Virtually dissect the basal cerebral arteries (Circle of Willis) with the function VolumeEdit by surrounding the vessels with the cursor. Then virtually dissect the vessel segment to be analyzed. Therefore, rotate the 3-dimensional model of the vascular tree in order to precisely separate all small branches from the main artery. It is essential for the further analysis to delete all vessels except the vessel segment to be analyzed.
4. Apply the function Autoskeleton with threshold set to visualization threshold, which generates a center-line-based SpatialGraph.
5. Then, apply the function SpatialGraphToLineSet to create a line set. Divide the line set into its single subsegments by manually choosing the single subsegments with the cursor and clicking on "split". This step is crucial in order to calculate the volumes of the single subsegments.
6. Use the function LineSetToSpatialGraph to create a Spatial Graph again.
7. Use the function SpatialGraphStatistics to determine length, volume, and diameter of each subsegment. For color-coded visualization representing the course of the vessel diameter, use the function SpatialGraphView. Set segment coloring to "thickness", which correlates with the vessel diameter. It is important to choose the same color map for all datasets belonging to the experimental series.
8. Add the lengths of the subsegments to determine which subsegments are to be included in the further analysis. In the present study we evaluated a vessel segment consisting of 1 mm of the internal carotid artery proximal of the carotid T and 2.5 mm of the middle cerebral artery distal of the carotid T. Then add the volumes to determine the vessel volume of the defined vessel segment.

Results

Virtual reconstruction of the 3-dimensional intracranial vascular tree

The 3-dimensionally reconstructed intracranial vascular trees provided a highly accurate vascular anatomy (Figure 1). To assess the accuracy, we performed a diameter-based comparison between vessel diameters determined microscopically and from the 3-dimensional virtual reconstructions at 2 anatomically defined points (1: left middle cerebral artery (MCA) 1 mm distal of the carotid T; 2: right MCA 1 mm distal of the carotid T). For microscopic determination of vessel diameters, a high-resolution camera (Infinity X-21, Deltapix) with DeltaPix Insight software version 2.0.1 calibrated to a micrometer scale was used. For this evaluation, 10 brain samples (5 SAH, 5 sham) were analyzed. These were from a series of 12 mice, in 7 of which an SAH was induced, while 5 underwent sham surgery (2 animals of the SAH group died on postoperative days 1 and 2, respectively). There were no significant differences between the diameters determined microscopically and virtually, indicating an accurate virtual reconstruction of the intracranial vascular anatomy (microscopic determination vs. virtual reconstruction, mean ± standard error: left MCA 150 ± 9 µm vs. left MCA 150 ± 8 µm; right MCA 153 ± 8 µm vs.154 ± 9 µm, see Figure 2).

Quantification of cerebral vasospasm in mice with SAH

To quantify cerebral vasospasm, (i) the volume of a predefined representative 3.5 mm vessel segment, consisting of 1 mm internal carotid artery (ICA) and 2.5 mm MCA on the left, and (ii) the vessel diameters at 2 anatomically defined points (left and right MCA) were determined in brain samples from the SAH and sham-operated animals (n = 5). The vessel volume was significantly lower in SAH compared to sham (36 ± 4 nL vs. 71 ± 9 nL, p < 0.05). Vessel diameters were lower in SAH compared to sham (left MCA: 140 ± 11 µm vs. 160 ± 10 µm, p = 0.11; right MCA: 130 ± 16 µm vs. 158 ± 13 µm, p < 0.05; see Figure 3), while the level of significance was only reached for analysis of the right MCA.

Figure 1. Virtual reconstruction of the intracranial vascular tree. (A) shows a representative brain sample; (B) shows the corresponding virtually reconstructed vascular tree. (C) Color-coded visualization of the diameter of the left MCA. Please click here to view a larger version of this figure.

Figure 2. Accuracy of the digital reconstruction of the vasculature. Mean diameters measured from 3D reconstructed brain samples compared with those determined microscopically. Data are shown as mean ± standard error of the mean. Please click here to view a larger version of this figure.

Figure 3. Vessel volume and vessel diameter after SAH. (A) Comparison of MCA diameters measured from 3D reconstructed vasculature in SAH and sham mice. (B) Vessel volume in SAH and sham mice. Data are shown as mean ± standard error of the mean. *p < 0.05. Please click here to view a larger version of this figure.

Discussion

Murine SAH models are an important tool for basic SAH research. Cerebral vasospasm is frequently used as an endpoint in experimental studies investigating the mechanisms leading to DCI after SAH^9,11. However, quantification of cerebral vasospasm in mice or other small animal models of SAH is challenging. Commonly, vasospasm is quantified by ex vivo determination of vessel diameters at defined anatomical points after endovascular perfusion and casting^7,9,12 or by determination of the circumference of defined vessels on histological sections^10,11. However, these methods have some disadvantages: Vasospasm is evaluated only at defined anatomical points; vasospasm of neighboring vessel segments may escape evaluation. Histological artifacts present another source of errors. Furthermore, the evaluation can be rather subjective, because the exact position where the vessel diameter is measured is determined by the investigator.

The goal was therefore to establish a method that quantifies cerebral vasospasm by calculating the vessel volume of entire cerebral vessel segments from cross-sectional imaging data⁷. The most important advantage of the volumetric method presented here is that entire vessel segments can be examined. This avoids the need of definition of a point where the vessel diameter is measured. A further advantage of the evaluation of whole vessel segments is that it presumably presents a more objective parameter to quantify vasospasm than determination of vessel diameters at defined points where vasospasm of the more proximal or distal vessel may escape evaluation. Digital representation of vessel diameters using a color code allows an intuitive estimation of the degree of vasospasm. Furthermore, volumetric evaluation leads to larger differences between vasospastic vessels compared to evaluation of vessel diameters as shown in the representative results. The virtual reconstruction achieved with the method presented here reflects the vascular anatomy accurately. This is shown by the evaluation of the representative series, in which vessel diameters measured microscopically and from the digital reconstructions were similar, reproducing the observations of a previous study⁷. However, despite its advantages, further studies are needed to evaluate whether or not the method presented here is superior to conventional methods of vasospasm analysis.

A limitation of the method presented here is that it affords more time compared to microscopic analysis of casted brain samples or histological analysis (micro CT scanning time 90 minutes per brain sample, data processing 45 min per brain sample). Furthermore, the availability of micro CT scanners may limit its application. The number of animals examined here was sufficient to demonstrate the feasibility of the protocol described in this manuscript. However, if the protocol should be used in treatment studies, animal numbers would have to be calculated based on the expected effects on vessel volumes and diameters. Another limitation of this and other studies using murine SAH models is that vasospasm is determined ex vivo. This makes longitudinal studies impossible that investigate baseline values before SAH induction and vasospasm at different time points. Although studies have demonstrated that it is possible to depict the anatomy of the large intracranial vessels of mice in vivo using magnetic resonance tomography¹⁸, computed tomography angiography¹⁹, or digital subtraction angiography²⁰, these methods, to our knowledge, have not yet been used to analyze cerebral vasospasm in murine SAH models in vivo. Of note, the digital reconstruction of the cerebral vasculature with subsequent volumetric evaluation of cerebral vasospasm presented here is not limited to the use on ex vivo micro CT data. If high resolution vascular cross sectional brain imaging in mice should become available in the future, it could be used to perform a volumetric analysis of vasospasm in vivo.

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
Medetomidin	Pfizer, Karlsruhe, Germany	n.a.
Midazolam	Ratiopharm, Ulm, Germany	n.a.
Fentanyl	Curamed, Karlsruhe, Germany	n.a.
Venofix 21G	B Braun Melsungen AG, Melsungen, Germany	n.a.	21G cannula
Dulbecco’s Phosphate Buffered Saline containing MgCl2 and CaCl2, pH 7.4	Sigma-Aldrich, Hamburg, Germany	D8662
4% paraformaldehyde solution	Sigma-Aldrich, Hamburg, Germany	100496
Microfil MV-122	Flowtech Inc., Carver, MA, USA	n.a.	Radiopaque
Micro-CT system Y.Fox	Yxlon, Garbsen, Germany	n.a.
Reconstruction Studio software version 1.2.8.1	TeraRecon, Frankfurt am Main, Germany	n.a.	Reconstruction software
Amira software version 5.4.2	FEI Visualization Sciences Group, Hillsboro, OR, USA	n.a.	Visualization software
PHD ultra syringe pump	Harvard Apparatus	70-3	Pressure controlled pump
anatomical forceps (blunt)	B Braun Melsungen AG, Melsungen, Germany	160323_v
Infinity X-21	Deltapix, Maalov, Denmark	n.a.	high resolution camera
DeltaPix Insight software version 2.0.1	Deltapix, Maalov, Denmark	n.a.
C57BL6 mice	Charles River, Cologne, Germany