
==== Front
EJNMMI Phys
EJNMMI Phys
EJNMMI Physics
2197-7364
Springer International Publishing Cham

39276263
682
10.1186/s40658-024-00682-8
Original Research
Quantitative SPECT imaging of 155Tb and 161Tb for preclinical theranostic radiopharmaceutical development
http://orcid.org/0000-0003-0794-9115
Koniar Helena hkoniar@triumf.ca

12
McNeil Scott 1
Wharton Luke 1
Ingham Aidan 1
Van de Voorde Michiel 3
Ooms Maarten 3
Sekar Sathiya 4
Rodríguez-Rodríguez Cristina 24
Kunz Peter 56
Radchenko Valery 17
Rahmim Arman 289
Uribe Carlos 8910
Yang Hua 16
Schaffer Paul 168
1 https://ror.org/03kgj4539 grid.232474.4 0000 0001 0705 9791 TRIUMF, Life Sciences Division, 4004 Wesbrook Mall, Vancouver, BC V6T 2A3 Canada
2 https://ror.org/03rmrcq20 grid.17091.3e 0000 0001 2288 9830 Department of Physics and Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, BC CAN V6T Canada
3 grid.8953.7 0000 0000 9332 3503 Institute for Nuclear Medical Applications Science, SCK CEN, Boeretang 200, Mol, BE 2400 Canada
4 https://ror.org/03rmrcq20 grid.17091.3e 0000 0001 2288 9830 Faculty of Pharmaceutical Sciences, University of British Columbia, 2405 Wesbrook Mall, Vancouver, BC V6T 1Z3 Canada
5 https://ror.org/03kgj4539 grid.232474.4 0000 0001 0705 9791 TRIUMF, Accelerator Division, 4004 Wesbrook Mall, Vancouver, BC V6T 2A3 Canada
6 https://ror.org/0213rcc28 grid.61971.38 0000 0004 1936 7494 Department of Chemistry, Simon Fraser University, 8888 University Drive, Burnaby, BC V5A 1S6 Canada
7 https://ror.org/03rmrcq20 grid.17091.3e 0000 0001 2288 9830 Department of Chemistry, University of British Columbia, 2036 Main Mall, Vancouver, BC V6T 1Z1 Canada
8 https://ror.org/03rmrcq20 grid.17091.3e 0000 0001 2288 9830 Department of Radiology, University of British Columbia, 2775 Laurel Street, Vancouver, BC V5Z 1M9 Canada
9 grid.248762.d 0000 0001 0702 3000 BC Cancer Research Centre, Department of Integrative Oncology, 675 W 10th Ave, Vancouver, BC V5Z 1L3 Canada
10 https://ror.org/03sfybe47 grid.248762.d 0000 0001 0702 3000 Functional Imaging, BC Cancer Agency, 600 West 10th Avenue, Vancouver, BC V5Z 4E6 Canada
14 9 2024
14 9 2024
12 2024
11 7722 4 2024
4 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Background

Element-equivalent matched theranostic pairs facilitate quantitative in vivo imaging to establish pharmacokinetics and dosimetry estimates in the development of preclinical radiopharmaceuticals. Terbium radionuclides have significant potential as matched theranostic pairs for multipurpose applications in nuclear medicine. In particular, 155Tb (t1/2 = 5.32 d) and 161Tb (t1/2 = 6.89 d) have been proposed as a theranostic pair for their respective applications in single photon emission computed tomography (SPECT) imaging and targeted beta therapy. Our study assessed the performance of preclinical quantitative SPECT imaging with 155Tb and 161Tb. A hot rod resolution phantom with rod diameters ranging between 0.85 and 1.70 mm was filled with either 155Tb (21.8 ± 1.7 MBq/mL) or 161Tb (23.6 ± 1.9 MBq/mL) and scanned with the VECTor preclinical SPECT/CT scanner. Image performance was evaluated with two collimators: a high energy ultra high resolution (HEUHR) collimator and an extra ultra high sensitivity (UHS) collimator. SPECT images were reconstructed from photopeaks at 43.0 keV, 86.6 keV, and 105.3 keV for 155Tb and 48.9 keV and 74.6 keV for 161Tb. Quantitative SPECT images of the resolution phantoms were analyzed to report inter-rod contrast, recovery coefficients, and contrast-to-noise metrics.

Results

Quantitative SPECT images of the resolution phantom established that the HEUHR collimator resolved all rods for 155Tb and 161Tb, and the UHS collimator resolved rods ≥ 1.10 mm for 161Tb and ≥ 1.30 mm for 155Tb. The HEUHR collimator maintained better quantitative accuracy than the UHS collimator with recovery coefficients up to 92%. Contrast-to-noise metrics were also superior with the HEUHR collimator.

Conclusions

Both 155Tb and 161Tb demonstrated potential for applications in preclinical quantitative SPECT imaging. The high-resolution collimator achieves < 0.85 mm resolution and maintains quantitative accuracy in small volumes which is advantageous for assessing sub organ activity distributions in small animals. This imaging method can provide critical quantitative information for assessing and optimizing preclinical Tb-radiopharmaceuticals.

Keywords

Terbium-155
Terbium-161
Radiopharmaceutical therapy
Theranostic pair
Preclinical imaging
Quantitative SPECT
National Science and Engineering Research Council of Canada (NSERC) Discovery GrantRGPIN-2022-03887 SAPIN2021–00030 RGPIN-2021–04093 RGPIN-2018–04997 Kunz Peter Radchenko Valery Yang Hua Schaffer Paul http://dx.doi.org/10.13039/501100000196 Canada Foundation for Innovation 25413 Government of Canada’s New Frontiers in Research Fund-ExplorationNFRFE-2019–00128 Kunz Peter National Research Council of CanadaNational Research Council of Canada issue-copyright-statement© Springer Nature Switzerland AG 2024
==== Body
pmcBackground

There are four terbium isotopes of clinical interest due to their versatile applications in nuclear medicine with suitable half-lives, diverse decay modes, and emission energies [1, 2]. 149Tb (t1/2 = 4.12 h), 152Tb (t1/2 = 17.5 h), 155Tb (t1/2 = 5.32 d), and 161Tb (t1/2 = 6.89 d) have applications in alpha therapy, PET, SPECT, and β-/Auger therapy, respectively [3–5]. Matched element-equivalent theranostic pairs have been investigated for many isotopes including 43Sc/44Sc/47Sc [6, 7], 64Cu/67Cu [8], 86Y/90Y [9], 123I/124I/131I [10], 132La/135La [11], 209At/211At [12], 203Pb/212Pb [13], and 225Ac/226Ac [14, 15]. Element-equivalent pairs benefit from identical radiolabelling conditions, pharmacokinetics, and biodistribution, facilitating quantitative assessment and optimization in radiopharmaceutical development and enabling accurate dosimetry estimates for clinical applications [16, 17].

Of the four terbium isotopes described above, 155Tb and 161Tb have been identified as a promising theranostic pair with 155Tb serving as a diagnostic isotope and 161Tb as a therapeutic isotope (Fig. 1). 155Tb decays purely via electron capture to stable 155Gd emitting photons with energies of 42.3 keV (30.7%), 43.0 keV (55.0%), 48.7 keV (17.4%), 49.3 keV (22.0%), 86.6 keV (32.0%), and 105.3 keV (25.1%) [18], making it a promising candidate for single photon emission computed tomography (SPECT) imaging. Quantitative SPECT imaging enables diagnostic information on disease spread and assists in pre-therapy dosimetry calculations [17, 19]. Furthermore, 155Tb emits an average of 13.9 Auger electrons and 0.77 conversion electrons per decay which could be leveraged for Auger electron therapy [20, 21]. 161Tb decays via β--decay to stable 161Dy with decay characteristics suitable for treating metastatic disease, emitting 154 keV mean energy β- particles and an average of 10.9 Auger electrons and 1.4 conversion electrons per decay [21, 22]. In addition to its therapeutic emissions, 161Tb also has photon emissions that enable SPECT imaging with energies at 25.7 keV (23.2%), 46.0 keV (11.2%), 48.9 keV (17.0%), and 74.6 keV (10.2%) [23]. The higher energy 48.9 keV and 74.6 keV gamma emissions are particularly useful for small animal preclinical SPECT imaging.

