An eNAMPT-neutralizing mAb reduces post-infarct myocardial fibrosis and left ventricular dysfunction

Myocardial infarction (MI) triggers adverse ventricular remodeling (VR), cardiac fibrosis, and subsequent heart failure. Extracellular nicotinamide phosphoribosyltransferase (eNAMPT) is postulated to play a significant role in VR processing via activation of the TLR4 inflammatory pathway. We hypothesized that an eNAMPT specific monoclonal antibody (mAb) could target and neutralize overexpressed eNAMPT post-MI and attenuate chronic cardiac inflammation and fibrosis. We investigated humanized ALT-100 and ALT-300 mAb with high eNAMPT-neutralizing capacity in an infarct rat model to test our hypothesis. ALT-300 was 99mTc-labeled to generate 99mTc-ALT-300 for imaging myocardial eNAMPT expression at 2 hours, 1 week, and 4 weeks post-IRI. The eNAMPT-neutralizing ALT-100 mAb (0.4 mg/kg) or saline was administered intraperitoneally at 1 hour and 24 hours post-reperfusion and twice a week for 4 weeks. Cardiac function changes were determined by echocardiography at 3 days and 4 weeks post-IRI. 99mTc-ALT-300 uptake was initially localized to the ischemic area at risk (IAR) of the left ventricle (LV) and subsequently extended to adjacent non-ischemic areas 2 hours to 4 weeks post-IRI. Radioactive uptake (%ID/g) of 99mTc-ALT-300 in the IAR increased from 1 week to 4 weeks (0.54 ± 0.16 vs. 0.78 ± 0.13, P < 0.01). Rats receiving ALT-100 mAb exhibited significantly improved myocardial histopathology and cardiac function at 4 weeks, with a significant reduction in the collagen volume fraction (%LV) compared to controls (21.5 ± 6.1% vs. 29.5 ± 9.9%, P < 0.05). Neutralization of the eNAMPT/TLR4 inflammatory cascade is a promising therapeutic strategy for MI by reducing chronic inflammation, fibrosis, and preserving cardiac function.

Resolution recovery on list mode MLEM reconstruction for dynamic cardiac SPECT system

The Dynamic Cardiac SPECT (DC-SPECT) system is being developed at the Massachusetts General Hospital, featuring a static cardio focus asymmetrical geometry enabling simultaneous high-resolution and high-sensitivity imaging. Among 14 design iterations of the DC-SPECT with varying number of detector heads, system sensitivity and resolution, the current version under development features 10 mm FWHM geometrical resolution (without resolution recovery) and 0.07% sensitivity at the center of the FOV, this is 1.5× resolution gain and 7× sensitivity gain compared to a conventional dual head gamma camera (0.01% sensitivity and 15-mm resolution). This work presents improvement in imaging resolution by implementing a spatially variant point spread function (SV-PSF) with list mode MLEM reconstruction. A resolution recovery method by PSF deconvolution is validated on list mode MLEM reconstruction for the DC-SPECT. A spatial invariant PSF is included as an additional test to show the influence of the PSF modelling accuracy on reconstructed image quality. We compare the MLEM reconstruction with and without PSF deconvolution; an analytic model is used for the calculation of system response, and the results are compared to the reconstruction with system modelling using Monte Carlo (MC) based methods. Results show that with PSF modelling applied, the quality of the reconstructed image is improved, and the DC-SPECT system can achieve a 4.5 mm central spatial resolution with average 795 counts/Mbq. Both the SV-PSF and the spatial-invariant PSF improve the image quality, and the reconstruction with SV-PSF generates line profiles closer to the ground truth. The results show substantial improvement over the GE Discovery 570c performance (7 mm spatial resolution with an average 460 counts/MBq, 5.8 mm resolution at the FOV center). The impact of PSF deconvolution is significant, improvement of the reconstructed image quality is evident in comparison to MC simulated system matrix with the same sampling size in the simulation.

Yttrium-90 Activity Quantification in PET/CT-Guided Biopsy Specimens from Colorectal Hepatic Metastases Immediately after Transarterial Radioembolization Using Micro-CT and Autoradiography

PURPOSE

To evaluate the yttrium-90 (90Y) activity distribution in biopsy tissue samples of the treated liver to quantify the dose with higher spatial resolution than positron emission tomography (PET) for accurate investigation of correlations with microscopic biological effects and to evaluate the radiation safety of this procedure.

MATERIALS AND METHODS

Eighty-six core biopsy specimens were obtained from 18 colorectal liver metastases (CLMs) immediately after 90Y transarterial radioembolization (TARE) with either resin or glass microspheres using real-time 90Y PET/CT guidance in 17 patients. A high-resolution micro-computed tomography (micro-CT) scanner was used to image the microspheres in part of the specimens and allow quantification of 90Y activity directly or by calibrating autoradiography (ARG) images. The mean doses to the specimens were derived from the measured specimens' activity concentrations and from the PET/CT scan at the location of the biopsy needle tip for all cases. Staff exposures were monitored.

RESULTS

The mean measured 90Y activity concentration in the CLM specimens at time of infusion was 2.4 ± 4.0 MBq/mL. The biopsies revealed higher activity heterogeneity than PET. Radiation exposure to the interventional radiologists during post-TARE biopsy procedures was minimal.

CONCLUSIONS

Counting the microspheres and measuring the activity in biopsy specimens obtained after TARE are safe and feasible and can be used to determine the administered activity and its distribution in the treated and biopsied liver tissue with high spatial resolution. Complementing 90Y PET/CT imaging with this approach promises to yield more accurate direct correlation of histopathological changes and absorbed dose in the examined specimens.

Mesh modeling of system geometry and anatomy phantoms for realistic GATE simulations and their inclusion in SPECT reconstruction

Objective.Monte-Carlo simulation studies have been essential for advancing various developments in single photon emission computed tomography (SPECT) imaging, such as system design and accurate image reconstruction. Among the simulation software available, Geant4 application for tomographic emission (GATE) is one of the most used simulation toolkits in nuclear medicine, which allows building systems and attenuation phantom geometries based on the combination of idealized volumes. However, these idealized volumes are inadequate for modeling free-form shape components of such geometries. Recent GATE versions alleviate these major limitations by allowing users to import triangulated surface meshes.Approach.In this study, we describe our mesh-based simulations of a next-generation multi-pinhole SPECT system dedicated to clinical brain imaging, called AdaptiSPECT-C. To simulate realistic imaging data, we incorporated in our simulation the XCAT phantom, which provides an advanced anatomical description of the human body. An additional challenge with the AdaptiSPECT-C geometry is that the default voxelized XCAT attenuation phantom was not usable in our simulation due to intersection of objects of dissimilar materials caused by overlap of the air containing regions of the XCAT beyond the surface of the phantom and the components of the imaging system.Main results.We validated our mesh-based modeling against the one constructed by idealized volumes for a simplified single vertex configuration of AdaptiSPECT-C through simulated projection data of123I-activity distributions. We resolved the overlap conflict by creating and incorporating a mesh-based attenuation phantom following a volume hierarchy. We then evaluated our reconstructions with attenuation and scatter correction for projections obtained from simulation consisting of mesh-based modeling of the system and the attenuation phantom for brain imaging. Our approach demonstrated similar performance as the reference scheme simulated in air for uniform and clinical-like123I-IMP brain perfusion source distributions.Significance.This work enables the simulation of complex SPECT acquisitions and reconstructions for emulating realistic imaging data close to those of actual patients.

PEGylated and Non-PEGylated TCP-1 Probes for Imaging of Colorectal Cancer

PURPOSE

Previous studies indicate that 99mTc- and fluorescent-labeled c[Cys-Thr-Pro-Ser-Pro-Phe-Ser-His-Cys]OH (TCP-1) peptides were able to detect colorectal cancer (CRC) and tumor-associated vasculature. This study was designed to characterize the targeting properties of PEGylated and non-PEGylated TCP-1 peptides for CRC imaging.

PROCEDURES

Cell uptake of cyanine 7 (Cy7)-labeled TCP-1 probes (Cy7-PEG4-TCP-1 and Cy7-TCP-1) was investigated in three CRC cell lines (human, HCT116 and HT29; mouse, CT26). Xenograft and orthotopic CRC tumor models with HCT116 and CT26 cells were used to characterize biodistribution and CRC tumor-targeting properties of TCP-1 fluorescence and radioligand with and without PEGylation, [99mTc]Tc-HYNIC-PEG4-TCP-1 vs. [99mTc]Tc-HYNIC-TCP-1.

RESULTS

