Evaluation of novel data-driven metrics of amyloid ß deposition for longitudinal PET studies.

PURPOSE: Positron emission tomography (PET) provides in vivo quantification of amyloid-ß (Aß) pathology. Established methods for assessing Aß burden can be affected by physiological and technical factors. Novel, data-driven metrics have been developed to account for these sources of variability. We aimed to evaluate the performance of four of these amyloid PET metrics against conventional techniques, using a common set of criteria. METHODS: Three cohorts were used for evaluation: Insight 46 (N=464, [18F]florbetapir), AIBL (N=277, [18F]flutemetamol), and an independent test-retest data (N=10, [18F]flutemetamol). Established metrics of amyloid tracer uptake included the Centiloid (CL) and where dynamic data was available, the non-displaceable binding potential (BPND). The four data-driven metrics computed were the amyloid load (Aß load), the Aß-PET pathology accumulation index (Aß index), the Centiloid derived from non-negative matrix factorisation (CLNMF), and the amyloid pattern similarity score (AMPSS). These metrics were evaluated using reliability and repeatability in test-retest data, associations with BPND and CL, variability of the rate of change and sample size estimates to detect a 25% slowing in Aß accumulation. RESULTS: All metrics showed good reliability. Aß load, Aß index and CLNMF were strong associated with the BPND. The associations with CL suggest that cross-sectional measures of CLNMF, Aß index and Aß load are robust across studies. Sample size estimates for secondary prevention trial scenarios were the lowest for CLNMF and Aß load compared to the CL. CONCLUSION: Among the novel data-driven metrics evaluated, the Aß load, the Aß index and the CLNMF can provide comparable performance to more established quantification methods of Aß PET tracer uptake. The CLNMF and Aß load could offer a more precise alternative to CL, although further studies in larger cohorts should be conducted.

Uncertainty analysis of MR-PET image registration for precision neuro-PET imaging.

Accurate regional brain quantitative PET measurements, particularly when using partial volume correction, rely on robust image registration between PET and MR images. We argue here that the precision, and hence the uncertainty, of MR-PET image registration is mainly driven by the registration implementation and the quality of PET images due to their lower resolution and higher noise compared to the structural MR images. We propose a dedicated uncertainty analysis for quantifying the precision of MR-PET registration, centred around the bootstrap resampling of PET list-mode events to generate multiple PET image realisations with different noise (count) levels. The effects of PET image reconstruction parameters, such as the use of attenuation and scatter corrections and different number of iterations, on the precision and accuracy of MR-PET registration were investigated. In addition, the performance of four software packages with their default settings for rigid inter-modality image registration were considered: NiftyReg, Vinci, FSL and SPM. Four distinct PET image distributions made of two early time frames (similar to cortical FDG) and two late frames using two amyloid PET dynamic acquisitions of one amyloid positive and one amyloid negative participants were investigated. For the investigated four PET frames, the biggest impact on the uncertainty was observed between registration software packages (up to 10-fold difference in precision) followed by the reconstruction parameters. On average, the lowest uncertainty for different PET frames and brain regions was observed with SPM and two iterations of fully quantitative image reconstruction. The observed uncertainty for the varying PET count-level (from 5% to 60%) was slightly lower than for the reconstruction parameters. We also observed that the registration uncertainty in quantitative PET analysis depends on amyloid status of the considered PET frames, with increased uncertainty (up to three times) when using post-reconstruction partial volume correction. This analysis is applicable for PET data obtained from either PET/MR or PET/CT scanners.

Heterogeneity in Preclinical Alzheimer's Disease Trial Cohort Identified by Image-based Data-Driven Disease Progression Modelling.

