Biomedical subjects

Richard M Leahy

Publications and source records attributed to Richard M Leahy.

At least 19 recordsLinked to original sources

PubMed2007

Digimouse: a 3D whole body mouse atlas from CT and cryosection data.

We have constructed a three-dimensional (3D) whole body mouse atlas from coregistered x-ray CT and cryosection data of a normal nude male mouse. High quality PET, x-ray CT and cryosection images were acquired post mortem from a single mouse placed in a stereotactic frame with fiducial markers visible in all three modalities. The image data were coregistered to a common coordinate system using the fiducials and resampled to an isotropic 0.1 mm voxel size. Using interactive editing tools we segmented and labelled whole brain, cerebrum, cerebellum, olfactory bulbs, striatum, medulla, masseter muscles, eyes, lachrymal glands, heart, lungs, liver, stomach, spleen, pancreas, adrenal glands, kidneys, testes, bladder, skeleton and skin surface. The final atlas consists of the 3D volume, in which the voxels are labelled to define the anatomical structures listed above, with coregistered PET, x-ray CT and cryosection images. To illustrate use of the atlas we include simulations of 3D bioluminescence and PET image reconstruction. Optical scatter and absorption values are assigned to each organ to simulate realistic photon transport within the animal for bioluminescence imaging. Similarly, 511 keV photon attenuation values are assigned to each structure in the atlas to simulate realistic photon attenuation in PET. The Digimouse atlas and data are available at http://neuroimage.usc.edu/Digimouse.html.

Anatomy, Cross-Sectional↗

PubMed2006

Iterative reconstruction techniques in emission computed tomography.

In emission tomography statistically based iterative methods can improve image quality relative to analytic image reconstruction through more accurate physical and statistical modelling of high-energy photon production and detection processes. Continued exponential improvements in computing power, coupled with the development of fast algorithms, have made routine use of iterative techniques practical, resulting in their increasing popularity in both clinical and research environments. Here we review recent progress in developing statistically based iterative techniques for emission computed tomography. We describe the different formulations of the emission image reconstruction problem and their properties. We then describe the numerical algorithms that are used for optimizing these functions and illustrate their behaviour using small scale simulations.

Algorithms↗

PubMed2006

A hyperspectral fluorescence system for 3D in vivo optical imaging.

In vivo optical instruments designed for small animal imaging generally measure the integrated light intensity across a broad band of wavelengths, or make measurements at a small number of selected wavelengths, and primarily use any spectral information to characterize and remove autofluorescence. We have developed a flexible hyperspectral imaging instrument to explore the use of spectral information to determine the 3D source location for in vivo fluorescence imaging applications. We hypothesize that the spectral distribution of the emitted fluorescence signal can be used to provide additional information to 3D reconstruction algorithms being developed for optical tomography. To test this hypothesis, we have designed and built an in vivo hyperspectral imaging system, which can acquire data from 400 to 1000 nm with 3 nm spectral resolution and which is flexible enough to allow the testing of a wide range of illumination and detection geometries. It also has the capability to generate a surface contour map of the animal for input into the reconstruction process. In this paper, we present the design of the system, demonstrate the depth dependence of the spectral signal in phantoms and show the ability to reconstruct 3D source locations using the spectral data in a simple phantom. We also characterize the basic performance of the imaging system.

Animals↗

PubMed2006

Quantitative evaluation of automated skull-stripping methods applied to contemporary and legacy images: effects of diagnosis, bias correction, and slice location.

