Trans-illuminated laser speckle imaging of collateral artery blood flow in ischemic mouse hindlimb

Abstract. The mouse ischemic hindlimb model is used widely for studying collateral artery growth (i.e., arteriogenesis) in response to increased shear stress. Nonetheless, precise measurements of regional shear stress changes along individual collateral arteries are lacking. Our goal is to develop and verify trans-illumination laser speckle flowmetry (LSF) for this purpose. Studies of defibrinated bovine blood flow through tubes embedded in tissue-mimicking phantoms indicate that trans-illumination LSF better maintains sensitivity with an increasing tissue depth when compared to epi-illumination, with an ∼50% reduction in the exponential decay of the speckle velocity signal. Applying trans-illuminated LSF to the gracilis muscle collateral artery network in vivo yields both improved sensitivity and reduced noise when compared to epi-illumination. Trans-illuminated LSF images reveal regional differences in collateral artery blood velocity after femoral artery ligation and are used to measure an ∼2-fold increase in the shear stress at the entrance regions to the muscle. We believe these represent the first direct measurements of regional shear stress changes in individual mouse collateral arteries. The ability to capture deeper vascular signals using a trans-illumination configuration for LSF may expand the current applications for LSF, which could have bearing on determining how shear stress magnitude and direction regulate arteriogenesis.


Introduction
The importance of adequate remodeling of pre-existing arterial interconnections to form endogenous collateral bypasses-i.e., arteriogenesis-is highlighted in the extensive link between adequate collateral development and enhanced patient outcomes with arterial occlusive disease. 1,2 The key initiating stimulus for arteriogenesis is increased shear stress. 3 Upon occlusion of a major artery, downstream pressure is reduced, causing an increase in the pressure gradient, blood flow, and shear stress along pre-existing collateral arteries that bypass the occlusion. Both the magnitude 4 and duration 5 of increased shear stress determine maximal collateral outgrowth and eventual resolution. Nonetheless, direct measurements of shear stress magnitude along individual collateral arteries have not been reported. We believe this paucity of hemodynamic information is due, at least in part, to the lack of a method capable of making such measurements. The primary challenges are the small caliber of the developing collateral vessels (often <100 μm in an initial diameter 2 ) and the need to integrate these data over a large area. Of the technologies with sufficient spatial resolution, such as multiphoton excitation fluorescence, 6,7 optical coherence tomography, 8,9 and photoacoustic tomography, 10,11 there are still multiple barriers to their adoption for studying arteriogenesis. These include cost and technical expertise, depth of imaging, and data integration over the observed vascular networks.
One potential solution is laser speckle flowmetry (LSF) microscopy. Both relative flow and functional microvascular density have been measured by LSF in cerebral, [12][13][14][15] retinal, 13,16,17 and dorsal skinfold window chamber microcirculation. From a theoretical standpoint, absolute LSF measurements are possible, 13 and we have recently shown LSF to be an efficient solution for measuring flow changes across large microvascular networks. 18 The use of LSF has, however, been limited by the depth of signal acquisition to the most superficial vascular structures. 19 In turn, this limits its application for analysis of collateral arteries embedded deeply in the tissue. To increase acquisition depth, a recent report 20 suggested that the reorientation of the illumination source to transmission through the tissue could project and capture the hemodynamic signals from deeper vascular structures. Here, we developed a trans-illuminationbased LSF system to measure, for the first time, the in vivo spatial distribution of collateral artery hemodynamics before and after femoral artery ligation (FAL) in the mouse ischemic hindlimb, which is the most widely used model of peripheral arteriogenesis.

