This article has Open Peer Review reports available.
Masked-Volume-Wise PCA and "reference Logan" illustrate similar regional differences in kinetic behavior in human brain PET study using [11C]-PIB
© Razifar et al; licensee BioMed Central Ltd. 2009
Received: 02 June 2008
Accepted: 07 January 2009
Published: 07 January 2009
Kinetic modeling using reference Logan is commonly used to analyze data obtained from dynamic Positron Emission Tomography (PET) studies on patients with Alzheimer's disease (AD) and healthy volunteers (HVs) using amyloid imaging agent N-methyl [11C]2-(4'-methylaminophenyl)-6-hydroxy-benzothiazole, [11C]-PIB. The aim of the present study was to explore whether results obtained using the newly introduced method, Masked Volume Wise Principal Component Analysis, MVW-PCA, were similar to the results obtained using reference Logan.
MVW-PCA and reference Logan were performed on dynamic PET images obtained from four Alzheimer's disease (AD) patients on two occasions (baseline and follow-up) and on four healthy volunteers (HVs). Regions of interest (ROIs) of similar sizes were positioned in different parts of the brain in both AD patients and HVs where the difference between AD patients and HVs is largest. Signal-to-noise ratio (SNR) and discrimination power (DP) were calculated for images generated by the different methods and the results were compared both qualitatively and quantitatively.
MVW-PCA generated images that illustrated similar regional binding patterns compared to reference Logan images and with slightly higher quality, enhanced contrast, improved SNR and DP, without being based on modeling assumptions. MVW-PCA also generated additional MVW-PC images by using the whole dataset, which illustrated regions with different and uncorrelated kinetic behaviors of the administered tracer. This additional information might improve the understanding of kinetic behavior of the administered tracer.
MVW-PCA is a potential multivariate method that without modeling assumptions generates high quality images, which illustrated similar regional changes compared to modeling methods such as reference Logan. In addition, MVW-PCA could be used as a new technique, applicable not only on dynamic human brain studies but also on dynamic cardiac studies when using PET.
Imaging modality such as dynamic Positron Emission Tomography, PET, when using amyloid imaging agent N-methyl [11C]2-(4'-methylaminophenyl)-6-hydroxy-benzothiazole, [11C]-PIB , is increasingly used as a diagnostic tool for providing early diagnosis of patients with Alzheimer's disease, AD [2–4]. The tracer principle is based on non-invasive estimation of beta-amyloid (Aβ) plaques in the brain of patients with AD vs. healthy volunteers (HVs). The technique is also valuable when exploring deposition of amyloid in different parts of the brain.
However, images generated using PET contain noise with different magnitude, spatial dependence and correlation, which impairs the visualization and affects the precision in quantification of the generated images. This is due to limitations of the amount of the radioactivity associated with the administered tracer, which is usually short-lived, to different technical limitations and to applied corrections and reconstruction algorithm used [5–7].
To extract quantitative values of a desired physiological information and to study the kinetic behavior of administered tracer in human PET studies, different methods can be used, such as kinetic modeling using e.g. the Patlak [4, 8, 9] and the Logan plots [3, 10]. These models need either the determination of the time course of the unchanged tracer concentration in arterial plasma or the definition of a reference region that is devoid of specific binding. Modeling can also be performed on a pixel-wise basis, generating parametric images. Prior knowledge and modeling experience about the kinetic behavior of the administered tracer in the human brain are needed to find the most appropriate method for a particular tracer. Furthermore, parametric images generally suffer from poor quality and a non-optimized signal-to-noise ratio (SNR). There are also many parameters and variables that need to be correctly considered to obtain proper results when applying these techniques.
In previous studies novel approaches for application of Principal Component Analysis (PCA) on dynamic PET images such as Masked Volume Wise PCA (MVW-PCA), were introduced [5, 6]. Images generated using MVW-PCA contained more detailed anatomical information with higher quality and precision as compared with images generated using other methods. We explored the performance of MVW-PCA as a way to both improve detection and visualization of significant changes in tracer kinetics by extracting different components and to enhance the discrimination between pathological and healthy regions in the brain. Several well-known tracers in clinical and research practice were used such as [11C]-L-deuterium-deprenyl ([11C]-DED), [11C]-GR205171 ([11C]-GLD), [11C]-5-hydroxytryptophan ([11C]-HTP), [11C]L-dihydroxyphenylalanine ([11C]-L-DOPA) and [11C]-PIB .
It was previously shown that MVW-PCA can be used as a multivariate image analysis technique that without modeling assumptions can extract and separate important kinetic information into different component images, MVW-PC images. There are no correlations between the kinetic information separated in these images .
The objectives of the current study were to explore the capabilities of MVW-PCA compared to kinetic modeling using the reference Logan method on the same dataset. The methods were used to extract signals which improve the discrimination between AD patients and HVs in a human brain study using PET and [11C]-PIB. In other words, the aim of the study was to explore whether images generated by MVW-PCA illustrate similar regional retention patterns compared to reference Logan images. Another aim was to study and compare these methods with respect to signal extraction based on calculating SNR and discrimination power (DP).
Dynamic PET data
Data were used from four arbitrarily chosen AD patients and four HVs recruited for different studies at Uppsala Imanet, GE Healthcare in collaboration with several groups such as Karolinska University Hospital Huddinge, Stockholm, Sweden. The four AD patients were also recruited in a follow-up (FU) study where they were scanned a second time at Uppsala Imanet, GE Healthcare . Data from both the baseline and FU studies were used in this study. Four HVs were selected randomly from a study performed at Uppsala Imanet in collaboration with Karolinska University Hospital Huddinge, Stockholm and University of Pittsburgh School of Medicine, Pittsburgh, USA .
Masked Volume Wise Principal Component Analysis
PCA [11, 12] is a data-driven and well-established multivariate analysis technique used to reduce dimensionality of multivariate input datasets such as dynamic PET images. PCA converts the projection of the original images into a new orthogonal coordinate system with lower dimensions in which the new axes explain the variance in the images in decreasing order of importance. This is done by calculating transformation vectors (principal components, i.e. PCs or weight factors), which define the directions of largest variance of the input multi-dimensional dataset in the multidimensional feature space. Each PC is orthogonal to all the others in multi-dimensional space; thus, the first PC (PC1) represents the linear transformation of the original variables which contain the largest variance. The second component, PC2, is the combination that contains the remaining variance as much as possible, orthogonal to the previous one and so on. "PC images" are generated by simply projecting all observations onto the PCs in the new multi-dimensional space. This gives new values along each PC.
If the matrix
X T = [X 1, X 2, X 3, ..., X p ]
where X T is transpose of the input matrix X, which contains pixel values of a single slice from each frame i, (i = (1, 2, 3, ..., p), as a column vector), has a variance-covariance matrix S with eigenvalues
λ = ⌊λ 1, λ 2, λ 3, ..., λ p ⌋
and corresponding eigenvectors
e = ⌊e 1, e 2, e 3, ..., e p ⌋
λ 1 ≥ λ 2 ≥ λ 3 ≥ ... ≥ λ p ≥ 0
and p corresponds to the number of the input column in the matrix X. If q = p in which q refers to number of principal component, then the q th PC is generated by using Eq. 1:
Y q = e'X = e q1 X 1 + e q2 X 2 + e q3 X 3 + ... + e qp X p
The condition Cov(Y q , Y i ) = 0 where i ≠ q is required, i.e. the components are uncorrelated. Principal components should explain the magnitude of variance in decreasing order. Here each element within the eigenvectors is used as weight-factor for creating images.
In a previous work, MVW-PCA was introduced as a new approach of application of PCA on dynamic PET images using various compounds, among others [11C]-PIB . This method is based on using noise prenormalized data that represents whole brain of each time sequence (frame) as a single variable after the background had been removed (masked out) before applying PCA. The result is MVW-PC images, which are separated into different principal components, MVW-PCs. The number of components is equal to the number of PET time frames. In each PET dataset, the whole sequence (all frames) was used to include all information in the dataset. Furthermore, the MVW-PC images can be seen as images representing different kinetic behaviors of the administered tracer and contain more detailed anatomical information with higher precision and quality. They have an improved SNR and visual contrast between the anatomical structures representing both affected and unaffected tissues compared with other methods .
However, in this work, only the first three MVW-PCs were explored since higher components contained only noise.
The software used for the application of MVW-PCA on dynamic PET images was developed in-house by one of the authors (PR) using Matlab 7.2 (The Mathworks, Natick, Massachusetts) with installed statistical and image processing toolboxes.
Reference Logan graphical analysis
The reference Logan graphical analysis is a regression method, which is appropriate for tracers with reversible kinetics. The method describes the kinetics of the tracer in a receptor-containing region and in a reference region devoid of specific binding, such that after a certain time there is a linear relationship between the reference Logan variables. The linear slope gives directly the distribution volume ratio (DVR), which corresponds to the ratio of the distribution volume (DV) of a receptor-containing region to the DV of a reference region. The DVR is widely used as a model parameter in PET studies since it is a linear function of the receptor availability. Parametric maps, in which the DVR is calculated for each pixel in the PET image, were generated using 25–60 min (last six frames) for 4 × 2 AD patients (each at both baseline and follow-up) and four HVs.
where ROI(t) and REF(t) are the radioactivity concentrations in a region of interest and a reference region, respectively. The term int' is the intercept and is the average tissue-to-plasma efflux constant. The slope DVR is obtained from the linear portion of the plot (T > t*), where T and t* are midframe scanning time and equilibrium, respectively. If the ratio ROI(t)/REF(t) is constant, the DVR can be obtained without the use of .
The following four ROIs of similar sizes representing different parts of the brain were drawn, in both AD patients and HVs: frontal cortex sinister (Fsin), frontal cortex dexter (Fdx), parietal cortex sinister (Psin) and parietal cortex dexter (Pdx). According to the study by Klunk  the difference in [11C]-PIB retention between AD patients and HVs is largest in these areas. SNR was calculated for all ROIs in images generated by the different methods and the mean results of the SNR were compared and plotted.
In this study DPs for both images generated by MVW-PCA and reference Logan were calculated using the same ROIs that were used for calculation of SNR (Fsin, Fdx, Psin and Pdx).
In the current study, Graph Pad Prism v. 4.03 (Graph Pad Software, Inc, San Diego, USA) was used to perform all statistical analysis and graphical illustrations.
To allow a visual comparison between the images generated using the two methods, the image color scale minimum- and maximum-levels were set to the image minimum and maximum intensity, respectively. The maximum and the minimum were calculated from a ROI positioned on the whole brain in the image.
Table 1 Statistical comparison between different areas in AD patients ( both in baseline and follow up) and HVs when using different methods.
Paired t-test, two-tailed
Mean of differences
SNR, MVW-PCA vs. Ref. Logan, base line
SNR, MVW-PCA vs. Ref. Logan, follow-up
SNR, MVW-PCA vs. Ref. Logan, HVs
Kinetic modeling can be used to estimate physiological parameters in vivo and to analyze dynamic images in human PET studies using various tracers. However, prior knowledge about kinetic behavior of the tracer under study is needed to select an appropriate modeling method to analyze the data. For instance, it is necessary to know whether the administered tracer binds reversibly or irreversibly to the target system and whether there is any proper reference region devoid of specific binding of administered tracer (if reference region is used as input function). To overcome these drawbacks we have previously introduced a new approach for application of PCA, and more specifically MVW-PCA, on dynamic PET images [5, 6]. The images generated using MVW-PCA, were shown to contain detailed high quality anatomical information. One important feature of MVW-PCA, as well as of other kinetic modeling, is the possibility to separate regions in different compartments, based on their different kinetic properties [5, 6].
The aim of the present study was to explore whether results obtained using MVW-PCA give similar results as using a kinetic modeling method such as reference Logan. We studied the detection and visualization of regional differences in [11C]-PIB retention between affected and unaffected regions in the human brain in patients with AD as compared with HVs. Moreover, we also compared the images generated using MVW-PCA and reference Logan with respect to signal extraction based on calculating SNR and differentiation between AD patients and HVs using DP.
Results obtained from this work show that MVW-PCA generates images which demonstrate similar regional changes in kinetic behavior of administered tracer in human brain in both AD patients and HVs as images generated using reference Logan. Moreover, several MVW-PC images are generated when applying MVW-PCA, images that provide additional information by discriminating and visualizing regions with different types of kinetics in the different MVW-PCs. It is thus possible to obtain separate images of e.g. regions with specific binding vs regions with nonspecific binding or vs. blood compartment. This is a property of great value that is not obtained when using reference Logan. This might improve the understanding of the kinetic behavior of administered tracer by illustrating additional differences in different parts of the brain. The MVW-PCA method is automated; faster and fewer numbers of clinical experiments are required to determine efficiency of the tracer or tracers of interest used in human PET studies. Furthermore, these images show slightly improved SNR and DP compared to the images generated using reference Logan.
So far, MVW-PCA does not yield absolute quantitative binding data, but only relative data. The reference Logan method is however a quantitative method that estimates the distribution volume ratio of the tracer, which is a linear function of receptor availability. MVW-PCA retains the relative proportions between different regions with different amyloid retention with higher qualitative values, but can so far only be used as a qualitative tool. Future studies will be performed to enable quantitative analysis of images generated using MVW-PCA.
To conclude, the results obtained from this study demonstrate the potential of MVW-PCA as a multivariate method that by using the whole dynamic dataset without modeling assumptions generates high quality images. Generated images using MVW-PCA demonstrate similar regional changes in kinetic behavior of administered tracer as compared with images generated using the reference Logan kinetic modeling method with the cerebellum as a reference region. In addition, MVW-PCA generates images with an enhanced contrast and higher SNR and DP. The present evaluation shows that MVW-PCA can be used in the study of beta amyloid distribution in AD patients as compared to HVs. This might improve the clinical interpretation and diagnosis of the PET studies on the human brain when using newly developed tracers or when using existing tracers in new clinical applications where kinetic behavior of employed tracer is not known. In addition, MVW-PCA could be used as a new technique, applicable not only on dynamic human brain studies but also on dynamic human cardiac studies when using PET.
The authors wish to thank Prof. Agneta Nordberg for providing data used in this study. The authors also wish to thank Dr. Sergio Estrada and Dr. Azita Monazzam for beneficial scientific discussions.
Written consent for publication was obtained from the patient or their relative.
- Mathis CA, Wang Y, Holt DP, Huang GF, Debnath ML, Klunk WE: Synthesis and evaluation of 11C-labeled 6-substituted 2-arylbenzothiazoles as amyloid imaging agents. J Med Chem. 2003, 46: 2740-2754. 10.1021/jm030026b.View ArticlePubMedGoogle Scholar
- Klunk WE, Engler H, Nordberg A, Wang Y, Blomqvist B, Holt DP, et al: Imaging Brain Amyloid in Alzheimer's Disease with Pittsburgh Compound-B. Ann Neurol. 2004, 55: 306-319. 10.1002/ana.20009.View ArticlePubMedGoogle Scholar
- Price JC, Klunk WE, Lopresti BJ, Lu X, Hoge JA, Ziolko SK, et al: Kinetic Modeling of Amyloid Binding in Humans using PET Imaging and Pittsburgh Compound-B. J Cereb Blood Flow Metab. 2005, 25: 1528-1547. 10.1038/sj.jcbfm.9600146.View ArticlePubMedGoogle Scholar
- Engler H, Forsberg A, Almkvist O, Blomquist G, Larsson E, Savitcheva I, et al: Two-year follow-up of amyloid deposition in patients with Alzheimer's disease. Brain. 2006, 2856-2866. 10.1093/brain/awl178. Pt10Google Scholar
- Razifar P, Axelsson J, Schneider H, Långström B, Bengtsson E, Bergström M: A new Application of Pre-normalized Principal Component Analysis for Improvement of Image Quality and Clinical Diagnosis in Human Brain PET Studies – Clinical brain studies using [11C]-GR20 [11C]-L-deuterium-deprenyl, [11C]-5-Hydroxy-L-Tryptophan, [11C]-L-DOPA and Pittsburgh Compound-B. Neuroimage. 5171, 33: 588-598. 10.1016/j.neuroimage.2006.05.060.View ArticleGoogle Scholar
- Razifar P, Axelsson J, Schneider H, Långström B, Bengtsson E, Bergström M: Volume-Wise Application of Principal Component Analysis on Masked Dynamic PET Data in Sinogram Domain. IEEE Trans Nucl Med. 2006, 53: 2759-2768. 10.1109/TNS.2006.878008.View ArticleGoogle Scholar
- Zaidi H, Montandon M-L, Meikle SR: Strategies for attenuation compensation in neurological PET studies. Neuroimage. 2007, 34: 518-541. 10.1016/j.neuroimage.2006.10.002.View ArticlePubMedGoogle Scholar
- Patlak CS, Blasberg RG, Fenstermacher JD: Graphical evaluation of blood-to-brain transfer constants from multiple-time uptake data. J Cereb Blood Flow Metab. 1983, 3: 1-7.View ArticlePubMedGoogle Scholar
- Gjedde A, Reith J, Dyve S, Léger G, Guttman M, Diksic M, et al: Dopa decarboxylase activity of the living human brain. Proc Natl Acad Sci Neurobiol. 1991, 88: 2721-2725. 10.1073/pnas.88.7.2721.View ArticleGoogle Scholar
- Logan J: Graphical analysis of PET data applied to reversible and irreversible tracers. Nucl Med Biol. 2000, 27: 661-670. 10.1016/S0969-8051(00)00137-2.View ArticlePubMedGoogle Scholar
- Pearson K: On lines and planes of closest fit to systems of points in space. Philosophical Magazine. 1901, 6: 559-572.View ArticleGoogle Scholar
- Gonzalez RC, Woods RE: Digital image processing. 1992, Addison Wesley Publishing Company, 148-159.Google Scholar
- Razifar P: Novel approaches for application of principal component analysis on dynamic PET images for improvement of image quality and clinical diagnosis. PhD thesis. 2005, Centre for image analysis, Uppsala University, SwedenGoogle Scholar
- Rosnow RL, Rosenthal R: Computing contrasts, effect sizes, and counternulls on other people's published data: General procedures for research consumers. Pyschological Methods. 1996, 1: 331-340. 10.1037/1082-989X.1.4.331.View ArticleGoogle Scholar
- Cohen J: Statistical power analysis for the behavioral sciences. 1998, Hillsdale, NJ: Lawrence Earlbaum Associates, 2Google Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2377/9/2/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.