Fluorescence images showed that TCP-1 probes were distributed in the cytoplasm and nucleus of CRC cells. When CT26 cells were treated with unlabeled TCP-1 peptide prior to the cell incubation with Cy7-PEG4-TCP-1, cell fluorescent signals were significantly reduced relative to the cells without blockade. Relative to Cy7-TCP-1, superior brilliance and visibility of fluorescence was observed in the tumor with Cy7-PEG4-TCP-1 and maintained up to 18 h post-injection. [99mTc]Tc-HYNIC-PEG4-TCP-1 images in xenograft and orthotopic CRC models demonstrated that TCP-1 PEGylation preserved tumor-targeting capability of TCP-1, but its distribution (%ID/g) in the liver and intestine was higher than that of [99mTc]Tc-HYNIC-TCP-1 (1.51 ± 0.29 vs 0.53 ± 0.12, P < 0.01). Better tumor visualization by [99mTc]Tc-HYNIC-TCP-1 was observed in the orthotopic CRC model due to lower intestinal radioactivity.

CONCLUSIONS

TCP-1-based probes undergo endocytosis and localize in the cytoplasm and nucleus of human and mouse CRC cells. Tumor detectability of fluorescent TCP-1 peptide with a PEG4 spacer is promising due to its enhanced tumor binding affinity and rapid clearance kinetics from nontumor tissues. Non-PEGylated [99mTc]Tc-HYNIC-TCP-1 exhibits lower nonspecific accumulation in the liver and gastrointestinal tract and might have better capability for detecting CRC lesions in clinical sites. TCP-1 may represent an innovative targeting molecule for detecting CRC noninvasively.

Use of radiolabeled hyaluronic acid for preclinical assessment of inflammatory injury and acute respiratory distress syndrome

Acute respiratory distress syndrome (ARDS) is accompanied by a dramatic increase in lung hyaluronic acid (HA), leading to a dose-dependent reduction of pulmonary oxygenation. This pattern is associated with severe infections, such as COVID-19, and other important lung injury etiologies. HA actively participates in molecular pathways involved in the cytokine storm of COVID-19-induced ARDS. The objective of this study was to evaluate an imaging approach of radiolabeled HA for assessment of dysregulated HA deposition in mouse models with skin inflammation and lipopolysaccharide (LPS)-induced ARDS using a novel portable intensified Quantum Imaging Detector (iQID) gamma camera system.

METHODS

HA of 10 kDa molecular weight (HA10) was radiolabeled with 125I and 99mTc respectively to produce [125I]I-HA10 and [99mTc]Tc-HA10, followed by comparative studies on stability, in vivo biodistribution, and uptake at inflammatory skin sites in mice with 12-O-tetradecanoylphorbol-13-acetate (TPA)-inflamed ears. [99mTc]Tc-HA10 was used for iQID in vivo dynamic imaging of mice with ARDS induced by intratracheal instillation of LPS.

RESULTS

[99mTc]Tc-HA10 and [125I]I-HA10 had similar biodistribution and localization at inflammatory sites. [99mTc]Tc-HA10 was shown to be feasible in measuring skin injury and monitoring skin wound healing. [99mTc]Tc-HA10 dynamic pulmonary images yielded good visualization of radioactive uptake in the lungs. There was significantly increased lung uptake and slower lung washout in mice with LPS-induced ARDS than in control mice. Postmortem biodistribution measurement of [99mTc]TcHA10 (%ID/g) was 11.0 ± 3.9 vs. 1.3 ± 0.3 in the ARDS mice (n = 6) and controls (n = 6) (P < 0.001), consistent with upregulated HA expression as determined by enzyme-linked immunosorbent assay (ELISA) and immunohistochemistry (IHC) staining.

CONCLUSIONS

[99mTc]Tc-HA10 is promising as a biomarker for evaluating HA dysregulation that contributes to pulmonary injury in ARDS. Rapid iQID imaging of [99mTc]Tc-HA10 clearance from injured lungs may provide a functional template for timely assessment and quantitative monitoring of pulmonary pathophysiology and intervention in ARDS.

Improved Performance of a Multipinhole SPECT for DAT Imaging by Increasing Number of Pinholes at the Expense of Increased Multiplexing

SPECT imaging of dopamine transporters (DAT) in the brain is a widely utilized study to improve the diagnosis of Parkinsonian syndromes, where conventional (parallel-hole and fan-beam) collimators on dual-head scanners are commonly employed. We have designed a multi-pinhole (MPH) collimator to improve the performance of DAT imaging. The MPH collimator focuses on the striatum and hence offers a better trade-off for sensitivity and spatial resolution than the conventional collimators within this clinically most relevant region for DAT imaging. Our original MPH design consisted of 9 pinholes with a background-to-striatal (Bkg/Str) projection multiplexing of 1% only. In this simulation study, we investigated whether further improvements in the performance of MPH imaging could be obtained by increasing the number of pinholes, hence by enhancing the sensitivity and sampling, despite the ambiguity in reconstructing images due to increased multiplexing. We performed analytic simulations of the MPH configurations with 9, 13, and 16 pinholes (aperture diameters: 4-6mm) using a digital phantom modeling DAT imaging. Our quantitative analyses indicated that using 13 (Bkg/Str: 12%) and 16 (Bkg/Str: 22%) pinholes provided better performance than the original 9-pinhole configuration for the acquisition with 2 or 4 angular views, but a similar performance with 8 and 16 views.

Cerebral SPECT imaging with different acquisition schemes using varying levels of multiplexing versus sensitivity in an adaptive multi-pinhole brain-dedicated scanner

Application of multi-pinhole collimator in pinhole-based SPECT increases detection sensitivity. The presence of multiplexing in projection images due to the usage of multiple pinholes can further improve the sensitivity at the cost of adding data ambiguity. We are developing a next-generation adaptive brain-dedicated SPECT system -AdaptiSPECT-C. The AdaptiSPECT-C can adapt the multiplexing level and system sensitivity using adaptable pinhole modules. In this study, we investigated the performance of 4 data acquisition schemes with different multiplexing levels and sensitivities on cerebral SPECT imaging. Schemes #1, #2, and #3 have <1%, 67%, and 31% overall multiplexing, respectively, while the 4th scheme without multiplexing is considered as ground truth. The ground-truth and schemes #1-3 have 1.0, 1.7, 5.1, and 4.0 times higher sensitivity, respectively, compared to a dual-headed parallel-hole SPECT system at matched spatial resolution. A customized XCAT brain perfusion digital phantom emulating the distribution of I-123 N-isopropyl iodoamphetamine (IMP) in a 99th percentile size male was used for simulations. Data acquisition for each scheme was performed at two count levels (low-count and high-count relative to the recommended clinical count level). The normalized root-mean-square error (NRMSE) for schemes #1, #2, and #3 with the low-count (high-count) scenario showed 11%, 4%, and 5% (10%, 5%, and 6%) deviation, respectively, from that of the multiplex-free ground truth. For both the low-count and high-count scenarios, scheme #1 resulted in the least accurate activity ratio (AR) for almost all the analyzed gray-matter brain regions. Further schemes #2 or #3 led to the most accurate AR values with both low-count and high-count scenarios for all the analyzed gray-matter regions. It was thus observed that even with this large head size which leads to significant multiplexing levels, the higher sensitivity from multiplexing could to some extent mitigate the data ambiguity and be translated into reconstructed images of higher quality.

Investigation of Axial and Angular Sampling in Multi-Detector Pinhole-SPECT Brain Imaging

We designed a dedicated multi-detector multi-pinhole brain SPECT scanner to generate images of higher quality compared to general-purpose systems. The system, AdaptiSPECT-C, is intended to adapt its sensitivity-resolution trade-off by varying its aperture configurations allowing both high-sensitivity dynamic and high-spatial-resolution static imaging. The current system design consists of 23 detector heads arranged in a truncated spherical geometry. In this work, we investigated the axial and angular sampling capability of the current stationary system design. Two data acquisition schemes using limited rotation of the gantry and two others using axial translation of the imaging bed were also evaluated concerning their impact on image quality through improved sampling. Increasing both angular and axial sampling in the current prototype system resulted in quantitative improvements in image quality metrics and qualitative appearance of the images as determined in studies with specifically selected phantoms. Visual improvements for the brain phantoms with clinical distributions were less pronounced but presented quantitative improvements in the fidelity (normalized root-mean-square error (NRMSE)) and striatal specific binding ratio (SBR) for a dopamine transporter (DAT) distribution, and in NRMSE and activity recovery for a brain perfusion distribution. More pronounced improvements with increased sampling were seen in contrast recovery coefficient, bias, and coefficient of variation for a lesion in the brain perfusion distribution. The negligible impact of the most cranial ring of detectors on axial sampling, but its significant impact on sensitivity and angular sampling in the cranial portion of the imaging volume-of-interest were also determined.

