In vivo photoacoustic lipid imaging in mice using the second near-infrared window

: Photoacoustic imaging has emerged as a promising technique to improve preclinical and clinical imaging by providing users with label-free optical contrast of tissue. Here, we present a proof-of-concept study for noninvasive in vivo murine lipid imaging using 1210 nm light to investigate differences in periaortic fat among mice of different gender, genotypes, and maturation. Acquired lipid signals suggest that adult male apoE − / − mice have greater periaortic fat accumulation compared to adolescent males, apoE − / − females, and wild-type mice. These results demonstrate the potential of photoacoustic tomography for studying vascular pathophysiology and improving the diagnosis of lipid-based diseases.


Introduction
Lipid-based diseases account for a significant portion of Western world mortality [1][2][3]. Clinical manifestations of lipid-based conditions appear in a multitude of debilitating disease states that include diabetes, obesity, fatty liver disease, demyelination, and vascular diseases such as coronary artery disease and peripheral artery disease. As a result, many imaging techniques including ultrasonography, magnetic resonance imaging, positron emission tomography, x-ray angiography, and optical imaging have been developed to better characterize these conditions. While these imaging modalities can be effective in studying a wide variety of disease states, each has its own inherent limitations such as high cost, limited penetration depth, and/or long imaging times [4]. Therefore, there is still a need to develop techniques that can better characterize compositional information in a rapid, label-free, noninvasive manner that complements currently used imaging methods.
Photoacoustic tomography (PAT) is an emerging modality that can improve in vivo disease characterization in humans and animals [4,5]. Compared to traditional photoacoustic techniques, PAT, using the second near-infrared window, provides optical-dependent contrast that allows users to differentiate biological components such as blood and fat [6,7]. This contrast is produced by using pulsed laser light to induce a transient temperature rise in tissue that causes acoustic wave propagation. These acoustic waves are acquired to reconstruct an image that displays the relative density and location of biological components. PAT also allows the user to tune the wavelength of light to monitor the excitation of distinct molecular overtones. For instance, the second overtone vibrations of the CH 2 bonds can be detected by using 1210 nm light to identify lipid deposition [4]. When overlaid with ultrasound images, we obtain anatomical and compositional information that other modalities cannot easily provide. In fact, previous ex vivo and intravascular studies have shown that PAT can identify mild fatty streaks and discriminate between fibrous and lipid-rich plaques [5]. Despite the exciting potential of PAT, its use in noninvasive in vivo imaging of lipid-based murine disease models to better understand disease initiation and progression has not been explored.
Here we present our findings on how PAT can be utilized to study lipid accumulation in hyperlipidemic mice. We focused on the periaortic fat in mice, as it is ample, metabolically active, and accumulates throughout the vasculature. Studying excessive periaortic fat may also have clinical implication as it is involved in cardiovascular and metabolic dysfunction [8]. We hypothesized that PAT can be utilized to image and detect significant differences in lipid accumulation between apolipoprotein E-deficient (apoE −/− ) and wild-type (WT) mice. The methods and results in this paper are significant as this is the first report to our knowledge that utilizes long wavelength PAT for noninvasive in vivo lipid detection in mice.