While 155Tb and 161Tb have been proposed as matched theranostic pair, given their individual decay characteristics they could also act as standalone theranostic agents with 155Tb for SPECT imaging and Auger electron therapy and 161Tb for SPECT imaging and β--therapy. Several routes for 155Tb and 161Tb production are currently under investigation for their feasibility in larger scale clinical quantities while maintaining high yield and radionuclidic purity [24, 25]. Although there are standing challenges in optimizing large scale routine production, this pair of Tb isotopes remain very promising for clinical applications in nuclear medicine [3].

Fig. 1 Decay characteristics of 155Tb (left) and 161Tb (right). X-ray and gamma photon emissionswith absolute intensity > 10% are shown. Half-lives and decay radiation properties were extracted from Livechart (IAEA Nuclear Data Section) based on the Nov 2023 ENSDF snapshot

Previous studies have demonstrated the in vivo potential and capabilities of Tb isotopes [1]. SPECT imaging of Derenzo hot rod resolution phantoms with both 155Tb and 161Tb has been previously conducted and revealed excellent spatial resolution [1, 22, 26, 27]. However, metrics on the limits of resolution, quantitative accuracy, and optimal image reconstruction methods for these Tb-based SPECT imaging systems have not yet been well characterized and require further investigation. Quantitative imaging at both the preclinical and clinical stages can assist in estimating dosimetry and optimizing therapy outcomes.

Preclinical studies have explored in vivo imaging with several promising radiopharmaceuticals including [155/161Tb]Tb-cm09 [1], [161Tb]Tb-PSMA-617 [28], [161Tb]Tb-DOTA-TOC [29], [155Tb]Tb-DOTA-TOC [27], [155/161Tb]Tb-crown-TATE [20], and [155/161Tb]Tb-crown-αMSH [30]. Further preclinical work involving 161Tb-radiopharmaceuticals has assessed the in vivo stability [31] and biodistribution profiles [32] but often lacks corresponding quantitative SPECT imaging data.

In-human clinical invesitagtions with 161Tb have remained limited relative to the β- emitter 177Lu, in part, due to production challenges however, there have been promising initial outcomes. Clinical protocols for 161Tb SPECT/CT imaging have been established with high resolution capabilities [33]. In two patients with metastatic neuroendocrine neoplasms, [161Tb]Tb-DOTA-TOC SPECT/CT imaging demonstrated the feasibility of visualizing even small metastases with relatively low injected activities [34]. Following a single treatment with [161Tb]Tb-PSMA-617, one patient with advanced metastatic castration-resistant prostate cancer exhibited impressive partial remission, even with disease progression following extensive [177Lu]Lu-PSMA-617 therapy [35]. Additionally, post-therapy SPECT/CT imaging of [161Tb]Tb-PSMA-617 is feasible for visualizing its distribution in the targeted lesions and non-target organs [36]. The ongoing VIOLET Phase I/II clinical trial aims to evaluate the therapeutic efficacy and safety of [161Tb]Tb-PSMA-I&T for metastatic prostate cancer [37]. Post-therapy quantitative SPECT/CT imaging was also evaluated over a range of activities and imaging times to produce suitable imaging for clinical review [38]. Further, in silico dosimetry studies have indicated that 161Tb may be a better candidate than 177Lu for irradiating single cancer cells and micrometastasis [39].

Evidently there has been tremendous interest in the development of theranostic pairs with 155Tb and 161Tb. As production and purification methods improve Tb availability for preclinical and clinical studies, quantitative SPECT imaging protocols are necessary to further optimize the development pipeline of novel Tb-radiopharmaceuticals. In this work, we assess and evaluate the imaging performance of 155Tb and 161Tb with a preclinical SPECT/CT scanner and establish a SPECT/CT protocol for quantitative imaging. Element-equivalent matched theranostic pairs provide accurate biodistribution and dosimetry estimates, and this quantitative pharmacokinetic information is essential in optimizing the efficacy of radiopharmaceuticals in preclinical settings.

Methods

155Tb and 161Tb production and activity quantification

155Tb was produced at TRIUMF (Vancouver, Canada) with its Isotope Separation and Acceleration (ISAC) facility. 480 MeV protons were used to irradiate a tantalum target, resulting in a heterogeneous ion beam of spallation products that are mass separated by the principles of isotope separation on-line (ISOL) [40, 41]. The A/q = 155 beam was implanted onto an aluminum implantation target with a thin ammonium chloride salt layer. 155Er (t1/2 = 5.3 m), 155Ho (t1/2 = 48 m), and 155Dy (t1/2 = 9.9 h) were co-implanted with the 155Tb activity, however, all these isotopes decay into 155Tb via β+/EC-decay during a 5-day cool-down period after end of beam delivery. The salt layer on the implantation target was dissolved with small volumes of water (< 500 µL) prior to transfer to a sample vial for further dilution into imaging phantoms. Production methods are described only briefly here and are fully described in Fiaccabrino et al. (2021) [41]. Yield measurements from all production runs are documented in the ISAC Yield Database [42].

161Tb was produced at SCK CEN (Mol, Belgium) with its BR2 reactor [22, 43]. A highly enriched 160Gd target (97.5% enrichment, as 160Gd2O3) was irradiated for 7 days using a high thermal neutron flux of 3 × 1014 neutrons/cm2/s. Neutron activation of the enriched gadolinium (160Gd) target yields the short-lived 161Gd (t1/2 = 3.66 m) which decays into 161Tb via β--decay. The target was dissolved with high purity 1 M HNO3. 161Tb was purified from the 160Gd target material and 161Dy decay product with an automated system of solid-phase extraction columns. Purification methods are described briefly here and are fully described in McNeil et al. (2022) [32].

Small samples of the 155Tb and 161Tb solutions were collected for activity quantifications. The activity was measured by gamma spectroscopy with a high-purity geranium (HPGe) detector (Canberra Industries, Meriden, CT) and analyzed with the Genie 2000 software package (Canberra Industries, Meriden, CT). Activity measurements have 8% estimated uncertainties, compounded from volume measurement and calibration source uncertainties [44]. For this study, activity concentrations of 21.8 ± 1.7 MBq/mL of 155Tb and 23.6 ± 1.9 MBq/mL of 161Tb were used for point source calibration and resolution phantom imaging.

SPECT image acquisition and reconstruction

The performance of 155/161Tb-based imaging was assessed with the VECTor (Versatile Emission Computed Tomography) microSPECT/CT (MILabs, Utretch, Netherlands). Its imaging capacities have been well described by Goorden et al. (2013) [45]. Briefly, the VECTor scanner consist of three stationary NaI gamma camera detectors placed in a triangle with pinhole collimators designed to collect data in a central field of view (CFOV) [46]. With a translation stage bed, phantoms and animals can be moved throughout imaging acquisition to collect scintillation events in the pinholes from different regions in the scanned subject [47].

While MILabs produces many collimators compatible with the VECTor scanner, two collimators were evaluated for their scanning performance: a high energy ultra high resolution (HEUHR) collimator and an extra ultra high sensitivity (UHS) collimator. The HEUHR collimator is ideal for high energy photons (> 350 keV) and high spatial resolution [48] while the UHS collimator is optimized for low energy photons (< 350 keV) and high sensitivity [49]. The HEUHR collimator is made of tungsten with an inner diameter of 44 mm and consists of 48 four-pinhole clusters with 0.35 mm diameter pinholes. When assessed with 99mTc and 18F, its peak sensitivity was 2,800 cps/MBq and 2,899 cps/MBq and the reconstructed spatial resolution was 0.5 mm and 0.8 mm, respectively. The UHS collimator is made of lead with an inner diameter of 46 mm and consists of 54 conical 2.0 mm diameter pinholes. When assessed with 99mTc, the peak sensitivity was 13,080 cps/MBq and reconstructed spatial resolution was 0.85 mm. Both collimators were evaluated to characterize 155Tb and 161Tb image quality and its applications for in vivo quantitative SPECT imaging. Robertson et al. [50] and Koniar et al. [14] serve as a basis for this work in our assessment of Tb-based imaging [14, 50].