Deep learning based methods for gamma ray interaction location estimation in monolithic scintillation crystal detectors

In this work, we explore deep learning based techniques using the information from mean detector response functions (MDRFs) as a new method to estimate gamma ray interaction location in monolithic scintillation crystal detectors. Compared with searching based methods, deep learning techniques do not require recording all the MDRF information once the prediction networks are trained, which means the memory cost could be significantly reduced. In addition, the event positioning process using deep learning techniques only requires running through the network once, without the need to do searching in the reference dataset. This could greatly speed up the positioning process for each event. We have designed and trained four different neural networks to estimate the gamma ray interaction location given the MDRF data. We have studied network structures consisting only of fully connected (FC) layers, as well as Conv neural networks (CNNs). In addition, we tried to use both regression and classification to generate the final prediction of the gamma ray interaction position. We evaluated the estimation accuracy, testing speed and memory cost (numbers of parameters) of different network architectures, and also compared them with the exhaustive search method. Our results indicate that deep learning based estimation methods with a well designed network structure can achieve a relative positioning error with respect to the ground truth determined by the exhaustive search method of below 1 mm in both x and y directions (depth information is not considered in this work), which would imply a very high performance positioning algorithm for practical monolithic scintillation crystal detectors. The deep learning network also achieves a testing speed that is more than 400 times faster than the exhaustive search method. With proper design of the network structure, the deep learning based positioning methods have the potential to save memory cost by a factor of up to 100.

Primary, scatter, and penetration characterizations of parallel-hole and pinhole collimators for I-123 SPECT

Multi-pinhole (MPH) collimators are known to provide better trade-off between sensitivity and resolution for preclinical, as well as for smaller regions in clinical SPECT imaging compared to conventional collimators. In addition to this geometric advantage, MPH plates typically offer better stopping power for penetration than the conventional collimators, which is especially relevant for I-123 imaging. The I-123 emits a series of high-energy (>300 keV, ~2.5% abundance) gamma photons in addition to the primary emission (159 keV, 83% abundance). Despite their low abundance, high-energy photons penetrate through a low-energy parallel-hole (LEHR) collimator much more readily than the 159 keV photons, resulting in large downscatter in the photopeak window. In this work, we investigate the primary, scatter, and penetration characteristics of a single pinhole collimator that is commonly used for I-123 thyroid imaging and our two MPH collimators designed for I-123 DaTscan imaging for Parkinson's Disease, in comparison to three different parallel-hole collimators through a series of experiments and Monte Carlo simulations. The simulations of a point source and a digital human phantom with DaTscan activity distribution showed that our MPH collimators provide superior count performance in terms of high primary counts, low penetration, and low scatter counts compared to the parallel-hole and single pinhole collimators. For example, total scatter, multiple scatter, and collimator penetration events for the LEHR were 2.5, 7.6 and 14 times more than that of MPH within the 15% photopeak window. The total scatter fraction for LEHR was 56% where the largest contribution came from the high-energy scatter from the back compartments (31%). For the same energy window, the total scatter for MPH was 21% with only 1% scatter from the back compartments. We therefore anticipate that using MPH collimators, higher quality reconstructions can be obtained in a substantially shorter acquisition time for I-123 DaTscan and thyroid imaging.

Towards continuous-to-continuous 3D imaging in the real world

Imaging systems are often modeled as continuous-to-discrete mappings that map the object (i.e. a function of continuous variables such as space, time, energy, wavelength, etc) to a finite set of measurements. When it comes to reconstruction, some discretized version of the object is almost always assumed, leading to a discrete-to-discrete representation of the imaging system. In this paper, we discuss a method for single-photon emission computed tomography (SPECT) imaging that avoids discrete representations of the object or the imaging system, thus allowing reconstruction on an arbitrarily fine set of points.

Fast gamma-ray interaction-position estimation using k-d tree search

We have developed a fast gamma-ray interaction-position estimation method using k-d tree search, which can be combined with various kinds of closeness metrics such as Euclidean distance, maximum-likelihood estimation, etc. Compared with traditional search strategies, this method can achieve both speed and accuracy at the same time using the k-d tree data structure. The k-d tree search method has a time complexity of [Formula: see text], where N is the number of entries in the reference data set, which means large reference datasets can be used to efficiently estimate each event's interaction position. This method's accuracy was found to be equal to that of the exhaustive search method, yielding the highest achievable accuracy. Most importantly, this method has no restriction on the data structure of the reference dataset and can still work with complicated mean-detector-response functions (MDRFs), meaning that it is more robust than other popular methods such as contracting-grid-search (CG) or vector-search (VS) methods that could yield locally optimal instead of globally optimal results.

Imaging assessment of cardioprotection mediated by a dodecafluoropentane oxygen-carrier administered during myocardial infarction

INTRODUCTION

The objective of this study was to investigate the cardioprotective effects of a dodecafluoropentane (DDFP)-based perfluorocarbon emulsion (DDFPe) as an artificial carrier for oxygen delivery to ischemic myocardium, using ^99mTc-duramycin SPECT imaging.

METHODS

Rat hearts with Ischemia-reperfusion (I/R) was prepared by coronary ligation for 45-min followed by reperfusion. The feasibility of ^99mTc-duramycin in detecting myocardial I/R injury and its kinetic profile were first verified in the ischemic hearts with 2-h reperfusion (n = 6). DDFPe (0.6 mL/kg) was intravenously administered at 10 min after coronary ligation in fifteen rats and saline was given in thirteen rats as controls. ^99mTc-duramycin SPECT images were acquired in the DDFPe-treated hearts and saline controls at 2-h (DDFPe-2 h, n = 7 and Saline-2 h, n = 6) or 24-h (DDFPe-24 h, n = 8 and Saline-24 h, n = 7) of reperfusion.

RESULTS

SPECT images, showing "hot-spot" ^99mTc-duramycin uptake in the ischemic myocardium, exhibited significantly lower radioactive retention and smaller hot-spot size in the DDFPe-2 h and DDFPe-24 h hearts compared to controls. The infarcts in the Saline-24 h hearts extended significantly relative to measurements in the Saline-2 h. The extension of infarct size did not reach a statistical difference between the DDFPe-2 h and DDFPe-24 h hearts. Ex vivo measurement of ^99mTc-duramycin activity (%ID/g) was lower in the ischemic area of DDFPe-2 h and DDFPe-24 h than that of the Saline-2 h and Saline-24 h hearts (P < 0.05). The area of injured myocardium, delineated by the uptake of ^99mTc-duramycin, extended more substantially outside the infarct zone in the controls.

CONCLUSIONS

Significant reduction in myocardial I/R injury, as assessed by ^99mTc-duramycin cell death imaging and histopathological analysis, was induced by DDFPe treatment after acute myocardial ischemia. ^99mTc-duramycin imaging can reveal myocardial cell death in ischemic hearts and may provide a tool for the non-invasive assessment of cardioprotective interventions.

An edge-readout, multilayer detector for positron emission tomography

PURPOSE

We present a novel gamma-ray-detector design based on total internal reflection (TIR) of scintillation photons within a crystal that addresses many limitations of traditional PET detectors. Our approach has appealing features, including submillimeter lateral resolution, DOI positioning from layer thickness, and excellent energy resolution. The design places light sensors on the edges of a stack of scintillator slabs separated by small air gaps and exploits the phenomenon that more than 80% of scintillation light emitted during a gamma-ray event reaches the edges of a thin crystal with polished faces due to TIR. Gamma-ray stopping power is achieved by stacking multiple layers, and DOI is determined by which layer the gamma ray interacts in.

METHOD

The concept of edge readouts of a thin slab was verified by Monte Carlo simulation of scintillation light transport. An LYSO crystal of dimensions 50.8 mm × 50.8 mm × 3.0 mm was modeled with five rectangular SiPMs placed along each edge face. The mean-detector-response functions (MDRFs) were calculated by simulating signals from 511 keV gamma-ray interactions in a grid of locations. Simulations were carried out to study the influence of choice of scintillator material and dimensions, gamma-ray photon energies, introduction of laser or mechanically induced optical barriers (LIOBs, MIOBs), and refractive indices of optical-coupling media and SiPM windows. We also analyzed timing performance including influence of gamma-ray interaction position and presence of optical barriers. We also modeled and built a prototype detector, a 27.4 mm × 27.4 mm × 3.0 mm CsI(Tl) crystal with 4 SiPMs per edge to experimentally validate the results predicted by the simulations. The prototype detector used CsI(Tl) crystals from Proteus outfitted with 16 Hamamatsu model S13360-6050PE MPPCs read out by an AiT-16-channel readout. The MDRFs were measured by scanning the detector with a collimated beam of 662-keV photons from a ¹³⁷ Cs source. The spatial resolution was experimentally determined by imaging a tungsten slit that created a beam of 0.44 mm (FWHM) width normal to the detector surface. The energy resolution was evaluated by analyzing list-mode data from flood illumination by the ¹³⁷ Cs source.