Materials and methods
ApoE −/− (n = 17) and WT (n = 10) mice were obtained from Jackson Laboratory (Bar Harbor, ME) and fed a standard chow diet. These mice were separated into six distinct groups based on gender and maturation as shown in Table 1 [9]. A small animal anesthesia system (SomnoSuite, Kent Scientific [10]) was used to keep the mouse under anesthesia using 1.5-3% isoflurane and 225 mL/min room air. Prior to imaging, eye lubricant was applied to each animal to prevent corneal desiccation. We kept rectal temperature at roughly 34-36°C with a heated stage and monitored this using a rectal probe. Furthermore, heart rate was maintained between 450 and 600 beats per minute while respiration rate was kept between 40 and 80 breaths per minute and monitored using electrodes built into the heated stage. We used a high-frequency small animal ultrasound system (Vevo2100, FUJIFILM Visual Sonics) equipped with a 40 MHz center frequency transducer (MS550D) to obtain in vivo long-axis B-mode and PAT images of the infrarenal aorta. A Nd:YAG pulsed optical parametric oscillator (OPO) laser (NT352C, Ekspla) capable of producing 670 to 2300 nm light was used to generate 5 ns pulses at 10 Hz. Our system design has been previously summarized [7] and has been recapitulated in Fig. 1. A function generator (4013B, B&K Precision) was used to synchronize pulse excitation and image acquisition between the laser and the ultrasound, and a pulse generator (9200, Quantum Composers) was used to coordinate the trigger of the laser with the ultrasound system. The function generator specifically produced a 10 Hz, 5 V square waveform. This setup provided us an ultrasound axial spatial resolution of 40 μm with pixel size of 13.2x55.0 μm (axial x lateral), and PAT axial spatial resolution of 124 μm with pixel size of 28.4x26.6 μm [7]. Fig. 1. Schematic depicting PAT setup. Nd:YAG OPO laser delivered light through a fiber optic bundle to an anesthetized mouse. A 40 MHz center frequency transducer was used to acquire the resulting PAT signal. A function generator coordinated pulse excitation and image acquisition between the laser and the Vevo2100 ultrasound system. A pulse generator was used to synchronize the trigger of the laser with the ultrasound system. The ultrasound transducer face is 20x5 mm and fiber optic bundle face is 12x2 mm.
In this study, 1210 nm light was used to detect periaortic fat signal along the anterior wall of the infrarenal aorta ( Fig. 2 [6]) of apoE −/− and WT mice (Table 1). We used 1100 nm light to detect blood signal [7] within the lumen of the aorta and 1250 nm light as a negative offresonance control [6]. Photoacoustic images were acquired with a priority of 80%, gain of 30 dB, power of 100%, brightness of 50, and contrast of 50 in the Vevo2100 system. Priority was decreased to 80% to shift the image intensity histogram leftward by 20%. This eliminated many of the background pixels and provided us with better lipid contrast. We chose to remove 20% as it ensured that no part of the histogram was truncated. To minimize minor PAT signal fluctuations and background noise we acquired ten PAT images for each wavelength and computed the pixel-wise arithmetic mean (Fig. 3(A)). We then performed a thresholding technique to discriminate lipid and background signal. A new histogram of each image was obtained, which produced two peaks that represented the background and lipid signal ( Fig. 3(B)). A threshold point, in the trough of these two peaks, was chosen to assign pixels a value of 255 if their value was greater than or equal to that of the threshold and assigned 0 otherwise. This PAT mask was overlaid on a grayscale PAT image to confirm the accuracy of our thresholding technique (Fig. 3(C)). If the mask did not properly align on the grayscale PAT image, the threshold was either increased or decreased by a value of 10 until the mask aligned with the grayscale image. We then plotted the profile of pixel intensities along multiple lines of interest extending anterior to posterior (Fig. 3(D)). The width of the periaortic fat peak was measured to quantify lipid accumulation on the aorta (Fig. 3(E)). This thresholding method was repeated 10 times for each mouse (270 total measurements) to take into account measurement variability. Oil Red O (ORO) lipophilic staining was performed to confirm presence of periaortic fat on an excised mouse aorta. Finally, an ANOVA with a Tukey post-hoc statistical test was used to determine significance between groups (p<0.05). Images are shown such that the animals are orientated in supine position with the head located to the left and tail is located to the right. Greater periaortic fat signal is evident in the apoE −/− mouse compared to the WT mouse. Using these images we can clearly discriminate the infrarenal aorta (red dotted outline), subcutaneous fat (white arrows), periaortic fat (orange arrows), and blood (red arrows).

Results
Using our PAT setup we were able to resolve subcutaneous and periaortic fat buildup in both apoE −/− and WT mice. We clearly observed lipid specific PAT signal that followed both the geometry of the skin due to the subcutaneous fat and infrarenal aorta due to periaortic fat accumulation. The adult male apoE −/− mice also showed more subcutaneous and periaortic fat accumulation compared to the adult male WT mice (Fig. 2). To confirm that our lipid signal is coming from periaortic fat, we performed ORO staining on excised aortae of apoE −/− mice. We observed that the entire vessel stained red, confirming the presence of periaortic fat (Fig.  3(F)). This suggests that our lipid signal is indeed coming from fat buildup on top of the vessel. All vessels that underwent this ORO process stained red, making it challenging to detect differences in periaortic fat between mice, thus highlighting the benefits of PAT. Prior to quantifying periaortic fat accumulation, we validated our analysis technique to ensure accurate quantification of lipid signal. First, we confirmed that our thresholding technique provides reliable measurements by measuring the inner diameter of polyethylene 50 (PE50) tubing. This was accomplished by acquiring 1100 nm PAT images of a murine blood filled PE50 tube. The thickness of the blood signal was measured at three locations along the tube and compared to the known inner diameter (0.58 mm). Using our threshold algorithm we obtained an experimental inner diameter of 0.553 ± 0.012 mm, resulting in a percent error of 4.7%. This variation is possibly due to differences in acoustic speeds through water and PE50 tubing. We also showed that there is no measurement variation due to respiratory or cardiac motion by quantifying lipid thickness in both the raw PAT and the pixel-wise arithmetic mean images. To prevent measurement bias we quantified thickness at the same location throughout the aorta (Fig. 4). Our analysis showed respiratory and cardiac motion does not affect our results, as there are no statistically significant differences in fat thickness between the raw and mean images (p>0.05). Furthermore, we observed greater lipid signal towards the surface of the mouse due to progressively fewer photons reaching deeper structures. To overcome this limitation we used lipid thickness as our quantification marker rather than signal intensity. We then used our thresholding technique to quantify lipid thickness in adolescent male (0.24 ± 0.06 mm), adult male (0.46 ± 0.12 mm), adolescent female (0.24 ± 0.08 mm), and adult female apoE −/− (0.21 ± 0.03 mm) mice, as well as in adolescent (0.19 ± 0.04 mm) and adult (0.25 ± 0.06 mm) male WT mice (Fig. 5). Adult male apoE −/− mice have significantly greater periaortic fat compared to other apoE −/− and WT groups (p<0.01).