LSF Setup and Data Analysis
When a coherent light source passes through a static tissue, it creates a random interference pattern known as speckle. The speckle pattern produced is characteristic for a given tissue structure. When particles move within the given tissue (i.e., red blood cells), it causes the speckle pattern to fluctuate, such that when these fluctuations are integrated over time it is possible to use spatial statistics to determine tissue motion. 21 We have previously shown that, when using a simplified algorithm developed to remove the assumption of Gaussian or Lorentzian distributions, 22 it is possible to find relative changes in blood flow in the microcirculation at an individual vessel level. 18 The approach relates the change in speckle contrast [C, Eq. (1)] to autocorrelation time [τ, Eq. (2)], which is inversely proportional to blood velocity [speckle intensity (SI), Eq. (2)] (1) Here, C is derived on a per-pixel basis from the standard deviation [σ, Eq. (1)] and mean pixel intensity [hIi, Eq. (1)] values in a 7 × 7 neighborhood. The 7 × 7 neighborhood size was chosen based on a balance between improved signal-to-noise ratio and decreased spatial resolution with an increasing window size. We have previously demonstrated this to be applicable for mouse arterioles of the diameter of gracilis collateral arteries. 18 For our calculations, we used a simplified algorithm [Eq.
(2)] for relating speckle contrast (C) and exposure time (T) to autocorrelation time [τ, Eq. (2)]. Vessel speckle velocity was normalized speckle index (NSI) to the background tissue according to Eq. (3). As previously described, 18 this process removes background artifact from general tissue motion and variations in speckle illumination to allow for comparison of velocity change across experiments In our arrangement, a 30-mW, 658-nm laser diode (LPM658-30, Newport Corporation, Irvine, CA) was coupled to a fiber optic cable and used to illuminate the tissue or tissue phantom from either an epi-illumination or trans-illumination orientation, where the fiber output was placed at the base of the (phantom) tissue. An illustration of the laser speckle imaging system, as applied to the mouse gracilis adductor collateral arteries in the trans-illumination orientation, is shown in Fig. 1. For both in vitro and in vivo testing, a cooled, monochrome chargecoupled device camera (Optronics Quantifier, Goleta, CA) (7.4 × 7.4 μm 2 pixels) was used to acquire the raw speckle images through an intravital microscope (Zeiss Axioskop, Thornwood, NY) using a 4× air objective (Zeiss Acroplan LD NA ¼ 0.1). As previously described, 18 the objective and camera were chosen to balance the resolution of the scope with the size of the speckle to satisfy the Nyquist sampling criteria of at least 2 pixels per individual speckle. 23 For each field of view, laser position, and flow setting, a sequence of ≥15 12-bit raw speckle images was acquired with 5-ms exposure time to capture average velocity over multiple cardiac cycles. Camera gain and light path intensity were varied to maintain mean pixel intensity of the raw speckle image. To minimize the influence of whole background tissue motion from mouse movement, the processed SI images were normalized to median background intensity. Code for converting raw speckle images into speckle flow maps was developed, implemented, and analyzed in Fiji, 24 and is freely available upon request. Individual flow images were then merged into larger two-dimensional maps using Adobe Photoshop (CS2, Adobe Systems Inc., San Jose, CA).

Tissue Phantom Studies
While initial studies suggested that speckle signal analysis from transmitted laser light yields dominant vascular signals from deeper tissue structures 20 and that speckle signals remain viable through up to 2 cm of biological tissue, 25 the applicability of signal processing used for epi-illumination versus trans-illumination was unknown. Therefore, light scattering tissue phantoms were used to test speckle velocity linearity between transillumination and epi-illumination speckle arrangements. Phantoms were generated as previously described 18 by adding titanium dioxide nanopowder (677469, Sigma-Aldrich, St. Louis, MO) to 184 Sylgard PDMS gel (Dow Corning, Midland, MI) at 10 mg∕mL and pouring the mixture over optical glass-bottomed molds. Phantoms were then allowed to cure overnight at 37°C. Molds 5-mm thick were composed of polyethylene tubing connected to glass capillary tubes with 142-μm inner diameter (Drummond Microcap, P1799 Sigma-Aldrich, St. Louis, MO) placed at 0, 120, 360, and 720 μm from the surface. These values were chosen to approximately represent the range of depths at which the collateral arteries reside under the top surface of the mouse gracilis muscle, which has a thickness of ∼750 μm.
For validating the linear operation of the speckle processing algorithm for both epi-illumination and trans-illumination arrangements, a sequence of known flows of defibrinated bovine blood (Hemostat Labs, Dixon, CA) was driven via a syringe pump (PHD2000, Harvard Apparatus, Holliston, MA) through the phantom. The randomized sequence and image acquisition were controlled by a custom-designed LabVIEW program (National Instruments, Austin, TX). Only data from the most Fig. 1 In vivo intravital imaging of collateral artery hemodynamics. The fiber optic cable was positioned for either trans-illumination or epiillumination. For the trans-illumination configuration, the fiber optic outlet was placed on the opposite side of the leg from the exposed muscle and against the depilated skin to yield an evenly illuminated field of view. For the epi-illumination configuration, the fiber optic outlet was positioned to produce an unobstructed and even sheet of light on the surface of the exposed muscle. The given field of view was shifted using an xy-stage.
superficial capillary tube was analyzed. Each complete program of flow rates was repeated five times (i.e., n ¼ 5) for both epiillumination and trans-illumination configurations.
For testing the effect of increasing depth on speckle signal attenuation, a single flow rate was applied to the phantom, and raw speckle images were captured at each capillary tube depth (0, 120, 360, and 720 μm). The image capture sequence at each depth setting was repeated five times for both epiillumination and trans-illumination configurations.

