Imaging of fluorescent yield and lifetime from multiply scattered light reemitted from random media

Near-Infrared Fluorescence-Enhanced Optical Tomography

Fluorescence-enhanced optical imaging using near-infrared (NIR) light developed for in vivo molecular targeting and reporting of cancer provides promising opportunities for diagnostic imaging. The current state of the art of NIR fluorescence-enhanced optical tomography is reviewed in the context of the principle of fluorescence, the different measurement schemes employed, and the mathematical tools established to tomographically reconstruct the fluorescence optical properties in various tissue domains. Finally, we discuss the recent advances in forward modeling and distributed memory parallel computation to provide robust, accurate, and fast fluorescence-enhanced optical tomography.

Experimental Comparison of Continuous-Wave and Frequency-Domain Fluorescence Tomography in a Commercial Multi-Modal Scanner

The performance evaluation of a variety of small animal tomography measurement approaches and algorithms for recovery of fluorescent absorption cross section has not been conducted. Herein, we employed an intensified CCD system installed in a commercial small animal CT (Computed Tomography) scanner to compare image reconstructions from time-independent, continuous wave (CW) measurements and from time-dependent, frequency domain (FD) measurements in a series of physical phantoms specifically designed for evaluation. Comparisons were performed as a function of (1) number of projections, (2) the level of preprocessing filters used to improve the signal-to-noise ratio (SNR), (3) endogenous heterogeneity of optical properties, as well as in the cases of (4) two fluorescent targets and (5) a mouse-shaped phantom. Assessment of quantitative recovery of fluorescence absorption cross section was performed using a fully parallel, regularization-free, linear reconstruction algorithm with diffusion approximation (DA) and high order simplified spherical harmonics ( SPN) approximation to the radiative transport equation (RTE). The results show that while FD measurements may result in superior image reconstructions over CW measurements, data acquisition times are significantly longer, necessitating further development of multiple detector/source configurations, improved data read-out rates, and detector technology. FD measurements with SP3 reconstructions enabled better quantitative recovery of fluorescent target strength, but required increased computational expense. Despite the developed parallel reconstruction framework being able to achieve more than 60 times speed increase over sequential implementation, further development in faster parallel acceleration strategies for near-real time and real-time image recovery and more precise forward solution is necessary.

Reconstruction of localized fluorescent target from multi-view continuous-wave surface images of small animal with lp sparsity regularization

Fluorescence diffuse optical tomography using a multi-view continuous-wave and non-contact measurement system and an algorithm incorporating the lp (0 < p ≤ 1) sparsity regularization reconstructs a localized fluorescent target in a small animal. The measurement system provides a total of 25 fluorescence surface 2D-images of an object, which are acquired by a CCD camera from five different angles of view with excitation from five different angles. Fluorescence surface emissions from five different angles of view are simultaneously imaged on the CCD sensor, thus leading to fast acquisition of the 25 images within three minutes. The distributions of the fluorophore are reconstructed by solving the inverse problem based on the photon diffusion equations. In the reconstruction process incorporating the lp sparsity regularization, the regularization term is reformulated as a differentiable function for gradient-based non-linear optimization. Numerical simulations and phantom experiments show that the use of the lp sparsity regularization improves the localization of the target and quantitativeness of the fluorophore concentration. A mouse experiment demonstrates that a localized fluorescent target in a mouse is successfully reconstructed.

Fusion of [(18)F]FDG PET with fluorescence diffuse optical tomography to improve validation of probes and tumor imaging

Purpose

Given the progress of fluorescence diffuse optical tomography (fDOT) technology, here, we study the additional benefits provided by multimodal PET/fDOT imaging by comparing the biodistribution of 2-deoxy-2-[¹⁸F]fluoro-d-glucose ([¹⁸F]FDG) in tumors with three fluorescent probes: a glucose analog, a protease activatable optical probe, and a ligand of αvβ3 integrin.

Procedures

Sequential fDOT/PET/computed tomography (CT) imaging of mice was performed with a custom multimodal mouse support that allows the subject to be transferred between the fDOT and the PET/CT scanners. Experiments were performed in xenografted tumor models derived from the human breast cancer line MDA-MB 231 and compared to ex vivo analysis.

Results

The three-dimensional signals showed that the fluorescent glucose analog is not colocalized with [¹⁸F]FDG, raising questions about its use as a surrogate probe of the PET tracer. Fusion of [¹⁸F]FDG with the other fluorescent probes showed evidence of high variability both for the protease activity and the αvβ3 integrin expression during tumor growth.

Conclusion

The added value of hybrid PET/fDOT over the two modalities was demonstrated for cross-validation of probes and for better characterization of tumor models.

Small animal fluorescence and bioluminescence tomography: a review of approaches, algorithms and technology update

Emerging fluorescence and bioluminescence tomography approaches have several common, yet several distinct features from established emission tomographies of PET and SPECT. Although both nuclear and optical imaging modalities involve counting of photons, nuclear imaging techniques collect the emitted high energy (100-511 keV) photons after radioactive decay of radionuclides while optical techniques count low-energy (1.5-4.1 eV) photons that are scattered and absorbed by tissues requiring models of light transport for quantitative image reconstruction. Fluorescence imaging has been recently translated into clinic demonstrating high sensitivity, modest tissue penetration depth, and fast, millisecond image acquisition times. As a consequence, the promise of quantitative optical tomography as a complement of small animal PET and SPECT remains high. In this review, we summarize the different instrumentation, methodological approaches and schema for inverse image reconstructions for optical tomography, including luminescence and fluorescence modalities, and comment on limitations and key technological advances needed for further discovery research and translation.

