20 September 2013 Microscope-based near-infrared stereo-imaging system for quantifying the motion of the murine epicardial coronary arteries in vivo
Author Affiliations +
J. of Biomedical Optics, 18(9), 096013 (2013). doi:10.1117/1.JBO.18.9.096013
Atherosclerosis is a leading cause of mortality in industrialized countries. In addition to “traditional” systemic risk factors for atherosclerosis, the geometry and motion of coronary arteries may contribute to individual susceptibility to the development and progression of disease in these vessels. To be able to test this, we have developed a high-speed (∼40 frames per second) microscope-based stereo-imaging system to quantify the motion of epicardial coronary arteries of mice. Using near-infrared nontargeted quantum dots as an imaging contrast agent, we synchronously acquired paired images of a surgically exposed murine heart, from which the three-dimensional geometry of the coronary arteries was reconstructed. The reconstructed geometry was tracked frame by frame through the cardiac cycle to quantify the in vivo motion of the vessel, from which displacements, curvature, and torsion parameters were derived. Illustrative results for a C57BL/6J mouse are presented.
Long, Zhu, and Friedman: Microscope-based near-infrared stereo-imaging system for quantifying the motion of the murine epicardial coronary arteries in vivo



Coronary artery atherosclerosis is a leading cause of morbidity and mortality in Western societies. Atherosclerosis is a progressive fibroinflammatory disease identified by intimal thickening, the focal accumulation of lipids, fibrous elements, and cellular elements within the walls of large arteries, known as atherosclerotic lesions.1 Lesions preferentially develop at arterial branches, the outer walls of bifurcations, and the inner wall of curved sections; the cause of this focal vasculopathy is not fully understood. It is, however, understood from epidemiological and clinical studies that individual susceptibility to the development and localization of atherosclerotic lesions is influenced by systemic risk factors, including smoking, diabetes mellitus, obesity, hypertension, and high cholesterol. However, these risk factors cannot account for half of the variability in occurrence of this disease;2,3 this indicates additional risk factors have not been identified. The epicardial coronary arteries, which are attached to the surface of the heart, undergo substantial motion during a cardiac cycle. Thus, one prevalent explanation of the focal nature of the disease is that the local fluid mechanical stresses at the walls of coronary arteries,4 as well as mechanical stresses within the vessel wall, may mediate the phenotype of endothelial cells, thereby producing atherosusceptible sites.5,6,7,8,9 Therefore, it has been speculated that certain aspects of arterial geometry and motion, which vary substantially among individuals, may increase an individual’s susceptibility to developing atherosclerosis—geometric risk factors.9,10,11

To study the pathogenesis and evolution of disease, including atherosclerosis and the geometric risk factor, longitudinal studies are effective. However, in addition to limited patient availability and associated risks and costs, the slow progression and long asymptomatic phase of atherosclerosis in humans make longitudinal studies of early human atherosclerosis impractical. As a consequence, animals, in particular mice, have become important models for experimental atherosclerosis research.12,13 Advantages of using the mouse are its short gestation period (21days), lower relative husbandry cost, and, in particular, the ability to knockout (or disrupt) endogenous genes associated with lipid metabolism and regulation.14,15 The two most commonly used murine models of atherosclerosis are the low-density lipoprotein receptor deficient mouse (LDLR/)16 and the apolipoproteinE-deficient mouse (apoE/). Even when fed normal chow, mice deficient in apoE develop hypercholesterolemia, with total plasma cholesterol levels nearly five times that of nonmutant mice, and spontaneously develop lesions that are phenotypically consistent, in composition and location, with those of human atherosclerosis.17,18,19 Hence, we can advantageously use genetically modified mice to explore the geometric risk factor concept, if we can image the in vivo dynamics of the coronary arteries. Not withstanding these benefits and the availability of knockout models, cardiac imaging of mice is very challenging due to the small size of the heart (8-mm major axis) and rapid heart rate (500600beats/min). For those reasons, cardiac imaging of mice requires a comparable increase in both spatial and temporal resolution relative to that used in the clinical setting.

Small-animal noninvasive imaging system resolution has increased with the use of respiratory and cardiac gating. Even so, current noninvasive imaging techniques such as microcomputed tomography (micro-CT),20,21 magnetic resonance imaging (MRI),22,23,24 and coronary microangiography25 still lack sufficient resolution to dynamically image the motion of murine epicardial coronary arteries over a cardiac cycle. For example, to image cardiopulmonary structures in mice, Badea and colleagues20 developed a micro-CT imaging system that uses cardiac and respiratory gating, and achieved an isotropic spatial resolution of 100 μm. However, since the diameter of the major murine coronary artery is on the order of 100 μm, this is insufficient spatial resolution and these images cannot be acquired real-time during a single heartbeat.26 MRI has also been applied to small-animal imaging. For example, using cardiac-gated MRI to achieve a spatial resolution of 100 μm, Ref. 24 acquired images of the left coronary artery of a living mouse. However, these images were acquired only during mid-diastole, when cardiac motion is at a minimum. Reference 25 used in vivo microangiography to acquire digital images of murine coronary vessels less than 500μm in diameter at a pixel resolution of 7 μm. Although the spatial resolution of microangiography is significantly higher than that of gated micro-CT or MRI, it lacks sufficient temporal resolution since the maximum frame rate was 30 frames per second (fps). To acquire murine arterial motion data over a cardiac cycle, either an ex vivo heart preparation must be used, or the chest of the mouse must be opened to expose the beating heart; coronary arteries on a beating heart in vivo are more natural than those on an isolated heart.

