Spatially resolved detection of fluorescence-labeled agent distribution in mice has been demonstrated to hold great potential in a broad spectrum of preclinical molecular imaging studies.1 Because of the mathematical difficulties in resolving the highly scattering nature of light photon propagation in tissue, the majority of in vivo small animal fluorescence imaging applications in the laboratory settings is mostly performed to generate two-dimensional (2-D) planar images following epi-illumination.2 Nonetheless, to improve localization as well as to enable quantification of fluorescent sources, efforts have been taken toward fully three-dimensional (3-D) tomographic reconstruction of fluorescence light distribution within heterogeneously turbid and arbitrarily shaped objects such as mice and, in some cases, up to the size of rats.3,4
Noncontact data acquisition approaches have replaced fiber-based contact imaging implementations for the most part mainly because of animal handling and efficiency shortcomings of the latter.5 Some noncontact imaging systems perform fluorescence optical tomography (FOT) reconstruction by slightly squeezing the imaged animal into, or between, a purposely shaped chamber of simple geometry (e.g., translucent vices),6,7 or within translucent compartments of (generally) cylindrical geometry with lesser spacing constraints forced onto the animal.8,9 Planar spacing by the use of opposite vices yields a much reduced boundary complexity involved in solving the inverse problem. However, this method substantially restricts freedom of excitation and detection geometry. Cylindrical animal burrow enclosures maintain a “normal” body shape and allow for arbitrary excitation and detection from the full set of 360 deg source angles and imaging projections.8,9 However, this geometry involves a rather complex problem space for the inverse problem. The imaging system embodied herein uses the latter approach incorporating an animal enclosure with a bore holder of 30-mm inner diameter. Fluorescence surface imaging (FSI) is performed by illuminating the entire imaged object, or local regions thereof, with a large, preferably homogeneous field of light to excite fluorochromes in a large area, generally the whole animal.
In almost all noncontact FOT approaches, though, laser beam generated point illumination is used to sequentially form an illumination pattern on the surface of the animal to create a sufficiently high number of source points on the surface, which are paired with some data projection schema, respectively.10 Such scanning strategy, however, can be time consuming. Adoption of galvanometer mirrors as suggested by Patwardhan et al. can accelerate the scanning process.11 This strategy, however, still induces a large amount of acquired data and, hence, yields a large space problem for the FOT image reconstruction. Following recent studies, wide field illumination using structured light or spatially modulated light has been extensively investigated for turbid media as well as for in vivo imaging.1213.14.15.–16 Those interesting approaches are intended for lessening the inherently ill-posed nature of FOT1213.–14 or lowering the number of excitation points.15,16
In the present study, we seek to investigate the use of single excitation line illumination as generated by a line laser source and conduct comparative experiments for FSI and FOT. The line source is fit and aligned with the axial center line of the imaging bore and rotationally mounted at some distance from the imaged object. Hence, a single 360-deg revolution (at certain angular stepping offsets) of the line source forms an illumination pattern along with a set of projection data that could be used for both FSI and FOT.
Combining in vivo optical imaging with further imaging modalities such as CT is being adopted more and more to improve reconstruction accuracy and, hence, signal localization.1718.–19 To conduct this study, the optical imaging setup, including all light sources, is for that reason placed into a trimodal small animal imaging system, incorporating also CT and SPECT.20 This enables simultaneous, coregistered FSI and FOT imaging with CT and SPECT.
Materials and Methods
In the given trimodal embodiment, cf. Fig. 1(a) and 1(b), the optical imaging part is composed of a light camera (ORCA-AG, Hamamatsu Photonics, Shizuoka, Japan) and two fitted laser diode sources for probe excitation: (1) a collimated point-laser with 20 mW power at 660 nm (FP-66/20AAF-AV, Laser Components GmbH, Olching, Germany) and (2) a collimated line-laser with 70 mW power at 660 nm (FP-66/70LF-Ö90-HOM, Laser Components GmbH, Olching, Germany). The point laser diode is aligned and can be shifted parallel to the long axis of the imaged object. Hence, any point radial along the long axis can be illuminated [Fig. 1(c) left]. The line laser beam is aligned with the axial long axis of the bore, as well. Its light field has a minimum length of 100 mm within a radius of 25 mm, hence illuminating a whole mouse [Fig. 1(c) right]. For all data acquisition studies conducted herein, a filter set combining a long-pass filter (RG695, Schott AG, Mainz, Germany) and a bandpass filter (FF01-716/40-25, Semrock Inc., Rochester, NY) is attached to the camera, and two bandpass filters (FF01-655/40-25, Semrock Inc., Rochester, NY) are fixed to the two laser sources, therefore enabling the acquisition of probes such as Cy5.5, DyLight 680, or Alexa Fluor 680.
The x-ray CT imaging part is composed of an x-ray tube having a 47.8 μm beam focus (Apogee Series 500, Oxford Instruments, Oxfordshire, UK) that is aligned to a flat panel x-ray detector (Shad-o-Box 2048, Rad-icon, Sunnyvale, CA) featuring a pixel array with 48 μm pixel size. Three-dimensional images are reconstructed employing the Feldkamp–Davis–Kress (FDK) algorithm.21
To accomplish SPECT imaging a self-build gamma camera is used incorporating pixelated NaI(Tl) crystals of in size, arranged on a array with 1.5 mm pixel spacing, the crystal array being fitted to a multianode PMT array (Hamamatsu H8500, Hamamatsu Photonics, Shizuoka, Japan). A slit-slat collimator as described in Ref. 22 is attached to the gamma camera enabling whole-body mouse imaging [axial field-of-view ]. Three-dimensional images are reconstructed employing the maximum-likelihood expectation-maximization algorithm.23
All imaging components are mounted on a common rotatable and translatable gantry, and share fully overlapping fields of view [Fig. 1(b)]. All subsystems are fully calibrated yielding intrinsically fused intermodal and FSI-FOT data.24 Careful calibration is particularly essential not only for fused imaging representation and derived diagnostic interpretation, but constitutes a crucial precondition for both FSI mapping and FOT reconstruction algorithms as laid out in the following.
3-D Fluorescence Surface Imaging
Figure 2 shows the flow diagram of the algorithm developed and used to obtain FSI data at each surface point of the 3-D object. In short, reconstructed x-ray CT data of the imaged object are segmented first by regional growing to provide the object’s boundary shape. Detected fluorescence signals are then mapped onto the surface for every imaging angle over a 360-deg rotation according to the specific intermodal system geometry. Finally, cumulated FSI data are normalized with respect to the (possibly varying) number of contributing light rays as illustrated in Fig. 2. This strategy is basically the same for either illumination pattern (point or line).
In order to visualize FSI data in the form of projectional views—which are similar and comparable to views from 2-D epi-illumination acquisitions—2-D projection images are digitally generated using the fused 3-D CT (printed in grayscale) and FSI (printed in pseudocolor) images at any arbitrarily chosen view angle (figures in Secs. 3.1.1 and 3.2.1 for the visualization results).
Fluorescence Optical Tomography
If point-laser source driven FOT is performed, then the 3-D image reconstruction algorithm as proposed and presented in Ref. 19 is being used. Three-dimensional tetrahedral meshes of the imaged object are generated from the boundary segmented CT data using the TetGen library.25
If, however, the line source illumination pattern is used for probe excitation, then a modified 2-D image reconstruction algorithm is used as described below to solve the inverse problem at any transversal plane of interest. Assuming that the line source excitation is homogeneous along its entire length, a stack of 2-D image reconstructions yields a volumetric (3-D) map of estimated fluorochrome distribution. Two-dimensional triangular meshes, generated at any transversal slice of interest, serve as the object’s geometric representation. The 2-D image reconstruction algorithm has been derived from the 3-D algorithm as given in Ref. 19 whereby the propagation of photons inside the imaged object follows the solution of coupled 2-D diffusion equations at any transversal slice as illuminated by the line source which is expressed as
Assuming that is the optical measurements for detection positions and refers to the estimated fluorescence distribution at mesh elements, the forward system model can be formulated as , in which is the system matrix calculated based on the numerical solution to Eq. (1) by the finite element method (FEM). The geometry used in FEM is based on 2-D meshes extracted from the corresponding CT data. Employing Tikhonov regularization, the corresponding reconstruction problem is formulated as
To validate the proposed method, a mouse-shaped phantom was adopted (XFM-2 Fluorescent Phantom, Caliper Life Science Inc., Hopkinton, MA; and ). Two point sources filled with water solution of Cy5.5 Mono NHS Ester (approximately 40 pmol each; GE Healthcare, Buckinghamshire, UK) at the tips of two cylindrical rods were inserted into the two holes of the phantom [Fig. 3(a)]. Figure 3(b) shows the exact positions of the two point sources in the phantom. CT data acquisition scans are first performed at 40 kV anode voltage, 0.4 mA tube current, 1 s acquisition time per projection, 240 projections per 360-deg rotation, whereby two axially shifted rotations are performed to cover an axial FOV of 100 mm. Following CT data acquisition, fluorescence data (60 projections per 360-deg rotation, 0.5 s exposure time per projection) were collected firstly employing point-laser excitation at seven axial positions (with axial intervals of 5 mm) with seven rotations, and thereafter employing line-laser excitation with a single rotation. For comparison, this phantom was also acquired by a self-built epi-illumination 2-D fluorescence imaging system which incorporates a CCD camera (ORCAII-BT-512G, Hamamatsu Photonics K.K., Hamamatsu City, Japan) that is attached with an emission bandpass filter (FF01-716/40-25, Semrock Inc., Rochester, NY) as well as 4 LEDs (LR W5SM, OSRAM Opto Semiconductors GmbH, Regensburg, Germany) equipped with appropriate excitation filters (FF01-655/40-25, Semrock Inc., Rochester, NY).
In Vivo Mouse Experiment
To further demonstrate the preclinical imaging ability of the method, an in vivo mouse study was performed. This study was carried out according to the policies and principles established by the German animal protection laws of district government Karlsruhe with the animal number G-61/10. A healthy nude mouse was fed with chlorophyll free diets (CRD FLI, Harlan Laboratories, Inc., Indiana, IN) for three consecutive days prior to imaging to reduce autofluorescence from abdomen and skin. One day prior to imaging, the mouse was injected intravenously with 150 μl RediJect Bone Probe 680 (Caliper Life Science, Inc., Hopkinton, MA, concentration: ). For cross-validation, SPECT data (120 projections in a 360-deg rotation, 30 s per projection) have been collected following the administration of 104 MBq tagged to methylene diphosphonate molecules (MDP), injected 50 min prior to imaging. CT imaging was performed with the same protocol as described in the phantom experiment. Following SPECT and CT scans, fluorescence data were acquired (60 projections per 360-deg rotation, 0.5 s exposure time per projection), first employing point-laser excitation and, thereafter, employing line-laser excitation. Image reconstruction was performed in reference to every modality as previously described.
Phantom Experiment Results
Figure 3(c) shows results of 2-D epi-illumination acquisition of the mouse phantom at supine and prone positions, as well as digital 2-D projection views of fused 3-D CT and FSI data using both line-laser and point-laser illumination patterns, respectively. The view angles of the fused CT-FSI data are selected to be identical to those for the 2-D acquisition with epi-illumination. As can be seen in Fig. 3(c), when 3-D FSI data are generated using point-laser illumination with only one axial position, the axial FOV is primarily covering one fluorescent point source. Using multiple point-laser or line-laser illumination, however, the whole phantom is exposed with one full rotation. The fluorescence signal intensities are normalized individually for each illumination strategy toward its maximum signal intensity. Figure 3(d) and 3(e) shows profiles of the signal intensities at the indicated positions providing quantitative comparison among different illumination types.
Figure 4(a) depicts part of the 3-D mesh as calculated from the acquired CT data. This data set includes 18,987 nodes and 73,339 tetrahedral elements. Figure 4(b) shows one of the derived 2-D meshes which includes 1017 nodes and 2593 triangular elements. Figure 5 shows results of 3-D FOT image reconstruction following point-laser excitation, alongside results of FOT 2-D image reconstruction following line-laser illumination; CT images are also provided as reference.
In Vivo Mouse Experiment Results
The visualization of FSI results of the in vivo mouse experiment using digital projection data is shown in Fig. 6(a). The upper row demonstrates the results following point-laser illumination at a single axial position, while the lower row holds the results following line-laser illumination. Profiles of the normalized fluorescence signal intensities are shown in Fig. 6(b) for a more quantitative comparison between the results following both illumination patterns.
With regard to FOT image reconstruction, the optical properties of mouse tissue were assumed to be homogeneous with and . In order to fuse the FOT data with the SPECT and CT results for diagnostic evaluation, the derived mesh geometries were reallocated into voxel-based representations at corresponding voxel resolutions. Results of 2-D FOT combining all slices are shown in Fig. 7 in the form of a maximum intensity projection view which is fused with SPECT–CT (left). CT, SPECT, 2-D FOT, and 3-D FOT data are also shown as transversal slices at the indicated axial positions (right).
Mapping the fluorescence signals onto the surface and normalization over a 360-deg rotation (FSI mode) provides the opportunity to view the imaging results similar to a planar acquisition with epi-illumination, however, at any arbitrary view angle [see Figs. 3(c) and 6]. It has been shown that changing the view angle in 2-D optical data acquisition can affect the measured signal intensity. Three-dimensional FSI using line-laser illumination is, therefore, advantageous to 2-D epi-illumination acquisition by providing digital projection data at any required angular view.
Whole-body fluorescence imaging, especially FOT, is challenging because the probe often does not form unique spots but distributes widely throughout the body, particularly when systemic assessment is required, e.g., in a skeleton scan. Employing line-laser excitation enables one to perform whole-body fluorescence imaging at a much shorter scan time as compared to data acquisition strategies employing a point-laser source. Moreover, because of the extended excitation pattern throughout the axial direction generated by the line laser, the approximation of Eq. (2) could be locally held in practice. It is, however, not the case for point-laser illumination. Two-dimensional image reconstruction of tomographic slices becomes, therefore, practically available over the entire animal following line-laser excitation. In comparison to 3-D FOT reconstruction with point-laser excitation, this method allows a much larger FOV given the same acquisition time. In the mouse phantom study with two point-like fluorescence inclusions, the 2-D FOT result from the line-laser excitation shows broader distributions of the resolved point sources, especially in the axial direction, as compared to the 3-D reconstruction result from the point-laser excitation (cf. Fig. 5). Although, in this phantom study, 3-D FOT with point-laser outperforms 2-D FOT with line-laser excitation. Note, however, that the data acquisition time of the latter case was seven times faster. The in vivo experiment performed in this study is in fact an extremely challenging example for FOT as the hydroxyapatite-targeting probe has an affinity to the entire skeleton system in the animal which is very different to the case of spatially distinct point sources as are being used in the mouse phantom. The 2-D FOT reconstruction is applied slice by slice along the whole axial FOV. Intrinsically coregistered SPECT results with , a probe with high bone-specificity as well, provides additional information for cross-validation. Two-dimensional FOT results are shown exemplarily in three different transaxial slices (cf. Fig. 7). It can be seen from “Slice 1” that the small bone structures close to surface are well resolved, while the spine, which shows a high activity uptake in SPECT, is not visible in FOT. In “Slice 2,” the 2-D FOT image shows not only the spine structure near to the surface but also probe uptake in the kidneys, while no tracer activity can be seen in SPECT. “Slice 3” includes spine and knee joints, showing similar information in FOT and SPECT. Employing point-laser excitation, however, only a small portion of the animal is within the active FOV. The reconstruction results show a more pronounced convergence at the surface area (see the last row of Fig. 7).
This study has been particularly conducted to lay out a frame for future investigations. While exemplary studies have been provided illustratively for each illumination pattern, the results are not conclusive with respect to the quantitative accuracy of either of the involved method.
In conclusion, a trimodal SPECT–CT–OT system was investigated in this study, specifically with the aim of investigating the use of a line-laser as the excitation source for 3-D fluorescence imaging and tomography. Phantom and in vivo experiments demonstrate the feasibility of line-laser excitation yielding an axially enlarged FOV in both 3-D FSI and FOT with short acquisition time.
The authors would like to thank Bartlomiej Grychtol for the support of mesh generation and Karin Leotta for excellent technical assistance during the animal experiment. Disclosure statement: The authors have no conflict of interest to declare.
S. A. HilderbrandR. Weissleder, “Near-infrared fluorescence: application to in vivo molecular imaging,” Curr. Opin. Chem. Biol. 14(1), 71–79 (2010).COCBF41367-5931http://dx.doi.org/10.1016/j.cbpa.2009.09.029Google Scholar
F. Leblondet al., “Pre-clinical whole-body fluorescence imaging: review of instruments, methods and applications,” J. Photochem. Photobiol. B 98(1), 77–94 (2010).JPPBEG1011-1344http://dx.doi.org/10.1016/j.jphotobiol.2009.11.007Google Scholar
A. CongG. Wang, “A finite-element-based reconstruction method for 3D fluorescence tomography,” Opt. Express 13(24), 9847–9857 (2005).OPEXFF1094-4087http://dx.doi.org/10.1364/OPEX.13.009847Google Scholar
H. Meyeret al., “Noncontact optical imaging in mice with full angular coverage and automatic surface extraction,” Appl. Opt. 46(17), 3617–3627 (2007).APOPAI0003-6935http://dx.doi.org/10.1364/AO.46.003617Google Scholar
P. Mohajeraniet al., “Compensation of optical heterogeneity-induced artifacts in fluorescence molecular tomography: theory and in vivo validation,” J. Biomed. Opt. 14(3), 034021 (2009).JBOPFO1083-3668http://dx.doi.org/10.1117/1.3149855Google Scholar
K. O. VasquezC. CasavantJ. D. Peterson, “Quantitative whole body biodistribution of fluorescent-labeled agents by non-invasive tomographic imaging,” PLoS One 6(6), e20594 (2011).1932-6203http://dx.doi.org/10.1371/journal.pone.0020594Google Scholar
N. Deliolaniset al, “Free-space fluorescence molecular tomography utilizing 360 degrees geometry projections,” Opt. Lett. 32(4), 382–384 (2007).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.32.000382Google Scholar
X. Songet al., “Reconstruction for free-space fluorescence tomography using a novel hybrid adaptive finite element algorithm,” Opt. Express 15(26), 18300–18317 (2007).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.15.018300Google Scholar
R. B. SchulzJ. RipollV. Ntziachristos, “Experimental fluorescence tomography of tissues with noncontact measurements,” IEEE Trans. Med. Imaging 23(4), 492–500 (2004).ITMID40278-0062http://dx.doi.org/10.1109/TMI.2004.825633Google Scholar
S. Patwardhanet al., “Time-dependent whole-body fluorescence tomography of probe bio-distributions in mice,” Opt. Express 13(7), 2564–2577 (2005).OPEXFF1094-4087http://dx.doi.org/10.1364/OPEX.13.002564Google Scholar
J. Duttaet al., “Illumination pattern optimization for fluorescence tomography: theory and simulation studies,” Phys. Med. Biol. 55(10), 2961–2982 (2010).PHMBA70031-9155http://dx.doi.org/10.1088/0031-9155/55/10/011Google Scholar
A. Bassiet al., “Spatial shift of spatially modulated light projected on turbid media,” J. Opt. Soc. Am. A 25(11), 2833–2839 (2008).JOAOD60740-3232http://dx.doi.org/10.1364/JOSAA.25.002833Google Scholar
C. D’Andreaet al., “Fast 3D optical reconstruction in turbid media using spatially modulated light,” Biomed. Opt. Express 1(2), 471–481 (2010).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.1.000471Google Scholar
N. Ducroset al., “Fluorescence molecular tomography of an animal model using structured light rotating view acquisition,” J. Biomed. Opt. 18(2), 020503 (2013).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.18.2.020503Google Scholar
C. Liet al., “Three-dimensional fluorescence optical tomography in small-animal imaging using simultaneous positron-emission-tomography priors,” Opt. Lett. 34(19), 2933–2935 (2009).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.34.002933Google Scholar
R. B. Schulzet al., “Hybrid system for simultaneous fluorescence and x-ray computed tomography,” IEEE Trans. Med. Imaging 29(2), 465–473 (2010).ITMID40278-0062http://dx.doi.org/10.1109/TMI.2009.2035310Google Scholar
L. CaoJ. Peter, “Bayesian reconstruction strategy of fluorescence-mediated tomography using an integrated SPECT–CT–OT system,” Phys. Med. Biol. 55(9), 2693–2708 (2010).PHMBA70031-9155http://dx.doi.org/10.1088/0031-9155/55/9/018Google Scholar
J. Peter, “Triple-modality imaging system,” U.S. Patent 2010/0030069 A1 (February 4 2010).Google Scholar
L. CaoJ. Peter, “Slit-slat collimator equipped gamma camera for whole-mouse SPECT–CT imaging,” IEEE Trans. Nucl. Sci. 59(3), 530–536 (2012).IETNAE0018-9499http://dx.doi.org/10.1109/TNS.2012.2192289Google Scholar
L. CaoM. BreithauptJ. Peter, “Geometrical co-calibration of a tomographic optical system with CT for intrinsically co-registered imaging,” Phys. Med. Biol. 55(6), 1591–1606 (2010).PHMBA70031-9155http://dx.doi.org/10.1088/0031-9155/55/6/004Google Scholar
H. SiK. Gärtner, “Meshing piecewise linear complexes by constrained Delaunay tetrahedralizations,” in Proc. 14th Int. Meshing Roundtable, pp. 147–163, Springer (2005).Google Scholar