RESULT

We find that in a block-detector-sized LYSO layer read out by five SiPMs per edge, illuminated by 511-keV photons, the average resolution is 1.49 mm (FWHM). With the introduction of optical barriers, average spatial resolution improves to 0.56 mm (FWHM). The DOI resolution is the layer thickness of 3.0 mm. We also find that optical-coupling media and SiPM-window materials have an impact on spatial resolution. The timing simulation with LYSO crystal yields a coincidence resolving time (CRT) of 200-400 ps, which is slightly position dependent. And the introduction of optical barriers has minimum influence. The prototype CsI(Tl) detector, with a smaller area and fewer SiPMs, was measured to have central-area spatial resolutions of 0.70 and 0.39 mm without and with optical barriers, respectively. These results match well with our simulations. An energy resolution of 6.4% was achieved at 662 keV.

CONCLUSION

A detector design based on a stack of monolithic scintillator layers that uses edge readouts offers several advantages over current block detectors for PET. For example, there is no tradeoff between spatial resolution and detection sensitivity since no reflector material displaces scintillator crystal, and submillimeter resolution can be achieved. DOI information is readily available, and excellent timing and energy resolutions are possible.

Radiolabeled methotrexate as a diagnostic agent of inflammatory target sites: A proof-of-concept study

Methotrexate (MTX), as a pharmaceutical, is frequently used in tumor chemotherapy and is also a part of the established treatment of a number of autoimmune inflammatory disorders. Radiolabeled MTX has been studied as a tumor-diagnostic agent in a number of published studies. In the present study, the potential use of technetium-99m-labelled MTX (^99mTc-MTX) as a radiotracer was investigated for the identification of inflammatory target sites. The labelling of MTX was carried out via a ^99mTc-gluconate precursor. Evaluation studies included in vitro stability, plasma protein binding assessment, partition-coefficient estimation, in vivo scintigraphic imaging and ex vivo animal experiments in an animal inflammation model. MTX was successfully labelled with ^99mTc, with a radiochemical purity of >95%. Stability was assessed in plasma, where it remained intact up to 85% at 4 h post-incubation, while protein binding of the radiotracer was observed to be ~50% at 4 h. These preclinical ex vivo and in vivo studies indicated that ^99mTc-MTX accumulates in inflamed tissue, as well as in the spinal cord, joints and bones; all areas with relatively high remodeling activity. The results are promising, and set the stage for further work on the development and application of ^99mTc-MTX as a radiotracer for inflammation associated with rheumatoid arthritis.

Maximum-Likelihood Estimation of Scintillation Pulse Timing

Including time-of-flight information in positron emission tomography (PET) reconstruction increases the signal-to-noise ratio if the timing information is sufficiently accurate. We estimate timing information by analyzing sampled waveforms, where the sampling frequency and number of samples acquired affect the accuracy of timing estimation. An efficient data-acquisition system acquires the minimum number of samples that contains the most timing information for a desired resolution. We describe a maximum-likelihood (ML) estimation algorithm to assign a time stamp to digital pulses. The method is based on a contracting-grid search algorithm that can be implemented in a field-programmable gate array and in graphics processing units. The Fisher-information (FI) matrix quantifies the amount of timing information that can be extracted from the waveforms. FI analyses on different segments of the waveform allow us to determine the smallest amount of data that we need to acquire in order to obtain a desired timing resolution. We describe the model and the procedure used to simulate waveforms for ML estimation and FI analysis, the ML-estimation algorithm and the timing resolution obtained from experimental data using a LaBr₃:Ce crystal and two photomultiplier tubes. The results show that for lengthening segments of the pulse, timing resolution approaches a limit. We explored the method as a function of sampling frequency and compared the results to other digital time pickoff methods. This information will be used to build an efficient data-acquisition system with reduced complexity and cost that nonetheless preserves full timing performance.

Characterization of TCP-1 probes for molecular imaging of colon cancer

Molecular probes capable of detecting colorectal cancer (CRC) are needed for early CRC diagnosis. The objective of this study was to characterize c[CTPSPFSHC]OH (TCP-1), a small peptide derived from phage display selection, for targeting human CRC xenografts using technetium-99m ((99m)Tc)-labeled TCP-1 and fluorescent cyanine-7 (Cy7)-labeled form of the peptide (Cy7-TCP-1). (99m)Tc-TCP-1 was generated by modifying TCP-1 with succinimidyl-6-hydrazino-nicotinamide (S-HYNIC) followed by radiolabeling. In vitro saturation binding experiments were performed for (99m)Tc-TCP-1 in human HCT116 colon cancer cells. SCID mice with human HCT116 cancer xenografts were imaged with (99m)Tc-TCP-1 or control peptide using a small-animal SPECT imager: Group I (n=5) received no blockade; Group II (n=5) received a blocking dose of non-radiolabeled TCP-1. Group III (n=5) were imaged with (99m)Tc-labeled control peptide (inactive peptide). SCID mice with human PC3 prostate cancer xenografts (Group IV, n=5) were also imaged with (99m)Tc-TCP-1. Eight additional SCID mice bearing HCT116 xenografts in dorsal skinfold window chambers (DSWC) were imaged by direct positron imaging of (18)F-fluorodeoxyglucose ((18)F-FDG) and fluorescence microscopy of Cy7-TCP-1. In vitro(99m)Tc-HYNIC-TCP-1 binding assays on HCT 116 cells indicated a mean Kd of 3.04±0.52nM. In cancer xenografts, (99m)Tc-TCP-1 radioactivity (%ID/g) was 1.01±0.15 in the absence of blockade and was reduced to 0.26±0.04 (P<0.01) with blockade. No radioactive uptake was observed in the PC3 tumors with (99m)Tc-TCP-1 or HCT116 tumors with inactive peptide. Cy7-TCP-1 activity localized not only in metabolically active tumors, as defined by (18)F-FDG imaging, but also in peritumoral microvasculature. In conclusion, TCP-1 probes may have a distinct targeting mechanism with high selectivity for CRC and tumor-associated vasculature. Molecular imaging with TCP-1 probes appears promising to detect malignant colorectal lesions.

Task-based optimization of flip angle for fibrosis detection in T1-weighted MRI of liver

Chronic liver disease is a worldwide health problem, and hepatic fibrosis (HF) is one of the hallmarks of the disease. The current reference standard for diagnosing HF is biopsy followed by pathologist examination; however, this is limited by sampling error and carries a risk of complications. Pathology diagnosis of HF is based on textural change in the liver as a lobular collagen network that develops within portal triads. The scale of collagen lobules is characteristically in the order of 1 to 5 mm, which approximates the resolution limit of in vivo gadolinium-enhanced magnetic resonance imaging in the delayed phase. We use MRI of formalin-fixed human ex vivo liver samples as phantoms that mimic the textural contrast of in vivo Gd-MRI. We have developed a local texture analysis that is applied to phantom images, and the results are used to train model observers to detect HF. The performance of the observer is assessed with the area-under-the-receiver-operator-characteristic curve (AUROC) as the figure-of-merit. To optimize the MRI pulse sequence, phantoms were scanned with multiple times at a range of flip angles. The flip angle that was associated with the highest AUROC was chosen as optimal for the task of detecting HF.

Detection of atherosclerotic plaques in ApoE-deficient mice using (99m)Tc-duramycin

Apoptosis of macrophages and smooth muscle cells is linked to atherosclerotic plaque destabilization. The apoptotic cascade leads to exposure of phosphatidylethanolamine (PE) on the outer leaflet of the cell membrane, thereby making apoptosis detectable using probes targeting PE. The objective of this study was to exploit capabilities of a PE-specific imaging probe, (99m)Tc-duramycin, in localizing atherosclerotic plaque and assessing plaque evolution in apolipoprotein-E knockout (ApoE(-/-)) mice.

METHODS

