Alteration of Blood Flow in a Venular Network by Infusion of Dextran 500: Evaluation with a Laser Speckle Contrast Imaging System

Abstract

This study examined the effect of dextran-induced RBC aggregation on the venular flow in microvasculature. We utilized the laser speckle contrast imaging (LSCI) as a wide-field imaging technique to visualize the flow distribution in venules influenced by abnormally elevated levels of RBC aggregation at a network-scale level, which was unprecedented in previous studies. RBC aggregation in rats was induced by infusing Dextran 500. To elucidate the impact of RBC aggregation on microvascular perfusion, blood flow in the venular network of a rat cremaster muscle was analyzed with a stepwise reduction of the arterial pressure (100 → 30 mmHg). The LSCI analysis revealed a substantial decrease in the functional vascular density after the infusion of dextran. The relative decrease in flow velocity after dextran infusion was notably pronounced at low arterial pressures. Whole blood viscosity measurements implied that the reduction in venular flow with dextran infusion could be due to the elevation of medium viscosity in high shear conditions (> 45 s^-1). In contrast, further augmentation to the flow reduction at low arterial pressures could be attributed to the formation of RBC aggregates (< 45 s^-1). This study confirmed that RBC aggregation could play a dominant role in modulating microvascular perfusion, particularly in the venular networks.

About the Authors:

Bumseok Namgung

Contributed equally to this work with: Bumseok Namgung, Yan Cheng Ng

Affiliation: Department of Biomedical Engineering, National University of Singapore, Singapore

Yan Cheng Ng

Contributed equally to this work with: Bumseok Namgung, Yan Cheng Ng

Affiliations Department of Biomedical Engineering, National University of Singapore, Singapore, NUS Graduate School for Integrative Sciences and Engineering, National University of Singapore, Singapore

Jeonghun Nam

Affiliation: Department of Biomedical Engineering, National University of Singapore, Singapore

Hwa Liang Leo

Affiliations Department of Biomedical Engineering, National University of Singapore, Singapore, NUS Graduate School for Integrative Sciences and Engineering, National University of Singapore, Singapore

Sangho Kim

* E-mail: [email protected]

Affiliations Department of Biomedical Engineering, National University of Singapore, Singapore, NUS Graduate School for Integrative Sciences and Engineering, National University of Singapore, Singapore, Department of Surgery, National University of Singapore, Singapore

Introduction

Microvasculature provides a large area that allows the exchange of substances between the blood stream and surrounding tissues. Thus, microcirculation plays a fundamental role in regulating tissue oxygenation, nutrients supply and eventually organ survival. As blood continuously circulates the entire human body while virtually reaching to all peripheral regions, tissue perfusion becomes an important indicator for evaluating the hemodynamic microvascular function. Accordingly, impaired tissue perfusion can potentially result from abnormal alteration in the rheological properties including red blood cell (RBC) aggregation, deformability and hematocrit. Among these factors, RBC aggregation has gathered great attention [1–5] due to its clinical relevance to various disease conditions [6–9]. Thus, it is important to understand the relation between RBC aggregation and flow resistance in the microvascular network since the latter is inversely proportional to the perfusion.

In many previous in vivo and in vitro studies, Dextran 500 has been used to induce RBC aggregation. Previous in vitro studies performed in a vertical glass tube reported that the formation of near-wall cell-free layer resulted in a reduction in the apparent viscosity [10,11]. In contrast, flow resistance could increase with RBC aggregation at low shear rates in a horizontal tube due mainly to sedimentation of RBCs [12]. Unlike the in vitro micro-tube experiments, the microvascular network has a complex geometry with multiple branches. Therefore, the effect of RBC aggregation on vascular resistance may not be adequately examined from the findings obtained in in vitro studies [13,14]. Moreover, a previous in vivo study highlighted that the tissue perfusion could decrease with RBC aggregation [3]. The increase in flow resistance due to RBC aggregation is prominent in venules due to their low shear condition which is favorable for the formation of RBC aggregates. Tissue perfusion was found to be significantly modulated at reduced arterial pressures (50 and 25 mmHg) with RBC aggregation when evaluated by the functional capillary density (FCD) [3]. However, most previous studies relied only on data obtained from a single-vessel or capillary-based analysis, thus greatly limiting the ability to highlight the effect of RBC aggregation on microvascular network flows.

