This work demonstrates the use of a multimodal ultrasound-based imaging platform for noninvasive imaging of ischemic stroke. This system allows for the quantification of blood oxygenation through photoacoustic imaging and impaired perfusion in the brain through acoustic angiography.
Presented here is an experimental ischemic stroke study using our newly developed noninvasive imaging system that integrates three acoustic-based imaging technologies: photoacoustic, ultrasound, and angiographic tomography (PAUSAT). Combining these three modalities helps acquire multi-spectral photoacoustic tomography (PAT) of the brain blood oxygenation, high-frequency ultrasound imaging of the brain tissue, and acoustic angiography of the cerebral blood perfusion. The multi-modal imaging platform allows the study of cerebral perfusion and oxygenation changes in the whole mouse brain after stroke. Two commonly used ischemic stroke models were evaluated: the permanent middle cerebral artery occlusion (pMCAO) model and the photothrombotic (PT) model. PAUSAT was used to image the same mouse brains before and after a stroke and quantitatively analyze both stroke models. This imaging system was able to clearly show the brain vascular changes after ischemic stroke, including significantly reduced blood perfusion and oxygenation in the stroke infarct region (ipsilateral) compared to the uninjured tissue (contralateral). The results were confirmed by both laser speckle contrast imaging and triphenyltetrazolium chloride (TTC) staining. Furthermore, stroke infarct volume in both stroke models was measured and validated by TTC staining as the ground truth. Through this study, we have demonstrated that PAUSAT can be a powerful tool in noninvasive and longitudinal preclinical studies of ischemic stroke.
Blood transports oxygen (via the hemoglobin protein) and other important nutrients to tissues in our bodies. When the flow of blood through tissues is interrupted (ischemia), severe damage to the tissues can occur, the most immediate effects of which are due to a lack of oxygen (hypoxia). Ischemic stroke is the result of interrupted blood flow to a certain region of the brain. The brain damage resulting from an ischemic stroke can occur within minutes of a vessel blockage, and can often have debilitating and lasting effects1,2. A highly valuable strategy to evaluate the physiopathology after ischemic stroke and identify and test new treatments is the use of small-animal models in the lab. Treatments discovered in the lab aim to be translated to clinical use and improve patients' lives. However, the use of animals in biomedical research needs to be carefully evaluated according to Russell and Burch’s 3Rs principles: replacement, reduction, and refinement3. The objective of the reduction component is to reduce the number of animals without compromising data collection. With this in mind, being able to longitudinally evaluate the lesion evolution via noninvasive imaging allows a great advantage in reducing the number of animals required, as well as maximizing the information obtained from each animal4.
Photoacoustic tomography (PAT) is a hybrid imaging modality that combines optical absorption contrast with ultrasound imaging spatial resolution5. The imaging mechanism of PAT is as follows. An excitation laser pulse is illuminated on the target being imaged. Assuming the target absorbs light at the wavelength of the excitation laser, it will increase in temperature. This quick increase in temperature results in a thermoelastic expansion of the target. The expansion causes an ultrasound wave to propagate out from the target. By detecting the ultrasound wave at many positions, the time required for the wave to propagate from the target to the detectors can be used to create an image through a reconstruction algorithm. The ability of PAT to detect optical absorption in deep tissue regions differentiates PAT from ultrasound imaging, which detects boundaries of differing acoustic impedances of tissues5. In the visible and near-infrared spectra, the primary highly absorbing biomolecules that are abundant in organisms are hemoglobin, lipids, melanin, and water7. Of particular interest in the study of stroke is hemoglobin. Since oxyhemoglobin and deoxyhemoglobin have different optical absorption spectra, PAT can be used with multiple excitation laser wavelengths to determine the relative concentration of the two states of the protein. This allows the oxygen saturation of hemoglobin (sO2), or blood oxygenation, to be quantified in and outside of the infarct region8,9. This is an important measure in ischemic stroke, as it can indicate the level of oxygen in the damaged brain tissue following ischemia.
Acoustic angiography (AA) is a contrast-enhanced ultrasound imaging method that is particularly useful for imaging the morphology of vasculature in vivo10. The method relies on the use of a dual-element wobbler transducer (a low frequency element and a high frequency element) in conjunction with microbubbles injected into the circulatory system of the imaging subject. The low-frequency element of the transducer is used for transmitting at the resonant frequency of the microbubbles (e.g., 2 MHz), while the high-frequency element is used to receive the super harmonic signals of the microbubbles (e.g., 26 MHz). When excited at a resonant frequency, the microbubbles have a strong nonlinear response, resulting in the production of super harmonic signals that surrounding body tissues do not produce11. By receiving with a high-frequency element, this ensures that only the microbubble signals are detected. Since the microbubbles are confined to the blood vessels, the result is an angiographic image of blood vessel morphology. AA is a powerful method for imaging ischemic stroke, as the microbubbles that flow through the circulatory system are not be able to flow through blocked vessels. This allows AA to detect regions of the brain that are not perfused due to ischemic stroke, which indicates the infarct region.
Preclinical ischemic stroke research generally relies on the use of histology and behavioral testing to assess the location and severity of the stroke. Triphenyltetrazolium chloride (TTC) staining is a common histological analysis used to determine the stroke infarct volume. However, it can only be used at an endpoint, since it requires the animal to be euthanized12. Behavioral tests can be used to determine motor function impairment at multiple time points, but they cannot provide quantitative anatomical or physiological values13. Biomedical imaging provides a more quantitative approach to studying the effects of ischemic stroke noninvasively and longitudinally9,14,15. However, existing imaging technologies (such as small-animal magnetic resonance imaging [MRI]) can come at a high cost, be unable to provide concurrent structural and functional information, or have limited penetration depth (as most optical imaging techniques).
Here, we combine photoacoustic, ultrasound, and angiographic tomography (PAUSAT; see system diagram in Figure 1), which allows for complementary structural and functional information of blood perfusion and oxygenation following ischemic stroke16. These are two important aspects in assessing the severity of injury and monitoring the recovery or response to treatments. Using these integrated imaging methods can increase the amount of information obtained by each animal, reducing the number of animals required and providing more information in the study of potential treatments for ischemic stroke.
Figure 1: PAUSAT diagram. (A) Complete schematic of the PAUSAT system, including the laser and OPO used for PAT. (B) Inside view of the PAUSAT system, including two ultrasound transducers. The dual-element wobbler transducer is used for both B-mode ultrasound and AA, and the linear-array transducer is used for PAT. Both transducers are mounted on the same 2D motorized stage, allowing for scanning to generate volumetric data. This figure has been modified from16. Please click here to view a larger version of this figure.
All animal procedures were approved by the Duke University Medical Center Animal Care and Use Committee and were conducted in accordance with the United States Public Health Service's Policy on Humane Care and Use of Laboratory Animals. Male and female C57BL/6J mice (see Table of Materials) were used for these studies. A minimum of three animals were imaged per stroke model group. See Figure 2 for the workflow followed in this protocol.
Figure 2: Summary of the experimental procedure for PAUSAT imaging applied to stroke. Created with Biorender.com. The figure shows the workflow of the imaging procedure starting from (A) the two main stroke models (pMCAO and PT stroke). (B) A retro-orbital injection of the microbubbles must be performed prior to positioning the animal on the PAUSAT membrane. (C) A mask providing continuous anesthesia and a heating pad to keep the animal's body temperature stable are required in this setup. The animal's body is placed on the heating pad while the head rests on the membrane of the system. (D) The order of image acquisition is also displayed in the figure. (E) TTC staining is performed to validate our results in this study. DPI: days post-injury. Please click here to view a larger version of this figure.
1. Inducing the stroke mouse model
2. Preparing PAUSAT for imaging
3. Preparing the animal for PAUSAT
NOTE: PAUSAT is performed 1 day following PT stroke surgery or 3 days following pMCAO surgery. Preparing PAUSAT for imaging (step 2) takes about 20 min and should be done immediately before preparing the animal for PAUSAT.
4. PAUSAT Imaging
NOTE: This is done to image the contra- and ipsi-lateral regions of the brain after stroke
5. Optional: Results validation
Imaging of blood vessel morphology in the brain
AA generates blood vessel morphology images by exciting microbubbles in the circulatory system at their resonant frequency and receiving the super harmonic response of the microbubbles. By using the customized ramp (Figure 2C) attached to a manually adjustable stage, we can image the mouse brain with AA mode at two different focal depths. When deeper regions are targeted, more superficial regions (such as the cerebral cortex) show a poorer resolution and signal strength (Figure 3A), and vice versa (Figure 3B). However, by acquiring two focal depths and combining them, AA images can provide information about a whole coronal section (Figure 3C,D). Additionally, by using the motorized stage to scan along a third dimension, PAUSAT can register a series of coronal images covering the region of interest (ROI) defined by the user. These series of images can be aligned and used to visualize a 3D representation of the whole brain, or the ROI defined by the user. However, additional results demonstrating that other quantitative measures (such as the stroke volume) can be analyzed by using the 3D information from the compilation of images acquired by PAUSAT are also provided here.
Anatomical and functional evaluation before and after ischemic stroke
To illustrate the potential of the PAUSAT system in preclinical studies, we studied two groups of animals subjected to one of the models of ischemic stroke analyzed here: pMCAO or PT stroke. These two stroke models differ in the principles by which the ischemic region is created. Briefly, in the pMCAO model, the middle cerebral artery is electro-cauterized, stopping the blood supply from this artery to the brain. This injury triggers a secondary injury, in which surrounding tissue becomes ischemic, enlarging the area affected by the stroke. We decided to image the pMCAO stroke brain at day three after the stroke, because this is the time at which the maximum area expected to be affected by the stroke is achieved. In the PT stroke, however, the maximum area of tissue affected by the stroke is achieved after the first day, so we decided to image PT stroke at day one after the stroke was performed. Although we chose these time points in our study, PAUSAT can be used for longitudinal monitoring of stroke at any desired time point.
First, a B-mode scan was acquired to ensure the correct position of the animal's head and identify the spatial limits between the skull and the brain (Figure 4A and Figure 5A). The most anterior part of the brain imaged includes the frontonasal suture of the skull as a reference point. An ROI of 20 mm (anterior to posterior) x 17.15 mm (lateral) covering the left (contralateral; CL) and the right (ipsilateral; IL) hemispheres was defined in the following studies to locate the stroke region. AA images provided information regarding the blood vessel structures by targeting the signal from the circulating microbubbles. Our results show that in the uninjured brain (baseline), both hemispheres present a similar distribution of blood vessels (Figure 4B and Figure 5B), as expected in the absence of injury. A similar result is observed for the tissue oxygenation map distribution based on PA images obtained at two different wavelengths (Figure 4C and Figure 5C). We decided to evaluate two different wavelengths based on the local maximum optical absorption for deoxygenated hemoglobin (756 nm) and the wavelength at which the deoxygenated and oxygenated hemoglobin presented equal optical absorption (798 nm)20. By capturing these two states of hemoglobin, we can accurately estimate tissue oxygenation (Figure 4D and Figure 5D).
The day after acquiring the baseline image from the brain, we performed surgery on the animal (pMCAO or PT stroke), as described in previous sections. A new set of images were acquired at a specific time point after the stroke (see Figure 2A) to evaluate the location and size of the stroke. The AA images acquired from mice subjected to pMCAO show a notable reduction in the intensity of the signal in the right lateral of the cortex (Figure 4B). The same region shows a reduction in the tissue oxygenation map from PA images, suggesting an ischemic area (Figure 4D). To validate our results, we decided to harvest the brain and perform TTC staining immediately after acquiring the images from the stroke animal. Our results show that the area identified as stroke is similar when compared with our results by PAUSAT and the well-established TTC staining (Figure 4B-E). We also demonstrate here that PAUSAT can identify stroke areas in the brain using a second model of stroke, PT stroke, evaluated at an earlier time point after stroke induction (1 day). As shown in Figure 5B, we were able to identify a region in the upper part of the cortex with a reduced blood flow supply and a concomitant decrease in oxygen saturation (Figure 5D). The results obtained from TTC staining match the location and size of the stroke previously identified by PAUSAT (Figure 5B–E).
The age of the mouse being imaged can heavily affect the image quality. The skull of a mouse becomes thicker with age, and since the skull has a significantly different acoustic impedance relative to soft tissue, a large percentage of the ultrasonic waves are reflected at the skull boundary. Because these integrated imaging systems are all acoustic-based, this leads to a decrease in imaging depth for older mice. This can be clearly seen in Figure 6, where mice of three different ages were imaged with AA under the same conditions. Similar low-signal results are expected if the imaging procedure is not followed as previously described.
Quantitative analysis of the stroke volume by PAUSAT
As described in the above-mentioned results, PAUSAT captures a series of coronal images targeting an ROI described by the user. We analyzed the volume of stroke by calculating the stroke area within the different coronal images and the distance between images. The stroke volume calculated by PAUSAT shows no statistical difference (p < 0.05) compared to the stroke volume calculated based on a similar approach using TTC staining images (Figure 7).
Figure 3: Effect of focal depth of wobbler transducer on the quality of AA images. Images at multiple focal depths can be acquired and later combined to produce the best whole-brain imaging results. (A) AA coronal section acquired at a deeper focal depth. (B) AA coronal section acquired at a more superficial focal depth. (C) Result of fusing deeper focus (green) and more superficially focused (magenta) images. (D) Result of combining deeper focus and more superficially focused images with a pixel scale value comparable to (A) and (B). This figure has been modified from16. Please click here to view a larger version of this figure.
Figure 4: Representative PAUSAT results of a mouse brain (#1) before and after pMCAO stroke. (A) Baseline and post-stroke B-mode ultrasound coronal images. (B) Baseline and post-stroke AA coronal images. (C) Baseline and post-stroke PA coronal images at 756 nm excitation. (D) Oxygen saturation map based on 756 nm and 798 nm excitation images. (E) TTC stained sections of a mouse brain, showing the same region of stroke. The distance between the lines on the bottom of the image represents 1 mm. This figure has been modified from16. Please click here to view a larger version of this figure.
Figure 5: Representative PAUSAT results of a mouse (#2) before and after PT stroke. (A) Baseline and post-stroke B-mode ultrasound coronal images. (B) Baseline and post-stroke AA coronal images. (C) Baseline and post-stroke PA coronal images at 756 nm excitation. (D) Oxygen saturation map based on 756 nm and 798 nm excitation images. (E) TTC stained sections of a mouse brain, showing the same region of stroke. The distance between the lines on the bottom of the image represents 1 mm. This figure has been modified from16. Please click here to view a larger version of this figure.
Figure 6: Effect of mouse age on the quality of AA images. This figure shows the difference in AA signal intensity and imaging depth when animals from different ages are imaged under comparable conditions. As seen in the figure, older mice (18 months) imaging produces lower quality images due to the size and thickness of the skull compared to younger mice (1.5 months), while adult mice show an intermediate signal (6 months). Due to the large acoustic impedance mismatch between the skull and brain tissue, ultrasound waves propagating through the skull are reflected and refracted, leading to signal loss and resolution degradation. Please click here to view a larger version of this figure.
Figure 7: Analysis of the volume of the ischemic stroke. (A) Example of TTC sections from PT stroke and corresponding TTC segmented images, AA images, and AA segmented images at comparable coronal locations. (B) The stroke volume calculated based on TTC and AA images for pMCAO and PT stroke, showing no significant differences (p > 0.05). The graph shows the mean ± standard deviation (S.D.) (n = 3 animals per group). This figure has been modified from16. Please click here to view a larger version of this figure.
There are a few vital aspects of this method that, if done incorrectly, can lead to significantly decreased image quality and quantitative analysis. The most commonly occurring result of user-error in PAUSAT images is either a lack of signal or very low signal strength, both of which can occur for a variety of reasons. One such reason is a problem with the acoustic coupling. Large air bubbles in the water surrounding the mouse's head during imaging can often block the ultrasound from travelling to or from the transducer, causing a shadow region in the resulting image for all three modalities of the system. This can be prevented by ensuring that enough water is present between the system membrane and the sample to be imaged. It could also occur that there is a lack of AA signal, but a surprisingly high PA signal, independent of the wavelengths used. This could be due to the presence of hair-especially dark hair-interfering with the absorption of the light. To avoid this, it is necessary to have previously shaved the animal's head until no hair is visually identifiable.
Regarding AA images, another problem that might occur and would lead to the lack of or very low signal is a low concentration of microbubbles present in the circulatory system. If the microbubbles are too diluted or they were injected improperly, the resulting signal will be very weak. The retro-orbital injection must be performed only by well-trained personnel. Another reason could be a long period of time between microbubble injection and the start of imaging, which can lead to a reduction of the microbubbles in the bloodstream. To avoid this, it is recommended to prepare the system for imaging before injecting the microbubbles, so the animal can be transferred to the PAUSAT system membrane immediately after injection. It is worth mentioning that alternative intravenous routes (such as tail vein injection) can also be used once the animal is already positioned on the PAUSAT system membrane, shortening the time between injection and imaging. Additionally, in the interest of maintaining the maximum concentration of microbubbles possible, it is also recommended to perform the AA imaging before the PA imaging.
The PA signal can also be affected if the necessary steps are not properly performed. One of the issues involves the quality of the excitation laser. The two main components that describe the laser source are the pulse energy and the wavelength. It is highly recommended to use an independent power meter and spectrometer to measure these quantities and ensure they are set to the desired values. If an incorrect value is assumed, the functional calculations will provide misleading results.
It is also possible that when AA and PA images are combined, they are out of alignment. The main reason for this is that the coordinates established between the transducers were not accurate. To prevent this, it is critical to perform the required experiments using a phantom grid beforehand to determine the exact coordinates for a successful co-registration of AA and PA images. Another potential source of inaccuracy is due to incorrect calibration of the OPO. To mitigate this, the OPO must be correctly calibrated using an independent spectrometer.
There remain significant areas of improvement, particularly regarding the image quality of the PAT component of this integrated system. The current PAT system is based on a scanning linear-array configuration. We are able to observe large-scale differences in the oxygenation of infarct versus healthy regions of the brain using several stroke mouse models, however, the detailed vasculature (microvessels) is not able to be seen. There are two main problems that contribute to this low image quality. The first is the limited-view problem. In order to create a full-view reconstructed image of a sample in PAT, the detector must completely surround the object (or have a solid angle of 4π). However, in the experimental setting, this is difficult. This gives rise to the limited view problem, leading to blood vessels that are orthogonal to the transducer array, being undetectable21. Solutions to reduce the effects of the limited-view problem in linear-array PAT do exist, the most promising of which is the use of microbubbles as virtual point sources22. The second problem of linear-array PAT is poor elevational resolution and sensitivity, due to the weak focus of the acoustic lens. However, the use of single-slit diffraction in the hardware design of a linear-array system has been shown to create isotropic resolution and sensitivity for linear-array PAT23,24. Deep learning approaches have also been shown to partially resolve both the limited-view problem and the poor elevational resolution of linear-array PAT25,26,27. The combination of these solutions would significantly enhance the image quality of the PAT component of our integrated imaging system.
Here, we have presented a novel noninvasive multi-modal imaging method for the structural and functional quantification of ischemic stroke in the preclinical setting. Through AA, clear morphology and perfusion mapping of blood vessels in the mouse brain is possible. A locally defined lack of signal can indicate an infarct region where there is decreased blood perfusion, allowing for the estimation of infarct volume noninvasively and longitudinally. Through PAT, the oxygen saturation of hemoglobin can be measured in- and outside of the stroke region, showing the hypoxic state of brain tissues in the infarct region. This method of whole-brain imaging is low-cost with respect to some alternatives, such as small-animal MRI. Furthermore, it allows for the combination of functional and structural information in deep tissue that other preclinical imaging devices are unable to achieve. Compared to histological analyses, this method allows for measurement of these metrics longitudinally and at multiple time points, as opposed to at a single endpoint. The ability to do so will provide unprecedented quantitative detail in the preclinical study of ischemic stroke. Additionally, the structural and functional imaging capability of this system can be applied to many preclinical studies beyond ischemic stroke28. Any disease state or biological phenomenon affecting the vascular system could be studied in depth using PAUSAT, making it a powerful preclinical imaging tool translatable to many other preclinical fields (e.g., cancer).
The authors have nothing to disclose.
The authors would like to acknowledge the engineering team at SonoVol Inc. for their technical support. This work was partially sponsored by the American Heart Association Collaborative Sciences Award (18CSA34080277), to J. Yao and W. Yang; The United States National Institutes of Health (NIH) grants R21EB027981, R21 EB027304, RF1 NS115581 (BRAIN Initiative), R01 NS111039, R01 EB028143; The United States National Science Foundation (NSF) CAREER award 2144788; the Chan Zuckerberg Initiative Grant (2020-226178), to J. Yao; and NIH grants R21NS127163 and R01NS099590 to W. Yang.
20 GA catheter | BD Insyte Autoguard Winged | 381534 | For mouse intubation |
2,3,5-Triphenyltetrazolium chloride | Sigma | T8877 | Necessary for TTC-staining brain for validation |
532nm Laser | Quantel | Q-smart 850 | Laser used to pump the OPO for PAT |
Automatic Ventilator Rovent Jr. | Kent Scientific | RV-JR | To keep mice under anesthesia during surgical procedure |
Black braided silk 4-0 USP | Surgical Specialties | SP116 | Used for sutures on the neck for pMCAO surgery |
Bupivacaine | Hospira | 0409-1159-18 | Used prior to closing wounds during surgical procedure |
C57BL/6 Mice | Jackson Lab | #000664 | Mice used for studying ischemic stroke (2-6 month old male/female) |
Clear suture | Ethicon | 8606 | Used for closing wound (PT stroke and pMCAO). A clear suture won't interfere with PAT |
Cold Light LED | Schott | KL 1600 | Needed to create PT stroke |
Disposable Razor Blade | Accutec Blades | 74-0002 | For sectioning mouse brain |
Electric drill | JSDA | JD-700 | Used to expose MCA during pMCAO procedure |
Electrocauterization tool | Wet-Field | Wet-Field Bipolar-RG | Stops blood flow after drilling during pMCAO procedure |
Hair removal gel | Veet | 8282651 | Used to remove hair from mouse prior to imaging |
High Temperature Cautery Loop Tip | BOVIE Medical Corporation | REF AA03 | Used to avoid bleeding when separating the temporal muscle from the skull |
IR Detector Card | Thorlabs | VRC5 | Used to ensure light path is aligned |
Laser Power Meter | Ophir | StarBright, P/N 7Z01580 | Can be used to calibrate the laser energy prior to imaging |
Laser Speckle Imaging System | RWD Life Science Co. | RFLSI-III | Can be used to validate stroke surgery success |
Lubricant Eye Ointment | Soothe | AB31336 | Can be used to avoid drying of the eyes |
Manually adjustable stage | Thorlabs | L490 | Used with custom ramp for multiple focal depth AA imaging |
Modified Vega Imaging System | Perkin Elmer | LLA00061 | System containing both B-mode/AA and PAT transducers |
Optical Parametric Oscillator | Quantel | versaScan-L532 | Allows for tuning of excitation wavelength in a large range |
Programmable Ultrasound System | Verasonics | Vantage 256 | Used for PAT part of system |
Rose Bengal | Sigma | 330000 | Necessary to induce PT stroke |
Suture | LOOK | SP116 | Used for permanent ligation of CCA |
Temperature Contoller | Physitemp | TCAT-2 | Used to maintain stable body temperature of mice during procedures |
VesselVue Microbubbles | Perkin Elmer | P-4007001 | Used for acoustic angiography (2.43 × 10^9 microbubbles/mL) |