Partial volume correction of brain PET studies using iterative deconvolution in combination with HYPR denoising

Application of the Van Cittert Algorithm for Deconvolving Loss Features in X-ray Photoelectron Spectroscopy Spectra

The convolution of two physical entities, denoted as f and g, delineates the manner in which one entity undergoes modification in response to the other. This transformative process is mathematically represented by the expression f ⨂ g, symbolizing the convolution of the two entities in a resultant function h. Frequently, it becomes imperative to comprehend the magnitude of the induced modifications. From the derived function h, a crucial step involves the separation of the two original signals, a process commonly referred to as deconvolution. Various techniques have been proposed to facilitate the calculation of the deconvolution, with one notable approach originating in 1931 by van Cittert. The algorithm, based on an iterative method, has been scrutinized over time, notably by Bracewell and, more recently, by Jansson. This work represents the current state-of-the-art, focusing specifically on the analysis of Auger spectra obtained through XPS. Emphasis is placed on delineating the procedural aspects of the analysis, and the algorithm utilized in the open-source software RxpsG is comprehensively described.

Longitudinal Tau PET Using 18F-Flortaucipir: The Effect of Relative Cerebral Blood Flow on Quantitative and Semiquantitative Parameters

Semiquantitative PET measures such as SUV ratio (SUVr) have several advantages over quantitative measures, such as practical applicability and relative computational simplicity. However, SUVr may potentially be affected by changes in blood flow, whereas quantitative measures such as nondisplaceable binding potential (BP_ND) are not. For ¹⁸F-flortaucipir PET, the sensitivity of SUVr for changes in blood flow is currently unknown. Therefore, we compared semiquantitative (SUVr) and quantitative (BP_ND) parameters of longitudinal ¹⁸F-flortaucipir PET scans and assessed their vulnerability to changes in blood flow. Methods: Subjects with subjective cognitive decline (n = 38) and Alzheimer disease patients (n = 24) underwent baseline and 2-y follow-up dynamic ¹⁸F-flortaucipir PET scans. BP_ND and relative tracer delivery were estimated using receptor parametric mapping, and SUVr at 80-100 min was calculated. Regional SUVrs were compared with corresponding distribution volume ratio (BP_ND + 1) using paired t tests. Additionally, simulations were performed to model effects of larger flow changes in different binding categories. Results: Results in subjective cognitive decline and Alzheimer disease showed only minor differences between SUVr and BP_ND changes over time. Relative tracer delivery changes were small in all groups. Simulations illustrated a variable bias for SUVr depending on the amount of binding. Conclusion: SUVr provided an accurate estimate of changes in specific binding for ¹⁸F-flortaucipir over a 2-y follow-up during which changes in flow were small. Notwithstanding, simulations showed that large(r) flow changes may affect ¹⁸F-flortaucipir SUVr. Given that it is currently unknown to what order of magnitude pharmacotherapeutic interventions may induce changes in cerebral blood flow, caution may be warranted when changes in flow are potentially large(r), as in clinical trials.

Impact of Aggregation Methods for Texture Features on Their Robustness Performance: Application to Nasopharyngeal 18F-FDG PET/CT

PURPOSE

This study aims to investigate the impact of aggregation methods used for the generation of texture features on their robustness of nasopharyngeal carcinoma (NPC) based on ¹⁸F-FDG PET/CT images.

METHODS

128 NPC patients were enrolled and 95 texture features were extracted for each patient including six feature families under different aggregation methods. For GLCM and GLRLM features, six aggregation methods were considered. For GLSZM, GLDZM, NGTDM and NGLDM features, three aggregation methods were considered. The robustness of the features affected by aggregation methods was assessed by the pair-wise intra-class correlation coefficient (ICC). Furthermore, the effects of discretization and partial volume correction (PVC) on the percent of ICC categories of all texture features were evaluated by overall ICC instead of the pair-wise ICC.

RESULTS

There were 12 features with excellent pair-wise ICCs varying aggregation methods, namely joint average, sum average, autocorrelation, long run emphasis, high grey level run emphasis, short run high grey level emphasis, long run high grey level emphasis, run length variance, SZM high grey level emphasis, DZM high grey level emphasis, high grey level count emphasis and dependence count percentage. For GLCM and GLRLM features, 19/25 and 14/16 features showed excellent pair-wise ICCs varying aggregation methods (averaged and merged) on the same dimensional features (2D, 2.5D or 3D). Different discretization levels and partial volume corrections lead to consistent robustness of textural features affected by aggregation methods.

CONCLUSION

Different dimensional features with the same aggregation methods showed worse robustness compared with the same dimensional features with different aggregation methods. Different discretization levels and PVC algorithms had a negligible effect on the percent of ICC categories of all texture features.

Cognitive performance in multiple sclerosis: what is the role of the gamma-aminobutyric acid system?