Recently, in vivo studies have adopted the laser speckle contrast imaging (LSCI) technique to visualize the microvascular remodeling and hemodynamic changes [15,16]. LSCI provides high-resolution perfusion images in a relatively wide field by analyzing laser speckle patterns with time. Owing to the great advantages of LSCI including low-cost, ease of setup, and noninvasive measurement, it has been widely used in various studies on cerebral blood flow [17,18], retinal blood flow [19,20], liver microcirculation [21], and wound healing angiogenesis [15]. In addition, a recent study proposed the use of speckle contrast imaging to quantify the functional vascular density (FVD) [22] which has been used to examine changes in blood flow. In the present study, we adopted the LSCI to visualize the microvasculature of the rat cremaster muscle to obtain quantitative flow information at a network scale. Particularly, this study aimed to elucidate how the blood flow in venules at the network level is altered by the arterial pressure reduction and/or the aggregation induced by Dextran 500.

Materials and Methods

Animal preparation

All animal handling procedures in this study have been approved by the National University of Singapore Institutional Animal Care and Use Committee (approved protocol no. 098/12). A total of ten male rats (Sprague-Dawley) weighing 197 ± 19 g were used in this study. The animals were initially anesthetized with ketamine (37.5 mg/ml) and xylazine (5 mg/ml) cocktail through intraperitoneal injection (2 ml/kg). The surgery was performed on a heating pad to maintain the animal body temperature at 37°C. The animal was tracheotomized and the jugular vein was catheterized for the administration of additional anesthetic and dextran solutions during the experiment. The femoral artery was catheterized for blood sample withdrawal and real-time pressure monitoring. All catheters were heparinized with saline (30 IU/mL) solution to prevent blood clotting. The rat cremaster muscle was surgically exposed, and stretched on a customized animal platform with two heating elements to maintain the muscle temperature at 35°C. At the end of experiments, the animal was euthanized with an overdose of pentobarbital sodium.

Induction of RBC aggregation

Dextran 500 (Avg. MW = 450–550 kDa, Pharmacosmos A/S, Denmark) was first dissolved in saline (60 mg/ml). A total of 250 mg/kg was infused into rats over the course of 1–2 min through the jugular vein to achieve a plasma-dextran concentration of ~0.78% in the rat blood [4]. There was no discernable adverse reaction (e.g., visible swelling of the limbs) to the Dextran 500 infusion in the rats used in this study.

Hematocrit and viscosity measurement

A blood sample (~0.1 ml) was withdrawn from the femoral artery before and 15-min after the dextran infusion for the hematocrit measurement (Sigma 1–14 Microcentrifuge, Sigma, Germany). The whole blood viscosity was measured at 35°C over a range of shear rates ( = 11.25 − 90 s^-1) before and after the dextran infusion with a cone-and-plate viscometer (DV-II+, Brookfield Engineering Laboratories, Inc., Middleboro, MA, USA). Similarly, blood plasma viscosity was measured after extraction, from a blood sample by centrifugation (Sigma 2–6, Goettingen, Germany).

Laser speckle contrast imaging

Fig 1 shows the schematic of our intravital microscopic laser speckle contrast imaging (LSCI) system used to visualize the microvascular flow of the rat cremaster muscle. The rat was placed on the microscopic stage and allowed to stabilize over a period of ~10 min. The arterial cannulation was connected to a physiological data-acquisition system (MP 100 System with TSD 104A, BIOPAC Systems, USA) for continuous pressure monitoring during the experiments. Raw speckle images of microcirculation were acquired with an adjustable focus laser module (CPS196, 635 nm, 4.5mW, Thorlabs Inc., New Jersey, USA) and a high-sensitive sCMOS camera (Pco.edge, PCO AG, Kelheim, Germany) mounted on an intravital microscope (BX51, Olympus, Japan) with a 4X objective (UPlanFl, NA = 0.13, Olympus, Japan). The laser diode was equipped with a micromanipulator (NMN-21, Narishige, Japan) on the fixed microscopic stage platform for constant illumination.

[Figure omitted. See PDF.]

Fig 1. Schematic diagram of intravital microscopic laser speckle contrast imaging system (LSCI).

https://doi.org/10.1371/journal.pone.0140038.g001

The microscopic system is capable of providing images at 1280 x 1080 pixels with a pixel resolution (ρ_pixel) of 3.22 μm/pixel. To satisfy the Nyquist criterion [23], the speckle size (ρ_speckle) was ensured to be at least two times greater than the image pixel size (ρ_speckle ≥ 2ρ_pixel). The speckle size can be given by ρ_speckle = 2.44λ(1+M)f/#, where λ is the wavelength of the laser light, M is the magnification of the imaging system, and f/# is the f number of the system [24]. The speckle size relative to the image pixel size in the present study was ~9 pixels per speckle. The light intensity was carefully adjusted to match the dynamic range of the camera during the image recording for consistent image quality. The speckle contrast images were recorded at 10 frame/s to obtain a total of 80 images with a fixed exposure time of 10 ms.