Phantom and mouse experiments of time-domain fluorescence tomography using total light approach

Phantom and mouse experiments of time-domain fluorescence tomography were conducted to demonstrate the total light approach which was previously proposed by authors. The total light approach reduces the computation time to solve the forward model for light propagation. Time-resolved temporal profiles were acquired for cylindrical phantoms having single or double targets containing indocyanine green (ICG) solutions. The reconstructed images of ICG concentration reflected the true distributions of ICG concentration with a spatial resolution of about 10 mm. In vivo experiments were conducted using a mouse in which an ICG capsule was embedded beneath the skin in the abdomen. The reconstructed image of the ICG concentration again reflected the true distribution of ICG although artifacts due to autofluorescence appeared in the vicinity of the skin. The effectiveness of the total light approach was demonstrated by the phantom and mouse experiments.

Accelerated image reconstruction in fluorescence molecular tomography using dimension reduction

With the development of charge-coupled device (CCD) camera based non-contact fluorescence molecular tomography (FMT) imaging systems, multi projections and densely sampled fluorescent measurements used in subsequent image reconstruction can be easily obtained. However, challenges still remain in fast image reconstruction because of the large computational burden and memory requirement in the inverse problem. In this work, an accelerated image reconstruction method in FMT using principal components analysis (PCA) is presented to reduce the dimension of the inverse problem. Phantom experiments are performed to verify the feasibility of the proposed method. The results demonstrate that the proposed method can accelerate image reconstruction in FMT almost without quality degradation.

In vivo validation of quantitative frequency domain fluorescence tomography

We have developed a hybrid frequency domain fluorescence tomography and magnetic resonance imaging system (MRI) for small animal imaging. The main purpose of this system is to obtain quantitatively accurate fluorescence concentration and lifetime images using a multi-modality approach. In vivo experiments are undertaken to evaluate the system. We compare the recovered fluorescence parameters with and without MRI structural a priori information. In addition, we compare two optical background heterogeneity correction methods: Born normalization and utilizing diffuse optical tomography (DOT) functional a priori information. The results show that the concentration and lifetime of a 4.2-mm diameter indocyanine green inclusion located 15 mm deep inside a rat can be recovered with less than a 5% error when functional a priori information from DOT and structural a priori information from MRI are utilized.

Information loss and reconstruction in diffuse fluorescence tomography

This paper is a theoretical exploration of spatial resolution in diffuse fluorescence tomography. It is demonstrated that, given a fixed imaging geometry, one cannot-relative to standard techniques such as Tikhonov regularization and truncated singular value decomposition-improve the spatial resolution of the optical reconstructions via increasing the node density of the mesh considered for modeling light transport. Using techniques from linear algebra, it is shown that, as one increases the number of nodes beyond the number of measurements, information is lost by the forward model. It is demonstrated that this information cannot be recovered using various common reconstruction techniques. Evidence is provided showing that this phenomenon is related to the smoothing properties of the elliptic forward model that is used in the diffusion approximation to light transport in tissue. This argues for reconstruction techniques that are sensitive to boundaries, such as L1-reconstruction and the use of priors, as well as the natural approach of building a measurement geometry that reflects the desired image resolution.

In vivo fluorescence lifetime imaging monitors binding of specific probes to cancer biomarkers

One of the most important factors in choosing a treatment strategy for cancer is characterization of biomarkers in cancer cells. Particularly, recent advances in Monoclonal Antibodies (MAB) as primary-specific drugs targeting tumor receptors show that their efficacy depends strongly on characterization of tumor biomarkers. Assessment of their status in individual patients would facilitate selection of an optimal treatment strategy, and the continuous monitoring of those biomarkers and their binding process to the therapy would provide a means for early evaluation of the efficacy of therapeutic intervention. In this study we have demonstrated for the first time in live animals that the fluorescence lifetime can be used to detect the binding of targeted optical probes to the extracellular receptors on tumor cells in vivo. The rationale was that fluorescence lifetime of a specific probe is sensitive to local environment and/or affinity to other molecules. We attached Near-InfraRed (NIR) fluorescent probes to Human Epidermal Growth Factor 2 (HER2/neu)-specific Affibody molecules and used our time-resolved optical system to compare the fluorescence lifetime of the optical probes that were bound and unbound to tumor cells in live mice. Our results show that the fluorescence lifetime changes in our model system delineate HER2 receptor bound from the unbound probe in vivo. Thus, this method is useful as a specific marker of the receptor binding process, which can open a new paradigm in the “image and treat” concept, especially for early evaluation of the efficacy of the therapy.

Noise characteristics of heterodyne/homodyne frequency-domain measurements

We theoretically develop and experimentally validate the noise characteristics of heterodyne and/or homodyne measurements that are widely used in frequency-domain diffusive imaging. The mean and covariance of the modulated heterodyne output are derived by adapting the random amplification of a temporal point process. A multinomial selection rule is applied to the result of the temporal noise analysis to additionally model the spatial distribution of intensified photons measured by a charge-coupled device (CCD), which shows that the photon detection efficiency of CCD pixels plays an important role in the noise property of detected photons. The approach of using a multinomial probability law is validated from experimental results. Also, experimentally measured characteristics of means and variances of homodyne outputs are in agreement with the developed theory. The developed noise model can be applied to all photon amplification processes.

Fluorescence imaging to quantify the fluorescent microspheres in cardiac tissue