Importance: Undetected biological heterogeneity adversely impacts trials in Alzheimer's disease because rate of cognitive decline - and perhaps response to treatment - differs in subgroups. Recent results show that data-driven approaches can unravel the heterogeneity of Alzheimer's disease progression. The resulting stratification is yet to be leveraged in clinical trials. Objective: Investigate whether image-based data-driven disease progression modelling could identify baseline biological heterogeneity in a clinical trial, and whether these subgroups have prognostic or predictive value. Design: Screening data from the Anti-Amyloid Treatment in Asymptomatic Alzheimer Disease (A4) Study collected between April 2014 and December 2017, and longitudinal data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) observational study downloaded in February 2022 were used. Setting: The A4 Study is an interventional trial involving 67 sites in the US, Canada, Australia, and Japan. ADNI is a multi-center observational study in North America. Participants: Cognitively unimpaired amyloid-positive participants with a 3-Tesla T1-weighted MRI scan. Amyloid positivity was determined using florbetapir PET imaging (in A4) and CSF Aß(1-42) (in ADNI). Main Outcomes and Measures: Regional volumes estimated from MRI scans were used as input to the Subtype and Stage Inference (SuStaIn) algorithm. Outcomes included cognitive test scores and SUVr values from florbetapir and flortaucipir PET. Results: We included 1,240 Aß+ participants (and 407 Aß- controls) from the A4 Study, and 731 A4-eligible ADNI participants. SuStaIn identified three neurodegeneration subtypes - Typical, Cortical, Subcortical - comprising 523 (42%) individuals. The remainder are designated subtype zero (insufficient atrophy). Baseline PACC scores (A4 primary outcome) were significantly worse in the Cortical subtype (median = -1.27, IQR=[-3.34,0.83]) relative to both subtype zero (median=-0.013, IQR=[-1.85,1.67], P<.0001) and the Subcortical subtype (median=0.03, IQR=[-1.78,1.61], P=.0006). In ADNI, over a four-year period (comparable to A4), greater cognitive decline in the mPACC was observed in both the Typical (-0.23/yr; 95% CI, [-0.41,-0.05]; P=.01) and Cortical (-0.24/yr; [-0.42,-0.06]; P=.009) subtypes, as well as the CDR-SB (Typical: +0.09/yr, [0.06,0.12], P<.0001; and Cortical: +0.07/yr, [0.04,0.10], P<.0001). Conclusions and Relevance: In a large secondary prevention trial, our image-based model detected neurodegenerative heterogeneity predictive of cognitive heterogeneity. We argue that such a model is a valuable tool to be considered in future trial design to control for previously undetected variance.

Tau accumulation in autosomal dominant Alzheimer's disease: a longitudinal [18F]flortaucipir study.

Cortical tau accumulation is a key pathological event that partly defines Alzheimer's disease (AD) onset and is associated with cognitive decline and future disease progression. However, an improved understanding of the timing and pattern of early tau deposition in AD and how this may be tracked in vivo is needed. Data from 59 participants involved in two longitudinal cohort studies of autosomal dominant AD (ADAD) were used to investigate whether tau PET can detect and track presymptomatic change; seven participants were symptomatic, and 52 were asymptomatic but at a 50% risk of carrying a pathogenic mutation. All had baseline flortaucipir (FTP) PET, MRI and clinical assessments; 26 individuals had more than one FTP PET scan. Standardised uptake value ratios (SUVRs) in prespecified regions of interest (ROIs) were obtained using inferior cerebellar grey matter as the reference region. We compared the changes in FTP SUVRs between presymptomatic carriers, symptomatic carriers and non-carriers, adjusting for age, sex and study site. We also investigated the relationship between regional FTP SUVRs and estimated years to/from symptom onset (EYO). Compared to both non-carriers and presymptomatic carriers, FTP SUVRs were significantly higher in symptomatic carriers in all ROIs tested (p < 0.001). There were no significant regional differences between presymptomatic carriers and non-carriers in FTP SUVRs, or their rates of change (p > 0.05), although increased FTP signal uptake was seen posteriorly in some individuals around the time of expected symptom onset. When we examined the relationship of FTP SUVR with respect to EYO, the earliest significant regional difference between mutation carriers and non-carriers was detected within the precuneus prior to estimated symptom onset in some cases. This study supports preliminary studies suggesting that presymptomatic tau tracer uptake is rare in ADAD. In cases where early uptake was seen, there was often a predilection for posterior regions (the precuneus and post-cingulate) as opposed to the medial temporal lobe, highlighting the importance of examining in vivo tau uptake beyond the confines of traditional Braak staging.

Operationalizing the centiloid scale for [18F]florbetapir PET studies on PET/MRI.

INTRODUCTION: The Centiloid scale aims to harmonize amyloid beta (Aß) positron emission tomography (PET) measures across different analysis methods. As Centiloids were created using PET/computerized tomography (CT) data and are influenced by scanner differences, we investigated the Centiloid transformation with data from Insight 46 acquired with PET/magnetic resonanceimaging (MRI). METHODS: We transformed standardized uptake value ratios (SUVRs) from 432 florbetapir PET/MRI scans processed using whole cerebellum (WC) and white matter (WM) references, with and without partial volume correction. Gaussian-mixture-modelling-derived cutpoints for Aß PET positivity were converted. RESULTS: The Centiloid cutpoint was 14.2 for WC SUVRs. The relationship between WM and WC uptake differed between the calibration and testing datasets, producing implausibly low WM-based Centiloids. Linear adjustment produced a WM-based cutpoint of 18.1. DISCUSSION: Transformation of PET/MRI florbetapir data to Centiloids is valid. However, further understanding of the effects of acquisition or biological factors on the transformation using a WM reference is needed. HIGHLIGHTS: Centiloid conversion of amyloid beta positron emission tomography (PET) data aims to standardize results.Centiloid values can be influenced by differences in acquisition.We converted florbetapir PET/magnetic resonance imaging data from a large birth cohort.Whole cerebellum referenced values could be reliably transformed to Centiloids.White matter referenced values may be less generalizable between datasets.