Beyond the small heart size and rapid heart rate, optically imaging the murine coronary arteries has an additional challenge. As with other rodents, the murine coronary artery becomes intramyocardial just distal of its ostium,27 so it is not easily distinguished from the myocardium, even under direct microscopic observation.28 In order to image the dynamic motion of the artery, a fluorescence vascular contrast agent is then needed with four specific requirements: (1) minimal to no extravasation after introduction so that multiple successive cardiac cycles can be continuously imaged; (2) exceptionally bright and stable fluorescence signal so that paired images can be acquired at a sufficient frame rate (40fps); (3) red-shifted (>650nm) for increased tissue penetration to distinguish artery from myocardium; and (4) a small required intravenous injectate volume (<0.2ml). We selected nontargeted near-infrared quantum dots (Qdots) that have an emission maximum wavelength of 705nm. Reference 29 reviewed the fabrication, properties, and use of Qdots for live-cell and in vivo imaging. Qdots are roughly protein-sized (10 to 20 nm) bright fluorophores with a high quantum yield and are highly resistant to photobleaching. Qdots have been used in numerous in vivo imaging studies from whole animal,30,31 to angiography of murine squamous cell carcinoma,32 to microcirculation angiography,33 and as intravascular probes.34 Qdots are an ideal vascular contrast agent for imaging the in vivo dynamics of the murine coronary arteries.

To optically image the in vivo dynamics of murine epicardial coronary arterial segments over a single cardiac cycle, we have developed an imaging system. The system combines high-speed (40fps) near-infrared stereo-imaging into a standard microscope. The microscope is used as both the surgical microscope and the optics for the imaging system. Using this system, paired coronary artery images are synchronously acquired; the two images are calibrated and, under the constraints of epipolar geometry, the three-dimensional (3-D) vessel geometry is reconstructed. The reconstructed axis is tracked frame-by-frame to quantify the in vivo motion of the vessel. Illustrative results show that this system is a promising tool for laboratory investigations of vascular mechanics and experimental atherosclerosis.


Materials and Methods


Animals and Surgery

All animal experiments were conducted under a protocol approved by the Institutional Animal Care and Use Committee of Duke University and in accordance with the Guide for the Care and Use of Laboratory Animals (NIH Publication No. 85–23). Experiments were performed using a healthy male wild-type mouse (C57BL/6J, The Jackson Laboratory, Bar Harbor, ME), 15 weeks of age, weighing 29.2 g. The mouse was presedated with an intraperitoneal injection of ketamine (100mg/kg, Fort Dodge Animal Health, Fort Dodge, IA) and xylazine (5mg/kg, Lloyd Laboratories, Shenandoah, IA) and then placed in the supine position on an operating board and 38°C heating pad (Gaymar Industries Inc., Orchard Park, NY). The animal’s body temperature was monitored throughout the entire procedure with a rectal temperature probe (RET-3, Physitemp Instruments, Clifton, NJ). An endotracheal tube was inserted perorally and a small-animal ventilator (1.5kPa ventilator pressure, 100breaths/min respiratory rate, SAR-830A/P, CWE Inc., Ardmore, PA) and an Isoflurane (1 to 2 vol.%) medical-oxygen mixture was used to maintain anesthesia during the procedure. A PE-10 polyethylene catheter (0.28 mm inside diameter, 0.61 mm outside diameter) was placed in the right external jugular vein for contrast agent infusion. To expose the left epicardial coronary artery for imaging, an oblique skin incision (10mm) was made just below (2 to 3 mm) the animal’s left axilla toward its xiphoid cartilage. The skin and pectoralis muscle were dissected away by blunt dissection so that the heart and left lung could be visualized through the intercostal muscles. A small incision (<10mm) was made in the intercostal space between the fourth and fifth ribs counting from the xiphoid cartilage. Next, the ribs were gently retracted using a small-animal retraction system (Fine Science Tools, Foster City, CA). Using fine forceps, the pericardium was gently picked up and opened. The pericardium was attached to the animal’s chest wall to provide a “natural cradle” for the heart. All mice were euthanized via exsanguination under anesthesia.


Imaging System