Data was acquired in list mode and sorted into 512 energy bins with 2.34 keV width. The FOV was set to 54 mm × 54 mm × 140 mm for both the UHS and HEUHR collimator for all image acquisitions. The images were reconstructed from photopeaks at 43.0 keV, 86.6 keV, and 105.3 keV for 155Tb and at 48.9 keV and 74.6 keV for 161Tb (see Fig. 2). The 43.0 keV photopeak for 155Tb and the 48.9 keV photopeak for 161Tb are the result of multiple photon emissions within a close energy range. More precisely, the 43.0 keV photopeak is a combination of 42.3 keV (30.7%), 43.0 keV (55.0%), 48.7 keV (17.4%), 49.3 keV (22.0%) emissions and the 48.9 keV photopeak is a combination of 46.0 keV (11.2%) and 48.9 keV (17.0%) emissions. For simplicity and clarity, reconstructions from these photopeak are referred to by their most probable photon emission, namely the 43.0 keV and 48.9 keV emissions. Figure 2. demonstrates the relative count rate, or sensitivity normalized to the activity concentration and scan time, for both the HEUHR and UHS collimators as a function of photon energy and corresponds well with published photon intensities seen in Fig. 1.

Fig. 2 Photon energy spectrum of 155Tb (left) and 161Tb (right)

A pixel-based ordered subset expectation maximization (OSEM) iterative reconstruction algorithm (16 subsets, 6 iterations) was used to reconstruct the images with 0.4 mm3 voxel size. This number of subsets and iterations (96 MLEM-equivalent iterations) was determined to optimize the contrast-to-noise ratio for most resolution phantom images when considering different combinations of isotopes (155Tb and 161Tb) and collimators (UHS and HEUHR). The contrast-to-noise ratio calculation is presented in further detail below. The reconstruction algorithm was maintained for all images to investigate the effect of photopeak and collimator choice on resulting SPECT image quality. SPECT images were corrected for background and scatter by the triple energy window method and for attenuation by CT image co-registration (Table 1). The width of energy windows was selected to include mostly primary photons counted by the detectors. The width of the upper and lower energy background windows was set with 2.5 background weight as specified by the manufacturer. SPECT images were re-sampled to CT images with 0.169 mm3 voxels. Non-uniform attenuation correction was employed post-reconstruction with a CT-based attenuation map [51].

Table 1 Energy window parameters for SPECT image reconstruction

Isotope	Photopeak energy	Photopeak window	Upper/Lower window	
155Tb	43.0 keV	50%	4.3%	
86.6 keV	18%	3.1%	
105.3 keV	15%	3.2%	
161Tb	48.9 keV	50%	4.9%	
74.6 keV	20%	3.6%	

Point source calibration

Quantitative SPECT images were obtained via calibration factors determined through SPECT scans of activity point source phantoms for each collimator and isotope combination as described by the manufacturer. Activity concentration was also determined independently via gamma spectroscopy. Point source phantoms consisted of Eppendorf tubes filled with decay corrected 155Tb or 161Tb activity as listed in Table 2. Scan time, and total counts in the photopeak and background windows for all combinations of reconstruction photopeak and collimator are presented in Table 2. Calibration factors were determined with volumes of interest centered inside the point source phantom SPECT reconstructions for all collimator and photopeak reconstruction combinations. All reconstructed SPECT/CT images were analyzed with Python (v3.8.8) scripts.

Table 2 Point source phantom scan time, activity decay-corrected to time of imaging, and counts in each reconstruction window for each isotope and collimator combination

		UHS	HEUHR	
Isotope	Photopeak energy	Scan time	Activity concentration (MBq/mL)	Photopeak counts
(millions)	Background counts
(millions)	Scan time	Activity concentration (MBq/mL)	Photopeak counts
(millions)	Background counts
(millions)	
155Tb	43.0 keV	15 min	47.1 ± 3.8	35.83	2.61	30 min	46.6 ± 3.7	14.00	0.96	
	86.6 keV	13.80	2.06	5.00	0.75	
	105.3 keV	10.61	2.15	4.17	0.74	
161Tb	48.9 keV	30 min	23.2 ± 1.9	32.56	1.35	30 min	23.5 ± 1.9	6.01	0.29	
	74.6 keV	8.99	0.82	1.57	0.19	

Resolution phantom

A hot rod resolution phantom with thin rod clusters of diameters 0.85 mm, 0.95 mm, 1.10 mm, 1.30 mm, 1.50 mm, and 1.70 mm was used to assess the resolution of the VECTor imaging system (see Fig. 3). The resolution phantom was filled with decay corrected activities of 155Tb (21.4 ± 1.7 MBq/mL for the UHS scan; 20.9 ± 1.7 MBq/mL for the HEUHR scan) or with 161Tb (23.4 ± 1.9 MBq/mL for the UHS scan; 23.5 ± 1.9 MBq/mL for the HEUHR scan). Due to the high activity concentrations of 155Tb and 161Tb, images were acquired rapidly, with 60 min total acquisition time for all isotope and collimator combinations. Quantitative SPECT images of the resolution phantoms were analyzed to report spatial resolution, quantitative accuracy, and noise metrics as defined by Walker et al. (2014) [52]. Briefly here, circular ROIs with 0.8 times the physical diameter were placed inside the rods and in between the rods to measure the mean value and repeated across 30 planes for an axial thickness of 5.1 mm.

Fig. 3 Sample of a cross section of the resolution phantom with 0.85–1.70 mm diameter rods. Red and green circles represent the hot (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{h}_{d}$$\end{document}) and background (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{b}_{d}$$\end{document}) ROIs, resepectively, for each rod diameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:d$$\end{document}

Inter-rod contrast (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{C}_{d}$$\end{document}) quantifies resolvability in small hot regions. Inter-rod contrast for rods with diameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:d$$\end{document} was defined as:\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{C}_{d}=\frac{\stackrel{-}{{h}_{d}}-\stackrel{-}{{b}_{d}}}{\stackrel{-}{{h}_{d}}}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\stackrel{-}{{h}_{d}}$$\end{document} is the mean value inside the rods and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\stackrel{-}{{b}_{d}}$$\end{document} is the mean value in between the rods. The ideal inter-rod contrast is 1, and rods with inter-rod contrast greater than 0.2 were visually resolvable.

The recovery coefficient (RC) quantifies the accuracy of the apparent activity concentrations. The RC was defined as:\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:RC=\frac{\stackrel{-}{{h}_{d}}}{{h}_{0}}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{h}_{0}$$\end{document} is the activity concentration determined independently via gamma spectroscopy. The ideal RC is also 1.

The measure for quantifying the variability between ROI mean values for rods with diameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:d$$\end{document} was defined as:\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{N}_{d}=\frac{\sqrt{{\sigma\:}_{{h}_{d}}^{2}+{\sigma\:}_{{b}_{d}}^{2}}}{\overline{RO{I}_{d}}}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{\sigma\:}_{{h}_{d}}^{2}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{\sigma\:}_{{h}_{d}}^{2}$$\end{document}is the variance in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{h}_{d}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{b}_{d}$$\end{document}, calculated over 6 of the 30 image slices to reduce covariance between adjacent planes, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\stackrel{-}{RO{I}_{d}}$$\end{document} is the mean value from all ROIs (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{h}_{d}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{b}_{d}$$\end{document}) with rod diameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:d$$\end{document}. The contrast-to-noise ratio (CNR) was defined as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{C}_{d}/{N}_{d}$$\end{document} as another measure for image quality.

Impact of activity concentration on image quality metrics

Spatial resolution and noise metrics were also assessed in SPECT images reconstructed with fewer counts, to simulate lower activity concentrations or shorter scan times more commonly seen in preclinical in vivo imaging. Count reductions randomly select a portion of the list mode acquired data to be ignored during image reconstruction and emulate images acquired with lower activity concentrations. Images were reconstructed with 20%, 5%, 1%, and 0.5% of the total counts to investigate the effects on image quality.

Results

Resolution phantom

Figure 4 shows quantitative SPECT images of the resolution phantom acquired with the HEUHR and UHS collimators with various photon emission photopeak reconstructions. Qualitatively, it is apparent that the HEUHR collimator resolves rods better than the UHS collimator for all 155Tb and 161Tb emission photopeak reconstructions. Further, it is visually evident that the HEUHR collimator resolves all rods, indicating a resolution of < 0.85 mm for all reconstructions.