Fig 2 shows a typical example of the microvasculature in the rat cremaster muscle visualized by LSCI. The speckle contrast (K) at each individual pixel (i, j) of the image was obtained by adopting the previously described temporal scheme [25] as follows:(1)where μ₈₀ and σ₈₀ represent the pixel-wise mean and standard deviation over 80 sequential images, respectively. The velocity (v) is inversely related to the speckle correlation time (τ_c) [26] which is related to K by:(2)where T represents the exposure time of the camera. This simplified relation was adopted to enhance the image processing performance [20,27]. All images utilized in this study were ensured to satisfy K < 0.6 in order for Eq 2 to be valid [27,28] (S1 Fig). All the image processing was performed with a custom-built image processing script using MATLAB (Mathworks, MA).

[Figure omitted. See PDF.]

Fig 2. Velocity map of venular network flow.

A: Speckle contrast image obtained at normal arterial pressure (MAP = 100 mmHg). B: Binary image for the vessel mask (Black: venules, Gray: arterioles). C: Venular vessel mask. D: Venular velocity map. Scale bar = 500 μm.

https://doi.org/10.1371/journal.pone.0140038.g002

Experimental protocol and vessel mask preparation

The mean arterial pressure was varied from 100 to 30 mmHg by partially occluding the abdominal aorta via a surgically inserted pneumatic cuff. The aortic occlusion was repeated for each pressure condition and maintained for 1 min before the recording of images. A vascular region, where the venular vessel branching was substantial and clearly focused, was selected for the recording. The field of view was then fixed during the course of the experiment after confirming the region of interest.

The functional vessel region was generally less than 10% of the entire field of view in the present study, thus an average value for the velocities over the entire field of view might not well represent the flow in the microvasculature. Hence, a vessel mask technique was utilized to isolate the venular network of interest for this study (Fig 2). A semi-automated image processing was performed on the speckle contrast image (K) acquired at the normal mean arterial pressure (MAP ~100 mmHg, Fig 2A) to obtain a representative venular network mask (Fig 2C). Accordingly, all the vessels were assumed to be functional at the normal MAP. The contrast image was subsequently converted into a binary image (Fig 2B) by utilizing an adaptive thresholding algorithm [29]. Small isolated objectives (number of pixels < 100) were considered as noise and removed from the binary image. Additionally, discontinuities in the vessels were manually compensated by collating with the original speckle contrast image in Fig 2A. Arteriolar vessels were then identified and removed from the venular network mask (Fig 2B and 2C). Concomitantly, the physical velocities were estimated by calibrating our LSCI measurements with in vitro experiments using predetermined flows in microfluidic channels. The details of the calibration experiments can be found in S1 Text. Finally, the venular network mask was overlaid onto the original velocity map to acquire the velocity map for the venular network (Fig 2D). As, also reported previously [30], there was no significant change in the diameter of the venules with a reduction in the MAP. Hence, the representative venular vessel mask obtained at the normal MAP (Fig 2C) was then applied to images obtained in all other pressure conditions to acquire the respective velocity maps.

Functional vascular density

To evaluate the impairment of vascular perfusion by a reduction in MAP and/or infusion of Dextran 500, we employed the LSCI to quantify the functional vascular density (FVD) in the rat muscle as reported in an earlier study [22]. The FVD is an indicator of identifying the functionality of vessels with blood flow since the flow can be visualized by LSCI based on the movement of RBCs [22]. The speckle contrast image was smoothened using the Gaussian filter and converted into a binary image (Figure B in S3 Fig). The binary image was then skeletonized until a single pixel remained across the vessel diameter, and three pixels at each branching point were removed to segment each vessel (Figure C in S3 Fig). Consequently, the FVD was calculated by dividing the total pixels over all vessel segments (total length of vessels) with the total pixels in the field of view. It is important to note that previous studies reported no observable change in the venular diameter during a stepwise arterial pressure reduction to 20 mmHg in the cat sartorius skeletal muscle [31] and with arterial pressure reduction to 40 mmHg in the rat spinotrapezius muscle [30]. In addition, no change in the venular diameter was observed with Dextran 500 infusion at a plasma-dextran concentration of ~0.6%. [30]. These findings are in consensus with ours when using a dextran concentration of ~0.78% (S4 Fig). There was no significant change in the venular diameter due to the dextran infusion or flow reduction.

Statistical analysis

A statistical software package (Prism 6.0, GraphPad) was used for all statistical analyses. Paired t-test with a two-tailed test was performed for comparison of two groups (before and after dextran infusion) for FVD and relative velocity decrease at each MAP, whereas unpaired t-test was performed for the viscosity comparison that was made from two separated groups of animals. P < 0.05 was considered statistically significant. Linear regression was performed to examine the statistical significance in the relation between two variables. The relation was considered to be statistically significant when the slope of regression line was significantly deviated from zero (P < 0.05). In addition, the analysis of covariance (ANCOVA) was performed to evaluate whether the difference between the slopes is significant.