The imaging system shown in Fig. 1(a), which also served as the surgical microscope, consisted of a stereo-microscope (SMZ-1000) with a 0.5× objective (P-Plan Apo, numerical aperture of 0.012 at 0.8× zoom magnification and 0.05 at 8× zoom magnification), fluorescent filter attachment, and dual beam splitter (Nikon Instruments, Melville, NY), an excitation light source, a computer, a calibration pyramid, and two high-performance synchronized digital CCD cameras with 2 GB of internal memory (pco.1600, PCO-TECH, Romulus, MI). The angle of view difference between the two cameras was 6.3 deg. During image acquisition, the arterial images were stored on the camera’s internal memory and then, after each experiment, downloaded to a computer for image processing. The near-infrared emission light of the contrast agent was filtered using 18-mm diameter 690- to 730-nm bandpass emission filters (710AF40, Omega Optical, Brattleboro, VT). To excite the fluorescent contrast agent, a 150-W fiber-optic illuminator (Schott-Fostec, LLC, Auburn, NY) was used; it was outfitted with an EKE high-output halogen bulb (USHIO, Cypress, CA) and depleted of infrared wavelengths using a hot mirror (0 angle of incidence, Edmund Optics, Barrington, NJ). Attached to the illuminator was a liquid light guide (8 mm core diameter, 1 m long, Newport, Irvine, CA) and a fiber-optic focusing lens unit (Edmund Optics), which was modified to hold the 25-mm-diameter Qdot excitation filter (shortpass 650 nm cutoff, Edmund Optics) [see Fig. 1(b)]; the focusing unit was attached to an articulating arm mounting system to position the focusing unit during image acquisition. The calibration pyramid [see Fig. 1(c)] used in this study to define the imaging geometry was analogous to that used in previous biplane angiographic studies of human coronary arterial motion.10,35 The pyramid was a black Delrin pyramid (11.5mm×11.5mm base, 8.0 mm height) embedded with 88 0.5-mm-diameter stainless steel spheres (Salem Specialty Ball, Canton, CT), with the centroid of each sphere’s relative location known.

Fig. 1

(a) The microscope-based imaging system developed to quantify the motion of the murine epicardial coronary arteries consists of the following: (a) stereomicroscope, (b) dual beam splitter, (c) fluorescent filter attachment, (d) two high-performance synchronized CCD cameras, (e) computer monitor to observe the image acquisition, (f) an excitation light source, and a computer (not shown) to control the cameras. (B) The transmission spectra versus wavelength for the excitation filter (dashed line) and emission filter (solid line). (C) Photomicrographs of the top view (left image) and side view (right image) of the calibration pyramid used to calibrate the acquired images for image reconstruction. Note: The horizontal distance between marks corresponds to 1 mm.



Imaging Contrast Agent

In order to enhance the contrast between the left epicardial coronary artery and the epicardium for imaging, each mouse was intravenously administered nontargeted near-infrared Qdots (Qtracker 705, 2 μM, Invitrogen/Molecular Probes Inc., Eugene, OR, maximum emission wavelength of 705 nm) via the jugular vein catheter just before image acquisition. The Qdot solution was administered at a volume flow rate of 5μl/s using a syringe pump (Kent Scientific Corporation, Torrington, CT) for a final injected concentration of 8.5pmol-Qdot/g animal weight, which for a 30-g mouse corresponds to a 100 μl injectate volume. Image was acquired immediately after injecting the contrast agent, and using the pco.camware software (PCO-TECH), a sequence of paired fluorescent arterial images (400pixels×400pixels) was synchronously acquired at a frame rate of 40fps for two to three successive cardiac cycles, which corresponded to approximately eight frames per cardiac cycle. The paired images were acquired during a constant-pressure breath hold (pressure plateau). To achieve the breath-hold, the ventilator was switched to constant inhalation, and an external air pressure regulator (Type 70, 0 to 2 psi, Bellofram Corporation, Newell, WV) attached to the ventilator created a pressure plateau of 15cm H2O, which was preset before experiments. After imaging the beating heart, the ventilator was reset to normal operation, and without disturbing the microscope, the animal was removed and replaced with the calibration pyramid in approximately the location of the animal’s heart during image acquisition. Paired images of the calibration pyramid were acquired.


Image Processing

The image processing method used to reconstruct and track the 3-D murine coronary artery axis using a sequence of paired images was based on the technique developed for tracking human coronary arteries in biplane angiograms.35,36 It involves three main steps: (1) calibrating the imaging system; (2) reconstructing the 3-D vessel axis in the first frame; and (3) tracking the 3-D vessel axis through the entire sequence of images.


Imaging System Calibration

At the time of image acquisition, the imaging geometry of the system was not known; therefore, a calibration process was needed to reconstruct a 3-D vessel axis from two two-dimensional (2-D) projection images of the vessel. The calibration process used in this study was described in greater detail elsewhere35,36 and will only be summarized here. The 2-D position of each 0.5-mm-diameter stainless steel sphere in the calibration pyramid image was located by detecting circles using Hough transform techniques. A projective transformation of a point from 3-D to a 2-D plane can be represented by


where (u,v) is the 2-D coordinate of the point in the projection plane (image), and (x,y,z) is the coordinate of the point in the 3-D space. The coefficients c11,c12,,c34 can be solved when six or more spheres are located in the image. For the stereo-imaging system, there were two projection transformations, T1 and T2, which represent, respectively, the first and second views obtained from each camera. Once the coefficients were determined, the projection transformations T1 and T2 were used to reconstruct the 3-D vessel axis.


3-D Vessel Axis Reconstruction