Cognitive impairment occurs in 40-65% of persons with multiple sclerosis and may be related to alterations in glutamatergic and GABAergic neurotransmission. Therefore, the aim of this study was to determine how glutamatergic and GABAergic changes relate to cognitive functioning in multiple sclerosis in vivo. Sixty persons with multiple sclerosis (mean age 45.5 ± 9.6 years, 48 females, 51 relapsing-remitting multiple sclerosis) and 22 age-matched healthy controls (45.6 ± 22.0 years, 17 females) underwent neuropsychological testing and MRI. Persons with multiple sclerosis were classified as cognitively impaired when scoring at least 1.5 standard deviations below normative scores on ≥30% of tests. Glutamate and GABA concentrations were determined in the right hippocampus and bilateral thalamus using magnetic resonance spectroscopy. GABA-receptor density was assessed using quantitative [¹¹C]flumazenil positron emission tomography in a subset of participants. Positron emission tomography outcome measures were the influx rate constant (a measure predominantly reflecting perfusion) and volume of distribution, which is a measure of GABA-receptor density. Twenty persons with multiple sclerosis (33%) fulfilled the criteria for cognitive impairment. No differences were observed in glutamate or GABA concentrations between persons with multiple sclerosis and healthy controls, or between cognitively preserved, impaired and healthy control groups. Twenty-two persons with multiple sclerosis (12 cognitively preserved and 10 impaired) and 10 healthy controls successfully underwent [¹¹C]flumazenil positron emission tomography. Persons with multiple sclerosis showed a lower influx rate constant in the thalamus, indicating lower perfusion. For the volume of distribution, persons with multiple sclerosis showed higher values than controls in deep grey matter, reflecting increased GABA-receptor density. When comparing cognitively impaired and preserved patients to controls, the preserved group showed a significantly higher volume of distribution in cortical and deep grey matter and hippocampus. Positive correlations were observed between both positron emission tomography measures and information processing speed in the multiple sclerosis group only. Whereas concentrations of glutamate and GABA did not differ between multiple sclerosis and control nor between cognitively impaired, preserved and control groups, increased GABA-receptor density was observed in preserved persons with multiple sclerosis that was not seen in cognitively impaired patients. In addition, GABA-receptor density correlated to cognition, in particular with information processing speed. This could indicate that GABA-receptor density is upregulated in the cognitively preserved phase of multiple sclerosis as a means to regulate neurotransmission and potentially preserve cognitive functioning.

Implementing the Point Spread Function Deconvolution for Better Molecular Characterization of Newly Diagnosed Gliomas: A Dynamic 18F-FDOPA PET Radiomics Study

Purpose: This study aims to investigate the effects of applying the point spread function deconvolution (PSFd) to the radiomics analysis of dynamic L-3,4-dihydroxy-6-[18F]-fluoro-phenyl-alanine (18F-FDOPA) positron emission tomography (PET) images, to non-invasively identify isocitrate dehydrogenase (IDH) mutated and/or 1p/19q codeleted gliomas. Methods: Fifty-seven newly diagnosed glioma patients underwent dynamic 18F-FDOPA imaging on the same digital PET system. All images were reconstructed with and without PSFd. An L1-penalized (Lasso) logistic regression model, with 5-fold cross-validation and 20 repetitions, was trained with radiomics features extracted from the static tumor-to-background-ratio (TBR) and dynamic time-to-peak (TTP) parametric images, as well as a combination of both. Feature importance was assessed using Shapley additive explanation values. Results: The PSFd significantly modified 95% of TBR, but only 79% of TTP radiomics features. Applying the PSFd significantly improved the ability to identify IDH-mutated and/or 1p/19q codeleted gliomas, compared to PET images not processed with PSFd, with respective areas under the curve of 0.83 versus 0.79 and 0.75 versus 0.68 for a combination of static and dynamic radiomics features (p < 0.001). Without the PSFd, four and eight radiomics features contributed to 50% of the model for detecting IDH-mutated and/or 1p/19q codeleted gliomas, respectively. Application of the PSFd reduced this to three and seven contributive radiomics features. Conclusion: Application of the PSFd to dynamic 18F-FDOPA PET imaging significantly improves the detection of molecular parameters in newly diagnosed gliomas, most notably by modifying TBR radiomics features.

PET image enhancement using artificial intelligence for better characterization of epilepsy lesions

INTRODUCTION

[18F]fluorodeoxyglucose ([18F]FDG) brain PET is used clinically to detect small areas of decreased uptake associated with epileptogenic lesions, e.g., Focal Cortical Dysplasias (FCD) but its performance is limited due to spatial resolution and low contrast. We aimed to develop a deep learning-based PET image enhancement method using simulated PET to improve lesion visualization.

METHODS

We created 210 numerical brain phantoms (MRI segmented into 9 regions) and assigned 10 different plausible activity values (e.g., GM/WM ratios) resulting in 2100 ground truth high quality (GT-HQ) PET phantoms. With a validated Monte-Carlo PET simulator, we then created 2100 simulated standard quality (S-SQ) [18F]FDG scans. We trained a ResNet on 80% of this dataset (10% used for validation) to learn the mapping between S-SQ and GT-HQ PET, outputting a predicted HQ (P-HQ) PET. For the remaining 10%, we assessed Peak Signal-to-Noise Ratio (PSNR), Structural Similarity Index Measure (SSIM), and Root Mean Squared Error (RMSE) against GT-HQ PET. For GM and WM, we computed recovery coefficients (RC) and coefficient of variation (COV). We also created lesioned GT-HQ phantoms, S-SQ PET and P-HQ PET with simulated small hypometabolic lesions characteristic of FCDs. We evaluated lesion detectability on S-SQ and P-HQ PET both visually and measuring the Relative Lesion Activity (RLA, measured activity in the reduced-activity ROI over the standard-activity ROI). Lastly, we applied our previously trained ResNet on 10 clinical epilepsy PETs to predict the corresponding HQ-PET and assessed image quality and confidence metrics.