Results and Discussion

Systemic conditions

The mean arterial pressures (MAP) and hematocrits are listed in Table 1. There was no statistical difference in the systemic parameters before and after dextran infusion.

[Figure omitted. See PDF.]

Table 1. Systemic parameters.

https://doi.org/10.1371/journal.pone.0140038.t001

Calibration and validation of velocity measurement

Typical examples of venular velocity profiles before and after dextran infusion at different MAPs are shown in Fig 3A. The centerline velocities (V_c) determined in this study were in agreement with the range observed in a 50-μm venule by Bishop et al. [32] at 120 mmHg (V_c = 4.5–5.8 mm/s) and 50 mmHg (V_c = 0.41–0.43 mm/s). In addition, Durussel et al. [5] reported a positive relation between RBC velocity and vessel diameter for venules up to 150 μm with corresponding velocities of up to 4 mm/s which was also observed in this study (Fig 3B), thus validating the measurements of venular velocity using our LSCI system.

[Figure omitted. See PDF.]

Fig 3. Typical example of velocity profile and distribution in the venular network.

A: Representative velocity profiles at various MAPs in a venule (ID ~95 μm). B: Centerline velocity as a function of vessel diameter at various MAPs (total 55 vessel segments). Data were obtained from the single venular network shown in A. Solid lines represent the linear regression fits. P-values represent the statistical significance in the difference between the slopes of regression lines before and after dextran infusion.

https://doi.org/10.1371/journal.pone.0140038.g003

The decrease in the V_c was apparent after dextran infusion when compared to that before dextran infusion in all pressure conditions (Fig 3A and 3B). Furthermore, the velocity profiles seemed to be more blunted with dextran infusion and decreasing MAP, which could be due mainly to aggregation. A similar trend in the flow velocity following dextran infusion was also reported in a previous study [1] which examined venular flows at a single-vessel level. Interestingly, as shown in Fig 3B, the centerline velocity in the venules appeared to decrease after dextran infusion for the entire range of vessel size in this study (D = 24–120 μm). In particular, the maximum centerline velocity in the venular network declined from 3.50 mm/s to 3.02 mm/s (D = 120 μm at MAP = 100 mmHg).

Effect of dextran infusion on functional vascular density

Fig 4A and 4B show a typical example of the speckle contrast images and the functional vascular density respectively, over a range of MAPs from 30 to 100 mmHg before and after dextran infusion. As illustrated in Fig 4A, the contrast value (K) was lower in the reduced MAP condition as compared to that at 100 mmHg. The decline in K was due to a decrease in the standard deviation of speckles (σ) attributed to the reduction in RBC flow. As expected, there was a positive relation between the functional vascular density (FVD) and MAP regardless of dextran infusion. The FVD after dextran infusion was significantly lower than that before dextran infusion at MAPs below 70 mmHg. Although the difference in FVD before and after dextran infusion was not statistically significant, the mean FVD after dextran infusion was consistently lower than that before dextran infusion over the entire range of MAPs (Fig 4B).

[Figure omitted. See PDF.]

Fig 4. Effect of arterial pressure reduction and dextran infusion on functional vascular density (FVD).

A: Typical examples of speckle contrast images at 30 and 100 mmHg before and after dextran infusion. B: FVD corresponding to the MAP reduction before and after dextran infusion. Scale bar = 500 μm. (* P < 0.5, ** P < 0.005; significant decrease due to the dextran infusion.)

https://doi.org/10.1371/journal.pone.0140038.g004

A previous in vivo study reported similar findings in the alteration of functional capillary density (FCD) by RBC aggregation [3]. It was reported that even at a lower plasma-dextran concentration (0.6% wt/vol), RBC aggregation could impede the FCD at normal MAP. Moreover, this reduction was more pronounced at a reduced MAP. In the present study, there also seemed to be a trend of reduction in FVD after dextran infusion. In particular, the significant reduction in the FVD after dextran infusion at low MAPs (< 70 mmHg) indicates that RBC aggregation has an important role in reducing RBC flux in the small vessels. It should be noted that although the LSCI system used in this study could not take into account the capillary flow (FCD) due to the limited image resolution (3.22 µm/pixel), we herein demonstrate an overall reduction in microvascular flow by dextran-induced RBC aggregation, particularly in low shear conditions. In contrast, our results suggest that the decrease in the FVD after dextran infusion did not have any apparent relation with the MAP. Since there was no notable difference in reduction of FVD between the low and high shear conditions, RBC aggregation may not have a significant role in influencing the FVD. This could also be partly attributed to the elevated RBC flux in the capillaries resulting from a reduction in the FCD at low MAPs [3]. As such, this also implies that RBC flows in capillaries (and thus the FCD) could be more sensitive to rheological alterations such as aggregation and pressure as compared to the FVD.