The axis of a coronary artery segment was modeled parametrically by a 3-D space curve, p(s)=[x(s),y(s),z(s)]T, where x, y, and z were the 3-D coordinate functions, and s[0,1] were the parametric space along the length of the vessel axis. For a vessel point in one image, its corresponding vessel point in the other image could be located under the constraints of epipolar geometry, and the 3-D coordinates of this vessel point could be reconstructed. An initial plan for the vessel axis was obtained in the form of a bead chain by simply selecting several points along the axis and connecting them with straight lines. Assuming that the 3-D vessel axis projected into the intensity ridge in the two images, an energy function that was a combination of the intensities of the two images was minimized to refine the 3-D initial plan. The 3-D refined plan was then fit by a cubic spline curve and resampled at intervals of 0.2mm. This was the reconstructed axis of the 3-D vessel segment in the first frame of the image sequence.


Vessel Tracking in 3-D

After the 3-D vessel axis was located in the first pair of frames, it was automatically tracked for the remaining sequence of images. The frame-by-frame motion of the vessel was decomposed into two parts: (1) a global motion modeled by 3-D affine transformation and (2) a local motion of each point along the vessel axis, which makes the whole tracking model deformable. Similar to locating the initial plan in the first frame, the vessel was tracked by regarding its position in the previous frame as an initial guess. In addition to attracting the curve to the brightest position in the images, a string of small neighborhood windows along the vessel axis was set to be matched to the next frame, so that each point on the vessel could find its corresponding position at the next time point as accurately as possible.


Quantification of Motion Parameters

Once a vessel axis segment had been reconstructed and tracked from a sequence of paired images, it could be used to characterize coronary arterial motion, including torsion, τ, and curvature, κ. Curvature is a measure of how rapidly the curve changes direction at a point. Torsion is a measure of the rate at which the oscillating plane is turning along the vessel axis. These quantities could be determined using a cubic polynomial fitting method previously described.36 For a given point, P, on the reconstructed 3-D vessel axis, a 2-mm window centered at P was selected. A half window typically contained five points. Within a window, the x, y, and z coordinates were individually fitted in a least squares sense to cubic polynomials:


where η[1,1]. Derivatives of x, y, and z with respect of η are taken and evaluated at η=0 to determine torsion and curvature at P given by




where r(η)=x(η)e^x+y(η)e^y+z(η)e^z.


Validation with a Coronary Arterial Phantom

To assess the capabilities of this stereo-imaging system and the 3-D reconstruction process, we constructed, imaged, and then reconstructed the geometry of a coronary artery phantom whose shape was known. The phantom was constructed to closely represent both the size and 3-D path of the main murine coronary artery.37 The phantom artery was a 5.25-mm space curve engraved into an 8-mm-diameter acrylic cylinder [see Fig. 2(a)]. The path of the phantom was described parametrically by the following equations:


where R is the radius of the cylinder, h=4mm, and the parameter θ ranges from π/8 to π/8. [Note: the x-y-z coordinate system is defined in Fig. 2(a).] The phantom artery was placed under the imaging system, and under epi-illumination, paired images of the phantom were acquired. Without disturbing the microscope, the phantom was replaced with the calibration pyramid and paired images of it were acquired. Using these images, the imaging system was calibrated and the 3-D phantom axis was reconstructed.

Fig. 2

(a) Paired images of the 5.25-mm-long phantom coronary artery (black curve) with the reconstructed three-dimensional (3-D) axis of the phantom (dotted white curve) superimposed on each image and, on the left image, the coordinate system. The phantom geometry is described parametrically by x=Rsinθ, y=0.5hsin(4θ), and z=R(1cosθ), where R is the radius of the cylinder (R=4mm), h=4mm, and θ ranges between 22.5 and 22.5 deg. Comparison between the phantom axis given by the parametric equation and the reconstructed 3-D phantom axis in (b) the x-y plane, (c) the z-y plane, and (d) the z-x plane. Comparison between the actual values (solid curves) of the phantom axis and the values calculated from the reconstructed phantom axis (dotted curves) of (e) torsion and (f) curvature along the phantom axis.





Results of Coronary Arterial Phantom Validation

Paired images were acquired of the 5.25-mm-long phantom coronary artery [Fig. 2 (a), black curve]. The phantom coronary artery axis was reconstructed. Superimposed on each image in Fig. 2(a) is the reconstructed 3-D axis of the phantom [Fig. 2(a), dotted white curve] and, on the left image, the coordinate system. The actual phantom axis, defined by the parametric equations in Eq. (5), was compared against the reconstructed 3-D phantom axis and shown in Fig. 2, (b) the x-y plane, (c) the z-y plane, and (d) the z-x plane. The root-mean-square (RMS) error for the x, y, z coordinate was 0.010, 0.011, and 0.014 mm, respectively. In this figure, the solid line represents the value of interest for the actual phantom axis and the dotted line represents the value of interest for the reconstructed phantom axis. Point-by-point comparison was made between the actual values (solid curves) of the phantom axis and the values calculated from the reconstructed phantom axis (dotted curves) of (e) torsion and (f) curvature along the phantom axis. The RMS error in torsion was 0.098mm1, whereas the RMS error in curvature was 0.74mm1.


Results of In Vivo Coronary Artery Motion


3-D vessel axis reconstruction