RESULTS

Compared to S-SQ PET, P-HQ PET improved PNSR, SSIM and RMSE; significatively improved GM RCs (from 0.29 ± 0.03 to 0.79 ± 0.04) and WM RCs (from 0.49 ± 0.03 to 1 ± 0.05); mean COVs were not statistically different. Visual lesion detection improved from 38 to 75%, with average RLA decreasing from 0.83 ± 0.08 to 0.67 ± 0.14. Visual quality of P-HQ clinical PET improved as well as reader confidence.

CONCLUSION

P-HQ PET showed improved image quality compared to S-SQ PET across several objective quantitative metrics and increased detectability of simulated lesions. In addition, the model generalized to clinical data. Further evaluation is required to study generalization of our method and to assess clinical performance in larger cohorts.

Monash DaCRA fPET-fMRI: A dataset for comparison of radiotracer administration for high temporal resolution functional FDG-PET

BACKGROUND

"Functional" [18F]-fluorodeoxyglucose positron emission tomography (FDG-fPET) is a new approach for measuring glucose uptake in the human brain. The goal of FDG-fPET is to maintain a constant plasma supply of radioactive FDG in order to track, with high temporal resolution, the dynamic uptake of glucose during neuronal activity that occurs in response to a task or at rest. FDG-fPET has most often been applied in simultaneous BOLD-fMRI/FDG-fPET (blood oxygenation level-dependent functional MRI fluorodeoxyglucose functional positron emission tomography) imaging. BOLD-fMRI/FDG-fPET provides the capability to image the 2 primary sources of energetic dynamics in the brain, the cerebrovascular haemodynamic response and cerebral glucose uptake.

FINDINGS

In this Data Note, we describe an open access dataset, Monash DaCRA fPET-fMRI, which contrasts 3 radiotracer administration protocols for FDG-fPET: bolus, constant infusion, and hybrid bolus/infusion. Participants (n = 5 in each group) were randomly assigned to each radiotracer administration protocol and underwent simultaneous BOLD-fMRI/FDG-fPET scanning while viewing a flickering checkerboard. The bolus group received the full FDG dose in a standard bolus administration, the infusion group received the full FDG dose as a slow infusion over the duration of the scan, and the bolus-infusion group received 50% of the FDG dose as bolus and 50% as constant infusion. We validate the dataset by contrasting plasma radioactivity, grey matter mean uptake, and task-related activity in the visual cortex.

CONCLUSIONS

The Monash DaCRA fPET-fMRI dataset provides significant reuse value for researchers interested in the comparison of signal dynamics in fPET, and its relationship with fMRI task-evoked activity.

Deep learning-based image reconstruction and post-processing methods in positron emission tomography for low-dose imaging and resolution enhancement

Image processing plays a crucial role in maximising diagnostic quality of positron emission tomography (PET) images. Recently, deep learning methods developed across many fields have shown tremendous potential when applied to medical image enhancement, resulting in a rich and rapidly advancing literature surrounding this subject. This review encapsulates methods for integrating deep learning into PET image reconstruction and post-processing for low-dose imaging and resolution enhancement. A brief introduction to conventional image processing techniques in PET is firstly presented. We then review methods which integrate deep learning into the image reconstruction framework as either deep learning-based regularisation or as a fully data-driven mapping from measured signal to images. Deep learning-based post-processing methods for low-dose imaging, temporal resolution enhancement and spatial resolution enhancement are also reviewed. Finally, the challenges associated with applying deep learning to enhance PET images in the clinical setting are discussed and future research directions to address these challenges are presented.

Genetically identical twins show comparable tau PET load and spatial distribution

Tau accumulation starts during the preclinical phase of Alzheimer’s disease and is closely associated with cognitive decline. For preventive purposes, it is important to identify factors associated with tau accumulation and spread. Studying genetically identical twin-pairs may give insight into genetic and environmental contributions to tau pathology, as similarities in identical twin-pairs largely result from genetic factors, while differences in identical twin-pairs can largely be attributed to non-shared, environmental factors. This study aimed to examine similarities and dissimilarities in a cohort of genetically identical older twin-pairs in (i) tau load; and (ii) spatial distribution of tau, measured with ¹⁸F-flortaucipir PET.

We selected 78 genetically identical twins (39 pairs; average age 73 ± 6 years), enriched for amyloid-β pathology and APOE ε4 carriership, who underwent dynamic ¹⁸F-flortaucipir PET. We extracted binding potentials (BP_ND) in entorhinal, temporal, widespread neocortical and global regions, and examined within-pair similarities in BP_ND using age and sex corrected intra-class correlations. Furthermore, we tested whether twin-pairs showed a more similar spatial ¹⁸F-flortaucipir distribution compared to non-twin pairs, and whether the participant’s co-twin could be identified solely based on the spatial ¹⁸F-flortaucipir distribution. Last, we explored whether environmental (e.g. physical activity, obesity) factors could explain observed differences in twins of a pair in ¹⁸F-flortaucipir BP_ND.