To quantify the fluorescent microsphere (FM) content in cardiac tissue, which is an indicative of blood flow, fluorescence imaging of both sides of the pig heart slice was employed. Despite the light scattering inside the tissue and contributions from multiple tissue layers to the total emission, it is shown that the fluorescence intensity at any pixel is proportional to the FM content and the fluorescence image may be transformed to the image of the FM concentration. A convenient standard for the emission-FM concentration transformation is proposed. The approach has several advantages in comparison with the traditional "digestion & extraction" method such as: non-destructiveness, high spatial resolution, high throughput, repeatability and simplicity of operation.

Excitation-resolved fluorescence tomography with simplified spherical harmonics equations

Fluorescence tomography (FT) reconstructs the three-dimensional (3D) fluorescent reporter probe distribution inside biological tissue. These probes target molecules of biological function, e.g. cell surface receptors or enzymes, and emit fluorescence light upon illumination with an external light source. The fluorescence light is detected on the tissue surface and a source reconstruction algorithm based on the simplified spherical harmonics (SP(N)) equations calculates the unknown 3D probe distribution inside tissue. While current FT approaches require multiple external sources at a defined wavelength range, the proposed FT method uses only a white light source with tunable wavelength selection for fluorescence stimulation and further exploits the spectral dependence of tissue absorption for the purpose of 3D tomographic reconstruction. We will show the feasibility of the proposed hyperspectral excitation-resolved fluorescence tomography method with experimental data. In addition, we will demonstrate the performance and limitations of such a method under ideal and controlled conditions by means of a digital mouse model and synthetic measurement data. Moreover, we will address issues regarding the required amount of wavelength intervals for fluorescent source reconstruction. We will explore the impact of assumed spatially uniform and nonuniform optical parameter maps on the accuracy of the fluorescence source reconstruction. Last, we propose a spectral re-scaling method for overcoming the observed limitations in reconstructing accurate source distributions in optically non-uniform tissue when assuming only uniform optical property maps for the source reconstruction process.

Early detection of breast cancer: synthesis and characterization of novel target specific NIR-fluorescent estrogen conjugate for molecular optical imaging

Estrogen induced proliferation of existing mutant cells is widely understood to be the major risk determining factor in the development of breast cancer. Hence determination of the Estrogen Receptor[ER] status is of paramount importance. We have carried out the synthesis and characterization of a novel NIR fluorescent dye conjugate aimed at measuring ER+ve status in-vivo. The conjugate was synthesized by ester formation between 17-β estradiol and a cyanine dye namely: bis-1, 1-(4-sulfobutyl) indotricarbocyanine-5-carboxylic acid, sodium salt. The replacement of the sodium ion in the ester by a larger glucosammonium ion was found to enhance the hydrophilicity and reduce the toxic effect on cell lines. The excitation and emission peaks for the dye were recorded as 750 and 788 nm respectively; ideal for non-invasive optical imaging owing to minimal tissue attenuation and auto-fluorescence at these wavelengths. The dye (NIRDC1) has a significant drop in plasma-protein binding therefore leading to marked improvement in pharmacokinetic profile such as dye evacuation in comparison to ICG. In addition the dye showed enhanced fluorescence quantum yield, molar extinction coefficient and linearity in fluorescence relative to ICG. This dye can be potentially used as a target specific exogenous contrast agent in molecular optical imaging for early detection of breast cancer.

Angular domain fluorescence lifetime imaging: a tissue-like phantom study

We describe a fluorescence lifetime imaging technique employing the collimation detection capabilities of an angular filter array (AFA). The AFA accepts minimally scattered photons emitted from fluorophores up to 2 mm deep within turbid media. The technique, referred to as Angular Domain Fluorescence Lifetime Imaging (ADFLI), is described and its performance evaluated in comparison to a conventional (lens and pinhole) system. Results from a tissue-mimicking phantom demonstrated that ADFLI provides better spatial resolution and image contrast for fluorescent probes at greater depths compared to a lens and pinhole system.

Diffuse Optics for Tissue Monitoring and Tomography

This review describes the diffusion model for light transport in tissues and the medical applications of diffuse light. Diffuse optics is particularly useful for measurement of tissue hemodynamics, wherein quantitative assessment of oxy- and deoxy-hemoglobin concentrations and blood flow are desired. The theoretical basis for near-infrared or diffuse optical spectroscopy (NIRS or DOS, respectively) is developed, and the basic elements of diffuse optical tomography (DOT) are outlined. We also discuss diffuse correlation spectroscopy (DCS), a technique whereby temporal correlation functions of diffusing light are transported through tissue and are used to measure blood flow. Essential instrumentation is described, and representative brain and breast functional imaging and monitoring results illustrate the workings of these new tissue diagnostics.

Reconstruction of fluorescence/bioluminescence sources in biological medium with spatial filter

We propose a new method for reconstruction of emitting source distributions by use of a spatial filter and a successive updating process of the forward model for fluorescence/bioluminescence diffuse optical tomography. The spatial filter transforms a set of the measurement data to a single source strength at a position of interest, and the forward model is updated by use of the estimated source strengths. This updating process ignores the dispensable source positions from reconstruction according to the reconstructed source distribution, and the spatial resolution of the reconstructed image is improved. The estimated sources are also used for the reduction of artifacts induced by noises based on the singular value decomposition. Some numerical experiments show the advantages of the proposed method by comparing the present results with those obtained by the conventional methods of the least squares method nd Algebraic Reconstruction Technique. Finally the criteria for practical use of the method are quantitatively presented by the simulations for 2D and 3D geometries.

Image reconstruction of fluorescent molecular tomography based on the tree structured Schur complement decomposition

Background

The inverse problem of fluorescent molecular tomography (FMT) often involves complex large-scale matrix operations, which may lead to unacceptable computational errors and complexity. In this research, a tree structured Schur complement decomposition strategy is proposed to accelerate the reconstruction process and reduce the computational complexity. Additionally, an adaptive regularization scheme is developed to improve the ill-posedness of the inverse problem.