The imaging system was used to quantify the motion of a left epicardial coronary artery segment of a 29.2-g male C57BL/6J mouse. The mouse was intravenously injected with 8.5pmol/g of nontargeted, near-infrared Qdots and then imaged at a frame rate of 40 fps. Illustrative paired fluorescence images are shown in Fig. 3(a) for one instant in a cardiac cycle. For that instant, a 4.3-mm-long left coronary artery segment was reconstructed [Fig. 3(b)]. Using the reconstructed vessel axis, the torsion and curvature of the vessel centerline was determined. For this instant in the cardiac cycle, the vessel centerline torsion ranged from a minimum of 0.02mm1 to a maximum of 4.62mm1 with a mean torsion of 0.95mm1, whereas the vessel centerline curvature ranged from a minimum of 0.13mm1 to a maximum of 1.22mm1 with a mean curvature of 0.48mm1. Once the vessel axis was reconstructed for one instant in a cardiac cycle, the motion of the vessel axis was tracked and the frame-by-frame motion determined.

Fig. 3

(a) Paired fluorescence images of the main branch of the left epicardial coronary artery in vivo (denoted by arrow) on the beating heart of a 29.2-g C57BL/6J mouse imaged, at a frame rate of 40fps, with the stereo-imaging system shown in Fig. 1. The mouse was intravenously injected with 8.5pmol/g of nontargeted, near-infrared quantum dots imaging contrast agent. Superimposed on each image in (a) is the 3-D reconstructed vessel axis segment (dotted white line). (b) For the paired images in (a), the 4.3-mm-long reconstructed segment of the vessel axis is calculated from 22 locations spaced 0.2 mm along the vessel axis and displayed in approximately the same orientation of the heart in (a). Note: The distance between axis tick marks is 0.2 mm.



Vessel tracking in 3-D

Over two successive cardiac cycles, paired fluorescence images of the left coronary of a C57BL/6J mouse were acquired at 40fps (eight frames per cardiac cycle). Figures 4(a) to 4(c) show paired fluorescence images at (a) frame 1, (b) frame 4, and (c) frame 7 of a cardiac cycle. As illustrative results, trajectories at four axial locations are shown in Fig. 4(d). Over the two successive cardiac cycles, the total vessel axis arc length ranged from a minimum of 3.6mm to a maximum of 4.3mm.

Fig. 4