In Vivo Analysis of Deep Collateral Blood Flow in Trans-illuminated Mouse Hindlimb
All animal protocols were approved by the Institutional Animal Care and Use Committee at the University of Virginia and conformed to all regulations for animal use outlined in the American Heart Association Guidelines for the Use of Animals in Research. Male C57BL/6 mice (n ¼ 10) were purchased from Charles River Laboratory (Wilmington, MA) and housed in the animal facilities at the University of Virginia. To produce consistent flow increases in the collateral arteries in a superficial adductor muscle, the FAL scheme in Fig. 2 was chosen. This particular ligation pattern has been previously shown to produce consistent arteriogenesis in the collateral arteries of the gracilis adductor muscles [26][27][28][29] with minimal heterogeneity in the baseline collateral structure and with the predicted changes in flow direction from baseline. Mice, 9 to 13 weeks of age, were anesthetized (i.p. 120 mg∕kg ketamine, 12 mg∕kg xylazine, and 0.08 mg∕kg atropine), depilated, and prepped for aseptic surgery. On the left leg, an incision was made directly above and along the femoral artery, and a window of skin was dissected free and retracted directly above the superficial adductor muscles. The femoral artery was gently dissected from the femoral vein and nerve between the bifurcation of the superior epigastric artery and popliteal artery. Two 6.0 silk sutures were placed underneath the femoral artery, but not yet tied, immediately distal to the epigastric artery, which served as the origin of the muscular branch artery in all mice.
Upon placement of the ligature (untied), the mouse was transferred to the intravital microscopy station, and the exposed tissue was maintained under a constant drip of warmed Ringer's solution (137.9 mmol∕L NaCl, 4.7 mmol∕L KCl, 1.2 mmol∕L MgSO 4 , 1.9 mmol∕L CaCl 2 , 23 mmol∕L NaHCO 3 ) that was prepared with 5% CO 2 balance N 2 , and set to maintain 37°C fluid temperature. Mouse body temperature was maintained by a surgical heating pad. For imaging of gracilis collateral hemodynamics, we used the intravital microscopy setup shown in Fig. 2. Before baseline imaging, muscles were allowed 15 min to equilibrate while being superfused with warmed Ringer's solution. After the equilibration period, a baseline (preligation) series of image acquisitions using epi-illumination and then trans-illumination laser orientation were taken at multiple fields of view to encompass the gracilis collateral arteries between the saphenous and the muscular branch arteries. Between each series of raw speckle images (epi-and transillumination), images of the collateral arteries were taken using brightfield epi-illumination through an objective-mounted fiber optic light guide to allow for lumenal diameter measurements (A08650, Schott Inc., Elmsford, NY). After baseline acquisition, the ligature was tightened, and stoppage of blood flow through the femoral artery was visually confirmed in all mice. After 30 min post-ligation, trans-illumination speckle and brightfield then epi-illumination and brightfield imaging sequences were repeated. Predicted baseline and post-ligation flow direction were confirmed in an identical preparation on a separate set of mice through fluorescent imaging of intravascularly injected fluorescent microspheres (FP-3070-2, Spherotech Inc., Forest Lake, IL).
Blood velocity analysis of laser speckle images was limited to well-defined collateral artery regions, such that bulk microvascular flow in the surrounding gracilis muscle flow was not considered. The assumption of Poiseuille flow was made, where wall shear rate is proportional to average blood velocity divided by tube diameter. This assumption is valid for most cases in the microcirculation. 30 Therefore, proportional wall shear rate [speckle shear stress (SSR)] was calculated using the NSI measurement of average blood velocity and vessel diameter (D) at each designated location [Eq. (4)] 3 Results

Influence of Illumination Source on Sensitivity as a Function of Depth in Tissue Phantoms
The results from the tissue phantom experiments are shown in Fig. 3. Both trans-and epi-illumination exhibited a strong linear relationship between NSI and absolute velocity [ Fig. 3(a)]. However, the trans-illumination signal was not as significantly diminished with an increasing tissue depth, showing an ∼50% reduction of the exponential decay of the speckle velocity signal with an increasing depth compared to epi-illumination [1.20 × 10 −3 versus 2.70 × 10 −3 , Fig. 3(b) at a constant velocity of 5 mm∕s]. As such, trans-illumination maintained 40% of the initial superficial signal out to 720μm deep within the tissue phantom versus 15% with epiillumination.

