Identification of cholesterol crystals in plaques of atherosclerotic mice using hyperspectral CARS imaging.

lipid-rich regions and characterize general plaque morphology. Noninvasive Abstract The accumulation of lipids, including cholesterol, in the arterial wall plays a key role in the pathogenesis of atherosclerosis. Although several advances have been made in the detection and imaging of these lipid structures in plaque lesions, their morphology and composition have yet to be fully elucidated, particularly in different animal models of disease. To address this issue, we analyzed lipid morphology and composition in the atherosclerotic plaques of two animal models of disease, the low density lipoprotein receptor-deﬁ cient (LDLR (cid:2) / (cid:2) ) mouse and the ApoE lipopro-tein-deﬁ cient (ApoE (cid:2) / (cid:2) ) mouse, utilizing hyperspectral coherent anti-Stokes Raman scattering (CARS) microscopy in combination with principal component analysis (PCA). Hyperspectral CARS imaging revealed lipid-rich macrophage cells and condensed needle-shaped and plate-shaped lipid crystal structures in both mice. Spectral analysis with PCA and comparison to spectra of pure cholesterol and cholesteryl ester derivatives further revealed these lipid structures to be pure cholesterol crystals, which were predominantly observed in the ApoE (cid:2) / (cid:2) mouse model. These results illustrate the ability of hyperspectral CARS imaging in combination with multivariate analysis to characterize atherosclerotic lipid morphology and composition with chemical speciﬁ city, and consequently, provide new insight into the formation of cholesterol crystal structures in atherosclerotic plaque lesions. —Lim, R. S., • LDL receptor • ApoE • hyperspectral coherent anti-Stokes Raman scattering • nonlinear optical microscopy • principal component analysis