On visual inspection, Alzheimer’s disease-like ¹⁸F-flortaucipir PET patterns were observed, and although we mainly identified similarities in twin-pairs, some pairs showed strong dissimilarities. ¹⁸F-flortaucipir BP_ND was correlated in twins in the entorhinal (r = 0.40; P = 0.01), neocortical (r = 0.59; P < 0.01) and global (r = 0.56; P < 0.01) regions, but not in the temporal region (r = 0.20; P = 0.10). The ¹⁸F-flortaucipir distribution pattern was significantly more similar between twins of the same pair [mean r = 0.27; standard deviation (SD) = 0.09] than between non-twin pairings of participants (mean r = 0.01; SD = 0.10) (P < 0.01), also after correcting for proxies of off-target binding. Based on the spatial ¹⁸F-flortaucipir distribution, we could identify with an accuracy of 86% which twins belonged to the same pair. Finally, within-pair differences in ¹⁸F-flortaucipir BP_ND were associated with within-pair differences in depressive symptoms (0.37 < β < 0.56), physical activity (−0.41 < β < −0.42) and social activity (−0.32 < β < −0.36) (all P < 0.05).

Overall, identical twin-pairs were comparable in tau load and spatial distribution, highlighting the important role of genetic factors in the accumulation and spreading of tau pathology. Considering also the presence of dissimilarities in tau pathology in identical twin-pairs, our results additionally support a role for (potentially modifiable) environmental factors in the onset of Alzheimer’s disease pathological processes, which may be of interest for future prevention strategies.

Differential associations between neocortical tau pathology and blood flow with cognitive deficits in early-onset vs late-onset Alzheimer's disease

Purpose

Early-onset Alzheimer’s disease (EOAD) and late-onset Alzheimer’s disease (LOAD) differ in neuropathological burden and type of cognitive deficits. Assessing tau pathology and relative cerebral blood flow (rCBF) measured with [¹⁸F]flortaucipir PET in relation to cognition may help explain these differences between EOAD and LOAD.

Methods

Seventy-nine amyloid-positive individuals with a clinical diagnosis of AD (EOAD: n = 35, age-at-PET = 59 ± 5, MMSE = 23 ± 4; LOAD: n = 44, age-at-PET = 71 ± 5, MMSE = 23 ± 4) underwent a 130-min dynamic [¹⁸F]flortaucipir PET scan and extensive neuropsychological assessment. We extracted binding potentials (BP_ND) and R₁ (proxy of rCBF) from parametric images using receptor parametric mapping, in medial and lateral temporal, parietal, occipital, and frontal regions-of-interest and used nine neuropsychological tests covering memory, attention, language, and executive functioning. We first examined differences between EOAD and LOAD in BP_ND or R₁ using ANOVA (region-of-interest analysis) and voxel-wise contrasts. Next, we performed linear regression models to test for potential interaction effects between age-at-onset and BP_ND/R₁ on cognition.

Results

Both region-of-interest and voxel-wise contrasts showed higher [¹⁸F]flortaucipir BP_ND values across all neocortical regions in EOAD. By contrast, LOAD patients had lower R₁ values (indicative of more reduced rCBF) in medial temporal regions. For both tau and flow in lateral temporal, and occipitoparietal regions, associations with cognitive impairment were stronger in EOAD than in LOAD (EOAD BP_ND − 0.76 ≤ stβ ≤  − 0.48 vs LOAD − 0.18 ≤ stβ ≤  − 0.02; EOAD R₁ 0.37 ≤ stβ ≤ 0.84 vs LOAD − 0.25 ≤ stβ ≤ 0.16).

Conclusions

Compared to LOAD, the degree of lateral temporal and occipitoparietal tau pathology and relative cerebral blood-flow is more strongly associated with cognition in EOAD.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00259-021-05669-6.

Voxel-based partial volume correction of PET images via subtle MRI guided non-local means regularization

PURPOSE

Positron emission tomography (PET) images tend to be significantly degraded by the partial volume effect (PVE) resulting from the limited spatial resolution of the reconstructed images. Our purpose is to propose a partial volume correction (PVC) method to tackle this issue.

METHODS

In the present work, we explore a voxel-based PVC method under the least squares framework (LS) employing anatomical non-local means (NLMA) regularization. The well-known non-local means (NLM) filter utilizes the high degree of information redundancy that typically exists in images, and is typically used to directly reduce image noise by replacing each voxel intensity with a weighted average of its non-local neighbors. Here we explore NLM as a regularization term within iterative-deconvolution model to perform PVC. Further, an anatomical-guided version of NLM was proposed that incorporates MRI information into NLM to improve resolution and suppress image noise. The proposed approach makes subtle usage of the accompanying MRI information to define a more appropriate search space within the prior model. To optimize the regularized LS objective function, we used the Gauss-Seidel (GS) algorithm with the one-step-late (OSL) technique.