Effect of dextran infusion on venular network flow

Fig 5A shows a typical relative velocity map at 30 mmHg which was obtained by performing a pixel-wise comparison with that at the normal MAP (100 mmHg) (, where P = 30, 40… 90 mmHg, i and j = pixel coordinates). In addition, the relative velocities were averaged by the total number of pixels constituting the venules to evaluate the variations of the overall relative velocity in the venular network as shown in Fig 5B. Regardless of dextran infusion, the relative velocity in venules exhibited a significant decreasing trend with declining MAP. The mean relative velocity of the venular network appeared to decrease after dextran infusion. Furthermore, the relative velocity decrease was significantly augmented at low MAPs of 30 (P < 0.05) and 40 mmHg (P < 0.01). When the MAP decreased to 30 mmHg, the reduction of mean relative velocity after dextran infusion was ~20% greater than that before the infusion. The slopes of regression lines before (y = 1.075x + 0.136, R² = 0.58) and after (y = 1.422x − 0.196, R² = 0.70) dextran infusion (Fig 5B) correspondingly exhibited a significant difference (P < 0.05).

[Figure omitted. See PDF.]

Fig 5. Relative velocity in the venular network.

A: Typical example of the relative velocity decrease map at 30 mmHg. B: Relative velocity decrease at each MAP with respect to the normal arterial pressure (MAP = 100 mmHg). The solid lines represent linear regression fits before (y = 1.075x + 0.136, R² = 0.58, significance of slope: P < 0.0001) and after (y = 1.422x − 0.196, R² = 0.70, significance of slope: P < 0.0001) dextran infusion, respectively. Scale bar = 500 μm.

https://doi.org/10.1371/journal.pone.0140038.g005

Previous studies which utilized the cat gastrocnemius muscle showed a similar result in which the venous vascular resistance was found to be inversely related to the blood flow rate [33]. Conversely, this phenomenon was mainly due to the shear-dependent RBC aggregation that is inherent in the cat blood. Another study [1] that had used the rat spinotrapezius muscle with dextran induced RBC aggregation attributed this inverse relation to the enhanced bluntness of the venular velocity profile at low flow rates, which is also highlighted in the present study (Fig 3).

Effect of dextran infusion on blood viscosity

The overall increase in whole blood viscosity (μ_wb) due to the elevated plasma viscosity (μ_pl) and/or the formation of RBC aggregates after dextran infusion could probably contribute to the reduction in venular flow. Fig 6A shows the shear-thinning behavior of rat blood before (μ_{wb, before}) and after (μ_{wb, after}) dextran infusion over a range of shear rates ( = 11.25 − 90 s^-1) at 35°C. The curve fitting by the power-law model () clearly showed that μ_{wb, after} was consistently greater than μ_{wb, before}, and the difference between the two viscosities progressively increased with decreasing . To better depict the effect of RBC aggregation on the increase in blood viscosity, the values of μ_wb were normalized by μ_pl and represented as relative viscosity (μ_rel) in Fig 6B. μ_{rel, after} was significantly greater (P < 0.001) than μ_{rel, before} in low shear conditions (11.25–22.5 s^-1) whereas no significant difference was observed between them at shear rates greater than 45 s^-1. Concomitantly, the infusion of dextran increased μ_pl from 1.20 ± 0.03 to 1.52 ± 0.01 cP. This implied that the difference in μ_wb in high shear conditions ( > 45 s^-1) could be due mainly to the elevation in μ_pl with dextran infusion, whereas the elevation in μ_wb in low shear conditions ( < 45 s^-1) reflected the influence of RBC aggregation. It has been reported that arterial blood viscosity is lower than venous blood viscosity [34]. Since whole blood samples were withdrawn from the femoral artery for viscosity measurement in the present study, the apparent blood viscosity in the venular network could be higher than the viscosity shown in Fig 6.

[Figure omitted. See PDF.]

Fig 6. Shear-dependent blood viscosity.

A: Whole blood viscosity (μ_wb). B: Relative viscosity (μ_rel) to plasma viscosity (μ_pl) before and after dextran infusion at 35°C (n = 10). The solid lines represent curve fittings by the power-law model (). (μ_{wb, before}: m = 18.87, n = 0.75, R² = 0.79; μ_{wb, after}: m = 43.88, n = 0.60, R² = 0.91; μ_{rel, before}: m = 15.67 n = 0.75, R² = 0.79; μ_{rel, after}: m = 28.86, n = 0.60, R² = 0.91) (* P < 0.001, ** P < 0.0001).

https://doi.org/10.1371/journal.pone.0140038.g006

