Abstract
Objective:
The aim of our study was to quantify the effect of doses delivered by a He:Ne laser on individual macrophage kinetics, tissue oxidative stress, and wound closure using real-time in vivo imaging.
Background:
Photobiomodulation has been reported to reduce tissue inflammation and accelerate wound closure; however, precise parameters of laser settings to optimize macrophage behavior have not been established. We hypothesized that quantitative and real-time in vivo imaging could identify optimal fluence for macrophage migration, reduction of reactive oxygen species, and acceleration of wound closure.
Methods:
Larval zebrafish Tg(mpeg-dendra2) were loaded with dihydroethidium for oxidative stress detection. Fish were caudal fin injured, treated with 635 nm continuous 5 mW He:Ne laser irradiation at 3, 9, or 18 J/cm2 and time-lapsed imaged within the first 120 min postinjury. Images taken 1 and 24-h postinjury were compared for percentage wound closure.
Results:
A fluence of 3 J/cm2 demonstrated significant increases in macrophage migration speed, fewer stops along the way, and greatest directed migration toward the wound. These findings were associated with a significant reduction in wound content reactive oxygen species when compared with control wounded fins. Both 3 and 9 J/cm2 significantly accelerated wound closure when compared with nonirradiated control fish.
Conclusions:
Wound macrophage activity could be manipulated by applied fluence, leading to reduced levels of wound reactive oxygen species and accelerated wound closure. The zebrafish model provides a means to quantitatively compare wound macrophage behavior in response to a variety of laser treatment parameters in real time.
Introduction
Resolution of the inflammatory phase of wound repair is highly dependent on the recruitment and function of macrophages and removal of excessive oxidative stress. 1,2 Excessive recruitment or prolonged retention results in a persistent proinflammatory state leading to deficient wound healing. 3 In diseased states, macrophage dysfunction can result in the failure to adequately remove neutrophils and large number of free radicals, such as superoxide anion (O2 −.) corresponding to oxidative stress. 4 Excess ROS may impair macrophages, delaying the initiation of tissue repair while low threshold levels may also impair a prohealing macrophage alternative activation. 5 As a result, augmenting the ability of macrophages to mitigate oxidative stress and facilitate reparative tasks is a sought-after target for therapeutic intervention 1,6 _ENREF_2.
Photobiomodulation (PBM), most commonly delivered by red laser light near the 635 nm wavelength, has been reported to enhance immune cell phagocytosis and secretion of chemical cues, crucial to oxidative stress clearance. 1,6,7 –13 Studies have further suggested that red laser light treatment promotes inflammatory resolution and wound closure. 14 –18 Extension of in vitro studies to clinical application has been promising; 19 –22 however, guidelines defining laser treatment settings of wavelength and duration are not yet established. 23 The wide range of energy fluence studied both ex vivo (1–60 J/cm2) and in vivo to enhance macrophage function suggests a need for real-time quantification of the therapeutic effects of PBM. Without such granular and in vivo contextual analyses, a consensus in defining precise parameters for low-level red laser may lack precision and reproducibility. 24,25
The motivation for our work was to establish an in vivo model to interrogate the effects of energy fluence on real-time measures of macrophage kinetics and tissue oxidative stress in the context of tissue injury. The zebrafish model presents powerful transgenic phenotypes to facilitate in vivo real-time characterization of macrophage recruitment and oxidative stress after injury. 26 –33 Insights into leukocyte migration behavior and reactive oxygen species generation at zebrafish wound sites have yielded high impact discoveries, providing novel insights translatable to mammalian systems. 31,32,34 –36 To simultaneously track individual macrophage trajectories and tissue oxidative stress levels, zebrafish transgenic for a fluorescent macrophage marker, Tg(mpeg-dendra2), underwent dihydroethidium (DHE) loading. Fish were injured at the caudal fin and treated with red 635 nm, 5 mW output power, at treatment doses supported by the literature 37 –39 0, 3, 9, and 18 J/cm2. This model enabled quantitative studies of the effects of energy fluence from a constant wave source on macrophage responses and tissue oxidative stress.
Materials and Methods
Preparation of zebrafish
The Tg(mpeg:dendra2) zebrafish line, in which the promoter mpeg was used to drive dendra2 expression in macrophages, 27 was purchased through the Zebrafish International Resource Center (Cat.#ZL10389). Housing, water quality, and spawning were performed according to standard protocols. 40,41 E3 embryo medium was reconstituted in distilled water for a 1x solution (5 mM NaCl, 0.17 mM KCl, 0.33 mM CaCl2, 0.33 mM MgSO4) with 0.05% methylene blue for standard use with larval zebrafish for incubation at 28°C in 100 mm glass dishes. As per AVMA guidelines, larvae were anesthetized using buffered Tricaine, 3-aminobenzoic acid ethyl ester at certain time points as described hereunder (168 mg/L, adjusted to pH 7.0 using NaHCO3; Sigma). Immediately following experimental procedures, larvae were humanely euthanized first with 400 mg/mL Tricaine solution and then 10% sodium hypochlorite solution as per 2013 AVMA guidelines. All animal procedures, conditions, and ethics described herein were approved by the Animal Care Committee at the University of Illinois at Chicago.
ROS probe loading
Dechorionated zebrafish larvae at 3 days postfertilization (dpf) were placed in 75 mm glass dishes with 5 mL of E3 medium. We used the following inclusion criterion for fish for study entry: demonstrated fluorescence, normal physiology, and no visible caudal fin nicks or wounds using a Zeiss inverted epifluorescent microscope. DHE (Sigma) dye was reconstituted in 1 mL of anhydrous dimethyl sulfoxide (Sigma-Aldrich) at 10 mM and stored in aliquots. Zebrafish, in groups of 8 to 15, were transferred to a 35 mm diameter glass bottom dish (MATEK) supplemented with 35 μM DHE solution in 1x Dulbecco's phosphate-buffered saline and incubated for 40 min in the absence of light, then transferred to a separate 75 mm glass dish, and washed with prewarmed E3 medium to remove excess DHE.
Wound placement
After ROS loading, zebrafish larvae were anesthetized (168 mg/L, Tricaine), immobilized (1% low melting point agarose solution), and underwent single caudal fin incisions (through 28 gauge sterile needle). Animals that did not meet uniformity in wound perimeter size (320 ± 50 μm), distance from notochord, or demonstrated small breaks or nicks in the embryonic outline were euthanized and not included in subsequent stages of the study.
Laser irradiation for time-lapse imaging
PBM laser treatment was administered through a research grade 5 mW (laser Power Source box: Zerona Medical Laser, Model #: EML-2245; Erchonia Corporation, McKinney, TX) emitting 635 nm continuous wave. Laser spot was elliptical at 0.25 cm2 (4 × 2 mm) at the delivery site (Fig. 1). Laser dosing was quantified in energy fluency (unit Joules/cm2) equivalent, Table 1. Larvae were randomly assigned to five treatment groups: (1) no laser, no injury, (2) 0 J/cm2 with injury, (3) 3 J/cm2 laser dose with injury, (4) 9 J/cm2 laser dose with injury, and (5) 18 J/cm2 laser dose with injury. Laser irradiation was initiated 5–15 min postinjury.

Photobiomodulation therapy schematic shows laser treatment with a continuous laser of 635 nm wavelength at 5 mW output power. Laser is applied directly to the glass of 35 mm glass bottom dish to treat fish. The entirety of the fish with a spot size of 0.25 cm2 is treated with either 150, 450, or 900 sec of irradiation for an energy fluence of 3, 9, and 18 J/cm2, respectively.
Laser Treatment Parameters at Target
n/a, not applicable.
Time-lapse imaging
All time-lapse images were acquired using ZEN® software with a Zeiss LSM 510 confocal microscope, as described previously. 33 Fluorescent excitations 488 nm for dendra 2 and 561 nm for ethidium were targeted at 0.80 and 0.08 mW, respectively, through ZEN software. Time-lapse images were taken in 60–75 sec intervals from up to 25–150 min postinjury. Up to eight fish were imaged per session with at least one nonwound for group baseline comparative analysis.
Analysis of time-lapse images
All time-lapse images were analyzed using open source Zirmi software in MATLAB R2015a, as described previously. 33 Quantification of cell kinetics was based on centroid positions, not whole cell morphology. Macrophage centroid positions were tracked using PhagoSight software's keyhole algorithm and visualization tools for proofing. 42 Zirmi software was used to select macrophage tracks that matched a selection criterion or of distinguishable centroid positions in >70% of all frames and categorized as static when movement did not exceed 0.9 μm. Error-prone tracks of macrophages near wound region were eliminated from analysis if they aggregated together. Macrophage migration was analyzed 30–120 min postinjury. Time intervals were defined as T1 [30–60 min postinjury], T2 [30–90 min postinjury], and, an overview time interval, T3 [30–120 min postinjury]. To analyze the extracted data more efficiently, wound regions were normalized, orienting positions with respect to the epicenter or zero position (So) of the wound margin for discrete distance to wound determinations. 33 In unwounded fish and for baseline measurements, the most distal position of the caudal fin was selected as the zero position. Zirmi software was used to manually trace wound region to isolate background subtracted DHE mean fluorescent intensity values ∼65 μm radially extended from the wound margin. 33
Fin regeneration assay
Fish were imaged using a Zeiss inverted epifluorescent microscope in the first 60 min postinjury. Immediately after, fish were treated with laser 0 (n = 8), 3 (n = 12), 9 (n = 9), and 18 (n = 10) J/cm2 doses. After treatment, fish were maintained at 28°C in E3 medium and, 24 h after injury, were anesthetized and imaged again. Total wound closure was reported as percentage of wound closure. The wound perimeter was measured using ImageJ (NIH) for 1- and 24-h fins per fish. Percent closure (
Statistical analysis
Statistical analyses were performed in GraphPad Prism and Microsoft Excel. Linear correlations were performed to derive Pearson r coefficients. p values were used for testing the hypothesis that there is no relationship between the observed phenomenon (null hypothesis). p values <0.05 were considered significant.
Results
Laser dose–response of oxidative stress after wounding
Control (Fig. 2A) untreated fish demonstrated higher pixel (red) intensity values around the margin of injury (Fig. 2B) than red laser-treated 3 J/cm2 fish (Fig. 2C, D). DHE mean intensity values were used to calculate geometric means and 95% confidence intervals at 60 and 120 min postinjury of control (n = 21) and dose groups 3 J/cm2 (n = 11), 9 J/cm2 (n = 8), and 18 J/cm2 (n = 8; Fig. 2E). At 60 min postinjury, a reverse dose trend was noted as 3 J/cm2 exhibited an overall 29% ROS reduction in wound region (11,013.4 ± 2,000.59, p = 0.03) compared with control fish (13,599.55 ± 1,563.5). In contrast, 9 J/cm2 exhibited values that were near control [14,108.78 ± 2,456.05, p = not significant (n.s.)] and 18 J/cm2 exhibited an overall 19% increase in ROS values (16,337.55 ± 2,893.64, p = n.s.) at near wound region compared with control. At 120 min postinjury, the 3 J/cm2 laser-treated group mean ROS remained significantly reduced when compared with the control group (p = 0.01), whereas the 9 J/cm2 laser-treated group had no significant effect (p = n.s.), and the 18 J/cm2laser-treated group showed significantly increased ROS (p = 0.02). In addition, at 120 min postinjury, a more significant reverse dose–response was noted as mean oxidative stress in the 3 J/cm2 dose group (10,190.45 ± 1,674.98) was significantly less than that in the 9 J/cm2 dose group (26%, 13,753 ± 2,527.06, p = 0.02) and was significantly less than that in the 18 J/cm2 dose group (37%, 16,161.02 ± 2,416.72, p < 0.001).

Laser dosing exhibits different levels of macrophage kinetics in response to wounding
We investigated the following number of available individual macrophages and their trajectories through three time periods [T1, T2, and T3]: 0 J/cm2 (n = 19; 227 macrophages tracked), 3 J/cm2 (n = 12; 196 macrophages tracked), 9 J/cm2 (n = 10; 96 macrophages tracked), and 18 J/cm2 (n = 8; 118 macrophages tracked). When compared with control, 3 J/cm2 macrophages demonstrated significantly more wound-directed migration (Fig. 3A, B). Speed (Fig. 3C) and meandering index (Fig. 3D) were used for determining persistent travel. Cumulatively, 3 J/cm2 macrophages exhibited significantly faster speeds (T3: 4.5 ± 0.7 vs. 3.9 ± 0.5 μm/min, p = 0.01) and more “straightness” or higher meandering index (T3: 0.44 ± 0.09 vs. 0.37 ± 0.09, p = 0.04) than nontreated control fish. Higher static ratio (Fig. 3E) was used to provide more insight into the pattern of travel, to provide insights into “pausing” and henceforth efficiency in travel. Compared with controls (T3: 0.177 ± 0.05, p = 0.03), 3 J/cm2 exhibited significantly less “pausing” overall (T3: 0.135 ± 0.04, p = 0.03), whereas 9 J/cm2 (T3: 0.139 ± 0.04, p = 0.1) and 18 J/cm2 (T3: 0.154 ± 0.02, p = 0.25) were only slightly less than control in a dose-dependent manner. Because persistence through speed, meandering index, and static ratio does not demonstrate a migration bias but instead reflects a consistent direction orientation, we used vector distances from macrophage centroid positions to the centroid position of the wound gap to identify and compare movement orientation to the site of injury. Higher net distances traveled toward the wound were used to indicate a preferential wound travel behavior or migration bias (Fig. 3F). Laser demonstrated significantly higher net distances toward wound in 3 J/cm2 overall (T3: 125 ± 50 μm) when compared with 9 J/cm2 (T3: 92 ± 52 μm, p = 0.8) and when compared with nontreated fish (T3: 89 ± 38 μm, p = 0.02). The 18 J/cm2 laser-treated group demonstrated an overall decrease in macrophage wound-oriented net distances (T3: 75 ± 38, p = 0.2).

Wound recruited macrophages of 0, 3, 9, 18 J/cm2 treatment groups of fish (n = 20, 12, 10, and 8, respectively) were tracked using software analysis. Representative control fish
Laser dose–response of 24-h wound closure
To examine the 24-h regenerative capacity of the caudal fin, wounds were measured 1-h (Fig. 4A, B) and 24-h postinjury (Fig. 4C, D). Single wounds were induced and randomly assigned to laser treatment. The following mean baseline wound perimeters were measured at 1-h postinjury: control 0 J/cm2 wound margin 343 ± 102 μm (n = 8), 3 J/cm2 wound margin 309 ± 64 μm (n = 12), 9 J/cm2 wound margin 318 ± 81 μm (n = 9), and 18 J/cm2 wound margin 318 ± 88 μm (n = 10). To ensure 1 h wound margin size did not influence 24 h wound size outcomes, a linear regression was performed indicating no correlation with 1 h wound size and 24 wound size per group: control (r = 0.22, p = n.s.), 3 J/cm2 (r = 0.14, p = n.s.), 9 J/cm2 (r = 0.37, p = n.s.), and 18 J/cm2 (r = 0.09, p = n.s.). All laser-treated fish exhibited accelerated wound closure outcomes 24 h postinjury when compared with nonirradiated controls. Percentage wound closure (geometric mean ± 95% confidence interval) of unirradiated controls (38% ± 8%) was compared with mean closure for each treatment group: 3 J/cm2 (51% ± 12%, p = 0.02), 9 J/cm2 (48% ± 11%, p = 0.04), and 18 J/cm2 (45% ± 12%, p = n.s.; Fig. 4E). These results indicate that a fluence of <9 J/cm2 may be an ideal therapeutic threshold for zebrafish caudal fin-accelerated wound closure.

Discussion
The zebrafish model facilitated the real-time detection of a dose–response in wound-oriented macrophage trajectories and oxidative stress. A dose of 3 J/cm2 decreased wound oxidative stress and elevated macrophage speeds and meandering index values for a resultant contribution for more “straightness” and less “pausing” in migration when compared with either nontreated or higher dose-treated groups. Significantly, a dose of 18 J/cm2 resulted in slower macrophage responses and highest oxidative stress values at the wound region. The 9 J/cm2 cell kinetic measures consistently fell between both high and low energy fluence. Both 3 and 9 J/cm2 demonstrated significantly greater wound closure than nontreated controls. Given the deleterious effects of excessive reactive oxygen species on wound closure rate observed in clinical practice, 43 our results support the need for a selection of laser dose to quantitatively optimize macrophage activity for the reduction of ROS and facilitation of wound closure.
Based on the large data sets involved in time-lapse imaging, this study focused on macrophage initial recruitment; selection of the quantitative measures of migration was important for quantifying an elicited effect on cellular kinetics. Our study outlines macrophage kinetics speeds, 3–6 μm/min on average, and meandering index, 0.4 to 0.6, that fall within the range in zebrafish literature. 27,44 However, the time postinjury is important in considering the average macrophage migration behavior. Initial recruitment 30–60 min postinjury demonstrated elevated speeds and less meandering for all zebrafish groups. These values lessened as time progressed. Laser-treated zebrafish did not appear to experience supernatural kinetic behavior; speeds, for instance, did not exceed the range of macrophage ranges previously reported. Macrophage static ratios were significantly correlated to respective speeds (for 3 J/cm2 group, r = −0.56, p < 0.01). Laser-treated macrophages with 3 J/cm2 demonstrated more movement, shown by decreased static ratio, implicating macrophages were simply more active rather than traveling at abnormally faster speeds frame-by-frame.
It is well accepted that speed and meandering index can be used for determining persistence, also defined as the period by which a cell moves in the same direction. 45 The 3 J/cm2 laser-treated group demonstrated considerably higher net distances toward wound and higher meandering index value. These macrophages demonstrated a preferential wound migration bias with less “stopping” and lower static ratios than the untreated controls, indicating more persistent or efficient wound-oriented tracks. Recently, reports have highlighted the role of enhanced levels of ATP, as a means of stimulating changes in speed and migration behaviors of immune cells. 46 –48 Since laser treatment has been linked to enhanced ATP production, 49 –52 it is possible that increased available ATP spurred increased macrophage persistence. 53 Together our findings coupled with previous reports support the notion that application of 635 nm energy at a fluence of 3 J/cm2 provides a threshold of utilizable energy that promotes greater macrophage activity, increased efficiencies of action, and a reduction in the time required for wound closure. The dose-dependent relationship may suggest that higher dose ranges are not as optimal as lower dose ranges for facilitating the resolution of early inflammatory events.
It is possible that increased macrophage cellular kinetics may play a direct role in the reduction of oxidative stress by hastening the pace of neutrophil clearance. Karu et al. reported that red laser (638.2 nm) at low energy fluence (0.01 to 0.03 J/cm2) could induce respiratory burst, important for clearing debris, in murine phagocytes in vitro measured by chemiluminescence. 10 Lu et al. demonstrated that red laser (632.8 nm) at lower energy fluence (1 and 2 J/cm2) increased mouse peritoneal macrophage phagocytosis of microspheres in vitro when compared with nonirradiated controls. 11 Additional findings demonstrate a reduction in proinflammatory macrophage behavior after red laser (660 nm at 4.5 J/cm2), 12 commonly tied to reactive oxygen species generation. It is possible that in vivo reduction of oxidative stress is a consequence of decreased neutrophil persistence and increased clearance.
Our findings agree with previous reports demonstrating enhanced wound closure after red laser treatment; however, wavelength and fluence varied. Rats treated with 660 nm at 4 J/cm2 demonstrated accelerated dermal healing at 7 and 14 days postinjury. 14 Alternatively, a review of laser treatments in rats demonstrated reductions in wound size and inflammatory infiltrates 54 ; however, these findings varied in dosing parameters (600–1000 nm; 1–60 J/cm2). Such widened parameters may be optimized for improved and perhaps more consistent therapeutic results.
We conclude that the zebrafish model provides a means to quantitatively test real-time macrophage responses to a variety of laser settings, using ROS and wound closure as clinically relevant outcome measures. Investigations examining each parameter of PBMT on macrophage and ROS behavior individually and in concert with each other in zebrafish provide a new means of defining this therapy. Previously, laser parameters tested comparatively in vivo using real-time quantitative methods have been difficult to accomplish. Future zebrafish investigations could provide direct comparisons of effects per laser-specific attribute such as constant versus intermittent energy delivery and differences in wavelength. The model can be interrogated to characterize macrophage polarization, neutrophil migration, and neutrophil and ROS clearance for optimization of the therapeutic in inflammatory conditions. Reduction of oxidative stress and facilitated tissue repair with PBMT would be clinically relevant for the reduction of ischemia–reperfusion injury after organ preservation or reperfusion injury after myocardial infarction, as well as a variety of additional macrophage-driven inflammatory conditions.
Footnotes
Acknowledgments
A.D.P. was funded by Diversifying Higher Education Faculty in Illinois Program Fellowship. This study was funded by a grant from the UIC Department of Surgery. We would like to thank Steve Shanks, Erchonia, for the use of the laser.
Author Disclosure Statement
No competing financial interests exist.