Fig. 4 Quantitative SPECT images of the resolution phantom reconstructed with various photopeak windows after acquisition with the HEUHR and UHS collimators. No post-reconstruction smoothing was applied

Inter-rod contrast and resolution limit

Inter-rod contrast measurements (Fig. 5) determine the minimum resolvability of SPECT images of the resolution phantom. Images acquired with the HEUHR collimator had inter-rod contrast values > 40% for all rods and photopeak reconstructions, indicating that the imaging system’s minimum resolvability is < 0.85 mm. Images reconstructed from higher energy photopeaks (105 keV for 155Tb and 75 keV for 161Tb) had higher inter-rod contrast values than those reconstructed from lower energy photopeaks (45 keV and 87 keV for 155Tb and 49 keV for 161Tb). This is seen qualitatively in Fig. 4 by less blurring, or lower apparent activity, in between the rods for SPECT images reconstructed from higher energy photopeaks. However, in images acquired with the UHS collimator, only rods with diameters ≥ 1.10 mm for 161Tb and ≥ 1.30 mm for 155Tb were resolvable. Unlike with the HEUHR collimator, the inter-rod contrast metrics from images acquired with the UHS collimator did not vary significantly with different photopeak energy reconstructions.

Fig. 5 Inter-rod contrast (Cd) measurements from SPECT images of the resolution phantom acquired with the HEUHR (left) and UHS (right) collimators with 155Tb in red and 161Tb in blue

At lower activity concentrations, the inter-rod contrast with the HEUHR collimator remained consistent down to 1 MBq/mL (5% total counts) maintaining the same overall resolvability seen in the full count SPECT images (see Fig. 6). Inter-rod contrast with the UHS collimator was less affected by lower total counts than the HEUHR collimator as system resolution remained consistent for all activity concentrations ≥ 0.1 MBq/mL (0.5% total counts).

Fig. 6 Inter-rod contrast measurements from SPECT images of the resolution phantom with reduced activity concentrations, simulated by reconstructing with 20%, 5%, 1%, and 0.5% of the total acquired counts. SPECT images were acquired with the HEUHR (left) and UHS (right) collimators with 155Tb in red and 161Tb in blue

Recover coefficient and quantitative accuracy

Figure 7 shows the RC metrics for evaluating the quantitative accuracy of SPECT images of the resolution phantom. For all photopeak reconstructions of 155Tb and 161Tb, the HEUHR collimator outperformed the UHS collimator with higher RC values indicative of better quantitative accuracy. For the 1.70 mm diameter rods, RC values ranged between 71 and 91% with the HEUHR collimator and between 69 and 75% with the UHS collimator. Like the inter-rod contrast value results presented above, SPECT images acquired with the HEUHR collimator and reconstructed from higher energy photopeaks had better quantitative accuracy. The best overall quantitative accuracy was seen in the SPECT images of the 155Tb-filled resolution phantom reconstructed from the high energy 105 keV photopeak; with RC values > 75% for rods with diameters ≥ 1.1 mm and maximum value of 92% quantitative accuracy. Similarily, SPECT images of the 161Tb-filled resolution phantom showed better quantitative accuracy with its 75 keV reconstruction than the lower energy 49 keV reconstruction, with higher RC values demonstrated for all rod diameters. The RC values for images acquired with the UHS collimator were relatively independent from the photopeak energy.

Fig. 7 Recovery coefficients from SPECT images of the resolution phantom acquired with the HEUHR (left) and UHS (right) collimators with 155Tb in red and 161Tb in blue

At lower activity concentrations, RC values become more variable and are generally lower than those seen in the full count SPECT images (see Fig. 8). The HEUHR collimator maintains consistent quantitative accuracy for activity concentrations ≥ 1MBq/mL while the UHS collimator maintains performance with activity concentrations ≥ 0.1 MBq/mL.

Fig. 8 Recovery coefficients from SPECT images of the resolution phantom with reduced activity concentrations, simulated by reconstructing with 20%, 5%, 1%, and 0.5% of the total acquired counts. SPECT images were acquired with the HEUHR (left) and UHS (right) collimators with 155Tb in red and 161Tb in blue

Contrast-to-noise metrics

CNR values are presented in Fig. 9. Images acquired with the HEUHR collimator generally had a higher CNR with increasing rod diameter, although the CNR plateaued for 155Tb images. The lower energy photopeak reconstructions (43 keV for 155Tb and 49 keV for 161Tb) slightly outperformed the higher energy photopeak reconstructions. Images acquired with the UHS collimator demonstrated notably better CNRs with the 155Tb-filled phantom than the 161Tb-filled phantom. The reconstructed photopeak energy did not have a meaningful effect on CNR values for images acquired with the UHS collimator. For both collimators, CNRs are positively correlated to the activity concentration as seen by decreasing CNRs with lower total counts (see Fig. 10).

Fig. 9 Contrast-to-noise ratios from SPECT images of the resolution phantom acquired with the HEUHR (left) and UHS (right) collimators with 155Tb in red and 161Tb in blue

Fig. 10 Contrast-to-noise ratios from SPECT images of the resolution phantom with reduced activity concentrations, simulated by reconstructing with 20%, 5%, 1%, and 0.5% of the total acquired counts. SPECT images were acquired with the HEUHR (left) and UHS (right) collimators with 155Tb in red and 161Tb in blue

Discussion

In this work, we have evaluated the performance of two different collimators for preclinical SPECT imaging of 155Tb and 161Tb by quantifying their image characteristics and limitations with a hot rod resolution phantom.

The HEUHR collimator universally achieves a resolution of < 0.85 mm, making it suitable for in vivo imaging in mice. At this scale, sub-organ localization of Tb-labelled radiopharmaceuticals can be discerned. This level of resolution is consistent with findings from previous studies that evaluated this collimator [48]. With the HEUHR collimator, SPECT imaging with of 225Ac’s progeny isotopes 221Fr (218 keV) and 213Bi (440 keV) achieved comparable resolution limits below 0.85 mm [50]. Similarly, SPECT imaging of 226Ac with its 158 keV and 230 keV gamma emissions also had resolvability < 0.85 mm with the same collimator [14].

Images acquired with the HEUHR collimator exhibited better resolved rods and noticeably higher inter-rod contrast when reconstructed from higher energy photopeaks. This is likely due to lower energy photons are having a higher probability of scattering and thus requiring wider photopeak energy windows (as shown in Table 1). The addition of scattered photons in the primary photopeak window results in blurring and degraded rod profiles leading to less contrast between rods [53]. Considering the low energy (< 140 keV) of photons emitted from 155Tb and 161Tb, this result is consistent with the findings of previous work with this imaging system. The improved inter-rod contrast could be attributed to the image reconstruction system matrix provided by the manufacturer which is optimized for 140 keV energy photons. We have observed that as the photon energy approaches 140 keV, the inter-rod contrast improves. The effect of photon energy on inter-rod contrast has also been seen with higher energy (> 140 keV) photons from the 225Ac decay chain, where the 218 keV (221Fr) image was better than the 440 keV (213Bi) image [50] and for 226Ac where the 158 keV image was better than the 230 keV image [14]. Despite being designed for collimating high energy photons, this work indicates that the HEUHR collimator is well suited to maintain its resolution for photon emissions with energies ≥ 40 keV with sufficient activity concentrations (see Fig. 6). The resolution and quantitative accuracy demonstrated in this preclinical imaging system benefits from the relatively small proportion of attenuated photons in small animals. These imaging cabapilities might not translate to clinical scanners since there will be significantly more scatter and attenuation of low energy photon emissions in humans. However, the ability to quantify preclinical images should not be understated in the optimization and development of preclinical Tb-radiopharmaceuticals.

The UHS collimator yielded resolvability of rods ≥ 1.10 mm for 161Tb and ≥ 1.30 mm for 155Tb, consistent with previous work on this imaging system. With the UHS collimator, SPECT images of 225Ac’s progeny 221Fr (218 keV) resolved rods ≥ 1.30 mm [50] and images of 226Ac resolved rods ≥ 1.30 mm and ≥ 1.50 mm for the 158 keV and 230 keV images, respectively [14]. Inter-rod contrast with the UHS collimator was less affected by photopeak energy than the HEUHR collimator which could be attributed to the higher sensitivity of the UHS collimator which permits a higher relative count rate than the HEUHR collimator (Fig. 2) and lessens the effect of including scattered photons in the primary photopeak window. While the UHS collimator does not have as powerful resolution capabilities as the HEUHR collimator, it is less affected by lower activity concentrations. This is especially important in considerations for in vivo imaging applications with generally lower activity concentrations and shorter scan times than those seen in a phantom imaging study.