Atherosclerosis was induced in ApoE(-/-) mice by feeding an atherogenic diet. (99m)Tc-duramycin images were acquired using a small-animal SPECT imager. Six ApoE(-/-) mice at 20weeks of age (Group I) were imaged and then sacrificed for ex vivo analyses. Six additional ApoE(-/-) mice (Group II) were imaged at 20 and 40weeks of age before sacrifice. Six ApoE wild-type (ApoE(+/+)) mice (Group III) were imaged at 40weeks as controls. Five additional ApoE(-/-) mice (40weeks of age) (Group IV) were imaged with a (99m)Tc-labeled inactive peptide, (99m)Tc-LinDUR, to assess (99m)Tc-duramycin targeting specificity.

RESULTS

Focal (99m)Tc-duramycin uptake in the ascending aorta and aortic arch was detected at 20 and 40weeks in the ApoE(-/-) mice but not in ApoE(+/+) mice. (99m)Tc-duramycin uptake in the aortic lesions increased 2.2-fold on quantitative imaging in the ApoE(-/-) mice between 20 and 40weeks. Autoradiographic and histological data indicated significantly increased (99m)Tc-duramycin uptake in the ascending aorta and aortic arch associated with advanced plaques. Quantitative autoradiography showed that the ratio of activity in the aortic arch to descending thoracic aorta, which had no plaques or radioactive uptake, was 2.1 times higher at 40weeks than at 20weeks (6.62±0.89 vs. 3.18±0.29, P<0.01). There was barely detectable focal uptake of (99m)Tc-duramycin in the aortic arch of ApoE(+/+) mice. No detectable (99m)Tc-LinDUR uptake was observed in the aortas of ApoE(-/-) mice.

CONCLUSIONS

PE-targeting properties of (99m)Tc-duramycin in the atherosclerotic mouse aortas were noninvasively characterized. (99m)Tc-duramycin is promising in localizing advanced atherosclerotic plaques.

Molecular Imaging in the College of Optical Sciences - An Overview of Two Decades of Instrumentation Development

During the past two decades, researchers at the University of Arizona's Center for Gamma-Ray Imaging (CGRI) have explored a variety of approaches to gamma-ray detection, including scintillation cameras, solid-state detectors, and hybrids such as the intensified Quantum Imaging Device (iQID) configuration where a scintillator is followed by optical gain and a fast CCD or CMOS camera. We have combined these detectors with a variety of collimation schemes, including single and multiple pinholes, parallel-hole collimators, synthetic apertures, and anamorphic crossed slits, to build a large number of preclinical molecular-imaging systems that perform Single-Photon Emission Computed Tomography (SPECT), Positron Emission Tomography (PET), and X-Ray Computed Tomography (CT). In this paper, we discuss the themes and methods we have developed over the years to record and fully use the information content carried by every detected gamma-ray photon.

A New Columnar CsI(Tl) Scintillator for iQID detectors

A 1650 μm thick columnar CsI(Tl) scintillator for upgrading iQID detectors, which is a high-resolution photon-counting gamma-ray and x-ray detector recently developed at the Center for Gamma-Ray Imaging (CGRI), has been studied in terms of sensitivity, spatial resolution and depth-of-interaction effects. To facilitate these studies, a new frame-parsing algorithm for processing raw event data is also proposed that has more degrees of freedom in data processing and can discriminate against a special kind of noise present in some low-cost intensifiers. The results show that in comparison with a 450 μm-thickness columnar CsI(Tl) scintillator, the 1650 μm thick CsI(Tl) scintillator provides more than twice the sensitivity at the expense of some spatial resolution degradation. The depth-of-interaction study also shows that event size and amplitude vary with scintillator thickness, which can assist in future detector simulations and 3D-interaction-position estimation.

A SPECT system simulator built on the SolidWorks TM 3D-Design package

We have developed a GPU-accelerated SPECT system simulator that integrates into instrument-design workflow [1]. This simulator includes a gamma-ray tracing module that can rapidly propagate gamma-ray photons through arbitrary apertures modeled by SolidWorks ^TM -created stereolithography (.STL) representations with a full complement of physics cross sections [2, 3]. This software also contains a scintillation detector simulation module that can model a scintillation detector with arbitrary scintillation crystal shape and light-sensor arrangement. The gamma-ray tracing module enables us to efficiently model aperture and detector crystals in SolidWorks ^TM and save them as STL file format, then load the STL-format model into this module to generate list-mode results of interacted gamma-ray photon information (interaction positions and energies) inside the detector crystals. The Monte-Carlo scintillation detector simulation module enables us to simulate how scintillation photons get reflected, refracted and absorbed inside a scintillation detector, which contributes to more accurate simulation of a SPECT system.

Fabrication of the pinhole aperture for AdaptiSPECT

AdaptiSPECT is a pre-clinical pinhole SPECT imaging system under final construction at the Center for Gamma-Ray Imaging. The system is designed to be able to autonomously change its imaging configuration. The system comprises 16 detectors mounted on translational stages to move radially away and towards the center of the field-of-view. The system also possesses an adaptive pinhole aperture with multiple collimator diameters and pinhole sizes, as well as the possibility to switch between multiplexed and non-multiplexed imaging configurations. In this paper, we describe the fabrication of the AdaptiSPECT pinhole aperture and its controllers.

Adaptive SPECT imaging with crossed-slit apertures

Preclinical single-photon emission computed tomography (SPECT) is an essential tool for studying the progression, response to treatment, and physiological changes in small animal models of human disease. The wide range of imaging applications is often limited by the static design of many preclinical SPECT systems. We have developed a prototype imaging system that replaces the standard static pinhole aperture with two sets of movable, keel-edged copper-tungsten blades configured as crossed (skewed) slits. These apertures can be positioned independently between the object and detector, producing a continuum of imaging configurations in which the axial and transaxial magnifications are not constrained to be equal. We incorporated a megapixel silicon double-sided strip detector to permit ultrahigh-resolution imaging. We describe the configuration of the adjustable slit aperture imaging system and discuss its application toward adaptive imaging, and reconstruction techniques using an accurate imaging forward model, a novel geometric calibration technique, and a GPU-based ultra-high-resolution reconstruction code.

Imaging properties of pixellated scintillators with deep pixels

We have investigated the light-transport properties of scintillator arrays with long, thin pixels (deep pixels) for use in high-energy gamma-ray imaging. We compared 10×10 pixel arrays of YSO:Ce, LYSO:Ce and BGO (1mm × 1mm × 20 mm pixels) made by Proteus, Inc. with similar 10×10 arrays of LSO:Ce and BGO (1mm × 1mm × 15mm pixels) loaned to us by Saint-Gobain. The imaging and spectroscopic behaviors of these scintillator arrays are strongly affected by the choice of a reflector used as an inter-pixel spacer (3M ESR in the case of the Proteus arrays and white, diffuse-reflector for the Saint-Gobain arrays). We have constructed a 3700-pixel LYSO:Ce Prototype NIF Gamma-Ray Imager for use in diagnosing target compression in inertial confinement fusion. This system was tested at the OMEGA Laser and exhibited significant optical, inter-pixel cross-talk that was traced to the use of a single-layer of ESR film as an inter-pixel spacer. We show how the optical cross-talk can be mapped, and discuss correction procedures. We demonstrate a 10×10 YSO:Ce array as part of an iQID (formerly BazookaSPECT) imager and discuss issues related to the internal activity of ¹⁷⁶Lu in LSO:Ce and LYSO:Ce detectors.

Inflammation imaging of atherosclerosis in Apo-E-deficient mice using a (99m)Tc-labeled dual-domain cytokine ligand

Interleukin-1β (IL-1β) and tumor necrosis factor-α (TNF-α) play a critical role in initiating and accelerating atherosclerosis. This study evaluated the imaging properties of (99m)Tc-TNFR2-Fc-IL-1RA ((99m)Tc-TFI), a dual-domain cytokine radioligand that targets TNF-α and IL-1β pathways, in assessing atherosclerosis development in apolipoprotein-E-deficient (ApoE(-)(/)(-)) mice.

METHODS

The feasibility and specificity of detecting atherosclerosis with (99m)Tc-TFI SPECT imaging were investigated in ApoE(-)(/)(-) and ApoE(+)(/)(+) mice. Fifty-four ApoE(-)(/)(-) mice were fed either an atherogenic diet (AGD) or a normal diet (ND) beginning at 5 weeks of age. Eighteen Apo-E wild-type (ApoE(+)(/)(+)) mice were fed an ND. Two groups of ApoE(-)(/)(-) mice (n=12 each group) on AGD and ND were imaged three times with (99m)Tc-TFI and a high-resolution SPECT system at 20-25, 30-40, and 48-52 weeks to study the evolution of atherosclerotic plaque.

RESULTS

