1.IntroductionDiffuse in vivo flow cytometry (DiFC) is an emerging optical technique that enables fluorescence detection of rare circulating cells in the bloodstream in the optically diffusive medium.1–4 Dual slope or dual ratio (DR) is a new diffuse optical technique that is designed to suppress noise in the optical signal and reduce sensitivity to superficial tissue regions.5–8 A challenge of DiFC is the contamination of the target fluorescence signal from noise, which may be associated with background autofluorescence (AF). In this work, we investigate the possibility of utilizing DR techniques to suppress this noise and AF, thus enabling better signal-to-noise ratio (SNR) of DiFC measurements. 1.1.Diffuse In Vivo Flow CytometryIn DiFC, the tissue surface is illuminated with laser light, typically delivered by an optical fiber bundle. As fluorescently labeled circulating cells pass through the DiFC field of view, a transient fluorescent peak may be detected at the surface using a collection fiber. Since DiFC uses diffuse light, it is possible to detect circulating cells from relatively deep-seated and large blood vessels several mm into tissue. Hence, the key advantage of DiFC is that it allows for the interrogation of relatively large circulating blood volumes enabling the detection of rare cells. For example, a major application of DiFC has been in mouse cancer research. Circulating tumor cells (CTCs) migrate from solid tumors, move through the circulatory system, and may form secondary tumor sites. As such, CTCs are instrumental in hematogenous cancer metastasis and are a major focus of medical research. However, CTCs are extremely rare and frequently number fewer than of blood in metastatic patients.9 In the clinic, the current gold standard method to study CTCs is via liquid biopsy (blood draws) using the Food and Drug Administration cleared system CellSearch.10,11 For example, in breast cancer patients, detected in 7.5 mL blood samples with CellSearch is associated with poor cancer prognosis.9 However, the small amount of blood analyzed (0.015% of the total peripheral blood volume) is known to yield poor sampling statistics, even assuming ideal Poisson statistics.12 Moreover, short-term temporal fluctuations of these cells may cause errors in estimating the CTC numbers based on when blood was drawn.13 Thus DiFC to enumerate CTCs continuously and non-invasively in large circulating blood volumes in vivo may help address these limitations. In mouse xenograft tumor models, we have previously performed DiFC on the mouse tail above the ventral caudal artery. These studies showed that it is possible to non-invasively sample the entire peripheral blood volume in , permitting the detection of very rare CTCs.2,14 Moreover, if paired with specific molecular contrast agents for CTCs (e.g., as those for fluorescence-guided surgery), translation of DiFC to humans may be possible.15–17 We recently demonstrated that CTCs could be labeled directly in the blood of mice using a folate-targeted fluorescent target (EC-17) and be detected externally with DiFC.18 However, two major technical challenges exist when applying DiFC to humans. First, the measured DiFC signal combines a fluorescence signal from the circulating cell and a non-specific background AF signal from the surrounding tissue. Although the non-specific background signal can be subtracted, noise contributed from the background cannot be removed. As such, the noise may obscure small-amplitude fluorescence signals. To illustrate this SNR consideration, example DiFC data showing signal peaks from detected flowing fluorescent microspheres embedded 0.75 and 1.00 mm deep in a phantom are shown in Figs. 1(e) and 1(f), respectively. Note that each peak represents individual mirco-spheres passing through DiFC field-of-view. The second technical challenge relates to the depth of blood vessels in humans. Suitable blood vessels, such as the radial artery in the wrist [Figs. 1(a)–1(c)], are expected to be about 2 to 4 mm deep (i.e., significantly deeper than in mice).19 Alignment to the radial artery can be potentially achieved in humans by placing DiFC probes in the “volar” wrist beneath the thumb, between the wrist bone and the tendon. Some people may be able to observe the vessel through the skin, facilitating alignment. Small misalignments () do not significantly affect the signal quality because of the wide sensitivity volume of diffuse light. We recently showed that individual cell fluorescence signals might be detectable up to 4 mm deep using near-infrared (NIR) fluorophores and a source–detector distance of , increased depth and larger relative non-specific background signals presents challenges for detection of weakly labeled CTCs.20 Hence methodology for reduction of the background signal and its contributed noise, such as dual slope or DR described in Sec. 1.2, would be extremely valuable for potential human translation of DiFC. 1.2.Dual RatioOne of the most basic measurements in diffuse optics is that of the single-distance (SD) reflectance (), which is measured by placing a point source and a point detector on the surface of a highly scattering medium. In practice, the detector does not measure the theoretical but instead, some optical intensity () that is proportional to the theoretical (Sec. 7.2.1). The proportionality of to is controlled by coupling coefficients (’s), which may be associated with the coupling efficiency of the optodes (sources or detectors) with the medium as well as any other multiplicative factor on . Recently, we introduced new measurement types that attempt to represent the theoretical by canceling any ’s associated with optodes. These are the dual slope7 or the DR8 data types, which in their calculation cancel any multiplicative factors associated with optodes.8,21 These measurement types essentially recover the spatial dependence of over multiple ’s. This is achieved from a symmetrical optode arrangement of sources and detectors, such as the one in Fig. 1(c). A further advantage of the dual-slope or DR data type is the reduced contribution to the signal from superficial parts of the diffuse medium.5–7 With the hypothesis that the unwanted AF contributors in DiFC measurements are mainly near the surface, we decided to explore if these data types could suppress the AF component of the DiFC measurement. We hypothesize that this AF suppression, together with the cancellation of noise through the cancellation of ’s, will improve the SNR of DR over traditional SD DiFC measurements using NIR optical wavelengths (’s). This improved SNR of DR DiFC could enable the detection of deep fluorescent targets, which is a current struggle of DiFC as discussed in Sec. 1.1. Therefore, this work investigates the potential of DR in the DiFC application and identifies the key parameters that control whether or not DR will be advantageous. 2.Methods2.1.Proposed Dual-Ratio SignalMost diffuse optical measurements recover an optical signal using a source and a detector, often each placed on a tissue surface. We consider each measurement of with a source and detector as a SD measurement (Sec 6.1). Previous DiFC work has considered these types of measurements;20 however, in this work, we introduce a DiFC measurement type based on a combination of SDs to form a DR (Sec. 6.3).7,8 DR is defined in detail in Sec. 6.3 and will be summarized here. The DR measurement is defined as follows: where the and subscripts represent measurements at long or short ’s, respectively. Additionally, the I and II subscripts show whether it is the first or second symmetric measurement in a DR. The exact geometric requirements for I and II are described in more detail in Sec. 6.3 and previous publications (where DR is referred to as dual slope, since geometric requirements for the two are the same).22 For this work, consider Fig. 1(c). In this case, , I and , I are detector A and source 1 (SD A1) and SD A2, respectively. And similarly, for detector B, , II and , II are SD B2 and SD B1, respectively.2.2.Monte-Carlo ModelsThis work presents simulation results based on a Monte-Carlo (MC) model run in Monte-Carlo Extreme (MCXLAB 2020 running in MATLAB 2023a on an NVIDIA RTX 4090).23 For all simulations, optical properties meant to represent of 810 nm were used.20 These were an absorption coefficient of , a scattering coefficient of , an anisotropy factor () of 0.9, and an index of refraction () of 1.37. The MC was run by launching photons into a volume, with grid spacing of 0.1 mm and 10 time gates ending at 10 ns. To calculate the continuous wave , the MC’s time-domain was integrated over the time gates. For the coordinate system, the surface of the medium was considered to be (positive pointing into the medium), and the center (with the edges 15 mm away in each direction) to be and . Using this coordinate system, detector A was placed at , detector B at , source 1 at , and source 2 at . In the case of the 0 mm simulation, both source and detector were placed at the origin. A separate MC was run for each of these optodes to find the fluence rate distribution for a MC source placed at the optode position. This follows the adjoint method to calculate the for a given SD set.24 In the case of sources, a MC pencil beam was used, and in the case of detectors, a MC cone beam with a numerical aperture of 0.5 was used to simulate the collection geometry of a fiber more accurately. This type of source distribution associated with the detectors is done to be closer to the condition of validity of the reciprocity theorem, which is crucial for applying the adjoint method. These ’s were used with Eq. (15) to yield the fluorescence Jacobian (),20 a key parameter used in all the simulation results as described in Appendix B. 2.3.Experimental Measurements with Diffuse Flow Cytometry2.3.1.Phantom experiments with fluorescent microspheres in-vitroTo estimate the potential performance of a DR DiFC instrument, we performed DiFC measurements in tissue-mimicking flow phantoms with NIR wavelengths (770 nm excitation and 810 nm emission).20 Briefly, we used Jade Green high intensity (Spherotech Inc., Lake Forest, Illinois, United States) cell-mimicking fluorescent microspheres (size to ) for NIR SD measurements.3 Microspheres were suspended at a concentration of in phosphate buffered saline and flowed through Tygon tubing (internal diameter of 0.25 mm; TGY-010-C, Small Parts Inc., Seattle, Washington, United States) embedded in a tissue-mimicking optical flow phantom made of high-density polyethylene. Microsphere suspensions were pumped with a syringe pump (0-2209, Harvard Apparatus, Holliston, Massachusetts, United States) at a flow rate of or average flow velocity of to approximate linear velocities in the mouse leg femoral artery.2,25 The tubing at a depth of 1.50 mm in the phantom. We performed SD DiFC measurements using a of 3 mm [Fig. 1(c) ]. Two fiber bundle probes were used, one as a source fiber and the other as a detector fiber. We collected DiFC data for the source (labeled as numbers) and detector (labeled as letters) pairs A1 and B2. To test the effect of the geometric arrangement of the probes relative to the flow tube, we performed measurements in the perpendicular or parallel directions as in Figs. 1(a) and 1(b). For illustration, a photograph of the A2 (at ) SD NIR DiFC measurement is included in Fig. 1(d), which shows a perpendicular placement of the probes to the tubing. We collected 15 min DiFC scans and processed the data as described previously.1 First, we performed background subtraction using a 5 s moving median window. Transit fluorescent peaks were measured as fluorescent microspheres were detected with the DiFC system. Here peaks were defined as having a minimum amplitude of five times the noise (standard deviation of the data in a 1 min moving window). We analyzed peak amplitude, width (i.e., shape), and noise properties. These parameters were used in Appendix C to find the important parameters of fluorescence efficiency and fluorescent needed to simulate SNR of DiFC measurements from the MC output. 2.3.2.Measurement of in vivo autofluorescence on miceAs discussed in more detail in Sec. 4.1 and Sec. 8.1.1, knowledge of the distribution of sources of AF is necessary for building a computational model of DR DiFC noise and background. To address this, we measured the AF background DiFC signal of the shaved inner thigh of recently euthanized mice: (1) at the surface of the skin and (2) at the surface of the exposed underlying muscle (without skin). We used female BALB/c mice on a low-fluorescence chow diet for this experiment. The hair of the inner left thigh of these euthanized mice was removed with depilatory cream (Nair, Church & Dwight Co., Inc., Ewing, New Jersey, United States). The skin of the right inner thigh was surgically removed, exposing the underlying muscle. We used our red DiFC instrument1 to measure the background tissue AF background SD signal in both cases. During these measurements, ultrasound gel was applied between the DiFC probe and the tissue surface to minimize the index of refraction mismatch. DiFC was performed for 15 min, and the mean background (i.e., with and without skin) was calculated from the resulting data. 3.Results3.1.Monte-Carlo Simulated Maps and SignalsIn Fig. 2, we show the simulated SNR from a fluorescent target placed anywhere on the plane. This shows that possible target positions would yield a measurable signal () and how strong the signal would be. One striking result that may be drawn from comparing the overall extent of the region with for the case with surface-weighted contributions of AF [Figs. 2(a)–2(d)] and the case with homogeneous contributions [Fig. 2(e)–2(h)]. In general, the case of homogeneous contributions of AF favors the shorter s. It also severely weakens the ability of DR to measure targets when they are deep within the medium. In fact, close examination of Figs. 2(f) and 2(g) suggests that, counter-intuitively, a shorter can measure deeper. However, this is only evident in the homogeneous AF contributors case and not in the surface-weighted case. Thus, telling us that the advantages of one measurement type over another depend on the spatial distribution of AF contributors. We may also compare the different s and the SD versus DR. For a 0 mm , we see a small bulb shape close to the surface, confirming what is expected for a co-localized source and detector. The non-zero s for SD show the typical banana shape expected in diffuse optics but with some differences. These differences arise from the small length scale used in these simulations, on the order of 1 mm so that the light does not act fully diffuse. For this reason, the MC source model matters much for the shape of the banana. To be more realistic, we modeled the source as a pencil beam and the detector as a cone with 0.5 numerical aperture. This resulted in one side of the banana (the one near the source) being deeper than the other since the pencil beam could more effectively launch photons along the direction into the medium. These observations show that the MC source condition simulated matters and should match reality. Focusing on DR, we see that it can measure deeper than every other measurement in the surface-weighted AF contributors case. However, in the homogeneous case, it does no better than the long ( of 4 mm) SD. This is primarily because DR focuses on canceling out signals from close to the surface. Therefore, it is advantageous when the confounding signal’s contributors, the AF’s chromophore concentration or efficiency, is surface weighted. This further reinforces the importance of understanding the distribution of AF in a realistic measurement case, like tissue, since it heavily affects which measurement type is preferable when detecting deep fluorescent targets. In Fig. 3, we present the expected signal profiles (in SNR) that one would measure if a fluorescent target flowed parallel to the source–detector line (in the plane). The axis shows the coordinate; in an actual DiFC measurement, this would be time with the position scaled by the velocity. Different colors represent targets flowing at different depths, and line type shows the type of distribution of contributors to AF. The gray region shows where the absolute SNR is , so if the curve drops into this region, we may say it is not measurable. Many of the conclusions drawn from Fig. 3 are similar to the ones that one may draw from Fig. 2 since these traces (Fig. 3) are simply line scans of Fig. 2. So we focus on more apparent features in these traces than the map. The first is the shape of a non-zero SD measurement. This is that of a double peak, which arises when the target first flows under the detector, making a strong signal, then under the source causing the target to make another strong signal. This is especially evident with shallow depths but is slightly present even in the broad traces simulated from deep s. Further, because of the different MC source conditions, the peak height when the target passes under the source is higher than when the target passes under the detector. Therefore, we may say that this model predicts that a non-zero SD measurement will produce a rather unique signal, with a double peak and a more significant peak height when the target passes under the source. Finally, lets examine the shape of the simulated DR signal [Figs. 3(d) and 3(h)]. In this case, the signal has positive and negative components. Here it is essential to understand that these traces are differences from a baseline measurement, so that Figs. 3(d) and 3(h) show the baseline DR subtracted from the current DR. This is detailed in Appendix A and Eq. (8). Therefore negative values in the DR signal mean that the current DR value is less than the baseline value. Since one wishes to identify the presence of a fluorescent target, it does not matter whether the signal is positive or negative, so even SNR of may be considered a signal which can identify the target. Knowing this, we see that the DR signal is unique and could further help identify if a target is genuinely detected. The DR signal is expected first to decrease as the target flows under the detector, then an increase when the target is under the source, followed by a decrease when the target finally flows under the second detector. Therefore, these results suggest that the DR (as well as the SD) signal has a unique shape that may be used to identify a true positive detection of the fluorescent target. 3.2.Comparison of Phantom Experiments and Simulated SignalsWe collected experimental DiFC data using a phantom and fluorescent microspheres as described in Sec. 2.3.1. Figure 4 shows sample normalized and smoothed (with the shaded region showing noise) fluorescent microsphere peaks (i.e., detections) flowing 1.50 mm deep in a phantom. DiFC measurements from two SD pairs, A1 and B2, are shown in perpendicular [Figs. 4(a)–4(c)] and parallel [Figs. 4(d)–4(f)] probe to tube configurations. Normalization was performed to set the mean peak maximum of Figs. 4(e) and 4(f) to one while using the same normalization factor for Figs. 4(b) and 4(c). This allows for comparing the relative amplitudes between all traces in Fig. 4. So that one may relate these normalized values to measured photomultiplier tube (PMT) currents, we provide the experimentally measured amplitudes for this data. For the parallel tube to probe configuration [Figs. 4(e)–4(f)], DiFC peak maximum amplitudes averaged 41 and 38 nA for A1 and B2, respectively. Meanwhile, for the perpendicular tube to probe configuration [Figs. 4(a) and 4(c)], the peak amplitudes averaged 16.64 and 17.22 nA for 1A and 2B, respectively. We have also overlaid normalized peaks simulated using MC (red dashed lines), using the same normalization method as the experimental data. This type of normalization ensures the mean amplitudes match experimental data for Figs. 4(e) and 4(f) but does not ensure a match for Figs. 4(b) and 4(c) or the individual amplitudes for Figs. 4(e) and 4(f). Additionally, the simulated velocity for the MC data was chosen to match the peak width with the experimental data. For this, a velocity of was used. However, since this velocity was assumed for all panels of Fig. 4, a match to the peak width for an individual panel is not ensured. Therefore, a comparison between experimental data and MC simulations should be done against features that were not made to match. These features were the amplitude of the peaks in Figs. 4(b) and 4(c), the amplitude of the lesser second peak in Figs. 4(e) and 4(f), and the overall shape of the signal. With this in mind, we see excellent agreement between the MC simulations and experimental measurements. This is particularly true regarding the doublet peak shown in Figs. 4(e) and 4(f), which is predicted by MC and observed experimentally. Note that the double peak in the parallel configuration is generated by the microsphere passing through high values of twice, once under the detector and once under the source. Additionally, the relative amplitudes of the doublet are influenced by the relationship between the numerical aperture of the source and detector, reinforcing the choice of a cone instead of a pencil beam for the detector in the MC simulation. In fact, if a pencil beam was chosen for both the source and detector MCs, the match to experimental data would be poor as the doublet peak would have the same maximum for both peaks in the doublet. 4.Discussion4.1.Expected DepthThe primary reason for this work is to explore if the DR technique may be applied to DiFC. With this in mind, we must consider what metric to use to compare measurement methods, such as SD versus DR. Since the aim of DiFC is to eventually non-invasively measure blood vessels of humans in vivo, the maximum depth that a method can detect a fluorescent target is the metric, which should be used for this comparison. One can determine this depth by examining Fig. 2 and finding the deepest colored region for each data type. However, a line scan in of Fig. 2, as shown in Fig. 5, makes the relationships more apparent. The values chosen in Fig. 5 were the centroid of the optodes (e.g., for 1A); thus the maximum measurable depth can be found by looking at where the curves in Fig. 5 drop below an SNR of one (the gray region). Table 1 summarizes these maximum measurable depths. Note that a suitable vessel to measure in humans, such as the radial artery, is located to 4 mm deep.19 Table 1Maximum measurable depths of fluorescent target for different measurement types.
Note: SNR, signal-to-noise ratio; SD, single distance; DR, dual ratio; ρ, source–detector distance; and AF, autofluorescence.Parameters: the same as Fig. 2 with assumed 5% NC noise. As is evident in Table 1, the depth of DR is heavily affected by the AF contributor distribution. Specifically, if the AF contributor distribution is surface-weighted, DR is the best measurement type, while this is not true for homogeneous. Evidence suggests that the AF contributor distribution is surface-weighted in vivo. To investigate the spatial distribution of the AF contributors, we measured the background signal of a mouse as described in Sec. 2.3.2. The mean background with skin was , whereas the mean background of the bare muscle (without skin) was , which is a reduction of 47%. These results imply that DiFC AF is weighted to superficial tissue layers, and as discussed below, from what we can find in the literature, this seems to be the case qualitatively. A literature search did not reveal a specific quantification of the spatial distribution of AF contributors. Nevertheless, the available literature on known endogenous fluorophores strongly suggests that these should be weighted more heavily in the skin. Specifically, reduced nicotinamide adenine dinucleotide, flavin adenine dinucleotide, collagen, and elastin fluoresce in the visible (violet and blue) wavelength regions26,27 and, to a lesser degree, in the red and NIR wavelength regions,28,29 the latter being most relevant for the potential use of DR DiFC in humans.20 Other significant contributors of tissue AF in the red and NIR wavelength regions are less obvious due to their rarity.30 Furthermore, AF in the red and NIR wavelength region may be more complex than visible AF because of the overall reduction in optical scattering and absorption at longer wavelengths.31 However, we found lipofuscin,32 melanin,33,34 and heme-metabolic compounds, such as porphyrins and bilirubin,35,36 to be reported as AF in the red and NIR wavelength region, all of which are present in the skin. Bilirubin and other heme-metabolic products are mostly known for contributing to liver AF.37 However, in the context of DiFC, these compounds may accumulate in the skin from hematomas in the case of injury. Thus the available literature and experimental DiFC measurements in mice (above) suggest that contributors to the background AF in the red and NIR wavelength region for potential DR DiFC are mostly in the skin as opposed to underlying tissue structures. 4.2.Consideration of NoiseThe noise and signal level used in these simulations was based on experimental data collected at 3 mm SD (Appendix C). However, one noise parameter that cannot be derived from the collected experimental data is the correlation of the noise measured by different optodes. This correlation between optode noises does not affect the SD results but does significantly affect the DR results since it affects how the DR can cancel noise. We modeled noise as Gaussian random variables (coupling coefficient random variable ) that multiplied the theoretical to yield the measured intensity random variable (Sec. 7.2). A separate independent random variable represents each optode (source or detector) and can be considered noise in the coupling, gain, power, or any other optode-specific noise. This concept, explained in detail in Appendix C, can be summarized by where using the random variable notation where the first argument is the mean and the second the variance. Since , the defined by the right-hand side of Eq. (2) is an excellent approximation of a random variable with mean and variance (). The DR cancels all ’s associated with specific optodes [Eq. (33)] but not the non-cancelable (NC) . NC is named such since it is the noise not canceled by DR. From this, we see that the parameter controlling the amount of noise that propagates into DR is , which is the fraction of the variance that is NC by DR.In all the above results in this work, was assumed to be 5%. Note that we do not know the physical origin of this noise if all SDs are acquired simultaneously, so we apply this noise to noise introduced from multiplexing. We expect that advanced multiplexing schemes may alleviate this noise (Sec. 4.3) but do not know what it is for the current system, so we assumed what we consider a reasonable value of 5%. To experiment with the effect of , we varied its value and determined the maximum measurable depth (where SNR is ) for each data type in Fig. 6. Note that all SD measurements experience the same noise regardless of ’s value, and only DR is affected. This is because the model was designed so that the total noise is the same regardless of the value. From Fig. 6, we can find the allowable , which makes DR sense deeper than the deepest SD modeled. For the simulation with surface-weighted AF contributions, this is 12%, while for the homogeneous case, it is 2%. This further emphasizes the advantage DR has in the surface-weighted AF medium. Additionally, this guides the amount of , which a system can have for DR to have an advantage over SD. 4.3.Future ImplementationThe experimental and computational data presented here were based on measurements using SD or single-ratio (SR) (not presented for brevity) measurements taken in series using our existing DiFC instrument.1,3 In practice, implementation of a DR for DiFC would require the construction of a new DiFC instrument capable of making simultaneous measurements with two sources and two detectors as in Fig. 1(c). Several optical designs would enable this. For example, we could frequency encode the two laser sources by modulating them at different frequencies and demodulating the fluorescence signals measured by the two detectors to separate the contributions from the two sources (frequency multiplexing). Likewise, we could alternately illuminate the two laser sources (S1 and S2) in an on-off pattern (time multiplexing). Assuming a peak width of 10 ms in vivo,1 this should be achievable by time multiplexing the laser output at about 1 kHz. Fast time multiplexing or frequency multiplexing would be desirable because the DR strategy would inherently cancel out signal drift and noise. DR cancelable artifacts are associated with a multiplicative factor applied to a source or detector that does not change within the multiplexing cycle (Sec. 7.2.1). This could include instrument drift or coupling of the sources and detectors to the skin surface. However, this approach cannot cancel some sources of random noise in the signal (Sec. 7.2.1 and Fig. 6). In principle, using modulated laser sources and lock-in detection for frequency multiplexing would improve the SNR of detected peaks in the presence of NC noise. Furthermore, as shown in Figs. 3 and 4, the DR DiFC measurement using a parallel configuration [Fig. 1(b)] yielded a unique double-peak shape. Therefore, we expect that using matched filters (or machine learning methods based on signal temporal shape) could improve peak detection. DR noise would also be suppressed when the two s are as different as possible.22 Therefore, custom-designed optical fiber bundles may be convenient for delivering and collecting light in a DR DiFC instrument. For instance, we can design a new DiFC optical fiber bundle that incorporates multiple source fibers with symmetric separations and large detector areas. Finally, we note that this work is the first case of dual slope/DR applied to fluorescence in general, not only to DiFC. Dual slope/DR was first developed for NIR spectroscopy (NIRS) and measurement of changes in .5,7 Now, with the methods presented in this work, we believe dual slope/DR may also be applied to diffuse fluorescence spectroscopy and tomography.27,35,38–41 Dual slope/DR should lend itself well to any application, which aims to measure changes in fluorescence (i.e., instead of as it was used for in NIRS). Thus, we open the door to future work regarding dual slope/DR fluorescence. 5.ConclusionThis work explored the feasibility of using DR for DiFC. We accomplished this by running MC models of the expected DR SNR based on noise and fluorescence parameters extracted from experimental DiFC data. From this exploration, we modify two key factors that control whether DR is beneficial to DiFC or not. The first is the distribution of AF contributors in vivo, with surface-weighted giving DR methods an advantage. Experiments on mice and a literature search suggest that the in vivo AF contributor distribution is concentrated in the skin, suggesting that DR will have an advantage over traditional SD methods. The second key parameter is the portion of NC (by DR) noise in the measurement. This noise is not associated with multiplicative optode factors (e.g., tissue coupling, source power, or detector efficiency), which are constant within a multiplexing cycle. A DiFC instrument designed for DR would need to be built to explore and address this. A future DiFC instrument may achieve the DR measurement by time- or frequency multiplexing leading to the required level of NC noise. Therefore, this work laid the groundwork to identify the key parameters of concern and aid in designing a DR DiFC instrument. 6.Appendix A: Definition of Measurement Types6.1.Single DistanceThe raw measurement from a DiFC system is the measured between a source and a detector that we also refer to as the SD here. Using this raw measurement, we define the change in the SD () as the background-subtracted as follows: where is expressed as a function of time () and is the average or over a baseline . could be found in various ways, for example: based on an initial measurement or an average over a preceding time period, regardless represents a change in over . Here the primary measurement parameter is the source–detector distance , in this work, SDs are evaluated at of either 3 or 4 mm.6.2.Single RatioNext, we define the SR, which is a ratio of s measured at long and short s: where the or subscripts represent measurements at either long or short s, respectively. Using Eq. (5), we next express the change in the SR (): which, as with , is also expressed as a function of time. The primary parameters for SR are the long and short s, which are 3 and 4 mm in this work.6.3.Dual RatioFinally, we define the DR, which is the geometric mean of two SRs related through geometric requirements. These requirements can be summarized by stating that optodes that contribute to the short for one SR contribute to the long for the other SR in a DR and vice versa. Suppose we have two SRs named and . Naming sources as numbers and detectors as letters, we can say that contains one detector (A) and two sources (1 and 2) such that the short is obtained between 1 and A and the long between 2 and A [Fig. 1(d)]. Now consider to be made of detector B and the same two sources 1 and 2 with the short being between B and 2 and the long between B and 1. Now notice that for 1 forms the short but for 1 forms the long . Similarly, for 2 forms the short but for 2 forms the long . Therefore, the geometric requirements are satisfied, and the geometric mean of and form a DR: and thus change in the DR () can we written asSimilar to SR, the parameters that are important for DR are the short and long s. However, in this case, there are two short and two long s. In this work, both short are considered to be 3 mm and both long 4 mm. 7.Appendix B: Theory7.1.Modeling Fluorescent ReflectanceTo model the detected fluorescent reflectance (), we must consider two processes: first, the delivery of power to the fluorophores and second, the transport of emitted light from the fluorophores to the detector. 7.1.1.Dependence on fluorophore positionThe expression for the power volume density () absorbed by a fluorophore [, unit of ] at the position vector of the fluorophore can be written as where is the absorption coefficient for the fluorophore (unit of ) and is the fluence rate (using the excitation optical wavelength optical properties) from a pencil beam at the source position to the fluorophore (unit of ). Next, we consider the emitted by the fluorophore [, unit of ]: where fluorescence efficiency represents the proportion of absorbed power converted to fluorescence. Thinking of the medium as voxelized with voxels with volume (), we volumetrically integrate Eq. (10) over (unit of ) to yield the power () emitted by the fluorophore at position inside the voxel [, unit of mW]:Finally, we find the fluorescent from the fluorophore within at , which is excited by a pencil beam at and detected at ( (unit of ). is obtained by multiplying by the Green’s function (using the emission optical properties) from to [, unit of ]: To simplify these expressions, we rewrite them to yield: Note that is rather difficult to calculate using methods, such as MC, so we may use the Green’s function of the (using the emission optical properties) with a source placed at [, unit of ] as an approximation according to the adjoint method.24 Additionally, we note that can be expressed as the Green’s function for [, unit of ] multiplied by the of the source (, unit of mW). These approximations and substitutions yield: We remind that and are the Green’s functions calculated at the excitation and emission wavelengths, respectively. In this work, however, we assume the optical properties at the two wavelengths to be the same. Last, we define the fluorescence Jacobian (),20 which is a helpful parameter in understanding how the depends on the spatial location of the fluorophore: In fact if and are both spatially independent, (unit of ) is directly proportional to and representative of how the spatial location of the fluorophore contributes to the detected fluorescence. 7.1.2.Fluorescent reflectance from multiple fluorophoresIt is assumed that the measured from each voxel containing a fluorophore may be linearly combined. Therefore: and if and (where the subscript indicates the ’th voxel) are spatially constant: whereThe background is a case where the sum of all the AF fluorophores from every voxel must be considered. In this case, Eq. (17) is considered for a heterogeneous distribution of AF or (i.e., AF contributions) and Eq. (18) for a homogeneous AF or . Any signal from a target fluorescent target is then added to this background using Eq. (16). This yields an expression for the fluorescent considering contributions from both the target fluorescent target at and all AF contributors: Here we have dropped the notation showing the source and detector for brevity. Notice that the background expressed as is the AF term, such that: andFor simulations these expressions, namely Eqs. (20) and (21), are used as the simulated s to obtained simulated measured using methods in Sec. 7.2. These simulated measured s can then be used with the expressions in Appendix A to obtain simulated measurement types SD and DR. 7.2.Simulating Noise and Coupling Coefficients7.2.1.Coupling coefficientsWe assume that the measured is related to the theoretical through multiplicative coupling coefficients (’s) typically associated with each optode. Notice that we distinguish between the measured and theoretical values by naming the former and the latter . As such, the equations in Appendix A show measurement types derived from measured (or simulated measured) s, whereas the expressions in Appendix B show calculation of theoretical s. The use of s in this section serves to connect and . Therefore, may take the proper units to convert the measured (measured in nA from a PMT for example) to the unit of (). Additionally, if noise is considered, may be a random variable representing the noise that confounds . Given the optode arrangement in Fig. 1(c), we have two sources (1 and 2) and two detectors (A and B). Therefore, four measurements of between a source and detector are possible: , , , and , which are related to through s as follows: where the s with 1, 2, A, or B subscripts are associated with the corresponding optode (source or detector). The s are NC coupling factors or noise not associated with a specific optode and from an unknown source. NC refers to these s not being canceled by the DR measurement, as is shown later in this section. Notice that there is also a unique (shown with Greek subscripts) for each measurement.7.2.2.NoiseNow we consider noise and thus model the s as Gaussian coupling coefficient random variables (’s). In the following, we assume and have the same units and scale (making unit-less and mean one); however, in reality, a coefficient converting to multiplies the ’s. We assume that the intensity noise is proportional to the such that relative noise is constant across different ’s, where Therefore, the intensity random variable from Eq. (2) has the mean value and variance (). This can be modeled by taking the mean of all ’s one and the sum of the three ’s variance who contribute to a single equal to . To model this, we introduce two parameters, the fraction of noise contributed from each optode () and the fraction of NC noise (). Thus we can model the noise as whereIn the simulations shown in this work, the ’s takes the form of a Gaussian distribution. 7.2.3.Propagation of coupling coefficients and noise to dual ratioThis work presents the results of SD and DR. Therefore, it is helpful to understand how these noise models affect each. For the case of SD, coupling and noise affect the measurement in the same way as is shown to be involved above. However, for DR, all coupling and noise are canceled except for the NC component. To show this, we rewrite Eq. (7) considering the optodes in Fig. 1(d) but replacing s with s and s: leading to where is the theoretical DR. Therefore, the measured and theoretical DR are equal if . Additionally, this demonstrates that DR eliminates s associated with multiplicative factors belonging to the optodes, and only the NC s remain in the DR measurement.8.Appendix C: Simulation Parameters from Experimental Data8.1.Determining Background and Peak Fluorescence CoefficientsThe experimental data from Sec. 2.3.1 was used to calculate and for both background and target ( and , respectively). This allowed for the simulated measurement which we present. The PMTs output measurement was in nA; therefore, it was assumed that the s in Eq. (23) together had the units of making have the unit of . We assumed this value to be based on detector gain (), detector area (), and (810 nm). However, since this value cancels whenever a unit-less quantity is calculated (such as SNR), we stress that what we assume for this value does not matter for the results; we simply assume one to match units between the experimental data and MC. As for MC parameters for calculation of , we assume a of , a scattering coefficient of , and a index of refraction () of 1.37. The experimental data collected for these calculations were measured at a of 3 mm and a source of 75 mW. Finally, the MC used a voxel size of and launched photons. 8.1.1.Background autofluorescenceTo find the background , the moving mean using a 1 s window (with a sample rate of 2 kHz) was found for all experimental data, with outliers removed (defining outliers as values more than three scaled mean absolute deviations from the median). Then the median of all moving mean values was found and taken as the background. Using this method, we found an average measured background at , of 170 nA or . When assuming the AF contributor distribution to be heterogeneous and surface-weighted, we model the as exponentially decaying as one goes deeper into the medium: Then rewrite Eq. (21) as follows: and solve for :With this method, we found that for the measured background AF. In the case of homogeneous AF contributors, we instead rewrite and solve Eq. (21) as which yields , assuming the optical properties used above for the heterogeneous case.These values of , whether homogeneous or heterogeneous, were used for the simulation results presented. Specifically, by implementing their values in Eqs. (21)–(22) to yield the simulated . 8.1.2.Fluorescence peak amplitudeFor the peak amplitude, we found peaks with an amplitude at least five times the noise (calculated in Sec. 8.2) and at least 1 s apart. Considering the peak amplitude as the background-subtracted measurement, we found a mean amplitude of 47.3 nA or . This was for and a target depth of 1.5 mm. Then rewriting Eq. (22) as we can find the by assuming that the peak maximum occurred when the target was at the location with the highest sensitivity (). Using this method, we found . This value, with the ones for the AF above, was used for the simulations in this paper using Eq. (22).8.2.Measurement of NoiseTo find the noise of the , the moving standard deviation using a 1 s window (with a sample rate of 2 kHz) was found for all experimental data, with outliers removed (defining outliers as values more than three scaled mean absolute deviations from the median), similar to the method used to find the background. Then the median of all moving standard deviation values was found and taken as the noise. Using this method, we found an average measured noise at of 5.29 nA or . Considering the peak amplitude found to be 47.3 nA or , the SNR at SD was found to be 8.95 on average. The [Eq. (27)] was the key parameter used to simulate noise for other distances. Given the background value of 170 nA or , this results in a of 0.031. This value was used for the noise simulations as described in Sec. 7.2 leading to the SNR results presented. Data, Materials, and Code AvailabilitySupporting code and data can be found at: github.com/DOIT-Lab/DOIT-Public/tree/master/DualRatio_FluorescenceFlowCytometry. Any additional supporting code and data are available from the authors upon reasonable request. AcknowledgmentsThis work was supported by the National Institutes of Health, grant nos. R01-EB029414 (to S.F.), R21-CA246413 (to M.N.), and R01-CA260202 (to M.N.). G.B. is also supported by grant no. K12-GM133314. The content is solely the authors’ responsibility and does not necessarily represent the official views of the awarding institutions. ReferencesX. Tan et al.,
“In vivo flow cytometry of extremely rare circulating cells,”
Sci. Rep., 9 3366 https://doi.org/10.1038/s41598-019-40143-2
(2019).
Google Scholar
R. A. Patil et al.,
“Fluorescence monitoring of rare circulating tumor cell and cluster dissemination in a multiple myeloma xenograft model in vivo,”
J. Biomed. Opt., 24 085004 https://doi.org/10.1117/1.JBO.24.8.085004 JBOPFO 1083-3668
(2019).
Google Scholar
J. Pace et al.,
“Near-infrared diffuse in vivo flow cytometry,”
J. Biomed. Opt., 27 097002 https://doi.org/10.1117/1.JBO.27.9.097002 JBOPFO 1083-3668
(2022).
Google Scholar
E. W. Zettergren et al.,
“Instrument for fluorescence sensing of circulating cells with diffuse light in mice in vivo,”
J. Biomed. Opt., 17 037001 https://doi.org/10.1117/1.JBO.17.3.037001 JBOPFO 1083-3668
(2012).
Google Scholar
A. Sassaroli, G. Blaney and S. Fantini,
“Dual-slope method for enhanced depth sensitivity in diffuse optical spectroscopy,”
J. Opt. Soc. Am. A, 36 1743
–1761 https://doi.org/10.1364/JOSAA.36.001743 JOAOD6 0740-3232
(2019).
Google Scholar
S. Fantini, G. Blaney and A. Sassaroli,
“Transformational change in the field of diffuse optics: from going bananas to going nuts,”
J. Innov. Opt. Health Sci., 13 1930013 https://doi.org/10.1142/S1793545819300131
(2019).
Google Scholar
G. Blaney et al.,
“Phase dual-slopes in frequency-domain near-infrared spectroscopy for enhanced sensitivity to brain tissue: first applications to human subjects,”
J. Biophotonics, 13 e201960018 https://doi.org/10.1002/jbio.201960018
(2020).
Google Scholar
G. Blaney, A. Sassaroli and S. Fantini,
“Method for measuring absolute optical properties of turbid samples in a standard cuvette,”
Appl. Sci., 12 10903 https://doi.org/10.3390/app122110903
(2022).
Google Scholar
S. Pang et al.,
“Circulating tumour cells at baseline and late phase of treatment provide prognostic value in breast cancer,”
Sci. Rep., 11 13441 https://doi.org/10.1038/s41598-021-92876-8
(2021).
Google Scholar
C. Alix-Panabières and K. Pantel,
“Circulating tumor cells: liquid biopsy of cancer,”
Clin. Chem., 59 110
–118 https://doi.org/10.1373/clinchem.2012.194258
(2013).
Google Scholar
S. Mader and K. Pantel,
“Liquid biopsy: current status and future perspectives,”
Oncol. Res. Treat., 40 404
–408 https://doi.org/10.1159/000478018
(2017).
Google Scholar
A. Mishra et al.,
“Ultrahigh-throughput magnetic sorting of large blood volumes for epitope-agnostic isolation of circulating tumor cells,”
Proc. Natl. Acad. Sci. U. S. A., 117 16839
–16847 https://doi.org/10.1073/pnas.2006388117
(2020).
Google Scholar
Z. Diamantopoulou et al.,
“The metastatic spread of breast cancer accelerates during sleep,”
Nature, 607 156
–162 https://doi.org/10.1038/s41586-022-04875-y
(2022).
Google Scholar
A. L. Williams et al.,
“Short-term circulating tumor cell dynamics in mouse xenograft models and implications for liquid biopsy,”
Front. Oncol., 10 601085 https://doi.org/10.3389/fonc.2020.601085
(2020).
Google Scholar
J. Y. K. Lee et al.,
“Review of clinical trials in intraoperative molecular imaging during cancer surgery,”
J. Biomed. Opt., 24 120901 https://doi.org/10.1117/1.JBO.24.12.120901 JBOPFO 1083-3668
(2019).
Google Scholar
S. Hernot et al.,
“Latest developments in molecular tracers for fluorescence image-guided cancer surgery,”
Lancet Oncol., 20 e354
–e367 https://doi.org/10.1016/S1470-2045(19)30317-1
(2019).
Google Scholar
M. Niedre,
“Prospects for fluorescence molecular in vivo liquid biopsy of circulating tumor cells in humans,”
Front. Photonics, 3 910035 https://doi.org/10.3389/fphot.2022.910035
(2022).
Google Scholar
R. A. Patil et al.,
“Fluorescence labeling of circulating tumor cells with a folate receptor-targeted molecular probe for diffuse in vivo flow cytometry,”
Mol. Imaging Biol., 22 1280
–1289 https://doi.org/10.1007/s11307-020-01505-9
(2020).
Google Scholar
V. Selvaraj and F. S. Buhari,
“Ultrasound evaluation of effect of different degree of wrist extension on radial artery dimension at the wrist joint,”
Ann. Cardiac Anaesth., 19 63 https://doi.org/10.4103/0971-9784.173022
(2016).
Google Scholar
F. Ivich et al.,
“Signal and measurement considerations for human translation of diffuse in vivo flow cytometry,”
J. Biomed. Opt., 27 067001 https://doi.org/10.1117/1.JBO.27.6.067001 JBOPFO 1083-3668
(2022).
Google Scholar
G. Blaney et al.,
“Functional brain mapping with dual-slope frequency-domain near-infrared spectroscopy,”
Proc. SPIE, 11946 1194602 https://doi.org/10.1117/12.2610664 PSISDG 0277-786X
(2022).
Google Scholar
G. Blaney, A. Sassaroli and S. Fantini,
“Design of a source–detector array for dual-slope diffuse optical imaging,”
Rev. Sci. Instrum., 91 093702 https://doi.org/10.1063/5.0015512 RSINAK 0034-6748
(2020).
Google Scholar
Q. Fang and D. A. Boas,
“Monte Carlo simulation of photon migration in 3D turbid media accelerated by graphics processing units,”
Opt. Express, 17 20178
–20190 https://doi.org/10.1364/OE.17.020178 OPEXFF 1094-4087
(2009).
Google Scholar
R. Yao, X. Intes and Q. Fang,
“Direct approach to compute Jacobians for diffuse optical tomography using perturbation Monte Carlo-based photon “replay”,”
Biomed. Opt. Express, 9 4588
–4603 https://doi.org/10.1364/BOE.9.004588 BOEICL 2156-7085
(2018).
Google Scholar
C. Hartmann et al.,
“Fluorescence detection, enumeration and characterization of single circulating cells in vivo: technology, applications and future prospects,”
Phys. Med. Biol., 63 01TR01 https://doi.org/10.1088/1361-6560/aa98f9 PHMBA7 0031-9155
(2017).
Google Scholar
M. Monici,
“Cell and tissue autofluorescence research and diagnostic applications,”
Biotechnol. Annu. Rev., 11 227
–256 https://doi.org/10.1016/S1387-2656(05)11007-2
(2005).
Google Scholar
A. C. Croce and G. Bottiroli,
“Autofluorescence spectroscopy and imaging: a tool for biomedical research and diagnosis,”
Eur. J. Histochem., 58 2461 https://doi.org/10.4081/ejh.2014.2461 EJHIE2 0391-7258
(2014).
Google Scholar
G. Hong, A. L. Antaris and H. Dai,
“Near-infrared fluorophores for biomedical imaging,”
Nat. Biomed. Eng., 1 0010 https://doi.org/10.1038/s41551-016-0010
(2017).
Google Scholar
S. L. Gibbs,
“Near infrared fluorescence for image-guided surgery,”
Quant. Imaging Med. Surg., 2 177
–187 https://doi.org/10.3978/j.issn.2223-4292.2012.09.04
(2012).
Google Scholar
G. Thomas et al.,
“Identifying the novel endogenous near-infrared fluorophore within parathyroid and other endocrine tissues,”
in Biomed. Opt.,
PTu3A.5
(2016). Google Scholar
A. N. Semenov et al.,
“The oxidation-induced autofluorescence hypothesis: red edge excitation and implications for metabolic imaging,”
Molecules, 25 1863 https://doi.org/10.3390/molecules25081863
(2020).
Google Scholar
A. D. Marmorstein et al.,
“Spectral profiling of autofluorescence associated with lipofuscin, Bruch’s membrane, and sub-RPE deposits in normal and AMD eyes,”
Investig. Ophthalmol. Vis. Sci., 43 2435
–2441
(2002).
Google Scholar
Z. Huang et al.,
“Cutaneous melanin exhibiting fluorescence emission under near-infrared light excitation,”
J. Biomed. Opt., 11 034010 https://doi.org/10.1117/1.2204007 JBOPFO 1083-3668
(2006).
Google Scholar
S. Wang et al.,
“In vivo near-infrared autofluorescence imaging of pigmented skin lesions: methods, technical improvements and preliminary clinical results,”
Skin Res. Technol., 19
(1), 20
–26 https://doi.org/10.1111/j.1600-0846.2012.00632.x
(2013).
Google Scholar
N. M. Htun et al.,
“Near-infrared autofluorescence induced by intraplaque hemorrhage and heme degradation as marker for high-risk atherosclerotic plaques,”
Nat. Commun., 8 75 https://doi.org/10.1038/s41467-017-00138-x NCAOBW 2041-1723
(2017).
Google Scholar
J. Lifante et al.,
“The role of tissue fluorescence in in vivo optical bioimaging,”
J. Appl. Phys., 128 171101 https://doi.org/10.1063/5.0021854 JAPIAU 0021-8979
(2020).
Google Scholar
S. G. Demos et al.,
“Near-infrared autofluorescence imaging for detection of cancer,”
J. Biomed. Opt., 9 587
–592 https://doi.org/10.1117/1.1688812 JBOPFO 1083-3668
(2004).
Google Scholar
D. Grosenick et al.,
“Review of optical breast imaging and spectroscopy,”
J. Biomed. Opt., 21 091311 https://doi.org/10.1117/1.JBO.21.9.091311 JBOPFO 1083-3668
(2016).
Google Scholar
Y. M. Wang et al.,
“Deep-tissue focal fluorescence imaging with digitally time-reversed ultrasound-encoded light,”
Nat. Commun., 3 928 https://doi.org/10.1038/ncomms1925 NCAOBW 2041-1723
(2012).
Google Scholar
A. D. Klose, V. Ntziachristos and A. H. Hielscher,
“The inverse source problem based on the radiative transfer equation in optical molecular imaging,”
J. Comput. Phys., 202 323
–345 https://doi.org/10.1016/j.jcp.2004.07.008 JCTPAH 0021-9991
(2005).
Google Scholar
I. Georgakoudi et al.,
“Trimodal spectroscopy for the detection and characterization of cervical precancers in vivo,”
Am. J. Obstetr. Gynecol., 186
(3), 374
–382 https://doi.org/10.1067/mob.2002.121075 AJOGAH 0002-9378
(2002).
Google Scholar
BiographyGiles Blaney is a National Institutes of Health Institutional Research and Academic Career Development Award Postdoctoral Scholar at the Diffuse Optical Imaging of Tissue Laboratory of Tufts University. He received his PhD from Tufts University, Medford, Massachusetts, in 2022, after working in the same laboratory with Prof. Sergio Fantini as his advisor. Before that, he received an undergraduate degree in mechanical engineering and physics from the Northeastern University, Boston, Massachusetts. His current research interests include diffuse optics and its possible applications within and outside of medical imaging. Fernando Ivich received his BS and MS degrees in biomedical engineering from the University of Arizona in 2017 and 2019, respectively. He is a PhD candidate at the Northeastern University Department of Bioengineering, Boston, Massachusetts, United States. Angelo Sassaroli received his PhD in physics from the University of Electro-Communications, Tokyo, Japan, in 2002. From July 2002 to August 2007, he was a research associate in the research group of Prof. Sergio Fantini at Tufts University. In September 2007, he was appointed by Tufts University as a research assistant professor. His fields of research are near-infrared spectroscopy and diffuse optical tomography. Mark Niedre received his PhD in medical physics from the University of Toronto, Toronto, Ontario, Canada, in 2004. He is a professor of bioengineering at the Northeastern University, Boston, Massachusetts, United States, and a senior member of SPIE and Optica. Sergio Fantini is a professor of biomedical engineering and principal investigator of the DOIT Lab at Tufts University. His research activities on applying diffuse optics to biological tissues resulted in about 120 peer-reviewed scientific publications and 12 patents. He co-authored with Prof. Irving Bigio, Boston University, Boston, Massachusetts, a textbook on “Quantitative Biomedical Optics” published by Cambridge University Press in 2016. He is a fellow of SPIE, Optica, and the American Institute for Medical and Biological Engineering. |
CITATIONS
Cited by 1 scholarly publication.
Signal to noise ratio
Fluorophores
Fluorescence
Tissues
Signal detection
Target detection
Near infrared