(a) to (c) Paired fluorescence image sequence of the left coronary artery of a C57BL/6J mouse acquired at a frame rate of 40 fps with the stereo-imaging system shown in Fig. 1. Superimposed on each image in (a) to (c) is the reconstructed vessel axis (dotted line) for that pair of images. The images are from the (a) first, (b) fourth, and (c) seventh frames during a cardiac cycle for which a total of eight frames were acquired. (d) A reconstructed axis (solid line) of the left coronary artery from the initial plan (i.e., first frame of cardiac cycle) (a), with trajectories, covering the first seven frames during a cardiac cycle [shown in (a) to (c), attached at four selected axial locations. Note: On each axis in (d), the distance between axis tick marks is 0.2 mm.



Quantification of motion parameters

Using the tracked and reconstructed sequence of images shown in Fig. 4, the temporal and spatial variations of motion parameters were determined. Figure 5 is a pseudocolor plot that shows the temporal and spatial variation of (a) torsion and (b) curvature for two successive cardiac cycles, with the first cardiac cycle corresponding to results shown in Fig. 4.

Fig. 5

For two sequential cardiac cycles, a pseudocolor plot of the variation with axial position and time of (a) torsion and (b) curvature of the left epicardial coronary artery segment shown in Fig. 4.


The variations in torsion and curvature were calculated from a reconstructed image sequence of 15 frames that was captured at 40 fps. Vessel axis point 0 corresponds to the proximal end of the reconstructed segment, whereas point 21 corresponds to the distal end of the reconstructed vessel segment. These quantities, as seen in Fig. 5, vary both spatially and temporally during each cardiac cycle. Over the two successive cardiac cycles, the vessel axis torsion ranged from a minimum of 0mm1 to a maximum of 6.59mm1 with a mean torsion of 1.11mm1. The vessel axis curvature ranged from a minimum of 0.05mm1 to a maximum of 1.58mm1 with a mean curvature of 0.57mm1.



In summary, we presented a microscope-based stereo-imaging system and image reconstruction and motion tracking process that is capable of quantifying the in vivo dynamics of the murine coronary arteries. Illustrative imaging and reconstruction results for C57BL/6J mouse were presented. The reconstructed vessel axis geometry was tracked frame by frame through two successive cardiac cycles to quantify the in vivo motion of the vessel segment, and displacements, curvature, and torsion parameters were determined. We have restricted the imaging to the left coronary artery. As reviewed by Ref. 38, it appears that human coronary artery atherosclerosis is more prevalent in the left coronary artery than the right coronary artery. Although the hemodynamics and anatomical differences between these two arteries may play a role, the pathophysiological mechanisms involved remain unknown.

This microscope-based stereo-imaging system complements small-animal noninvasive cardiac imaging systems such as MRI and micro-CT. Key features/benefits of the present system include (1) the use of near-infrared Qdots, which provide increased tissue penetration and are resistant to photobleaching and extravasation, allowing successive cardiac cycle imaging; (2) sufficient temporal resolution to permit measurement of murine epicardial coronary arterial dynamics during a cardiac cycle; (3) the use of stereoscopy for 3-D imaging; and (4) the ability to support serial studies.

All Qdots when continuously excited exhibit fluorescence intermittency (blinking) and, therefore, spontaneously alternate between emitting (on) and nonemitting (off) states. This alternating is most easily observed in single Qdot tracking/imaging;39,40 it can also be exploited to verify the presence of single Qdots due to their blinking, or clusters of Qdots due to their apparent nonblinking.41 Moreover, fluorescence intermittency does not substantially impact measurements/imaging of ensemble or clusters of Qdots, such as the case for vascular contrast agents. Previous studies have reported the use of Qdots as a vascular contrast agent with blinking not adversely affecting their measurements.30,42

Because optical imaging is used, the animal’s chest must be surgically opened to expose the beating heart. In a study of the differences in hemodynamics and force–frequency relationships between the open-chest and closed-chest mouse models, Ref. 43 determined that differences in these factors do, in fact, exist. In these earlier experiments, a bilateral thoracotomy was employed. The surgical technique used in the present work is based on the survival surgical procedure used in coronary artery ligation studies, which is much less invasive than a thoracotomy and may result in a smaller alteration of thoracic mechanics. Because a survival surgical procedure is used, the experimental system described here can be used in serial studies if tissue collection is not needed at intermediate time points.

One of the central goals of atherosclerosis research is to identify risk factors of the disease to promptly identify individuals at risk, or to improve the prognosis of diagnosed individuals. We believe that specific coronary arterial kinematics will be directly correlated to the initiation and progression of atherosclerotic lesions in these vessels. For example, in recent studies of the dynamics of human coronary arteries,36,35 it was found that the motion of the left anterior descending artery is much different than that of the right coronary artery: (1) less displacement and (2) greater axial variation in torsion. The motion parameters were also found to vary greatly between individuals and along an individual vessel axis. Motion parameters that exhibit large axial variations might be related to variations in susceptibility to atherosclerosis. No studies have integrated the analysis of motion parameters, histomorphometry, and genomics. We have developed a microscope-based stereo-imaging system capable of imaging the in vivo motion of murine epicardial coronary arteries. Using the system described in this paper, we can investigate whether certain arterial motions produce a pro-atherogenic environment, by imaging the coronary arteries of apoE-deficient mice with different stages of disease, followed by histological and genomic characterization. In future studies, this stereo-imaging system could be used to investigate the relationship between the motion of epicardial coronary arteries and the degree of atherosclerotic change in those vessels. Researchers have established18 that mice deficient in apoE and on normal chow spontaneously develop most of the pathological aspects (i.e., composition and sites of disease) of human atherosclerosis over a few months. This relatively short time-to-disease, compared to that of humans, and the imaging system described herein will enable a serial study of atherosclerosis that integrates arterial motion parameters, histomorphometry, and genomics using a mouse model that is relevant to the disease in humans.


This work was supported by National Institutes of Health grants HL058856 (M. H. F., H. Z.) and HL083715 (D. S. L.). The authors would like to thank Dr. Lan Mao (Department of Medicine, Duke University) for training in mouse cardiac surgery and Ellen Dixon-Tulloch (Department of Biomedical Engineering, Duke University) for assistance with animal studies.



R. Ross, “Atherosclerosis—an inflammatory disease,” N. Engl. J. Med. 340(2), 115–126 (1999).NEJMAG0028-4793http://dx.doi.org/10.1056/NEJM199901143400207Google Scholar


T. Pearson, “Coronary arteriography in the study of the epidemiology of coronary disease,” Epidemiol. Revs. 6(1) 1984).0193-936XGoogle Scholar


J. R. Crouse, “Progress in coronary artery disease risk-factor research: what remains to be done?,” Clin. Chem. 30(7), 1125–1127 (1984).CLCHAU0009-9147Google Scholar


W. C. Aird, “Endothelial cell heterogeneity and atherosclerosis,” Curr. Atheroscler. Rep. 8(1), 69–75 (2006).CARUCZ1523-3804http://dx.doi.org/10.1007/s11883-006-0067-zGoogle Scholar


C. G. CaroJ. M. Fitz-GeraldR. C. Schroter, “Atheroma and arterial wall shear, observation, correlation and proposal of a shear dependent mass transfer mechanism for atherogenesis,” Proc. R. Soc. Lond. B Biol. Sci. 177(1046), 109–159 (1971).PRLBA40080-4649http://dx.doi.org/10.1098/rspb.1971.0019Google Scholar


W. C. Aird, “Phenotypic heterogeneity of the endothelium: II. Representative vascular beds,” Circ. Res. 100(2), 174–190 (2007).CIRUAL0009-7330http://dx.doi.org/10.1161/01.RES.0000255690.03436.aeGoogle Scholar


S. G. FrangosV. GahtanB. Sumpio, “Localization of atherosclerosis—role of hemodynamics,” Arch. Surg. 134(10), 1142–1149 (1999).ARSUAX0004-0010http://dx.doi.org/10.1001/archsurg.134.10.1142Google Scholar