Focal radioactive accumulations in the aortic arch region were observed in the ApoE(-)(/)(-) mice (n=12) on AGD but not in the ApoE(+)(/)(+) mice on ND (n=10). Apo-E(-)(/)(-) mice on ND (n=11) exhibited lower radioactive uptake than ApoE(-)(/)(-) mice on AGD (P<0.05). Co-injection of an excess of cold ligand with (99m)Tc-TFI resulted in significant reduction of (99m)Tc-TFI uptake in the ApoE(-)(/)(-) mice on AGD. Longitudinal studies showed that (99m)Tc-TFI uptake in the aortas of ApoE(-)(/)(-) mice progressively increased from 20 to 48 weeks. Real-time PCR assays demonstrated that atherosclerotic aortas expressed significantly higher IL-1β and TNF-α than the aortas from wild-type controls.

CONCLUSIONS

Atherosclerotic plaques were detected by (99m)Tc-TFI imaging in ApoE(-)(/)(-) mice. (99m)Tc-TFI is promising for specific detection of inflammatory response in atherosclerotic plaques.

SPECT imaging of inflammatory response in ischemic-reperfused rat hearts using a 99mTc-labeled dual-domain cytokine ligand

Soluble tumor necrosis factor (TNF) receptor-2 (TNFR2) and interleukin-1 receptor antagonist (IL-1ra) were fused to the Fc portion of IgG1 using recombinant DNA technology. The resulting dual-domain cytokine ligand, TNFR2-Fc-IL-1ra, specifically binds to TNF and to the type I IL-1 receptor (IL-1RI). This study was designed to characterize the kinetic profile of (99m)Tc-labeled TNFR2-Fc-IL-1ra (TFI) for imaging inflammatory response in an ischemic-reperfused (IR) rat heart model.

METHODS

The IR model was created by ligating the left coronary artery for 45 min, followed by 2-h reperfusion. Cardiac SPECT images of TFI in the IR model (n = 6) were dynamically acquired for 3 h. Correlative data of myocardial TFI distribution versus microsphere-determined tissue blood flow were acquired in 3 extra IR hearts. Inflammation targeting affinity of TFI was compared with 2 individual cytokine radioligands, (99m)Tc-IL-1ra-Fc (IF) and (99m)Tc-TNFR2-Fc (TF) (n = 6 each group). Myocardial cytokine expression was evaluated by immunochemical assay.

RESULTS

Increased TFI uptake was found in the ischemic area and correlated with the severity of ischemia. At 3 h after injection, the ratio of hot-spot accumulation in the ischemic area to a remote viable zone was 5.39 ± 1.11 for TFI, which was greater than that for IF (3.28 ± 0.81) and TF (3.29 ± 0.75) (P < 0.05). The in vivo uptake profiles of TFI, TF, and IF were consistent with ex vivo radioactive measurements and correlated with upregulated IL-1 and TNF expression.

CONCLUSION

The dual-domain TFI is promising for noninvasive detection of inflammatory reactions in IR myocardium because of its more potent affinity to the inflammatory sites compared with TF and IF.

Image Science with Photon-Processing Detectors

We introduce and discuss photon-processing detectors and we compare them with photon-counting detectors. By estimating a relatively small number of attributes for each collected photon, photon-processing detectors may help understand and solve a fundamental theoretical problem of any imaging system based on photon-counting detectors, namely null functions. We argue that photon-processing detectors can improve task performance by estimating position, energy, and time of arrival for each collected photon. We consider a continuous-to-continuous linear operator to relate the object being imaged to the collected data, and discuss how this operator can be analyzed to derive properties of the imaging system. Finally, we derive an expression for the characteristic functional of an imaging system that produces list-mode data.

A SPECT imager with synthetic collimation

This work outlines the development of a multi-pinhole SPECT system designed to produce a synthetic-collimator image of a small field of view. The focused multi-pinhole collimator was constructed using rapid-prototyping and casting techniques. The collimator projects the field of view through forty-six pinholes when the detector is adjacent to the collimator. The detector is then moved further from the collimator to increase the magnification of the system. The amount of pinhole-projection overlap increases with the system magnification. There is no rotation in the system; a single tomographic angle is used in each system configuration. The maximum-likelihood expectation-maximization (MLEM) algorithm is implemented on graphics processing units to reconstruct the object in the field of view. Iterative reconstruction algorithms, such as MLEM, require an accurate model of the system response. For each system magnification, a sparsely-sampled system response is measured by translating a point source through a grid encompassing the field of view. The pinhole projections are individually identified and associated with their respective apertures. A 2D elliptical Gaussian model is applied to the pinhole projections on the detector. These coefficients are associated with the object-space location of the point source, and a finely-sampled system matrix is interpolated. Simulations with a hot-rod phantom demonstrate the efficacy of combining low-resolution non-multiplexed data with high-resolution multiplexed data to produce high-resolution reconstructions.

The effect of gain variation in micro-channel plates on gamma-ray energy resolution

A Monte Carlo simulation of micro-channel plate (MCP) with particular interest in its effect on energy resolution performance is presented. Important physical processes occurring in MCP channels are described and modeled, including secondary electron (SE) yield, SE emission, and primary electron reflection. The effects causing channel saturation are also introduced. A two dimensional Monte Carlo simulation is implemented under the assumption of unsaturated channel. Simulation results about basic MCP performances and especially gain and energy resolution performances are presented and analyzed. It's found that energy resolution as an intrinsic property of MCP cannot be improved simply by adjusting system parameters; however it can be improved by increasing input signal or number of photoelectrons (PEs) in the context of image intensifier. An initial experiment with BazookaSPECT detector and CsI(Tl) scintillator is performed to validate and correlate with the simulation results and good agreement is achieved.

Integration of AdaptiSPECT, a small-animal adaptive SPECT imaging system

AdaptiSPECT is a pre-clinical adaptive SPECT imaging system under final development at the Center for Gamma-ray Imaging. The system incorporates multiple adaptive features: an adaptive aperture, 16 detectors mounted on translational stages, and the ability to switch between a non-multiplexed and a multiplexed imaging configuration. In this paper, we review the design of AdaptiSPECT and its adaptive features. We then describe the on-going integration of the imaging system.

SmartCAM: An adaptive clinical SPECT camera

An adaptive pinhole aperture that fits a GE MaxiCam Single-Photon-Emission Computed Tomography (SPECT) system has been designed, built, and is undergoing testing. The purpose of an adaptive aperture is to allow the imaging system to make adjustments to the aperture while imaging data are being acquired. Our adaptive pinhole aperture can alter several imaging parameters, including field of view, resolution, sensitivity, and magnification. The dynamic nature of such an aperture allows for imaging of specific regions of interest based on initial measurements of the patient. Ideally, this mode of data collection will improve the understanding of a patient's condition, and will facilitate better diagnosis and treatment. The aperture was constructed using aluminum and a low melting point, high-stopping-power metal alloy called Cerrobend. The aperture utilizes a rotating disk for the selection of a pinhole configuration; as the aluminum disk rotates, different pinholes move into view of the camera face and allow the passage of gamma rays through that particular pinhole. By controlling the angular position of the disk, the optical characteristics of the aperture can be modified, allowing the system to acquire data from controlled regions of interest. First testing was performed with a small radioactive source to prove the functionality of the aperture.

Design and performance of a small-animal imaging system using synthetic collimation

This work outlines the design and construction of a single-photon emission computed tomography imaging system based on the concept of synthetic collimation. A focused multi-pinhole collimator is constructed using rapid-prototyping and casting techniques. The collimator projects the centre of the field of view (FOV) through 46 pinholes when the detector is adjacent to the collimator, with the number reducing towards the edge of the FOV. The detector is then moved further from the collimator to increase the magnification of the system. The object distance remains constant, and each new detector distance is a new system configuration. The level of overlap of the pinhole projections increases as the system magnification increases, but the number of projections subtended by the detector is reduced. There is no rotation in the system; a single tomographic angle is used in each system configuration. Image reconstruction is performed using maximum-likelihood expectation-maximization and an experimentally measured system matrix. The system matrix is measured for each configuration on a coarse grid, using a point source. The pinholes are individually identified and tracked, and a Gaussian fit is made to each projection. The coefficients of these fits are used to interpolate the system matrix. The system is validated experimentally with a hot-rod phantom. The Fourier crosstalk matrix is also measured to provide an estimate of the average spatial resolution along each axis over the entire FOV. The 3D synthetic-collimator image is formed by estimating the activity distribution within the FOV and summing the activities in the voxels along the axis perpendicular to the collimator face.

Imaging of rat cerebral ischemia-reperfusion injury using(99m)Tc-labeled duramycin

OBJECTIVES