RESULTS

After the import of NLMA, the visual and quality results are all improved. With a visual check, we notice that NLMA reduce the noise compared to other PVC methods. This is also validated in bias-noise curve compared to non-MRI-guided PVC framework. We can see NLMA gives better bias-noise trade-off compared to other PVC methods.

CONCLUSIONS

Our efforts were evaluated in the base of amyloid brain PET imaging using the BrainWeb phantom and in vivo human data. We also compared our method with other PVC methods. Overall, we demonstrated the value of introducing subtle MRI-guidance in the regularization process, the proposed NLMA method resulting in promising visual as well as quantitative performance improvements.

Deconvolution-based partial volume correction of PET images with parallel level set regularization

The partial volume effect (PVE), caused by the limited spatial resolution of positron emission tomography (PET), degrades images both qualitatively and quantitatively. Anatomical information provided by magnetic resonance (MR) images has the potential to play an important role in partial volume correction (PVC) methods. Post-reconstruction MR-guided PVC methods typically use segmented MR tissue maps, and further, assume that PET activity distribution is uniform in each region, imposing considerable constraints through anatomical guidance. In this work, we present a post-reconstruction PVC method based on deconvolution with parallel level set (PLS) regularization. We frame the problem as an iterative deconvolution task with PLS regularization that incorporates anatomical information without requiring MR segmentation or assuming uniformity of PET distributions within regions. An efficient algorithm for non-smooth optimization of the objective function (invoking split Bregman framework) is developed so that the proposed method can be feasibly applied to 3D images and produces sharper images compared to PLS method with smooth optimization. The proposed method was evaluated together with several other PVC methods using both realistic simulation experiments based on the BrainWeb phantom as well asin vivohuman data. Our proposed method showed enhanced quantitative performance when realistic MR guidance was provided. Further, the proposed method is able to reduce image noise while preserving structure details onin vivohuman data, and shows the potential to better differentiate amyloid positive and amyloid negative scans. Overall, our results demonstrate promise to provide superior performance in clinical imaging scenarios.

SimPET-An open online platform for the Monte Carlo simulation of realistic brain PET data. Validation for 18 F-FDG scans

Purpose

SimPET (www.sim‐pet.org) is a free cloud‐based platform for the generation of realistic brain positron emission tomography (PET) data. In this work, we introduce the key features of the platform. In addition, we validate the platform by performing a comparison between simulated healthy brain FDG‐PET images and real healthy subject data for three commercial scanners (GE Advance NXi, GE Discovery ST, and Siemens Biograph mCT).

Methods

The platform provides a graphical user interface to a set of automatic scripts taking care of the code execution for the phantom generation, simulation (SimSET), and tomographic image reconstruction (STIR). We characterize the performance using activity and attenuation maps derived from PET/CT and MRI data of 25 healthy subjects acquired with a GE Discovery ST. We then use the created maps to generate synthetic data for the GE Discovery ST, the GE Advance NXi, and the Siemens Biograph mCT. The validation was carried out by evaluating Bland‐Altman differences between real and simulated images for each scanner. In addition, SPM voxel‐wise comparison was performed to highlight regional differences. Examples for amyloid PET and for the generation of ground‐truth pathological patients are included.

Results

The platform can be efficiently used for generating realistic simulated FDG‐PET images in a reasonable amount of time. The validation showed small differences between SimPET and acquired FDG‐PET images, with errors below 10% for 98.09% (GE Discovery ST), 95.09% (GE Advance NXi), and 91.35% (Siemens Biograph mCT) of the voxels. Nevertheless, our SPM analysis showed significant regional differences between the simulated images and real healthy patients, and thus, the use of the platform for converting control subject databases between different scanners requires further investigation.

Conclusions

The presented platform can potentially allow scientists in clinical and research settings to perform MC simulation experiments without the need for high‐end hardware or advanced computing knowledge and in a reasonable amount of time.

In vivo tau pathology is associated with synaptic loss and altered synaptic function

BACKGROUND

The mechanism of synaptic loss in Alzheimer's disease is poorly understood and may be associated with tau pathology. In this combined positron emission tomography (PET) and magnetoencephalography (MEG) study, we aimed to investigate spatial associations between regional tau pathology ([18F]flortaucipir PET), synaptic density (synaptic vesicle 2A [11C]UCB-J PET) and synaptic function (MEG) in Alzheimer's disease.

METHODS

Seven amyloid-positive Alzheimer's disease subjects from the Amsterdam Dementia Cohort underwent dynamic 130-min [18F]flortaucipir PET, dynamic 60-min [11C]UCB-J PET with arterial sampling and 2 × 5-min resting-state MEG measurement. [18F]flortaucipir- and [11C]UCB-J-specific binding (binding potential, BPND) and MEG spectral measures (relative delta, theta and alpha power; broadband power; and peak frequency) were assessed in cortical brain regions of interest. Associations between regional [18F]flortaucipir BPND, [11C]UCB-J BPND and MEG spectral measures were assessed using Spearman correlations and generalized estimating equation models.

RESULTS