The present study showed that a substantial reduction in the relative venular flow was discernible by decreasing MAP down to 30 mmHg (Fig 5B). This systemic pressure-dependent reduction in the venular flow with dextran infusion implied that RBC aggregation could contribute largely to the velocity reduction at such low flow conditions. As seen in the Fig 6B, RBC aggregation seemed to be the dominant factor in elevating flow resistance in the venular network [5]. In a comprehensive review on the hemodynamic effects of RBC aggregation by Baskurt and Meiselman [35], they stated that the effect of RBC aggregation on venous flow resistance could vary with different levels of RBC aggregation. A previous study by Yalcin et al. [36] reported a triphasic response of flow resistance to the level of RBC aggregation in the isolated hind limb of guinea pigs. In their study, a significant increase in flow resistance was observed at a moderate aggregation level of 97% increase in erythrocyte sedimentation rate (ESR), followed by a reduction at enhanced ESR of 136% and 162%. Subsequently, increasing ESR by 200% resulted in an increase in flow resistance again. The plasma-dextran concentration achieved in this study (~0.78%) was similar to a pathological human level [4], and venular flows were substantially decreased by reducing the MAP down to 30 mmHg. While the augmented cell-free layer formation by RBC aggregation could lead to a reduction in the resistance by lowering the effective viscosity, this could be negated by the sedimentation of RBC aggregates and the disruption of cell-free layer formation due to frequent branching in the venular network [35].

Potential limitations

In the present study, a constant exposure time (10 ms) was used for all MAPs. This could decrease the sensitivity of our LSCI-based FVD measurement as compared to previous studies [22] which had used a longer exposure time (1000 ms). However, utilizing such a long exposure time was infeasible with our low-cost system as the pixel intensity of the raw speckle image was saturated at 1000 ms. In addition, fast flow generally requires a shorter exposure time for better spatial resolution and vice versa. Therefore, the spatial resolution of LSCI could be improved by adopting the multi-exposure LSCI [37,38], in particular under low flow conditions. In spite of the potential limitations, the flow measurement of our LSCI system, based on the conservation of volume flow at bifurcations, demonstrated a relatively high accuracy with errors of 4.0 ± 2.8% regardless of MAP.

Conclusion

In summary, using LSCI, we successfully visualized and quantified the blood perfusion in venules at a network level. Our results highlighted that dextran infusion could reduce the FVD, in particular at low MAPs (30–70 mmHg). Moreover, the reduction in venular flow velocity by decreasing MAP was further impaired by dextran infusion, in particular under low shear conditions (MAP ≤ 40 mmHg). Blood viscosity analyses demonstrated that the further impairment in the venular network perfusion could be dominantly attributed to the formation of RBC aggregates in low shear conditions whereas the dextran-induced elevation of μ_pl could partly contribute to the impairment in high shear conditions.

Supporting Information

[Figure omitted. See PDF.]

S1 Fig. Speckle contrast (K) as a function of T/τ_c.

T/τ_c represented by the asymptote (simplified algorithm, Eq 3), Lorentzian (Eq 4), Gaussian (Eq 5), and alternative Gaussian (Eq 6) velocity distribution assumptions. The asymptote provides an identical contrast value when K is smaller than 0.6 (Figure A). Probability distributions of K over the entire set of speckle contrast images used in the present study. All the analyzed speckle contrast images satisfied the criterion for the use of asymptote (K < 0.6) (Figure B).

(3)

(4)

(5)

(6)

https://doi.org/10.1371/journal.pone.0140038.s001

(TIFF)

S2 Fig. Flow phantom study for obtaining a convergent factor between mean 1/τ_c and mean velocity.

Cross-sectional view of the microchannel (Figure A). Typical example of 1/τ_c map obtained by LSCI with a mean flow velocity of 8.23 mm/s in the microchannel (Figure B). Mean 1/τ_c as a function of mean velocity (Figure C). The mean velocity was calculated based on the relation () between a given volumetric flow rate (Q) and cross-sectional area of the channel (A). The solid line indicates the curve fitting of the experimental data (y = 9.62x, R² = 0.94, n = 4). Scale bar = 500 μm.

https://doi.org/10.1371/journal.pone.0140038.s002

(TIFF)

S3 Fig. Determination of functional vascular density (FVD).

Functional vascular density was defined as the total length of all vessel segments over the area of image of microvasculature. The speckle contrast image (Figure A) was converted into binary image (Figure B) and subsequently skeletonized (Figure C) to obtain the FVD. Scale bar = 500 μm.

https://doi.org/10.1371/journal.pone.0140038.s003

(TIFF)

S4 Fig. Venular diameter before and after dextran infusion at normal and reduced MAPs.