Prompt identification of necrosis and apoptosis in the infarct core and penumbra region is critical in acute stroke for delineating the underlying ischemic/reperfusion molecular pathologic events and defining therapeutic alternatives. The objective of this study was to investigate the capability of (99m)Tc-labeled duramycin in detecting ischemia-reperfusion injury in rat brain after middle cerebral artery (MCA) occlusion.

METHODS

Ischemic cerebral injury was induced in ten rats by vascular insertion of a nylon suture in the left MCA for 3 hr followed by 21-24hr reperfusion. After i.v. injection of (99m)Tc-duramycin (1.0-3.5 mCi), dynamic cerebral images were acquired for 1 hr in six rats using a small-animal SPECT imager. Four other rats were imaged at 2 hr post-injection. Ex vivo images were obtained by autoradiography after sacrifice. Histologic analyses were performed to assess cerebral infarction and apoptosis.

RESULTS

SPECT images showed that (99m)Tc-duramycin uptake in the left cerebral hemisphere was significantly higher than that in the right at 1 and 2 hr post-injection. The level of radioactive uptake in the ischemic brain varied based on ischemic severity. The average ratio of left cerebral hot-spot uptake to right hemisphere radioactivity, as determined by computerized ROI analysis, was 4.92±0.79. Fractional washout at 1 hr was 38.2±4.5% of peak activity for left cerebral hot-spot areas and 80.9±2.0% for remote control areas (P<0.001). Based on triphenyltetrazolium chloride staining and autoradiograph image data, the hotspot uptake may be associated primarily with the ischemic penumbra, in which high apoptotic activity was observed by cleaved caspase-3 immunocytochemical staining.

CONCLUSIONS

(99m)Tc-duramycin SPECT imaging may be useful for detecting and quantifying ongoing apoptotic neuronal cell loss induced by ischemia-reperfusion injury.

The AdaptiSPECT Imaging Aperture

In this paper, we present the imaging aperture of an adaptive SPECT imaging system being developed at the Center for Gamma Ray Imaging (AdaptiSPECT). AdaptiSPECT is designed to automatically change its configuration in response to preliminary data, in order to improve image quality for a particular task. In a traditional pinhole SPECT imaging system, the characteristics (magnification, resolution, field of view) are set by the geometry of the system, and any modification can be accomplished only by manually changing the collimator and the distance of the detector to the center of the field of view. Optimization of the imaging system for a specific task on a specific individual is therefore difficult. In an adaptive SPECT imaging system, on the other hand, the configuration can be conveniently changed under computer control. A key component of an adaptive SPECT system is its aperture. In this paper, we present the design, specifications, and fabrication of the adaptive pinhole aperture that will be used for AdaptiSPECT, as well as the controls that enable autonomous adaptation.

A Prototype Detector for a Novel High-Resolution PET System: BazookaPET

We have designed and are developing a novel proof-of-concept PET system called BazookaPET. In order to complete the PET configuration, at least two detector elements are required to detect positron-electron annihilation events. Each detector element of the BazookaPET has two independent data acquisition channels. One side of the scintillation crystal is optically coupled to a 4×4 silicon photomultiplier (SiPM) array and the other side is a CCD-based gamma camera. Using these two separate channels, we can obtain data with high energy, temporal and spatial resolution data by associating the data outputs via several maximum-likelihood estimation (MLE) steps. In this work, we present the concept of the system and the prototype detector element. We focus on characterizing individual detector channels, and initial experimental calibration results are shown along with preliminary performance-evaluation results. We measured energy resolution and the integrated traces of the slit-beam images from both detector channel outputs. A photo-peak energy resolution of ~5.3% FWHM was obtained from the SiPM and ~48% FWHM from the CCD at 662 keV. We assumed SiPM signals follow Gaussian statistics and estimated the 2D interaction position using MLE. Based on our the calibration experiments, we computed the Cramér-Rao bound (CRB) for the SiPM detector channel and found that the CRB resolution is better than 1 mm in the center of the crystal.

Event Processing for Modular Gamma Cameras with Tiled Multi-Anode Photomultiplier Tubes

Multi-anode photomultiplier tubes (MAPMTs) are good candidates as light sensors for a new generation of modular scintillation cameras for Single-photon emission computed tomography (SPECT) and Positron emission tomography (PET) applications. MAPMTs can provide improved intrinsic spatial resolution (<1mm) compared to arrays of larger individual PMTs due to their small anode sizes, and the increased number of channels also allows accurate estimation of depth-of-interaction (DOI). However, the area of a single MAPMT module is small for a modular gamma camera, so we are designing read-out electronics that will allow multiple individual MAPMT modules to be optically coupled to a single monolithic scintillator crystal. In order to allow such flexibility, the read-out electronics, which we refer to as the event processor, must be compact and adaptable. In combining arrays of MAPMTs, which may each have 64 to 1024 anodes per unit, issues need to be overcome with amplifying, digitizing, and recording potentially very large numbers of channels per gamma-ray event. In this study, we have investigated different event-processor strategies for gamma cameras with multiple MAPMTs that will employ maximum-likelihood (ML) methods for estimation of 3D spatial location, deposited energy and time of occurrence of events. We simulated anode signals for hypothetical gamma-camera geometries based on models of the stochastic processes inherent in scintillation cameras. The comparison between different triggering and read-out schemes was carried out by quantifying the information content in the anode signals via the Fisher Information Matrix (FIM). We observed that a decline in spatial resolution at the edges of the individual MAPMTs could be improved by the inclusion of neighboring MAPMT anode signals for events near the tiling boundaries. Thus in order to maintain spatial resolution uniformity throughout the modular camera face, we propose dividing an MAPMT's array of anode signals into regions such to help determine when triggers from one MAPMT need to be passed to a neighboring MAPMT so that it can contribute anode information for events between them.

A System Calibration and Fast Iterative Reconstruction Method for Next-Generation SPECT Imagers

Recently, high-resolution gamma cameras have been developed with detectors containing> 10⁵-10⁶ elements. SPECT imagers based on these detectors usually also have a large number of voxel bins and therefore face memory storage issues for the system matrix when performing fast tomographic reconstructions using iterative algorithms. To address these issues, we have developed a method that parameterizes the detector response to a point source and generates the system matrix on the fly during MLEM or OSEM on graphics hardware. The calibration method, interpolation of coefficient data, and reconstruction results are presented in the context of a recently commissioned small-animal SPECT imager, called FastSPECT III.

High-resolution, anamorphic, adaptive small-animal SPECT imaging with silicon double-sided strip detectors

We are developing a prototype gamma-ray imaging system that consists of two sets of movable, keel-edged copper-tungsten blades configured as crossed slits. These apertures can be positioned independently between the object and detector, producing an anamorphic image in which the axial and transaxial magnifications are not constrained to be equal. The detector is a 60 mm × 60 mm, millimeter thick, one-megapixel silicon double-sided strip detector. The flexible nature of this system allows the application of adaptive imaging techniques. We will discuss system details, calibration and acquisition methods, and our progress towards biological imaging applications.

Investigation of the possibility of gamma-ray diagnostic imaging of target compression at NIF

The National Ignition Facility at Lawrence Livermore National Laboratory is the world's leading facility to study the physics of igniting plasmas. Plasmas of hot deuterium and tritium, undergo d(t,n)α reactions that produce a 14.1 MeV neutron and 3.5 MeV a particle, in the center of mass. As these neutrons pass through the materials surrounding the hot core, they may undergo subsequent (n,x) reactions. For example, (12)C(n,n'γ)(12)C reactions occur in remnant debris from the polymer ablator resulting in a significant fluence of 4.44 MeV gamma-rays. Imaging of these gammas will enable the determination of the volumetric size and symmetry of the ablation; large size and high asymmetry is expected to correlate with poor compression and lower fusion yield. Results from a gamma-ray imaging system are expected to be complimentary to a neutron imaging diagnostic system already in place at the NIF. This paper describes initial efforts to design a gamma-ray imaging system for the NIF using the existing neutron imaging system as a baseline for study. Due to the cross-section and expected range of ablator areal densities, the gamma flux should be approximately 10(-3) of the neutron flux. For this reason, care must be taken to maximize the efficiency of the gamma-ray imaging system because it will be gamma starved. As with the neutron imager, use of pinholes and/or coded apertures are anticipated. Along with aperture and detector design, the selection of an appropriate scintillator is discussed. The volume of energy deposition of the interacting 4.44 MeV gamma-rays is a critical parameter limiting the imaging system spatial resolution. The volume of energy deposition is simulated with GEANT4, and plans to measure the volume of energy deposition experimentally are described. Results of tests on a pixellated LYSO scintillator are also presented.