Across subjects, higher regional [18F]flortaucipir uptake was associated with lower [11C]UCB-J uptake. Within subjects, the association between [11C]UCB-J and [18F]flortaucipir depended on within-subject neocortical tau load; negative associations were observed when neocortical tau load was high, gradually changing into opposite patterns with decreasing neocortical tau burden. Both higher [18F]flortaucipir and lower [11C]UCB-J uptake were associated with altered synaptic function, indicative of slowing of oscillatory activity, most pronounced in the occipital lobe.

CONCLUSIONS

These results indicate that in Alzheimer's disease, tau pathology is closely associated with reduced synaptic density and synaptic dysfunction.

Comparison of MRI-based and PET-based image pre-processing for quantification of 11C-PBB3 uptake in human brain

PURPOSE

Quantification of tau load using ¹¹C-PBB3-PET has the potential to improve diagnosis of neurodegenerative diseases. Although MRI-based pre-processing is used as a reference method, not all patients have MRI. The feasibility of a PET-based pre-processing for the quantification of ¹¹C-PBB3 tracer was evaluated and compared with the MRI-based method.

MATERIALS AND METHODS

Fourteen patients with decreased recent memory were examined with ¹¹C-PBB3-PET and MRI. The PET scans were visually assessed and rated as either PBB3(+) or PBB3(-). The image processing based on the PET-based method was validated against the MRI-based approach. The regional uptakes were quantified using the Mesial-temporal/Temporoparietal/Rest of neocortex (MeTeR) regions. SUVR values were calculated by normalizing to the cerebellar reference region to compare both methods within the patient groups.

RESULTS