Methods

The global system is decomposed level by level with the Schur complement system along two paths in the tree structure. The resultant subsystems are solved in combination with the biconjugate gradient method. The mesh for the inverse problem is generated incorporating the prior information. During the reconstruction, the regularization parameters are adaptive not only to the spatial variations but also to the variations of the objective function to tackle the ill-posed nature of the inverse problem.

Results

Simulation results demonstrate that the strategy of the tree structured Schur complement decomposition obviously outperforms the previous methods, such as the conventional Conjugate-Gradient (CG) and the Schur CG methods, in both reconstruction accuracy and speed. As compared with the Tikhonov regularization method, the adaptive regularization scheme can significantly improve ill-posedness of the inverse problem.

Conclusions

The methods proposed in this paper can significantly improve the reconstructed image quality of FMT and accelerate the reconstruction process.

Near-infrared quantum dots for deep tissue imaging

Developments in nanotechnology have paved the way for the early detection, treatment, and prevention of several tumors which affect mankind. In the past few years, near-infrared (NIR) fluorescence imaging techniques have emerged that enable the in vivo imaging of physiological, metabolic, and molecular function. The NIR window, also known as the diagnostic window (700-900 nm), can be explored for sensitive detection techniques. Nanoparticles, particularly semiconductor quantum dots (QDs), can be utilized for the purpose of optical imaging. These semiconductor QDs possess novel electronic, optical, magnetic, and structural properties which are quite different from those of bulk materials. NIR QDs with these unique properties can be utilized as contrast agents for optical imaging, particularly for deep tissue imaging. Deep tissue imaging provides more information about the pathological status of the disease, which makes the treatment more effective and efficient. In this review we highlight the importance of NIR QDs as probes for optical imaging. We describe the different types of NIR QDs, their synthesis, and their application for deep tissue imaging along with recently developed self-illuminating NIR QDs.

Low-frequency wide-field fluorescence lifetime imaging using a high-power near-infrared light-emitting diode light source