Performance of automated methods to isolate brain from nonbrain tissues in magnetic resonance (MR) structural images may be influenced by MR signal inhomogeneities, type of MR image set, regional anatomy, and age and diagnosis of subjects studied. The present study compared the performance of four methods: Brain Extraction Tool (BET; Smith [2002]: Hum Brain Mapp 17:143-155); 3dIntracranial (Ward [1999] Milwaukee: Biophysics Research Institute, Medical College of Wisconsin; in AFNI); a Hybrid Watershed algorithm (HWA, Segonne et al. [2004] Neuroimage 22:1060-1075; in FreeSurfer); and Brain Surface Extractor (BSE, Sandor and Leahy [1997] IEEE Trans Med Imag 16:41-54; Shattuck et al. [2001] Neuroimage 13:856-876) to manually stripped images. The methods were applied to uncorrected and bias-corrected datasets; Legacy and Contemporary T1-weighted image sets; and four diagnostic groups (depressed, Alzheimer's, young and elderly control). To provide a criterion for outcome assessment, two experts manually stripped six sagittal sections for each dataset in locations where brain and nonbrain tissue are difficult to distinguish. Methods were compared on Jaccard similarity coefficients, Hausdorff distances, and an Expectation-Maximization algorithm. Methods tended to perform better on contemporary datasets; bias correction did not significantly improve method performance. Mesial sections were most difficult for all methods. Although AD image sets were most difficult to strip, HWA and BSE were more robust across diagnostic groups compared with 3dIntracranial and BET. With respect to specificity, BSE tended to perform best across all groups, whereas HWA was more sensitive than other methods. The results of this study may direct users towards a method appropriate to their T1-weighted datasets and improve the efficiency of processing for large, multisite neuroimaging studies.

Adult↗

PubMed2006

Generic head models for atlas-based EEG source analysis.

We describe a method for using a generic head model, in the form of an anatomical atlas, to produce EEG source localizations. The atlas is fitted to the subject by a nonrigid warp using a set of surface landmarks. The warped atlas is used to compute a finite element model (FEM) of the forward mapping or lead-fields between neural current generators and the EEG electrodes. These lead-fields are used to localize current sources from the subject's EEG data and the sources are then mapped back to the anatomical atlas. This approach provides a mechanism for comparing source localizations across subjects in an atlas-based coordinate system, which can be used in the large fraction of EEG studies in which MR images are not available. The Montreal brain atlas was used as the reference anatomical atlas and 10 individual MR volumes were used to evaluate the method. The atlas was fitted to each subject's head by a thin-plate-spline (TPS) warp. The spatial locations of a generic 155-electrode configuration were used to constrain the warp. For the purposes of evaluation, dipolar sources were placed on the inner cortical surface in the atlas geometry and transferred to each subject's brain space using a polynomial warp. The parameters of the warp were computed using an intensity-based matching of the atlas and subject brains, thus ensuring that the sources were placed at approximately the same anatomical location in each case. Data were simulated in the subject geometry and a dipole fit was performed on these data using an FEM of the TPS warped atlas. The source positions found in the warped atlas were transferred back to the original atlas and compared to the original position. Sources were simulated at 972 locations evenly distributed over the inner cortical surface of the atlas. The mean error over all 10 subjects was 8.1 mm in the subject space and 15.2 mm in the atlas space. In comparison, using an affine transformation of the electrodes into atlas space and an FEM model generated from the atlas produced mean errors of 22.3 mm in subject space and 19.6 mm in atlas space. With a standard three-shell spherical model the errors were 27.2 mm in the subject space and 34.7 mm when mapped to atlas space.

Brain↗

PubMed2006

Task-based comparison of inverse methods in magnetoencephalography.

Magnetoencephalography (MEG) provides unique insights into the spatio-temporal dynamics of neural activation in the human brain. Unfortunately, the accuracy with which neural sources can be localized is limited by the highly illposed nature of the inverse problem. A large number of inverse methods have been proposed that deal with this illposedness using a range of different modeling and regularization procedures. Here we describe an objective task-based framework for comparing different inverse methods. Using the free-response receiver operating characteristic (FROC) we compare the performance of matched filters, cortically constrained dipole scanning, and minimum norm imaging methods for the task of detecting focal cortical activation. Our results indicate that the scanning methods outperform matched filters and minimum norm imaging for the case of one and two 2 cm2 patches of cortical activity when the dynamics of the two patches are both strongly and weakly correlated and irrespective of the spacing of the two activated regions.

Algorithms↗

PubMed2006

Mean and covariance properties of dynamic PET reconstructions from list-mode data.

We derive computationally efficient methods for the estimation of the mean and variance properties of penalized likelihood dynamic positron emission tomography (PET) images. This allows us to predict the accuracy of reconstructed activity estimates and to compare reconstruction algorithms theoretically. We combine a bin-mode approach in which data is modeled as a collection of independent Poisson random variables at each spatiotemporal bin with the space-time separabilities in the imaging equation and penalties to derive rapidly computable analytic mean and variance approximations. We use these approximations to compare bias/variance properties of our dynamic PET image reconstruction algorithm with those of multiframe static PET reconstructions.

Algorithms↗

PubMed2006

Statistical modeling and reconstruction of randoms precorrected PET data.

Randoms precorrected positron emission tomography (PET) data is formed as the difference of two Poisson random variables. Its exact probability mass function (PMF) is inconvenient for use in likelihood-based iterative image reconstruction as it contains an infinite summation. The shifted Poisson model is a tractable approximation to this PMF but requires that negative values are truncated, resulting in positively biased reconstructions in low count studies. Here we analyze the properties of the exact PMF and propose a simple but accurate approximation that allows negative valued data. We investigate the properties of this approximation and demonstrate its application to penalized maximum likelihood image reconstruction.

Algorithms↗

PubMed2005

Hyperspectral and multispectral bioluminescence optical tomography for small animal imaging.

For bioluminescence imaging studies in small animals, it is important to be able to accurately localize the three-dimensional (3D) distribution of the underlying bioluminescent source. The spectrum of light produced by the source that escapes the subject varies with the depth of the emission source because of the wavelength-dependence of the optical properties of tissue. Consequently, multispectral or hyperspectral data acquisition should help in the 3D localization of deep sources. In this paper, we describe a framework for fully 3D bioluminescence tomographic image acquisition and reconstruction that exploits spectral information. We describe regularized tomographic reconstruction techniques that use semi-infinite slab or FEM-based diffusion approximations of photon transport through turbid media. Singular value decomposition analysis was used for data dimensionality reduction and to illustrate the advantage of using hyperspectral rather than achromatic data. Simulation studies in an atlas-mouse geometry indicated that sub-millimeter resolution may be attainable given accurate knowledge of the optical properties of the animal. A fixed arrangement of mirrors and a single CCD camera were used for simultaneous acquisition of multispectral imaging data over most of the surface of the animal. Phantom studies conducted using this system demonstrated our ability to accurately localize deep point-like sources and show that a resolution of 1.5 to 2.2 mm for depths up to 6 mm can be achieved. We also include an in vivo study of a mouse with a brain tumour expressing firefly luciferase. Co-registration of the reconstructed 3D bioluminescent image with magnetic resonance images indicated good anatomical localization of the tumour.

Algorithms↗

PubMed2005

A comparison of random field theory and permutation methods for the statistical analysis of MEG data.

We describe the use of random field and permutation methods to detect activation in cortically constrained maps of current density computed from MEG data. The methods are applicable to any inverse imaging method that maps event-related MEG to a coregistered cortical surface. These approaches also extend directly to images computed from event-related EEG data. We determine statistical thresholds that control the familywise error rate (FWER) across space or across both space and time. Both random field and permutation methods use the distribution of the maximum statistic under the null hypothesis to find FWER thresholds. The former methods make assumptions on the distribution and smoothness of the data and use approximate analytical solutions, the latter resample the data and rely on empirical distributions. Both methods account for spatial and temporal correlation in the cortical maps. Unlike previous nonparametric work in neuroimaging, we address the problem of nonuniform specificity that can arise without a Gaussianity assumption. We compare and evaluate the methods on simulated data and experimental data from a somatosensory-evoked response study. We find that the random field methods are conservative with or without smoothing, though with a 5 vertex FWHM smoothness, they are close to exact. Our permutation methods demonstrated exact specificity in simulation studies. In real data, the permutation method was not as sensitive as the RF method, although this could be due to violations of the random field theory assumptions.

Data Interpretation, Statistical↗

PubMed2005

Segmentation of skull and scalp in 3-D human MRI using mathematical morphology.

We present a new technique for segmentation of skull and scalp in T(1)-weighted magnetic resonance images (MRIs) of the human head. Our method uses mathematical morphological operations to generate realistic models of the skull, scalp, and brain that are suitable for electroencephalography (EEG) and magnetoencephalography (MEG) source modeling. We first segment the brain using our Brain Surface Extractor algorithm; using this, we can ensure that the brain does not intersect our skull segmentation. We next generate a scalp mask using a combination of thresholding and mathematical morphology. We use the scalp mask in our skull segmentation procedure, as it allows us to automatically exclude background voxels with intensities similar to those of the skull. We find the inner and outer skull boundaries using thresholding and morphological operations. Finally, we mask the results with the scalp and brain volumes to ensure closed and nonintersecting skull boundaries. Visual evaluation indicated accurate segmentations of the cranium at a gross anatomical level (other than small holes in the zygomatic bone in eight subjects) in all 44 MRI volumes processed when run using default settings. In a quantitative comparison with coregistered CT images as a gold standard, MRI skull segmentation accuracy, as measured using the Dice coefficient, was found to be similar to that which would be obtained using CT imagery with a registration error of 2-3 mm.

Algorithms↗

PubMed2004

Optimization and performance evaluation of the microPET II scanner for in vivo small-animal imaging.

MicroPET II is a newly developed PET (positron emission tomography) scanner designed for high-resolution imaging of small animals. It consists of 17,640 LSO crystals each measuring 0.975 x 0.975 x 12.5 mm3, which are arranged in 42 contiguous rings, with 420 crystals per ring. The scanner has an axial field of view (FOV) of 4.9 cm and a transaxial FOV of 8.5 cm. The purpose of this study was to carefully evaluate the performance of the system and to optimize settings for in vivo mouse and rat imaging studies. The volumetric image resolution was found to depend strongly on the reconstruction algorithm employed and averaged 1.1 mm (1.4 microl) across the central 3 cm of the transaxial FOV when using a statistical reconstruction algorithm with accurate system modelling. The sensitivity, scatter fraction and noise-equivalent count (NEC) rate for mouse- and rat-sized phantoms were measured for different energy and timing windows. Mouse imaging was optimized with a wide open energy window (150-750 keV) and a 10 ns timing window, leading to a sensitivity of 3.3% at the centre of the FOV and a peak NEC rate of 235,000 cps for a total activity of 80 MBq (2.2 mCi) in the phantom. Rat imaging, due to the higher scatter fraction, and the activity that lies outside of the field of view, achieved a maximum NEC rate of 24,600 cps for a total activity of 80 MBq (2.2 mCi) in the phantom, with an energy window of 250-750 keV and a 6 ns timing window. The sensitivity at the centre of the FOV for these settings is 2.1%. This work demonstrates that different scanner settings are necessary to optimize the NEC count rate for different-sized animals and different injected doses. Finally, phantom and in vivo animal studies are presented to demonstrate the capabilities of microPET II for small-animal imaging studies.

Animals↗

PubMed2004

Investigation of different transcript quantitation tools for high-throughput mapping of brain gene expression using voxelation.

Voxelation is a new approach for genome scale acquisition of brain gene expression patterns. The method employs high-throughput analysis of spatially registered voxels (cubes) to create multiple volumetric images of brain gene expression, similar to those obtained from biomedical imaging systems. The spatial resolution of voxelation depends on voxel size, with smaller voxels giving higher resolution. An important question is the applicability of different transcript profiling tools for the various levels of resolution that can be employed. Here, we describe the use of three methods to analyze voxel transcript abundance: real-time PCR, microarray analysis and linear amplification coupled with microarrays. We show statistically significant concordance between real-time PCR and microarray analysis for the myelin basic protein gene in human brain specimens at differing levels of spatial resolution. In addition, we also demonstrate the feasibility of using linear amplification coupled with microarray analysis to create voxelation maps from the mouse brain at high resolution, 1 microl. These data indicate the suitability of a number of transcript profiling tools for various levels of spatial resolution in voxelation.

Brain↗

PubMed2004

Accurate estimation of the Fisher information matrix for the PET image reconstruction problem.

The Fisher information matrix (FIM) plays a key role in the analysis and applications of statistical image reconstruction methods based on Poisson data models. The elements of the FIM are a function of the reciprocal of the mean values of sinogram elements. Conventional plug-in FIM estimation methods do not work well at low counts, where the FIM estimate is highly sensitive to the reciprocal mean estimates at individual detector pairs. A generalized error look-up table (GELT) method is developed to estimate the reciprocal of the mean of the sinogram data. This approach is also extended to randoms precorrected data. Based on these techniques, an accurate FIM estimate is obtained for both Poisson and randoms precorrected data. As an application, the new GELT method is used to improve resolution uniformity and achieve near-uniform image resolution in low count situations.

Algorithms↗

PubMed2003

High-resolution voxelation mapping of human and rodent brain gene expression.

Voxelation allows high-throughput acquisition of multiple volumetric images of brain gene expression, similar to those obtained from biomedical imaging systems. To obtain these images, the method employs analysis of spatially registered voxels (cubes). For creation of high-resolution maps using voxelation, relatively small voxel sizes are necessary and instruments will be required for semiautomated harvesting of such voxels. Here, we describe two devices that allow spatially registered harvesting of voxels from the human and rodent brain, giving linear resolutions of 3.3 and 1 mm, respectively. Gene expression patterns obtained using these devices showed good agreement with known expression patterns. The voxelation instruments and their future iterations represent a valuable approach to the genome scale acquisition of gene expression patterns in the human and rodent brain.

Animals↗

PubMed2003

Genes regulated by learning in the hippocampus.

The enduring changes in long-term memory probably depend on regulation of gene expression in the hippocampus. To seek genes regulated by learning, we used microarray technology to compare hippocampal gene expression in mice undergoing training in the Morris water maze and control mice forced to swim for the same period in the absence of a hidden platform. ANOVA was employed to prioritize genes for further study, and three genes were confirmed by real-time PCR as being regulated during learning. One of the genes was the alpha subunit of the platelet-derived growth factor receptor (Pdgfra); another showed homology to DnaJ and cAMP response element-binding protein 2 (CREB2); and a third was novel. These genes may provide useful insights into the molecular mechanisms of hippocampal learning.

Activating Transcription Factor 4↗

PubMed2003

Spatiotemporal localization of significant activation in MEG using permutation tests.

We describe the use of non-parametric permutation tests to detect activation in cortically-constrained maps of current density computed from MEG data. The methods are applicable to any inverse imaging method that maps event-related MEG to a coregistered cortical surface. To determine an appropriate threshold to apply to statistics computed from these maps, it is important to control for the multiple testing problem associated with testing 10's of thousands of hypotheses (one per surface element). By randomly permuting pre- and post-stimuius data from the collection of individual epochs in an event related study, we develop thresholds that control the familywise (type 1) error rate. These thresholds are based on the distribution of the maximum intensity, which implicitly accounts for spatial and temporal correlation in the cortical maps. We demonstrate the method in application to simulated data and experimental data from a somatosensory evoked response study.

Algorithms↗

PubMed2003

Error-correcting microarray design.

We describe a microarray design based on the concept of error-correcting codes from digital communication theory. Currently, microarrays are unable to efficiently deal with "drop-outs," when one or more spots on the array are corrupted. The resulting information loss may lead to decoding errors in which no quantitation of expression can be extracted for the corresponding genes. This issue is expected to become increasingly problematic as the number of spots on microarrays expands to accommodate the entire genome. The error-correcting approach employs multiplexing (encoding) of more than one gene onto each spot to efficiently provide robustness to drop-outs in the array. Decoding then allows fault-tolerant recovery of the expression information from individual genes. The error-correcting method is general and may have important implications for future array designs in research and diagnostics.

Data Interpretation, Statistical↗