Trans-illumination LSF Microscopy Better Demonstrates Regional Variation of Hemodynamic Changes Along Collateral Arteries After FAL
A representative comparison of in vivo images created by epiillumination and trans-illumination LSF is shown in Fig. 4. While epi-illumination LSF [ Fig. 4(a)] was somewhat capable of visualizing the most superficial elements of the gracilis collateral network (e.g., the saphenous artery and upstream region of the muscular branch), it generally yielded much poorer quality images. Epi-illumination LSF images exhibited significantly more noise when compared to the trans-illumination images [ Fig. 4(b)], with much of this noise being generated by superficial fat cells overlying the muscle [ Fig. 4(a), white arrows), especially in the collateral entry region near the saphenous artery. Epi-illumination images were also susceptible to errors created by light being directly reflected off the tissue surface [ Fig. 4(a), green arrows]. As seen in Fig. 4(b), these noise and error issues were entirely obviated by switching to transillumination.
Post-FAL, as expected, both the central [ Fig. 4(b), white boxes] and entrance regions [ Fig. 4(b), yellow arrows] of the gracilis collateral network experienced significant increases in the blood flow that were clearly evident with trans-illumination. Epi-illumination imaging was capable of capturing some of the characteristics of increased flow, but sensitivity within the collateral network was considerably poorer. The FAL also induced an expected drop in saphenous artery flow [ Fig. 4(b), yellow boxes], which was evident with both epi-and trans-illumination. However, the pre-FAL epi-illumination signal was often saturated, and both the pre-and post-FAL saphenous flow signals were prone to considerable noise [ Fig. 4(a), cyan arrows]. Finally, it is important to note that, while trans-illumination yielded better image quality, further improvements in spatial resolution could be realized by incorporating more sophisticated temporal laser speckle contrast analyses. 31,32 Limitations on camera control and frame rate precluded such analyses in the current study.

Velocity and Shear Rate Changes Generated by FAL
Because there were no changes in depth of the collateral artery network within the gracilis muscle from pre-to 30 min post-FAL, NSI could be used to determine relative velocity changes [ Fig. 5(a)]. When coupled with diameter measurements obtained from brightfield imaging, the velocity changes were used to calculate relative change in wall shear rate from pre-to 30 min post-FAL [ Fig. 5(b)]. As expected, there was an increase in the blood flow and velocity along the gracilis collateral network. Interestingly, the percent increase in the velocity [102.9% AE 9.5% and 99.7%AE 11.9% for the muscular branch (left) and saphenous artery (right) entrance regions, respectively, p ¼ 0.84] and shear stress (102% AE 12% and 97% AE 10% for muscular branch and saphenous artery entrance regions, respectively, p ¼ 0.73) were not different (n ¼ 10 mice). The central region of the gracilis collateral network showed the largest relative change in velocity after FAL, going from negligible flow (minimal signal prevented accurate baseline measurement) pre-FAL to high velocity post-FAL [ Fig. 4(b), white boxes]. The switch from convergent blood flow toward the central region of the gracilis muscle pre-FAL from both the saphenous and muscular branch arteries to high flow across the length of the collateral artery was present in all mice, matching the predicted hemodynamic patterns seen in Fig. 1.

Conclusions
We report here the development and application of an adapted LSF method for measuring blood velocity in collateral arteries in response to FAL in the widely used mouse ischemic hindlimb model. By adapting the LSF technique to incorporate a transillumination configuration, we substantially reduced the attenuation of the speckle signal with an increasing depth, which has previously limited the usefulness of the technique to very select instances, without impacting the relationship for determining relative velocity. In addition, noise and signal saturation issues were significantly improved by switching to trans-illumination. Ultimately, we were able to harness the power of intravital LSF microscopy to generate the first measurements of altered shear stress and velocity in individual collateral arterioles. Because the software and hardware components necessary for the technique are widely available and present a minimal financial barrier for integration into existing intravital imaging setups, this adapted technique could rapidly expand to additional models and provide similarly difficult-to-obtain hemodynamic data at the individual microvessel level. Going forward, the ability to capture these shear stress changes in vivo using our adapted LSF technology may have important bearing on determining how shear stress magnitude and direction regulate collateral artery remodeling. Our results indicate there are at least three distinct hemodynamic conditions: a nonreversing increase in the shear stress near the feeding entrance to the collateral loop, an increase in the shear stress from low/oscillating flow to sustained high shear stress at the central anastomotic region, and an increase in the shear stress but in a reversed direction at the downstream outlet back into the occluded arterial tree. Future studies may be aimed at determining whether these different flow regimes elicit differences in endothelial cell proliferation, planar polarity, and signaling, as well as differences in collateral arteriogenesis.