There was no significant difference in venular diameter (ID = 19.3–40.5 μm) before and after dextran infusion regardless of MAPs (n = 9).

https://doi.org/10.1371/journal.pone.0140038.s004

(TIFF)

S1 Text. Validation study.

https://doi.org/10.1371/journal.pone.0140038.s005

(DOCX)

Acknowledgments

The authors wish to thank Ms. Yacincha S Liem for her expert technical assistance. This work was supported by National Medical Research Council (NMRC)/Cooperative Basic Research Grant (CBRG)/0019/2012.

Author Contributions

Conceived and designed the experiments: BN SK. Performed the experiments: BN YCN JN. Analyzed the data: BN YCN. Contributed reagents/materials/analysis tools: BN YCN HLL. Wrote the paper: BN YCN SK.

Citation: Namgung B, Ng YC, Nam J, Leo HL, Kim S (2015) Alteration of Blood Flow in a Venular Network by Infusion of Dextran 500: Evaluation with a Laser Speckle Contrast Imaging System. PLoS ONE 10(10): e0140038. https://doi.org/10.1371/journal.pone.0140038

References

1. Bishop JJ, Nance PR, Popel AS, Intaglietta M, Johnson PC. Effect of erythrocyte aggregation on velocity profiles in venules. American journal of physiology Heart and circulatory physiology 2001; 280: H222–236. pmid:11123237

2. Bishop JJ, Popel AS, Intaglietta M, Johnson PC. Effects of erythrocyte aggregation and venous network geometry on red blood cell axial migration. American journal of physiology Heart and circulatory physiology 2001; 281: H939–950. pmid:11454601

3. Kim S, Popel AS, Intaglietta M, Johnson PC. Effect of erythrocyte aggregation at normal human levels on functional capillary density in rat spinotrapezius muscle. American journal of physiology Heart and circulatory physiology 2006; 290: H941–947. pmid:16183731

4. Ong PK, Jain S, Namgung B, Woo YI, Kim S. Cell-free layer formation in small arterioles at pathological levels of erythrocyte aggregation. Microcirculation (New York, NY: 1994) 2011; 18: 541–551.

5. Durussel JJ, Berthault MF, Guiffant G, Dufaux J. Effects of red blood cell hyperaggregation on the rat microcirculation blood flow. Acta Physiol Scand 1998; 163: 25–32. pmid:9648620

6. Lominadze D, Joshua IG, Schuschke DA. Increased erythrocyte aggregation in spontaneously hypertensive rats. American Journal of Hypertension 1998; 11: 784–789. pmid:9683038

7. Luquita A, Urli L, Svetaz MJ, Gennaro AM, Volpintesta R, Palatnik S, et al. Erythrocyte aggregation in rheumatoid arthritis: Cell and plasma factor's role. Clin Hemorheol Microcirc 2009; 41: 49–56. pmid:19136742

8. Ozanne P, Francis RB, Meiselman HJ. Red blood cell aggregation in nephrotic syndrome. Kidney international 1983; 23: 519–525. pmid:6842967

9. Ami RB, Barshtein G, Zeltser D, Goldberg Y, Shapira I, Roth A, et al. Parameters of red blood cell aggregation as correlates of the inflammatory state. American journal of physiology Heart and circulatory physiology 2001; 280: H1982–1988. pmid:11299197

10. Alonso C, Pries AR, Gaehtgens P. Time-dependent rheological behavior of blood at low shear in narrow vertical tubes. The American journal of physiology 1993; 265: H553–561. pmid:8368359

11. Cokelet GR, Goldsmith HL. Decreased hydrodynamic resistance in the two-phase flow of blood through small vertical tubes at low flow rates. Circulation research 1991; 68: 1–17. pmid:1984854

12. Alonso C, Pries AR, Gaehtgens P. Time-dependent rheological behaviour of blood flow at low shear in narrow horizontal tubes. Biorheology 1989; 26: 229–246. pmid:2605330

13. Vicaut E. Opposite effects of red blood cell aggregation on resistance to blood flow. J Cardiovasc Surg (Torino) 1995; 36: 361–368.

14. Reinke W, Johnson PC, Gaehtgens P. Effect of shear rate variation on apparent viscosity of human blood in tubes of 29 to 94 microns diameter. Circulation research 1986; 59: 124–132. pmid:3742742

15. Rege A, Thakor NV, Rhie K, Pathak AP. In vivo laser speckle imaging reveals microvascular remodeling and hemodynamic changes during wound healing angiogenesis. Angiogenesis 2012; 15: 87–98. pmid:22198198

16. Meisner JK, Sumer S, Murrell KP, Higgins TJ, Price RJ. Laser speckle flowmetry method for measuring spatial and temporal hemodynamic alterations throughout large microvascular networks. Microcirculation 2012; 19: 619–631. pmid:22591575