spectroscopy, is typically implemented only at selected pixels of interest. As compositional variation within the tissue cannot be assumed to simply be a function of morphology (i.e., spectra acquired from selected pixels do not necessarily refl ect the composition of the overall structure), spatially resolved molecular maps are required to better connect lipid composition and distribution to plaque progression.
In this work, we aim not only to identify atherosclerotic lipid structures in plaque lesions but also to address their composition. Specifi cally, we aim to track the distribution of molecular compounds that have been shown to play major roles in atherogenesis, particularly cholesterol and its derivatives. To selectively detect these compounds in plaque lesions, we have advanced the CARS imaging method such that spectrally resolved molecular maps of the lesion can be taken. Supplementing coherent Raman images with spectral information enables the visualization of several lipid classes selectively ( 31 ). Our method is based on hyperspectral CARS imaging, which was recently successfully used to unveil the compositional lipid variations in gland tissues ( 32 ).
Here we use hyperspectral CARS in combination with principal component analysis (PCA) to extract spectrally signifi cant variations of molecular distributions throughout the lesion. PCA is a statistical analysis tool that enables multidimensional datasets (i.e., hundreds of thousands of pixels taken at multiple wavelengths) to be reexpressed in terms of fewer variables ( 33 ). In this study, we demonstrate that, with PCA, major constituents can be discriminated based on the correlated spectral variations within the dataset itself, thereby providing an unsupervised and effi cient assessment of chemical composition. On this note, the application of PCA for dimensionality reduction, in lieu of other nonlinear techniques, has been recently highlighted ( 34 ), and its utility for spectral decomposition and image reconstruction of hyperspectral CARS datasets has been illustrated in several studies ( 35,36 ), including our own work ( 32 ).
The combined hyperspectral CARS imaging/PCA method has allowed us to examine lipid structures in plaque lesions with greatly improved chemical detail. Importantly, we have found that this technique can successfully discriminate lipophilic compounds, in both intra-and extracellular regions, in a label-free manner. We have applied this method to study the structure and chemical makeup of plaques in two atherosclerotic mouse models: the LDL receptor-defi cient (LDLR Ϫ / Ϫ ) mouse and the ApoE-defi cient (ApoE Ϫ / Ϫ ) mouse (37)(38)(39)(40)(41). Our label-free imaging results reveal clusters of macrophage cells as well as needle-shaped and plate-shaped lipid crystal structures within plaque lesions. Our CARS spectral analysis further establishes that a large fraction of these condensed crystal microstructures are, in fact, ChCs. These results provide important information about plaque morphology and composition, and consequently, carry strong implications on the structural development and progression of plaques during atherogenesis. methods, such as computed tomography (CT) and magnetic resonance imaging (MRI), provide general information about plaque structure, such as the degree of stenosis and the presence of lipid-rich tissue, respectively (13)(14)(15). These methods, although informative of general plaque morphology and features of the vessel wall, typically suffer from either poor contrast resolution (CT) or poor spatial resolution (MRI), and therefore preclude analysis at the cellular level. Invasive techniques, however, such as intravascular ultrasound (IVUS) and optical coherence tomography (OCT), have shown potential for yielding structural and chemical information at a highly detailed level (16)(17)(18). These methods are currently under investigation in large-scale clinical trials. Although these modalities provide insight into plaque structure, they typically require the use of exogenous dyes to provide contrast in the tissue. Furthermore, none of these existing imaging techniques permits chemically selective detection of ChCs. These limitations necessitate the need for labelfree, high-resolution imaging techniques that can simultaneously provide rich, spatially resolved chemical information.
Nonlinear optical microscopy provides a label-free means for imaging key components in biological tissue at submicron resolution and with highly specifi c chemical identifi cation (19)(20)(21). One nonlinear technique -coherent anti-Stokes Raman scattering (CARS) -has emerged as an ideal method for selective imaging of lipids (21)(22)(23). By probing the intrinsic vibrational signatures of carbonhydrogen (CH) bonds, CARS microscopy can easily detect CH-rich lipid structures in tissue. Consequently, CARS, in addition to other nonlinear microscopy techniques, has been applied to several studies on atherosclerosis, including the visualization of extracellular lipids and lipid-rich cells associated with plaque lesions ( 24,25 ) and the determination of lipid concentration levels in various lesions of the arterial intima ( 26,27 ). Previous work of our own utilized a multimodal CARS-based imaging technique to determine the impact of diet on macrophage infi ltration and lipid accumulation in early-stage atherosclerotic plaques ( 28 ).
Although these studies demonstrate the ability of CARS to detect lipids in atherosclerotic tissue, determining the exact composition of these lipid structures remains a diffi cult, yet important, endeavor. In fact, studies have shown that the composition of the plaque, rather than its size, determines plaque stability ( 29 ). Recently, Kim et al. used a multiplex CARS imaging system to identify four general morphologies of atherosclerotic lipids and characterize their chemical profi les ( 30 ). They were able to acquire CARS spectra in the 2,650-3,050 cm Ϫ 1 range from regions with distinct lipid morphology, and it was concluded that changes in the spectra were due to modifi cations of the lipids during the course of disease progression. Nonetheless, the exact makeups of the lipid structures and pools, which include cholesterol and cholesteryl esters, have yet to be fully determined. Furthermore, spectral acquisition using multiplex CARS, although undoubtedly faster than conventional Raman

Principal component analysis
PCA is a chemometric method that reduces the dimensionality of the data, based on the variance in the data itself, into a new coordinate system, also known as principal components ( 33 ). In this study, we utilized PCA to decompose the hyperspectral data in terms of a manageable set of principal components that best describes the spectral variations in the original data. Detailed descriptions on the implementation of PCA for image analysis can be found in several reviews ( 46,47 ), and its applications for hyperspectral CARS imaging has been recently demonstrated ( 32,36 ).
Briefl y, the hyperspectral data (x,y,) can be organized in a single matrix Y, where each column corresponds to intensities measured at different vibrational frequencies. Hence, for a 512by-512 image, Y typically consists of 262,144 rows and 64 columns, corresponding to the total number of pixels and total number of spectral points, respectively. In some cases, several datasets are combined to statistically improve the classifi cation of major components in the image. MATLAB's Singular Value Decomposition (SVD) function was used to compute eigenvectors (V) and eigenvalues (D). Subsequently, the principal components matrix (Z) is calculated based on Z = Y V , exclusively incorporating eigenvectors that contain the most signifi cant variance. Prior to decomposition, Y is centered by subtracting the mean intensity for each column from each pixel.
In addition to dimensionality reduction, the spectra from each pixel can also be generated by using selected principal components ( p ), thereby minimizing the unwanted contribution from noise in the data. For this purpose, the eigenvectors in matrix V are arranged in a decreasing fashion based on the magnitude of its corresponding eigenvalues. The spectra (S) are then computed as S = Z p V p T , where Z p and V p are the truncated principal components and eigenvectors, respectively. For most datasets, the fi rst three principal components account for more than 80% of the variances. In our previous work, we showed that the