Advanced quantitative evaluation of PET systems using the ACR phantom and NiftyPET software.

PURPOSE: A novel phantom-imaging platform, a set of software tools, for automated and high-precision imaging of the American College of Radiology (ACR) positron emission tomography (PET) phantom for PET/magnetic resonance (PET/MR) and PET/computed tomography (PET/CT) systems is proposed. METHODS: The key feature of this platform is the vector graphics design that facilitates the automated measurement of the knife-edge response function and hence image resolution, using composite volume of interest templates in a 0.5 mm resolution grid applied to all inserts of the phantom. Furthermore, the proposed platform enables the generation of an accurate µ $\mu$ -map for PET/MR systems with a robust alignment based on two-stage image registration using specifically designed PET templates. The proposed platform is based on the open-source NiftyPET software package used to generate multiple list-mode data bootstrap realizations and image reconstructions to determine the precision of the two-stage registration and any image-derived statistics. For all the analyses, iterative image reconstruction was employed with and without modeled shift-invariant point spread function and with varying iterations of the ordered subsets expectation maximization (OSEM) algorithm. The impact of the activity outside the field of view (FOV) was assessed using two acquisitions of 30 min each, with and without the activity outside the FOV. RESULTS: The utility of the platform has been demonstrated by providing a standard and an advanced phantom analysis including the estimation of spatial resolution using all cylindrical inserts. In the imaging planes close to the edge of the axial FOV, we observed deterioration in the quantitative accuracy, reduced resolution (FWHM increased by 1-2 mm), reduced contrast, and background uniformity due to the activity outside the FOV. Although it slows convergence, the PSF reconstruction had a positive impact on resolution and contrast recovery, but the degree of improvement depended on the regions. The uncertainty analysis based on bootstrap resampling of raw PET data indicated high precision of the two-stage registration. CONCLUSIONS: We demonstrated that phantom imaging using the proposed methodology with the metric of spatial resolution and multiple bootstrap realizations may be helpful in more accurate evaluation of PET systems as well as in facilitating fine tuning for optimal imaging parameters in PET/MR and PET/CT clinical research studies.

Reduced acquisition time PET pharmacokinetic modelling using simultaneous ASL-MRI: proof of concept.

Pharmacokinetic modelling on dynamic positron emission tomography (PET) data is a quantitative technique. However, the long acquisition time is prohibitive for routine clinical use. Instead, the semi-quantitative standardised uptake value ratio (SUVR) from a shorter static acquisition is used, despite its sensitivity to blood flow confounding longitudinal analysis. A method has been proposed to reduce the dynamic acquisition time for quantification by incorporating cerebral blood flow (CBF) information from arterial spin labelling (ASL) magnetic resonance imaging (MRI) into the pharmacokinetic modelling. In this work, we optimise and validate this framework for a study of ageing and preclinical Alzheimer's disease. This methodology adapts the simplified reference tissue model (SRTM) for a reduced acquisition time (RT-SRTM) and is applied to [18F]-florbetapir PET data for amyloid-ß quantification. Evaluation shows that the optimised RT-SRTM can achieve amyloid burden estimation from a 30-min PET/MR acquisition which is comparable with the gold standard SRTM applied to 60 min of PET data. Conversely, SUVR showed a significantly higher error and bias, and a statistically significant correlation with tracer delivery due to the influence of blood flow. The optimised RT-SRTM produced amyloid burden estimates which were uncorrelated with tracer delivery indicating its suitability for longitudinal studies.

NiftyPET: a High-throughput Software Platform for High Quantitative Accuracy and Precision PET Imaging and Analysis.

We present a standalone, scalable and high-throughput software platform for PET image reconstruction and analysis. We focus on high fidelity modelling of the acquisition processes to provide high accuracy and precision quantitative imaging, especially for large axial field of view scanners. All the core routines are implemented using parallel computing available from within the Python package NiftyPET, enabling easy access, manipulation and visualisation of data at any processing stage. The pipeline of the platform starts from MR and raw PET input data and is divided into the following processing stages: (1) list-mode data processing; (2) accurate attenuation coefficient map generation; (3) detector normalisation; (4) exact forward and back projection between sinogram and image space; (5) estimation of reduced-variance random events; (6) high accuracy fully 3D estimation of scatter events; (7) voxel-based partial volume correction; (8) region- and voxel-level image analysis. We demonstrate the advantages of this platform using an amyloid brain scan where all the processing is executed from a single and uniform computational environment in Python. The high accuracy acquisition modelling is achieved through span-1 (no axial compression) ray tracing for true, random and scatter events. Furthermore, the platform offers uncertainty estimation of any image derived statistic to facilitate robust tracking of subtle physiological changes in longitudinal studies. The platform also supports the development of new reconstruction and analysis algorithms through restricting the axial field of view to any set of rings covering a region of interest and thus performing fully 3D reconstruction and corrections using real data significantly faster. All the software is available as open source with the accompanying wiki-page and test data.