SPECT detectors: the Anger Camera and beyond

The development of radiation detectors capable of delivering spatial information about gamma-ray interactions was one of the key enabling technologies for nuclear medicine imaging and, eventually, single-photon emission computed tomography (SPECT). The continuous sodium iodide scintillator crystal coupled to an array of photomultiplier tubes, almost universally referred to as the Anger Camera after its inventor, has long been the dominant SPECT detector system. Nevertheless, many alternative materials and configurations have been investigated over the years. Technological advances as well as the emerging importance of specialized applications, such as cardiac and preclinical imaging, have spurred innovation such that alternatives to the Anger Camera are now part of commercial imaging systems. Increased computing power has made it practical to apply advanced signal processing and estimation schemes to make better use of the information contained in the detector signals. In this review we discuss the key performance properties of SPECT detectors and survey developments in both scintillator and semiconductor detectors and their readouts with an eye toward some of the practical issues at least in part responsible for the continuing prevalence of the Anger Camera in the clinic.

Design Considerations for the Next-Generation MAPMT-Based Monolithic Scintillation Camera

Multi-anode photomultiplier tubes (MAPMTs) offer high spatial resolution with their small size anodes that may range from 64 to 1024 in number per tube. In order to increase detector size, MAPMT modules can be arranged in arrays and combined in a single modular scintillation camera. However, then the large number of channels that require amplification and digitization become practically not feasible unless signals are combined or reduced in some manner. Conventional approaches use resistive charge division readouts with a centroid algorithm (or a variant of it) for simplicity in the electronic circuitry implementation and fast execution. However, coupling signals from many anodes may cause significant information loss and limit achievable resolution. In this study, a new approach for optimizing readout-electronics design for MAPMTs based on an analysis of information content in the signals is presented. An adaptive read-out scheme to be used with maximum-likelihood estimation methods is proposed. This scheme achieves precision in estimating event parameters that is close to what is achieved by retaining all signals.

Design and Validation of an Adaptive SPECT System: AdaptiSPECT

In order to obtain optimal image quality with respect to a particular task, adaptive imaging systems automatically change their acquisition parameters in response to preliminary data being recorded from the object under study. Currently, the adaptive aspect in Single Photon Emission Computed Tomography (SPECT) is limited to a manual collimator interchange and the choice of detector rotation radius. Furthermore, there is often no optimization of any kind with respect to a certain task. There is thus a need for more versatile SPECT systems that autonomously optimize their acquisition geometry for every task and every patient. Here we describe a pinhole SPECT imager, AdaptiSPECT, which is being developed at the Center for Gamma Ray Imaging (CGRI) to enable adaptive SPECT imaging in a pre-clinical context. Furthermore, ideas for an autonomous adaptation procedure are discussed and some preliminary results are reported upon.

List-mode MLEM Image Reconstruction from 3D ML Position Estimates

Current thick detectors used in medical imaging allow recording many attributes, such as the 3D location of interaction within the scintillation crystal and the amount of energy deposited. An efficient way of dealing with these data is by storing them in list-mode (LM). To reconstruct the data, maximum-likelihood expectation-maximization (MLEM) is efficiently applied to the list-mode data, resulting in the list-mode maximum-likelihood expectation-maximization (LMMLEM) reconstruction algorithm.In this work, we consider a PET system consisting of two thick detectors facing each other. PMT outputs are collected for each coincidence event and are used to perform 3D maximum-likelihood (ML) position estimation of location of interaction. The mathematical properties of the ML estimation allow accurate modeling of the detector blur and provide a theoretical framework for the subsequent estimation step, namely the LMMLEM reconstruction. Indeed, a rigorous statistical model for the detector output can be obtained from calibration data and used in the calculation of the conditional probability density functions for the interaction location estimates.Our implementation of the 3D ML position estimation takes advantage of graphics processing unit (GPU) hardware and permits accurate real-time estimates of position of interaction. The LMMLEM algorithm is then applied to the list of position estimates, and the 3D radiotracer distribution is reconstructed on a voxel grid.

Kinetic characterization of a novel cationic (99m)Tc(I)-tricarbonyl complex, (99m)Tc-15C5-PNP, for myocardial perfusion imaging

BACKGROUND

Intense liver uptake of (99m)Tc-sestamibi (MIBI) often interferes with visualization of myocardial perfusion in the inferior wall of the left ventricle. To develop improved myocardial perfusion agents, crown ether-containing dithiocarbamates and bisphosphines have been introduced in recent years. This study was designed to investigate the myocardial imaging properties and in vivo kinetics of a cationic (99m)Tc(I)-tricarbonyl complex, (99m)Tc-15C5-PNP, in comparison with MIBI.

METHODS

Dynamic cardiac images were acquired for 60 minutes after intravenous tracer injection using a small-animal SPECT system in healthy control rats and rats with myocardial infarcts. Myocardial and liver time-activity curves were generated for radiopharmaceutical kinetic analysis.

RESULTS

Good visualization of the left ventricular wall and perfusion defects could be achieved 20 minutes after (99m)Tc-15C5-PNP administration. (99m)Tc-15C5-PNP images in all hearts with infarcts showed perfusion defects, which were comparable to MIBI images. The kinetic curves plotted from 1 to 60 minutes demonstrated that (99m)Tc-15C5-PNP has a shorter washout half-life (6.4 ± 3.2 vs 124 ± 30.5 minutes, P < .01) in the liver, lower residual liver activity (14.5 ± 10.2% vs 36.5 ± 28.9%, P < .01), and higher heart/liver ratio than MIBI.

CONCLUSIONS

(99m)Tc-15C5-PNP has potential for rapid myocardial perfusion imaging with low liver uptake.

Maximum-Likelihood Calibration of an X-ray Computed Tomography System

We present a maximum-likelihood (ML) method for calibrating the geometrical parameters of an x-ray computed tomography (CT) system. This method makes use of the full image data and not a reduced set of data. This algorithm is particularly useful for CT systems that change their geometry during the CT acquisition, such as an adaptive CT scan. Our ML search method uses a contracting-grid algorithm that does not require initial starting values to perform its estimate, thus avoiding problems associated with choosing initialization values.

DR with a DSLR: Digital Radiography with a Digital Single-Lens Reflex camera

An inexpensive, portable digital radiography (DR) detector system for use in remote regions has been built and evaluated. The system utilizes a large-format digital single-lens reflex (DSLR) camera to capture the image from a standard fluorescent screen. The large sensor area allows relatively small demagnification factors and hence minimizes the light loss. The system has been used for initial phantom tests in urban hospitals and Himalayan clinics in Nepal, and it has been evaluated in the laboratory at the University of Arizona by additional phantom studies. Typical phantom images are presented in this paper, and a simplified discussion of the detective quantum efficiency of the detector is given.

System Integration of FastSPECT III, a Dedicated SPECT Rodent-Brain Imager Based on BazookaSPECT Detector Technology

FastSPECT III is a stationary, single-photon emission computed tomography (SPECT) imager designed specifically for imaging and studying neurological pathologies in rodent brain, including Alzheimer's and Parkinsons's disease. Twenty independent BazookaSPECT [1] gamma-ray detectors acquire projections of a spherical field of view with pinholes selected for desired resolution and sensitivity. Each BazookaSPECT detector comprises a columnar CsI(Tl) scintillator, image-intensifier, optical lens, and fast-frame-rate CCD camera. Data stream back to processing computers via firewire interfaces, and heavy use of graphics processing units (GPUs) ensures that each frame of data is processed in real time to extract the images of individual gamma-ray events. Details of the system design, imaging aperture fabrication methods, and preliminary projection images are presented.

Maximum Likelihood Event Estimation and List-mode Image Reconstruction on GPU Hardware

The scintillation detectors commonly used in SPECT and PET imaging and in Compton cameras require estimation of the position and energy of each gamma ray interaction. Ideally, this process would yield images with no spatial distortion and the best possible spatial resolution. In addition, especially for Compton cameras, the computation must yield the best possible estimate of the energy of each interacting gamma ray. These goals can be achieved by use of maximum-likelihood (ML) estimation of the event parameters, but in the past the search for an ML estimate has not been computationally feasible. Now, however, graphics processing units (GPUs) make it possible to produce optimal, real-time estimates of position and energy, even from scintillation cameras with a large number of photodetectors. In addition, the mathematical properties of ML estimates make them very attractive for use as list entries in list-mode ML image reconstruction. This two-step ML process-using ML estimation once to get the list data and again to reconstruct the object-allows accurate modeling of the detector blur and, potentially, considerable improvement in reconstructed spatial resolution.