Significant correlations were observed between the SUVRs of the MRI-based and the PET-based methods in the MeTeR regions (r_Me=0.91; r_Te=0.98; r_R=0.96; p<0.0001). However, the Bland-Altman plot showed a significant bias between both methods in the subcortical Me region (bias: -0.041; 95% CI: -0.061 to -0.024; p=0.003). As in the MRI-based method, the ¹¹C-PBB3 uptake obtained with the PET-based method was higher for the PBB3(+) group in each of the cortical regions and for the whole brain than for the PBB3(-) group (PET-based_Global: 1.11 vs. 0.96; Cliff's Delta (d)=0.68; p=0.04; MRI-based_Global: 1.11 vs. 0.97; d=0.70; p=0.03). To differentiate between positive and negative scans, Youden's index estimated the best cut-off of 0.99 from the ROC curve with good accuracy (AUC: 0.88±0.10; 95% CI: 0.67-1.00) and the same sensitivity (83%) and specificity (88%) for both methods.

CONCLUSION

The PET-based pre-processing method developed to quantify the tau burden with ¹¹C-PBB3 provided comparable SUVR values and effect sizes as the MRI-based reference method. Furthermore, both methods have a comparable discrimination accuracy between PBB3(+) and PBB3(-) groups as assessed by visual rating. Therefore, the presented PET-based method can be used for clinical diagnosis if no MRI image is available.

Tau pathology and relative cerebral blood flow are independently associated with cognition in Alzheimer's disease

Purpose

We aimed to investigate associations between tau pathology and relative cerebral blood flow (rCBF), and their relationship with cognition in Alzheimer’s disease (AD), by using a single dynamic [¹⁸F]flortaucipir positron emission tomography (PET) scan.

Methods

Seventy-one subjects with AD (66 ± 8 years, mini-mental state examination (MMSE) 23 ± 4) underwent a dynamic 130-min [¹⁸F]flortaucipir PET scan. Cognitive assessment consisted of composite scores of four cognitive domains. For tau pathology and rCBF, receptor parametric mapping (cerebellar gray matter reference region) was used to create uncorrected and partial volume-corrected parametric images of non-displaceable binding potential (BP_ND) and R₁, respectively. (Voxel-wise) linear regressions were used to investigate associations between BP_ND and/or R₁ and cognition_.

Results

Higher [¹⁸F]flortaucipir BP_ND was associated with lower R₁ in the lateral temporal, parietal and occipital regions. Higher medial temporal BP_ND was associated with worse memory, and higher lateral temporal BP_ND with worse executive functioning and language. Higher parietal BP_ND was associated with worse executive functioning, language and attention, and higher occipital BP_ND with lower cognitive scores across all domains. Higher frontal BP_ND was associated with worse executive function and attention. For [¹⁸F]flortaucipir R₁, lower values in the lateral temporal and parietal ROIs were associated with worse executive functioning, language and attention, and lower occipital R₁ with lower language and attention scores. When [¹⁸F]flortaucipir BP_ND and R₁ were modelled simultaneously, associations between lower R₁ in the lateral temporal ROI  and worse attention remained, as well as for lower parietal R₁ and worse executive functioning and attention.

Conclusion

Tau pathology was associated with locally reduced rCBF. Tau pathology and low rCBF were both independently associated with worse cognitive performance. For tau pathology, these associations spanned widespread neocortex, while for rCBF, independent associations were restricted to lateral temporal and parietal regions and the executive functioning and attention domains. These findings indicate that each biomarker may independently contribute to cognitive impairment in AD.

Electronic supplementary material

The online version of this article (10.1007/s00259-020-04831-w) contains supplementary material, which is available to authorized users.

Hippocampal [18F]flortaucipir BPND corrected for possible spill-in of the choroid plexus retains strong clinico-pathological relationships

BACKGROUND

Off-target [¹⁸F]flortaucipir (tau) PET binding in the choroid plexus causes spill-in into the nearby hippocampus, which may influence the correlation between [¹⁸F]flortaucipir binding and measures of cognition. Previously, we showed that partial volume correction (combination of Van Cittert iterative deconvolution and HYPR denoising; PVC HDH) and manually eroding the hippocampus resulted in a significant decrease of the choroid plexus spill-in. In this study, we compared three different approaches for the quantification of hippocampal [¹⁸F]flortaucipir signal using a semi-automated technique, and assessed correlations with cognitive performance across methods.

METHODS

Dynamic 130 min [¹⁸F]flortaucipir PET scans were performed in 109 subjects (45 cognitively normal subjects (CN) and 64 mild cognitive impairment/Alzheimer's disease (AD) dementia patients. We extracted hippocampal binding potential (BP_ND) using receptor parametric mapping with cerebellar grey matter as reference region. PVC HDH was performed. Based on our previous study in which we manually eroded 40% ± 10% of voxels of the hippocampus, three hippocampal volumes-of-interest (VOIs) were generated: a non-optimized 100% hippocampal VOI [100%], and combining HDH with eroding a percentage of the highest hippocampus BP_ND voxels (i.e. lowering spill-in) resulting in optimized 50%[50%HDH] and 40%[40%HDH] hippocampal VOIs. Cognitive performance was assessed with the Mini-Mental State Examination (MMSE) and Rey auditory verbal learning delayed recall. We performed receiver operating characteristic analyses to investigate which method could best discriminate MCI/AD from controls. Subsequently, we performed linear regressions to investigate associations between the hippocampal [¹⁸F]flortaucipir BP_ND VOIs and MMSE/delayed recall adjusted for age, sex and education.

RESULTS

We found higher hippocampal [¹⁸F]flortaucipir BP_ND in MCI/AD patients (BP_ND100%=0.27±0.15) compared to CN (BP_ND100%= 0.07±0.13) and all methods showed comparable discriminative effects (AUC_100%=0.85[CI=0.78-0.93]; AUC_50%HDH=0.84[CI=0.74-0.92]; AUC_40%HDH=0.83[CI=0.74-0.92]). Across groups, higher [¹⁸F]flortaucipir BP_ND was related to lower scores on MMSE (standardized β_100%=-0.38[CI=-0.57-0.20]; β_50%HDH= -0.37[CI=-0.54-0.19]; β_40%HDH=-0.35[CI=-0.53-0.17], all p<0.001) and delayed recall (standardized β_100%=-0.64[CI=-0.79-0.49]; β_50%HDH= -0.61[CI=-0.76-0.46]; β_40%HDH=-0.59[CI=-0.75-0.44]; all p<0.001), with comparable effect sizes for all hippocampal VOIs.

CONCLUSIONS

Hippocampal tau load measured with [¹⁸F]flortaucipir PET is strongly associated with cognitive function. Both discrimination between diagnostic groups and associations between hippocampal [¹⁸F]flortaucipir BP_ND and memory were comparable for all methods. The non-optimized 100% hippocampal VOI may be sufficient for clinical interpretation. However, proper correction for choroid plexus spillover and may be required in case of smaller effect sizes between subject groups or for longitudinal studies.

Associations between quantitative [18F]flortaucipir tau PET and atrophy across the Alzheimer's disease spectrum

Background

Neuropathological studies have linked tau aggregates to neuronal loss. To describe the spatial distribution of neurofibrillary tangle pathology in post-mortem tissue, Braak staging has been used. The aim of this study was to examine in vivo associations between tau pathology, quantified with [¹⁸F]flortaucipir PET in regions corresponding to Braak stages, and atrophy across the Alzheimer’s disease (AD) spectrum.

Methods

We included 100 subjects, including 58 amyloid-β positive patients with mild cognitive impairment (MCI, n = 6) or AD dementia (n = 52) and 42 controls with subjective cognitive decline (36% amyloid-β positive). All subjects underwent a dynamic [¹⁸F]flortaucipir PET to generate non-displaceable binding potential (BP_ND) maps. We extracted average [¹⁸F]flortaucipir BP_ND entorhinal, Braak III–IV (limbic) and Braak V–VI (neocortical) regions of interest (ROIs). T1-weighted MRI was used to assess gray matter (GM) volumes. We performed linear regression analyses using [¹⁸F]flortaucipir BP_ND ROIs as independent and GM density (ROI or voxelwise) as dependent variable.

Results

In MCI/AD subjects (age [mean ± SD] 65 ± 8 years, MMSE 23 ± 4), [¹⁸F]flortaucipir BP_ND was higher than in controls (age 65 ± 8, MMSE 29 ± 1) across all ROIs (entorhinal 0.06 ± 0.21 vs 0.46 ± 0.25 p < 0.001, Braak III–IV 0.11 ± 0.10 vs 0.46 ± 0.26, p < 0.001, Braak V–VI 0.07 ± 0.07 vs 0.38 ± 0.29, p < 0.001). In MCI/AD, greater [¹⁸F]flortaucipir BP_ND in entorhinal cortex was associated with lower GM density in medial temporal lobe (β − 0.40, p < 0.001). Greater [¹⁸F]flortaucipir BP_ND in ROI Braak III–IV and Braak V–VI was associated with smaller GM density in lateral and inferior temporal, parietal, occipital, and frontal lobes (range standardized βs − 0.30 to − 0.55, p < 0.01), but not in medial temporal lobe (β − 0.22, p 0.07). [¹⁸F]Flortaucipir BP_ND in ROI Braak I–II was not associated with GM density loss anywhere. When quantifying [¹⁸F]flortaucipir BP_ND across brain lobes, we observed both local and distant associations with GM atrophy. In controls, there were no significant associations between [¹⁸F]flortaucipir BP_ND and GM density (standardized βs ranging from − 0.24 to 0.02, all p > 0.05).

Conclusions

In MCI/AD patients, [¹⁸F]flortaucipir binding in entorhinal, limbic, and neocortical regions was associated with cortical atrophy.

Electronic supplementary material

The online version of this article (10.1186/s13195-019-0510-3) contains supplementary material, which is available to authorized users.

Partial-volume correction in dynamic PET-CT: effect on tumor kinetic parameter estimation and validation of simplified metrics

Background

Partial-volume effects generally result in an underestimation of tumor tracer uptake on PET-CT for small lesions, necessitating partial-volume correction (PVC) for accurate quantification. However, investigation of PVC in dynamic oncological PET studies to date is scarce. The aim of this study was to investigate PVC’s impact on tumor kinetic parameter estimation from dynamic PET-CT acquisitions and subsequent validation of simplified semi-quantitative metrics. Ten patients with EGFR-mutated non-small cell lung cancer underwent dynamic ¹⁸F-fluorothymidine PET-CT before, 7 days after, and 28 days after commencing treatment with a tyrosine kinase inhibitor. Parametric PVC was applied using iterative deconvolution without and with highly constrained backprojection (HYPR) denoising, respectively. Using an image-derived input function with venous parent plasma calibration, we estimated full kinetic parameters V_T, K₁, and k₃/k₄ (BP_ND) using a reversible two-tissue compartment model, and simplified metrics (SUV and tumor-to-blood ratio) at 50–60 min post-injection.

Results

PVC had a non-linear effect on measured activity concentrations per timeframe. PVC significantly changed each kinetic parameter, with a median increase in V_T of 11.8% (up to 25.1%) and 10.8% (up to 21.7%) without and with HYPR, respectively. Relative changes in kinetic parameter estimates vs. simplified metrics after applying PVC were poorly correlated (correlations 0.36–0.62; p < 0.01). PVC increased correlations between simplified metrics and V_T from 0.82 and 0.81 (p < 0.01) to 0.90 and 0.88 (p < 0.01) for SUV and TBR, respectively, albeit non-significantly. PVC also increased correlations between treatment-induced changes in simplified metrics vs. V_T at 7 (SUV) and 28 (SUV and TBR) days after treatment start non-significantly. Delineation on partial-volume corrected PET images resulted in a median decrease in metabolic tumor volume of 14.3% (IQR − 22.1 to − 7.5%), and increased the effect of PVC on kinetic parameter estimates.

Conclusion

PVC has a significant impact on tumor kinetic parameter estimation from dynamic PET-CT data, which differs from its effect on simplified metrics. However, it affected validation of these simplified metrics both as single measurements and as biomarkers of treatment response only to a small extent. Future dynamic PET studies should preferably incorporate PVC.

Trial registration

Dutch Trial Register, NTR3557.

Electronic supplementary material

The online version of this article (10.1186/s13550-019-0483-z) contains supplementary material, which is available to authorized users.

A novel partial volume correction method for accurate quantification of [18F] flortaucipir in the hippocampus

Background

Off-target binding in the choroid plexus (CP) may cause spill-in of the tau PET tracer [¹⁸F] flortaucipir into the adjacent hippocampus region. The impact of this spill-in on hippocampal uptake was assessed using a novel partial volume correction method (PVC).

Methods

PVC was performed on 20 [¹⁸F] flortaucipir dynamic PET scans (10 probable AD and 10 controls). Volumes of interest (VOIs) were defined for both hippocampus and CP. The correlation between hippocampal and CP distribution volume (V_T), with and without PVC, was determined. Both anatomically defined and eroded VOIs were used.

Results

For controls, the correlation between hippocampal and CP V_T was significantly reduced after using PVC along with an eroded VOI (r² = 0.59, slope = 0.80 versus r² = 0.15, slope = 0.15; difference: p < 0.05). The same was true for AD patients (p < 0.05).

Conclusion

PVC together with an optimized hippocampal VOI resulted in effective reduction of CP spill-in and improved accuracy of hippocampal V_T.

Electronic supplementary material

The online version of this article (10.1186/s13550-018-0432-2) contains supplementary material, which is available to authorized users.