A. G. Passeriniet al., “Coexisting proinflammatory and antioxidative endothelial transcription profiles in a disturbed flow region of the adult porcine aorta,” Proc. Natl. Acad. Sci. U S A 101(8), 2482–2487 (2004).PNASA60027-8424http://dx.doi.org/10.1073/pnas.0305938101Google Scholar


M. H. Friedmanet al., “Relationship between the geometry and quantitative morphology of the left anterior descending coronary artery,” Atherosclerosis 125(2), 183–192 (1996).ATHSBL0021-9150http://dx.doi.org/10.1016/0021-9150(96)05869-8Google Scholar


Z. H. DingM. H. Friedman, “Dynamics of human coronary arterial motion and its potential role in coronary atherogenesis,” J. Biomech. Eng. 122(5), 488–492 (2000).JBENDY0148-0731http://dx.doi.org/10.1115/1.1289989Google Scholar


L. Goubergritset al., “Geometry of the human common carotid artery. A vessel cast study of 86 specimens,” Pathol. Res. Pract. 198(8), 543–551 (2002).PARPDSGoogle Scholar


A. Daugherty, “Mouse models of atherosclerosis,” Am. J. Med. Sci. 323(1), 3–10 (2002).AJMSA90002-9629http://dx.doi.org/10.1097/00000441-200201000-00002Google Scholar


J. JawienP. NastalekR. Korbut, “Mouse models of experimental atherosclerosis,” J. Physiol. Pharmacol. 55(3), 503–517 (2004).JPHPEI0867-5910Google Scholar


A. S. Plumpet al., “Severe hypercholesterolemia and atherosclerosis in apolipoprotein E-deficient mice created by homologous recombination in ES cells,” Cell 71(2), 343–353 (1992).CELLB50092-8674http://dx.doi.org/10.1016/0092-8674(92)90362-GGoogle Scholar


J. A. Piedrahitaet al., “Generation of mice carrying a mutant apolipoprotein E gene inactivated by gene targeting in embryonic stem cells,” Proc. Natl. Acad. Sci. U S A 89(10), 4471–4475 (1992).PNASA60027-8424http://dx.doi.org/10.1073/pnas.89.10.4471Google Scholar


S. Ishibashiet al., “Hypercholesterolemia in low density lipoprotein receptor knockout mice and its reversal by adenovirus-mediated gene delivery,” J. Clin. Invest. 92(2), 883–893 (1993).JCINAO0021-9738http://dx.doi.org/10.1172/JCI116663Google Scholar


Y. Nakashimaet al., “ApoE-deficient mice develop lesions of all phases of atherosclerosis throughout the arterial tree,” Arterioscler. Thromb. 14(1), 133–140 (1994).ARTTE51049-8834http://dx.doi.org/10.1161/01.ATV.14.1.133Google Scholar


J. D. SmithJ. L. Breslow, “The emergence of mouse models of atherosclerosis and their relevance to clinical research,” J. Intern. Med. 242(2), 99–109 (1997).JINMEO1365-2796http://dx.doi.org/10.1046/j.1365-2796.1997.00197.xGoogle Scholar


J. L. JohnsonC. L. Jackson, “Atherosclerotic plaque rupture in the apolipoprotein E knockout mouse,” Atherosclerosis 154(2), 399–406 (2001).ATHSBL0021-9150http://dx.doi.org/10.1016/S0021-9150(00)00515-3Google Scholar


C. BadeaL. W. HedlundG. A. Johnson, “Micro-CT with respiratory and cardiac gating,” Med. Phys. 31(12), 3324–3329 (2004).MPHYA60094-2405http://dx.doi.org/10.1118/1.1812604Google Scholar


C. T. Badeaet al., “4-D micro-CT of the mouse heart,” Mol. Imaging 4(2), 110–116 (2005).MIOMBP1535-3508Google Scholar


L. W. Hedlundet al., “MR-compatible ventilator for small animals: computer-controlled ventilation for proton and noble gas imaging,” Magn. Reson. Imaging 18(6), 753–759 (2000).MRIMDQ0730-725Xhttp://dx.doi.org/10.1016/S0730-725X(00)00154-5Google Scholar


F. Wiesmannet al., “High-resolution MRI with cardiac and respiratory gating allows for accurate in vivo atherosclerotic plaque visualization in the murine aortic arch,” Magn. Reson. Med. 50(1), 69–74 (2003).MRMEEN0740-3194http://dx.doi.org/10.1002/(ISSN)1522-2594Google Scholar


J. Ruffet al., “Magnetic resonance imaging of coronary arteries and heart valves in a living mouse: techniques and preliminary results,” J. Magn. Reson. 146(2), 290–296 (2000).JOMRA40022-2364http://dx.doi.org/10.1006/jmre.2000.2147Google Scholar


T. Yamashitaet al., “In vivo angiographic detection of vascular lesions in apolipoprotein E-knockout mice using a synchrotron radiation microangiography system,” Circ. J. 66(11), 1057–1059 (2002).CJIOBY1346-9843http://dx.doi.org/10.1253/circj.66.1057Google Scholar


A. W. Wetzelet al., “Measurement and modeling of 4D live mouse heart volumes from CT time series,” Proc. SPIE 6491, 64910J (2007).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.704792Google Scholar