Quantitative accuracy in small volumes is important to assess in a resolution phantom to understand the impact of partial volume effects on in vivo activity measurements. The RC values from images acquired with the HEUHR collimator showed the best quantitative accuracy with the high energy 105 keV gamma from 155Tb demonstrating the highest RC values overall for all rod diameters. The UHS collimator does not perform as well in terms of maintaining quantitative accuracy in small volumes. This is expected given the results of the inter-rod contrast measurements which limit the inherent system resolution to ≥ 1.10 mm for 161Tb and ≥ 1.30 mm for 155Tb. In future in vivo applications, it is important to account for these RC values to estimate true activity concentrations more accurately from the apparent activity measured in small volumes from quantitative SPECT images.

In high-count SPECT images, the HEUHR collimator provides high CNR between small hot volumes with both 155Tb and 161Tb. The UHS collimator provides notably higher CNRs for 155Tb than 161Tb. As such, 155Tb is the stronger candidate for imaging in the development of preclinical Tb-radiopharmaceuticals. When considering in vivo imaging scenarios, it is important to note that low-count SPECT images will exhibit lower CNRs in small hot volumes such as the adrenal glands or gallbladder in mice.

Overall, the HEUHR collimator provides better spatial resolution and quantitative accuracy in SPECT images of the resolution phantom than the UHS collimator. However, the UHS collimator maintains better image quality in images reconstructed from simulated low activity concentrations ≥ 0.1 MBq/mL. For preclinical in vivo imaging applications, a collimator should be selected based on anticipated activity concentrations from injected activity levels and anticipated uptake in tumours or organs of interest. In cases where high-count SPECT images can be acquired the HEUHR collimator is superior and in cases where only low-count SPECT images can be obtained the UHS collimator is better suited. Images reconstructed from higher energy photopeaks demonstrated better overall image quality than those reconstructed from their counterpart lower energy photopeaks for some activity concentration and collimator combinations. The flexibility of image reconstruction from multiple different gamma emissions photopeaks is a strong advantage of SPECT imaging with 155Tb and 161Tb. It is quite feasible to reconstruct several images from different photon energies after an in vivo scan acquisition to optimize both quantitative and qualitative image quality metrics. This ability to adapt image reconstruction parameters is a strength to this set of isotopes, in addition to its element-equivalent in vivo behaviour as theranostic isotopes. These advantages can facilitate the preclinical development of radiopharmaceuticals in a more efficient manner.

Conclusion

Quantitative SPECT imaging of 155Tb and 161Tb can be performed with high spatial resolution using the preclinical small-animal SPECT/CT VECTor scanner. A high-resolution collimator is superior for quantitative accuracy in high count SPECT images; however, a high sensitivity collimator provides better qualitative images in low count SPECT images. Both 155Tb and 161Tb have multiple photon emissions that provide adaptability in image reconstruction protocols to optimize image quality for the development of novel preclinical Tb-based radiopharmaceuticals.

Abbreviations

CNR Contrast to noise ratio

CT Computed tomography

EC Electron capture

ENSDF Evaluated nuclear structure data file

HEUHR High energy ultra high resolution

HPGe High-purity geranium

ISAC Isotope separation and acceleration

ISOL Isotope separation on-line

MLEM Maximum likelihood expectation maximization

OSEM Ordered subset expectation maximization

RC Recovery coefficient

ROI Region of interest

SPECT Single photon emission computed tomography

UHS Extra ultra high sensitivity

VECTor Versatile emission computed tomography

Acknowledgements

Not applicable.

Author contributions

HK, SM, AR, CU, HY, and PS devised the study conceptualization. LW, AI, PK, VR, and HY produced, purified, and quantified 155Tb activity at TRIUMF (Vancouver, Canada). MV and MO produced 161Tb activity at SCK CEN (Mol, Belgium). SM purified 161Tb activity. HK, SS, and CR acquired SPECT/CT scans. HK reconstructed, quantified, and analyzed SPECT images of the resolution phantom. HK drafted the manuscript. All authors read and approved the final manuscript.

Data availability

The datasets generated and analysed in the current study are available from the corresponding author on reasonable request.

Declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