Full field spatially-variant image-based resolution modelling reconstruction for the HRRT.

Accurate characterisation of the scanner's point spread function across the entire field of view (FOV) is crucial in order to account for spatially dependent factors that degrade the resolution of the reconstructed images. The HRRT users' community resolution modelling reconstruction software includes a shift-invariant resolution kernel, which leads to transaxially non-uniform resolution in the reconstructed images. Unlike previous work to date in this field, this work is the first to model the spatially variant resolution across the entire FOV of the HRRT, which is the highest resolution human brain PET scanner in the world. In this paper we developed a spatially variant image-based resolution modelling reconstruction dedicated to the HRRT, using an experimentally measured shift-variant resolution kernel. Previously, the system response was measured and characterised in detail across the entire FOV of the HRRT, using a printed point source array. The newly developed resolution modelling reconstruction was applied on measured phantom, as well as clinical data and was compared against the HRRT users' community resolution modelling reconstruction, which is currently in use. Results demonstrated improvements both in contrast and resolution recovery, particularly for regions close to the edges of the FOV, with almost uniform resolution recovery across the entire transverse FOV. In addition, because the newly measured resolution kernel is slightly broader with wider tails, compared to the deliberately conservative kernel employed in the HRRT users' community software, the reconstructed images appear to have not only improved contrast recovery (up to 20% for small regions), but also better noise characteristics.

Evaluation of a direct 4D reconstruction method using generalised linear least squares for estimating nonlinear micro-parametric maps.

OBJECTIVE: Estimation of nonlinear micro-parameters is a computationally demanding and fairly challenging process, since it involves the use of rather slow iterative nonlinear fitting algorithms and it often results in very noisy voxel-wise parametric maps. Direct reconstruction algorithms can provide parametric maps with reduced variance, but usually the overall reconstruction is impractically time consuming with common nonlinear fitting algorithms. METHODS: In this work we employed a recently proposed direct parametric image reconstruction algorithm to estimate the parametric maps of all micro-parameters of a two-tissue compartment model, used to describe the kinetics of [[Formula: see text]F]FDG. The algorithm decouples the tomographic and the kinetic modelling problems, allowing the use of previously developed post-reconstruction methods, such as the generalised linear least squares (GLLS) algorithm. RESULTS: Results on both clinical and simulated data showed that the proposed direct reconstruction method provides considerable quantitative and qualitative improvements for all micro-parameters compared to the conventional post-reconstruction fitting method. Additionally, region-wise comparison of all parametric maps against the well-established filtered back projection followed by post-reconstruction non-linear fitting, as well as the direct Patlak method, showed substantial quantitative agreement in all regions. CONCLUSIONS: The proposed direct parametric reconstruction algorithm is a promising approach towards the estimation of all individual microparameters of any compartment model. In addition, due to the linearised nature of the GLLS algorithm, the fitting step can be very efficiently implemented and, therefore, it does not considerably affect the overall reconstruction time.

Optimized data preprocessing for multivariate analysis applied to 99mTc-ECD SPECT data sets of Alzheimer's patients and asymptomatic controls.

Multivariate image analysis has shown potential for classification between Alzheimer's disease (AD) patients and healthy controls with a high-diagnostic performance. As image analysis of positron emission tomography (PET) and single photon emission computed tomography (SPECT) data critically depends on appropriate data preprocessing, the focus of this work is to investigate the impact of data preprocessing on the outcome of the analysis, and to identify an optimal data preprocessing method. In this work, technetium-99methylcysteinatedimer ((99m)Tc-ECD) SPECT data sets of 28 AD patients and 28 asymptomatic controls were used for the analysis. For a series of different data preprocessing methods, which includes methods for spatial normalization, smoothing, and intensity normalization, multivariate image analysis based on principal component analysis (PCA) and Fisher discriminant analysis (FDA) was applied. Bootstrap resampling was used to investigate the robustness of the analysis and the classification accuracy, depending on the data preprocessing method. Depending on the combination of preprocessing methods, significant differences regarding the classification accuracy were observed. For (99m)Tc-ECD SPECT data, the optimal data preprocessing method in terms of robustness and classification accuracy is based on affine registration, smoothing with a Gaussian of 12 mm full width half maximum, and intensity normalization based on the 25% brightest voxels within the whole-brain region.