Atherosclerotic mouse models and sample handling
The 5/6 nephrectomized ApoE-defi cient (ApoE Ϫ / Ϫ ; stock #002052, n = 6) and LDL receptor-defi cient (LDLR Ϫ / Ϫ ; stock #002207, n = 6) mice were purchased from Jackson Laboratory (Bar Harbor, ME). Details of 5/6 nephrectomy were previously described ( 42 ). Mice were fi rst anesthetized with ketamine/Xylazine. The left kidney was then decapsulated and one-third of each kidney pole ligated with a silk suture. The kidney tissue from each pole was excised with scissors. The muscle and skin layers of the incision were then closed with sterile surgical staples. The animals recovered for one week before undergoing a second surgery, in which the entire right kidney was decapsulated and surgically removed by ligating with a silk suture at the bottom of the renal artery and then excising the tissue with scissors. At 8 weeks of age, animals were fed a Western diet containing 0.15% cholesterol and 21% milk fat (TD88137, Harlan Teklad, Indianapolis, IN)) for 16 weeks. Animals were euthanized and then perfused with 4% buffered paraformaldehyde. Serum cholesterol and triglyceride levels were determined by enzymatic methods (Thermo Fisher Scientifi c, Waltham, MA).
For microscopic examination, complete aortas were removed, opened longitudinally from the heart to the iliac arteries, and mounted on a glass coverslip, luminal side facedown for en face imaging on the inverted microscope. Phosphate buffered saline (PBS) solution was added to prevent the tissue from drying out, and a second coverslip was secured with nail polish on top of the artery to prevent motional artifacts. Animal experiments were approved by the Institutional Animal Care and Research Advisory Committee of the University of Colorado at Denver.

Cholesterol, cholesteryl esters, and triglyceride
Pure cholesterol, tristearin (a triglyceride), and three cholesteryl ester forms (palmitate, oleate, and linoleate) were purchased from Sigma-Aldrich and used without further purifi cation. Tristearin and the cholesteryl esters were solubilized in chloroform and allowed to recrystallize on glass coverslips for CARS spectroscopic examination. Cholesterol monohydrate crystals were grown by recrystallizing cholesterol in water. Bovine serum albumin (BSA) was obtained from Sigma-Aldrich. A saturated aqueous solution of BSA was placed adjacent to the crystal of the lipophilic compound on the coverslip to provide a generic protein reference spectrum within the same fi eld of view ( 43 ).

CARS imaging and spectral analysis
The pump and Stokes beams required for the CARS process were derived from an optical parametric oscillator (Levante Emerald OPO, Berlin, Germany) pumped by a 76 MHz mode-locked Nd:vanadate laser (Picotrain, High-Q, Hohenems, Austria). The fundamental 1,064 nm, 7 ps pulses of the Nd:vanadate laser source provided the Stokes radiation, while the second harmonic 532 nm light was used to pump the OPO (generating tunable wavelength from 780-830 nm). The collinearly overlapped pump and Stokes beams were passed through a laser scanner (Fluoview 300, Olympus, Center Valley, PA) and focused with a 20×, 0.75 NA objective lens (UplanS Apo, Olympus) onto the sample. Images were collected for vibrational frequencies in the range 2,760-3,060 cm Ϫ 1 in intervals of approximately 4 cm Ϫ 1 for 64 different wavelengths ( Fig. 1 ). This vibrational region has been previously identifi ed to contain the most spectral variation for various cholesterol subclasses ( 44 ). The range of vibrational frequencies was achieved by spectral tuning of the OPO, which has been recently discussed in detail ( 45 ). Average acquisition time for one hyperspectral data stack was 20 min.

RESULTS
Lipid-rich macrophage cells CARS images of lipid-rich structures within the aortic walls of LDLR Ϫ / Ϫ mice are shown in Fig. 2 . These structures, which are defi ned by rounded morphologies and distinct dark nuclei, are identifi ed as clusters of macrophage cells ( 28 ). Images were acquired at various regions of the aorta, including the descending/thoracic aorta ( Fig. 2A ), the ascending aorta ( Fig. 2B ), and the abdominal aorta ( Fig. 2C ), thus indicating a broad distribution of macrophage cells throughout the entire artery.
A scoremap ( Fig. 2D ) and corresponding CARS spectra ( Fig. 2E ) derived from the PCA of a representative macrophage CARS image ( Fig. 2A ) were also generated. Two major contributing CARS spectra are seen. The fi rst type of spectra (green spectrum) reveals a typical lipid profi le characterized by the strong peak of the CH 2 symmetric vibration at 2,845 cm Ϫ 1 . Lipid spectra of similar profi le have previously been associated with pools of neutral lipids, such as long-chain aliphatic triglycerides ( 23 ). This type of spectra is prevalent in all regions of the image that can be associated with macrophages. The second type of spectra (red spectrum) yields major peaks at 2,845 cm Ϫ 1 , 2,865 cm Ϫ 1 , 2,905 cm , and 2,965 cm Ϫ 1 . This spectral signature, although markedly different from the fi rst type of spectra, also appears to associate with some of the macrophage cells in the image, therefore suggesting either compositional variation within the macrophage population or the presence of other, spectrally unique lipid structures within the macrophage cells.
reconstructed spectra based on the selected principal components, S, are identical to the original CARS spectra, thus confi rming that the linear transformation and subsequent spectral reconstruction in PCA do not distort the intrinsic features contained in the raw dataset. Because of the undistorted reproduction of the original CARS spectra, the spectral comparison can be performed directly on the CARS dataset without the need for converting the data to linear Raman spectra with phase retrieval algorithms ( 32 ). To enhance the contrast among different constituents in the image, the RGB scoremap is generated based on new axes, Z p , where the spectra from the same colored pixels are averaged and normalized to the peak value. To further facilitate visualization, structures of interest are colored green, red, yellow, violet, and orange, with other compounds colored blue.

Quantitative analysis
Lipid-rich crystal structures, specifi cally needle-shaped and plate-shaped crystals, were quantitatively analyzed in both ApoE Ϫ / Ϫ and LDLR Ϫ / Ϫ mice. Specifi cally, the gross accumulation of cholesterol-associated needle and plate structures was calculated as a function of total tissue area (i.e., number of crystal structures per 125,000 m 2 ). To calculate the number of cholesterol crystal structures, only the fraction of each individual image with spectra similar to pure cholesterol monohydrate (obtained from the PCA-derived scoremaps and then compared with the pure lipid components) was considered in our analysis. Images were then cross-referenced with their original CARS images to match spectral features with morphologic features, thus allowing for an accurate assessment of only cholesterol-associated crystal structures. CARS images from all wavelengths within our range of vibrational frequencies (2,760-3,060 cm Ϫ 1 ) were analyzed to ensure an accurate spatial assessment of all crystals. All values were compared using the t test. All P values less than 0.05 were considered statistically signifi cant. by Kim et al. ( 30 ) These needles and plates appear to form at various locations along the artery, with some located in close proximity to other lipid-rich structures ( Fig. 4A, D ) and others in more isolated regions ( Fig. 4E ). Fig. 4C and F show electronically zoomed-in (8×) views of the needles and plates, respectively. These zoomed-in images, specifically Fig. 4F , reveal even further the complex nature of these lipid-based structures (e.g., smaller plates within the larger plate). In Fig. 5 , a CARS image of representative needles ( Fig.  5A ), its PCA-derived scoremap ( Fig. 5B ), and its corresponding CARS spectra ( Fig. 5C ) are shown. In Fig. 5B , the needle structures are characterized by several distinct spectra (green and red). This suggests that there exists compositional variation even within structures of apparent morphologic similarity. Neighboring lipid-rich cell clusters vary in compositional makeup as well, as some share the same composition as the needles (the green and red regions), whereas others maintain a slightly different constituency (yellow regions). The corresponding CARS spectra in Fig. 5C support these fi ndings. The green spectrum assigned to some of the needle structures and cell clusters exhibits a strong major peak at the lipid CH 2 symmetric stretching vibration ( ‫ف‬ 2,845 cm Ϫ 1 ). However, the red and yellow spectra assigned to other needles and cell clusters exhibit additional strong peaks at 2,845 cm Ϫ 1 , 2,865 cm Ϫ 1 , 2,905 cm Ϫ 1 , 2,945 cm Ϫ 1 , and 2,965 cm Ϫ 1 , in agreement with the CARS point-spectra previously taken at selected crystals ( 30 ). This second spectral signature is also similar to the secondary type of spectra observed in some of the CARS images of lipid-rich structures in the ApoE Ϫ / Ϫ mice are shown in Fig. 3 . Similar to the LDLR Ϫ / Ϫ mice, dense macrophage cell clusters are seen in lesions throughout the entire length of the aorta ( Fig. 3A-C ). A corresponding scoremap ( Fig. 3D ) and CARS spectra derived from PCA ( Fig. 3E ) were also generated. A signifi cant fraction of the macrophage cells exhibit intracellular lipids that spectrally resemble the lipids contained in the macrophage cells of the LDLR Ϫ / Ϫ mouse, as characterized by the lipid-like spectrum in Fig. 3E (green spectrum). We also fi nd additional spectral signatures among the intracellular lipids in the ApoE Ϫ / Ϫ model. As shown in Fig. 3D , some cells exhibit a distinct spectrum that features a prominent peak at 2,930 cm Ϫ 1 (orange spectrum). This spectrum, which was also evident in some regions of the LDLR Ϫ / Ϫ mouse (data not shown), corresponds to a composition of lipophilic compounds that is different from the typical long-chain aliphatic lipid spectrum. Some cells in the image show both types of spectra in distinct intracellular domains. It can be seen that, within the sensitivity of the experiment, the CARS spectrum in the macrophage cells is heterogeneous.

Needle-shaped and plate-shaped lipid crystals
In addition to macrophage cells, other well-defi ned lipid structures are evident in both LDLR Ϫ / Ϫ and ApoE Ϫ / Ϫ mouse models. In particular, as shown in Fig. 4 , we observe thin, needle-like structures ("needles") and fl at, rectangular plate-shaped structures ("plates") in the ApoE Ϫ / Ϫlesions. Similar structures were previously observed

Cholesterol, cholesteryl esters, and triglyceride
To facilitate a direct comparison between the spectral profi les of the LDLR Ϫ / Ϫ and ApoE Ϫ / Ϫ lipid crystal structures and the spectral profi les of known plaque components, hyperspectral CARS imaging of pure cholesterol monohydrate, triglyceride (tristearin), and three forms of cholesteryl esters (cholesteryl palmitate, oleate, and linoleate) was carried out. These cholesteryl esters were chosen because they form main contributions to the LDL lipid content ( 48 ). BSA was analyzed as a control with a known protein spectrum. The results are shown in Fig. 7 . The CARS spectrum of pure cholesterol monohydrate exhibits strong characteristic peaks at 2,905 cm Ϫ 1 and 2,945 cm Ϫ 1 . Cholesteryl linoleate yields a similar spectral profi le but without the 2,945 cm Ϫ 1 (instead, slightly shifted to 2,930 cm Ϫ 1 ) and 2,965 cm Ϫ 1 peaks. Cholesteryl oleate displays a much more pronounced contribution from the macrophage cells of the LDLR Ϫ / Ϫ mouse ( Fig. 2E , red spectrum). This further supports the theory that other spectrally unique lipid structures (in this case, needles) may reside within the macrophage cells. Fig. 6 provides a hyperspectral CARS analysis of the observed plate-like structures. The red and yellow regions associated with the plate structure in Fig. 6B yield spectra that are similar to the red and yellow regions of the needles, but they are different from the predominantly green spectrum of the surrounding lipid-rich regions (only one dominant peak at ‫ف‬ 2,845 cm Ϫ 1 ). Similar to the needle structures, plates with similar morphology also contain spectral variation. These spectral variations implicate the existence of compositional variations within the different lipid-based structures. This fi nding contrasts a prior study that correlated specifi c morphologic features with composition ( 30 ).

Fig. 4. CARS images of various lipid-rich structures in ApoE
Ϫ / Ϫ mice. Thin, needle-shaped structures, some indicated by circles in (A-B), and fl at, rectangular plate-shaped structures, some indicated by boxes in D-E, were observed. Electronically zoomed-in (8×) images were also obtained (C, F) to highlight the needle and plate structures, respectively. . Tristearin and cholesteryl palmitate both yield strong CARS peaks at ‫ف‬ 2,845 cm Ϫ 1 and ‫ف‬ 2,885 cm Ϫ 1 , whereas the protein reference (BSA) shows a distinct, typical protein profi le with a maximum at ‫ف‬ 2,930 cm Ϫ 1 .

Quantitative analysis of lipid-rich crystal structures
From the hyperspectral CARS analysis of the needleshaped and plate-shaped crystal structures illustrated in Figs. 5 and 6 , we observe two predominant spectral signatures (red and yellow spectra) that directly correspond to pure cholesterol monohydrate, as evident by the hyperspectral analysis of pure cholesterol and cholesteryl esters ( Fig. 7 , further elaborated in the "Discussion" section). To quantitatively determine the amount of cholesterol-rich crystals in both ApoE Ϫ / Ϫ and LDLR Ϫ / Ϫ mice, only the fraction of each individual image (from the PCA-derived scoremaps in Figs. 5B and 6B ) with spectra similar to pure cholesterol monohydrate was considered in our analysis. These images were then cross-referenced with their original CARS images to match spectral features with morphologic features, thus allowing for an accurate assessment of only cholesterol-associated crystal structures. The total number of crystals was calculated as a function of total tissue area (i.e., number of crystal structures per 125,000 m 2 ), thereby revealing the abundance of cholesterol-associated crystal structures within each transgenic model.

DISCUSSION
We have demonstrated that morphologic and compositional features of lipid-rich atherosclerotic plaques can be determined in two murine models of disease: the LDLR Ϫ / Ϫ mouse and the ApoE Ϫ / Ϫ mouse. Utilizing hyperspectral CARS imaging in combination with multivariate analysis, we specifi cally observed large, broad clusters of lipid-rich macrophage cells, as well as thin, needle-shaped lipid structures and fl at, rectangular-shaped lipid plates, within aortic plaques. We further determined the composition of  LDLR Ϫ / Ϫ and ApoE Ϫ / Ϫ models. This lipid type is also found in needle-and plate-like microsolids. The second class is characterized by a similar CARS spectrum but with a more prominent peak at 2,930 cm Ϫ 1 , which is also found within the macrophage population. The compositional origin of this spectrum is unclear, but it is likely a mixture of lipophilic compounds and protein-rich structures. The third class constitutes condensed needle-and plate-like lipid structures characterized by two strong peaks at 2,905 and 2,945 cm Ϫ 1 . These needles and plates exhibit rectangular shapes and are thus crystalline in nature. For these crystalline materials, a direct comparison with the in situ spectra in Fig. 7 is possible. Within experimental uncertainty, the spectrum associated with these crystal structures is identical to the cholesterol monohydrate spectrum. Therefore, we attribute this spectral class to the presence of ChCs in the tissue. Finally, the last class is formed by the protein-like spectrum of the extracellular matrix, which, as expected, is retrieved in both mouse models. Given the important role of ChCs in atherogenesis, the ability to detect ChCs with chemical selectivity in unstained arterial tissue is signifi cant. Previous optical microscopy studies have utilized confocal refl ection microscopy to identify ChCs in macrophage cells in tissue cultures and in prepared arterial tissue slices ( 8,9 ). Whereas refl ection microscopy is sensitive to the enhanced refl ectivity of ChCs, it is not chemically selective and thus cannot discriminate between ChCs and other condensed materials in the plaque. In our study, for instance, we have found that the microscopic solid materials found in plaques are not exclusively cholesterol. We observed microstructures that, whereas similar in morphology, vary in composition (mainly of pure cholesterol monohydrate and aliphatic lipids that were devoid of cholesterol). Rather than identifying crystalline material through their refl ective properties, the hyperspectral CARS method assigns a chemical signature to each microstructure to avoid an erroneous identifi cation of crystalline materials. This unique signature allows for the detailed detection of lipid composition at the most basic molecular level. Fig. 2 illustrates this point. The green and red pixels associated with the macrophages have distinct spectra, with the green spectrum exhibiting a typical neutral lipid profi le and the red spectrum identifying pure cholesterol monohydrate. This discrepancy indicates the presence of ChC material within the macrophage cells, an observation not made in previous studies ( 30 ). Hence, the generation of a chemical map with spectral information contained in every single pixel allows for the visualization of cellular processes that were previously inaccessible. Furthermore, hyperspectral CARS imaging enables the identifi cation of ChCs in nonsegmented, label-free arterial tissues with preserved three-dimensional morphology. In principle, this method can be extended to in vivo applications.

Detection of cholesterol crystals in the mouse models
Apolipoproteins, which bind to cholesterol and other fatty acids in the plasma, are responsible for transporting and assisting the clearance of lipoproteins, such as LDL, from the plasma via the LDLR pathway in the liver. ApoE

CARS spectral analysis of lipophilic components
In this study, we have used CARS hyperspectral imaging to gain insight into the distribution and composition of intra-and extracellular lipid-rich pools in plaque lesions. We have limited ourselves to the spectral region of the CH stretching vibrations. This region is known to produce very strong CARS signals, while maintaining a high degree of spectral differentiation among major classes of lipophilic compounds. The CARS signal is related to the Raman spectrum of the molecular compounds in focus, although the features in CARS spectra are generally broader and the amplitude of the higher energy band vibrations is generally lower due to interferences with the nonresonant background. In principle, the CARS spectra can be converted to Raman spectra with the aid of phase retrieval algorithms. To avoid conversion artifacts due to the spectral range examined, we have chosen to analyze the vibrational information contained in the CARS spectra by direct comparison with the CARS spectra in situ of the major lipophilic compounds found in plaque lesions. Such a comparison is warranted in the case of pure states of the compound, as is the case for crystalline materials. Furthermore, PCA is performed to facilitate spectral discrimination, and its unsupervised nature allows for effi cient and direct analysis without knowledge of the sample composition. We note that the principal components are assembled strictly based on linear combinations of mutually correlated variables in the dataset itself. It has been previously shown that spectral reconstructions using the principal components do not distort the intrinsic spectral features contained in the raw dataset ( 32 ).
The CARS spectral analysis reveals four major spectral classes in the plaques. The CARS spectrum of the fi rst class resembles the spectral signatures of long-chain fatty acids and triglycerides in the fl uid or crystalline fl uid phase, with a strong major peak at 2,845 cm Ϫ 1 . This lipid-rich pool is predominantly found in macrophages that have infi ltrated the lesions, which are evident in both the and LDLR are not the only ligand-receptor pair; however, they constitute the majority of receptor-mediated uptake functions. Hence, deleting ApoE or LDLR in mice, in tandem with feeding a high-cholesterol diet, causes severe hypercholesterolemia and atherosclerosis (37)(38)(39)(40)(41). In our study, both LDLR Ϫ / Ϫ and ApoE Ϫ / Ϫ mice were fed the same high-cholesterol diet, yet large clusters of ChC were mainly observed in the ApoE Ϫ / Ϫ mouse. The distribution of ChCs in the ApoE Ϫ / Ϫ mouse extends from the region spanning the ascending aorta to the abdominal aorta, although it is predominant in the ascending aorta (data not shown). Whereas our observation of ChCs in lesions of ApoE Ϫ / Ϫ mice on a high-cholesterol diet is in agreement with previous studies ( 8,9 ), the difference in ChC abundance and distribution in similar regions of the LDLR Ϫ / Ϫ aorta has not yet been reported, although studies have indicated ChC formation in the region of the aortic sinus (tip of the aorta) ( 49,50 ) and the skin ( 51 ). Interestingly, in our samples, most of the ChCs observed in the aorta of the LDLR Ϫ / Ϫ model mainly exist as smaller needles, perhaps refl ecting the less severe atherosclerotic state of the LDLR Ϫ / Ϫ model. This hypothesis, however, requires further investigation. The lowered abundance of ChCs in the aortas of our LDLR Ϫ / Ϫ mice suggests a possible correlation between the absence of the ApoE lipoprotein and the formation of lipid crystals. The presence of ChC in ApoE Ϫ / Ϫ , LDLR Ϫ / Ϫ double-knockout mice supports this hypothesis ( 52 ). Furthermore, an active role for ApoE in ChC formation is supported by the observation of a signifi cant decrease in ChC upon addition of ApoE to the culture media of J774 macrophage cells ( 3 ). Other studies have indicated an alternate pathway for LDL clearance via the LDLR-related receptor (LRP) ( 53,54 ) and have highlighted the importance of ApoE in the formation of HDL, a marker that is inversely correlated with the progression of disease ( 55,56 ). It has been shown that the addition of ApoE to LDLR Ϫ / Ϫ , ApoB100/100 mice does not necessarily correspond to lower cholesterol levels, but it does reduce atherosclerosis ( 57 ). Hence, it is possible that the ApoE lipoprotein may play a role in the complementary biochemical processes responsible for the inhibition of cholesterol crystals. In this regard, hyperspectral CARS imaging studies in double-knockout (ApoE Ϫ / Ϫ , LDLR

CONCLUSION
In this study, we have demonstrated the ability to study the microscopic distribution of various lipophilic compounds within plaque lesions of atherosclerotic tissue. We implemented hyperspectral CARS imaging in combination with multivariate analysis to discriminate various lipid-rich structures, including macrophage cells and cholesterol