17. Dunn AK. Laser speckle contrast imaging of cerebral blood flow. Annals of biomedical engineering 2012; 40: 367–377. pmid:22109805

18. Dunn AK, Bolay H, Moskowitz MA, Boas DA. Dynamic imaging of cerebral blood flow using laser speckle. Journal of Cerebral Blood Flow & Metabolism 2001; 21: 195–201. pmid:11295873

19. Cheng H, Duong TQ. Simplified laser-speckle-imaging analysis method and its application to retinal blood flow imaging. Optics letters 2007; 32: 2188–2190. pmid:17671579

20. Cheng H, Yan Y, Duong TQ. Temporal statistical analysis of laser speckle images and its application to retinal blood-flow imaging. Optics Express 2008; 16: 10214–10219. pmid:18607429

21. Sturesson C, Milstein DM, Post IC, Maas AM, van Gulik TM. Laser speckle contrast imaging for assessment of liver microcirculation. Microvascular research 2013; 87: 34–40. pmid:23403398

22. White SM, George SC, Choi B. Automated computation of functional vascular density using laser speckle imaging in a rodent window chamber model. Microvasc Res 2011; 82: 92–95. pmid:21419785

23. Kirkpatrick SJ, Duncan DD, Wells-Gray EM. Detrimental effects of speckle-pixel size matching in laser speckle contrast imaging. Optics letters 2008; 33: 2886–2888. pmid:19079481

24. Boas DA, Dunn AK. Laser speckle contrast imaging in biomedical optics. Journal of biomedical optics 2010; 15: 011109-011109-011112.

25. Rege A, Murari K, Li N, Thakor NV. Imaging Microvascular Flow Characteristics Using Laser Speckle Contrast Imaging. 2010 Annual International Conference of the Ieee Engineering in Medicine and Biology Society (Embc) 2010: 1978–1981.

26. Le T, Paul JS, Ong SH. Laser Speckle Imaging for Blood Flow Analysis. In: Pham T, editor. Computational Biology; 2010: Springer New York. pp. 243–271.

27. Ramirez-San-Juan JC, Ramos-GarcÝa R, Guizar-Iturbide I, MartÝnez-Niconoff G, Choi B. Impact of velocity distribution assumption on simplified laser speckle imaging equation. Optics express 2008; 16: 3197–3203. pmid:18542407

28. Tom WJ, Ponticorvo A, Dunn AK. Efficient processing of laser speckle contrast images. Medical Imaging, IEEE Transactions on 2008; 27: 1728–1738.

29. Gonzalez RC, Woods RE. Digital Image Processing, 2nd Edition. New Jersey: Prentice Hall; 2002.

30. Bishop JJ, Nance PR, Popel AS, Intaglietta M, Johnson PC. Diameter changes in skeletal muscle venules during arterial pressure reduction. American journal of physiology Heart and circulatory physiology 2000; 279: H47–57. pmid:10899040

31. House SD, Johnson PC. Diameter and blood flow of skeletal muscle venules during local flow regulation. The American journal of physiology 1986; 250: H828–837. pmid:3706555

32. Bishop JJ, Popel AS, Intaglietta M, Johnson PC. Effect of aggregation and shear rate on the dispersion of red blood cells flowing in venules. Am J Physiol Heart Circ Physiol 2002; 283: H1985–1996. pmid:12384477

33. Cabel M, Meiselman HJ, Popel AS, Johnson PC. Contribution of red blood cell aggregation to venous vascular resistance in skeletal muscle. The American journal of physiology 1997; 272: H1020–1032. pmid:9124410

34. Mokken FC, van der Waart F, Henny CP, Goedhart P, Gelb A. Differences in peripheral arterial and venous hemorheologic parameters. Annals of hematology 1996; 73: 135–137. pmid:8841101

35. Baskurt OK, Meiselman HJ. Hemodynamic effects of red blood cell aggregation. Indian journal of experimental biology 2007; 45: 25. pmid:17249324

36. Yalcin O, Uyuklu M, Armstrong JK, Meiselman HJ, Baskurt OK. Graded alterations of RBC aggregation influence in vivo blood flow resistance. American journal of physiology Heart and circulatory physiology 2004; 287: H2644–2650. pmid:15284061

37. Parthasarathy AB, Tom WJ, Gopal A, Zhang X, Dunn AK. Robust flow measurement with multi-exposure speckle imaging. Optics Express 2008; 16: 1975–1989. pmid:18542277

38. Rege A, Murari K, Seifert A, Pathak AP, Thakor NV. Multiexposure laser speckle contrast imaging of the angiogenic microenvironment. Journal of biomedical optics 2011; 16: 056006-056006-056010.