This study was financially supported by the following NSERC (National Science and Engineering Research Council of Canada) Discovery Grants: RGPIN-2022-03887 (HY), SAPIN2021–00030 (PK), RGPIN-2021–04093 (PS), RGPIN-2018–04997 (VR). Additionally, funding was obtained through the Canada Foundation for Innovation (CFI) project no. 25413 (CR) and the Government of Canada’s New Frontiers in Research Fund-Exploration NFRFE-2019–00128 (PK). TRIUMF receives federal funding via a contribution agreement with the NRC (National Research Council of Canada).

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Müller C Zhernosekov K Köster U Johnston K Dorrer H Hohn A A Unique Matched Quadruplet of Terbium Radioisotopes for PET and SPECT and for α- and β–-Radionuclide Therapy: an in vivo proof-of-Concept Study with a new receptor-targeted Folate Derivative J Nucl Med 2012 53 1951 10.2967/jnumed.112.107540 23139086
Müller C, Zhernosekov K, Köster U, Johnston K, Dorrer H, Hohn A, et al. A Unique Matched Quadruplet of Terbium Radioisotopes for PET and SPECT and for α- and β–-Radionuclide Therapy: an in vivo proof-of-Concept Study with a new receptor-targeted Folate Derivative. J Nucl Med. 2012;53:1951.23139086 10.2967/jnumed.112.107540
2. Van Laere C Koole M Deroose CM de Voorde MV Baete K Cocolios TE Terbium radionuclides for theranostic applications in nuclear medicine: from atom to bedside Theranostics 2024 14 1720 43 10.7150/thno.92775 38389843
Van Laere C, Koole M, Deroose CM, de Voorde MV, Baete K, Cocolios TE, et al. Terbium radionuclides for theranostic applications in nuclear medicine: from atom to bedside. Theranostics. 2024;14:1720–43.38389843 10.7150/thno.92775
3. Naskar N Lahiri S Theranostic Terbium Radioisotopes: challenges in production for clinical application Front Med 2021 8 675014 10.3389/fmed.2021.675014
Naskar N, Lahiri S. Theranostic Terbium Radioisotopes: challenges in production for clinical application. Front Med. 2021;8:675014.10.3389/fmed.2021.675014
4. Collins SM Robinson AP Ivanov P Köster U Cocolios TE Russell B Half-life determination of 155 tb from mass-separated samples produced at CERN-MEDICIS Appl Radiat Isot 2022 190 110480 10.1016/j.apradiso.2022.110480 36209648
Collins SM, Robinson AP, Ivanov P, Köster U, Cocolios TE, Russell B, et al. Half-life determination of 155 tb from mass-separated samples produced at CERN-MEDICIS. Appl Radiat Isot. 2022;190:110480.36209648 10.1016/j.apradiso.2022.110480
5. Collins SM Gilligan C Pierson B Ramirez N Goodwin M Pearce AK Determination of the 161Tb half-life Appl Radiat Isot 2022 182 110140 10.1016/j.apradiso.2022.110140 35158282
Collins SM, Gilligan C, Pierson B, Ramirez N, Goodwin M, Pearce AK, et al. Determination of the 161Tb half-life. Appl Radiat Isot. 2022;182:110140.35158282 10.1016/j.apradiso.2022.110140
6. Müller C Bunka M Haller S Köster U Groehn V Bernhardt P Promising prospects for 44Sc-/47Sc-Based theragnostics: application of 47Sc for Radionuclide Tumor Therapy in mice J Nucl Med 2014 55 1658 64 10.2967/jnumed.114.141614 25034091
Müller C, Bunka M, Haller S, Köster U, Groehn V, Bernhardt P, et al. Promising prospects for 44Sc-/47Sc-Based theragnostics: application of 47Sc for Radionuclide Tumor Therapy in mice. J Nucl Med. 2014;55:1658–64.25034091 10.2967/jnumed.114.141614
7. Domnanich KA Eichler R Müller C Jordi S Yakusheva V Braccini S Production and separation of 43Sc for radiopharmaceutical purposes EJNMMI Radiopharm Chem 2017 2 14 10.1186/s41181-017-0033-9 29503855
Domnanich KA, Eichler R, Müller C, Jordi S, Yakusheva V, Braccini S, et al. Production and separation of 43Sc for radiopharmaceutical purposes. EJNMMI Radiopharm Chem. 2017;2:14.29503855 10.1186/s41181-017-0033-9
8. Blower PJ Lewis JS Zweit J Copper radionuclides and radiopharmaceuticals in nuclear medicine Nucl Med Biol 1996 23 957 80 10.1016/S0969-8051(96)00130-8 9004284
Blower PJ, Lewis JS, Zweit J. Copper radionuclides and radiopharmaceuticals in nuclear medicine. Nucl Med Biol. 1996;23:957–80.9004284 10.1016/S0969-8051(96)00130-8
9. Rösch F Herzog H Qaim SM The beginning and development of the Theranostic Approach in Nuclear Medicine, as exemplified by the Radionuclide Pair 86Y and 90Y Pharmaceuticals 2017 10 56 10.3390/ph10020056 28632200
Rösch F, Herzog H, Qaim SM. The beginning and development of the Theranostic Approach in Nuclear Medicine, as exemplified by the Radionuclide Pair 86Y and 90Y. Pharmaceuticals. 2017;10:56.28632200 10.3390/ph10020056
10. Silberstein EB Radioiodine: the Classic Theranostic Agent Semin Nucl Med 2012 42 164 70 10.1053/j.semnuclmed.2011.12.002 22475425
Silberstein EB. Radioiodine: the Classic Theranostic Agent. Semin Nucl Med. 2012;42:164–70.22475425 10.1053/j.semnuclmed.2011.12.002
11. Aluicio-Sarduy E Hernandez R Olson AP Barnhart TE Cai W Ellison PA Production and in vivo PET/CT imaging of the theranostic pair 132/135La Sci Rep 2019 9 10658 10.1038/s41598-019-47137-0 31337833
Aluicio-Sarduy E, Hernandez R, Olson AP, Barnhart TE, Cai W, Ellison PA, et al. Production and in vivo PET/CT imaging of the theranostic pair 132/135La. Sci Rep. 2019;9:10658.31337833 10.1038/s41598-019-47137-0
12. Crawford JR Robertson AKH Yang H Rodríguez-Rodríguez C Esquinas PL Kunz P Evaluation of 209At as a theranostic isotope for 209At-radiopharmaceutical development using high-energy SPECT Phys Med Biol 2018 63 045025 10.1088/1361-6560/aaaa95 29466250
Crawford JR, Robertson AKH, Yang H, Rodríguez-Rodríguez C, Esquinas PL, Kunz P, et al. Evaluation of 209At as a theranostic isotope for 209At-radiopharmaceutical development using high-energy SPECT. Phys Med Biol. 2018;63:045025.29466250 10.1088/1361-6560/aaaa95
13. McNeil BL Robertson AKH Fu W Yang H Hoehr C Ramogida CF Production, purification, and radiolabeling of the 203Pb/212Pb theranostic pair EJNMMI Radiopharm Chem 2021 6 6 10.1186/s41181-021-00121-4 33527221
McNeil BL, Robertson AKH, Fu W, Yang H, Hoehr C, Ramogida CF, et al. Production, purification, and radiolabeling of the 203Pb/212Pb theranostic pair. EJNMMI Radiopharm Chem. 2021;6:6.33527221 10.1186/s41181-021-00121-4
14. Koniar H Rodríguez-Rodríguez C Radchenko V Yang H Kunz P Rahmim A SPECT imaging of 226Ac as a theranostic isotope for 225Ac radiopharmaceutical development Phys Med Biol 2022 67 185009 10.1088/1361-6560/ac8b5f
Koniar H, Rodríguez-Rodríguez C, Radchenko V, Yang H, Kunz P, Rahmim A, et al. SPECT imaging of 226Ac as a theranostic isotope for 225Ac radiopharmaceutical development. Phys Med Biol. 2022;67:185009.10.1088/1361-6560/ac8b5f
15. Koniar H Wharton L Ingham A Rodríguez-Rodríguez C Kunz P Radchenko V In vivo quantitative SPECT imaging of actinium-226: feasibility and proof-of-concept Phys Med Biol 2024 69 155003 10.1088/1361-6560/ad5c37
Koniar H, Wharton L, Ingham A, Rodríguez-Rodríguez C, Kunz P, Radchenko V, et al. In vivo quantitative SPECT imaging of actinium-226: feasibility and proof-of-concept. Phys Med Biol. 2024;69:155003.10.1088/1361-6560/ad5c37
16. Türler A Matched pair theranostics: Medicinal Chemistry and Chemical Biology highlights Chimia 2019 73 947 10.2533/chimia.2019.947 31753079
Türler A. Matched pair theranostics: Medicinal Chemistry and Chemical Biology highlights. Chimia. 2019;73:947.31753079 10.2533/chimia.2019.947
17. Miller C Rousseau J Ramogida CF Celler A Rahmim A Uribe CF Implications of physics, chemistry and biology for dosimetry calculations using theranostic pairs Theranostics 2022 12 232 59 10.7150/thno.62851 34987643
Miller C, Rousseau J, Ramogida CF, Celler A, Rahmim A, Uribe CF. Implications of physics, chemistry and biology for dosimetry calculations using theranostic pairs. Theranostics. 2022;12:232–59.34987643 10.7150/thno.62851
18. Nica N Nuclear Data sheets for a = 155 Nucl Data Sheets 2019 160 1 404 10.1016/j.nds.2019.100523
Nica N. Nuclear Data sheets for a = 155. Nucl Data Sheets. 2019;160:1–404.10.1016/j.nds.2019.100523
19. Brosch-Lenz J Yousefirizi F Zukotynski K Beauregard J-M Gaudet V Saboury B Role of Artificial Intelligence in Theranostics: toward Routine Personalized Radiopharmaceutical therapies PET Clin 2021 16 627 41 10.1016/j.cpet.2021.06.002 34537133
Brosch-Lenz J, Yousefirizi F, Zukotynski K, Beauregard J-M, Gaudet V, Saboury B, et al. Role of Artificial Intelligence in Theranostics: toward Routine Personalized Radiopharmaceutical therapies. PET Clin. 2021;16:627–41.34537133 10.1016/j.cpet.2021.06.002
20. Wharton L McNeil SW Merkens H Yuan Z Van de Voorde M Engudar G Preclinical evaluation of [155/161Tb]Tb-Crown-TATE—A Novel SPECT Imaging Theranostic Agent Targeting Neuroendocrine Tumours Molecules 2023 28 3155 10.3390/molecules28073155 37049918
Wharton L, McNeil SW, Merkens H, Yuan Z, Van de Voorde M, Engudar G, et al. Preclinical evaluation of [155/161Tb]Tb-Crown-TATE—A Novel SPECT Imaging Theranostic Agent Targeting Neuroendocrine Tumours. Molecules. 2023;28:3155.37049918 10.3390/molecules28073155
21. Bolcaen J Gizawy MA Terry SYA Paulo A Cornelissen B Korde A Marshalling the potential of Auger Electron Radiopharmaceutical Therapy J Nucl Med 2023 64 1344 51 10.2967/jnumed.122.265039 37591544
Bolcaen J, Gizawy MA, Terry SYA, Paulo A, Cornelissen B, Korde A, et al. Marshalling the potential of Auger Electron Radiopharmaceutical Therapy. J Nucl Med. 2023;64:1344–51.37591544 10.2967/jnumed.122.265039
22. Lehenberger S Barkhausen C Cohrs S Fischer E Grünberg J Hohn A The low-energy β – and electron emitter 161Tb as an alternative to 177Lu for targeted radionuclide therapy Nucl Med Biol 2011 38 917 24 10.1016/j.nucmedbio.2011.02.007 21843788
Lehenberger S, Barkhausen C, Cohrs S, Fischer E, Grünberg J, Hohn A, et al. The low-energy β – and electron emitter 161Tb as an alternative to 177Lu for targeted radionuclide therapy. Nucl Med Biol. 2011;38:917–24.21843788 10.1016/j.nucmedbio.2011.02.007
23. Reich CW Nuclear Data sheets for a = 161 Nucl Data Sheets 2011 112 2497 713 10.1016/j.nds.2011.09.001
Reich CW. Nuclear Data sheets for a = 161. Nucl Data Sheets. 2011;112:2497–713.10.1016/j.nds.2011.09.001
24. Müller C van der Meulen NP Schibli R Opportunities and potential challenges of using terbium-161 for targeted radionuclide therapy in clinics Eur J Nucl Med Mol Imaging 2023 50 3181 4 10.1007/s00259-023-06316-y 37436459
Müller C, van der Meulen NP, Schibli R. Opportunities and potential challenges of using terbium-161 for targeted radionuclide therapy in clinics. Eur J Nucl Med Mol Imaging. 2023;50:3181–4.37436459 10.1007/s00259-023-06316-y
25. Dellepiane G Casolaro P Favaretto C Grundler PV Mateu I Scampoli P Cross section measurement of terbium radioisotopes for an optimized 155Tb production with an 18 MeV medical PET cyclotron Appl Radiat Isot 2022 184 110175 10.1016/j.apradiso.2022.110175 35344829
Dellepiane G, Casolaro P, Favaretto C, Grundler PV, Mateu I, Scampoli P, et al. Cross section measurement of terbium radioisotopes for an optimized 155Tb production with an 18 MeV medical PET cyclotron. Appl Radiat Isot. 2022;184:110175.35344829 10.1016/j.apradiso.2022.110175
26. Müller C Fischer E Behe M Köster U Dorrer H Reber J Future prospects for SPECT imaging using the radiolanthanide terbium-155 — production and preclinical evaluation in tumor-bearing mice Nucl Med Biol 2014 41 e58 65 10.1016/j.nucmedbio.2013.11.002 24360901
Müller C, Fischer E, Behe M, Köster U, Dorrer H, Reber J, et al. Future prospects for SPECT imaging using the radiolanthanide terbium-155 — production and preclinical evaluation in tumor-bearing mice. Nucl Med Biol. 2014;41:e58–65.24360901 10.1016/j.nucmedbio.2013.11.002
27. Favaretto C Talip Z Borgna F Grundler PV Dellepiane G Sommerhalder A Cyclotron production and radiochemical purification of terbium-155 for SPECT imaging EJNMMI Radiopharm Chem 2021 6 37 10.1186/s41181-021-00153-w 34778932
Favaretto C, Talip Z, Borgna F, Grundler PV, Dellepiane G, Sommerhalder A, et al. Cyclotron production and radiochemical purification of terbium-155 for SPECT imaging. EJNMMI Radiopharm Chem. 2021;6:37.34778932 10.1186/s41181-021-00153-w
28. Müller C Umbricht CA Gracheva N Tschan VJ Pellegrini G Bernhardt P Terbium-161 for PSMA-targeted radionuclide therapy of prostate cancer Eur J Nucl Med Mol Imaging 2019 46 1919 30 10.1007/s00259-019-04345-0 31134301
Müller C, Umbricht CA, Gracheva N, Tschan VJ, Pellegrini G, Bernhardt P, et al. Terbium-161 for PSMA-targeted radionuclide therapy of prostate cancer. Eur J Nucl Med Mol Imaging. 2019;46:1919–30.31134301 10.1007/s00259-019-04345-0
29. Borgna F Barritt P Grundler PV Talip Z Cohrs S Zeevaart JR Simultaneous visualization of 161Tb- and 177Lu-Labeled somatostatin analogues using dual-isotope SPECT Imaging Pharmaceutics 2021 13 536 10.3390/pharmaceutics13040536 33921467
Borgna F, Barritt P, Grundler PV, Talip Z, Cohrs S, Zeevaart JR, et al. Simultaneous visualization of 161Tb- and 177Lu-Labeled somatostatin analogues using dual-isotope SPECT Imaging. Pharmaceutics. 2021;13:536.33921467 10.3390/pharmaceutics13040536
30. Wharton L, McNeil SW, Zhang C, Engudar G, Van de Voorde M, Zeisler J et al. Preclinical evaluation of MC1R targeting theranostic pair [155Tb]Tb-crown-αMSH and [161Tb]Tb-crown-αMSH. Nucl Med Biol. 2024;136–7:108925.
31. Cassells I Ahenkorah S Burgoyne AR Van de Voorde M Deroose CM Cardinaels T Radiolabeling of Human serum albumin with Terbium-161 using mild conditions and evaluation of in vivo Stability Front Med 2021 8 675122 10.3389/fmed.2021.675122
Cassells I, Ahenkorah S, Burgoyne AR, Van de Voorde M, Deroose CM, Cardinaels T, et al. Radiolabeling of Human serum albumin with Terbium-161 using mild conditions and evaluation of in vivo Stability. Front Med. 2021;8:675122.10.3389/fmed.2021.675122
32. McNeil SW Van de Voorde M Zhang C Ooms M Bénard F Radchenko V A simple and automated method for 161Tb purification and ICP-MS analysis of 161Tb EJNMMI Radiopharm Chem 2022 7 31 10.1186/s41181-022-00183-y 36459299
McNeil SW, Van de Voorde M, Zhang C, Ooms M, Bénard F, Radchenko V, et al. A simple and automated method for 161Tb purification and ICP-MS analysis of 161Tb. EJNMMI Radiopharm Chem. 2022;7:31.36459299 10.1186/s41181-022-00183-y
33. Marin I Rydèn T Van Essen M Svensson J Gracheva N Köster U Establishment of a clinical SPECT/CT protocol for imaging of 161Tb EJNMMI Phys 2020 7 45 10.1186/s40658-020-00314-x 32613587
Marin I, Rydèn T, Van Essen M, Svensson J, Gracheva N, Köster U, et al. Establishment of a clinical SPECT/CT protocol for imaging of 161Tb. EJNMMI Phys. 2020;7:45.32613587 10.1186/s40658-020-00314-x
34. Baum RP Singh A Kulkarni HR Bernhardt P Rydén T Schuchardt C First-in-humans application of 161Tb: a feasibility study using 161Tb-DOTATOC J Nucl Med 2021 62 1391 10.2967/jnumed.120.258376 33547209
Baum RP, Singh A, Kulkarni HR, Bernhardt P, Rydén T, Schuchardt C, et al. First-in-humans application of 161Tb: a feasibility study using 161Tb-DOTATOC. J Nucl Med. 2021;62:1391.33547209 10.2967/jnumed.120.258376
35. Rosar F Maus S Schaefer-Schuler A Burgard C Khreish F Ezziddin S New Horizons in Radioligand Therapy: 161Tb-PSMA-617 in Advanced mCRPC Clin Nucl Med 2023 48 433 4 10.1097/RLU.0000000000004589 36758550
Rosar F, Maus S, Schaefer-Schuler A, Burgard C, Khreish F, Ezziddin S. New Horizons in Radioligand Therapy: 161Tb-PSMA-617 in Advanced mCRPC. Clin Nucl Med. 2023;48:433–4.36758550 10.1097/RLU.0000000000004589
36. Al-Ibraheem A Doudeen RM Juaidi D Abufara A Maus S 161Tb-PSMA Radioligand Therapy: first-in-human SPECT/CT imaging J Nucl Med 2023 64 1322 3 10.2967/jnumed.122.265291 36759201
Al-Ibraheem A, Doudeen RM, Juaidi D, Abufara A, Maus S. 161Tb-PSMA Radioligand Therapy: first-in-human SPECT/CT imaging. J Nucl Med. 2023;64:1322–3.36759201 10.2967/jnumed.122.265291
37. Buteau JP Kostos LK Alipour R Jackson P McIntosh L Emmerson B VIOLET: a phase I/II trial evaluation of radioligand treatment in men with metastatic castration-resistant prostate cancer with [161Tb]Tb-PSMA-I&T J Clin Oncol 2023 41 TPS281 281 10.1200/JCO.2023.41.6_suppl.TPS281
Buteau JP, Kostos LK, Alipour R, Jackson P, McIntosh L, Emmerson B, et al. VIOLET: a phase I/II trial evaluation of radioligand treatment in men with metastatic castration-resistant prostate cancer with [161Tb]Tb-PSMA-I&T. J Clin Oncol. 2023;41:TPS281–281.10.1200/JCO.2023.41.6_suppl.TPS281
38. McIntosh L Jackson P Emmerson B Buteau JP Alipour R Kong G Quantitative calibration of Tb-161 SPECT/CT in view of personalised dosimetry assessment studies EJNMMI Phys 2024 11 18 10.1186/s40658-024-00611-9 38372952
McIntosh L, Jackson P, Emmerson B, Buteau JP, Alipour R, Kong G, et al. Quantitative calibration of Tb-161 SPECT/CT in view of personalised dosimetry assessment studies. EJNMMI Phys. 2024;11:18.38372952 10.1186/s40658-024-00611-9
39. Alcocer-Ávila ME Ferreira A Quinto MA Morgat C Hindié E Champion C Radiation doses from 161Tb and 177Lu in single tumour cells and micrometastases EJNMMI Phys 2020 7 33 10.1186/s40658-020-00301-2 32430671
Alcocer-Ávila ME, Ferreira A, Quinto MA, Morgat C, Hindié E, Champion C. Radiation doses from 161Tb and 177Lu in single tumour cells and micrometastases. EJNMMI Phys. 2020;7:33.32430671 10.1186/s40658-020-00301-2
40. Kunz P, Andreoiu C, Brown V, Cervantes M, Even J, Garcia FH et al. Medical isotope collection from ISAC targets. EPJ Web Conf. 2020. p. 06003.
41. Fiaccabrino DE Kunz P Radchenko V Potential for production of medical radionuclides with on-line isotope separation at the ISAC facility at TRIUMF and particular discussion of the examples of 165Er and 155Tb Nucl Med Biol 2021 94–95 81 91 10.1016/j.nucmedbio.2021.01.003 33607326
Fiaccabrino DE, Kunz P, Radchenko V. Potential for production of medical radionuclides with on-line isotope separation at the ISAC facility at TRIUMF and particular discussion of the examples of 165Er and 155Tb. Nucl Med Biol. 2021;94–95:81–91.33607326 10.1016/j.nucmedbio.2021.01.003
42. Kunz P. ISAC Yield Database [Internet]. 2023. https://yield.targets.triumf.ca/search/yield/data
43. Van de Voorde M Van Hecke K Cardinaels T Binnemans K Radiochemical processing of nuclear-reactor-produced radiolanthanides for medical applications Coord Chem Rev 2019 382 103 25 10.1016/j.ccr.2018.11.007
Van de Voorde M, Van Hecke K, Cardinaels T, Binnemans K. Radiochemical processing of nuclear-reactor-produced radiolanthanides for medical applications. Coord Chem Rev. 2019;382:103–25.10.1016/j.ccr.2018.11.007
44. Crawford JR Yang H Kunz P Wilbur DS Schaffer P Ruth TJ Development of a preclinical 211Rn/211At generator system for targeted alpha therapy research with 211At Nucl Med Biol 2017 48 31 5 10.1016/j.nucmedbio.2017.01.011 28193502
Crawford JR, Yang H, Kunz P, Wilbur DS, Schaffer P, Ruth TJ. Development of a preclinical 211Rn/211At generator system for targeted alpha therapy research with 211At. Nucl Med Biol. 2017;48:31–5.28193502 10.1016/j.nucmedbio.2017.01.011
45. Goorden MC van der Have F Kreuger R Ramakers RM Vastenhouw B Burbach JPH VECTor: a preclinical imaging system for simultaneous submillimeter SPECT and PET J Nucl Med 2013 54 306 10.2967/jnumed.112.109538 23077113
Goorden MC, van der Have F, Kreuger R, Ramakers RM, Vastenhouw B, Burbach JPH, et al. VECTor: a preclinical imaging system for simultaneous submillimeter SPECT and PET. J Nucl Med. 2013;54:306.23077113 10.2967/jnumed.112.109538
46. Branderhorst W Vastenhouw B Beekman FJ Pixel-based subsets for rapid multi-pinhole SPECT reconstruction Phys Med Biol 2010 55 2023 34 10.1088/0031-9155/55/7/015 20299722
Branderhorst W, Vastenhouw B, Beekman FJ. Pixel-based subsets for rapid multi-pinhole SPECT reconstruction. Phys Med Biol. 2010;55:2023–34.20299722 10.1088/0031-9155/55/7/015
47. Vastenhouw B Beekman F Submillimeter Total-Body murine imaging with U-SPECT-I J Nucl Med 2007 48 487 93 17332628
Vastenhouw B, Beekman F. Submillimeter Total-Body murine imaging with U-SPECT-I. J Nucl Med. 2007;48:487–93.17332628
48. Miwa K Inubushi M Takeuchi Y Katafuchi T Koizumi M Saga T Performance characteristics of a novel clustered multi-pinhole technology for simultaneous high-resolution SPECT/PET Ann Nucl Med 2015 29 460 6 10.1007/s12149-015-0966-6 25894055
Miwa K, Inubushi M, Takeuchi Y, Katafuchi T, Koizumi M, Saga T, et al. Performance characteristics of a novel clustered multi-pinhole technology for simultaneous high-resolution SPECT/PET. Ann Nucl Med. 2015;29:460–6.25894055 10.1007/s12149-015-0966-6
49. Ivashchenko O van der Have F Goorden MC Ramakers RM Beekman FJ Ultra-high-sensitivity Submillimeter Mouse SPECT J Nucl Med 2015 56 470 5 10.2967/jnumed.114.147140 25678487
Ivashchenko O, van der Have F, Goorden MC, Ramakers RM, Beekman FJ. Ultra-high-sensitivity Submillimeter Mouse SPECT. J Nucl Med. 2015;56:470–5.25678487 10.2967/jnumed.114.147140
50. Robertson AKH Ramogida CF Rodríguez-Rodríguez C Blinder S Kunz P Sossi V Multi-isotope SPECT imaging of the 225Ac decay chain: feasibility studies Phys Med Biol 2017 62 4406 20 10.1088/1361-6560/aa6a99 28362640
Robertson AKH, Ramogida CF, Rodríguez-Rodríguez C, Blinder S, Kunz P, Sossi V, et al. Multi-isotope SPECT imaging of the 225Ac decay chain: feasibility studies. Phys Med Biol. 2017;62:4406–20.28362640 10.1088/1361-6560/aa6a99
51. Wu C van der Have F Vastenhouw B Dierckx RAJO Paans AMJ Beekman FJ Absolute quantitative total-body small-animal SPECT with focusing pinholes Eur J Nucl Med Mol Imaging 2010 37 2127 35 10.1007/s00259-010-1519-9 20577738
Wu C, van der Have F, Vastenhouw B, Dierckx RAJO, Paans AMJ, Beekman FJ. Absolute quantitative total-body small-animal SPECT with focusing pinholes. Eur J Nucl Med Mol Imaging. 2010;37:2127–35.20577738 10.1007/s00259-010-1519-9
52. Walker MD Goorden MC Dinelle K Ramakers RM Blinder S Shirmohammad M Performance Assessment of a preclinical PET scanner with Pinhole Collimation by comparison to a coincidence-based small-animal PET scanner J Nucl Med 2014 55 1368 10.2967/jnumed.113.136663 24904110
Walker MD, Goorden MC, Dinelle K, Ramakers RM, Blinder S, Shirmohammad M, et al. Performance Assessment of a preclinical PET scanner with Pinhole Collimation by comparison to a coincidence-based small-animal PET scanner. J Nucl Med. 2014;55:1368.24904110 10.2967/jnumed.113.136663
53. Hutton BF Buvat I Beekman FJ Review and current status of SPECT scatter correction Phys Med Biol 2011 56 R85 10.1088/0031-9155/56/14/R01 21701055
Hutton BF, Buvat I, Beekman FJ. Review and current status of SPECT scatter correction. Phys Med Biol. 2011;56:R85.21701055 10.1088/0031-9155/56/14/R01