Fluorescence lifetime imaging (FLi) could potentially improve exogenous near-infrared (NIR) fluorescence imaging, because it offers the capability of discriminating a signal of interest from background, provides real-time monitoring of a chemical environment, and permits the use of several different fluorescent dyes having the same emission wavelength. We present a high-power, LED-based, NIR light source for the clinical translation of wide-field (larger than 5 cm in diameter) FLi at frequencies up to 35 MHz. Lifetime imaging of indocyanine green (ICG), IRDye 800-CW, and 3,3(')-diethylthiatricarbocyanine iodide (DTTCI) was performed over a large field of view (10 cm by 7.5 cm) using the LED light source. For comparison, a laser diode light source was employed as a gold standard. Experiments were performed both on the bench by diluting the fluorescent dyes in various chemical environments in Eppendorf tubes, and in vivo by injecting the fluorescent dyes mixed in Matrigel subcutaneously into CD-1 mice. Last, measured fluorescence lifetimes obtained using the LED and the laser diode sources were compared with those obtained using a state-of-the-art time-domain imaging system and with those previously described in the literature. On average, lifetime values obtained using the LED and the laser diode light sources were consistent, exhibiting a mean difference of 3% from the expected values and a coefficient of variation of 12%. Taken together, our study offers an alternative to laser diodes for clinical translation of FLi and explores the use of relatively low frequency modulation for in vivo imaging.

Automatic exposure control and estimation of effective system noise in diffuse fluorescence tomography

A diffuse fluorescence tomography system, based upon time-correlated single photon counting, is presented with an automated algorithm to allow dynamic range variation through exposure control. This automated exposure control allows the upper and lower detection levels of fluorophore to be extended by an order of magnitude beyond the previously published performance and benefits in a slight decrease in system effective noise. The effective noise level is used as a metric to characterize the system performance, integrating both model-mismatch and calibration bias errors into a single parameter. This effective error is near 7% of the reconstructed fluorescent yield value, when imaging in just few minutes. Quantifying protoporphyrin IX concentrations down to 50 ng/ml is possible, for tumor-sized regions. This fluorophore has very low fluorescence yield, but high biological relevance for tumor imaging, given that it is produced in the mitochondria, and upregulated in many tumor types.

Pre-clinical whole-body fluorescence imaging: Review of instruments, methods and applications

Fluorescence sampling of cellular function is widely used in all aspects of biology, allowing the visualization of cellular and sub-cellular biological processes with spatial resolutions in the range from nanometers up to centimeters. Imaging of fluorescence in vivo has become the most commonly used radiological tool in all pre-clinical work. In the last decade, full-body pre-clinical imaging systems have emerged with a wide range of utilities and niche application areas. The range of fluorescent probes that can be excited in the visible to near-infrared part of the electromagnetic spectrum continues to expand, with the most value for in vivo use being beyond the 630 nm wavelength, because the absorption of light sharply decreases. Whole-body in vivo fluorescence imaging has not yet reached a state of maturity that allows its routine use in the scope of large-scale pre-clinical studies. This is in part due to an incomplete understanding of what the actual fundamental capabilities and limitations of this imaging modality are. However, progress is continuously being made in research laboratories pushing the limits of the approach to consistently improve its performance in terms of spatial resolution, sensitivity and quantification. This paper reviews this imaging technology with a particular emphasis on its potential uses and limitations, the required instrumentation, and the possible imaging geometries and applications. A detailed account of the main commercially available systems is provided as well as some perspective relating to the future of the technology development. Although the vast majority of applications of in vivo small animal imaging are based on epi-illumination planar imaging, the future success of the method relies heavily on the design of novel imaging systems based on state-of-the-art optical technology used in conjunction with high spatial resolution structural modalities such as MRI, CT or ultrasound.

Nonlinear reconstruction of absorption and fluorescence contrast from measured diffuse transmittance and reflectance of a compressed-breast-simulating phantom

We report on the nonlinear reconstruction of local absorption and fluorescence contrast in tissuelike scattering media from measured time-domain diffuse reflectance and transmittance of laser as well as laser-excited fluorescence radiation. Measurements were taken at selected source-detector offsets using slablike diffusely scattering and fluorescent phantoms containing fluorescent heterogeneities. Such measurements simulate in vivo data that would be obtained employing a scanning, time-domain fluorescence mammograph, where the breast is gently compressed between two parallel glass plates, and source and detector optical fibers scan synchronously at various source-detector offsets, allowing the recording of laser and fluorescence mammograms. The diffusion equations modeling the propagation of the laser and fluorescence radiation were solved in frequency domain by the finite element method simultaneously for several modulation frequencies using Fourier transformation and preprocessed experimental data. To reconstruct the concentration of the fluorescent contrast agent, the Born approximation including higher-order reconstructed photon densities at the excitation wavelength was used. Axial resolution was determined that can be achieved by various detection schemes. We show that remission measurements increase the depth resolution significantly.

Excitation spectroscopy in multispectral optical fluorescence tomography: methodology, feasibility and computer simulation studies

Molecular probes used for in vivo optical fluorescence tomography (OFT) studies in small animals are typically chosen such that their emission spectra lie in the 680-850 nm wavelength range. This is because tissue attenuation in this spectral band is relatively low, allowing optical photons even from deep sites in tissue to reach the animal surface and consequently be detected by a CCD camera. The wavelength dependence of tissue optical properties within the 680-850 nm band can be exploited for emitted light by measuring fluorescent data via multispectral approaches and incorporating the spectral dependence of these optical properties into the OFT inverse problem-that of reconstructing underlying 3D fluorescent probe distributions from optical data collected on the animal surface. However, in the aforementioned spectral band, due to only small variations in the tissue optical properties, multispectral emission data, though superior for image reconstruction compared to achromatic data, tend to be somewhat redundant. A different spectral approach for OFT is to capitalize on the larger variations in the optical properties of tissue for excitation photons than for the emission photons by using excitation at multiple wavelengths as a means of decoding source depth in tissue. The full potential of spectral approaches in OFT can be realized by a synergistic combination of these two approaches, that is, exciting the underlying fluorescent probe at multiple wavelengths and measuring emission data multispectrally. In this paper, we describe a method that incorporates both excitation and emission spectral information into the OFT inverse problem. We describe a linear algebraic formulation of the multiple wavelength illumination-multispectral detection forward model for OFT and compare it to models that use only excitation at multiple wavelengths or those that use only multispectral detection techniques. This study is carried out in a realistic inhomogeneous mouse atlas using singular value decomposition and analysis of reconstructed spatial resolution versus noise. For simplicity, quantitative results have been shown for one representative fluorescent probe (Alexa 700) and effects due to tissue autofluorescence have not been taken into account. We also demonstrate the performance of our method for 3D reconstruction of tumors in a simulated mouse model of metastatic human hepatocellular carcinoma.

Fluorescence lifetime imaging from time resolved measurements using a shape-based approach

We present a novel fluorescent tomography algorithm to estimate the spatial distribution of fluorophores and the fluorescence lifetimes from surface time resolved measurements. The algorithm is a hybridization of the level set technique for recovering the distributions of distinct fluorescent markers with a gradient method for estimating their lifetimes. This imaging method offers several advantages compared to more traditional pixel-based techniques as, for example, well defined boundaries and a better resolution of the images. The numerical experiments show that our imaging method gives rise to accurate reconstructions in the presence of data noise and fluorescence background even for complicated fluorophore distributions in several-centimiter-thick biological tissue.

Spectral distortion in diffuse molecular luminescence tomography in turbid media

The influence of tissue optical properties on the shape of near-infrared (NIR) fluorescence emission spectra propagating through multiple centimeters of tissue-like media was investigated. Fluorescence emission spectra measured from 6 cm homogeneous tissue-simulating phantoms show dramatic spectral distortion which results in emission peak shifts of up to 60 nm in wavelength. Measured spectral shapes are highly dependent on the photon path length and the scattered photon field in the NIR amplifies the wavelength-dependent absorption of the fluorescence spectra. Simulations of the peak propagation using diffusion modeling describe the experimental observations and confirm the path length dependence of fluorescence emission spectra. Spectral changes are largest for long path length measurements and thus will be most important in human tomography studies in the NIR. Spectrally resolved detection strategies are required to detect and interpret these effects which may otherwise produce erroneous intensity measurements. This observed phenomenon is analogous to beam hardening in x-ray tomography, which can lead to image artifacts without appropriate compensation. The peak shift toward longer wavelengths, and therefore lower energy photons, observed for NIR luminescent signals propagating through tissue may readily be described as a beam softening phenomenon.

Tissue drug concentration determines whether fluorescence or absorption measurements are more sensitive in diffuse optical tomography of exogenous contrast agents

The measurement sensitivities of absorbing and fluorescing objects in tissue are compared to determine conditions for which fluorescence data are favorable over those derived from absorption. A simulated human breast volume was used to model the relative perturbation in boundary data caused by a deeply embedded anomaly containing elevated concentrations of theoretical exogenous contrast agents with absorption properties resembling lutetium texaphyrin (LuTex) and Indocyanine Green (ICG). Synthetic data were used to produce quantum yield values representing the transition between conditions favorable to fluorescence versus absorption imaging. The parameters explored include tumor-to-background contrast, background drug concentration, and excitation light filtering efficiency. Drug concentration in the background was the primary factor that determined which contrast mechanism provided the more sensitive measurements. Specifically, fluorescence measurements are favorable if background drug concentrations are below 135-200 nM for LuTex and 25-50 nM for ICG, while absorption measurements are more sensitive above these ranges.

In vivo fluorescence lifetime tomography

Local molecular and physiological processes can be imaged in vivo through perturbations in the fluorescence lifetime (FLT) of optical imaging agents. In addition to providing functional information, FLT methods can quantify specific molecular events and multiplex diagnostic and prognostic information. We have developed a fluorescence lifetime diffuse optical tomography (DOT) system for in vivo preclinical imaging. Data is captured using a time-resolved intensified charge coupled device (ICCD) system to measure fluorescence excitation and emission in the time domain. Data is then converted to the frequency domain, and we simultaneously reconstruct images of yield and lifetime using an extension to the normalized Born approach. By using differential phase measurements, we demonstrate DOT imaging of short lifetimes (from 350 ps) with high precision (+/-5 ps). Furthermore, this system retains the efficiency, speed, and flexibility of transmission geometry DOT. We demonstrate feasibility of FLT-DOT through a progressive series of experiments. Lifetime range and repeatability are first measured in phantoms. Imaging of subcutaneous implants then verifies the FLT-DOT approach in vivo in the presence of inhomogeneous optical properties. Use in a common research scenario is ultimately demonstrated by imaging accumulation of a targeted near-infrared (NIR) fluorescent-labeled peptide probe (cypate-RGD) in a mouse with a subcutaneous tumor.

Investigation of detection limits for diffuse optical tomography systems: I. Theory and experiment

We present a statistical test using simulated photon migration data and a noise model derived from the hardware of a particular diffuse optical tomography system to predict its detection limits. Our method allows us to assess the spatial distribution of the detection sensitivity of arbitrary geometries and noise without requiring phantom measurements and reconstructions. We determine the minimal detectable lesion size at selected lesion positions and compare the predicted results with phantom measurements carried out in a cup geometry.

Monte Carlo algorithm for efficient simulation of time-resolved fluorescence in layered turbid media

We present an efficient Monte Carlo algorithm for simulation of time-resolved fluorescence in a layered turbid medium. It is based on the propagation of excitation and fluorescence photon bundles and the assumption of equal reduced scattering coefficients at the excitation and emission wavelengths. In addition to distributions of times of arrival of fluorescence photons at the detector, 3-D spatial generation probabilities were calculated. The algorithm was validated by comparison with the analytical solution of the diffusion equation for time-resolved fluorescence from a homogeneous semi-infinite turbid medium. It was applied to a two-layered model mimicking intra- and extracerebral compartments of the adult human head.

Image-guided Raman spectroscopic recovery of canine cortical bone contrast in situ

Raman scattering provides valuable biochemical and molecular markers for studying bone tissue composition with use in predicting fracture risk in osteoporosis. Raman tomography can image through a few centimeters of tissue but is limited by low spatial resolution. X-ray computed tomography (CT) imaging can provide high-resolution image-guidance of the Raman spectroscopic characterization, which enhances the quantitative recovery of the Raman signals, and this technique provides additional information to standard imaging methods. This hypothesis was tested in data measured from Teflon tissue phantoms and from a canine limb. Image-guided Raman spectroscopy (IG-RS) of the canine limb using CT images of the tissue to guide the recovery recovered a contrast of 145:1 between the cortical bone and background. Considerably less contrast was found without the CT image to guide recovery. This study presents the first known IG-RS results from tissue and indicates that intrinsically high contrasts (on the order of a hundred fold) are available.

A time domain fluorescence tomography system for small animal imaging

We describe the application of a time domain diffuse fluorescence tomography system for whole body small animal imaging. The key features of the system are the use of point excitation in free space using ultrashort laser pulses and noncontact detection using a gated, intensified charge-coupled device (CCD) camera. Mouse shaped epoxy phantoms, with embedded fluorescent inclusions, were used to verify the performance of a recently developed asymptotic lifetime-based tomography algorithm. The asymptotic algorithm is based on a multiexponential analysis of the decay portion of the data. The multiexponential model is shown to enable the use of a global analysis approach for a robust recovery of the lifetime components present within the imaging medium. The surface boundaries of the imaging volume were acquired using a photogrammetric camera integrated with the imaging system, and implemented in a Monte-Carlo model of photon propagation in tissue. The tomography results show that the asymptotic approach is able to separate axially located fluorescent inclusions centered at depths of 4 and 10 mm from the surface of the mouse phantom. The fluorescent inclusions had distinct lifetimes of 0.5 and 0.95 ns. The inclusions were nearly overlapping along the measurement axis and shown to be not resolvable using continuous wave (CW) methods. These results suggest the practical feasibility and advantages of a time domain approach for whole body small animal fluorescence molecular imaging, particularly with the use of lifetime as a contrast mechanism.

Molecular imaging with optics: primer and case for near-infrared fluorescence techniques in personalized medicine

We compare and contrast the development of optical molecular imaging techniques with nuclear medicine with a didactic emphasis for initiating readers into the field of molecular imaging. The nuclear imaging techniques of gamma scintigraphy, single-photon emission computed tomography, and positron emission tomography are first briefly reviewed. The molecular optical imaging techniques of bioluminescence and fluorescence using gene reporter/probes and gene reporters are described prior to introducing the governing factors of autofluorescence and excitation light leakage. The use of dual-labeled, near-infrared excitable and radio-labeled agents are described with comparative measurements between planar fluorescence and nuclear molecular imaging. The concept of time-independent and -dependent measurements is described with emphasis on integrating time-dependent measurements made in the frequency domain for 3-D tomography. Finally, we comment on the challenges and progress for translating near-infrared (NIR) molecular imaging agents for personalized medicine.

Fluorescence diffuse tomography for detection of red fluorescent protein expressed tumors in small animals

A fluorescence diffuse tomography (FDT) setup for monitoring tumor growth in small animals has been created. In this setup an animal is scanned in the transilluminative configuration by a single source and detector pair. To remove stray light in the detection system, we used a combination of interferometric and absorption filters. To reduce the scanning time, an experimental animal was scanned using the following algorithm: (1) large-step scanning to obtain a general view of the animal (source and detector move synchronously); (2) selection of the fluorescing region; and (3) small-step scanning of the selected region and different relative shifts between the source and detector to obtain sufficient information for 3D reconstruction. We created a reconstruction algorithm based on the Holder norm to estimate the fluorophore distribution. This algorithm converges to the solution with a minimum number of fluorescing zones. The use of tumor cell lines transfected with fluorescent proteins allowed us to conduct intravital monitoring studies. Cell lines of human melanomas Mel-P, Mel-Ibr, Mel-Kor, and human embryonic kidney HEK293 Phoenix were transfected with DsRed-Express and Turbo-RFP genes. The emission of red fluorescent proteins (RFPs) in the long-wave optical range permits detection of deep-seated tumors. In vivo experiments were conducted immediately after subcutaneous injection of fluorescing cells into small animals.

Tutorial on diffuse light transport

A tutorial introduction to diffuse light transport is presented. The basic analytic equations of time-resolved, steady-state and modulated light transport are introduced. The perturbation method for handling slight heterogeneities in optical properties is outlined. The treatment of boundary conditions such as an air/tissue surface is described. Finite mesh-based numerical methods are introduced to calculate the diffuse light field in complex tissues with arbitrary boundaries. Applications in tissue spectroscopy and imaging illustrate these theoretical and computational tools.

Radiative transport-based frequency-domain fluorescence tomography

We report the development of radiative transport model-based fluorescence optical tomography from frequency-domain boundary measurements. The coupled radiative transport model for describing NIR fluorescence propagation in tissue is solved by a novel software based on the established Attila particle transport simulation platform. The proposed scheme enables the prediction of fluorescence measurements with non-contact sources and detectors at a minimal computational cost. An adjoint transport solution-based fluorescence tomography algorithm is implemented on dual grids to efficiently assemble the measurement sensitivity Jacobian matrix. Finally, we demonstrate fluorescence tomography on a realistic computational mouse model to locate nM to microM fluorophore concentration distributions in simulated mouse organs.

A penalized linear and nonlinear combined conjugate gradient method for the reconstruction of fluorescence molecular tomography

Conjugate gradient method is verified to be efficient for nonlinear optimization problems of large-dimension data. In this paper, a penalized linear and nonlinear combined conjugate gradient method for the reconstruction of fluorescence molecular tomography (FMT) is presented. The algorithm combines the linear conjugate gradient method and the nonlinear conjugate gradient method together based on a restart strategy, in order to take advantage of the two kinds of conjugate gradient methods and compensate for the disadvantages. A quadratic penalty method is adopted to gain a nonnegative constraint and reduce the illposedness of the problem. Simulation studies show that the presented algorithm is accurate, stable, and fast. It has a better performance than the conventional conjugate gradient-based reconstruction algorithms. It offers an effective approach to reconstruct fluorochrome information for FMT.

A fast reconstruction algorithm for fluorescence optical diffusion tomography based on preiteration

Fluorescence optical diffusion tomography in the near-infrared (NIR) bandwidth is considered to be one of the most promising ways for noninvasive molecular-based imaging. Many reconstructive approaches to it utilize iterative methods for data inversion. However, they are time-consuming and they are far from meeting the real-time imaging demands. In this work, a fast preiteration algorithm based on the generalized inverse matrix is proposed. This method needs only one step of matrix-vector multiplication online, by pushing the iteration process to be executed offline. In the preiteration process, the second-order iterative format is employed to exponentially accelerate the convergence. Simulations based on an analytical diffusion model show that the distribution of fluorescent yield can be well estimated by this algorithm and the reconstructed speed is remarkably increased.

Fluorescence diffuse optical tomography with functional and anatomicala prioriinformation: feasibility study

Fluorescence diffuse optical tomography (FT) is an emerging molecular imaging technique that can spatially resolve both fluorophore concentration and lifetime parameters. In this study, we investigate the performance of a frequency-domain FT system for small inclusions that are embedded in a heterogeneous background. The results demonstrate that functional and structural a priori information is crucial to be able to recover both parameters with high accuracy. The functional a priori information is defined by the absorption and scattering maps at both excitation and emission wavelengths. Similarly, the boundaries of the small inclusion and different regions in the background are utilized as the structural a priori information. Without a priori information, the fluorophore concentration of a 5 mm inclusion in a 40 mm medium is recovered with 50% error, while the lifetime cannot be recovered at all. On the other hand, when both functional and structural information are available, the true lifetime can be recovered and the fluorophore concentration can be estimated only with 5% error. This study shows that a hybrid system that can acquire diffuse optical absorption tomography (DOT), FT and anatomical images in the same setting is essential to be able to recover the fluorophore concentration and lifetime accurately in vivo.

Fluorescent molecular tomographic image reconstruction based on the Green's function

The rapid and accurate computation of the Jacobian matrix, which is usually computationally intensive, is of critical importance for the reconstruction problem of fluorescent molecular tomography (FMT). An extension of the Green's function method for the Jacobian matrix computation suitable for two coupled differential equations is proposed, in combination with the parallel forward computing strategy for FMT image reconstruction. For further acceleration of the reconstruction process without significant quality reduction of the results, we also propose to reconstruct the FMT image on an adaptively refined mesh generated with a priori information incorporated. Experimental results demonstrate that the speed of the reconstruction process can be significantly improved with the proposed overall algorithm.

Fluorescence-enhanced three-dimensional lifetime imaging: a phantom study

Near-infrared fluorescence optical imaging has the unique opportunity of differentiating diseased lesions from normal lesions based upon environmentally indicated changes in the lifetime of a fluorescent imaging agent. In this paper, we demonstrate three-dimensional lifetime tomography using the gradient-based penalty modified barrier function with simple bounds truncated Newton with trust region method to reconstruct lifetime maps in a clinically relevant, single breast-shaped ( approximately 1081 cm(3)) phantom from point-frequency-domain photon migration measurements at 100 MHz. A reverse differentiation technique is used to calculate the gradients. This algorithm is desirable because the storage benefit from the use of the truncated Newton method and the reverse differentiation technique increase the speed. Two fluorescent contrast agents, indocyanine green and 3-3'-diethylthiatricarbocyanine iodide which differed in their fluorescence lifetimes by 0.62 ns, were used. Images of targets at a depth of 2.0 cm and target-to-background ratios (T:B) of 212:1 and 70:1 in fluoroscence absorption and 1:2.1 and 2.1:1 in lifetimes are successfully reconstructed. Our results show that image reconstruction is possible when there is (i) a longer lifetime in a target than the background and (ii) a shorter lifetime in a target than the background.

Feasibility analysis of an epidermal glucose sensor based on time-resolved fluorescence

The goal of this study is to test the feasibility of using an embedded time-resolved fluorescence sensor for monitoring glucose concentration. Skin is modeled as a multilayer medium with each layer having its own optical properties and fluorophore absorption coefficients, lifetimes, and quantum yields obtained from the literature. It is assumed that the two main fluorophores contributing to the fluorescence at these excitation and emission wavelengths are nicotinamide adenine dinucleotide (NAD)H and collagen. The intensity distributions of excitation and fluorescent light in skin are determined by solving the transient radiative transfer equation by using the modified method of characteristics. The fluorophore lifetimes are then recovered from the simulated fluorescence decays and compared with the actual lifetimes used in the simulations. Furthermore, the effect of adding Poissonian noise to the simulated decays on recovering the lifetimes was studied. For all cases, it was found that the fluorescence lifetime of NADH could not be recovered because of its negligible contribution to the overall fluorescence signal. The other lifetimes could be recovered to within 1.3% of input values. Finally, the glucose concentrations within the skin were recovered to within 13.5% of their actual values, indicating a possibility of measuring glucose concentrations by using a time-resolved fluorescence sensor.

Subsurface diffuse optical tomography can localize absorber and fluorescent objects but recovered image sensitivity is nonlinear with depth

Subsurface tomography with diffuse light has been investigated with a noncontact approach to characterize the performance of absorption and fluorescence imaging. Using both simulations and experiments, the reconstruction of local subsurface heterogeneity is demonstrated, but the recovery of target size and fluorophore concentration is not linear when changes in depth occur, whereas the mean position of the object for experimental fluorescent and absorber targets is accurate to within 0.5 and 1.45 mm when located within the first 10 mm below the surface. Improvements in the linearity of the response with depth appear to remain challenging and may ultimately limit the approach to detection rather than characterization applications. However, increases in tissue curvature and/or the addition of prior information are expected to improve the linearity of the response. The potential for this type of imaging technique to serve as a surgical guide is highlighted.

A parallel reconstruction scheme in fluorescence tomography based on contrast of independent inversed absorption properties

Based on an independent forward model in fluorescent tomography, a parallel reconstructed scheme for inhomogeneous mediums with unknown absorption property is proposed in this paper. The method considers the two diffusion equations as separately describing the propagation of excited light in tissues with and without fluorescent probes inside. Then the concentration of fluorophores is obtained directly through the difference between two estimations of absorption coefficient which can be parallel inversed. In this way, the multiparameter estimation problem in fluorescent tomography is transformed into two independent single-coefficient determined schemes of diffusion optical tomography (DOT). Any algorithms proved to be efficient and effective in DOT can be directly applied here. In this study the absorption property is estimated from the independent diffusion equations by a gradient-based optimization method with finite element method (FEM) solving the forward model. Simulation results of three representative occasions show that the reconstructed method can well estimate fluorescent property and tissue absorption distribution.

Fluorescence molecular imaging

There is a wealth of new fluorescent reporter technologies for tagging of many cellular and subcellular processes in vivo. This imposed contrast is now captured with an increasing number of available imaging methods that offer new ways to visualize and quantify fluorescent markers distributed in tissues. This is an evolving field of imaging sciences that has already achieved major advances but is also facing important challenges. It is nevertheless well poised to significantly impact the ways of biological research, drug discovery, and clinical practice in the years to come. Herein, the most pertinent technologies associated with in vivo noninvasive or minimally invasive fluorescence imaging of tissues are summarized. Focus is given to small-animal imaging. However, while a broad spectrum of fluorescence reporter technologies and imaging methods are outlined, as necessary for biomedical research, and clinical translation as well.