Understanding cerebral microvascular flow regulation and capacity to deliver oxygen to cortical tissue is of critical importance in assessing brain physiology and pathophysiology.1 In addition, cerebral cortical microvascular networks play a vital role in delivering oxygen and nutrients to support brain cognitive functions.2 Recent studies suggest that heterogeneity in capillary flow either passively correlates with3 or actively regulates4 oxygen extraction.2,5 To further understand the properties of the cerebral microcirculation, optical techniques,—for example, two-photon microscopy (TPM) and more recently, optical coherence tomography (OCT)—have been advanced to investigate the highly heterogeneous and dynamic capillary properties at high spatiotemporal resolution, and preferably over many capillaries simultaneously for strengthening the statistical power of the analysis.6,7
The OCT-based RBC-passage technique detects reflected optical intensity changes due to RBCs passing through the optical focal volume.7 Because of the high spatial sampling rate of OCT, this technique has the ability of measuring absolute RBC-flux of multiple capillaries simultaneously,7,8 which makes it suitable for characterizing the distribution of capillary flux in microvascular networks and neurocapillary coupling.9,10 However, an insufficient temporal resolution () of this RBC-passage technique may fail in capturing RBCs with very high speed when passing the focal volume, inducing aliasing in the measurements, consequently confounding physiological interpretation of the results. The aim of this study is to investigate how impacts the accuracy in estimating the capillary flux and to provide general guidance for the requirement when measuring the distribution of flux and/or mean flux of cortical capillaries in rodents by the OCT RBC-passage technique. Our simulations show that underestimates the flux with an underestimation of , at the true flux (synthetized with ) of ; in contrast, and 4.5 ms cause an underestimation of and , respectively. At the true flux of , , 3, and 4.5 cause an underestimation of , , and , respectively. Furthermore, experimental data with either or 1.5 ms were acquired in two separate groups of mice for evaluating the sufficiency of . The experimental results show that with the measured flux was underestimated by at compared with the flux measured with , which was the maximum flux that was measured in our experiments. This observed 5% underestimation of flux at is very close to our simulation result. The downsampled experimental fluxes with and 4.5 ms showed the compressed distribution of flux and the reduction in mean flux, which we also observed in the simulations. Other studies reported that a resting capillary flux of in rodents was very rare,11–14 therefore, is expected to permit accurately characterizing the distribution of flux and/or the mean flux in cortical capillaries.
Materials and Methods
All animal experimental procedures were reviewed and approved by the Massachusetts General Hospital Subcommittee on Research Animal Care. We used CD1-Elite mice (Charles River Laboratories, , male, 3 to 4 months old, ) in this study. All surgical procedures, including the sealed cranial window preparation and the cannulation of femoral artery, were conducted under 1.5% to 2% isoflurane anesthesia in mixture. In experiments, mice () were separated into two groups: group 1 () and group 2 (). Acquisition was conducted under -chloralose administered intravenously.10 Body temperature was maintained at throughout all experiments. In experiments, the mean systemic arterial blood pressure and arterial partial pressure of and averaged over mice were , 96, and 35 mmHg, respectively; the pH of systemic arterial blood was .
Spectral-Domain Optical Coherence Tomography
Capillary RBC-passage B-scans were conducted using an optimized commercial spectral-domain OCT (SD-OCT) system (Thorlabs Inc.).7 The light source has a central wavelength of 1310 nm with a 170-nm bandwidth yielding an axial resolution of in brain tissue. The transverse resolution was when using the objective (). This SD-OCT system has a maximum acquisition rate of .
Determination of the Red Blood Cell Flux
The details of the RBC-passage technique have been described elsewhere.7 Figures 1(a) and 1(b) help explain the acquisition protocol of the RBC-passage B-scan. Figure 1(a) displays a picture of a typical cranial window over mouse cortex. As shown by the red lines in Fig. 1(a), three lines were typically selected in the cranial window in each mouse for conducting the OCT RBC-passage B-scans. Each RBC-passage B-scan consisted of 30 or 60 A-lines and spanned a distance of 50 or in the cross-sectional plane, which yields or , respectively. Each line was consecutively scanned for 2 min. Because the OCT speckle fluctuation in rodent cortex is mainly due to the motion of RBCs,15 capillaries [bright pixels in Fig. 1(b)] could be identified and manually selected based on the B-scan variance images. For each selected capillary, the pixel intensities at the same coordinate of the consecutive B-scans were extracted to reveal the RBC-passage time course. Here, we used a custom-developed MATLAB graphic user interface to make sure the selected points were in capillaries by selecting points that had the largest intensity fluctuations. Next, the RBC-passage time courses were denoised by a temporal 1-D Gaussian filter with standard deviation to reduce false-positive noise peaks when counting the RBC passages. The relates to the FWHM as , so, would yield a FWHM of . The RBC-flux of corresponds to an RBC passing every 5 ms in average. Thus, using a Gaussian filter with , we expect to diminish the noise but not mask the RBC passages when RBCs pass every 5 ms. Such a high flux of was observed, which was driven by high RBC linear density instead of high speed.14 The RBC-passage time courses were then normalized with the range from 0 to 1 (i.e., the minimum value was subtracted from the time course then the time course was normalized by the maximum value). A MATLAB function “findpeaks” was used to count the RBC-passage peaks by employing a threshold set at the half-maximum intensity (i.e., the threshold was 0.5 after the time courses were normalized by the maximum value) of each RBC-passage time course. The number of RBC-passage peaks was averaged over the 2-min scan for each capillary to obtain a mean flux.
The performance of the RBC-passage peak counting procedure was evaluated by simulations. Briefly, speckle noise (filtered from experimental data) was added to synthetized RBC passages (described in Sec. 2.4). Peak counting was performed on those RBC passages with and without filtering (). Without filtering, the ratio of the estimated flux to true flux was due to false detections associated with noise. With filtering, the ratio decreased to .
Numerical Modeling of Red Blood Cell Flux
We modeled the intensity change due to RBC passage in a virtual capillary by a Gaussian function: , where represents a time point in an RBC-passage time course; is the center time point of a single RBC passage. According to the previous study,7 the RBC-speed () relates to the RBC-diameter () and the Gaussian FWHM by . Here, we set in our simulations to .7 As reported,14 the flux () relates to and the distance between neighboring RBCs () as . Then we numerically investigated the RBC-passage time courses for a range of and . Specifically, 40 were evenly selected between 0.1 and , which covers the range of capillary reported in the literature.11–14,16–26 The mean value of was reported to be .14 In simulation, 20 different were evenly chosen between 12 and . For each RBC-passage time course, is fixed for all the RBCs in that time course; and the specific between neighboring RBCs was uniformly distributed around a certain mean . So, combinations of 40 and 20 result in RBC-passage time courses. Every RBC-passage time course was generated for 2 min of simulation time with . Each of the 800 RBC-passage time courses was regenerated 100 times with different random RBC distributions for averaging. In order to investigate the impact of , each of the RBC-passage time courses with was downsampled to , 3, and 4.5 ms. With these four different ’s, a total of RBC-passage time courses were obtained; each of them yielded a mean flux . The from the individual RBC-passage time courses was binned for comparing results with different ’s.
Representative RBC-passage time courses (0.5-s trace) of four different capillaries are shown in Figs. 1(c)–1(f). For the four representative capillaries, the mean fluxes averaged over the 2-min scan are 91, 60, 56, and , respectively.
Effect of B-Scan Temporal Resolution: Numerical Modeling
We first investigate the impact of on the estimation of flux by numerical modeling. Figures 2(a)–2(c) show three synthetized RBC-passage time courses (0.2-s trace). In these three time courses, the mean was fixed at ; but the individual between neighboring RBCs were randomly varied around . In Figs. 2(a)–2(c), the is 2, 1, and , and the resultant flux is 100, 50, and , respectively. The RBC-passage time courses with are in red; and the in blue. This illustration indicates that might be appropriate for accurate estimation of flux when the but underestimates flux when .
To quantify this ceiling effect with simulations, we synthetized RBC-passage time courses with different combinations of and (see details in Sec. 2.4). The comparison between the simulated fluxes with and the fluxes with , 1.5, 3, and 4.5 ms is presented in Fig. 2(d). Here, the -axis is the simulated flux with ; the -axis is the simulated fluxes with , 1.5, 3, and 4.5 ms. The green line shows flux with , which serves as a reference for comparison. Ideally, the estimated flux should match the true flux. We expect that as becomes longer that the estimated flux will underestimate the true value, particularly at large flux values. Comparing with the fluxes with , Fig. 2(d) shows that with , 3, and 4.5 ms, the fluxes are underestimated by 9%, 57%, and 74%, respectively, at the true flux of . At the true flux of , , 3, and 4.5 cause an underestimation of , , and , respectively.
Effect of B-Scan Temporal Resolution: Experimental Data
In experiments, the RBC-passage B-scans were performed in mice of group 1 () with , and group 2 () with . In both groups, the measurements were obtained at cortical depths up to but mainly around . The original RBC-passage time courses of group 1 that were acquired with were downsampled to get two additional datasets with and 4.5 ms. Similarly, the original RBC-passage time courses of group 2 that were acquired with were downsampled to get three additional datasets with , 3, and 4.5 ms. We then calculated the mean flux for each RBC-passage time course with different ’s to investigate the impact of . Figure 2(e) shows fluxes measured in group 2. In Fig. 2(e), the -axis presents the experimental fluxes with , which served as the “truth” in the experimental scenario; and the -axis is the fluxes with , 1.5, 3, and 4.5 ms. Comparing with the fluxes with , Fig. 2(e) showed that fluxes with are accurately estimated with an underestimation of at , which is the maximum flux measured in our experiments. In contrast, with and 4.5 ms, the fluxes are underestimated by 19% and 42%, respectively, at the flux of .
The histograms in Figs. 2(f)–2(h) present the distribution of flux measured in 313 capillaries over mice. As shown, when becomes longer, the distribution of flux becomes more compressed and left-shifted to lower fluxes, as would be expected due to aliasing introducing a ceiling effect. The mean flux is 57, 47, and , with , 3, and 4.5 ms, respectively, showing a reduction of the mean flux with longer ’s.
Discussion and Conclusion
We investigated the impact of temporal resolution on the accuracy of estimating RBC-flux based on the OCT RBC-passage technique. Simulations show that with , the flux is underestimated by 9% for fluxes up to and becomes worse for larger fluxes. For longer ’s, aliasing contaminates the measurements and introduces a significant underestimation of flux. Specifically, with and 4.5 ms, the flux is underestimated by 57% and 74%, respectively, at the true flux of .
Capillary fluxes were experimentally measured in two groups of mice (total ) with and 0.75 ms in groups 1 and 2, respectively. The results demonstrate that with , 3, and 4.5 ms, the flux is underestimated by 5%, 19%, and 42%, respectively, at the maximum flux of that was measured with . Simulation shows a similar trend in that , 3, and 4.5 cause an underestimation of , , and , at the true flux of , respectively. The similarity of the simulation and experimental results support our interpretation that the underestimation of flux with longer ’s arises from aliasing. This aliasing results in misinterpretation of the RBC-flux characteristics, as illustrated in the histograms presented in Figs. 2(f)–2(h), revealing a compression of the distribution of flux as well as a reduction in the mean flux with longer ’s.
One may argue that shorter than 1.5 ms, e.g., , will enable accurately measuring higher fluxes. However, shorter means smaller FOV, in turn, fewer measurements. On the contrary, longer will allow a larger FOV and more measurements, which benefits a more powerful statistical analysis. Therefore, not readily pushing the speed limit of acquisition, we suggested that could be necessarily fast to measure typical capillary flux in mice cortex.
Some limitations exist in this study. First, while we found qualitative agreement between the simulations and experiments, some discrepancies can be observed between Figs. 2(d) and 2(e). The simulation in Fig. 2(d) showed a more pronounced and earlier ceiling effect with longer ’s. This could be a result of the simulation parameters, e.g., the and , being different from the experimental values. For example, might be distributed differently than the range used in our simulations. Experimentally, there could be capillaries with low and small , thus resulting in a high flux that is still resolvable despite a long . As an example, capillary flux as high as 14 was reported; but such a high flux was driven by a high RBC linear density instead of high speed (), which would be detectable with . Indeed, many studies have reported that the RBC-speed within most cortical capillaries is ,11–14,16–26 and thus, passage of a RBC would be , indicating that associated fluxes can be measured with . Another possibility for this discrepancy is due to the statistics of the measurement noise in the experimental data. Practically, the noise level increases with depth because of optical attenuation. As a result, the filter parameter should be adjusted to adapt to the varying signal-to-noise ratio. Although a constant filter parameter () has been used to avoid bias in analysis, it might result in a depth-dependent noise reduction. Therefore, the flux would likely be measured more accurately in some superficial capillaries while overestimated in deeper capillaries because of false positives resulting from stronger noise. Second, with the experimental fluxes, we simply evaluated the underestimation of the estimated flux at a downsampled relative to the flux with the fastest , as we did not know the true flux for the experimental data.
Another limitation of our RBC-passage technique is that we could not distinguish whether or not a dual-peak RBC-passage was because of random orientation of a flowing RBC or two RBCs flowing together. Such a dual-peak RBC-passage might be counted as one in estimating the flux because of the potential undersampling induced by fast RBCs; further, given the filter parameter used (), RBCs need to be spaced by more than 1.5 ms to be resolved as separate RBCs.
Last, we acknowledge that some capillaries with very high speed, e.g., thoroughfare capillaries, can exist18,27 and that the associated fluxes may not be accurately measured with . Therefore, would be more appropriate for estimating the mean flux of the majority of capillaries with many branch orders between arterioles and venules, but possibly not sufficient for the relatively rare thoroughfare capillaries.16 Note that possible thoroughfare capillaries are rare and are often located near the brain surface (e.g., 50 to );16,18,27 we thus expect that those rare cases do not have a large impact in characterizing the distribution of flux and the mean flux when measuring a large ensemble of capillaries.
In summary, we investigated the impact of temporal resolution on the accuracy of estimating capillary RBC-flux. Simulations and experiments characterized this impact. Our findings support that permits accurate measurements of RBC-flux up to . Longer ’s will result in a ceiling effect that causes underestimation of physiologically reasonable flux values. These results should serve as a general rule in determining the temporal resolution that is necessary to measure the distribution of flux and the mean flux in cortical capillaries in mice. These results are not limited to OCT but also extend to other optical scanning techniques for estimating RBC passage, including confocal microscopy and TPM.
This study was supported by the National Institutes of Health grants (R01-EB021018, P01-NS055104, and R01-NS091230). This study was also supported by the Air Force Office of Sponsored Research (AFOSR FA-9550-15-1-0473).