B. Fernandezet al., “The coronary arteries of the C57BL/6 mouse strains: implications for comparison with mutant models,” J. Anat. 212(1), 12–18 (2008).JOANAY0021-8782http://dx.doi.org/10.1111/joa.2008.212.issue-1Google Scholar


D. Ahnet al., “Induction of myocardial infarcts of a predictable size and location by branch pattern probability-assisted coronary ligation in C57BL/6 mice,” Am. J. Physiol. Heart Circ. Physiol. 286(3), H1201–H1207 (2003).0363-6135http://dx.doi.org/10.1152/ajpheart.00862.2003Google Scholar


X. Michaletet al., “Quantum dots for live cells, in vivo imaging, and diagnostics,” Science 307(5709), 538–544 (2005).SCIEAS0036-8075http://dx.doi.org/10.1126/science.1104274Google Scholar


B. Ballouet al., “Noninvasive imaging of quantum dots in mice,” Bioconjug. Chem. 15(1), 79–86 (2004).BCCHES1043-1802http://dx.doi.org/10.1021/bc034153yGoogle Scholar


M. K. Soet al., “Self-illuminating quantum dot conjugates for in vivo imaging,” Nat. Biotechnol. 24(3), 339–343 (2006).NABIF91087-0156http://dx.doi.org/10.1038/nbt1188Google Scholar


N. Y. Morganet al., “Real time in vivo non-invasive optical imaging using near-infrared fluorescent quantum dots,” Acad. Radiol. 12(3), 313–323 (2005).1076-6332http://dx.doi.org/10.1016/j.acra.2004.04.023Google Scholar


D. R. Larsonet al., “Water-soluble quantum dots for multiphoton fluorescence imaging in vivo,” Science 300(5624), 1434–1436 (2003).SCIEAS0036-8075http://dx.doi.org/10.1126/science.1083780Google Scholar


M. E. Akermanet al., “Nanocrystal targeting in vivo,” Proc. Natl. Acad. Sci. U S A 99(20), 12,617–12,621 (2002).PNASA60027-8424http://dx.doi.org/10.1073/pnas.152463399Google Scholar


Z. H. DingH. ZhuM. H. Friedman, “Coronary artery dynamics in vivo,” Ann. Biomed. Eng. 30(4), 419–429 (2002).ABMECF0090-6964http://dx.doi.org/10.1114/1.1467925Google Scholar


Z. H. DingM. H. Friedman, “Quantification of 3-D coronary arterial motion using clinical biplane cineangiograms,” Int. J. Card. Imaging 16(5), 331–346 (2000).1569-5794http://dx.doi.org/10.1023/A:1026590417177Google Scholar


D. Kumaret al., “Distinct mouse coronary anatomy and myocardial infarction consequent to ligation,” Coron. Artery Dis. 16(1), 41–44 (2005).CADIEX0954-6928http://dx.doi.org/10.1097/00019501-200502000-00008Google Scholar


Y. S. Chatzizisiset al., “Is left coronary system more susceptible to atherosclerosis than right? A pathophysiological insight,” Int. J. Cardiol. 116(1), 7–13 (2007).IJCDD50167-5273http://dx.doi.org/10.1016/j.ijcard.2006.03.029Google Scholar


I. Chunget al., “Description of the fluorescence intensity time trace of collections of CdSe nanocrystal quantum dots based on single quantum dot fluorescence blinking statistics,” Phys. Rev. E, Stat. Nonlin. Soft Matter Phys. 73(1), 011106 (2006).PLEEE81539-3755http://dx.doi.org/10.1103/PhysRevE.73.011106Google Scholar


S. F. LeeM. A. Osborne, “Photodynamics of a single quantum dot: fluorescence activation, enhancement, intermittency, and decay,” J. Am. Chem. Soc. 129(29), 8936–8937 (2007).JACSAT0002-7863http://dx.doi.org/10.1021/ja071876+Google Scholar


M. Jonaset al., “Detecting single quantum dot motion with nanometer resolution for applications in cell biology,” IEEE Trans. Nanobiosci. 5(4), 246–250 (2006).ITMCEL1536-1241http://dx.doi.org/10.1109/TNB.2006.886559Google Scholar


I. TexierV. Josser, “In vivo imaging of quantum dots,” Methods Mol. Biol. 544, 393–406 (2009).MMBIED0097-0816http://dx.doi.org/10.1007/978-1-59745-483-4Google Scholar


B. D. HoitN. BallR. A. Walsh, “Invasive hemodynamics and force-frequency relationships in open- versus closed-chest mice,” Am. J. Physiol. 273(5 Pt 2), H2528–H2533 (1997).AJPHAP0002-9513Google Scholar

David S. Long, Hui Zhu, Morton H. Friedman, "Microscope-based near-infrared stereo-imaging system for quantifying the motion of the murine epicardial coronary arteries in vivo," Journal of Biomedical Optics 18(9), 096013 (20 September 2013). http://dx.doi.org/10.1117/1.JBO.18.9.096013


3D image processing

Imaging systems

In vivo imaging

3D image reconstruction


Image segmentation

Back to Top