Determination of phantom conditions
Here, we initially estimated the activity of [18F]flortaucipir, [18F]THK5351, and [18F]MK6240 tau PET tracers in the brain with reference to published scan parameters [34–36]. The conventional injected doses (MBq) and uptake durations (min) were 370 and 70, 185 and 40, and 185 and 90 for [18F]flortaucipir, [18F]THK5351, and [18F]MK6240, respectively [34–36], and the acquisition duration was 20 min for all three tracers. The estimated whole brain activity at the start of PET acquisition was 4.0, 1.0, 1.5 MBq for [18F]flortaucipir, [18F]THK5351, and [18F]MK6240, respectively. The standard brain volume was equivalent to 1,200 mL. The activity concentration in a normal brain was 3.33, 0.83, 1.25 kBq/mL for [18F]flortaucipir, [18F]THK5351, and [18F]MK6240, respectively.
We then calculated activity in the phantom based on estimated brain activity. The distribution of tau PET tracers in the brain corresponded to that of tau deposition defined in terms of Braak stages [3, 4]. Alzheimer disease can be detected early if PET imaging can detect local tau accumulation in the medial temporal lobe. Local accumulation is underestimated if a Hoffman 3D brain phantom (Data Spectrum Corporation, Hillsborough, NC, USA) contains whole brain activity that was described in previous paragraph. Therefore, we considered that the target (local accumulation) and reference regions could be mimicked by the Hoffman 3D brain phantom and a cylindrical phantom (Itoi Plastics Co. Ltd., Kobe, Japan), respectively, containing different amounts of activity. The average activity concentration of the three tau PET tracers in the whole brain was ~ 2.0 kBq/mL. The volume of the cylindrical phantom and the amount of activity in it were 6 L and 12.0 MBq, respectively, at the start PET image acquisition. The concentration of activity in the brain phantom was taken as the standardized uptake value ratio (SUVR) in the medial temporal lobe. The SUVR of the medial temporal lobe or hippocampus calculated from the activity concentrations of [18F]flortaucipir, [18F]THK5351, and [18F]MK6240 in the cerebellar cortex as a reference region was ~ 2.0 [37–44]. Thus, the concentration of activity in the brain phantom was twice that in the cylindrical phantom. The volume of the brain phantom and the activity in it were 1.2 L and 4.8 MBq, respectively, at the start of PET acquisition.
Characteristics of equipment at four sites
Radioactivity in the phantoms was quantified using two brands of dose calibrators (Nippon RayTech Co., Ltd., Tokyo, Japan and Capintec Inc., Ramsey, NJ, USA), and four PET/computed tomography (CT) scanners (one from GE Healthcare, Milwaukee, WI, USA and three from Siemens Healthineers, Erlangen, Germany). Table 1 shows the PET/CT systems and dose calibrators. Reconstruction conditions except for the iterations and the Gaussian filter (Table 1) proceeded under the clinical conditions for brain PET examinations at each site. The performance of the PET/CT scanners has been described elsewhere [45–48].
Table 1
Characteristics of equipment at four sites
Site | PET | Dose calibrator (Manufacture) |
Scanner (Manufacture) | Reconstruction conditions | |
NUH | Biograph 16 (SIEMENS) | 3D-OSEM; subset, 16; pixel size, 2 mm | IGC-7F (Aloka) |
QST | Biograph mCT Flow (SIEMENS) | 3D-OSEM + TOF; subset, 21; pixel size, 2 mm | IGC-3 (Aloka) |
NMS | Biograph Vision (SIEMENS) | 3D-OSEM + TOF; subset, 5; pixel size, 2 mm | CRC-55tR (Capintec) |
TMIG | Discovery MI (GE Healthcare) | 3D-OSEM + TOF; subset, 16; pixel size, 2 mm | CRC-55tR (Capintec) |
NMS, Nippon Medical School; NUH, Nagoya University Hospital; PET, positron emission tomography; QST, National Institutes for Quantum Science and Technology; TMIG, Tokyo Metropolitan Institute of Gerontology; TOF, time-of-flight; 3D-OSEM, three-dimensional ordered subset expectation maximization. |
Phantom experiment
Computed tomography images were acquired from all scanners to correct attenuation, scatter and other issues except for the point-spread function. Thereafter, images were acquired from the brain and cylindrical phantoms initially containing 4.8 and 12.0 MBq of [18F]FDG, respectively, for 20 min based on the protocols for [18F]flortaucipir, [18F]THK5351, and [18F]MK6240.
Creation of a VOI template for brain phantom images
We created a VOI template (Fig. 1) to optimize tau PET images based on Braak stages using PMOD v. 3.8 (PMOD Technologies LLC, Zurich, Switzerland). We then analyzed tau PET images acquired from the brain phantom. Each VOI was placed according to the distribution of tau pathology in patients with AD. The locations and amounts of voxels in the inferotemporal cortex, lateral temporal lobe, precuneus, white matter (WM), and cerebellar cortex were 357, 375, 365, 576, and 769 voxels, respectively. The inferotemporal cortex, lateral temporal lobe and precuneus comprised the grey matter (GM) VOI and the cerebellar cortex was the reference VOI.
Image analysis
Determination of iteration number
The data acquired from brain phantom images were reconstructed with 1–10 iterations and no post filter. Mean activity concentrations in the GM and WM were measured using the VOI template. The ratio of grey-to-white matter contrast (contrast [%]) and the recovery coefficient (RC) at the GM and WM were calculated as:
$$Contrast \left(\%\right)=\frac{\left({GM}_{p}/{WM}_{p}-1\right)}{\left({GM}_{d}/{WM}_{d}-1\right)}\times 100$$
,
where GMp and WMp in the brain, and GMd and WMd in the digital brain phantom PET images are GM and WM activities respectively, in VOIs. The GMd and WMd values provided a true gray-to-white ratio of 4 and were applied to the image co-registered to the digital phantom. Contrast was measured using PMOD v. 3.8.
The RC at GM and WM was defined as the image-derived mean activity concentration determined as contrast divided by the activity concentration of the stock solution in the brain phantom. The activity concentration in the brain phantom derived from net phantom activity (measured using a dose calibrator at each site), divided by the fillable volume (1.14 L) of the brain phantom. The activity concentration in WM was 25% of that in GM. The convergence of contrast or RC in 1–10 iterations was defined as optimal iteration for tau PET imaging.
Determination of Gaussian filter magnitude
Data acquired from the cylindrical phantom was reconstructed using the optimal number of iterations determined as contrast (%), RC, and Gaussian filter magnitudes of 0‒10 mm at full width at half maximum (FWHM). A large circular ROI (13 cm diameter; nROI) was placed on the center of the cylindrical phantom image to evaluate noise as a coefficient of variation (CV) calculated as:
$$CV \left(\%\right)=\frac{{SD}_{nROI}}{{nROI}_{mean}}\times 100$$
,
where SDnROI is the standard deviation of the voxel numbers within the nROI, and nROImean is mean nROI activity. The optimal magnitudes of Gaussian filters were determined from CVs < 15%.