Discussion
In this study we evaluated the hypothesis that PAT can be used to image and quantify periaortic fat accumulation in vivo. By coupling a high frequency ultrasound system and a pulsed laser, we were able to acquire images of subcutaneous and periaortic fat. We are certain that our PAT images provide lipid-specific signal for multiple reasons. We specifically utilized 1210 nm light to excite the second overtone vibrations of CH 2 bonds, which are abundantly found in lipids [5]. The in vivo subcutaneous fat signal also follows the geometry of the skin, while the periaortic fat signal follows the geometry of the aorta. Finally, we confirmed that our images are indeed detecting periaortic fat by performing ORO lipophilic staining on excised murine aortae (Fig. 3(F)). Our results support the application of PAT as a noninvasive in vivo method for preclinical studies of lipid-based conditions.
We then demonstrated the utility of this technique by detecting differences in periaortic lipid accumulation among groups of mice with different genetic background, age, and gender. Our results showed that adult male apoE −/− mice had significantly greater periaortic fat accumulation compared to adolescent males, female apoE −/− , and WT mice (p<0.01). These results are consistent with the literature, which shows that male apoE −/− mice are more prone to develop cardiovascular disease compared to female apoE −/− mice, as they are used for aneurysm, deep vein thrombosis, and atherosclerosis studies [11][12][13][14][15]. In fact, male apoE −/− mice are more prone to developing atherosclerotic lesions especially when fed a high fat diet [11]. We also investigated if female apoE −/− mice fed a high fat diet would accumulate more periaortic fat compared to female mice fed a normal chow diet. In a study with separate animals, we found that there were no differences in periaortic fat thickness between mice that were placed on a high fat diet for 3 months versus normal chow (data not shown), suggesting that diet may affect luminal and periaortic fat buildup differently. While further studies are needed to explore the various mechanisms affecting luminal and periaortic fat deposition, PAT can be a powerful tool for providing noninvasive information about lipid accumulation.
While PAT has the potential to improve preclinical and clinical imaging, there still remains room to optimize this technique to improve depth of penetration and minimize photoacoustic artifacts. One limitation of our current PAT technique is that much of the 1210 nm light is absorbed by subcutaneous fat and reflected off of the skin, therefore decreasing the depth of penetration. Light attenuation by skin, muscle, and blood also limit penetration of photons to the posterior periaortic fat. As a result, we were unable to quantify lipid accumulation on the posterior aortic wall. Moreover, our current PAT setup is susceptible to artifacts that match the geometry of the skin layer and sometimes appear within a region of interest. This artifact typically appears at two times the vertical distance from the probe-skin boundary and we hypothesize that it appears due to light reflecting off of the skin surface and causing a photoacoustic effect at the probe interface. This is clearly seen in the WT images as a thick skin geometry-matching PAT signal at depths of 5 to 8 mm (Fig. 2). Imaging atherosclerotic plaques in rodents with our current PAT system is also challenging and requires further optimization for in vivo characterization of murine atherosclerosis. Therefore, future work will be focused on improving our system by increasing the light fluence through the skin and to employing mechanisms to minimize light reflection towards the PAT probe.

Conclusion
The work described here illustrates the unexplored application of PAT for preclinical in vivo lipid detection. While further optimization is possible, PAT has a clear advantage for rapid, noninvasive, and label-free imaging. Moreover, PAT is easily coupled with ultrasound, providing both anatomical and compositional information of tissue with better depth of penetration compared to conventional optical techniques. Development of this technology has the strong potential to lead to a myriad of applications, such as studying lipid-based disease pathology, evaluating therapeutics, and aiding in the diagnosis of human disease.

Funding
Funding was provided by the American Heart Association through the Scientist Development Grant (14SDG18220010) and a Jim and Diann Robbers Research Grant to CJG.