<figure><img src="image_1.png" /></figure>

HAL
open science

Implementation of patient dosimetry in the clinical
practice after targeted radiotherapy using
[177Lu-[DOTA0, Tyr3]-octreotate

Lore Santoro, Erick Mora-Ramirez, Dorian Trauchessec, Soufiane Chouaf, Pierre Eustache, Jean-Pierre Pouget, Pierre-Olivier Kotzki, Manuel Bardiès, Emmanuel Deshayes

<figure><img src="image_2.png" /></figure>

To cite this version:

Lore Santoro, Erick Mora-Ramirez, Dorian Trauchessec, Soufiane Chouaf, Pierre Eustache, et al. Implementation of patient dosimetry in the clinical practice after targeted radiotherapy using [177Lu-[DOTA0, Tyr3]-octreotate. EJNMMI Research, 2018, 8, pp.103. 10.1186/s13550-018-0459-4. hal-02292273

HAL Id: hal-02292273

https://hal.umontpellier.fr/hal-02292273

Submitted on 19 Sep 2019

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.
ORIGINAL RESEARCH

Open Access

<figure><img src="image_3.png" /></figure>

CrossMark

# Implementation of patient dosimetry in the clinical practice after targeted radiotherapy using [¹⁷⁷Lu-[DOTA0, Tyr3]-octreotate

Lore Santoro<fnref n="1" />*, Erick Mora-Ramirez²,³,⁴, Dorian Trauchessec<fnref n="1" />, Soufiane Chouaf<fnref n="1" />, Pierre Eustache<fnref n="1" />, Jean-Pierre Pouget⁵, Pierre-Olivier Kotzki<fnref n="1" />,⁵, Manuel Bardiès²,³ and Emmanuel Deshayes<fnref n="1" />,⁵

**Abstract**

**Background:** This study's aim was to develop our dosimetric methodology using a commercial workstation for the routine evaluation of the organs at risk during peptide receptor radionuclide therapy (PRRT) with ¹⁷⁷Lu.

**Methods:** First, planar and SPECT sensitivity factors were determined on phantoms. The reconstruction parameters were optimized by SPECT/CT image acquisition using a NEMA IEC phantom containing a 500 ml bottle of ¹⁷⁷Lu, to simulate a kidney. The recovery coefficients were determined on various phantoms. For the red marrow, this was calculated using a NEMA IEC phantom that contained a centrally placed bottle of 80 ml of ¹⁷⁷Lu (to model the L2-L4 red marrow) flanked by two 200 ml bottles with ¹⁷⁷Lu to simulate the kidneys. Then, SPECT/CT images were acquired at 4, 24, 72, and 192 h after injection in 12 patients with neuroendocrine tumors who underwent PRRT with ¹⁷⁷Lu-DOTATATE. SPECT data were reconstructed using the iterative ordered subset expectation maximization (OSEM) method, with six iterations and ten subsets, attenuation, scatter, recovery resolution corrections, and a Gaussian post-filter of 0.11 cm. The liver, spleen, kidneys, and red marrow dose per administered activity (AD/A admin) values were calculated with the Medical Internal Radiation Dose (MIRD) formalism and the residence times (Dosimetry toolkit® application) using standard and CT imaging-based organ masses (OLINDA/EXM® V1.0 software).

**Results:** Sensitivity factors of 6.11 ± 0.01 and 5.67 ± 0.08 counts/s/MBq were obtained with planar and SPECT/CT acquisitions, respectively. A recovery coefficient of 0.78 was obtained for the modeled L2–L4 red marrow. The mean AD/A admin values were 0.43 ± 0.13 mGy/MBq [0.27–0.91] for kidneys, 0.54 ± 0.58 mGy/MBq [0.12–2.26] for liver, 0.61 ± 0.13 mGy/MBq [0.42–0.89] for spleen, and 0.04 ± 0.02 mGy/MBq [0.01–0.09] for red marrow. The AD/A admin values varied when calculated using the personalized and standard organ mass, particularly for kidneys (p = 1 × 10⁻⁷), spleen (p = 0.0069), and red marrow (p = 0.0027). Intra-patient differences were observed especially in organs close to or including tumor cells or metastases.

**Conclusions:** The obtained AD/A admin values were in agreement with the literature data. This study shows the technical feasibility of patient dosimetry in clinical practice and the need to obtain patient-specific information.

**Keywords:** Peptide receptor radionuclide therapy, [¹⁷⁷Lu- [DOTA0, Tyr3]-octreotate, Medical internal radiation dose, Patient-specific dosimetry

<figure><img src="image_4.png" /></figure>

<footnote>* Correspondence: lore.santoro@icm.unicancer.fr</footnote>

<footnote n="1">¹Nuclear Medicine Department, Montpellier Cancer Institute (ICM), University of Montpellier, 208 Avenue des Apothicaires, 34298 Montpellier Cedex5, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
## Background

In recent years, new radiopharmaceuticals, such as radio-peptides, have been developed for targeted radiotherapy. Consequently, in addition to red marrow, healthy liver and kidneys also are now identified as organs at risk (OAR) [1, 2]. Moreover, the increasing evidence that treatment outcomes correlate with the absorbed doses delivered to tumors and healthy organs [3, 4] makes individualized dosimetry all the more necessary [5].

Peptide receptor radionuclide therapy (PRRT) is a promising treatment for patients with neuroendocrine tumors (NET). Recently, the NETTER-1 trial [6] showed that treatment with ¹⁷⁷Lu-DOTATATE results in markedly longer progression-free survival and higher response rate than treatment with high-dose octreotide alone in patients with advanced midgut NET. In September 2017, ¹⁷⁷Lu-DOTATATE (LUTATHERA®) was approved for this indication by the European Medicines Agency. It has been shown that ¹⁷⁷Lu-DOTATATE uptake in kidneys and tumors greatly varies among patients [7], and bone marrow and kidneys are considered as dose-limiting organs. Although no clear cumulative absorbed dose cut-off has been identified in PRRT with ¹⁷⁷Lu to predict the risk of organ failure, which is probably multifactorial [4, 8], 23 Gy for kidneys and 2 Gy for bone marrow are sometimes proposed, although these values are based on studies with fractionated external beam radiotherapy (EBRT) that has different physical and radiobiological mechanisms. For example, based on fractionated EBRT, a 5% risk of renal dysfunction at 5 years has been described for a mean absorbed dose of 18–23 Gy and 0.5–1.25 Gy/fraction [9, 10].

In PRRT with ¹⁷⁷Lu, the important inter-patient variations in peptide pharmacokinetics require treatment individualization by tailoring the number of cycles or the administered activity [11]. In this context, therapy planning based on the maximum tolerable absorbed dose to non-target organs (“as high as safely attainable”, AHASA approach) [12] could be considered, instead of the “as low as reasonably achievable” (ALARA) approach (i.e., the dose to non-target tissues should be reasonably low). However, currently, a fixed activity of 7.4 GBq per cycle, as described in the NETTER-1 trial, is usually administered. Hence, personalized dosimetry is often performed mainly to ensure safety and evaluate the absorbed dose to the tumor rather than to optimize the administered activity and to assess the dose-response relationship.

Patient dosimetry requires the accurate estimation of the activity in the targeted organs at several time points [13, 14]. Therefore, the preliminary calibration and quantification steps are crucial [15, 16]. Although there are Medical Internal Radiation Dose (MIRD) guidelines [17, 18], a standardized dosimetry protocol to evaluate safety and toxicity, and to perform dosimetric

evaluations is crucially required. The need to determine the absorbed doses delivered to kidneys after ¹⁷⁷Lu-DOTATATE treatment was an evidence for nuclear medicine physicians in our department [11]. Therefore, before starting the first treatment with ¹⁷⁷Lu-DOTATATE, we implemented the imaging protocol based on the MIRD pamphlet n°26 [18], and performed a preliminary study on phantoms. Then, we defined the imaging schedule for patient dosimetry based on published data and our own institute logistics. We performed the dosimetry analyses using the MIRD formalism and the tools available in our nuclear medicine department. Using this protocol for dosimetry after ¹⁷⁷Lu-DOTATATE treatment and a commercial workstation, we could calculate the absorbed doses per unit of administered activity (AD/A admin) in the OARs (kidneys, liver, spleen, and red marrow) in patients with NET.

## Material and methods

### SPECT/CT equipment

All imaging acquisitions were performed using the SPECT/CT Discovery NM/CT 670 system (General Electric [GE] Healthcare) with a Bright Speed 16 CT scanner and 3/8-in NaI(Tl) crystal thickness. Nuclear medicine images were acquired using a medium-energy general purpose parallel-hole collimator. A 20% energy window centered on the 208 keV photopeak and a 10% scatter correction window centered on 177 keV were applied [18]. A 128 × 128 pixel matrix was used.

### Software

The Dosimetry Toolkit®, an application of the Xeleris® software provided with the SPECT/CT Discovery NM/CT 670 system (GE Healthcare) [19], was used to determine the radiopharmaceutical residence time in segmented organs with the multi-SPECT/CT scenario, as described by Kupitz et al. [20]. First, the application “Preparation for dosimetry toolkit” was used for SPECT/CT raw data reconstruction and CT data registration. Then, the “Dosimetry Toolkit” application was used to segment the different organs, to create the time activity curves fitted by a mono-exponential function, and to calculate the residence time.

The OLINDA/EXM® V1.0 software [21] was used to calculate the organ absorbed and effective doses. It contains a large pre-established database of radionuclides and S factors to calculate the absorbed dose per unit of administered activity. Standard or patient-adapted organ masses can be used with this software. The residence times obtained for each organ with Dosimetry Toolkit® were used as input data.

## Preliminary study on phantoms

### Correction maps

Energy and uniformity correction maps were computed using static data, acquired without collimator, of a
syringe containing 91.8 MBq of <sup>177</sup>Lu placed at a distance that corresponded to five times the field of view. A 10% energy window centered on the 208 keV photopeak was applied. A 256 × 256 pixel matrix and 10,000 kcounts were used for the energy and uniformity correction maps, and a 512 × 512 pixel matrix and 60,000 kcounts for the uniformity correction map.

## Sensitivity factors

The system sensitivity was evaluated with two methods (see below). The sensitivity factor (expressed in counts/ s/MBq) allowed converting the numbers of events (counts) detected by the gamma camera into activity values. For each method, the used activity was accurately measured with a CRC-25R from Capintec dose calibrator (Berthold Technologies). The calibration factor was calculated using a calibration vial of <sup>177</sup>Lu (Advanced Accelerator Applications, Saint Genis Pouilly, France), with a maximum activity measurement error of 5%. The clock was synchronized with the gamma camera clock to allow accurate decay correction using a half-life of 6.647 days.

**Method 1 (planar sensitivity factor)** According to the manufacturer's recommendations (GE Healthcare), a syringe of <sup>177</sup>Lu containing an activity of 59.1 ± 2.9 MBq was placed on the examination table of the gamma camera at 10 cm from the detectors. Planar imaging was performed for 5 min to obtain the number of events detected by the detectors. The geometric mean was calculated and reported to the real activity within the syringe and the acquisition time.

Then, a 16 ml hollow sphere filled with <sup>177</sup>Lu (75.8 ± 3.8 MBq activity) was placed in the air between the detectors. Planar acquisitions were performed at different distances from the detectors (8, 13, and 18 cm) for 5 min to determine the system sensitivity.

**Method 2 (SPECT sensitivity factor)** A NEMA IEC phantom (Body Phantom NU2-2001/2007) was filled with non-radioactive water. A 500 ml bottle to simulate a kidney was filled with 0.54 ± 0.03 MBq/ml of <sup>177</sup>Lu solution and fixed inside the phantom. SPECT/CT images were acquired using the body contour option, rotation of 180° per detector, total of 60 projections, and 120 s/each. For attenuation correction, CT images were acquired (120 kV, 200 mA, slice thickness: 1.25 mm, rotation time: 0.8 s, pitch: 1.375, 512 × 512 pixel matrix), with standard reconstruction. For SPECT acquisitions, the ordered subset expectation maximization (OSEM) iterative reconstruction algorithm was used, with optimized reconstruction parameters.

## Optimization of the reconstruction parameters

The reconstruction parameters were optimized using the SPECT/CT images acquired with the previously used NEMA IEC phantom. Using the “Preparation for Dosimetry toolkit” application, SPECT data were reconstructed with different numbers of iterations and subsets, with or without the manufacturer’s corrections (scatter, CT-based attenuation correction, resolution recovery), and different Gaussian post-filters. With the “Dosimetry toolkit” application, regions of interest were segmented in the CT slices. Information about the radionuclide and the previously computed planar sensitivity factor were entered in the appropriate interface. The segmented volume and internal activity were obtained for each evaluated reconstruction, and compared with those expected for the phantom.

## Recovery coefficients

The recovery coefficients (i.e., the ratio between the activity concentration estimated from the image and the true activity concentration in the object) [17] were evaluated using two SPECT/CT acquisitions performed using a Deluxe Jaszczak phantom. For both acquisitions, six hollow spheres of 0.5 ml, 1.0 ml, 2.0 ml, 4.0 ml, 8.0 ml, and 16.0 ml, containing 0.537 ± 0.028 MBq/ml of <sup>177</sup>Lu, were placed inside the phantom. For the first acquisition, the phantom was filled with non-radioactive water, and for the second one, with 474.3 ± 23.7 MBq of <sup>177</sup>Lu (i.e., a background concentration of 0.074 ± 0.004 MBq/ml). The SPECT/CT acquisition parameters were the same as those used with the NEMA IEC phantom.

The recovery coefficient of a large volume (500 ml) was estimated by SPECT/CT using the previously described NEMA IEC phantom.

A second NEMA IEC phantom (Body Phantom NU2-2001/2007) that contained a central bottle of 80 ml with 2.2 ± 0.11 MBq of <sup>177</sup>Lu (3.5 cm × 3.5 cm × 6.5 cm; to model the red marrow at the L2–L4 level) flanked by two 200 ml bottles (with 73.7 ± 3.7 MBq of <sup>177</sup>Lu/each, 30-fold higher activity to simulate the kidneys) was also used. The background was filled with non-radioactive water. SPECT/CT images were acquired using the body contour option, rotation of 180° per detector, 60 projections, and 45 s/each. For attenuation correction, CT images were acquired (120 kV, automatic mA regulation with a max at 200 mA, a noise index at 6.43, slice thickness of 5 mm, rotation time of 0.8 s, pitch 1.375, 512 × 512 pixels matrix), with standard reconstruction.

For all SPECT acquisitions, the OSEM iterative reconstruction method algorithm was used with the reconstruction parameters selected after optimization (see previous section).

## Clinical dosimetry method application

### Patients and treatment

Between May 2016 and February 2018, 12 patients (10 men and 2 women) with a NET who underwent PRRT with <sup>177</sup>Lu-[DOTA0, Tyr3]-octreotate (LUTATHERA®;
Advanced Accelerator Applications, Saint Genis Pouilly, France) were included in this dosimetric evaluation. Table 1 shows the patients' characteristics. The study was approved by the local ethics review board.

PRRT consisted in one intravenous injection of 7.4 GBq LUTATHERA® every 8 weeks for a total of four cycles. Lysine and arginine were administered concomitantly to ensure kidney protection by reducing the tubular reabsorption of the radiolabeled peptides. All patients were hospitalized in special radioprotection rooms for 24 h after injection. After injection, the residual activity in the vial was measured with the Capintec CRC-25R dose calibrator. By taking into account the ¹⁷⁷Lu physical decay, the real administered activity was determined by subtracting the residual activity in the vial from the activity before injection.

## Dosimetry imaging protocol

Dosimetry calculations were based on the imaging data collected after the first two cycles. Images were acquired with the SPECT/CT Discovery NM/CT 670, at 4 h, 24 h, 72 h, and 192 h after the first and second injection. In total, 60 projections (45 s per projection) were acquired with a 128 × 128 pixel matrix (pixel size: 4.416 mm). For medical consideration, a whole body (WB) scan was performed at 72 h with a scan speed of 15 cm/min. At 4 h post-injection, CT images were

acquired with a better image quality (120 kV, automatic mA regulation with max = 200 mA, noise index of 6.43, slice thickness of 5 mm, rotation time of 0.8 s, pitch 1.375, 512 × 512 pixel matrix, and standard reconstruction filter). For this first CT scan, patient set-up and immobilization devices were recorded to reproducibly position the patients for the next SPECT/CT image acquisitions. The other CT scans were acquired with the same parameters, except for the rotation time (0.6 s) and 80 mA fixed. After the third and fourth cycles, SPECT/CT and WB scan were performed only at 24 h.

## Dosimetry calculations

Following the last SPECT/CT acquisition, at 192 h after cycle 2, all SPECT/CT data were loaded on the “Preparation for dosimetry toolkit” application. Transversal slices were reconstructed using the OSEM algorithm with the reconstruction settings defined in the optimization study associated with correction of the patient's movements. Each SPECT acquisition was registered with the corresponding CT acquisition. A rigid registration between CT scans was performed. The reconstruction results were loaded on the “Dosimetry toolkit” application. The OARs (red marrow, kidneys, liver, and spleen) were segmented using the CT and SPECT images collected at 4 h post-injection. They were then replicated for the 24 h, 72 h, and 192 h images. The red

<table><caption>**Table 1** Patients' characteristics</caption><thead><tr><th rowspan="2">Patient</th><th rowspan="2">Sex</th><th rowspan="2">Age (years)</th><th rowspan="2">Weight (kg)</th><th rowspan="2">Primary tumor</th><th rowspan="2">Metastases</th><th colspan="2">Injected activity (MBq)</th><th rowspan="2">Number of treatment cycles</th><th colspan="2">SPECT/CT acquisition times</th></tr><tr><th>C1</th><th>C2</th><th>C1</th><th>C2</th></tr></thead><tbody><tr><td>1</td><td>M</td><td>75</td><td>74</td><td>Small intestine NET</td><td>Nodes, mesentery, liver</td><td>7167</td><td>7073</td><td>4</td><td>NA</td><td>4 h, 24 h, 72 h, 192 h</td></tr><tr><td>2</td><td>M</td><td>59</td><td>72</td><td>Small intestine NET</td><td>Nodes, liver, bone</td><td>7287</td><td>†</td><td>1</td><td>4 h, 24 h, 168 h</td><td>†</td></tr><tr><td>3</td><td>M</td><td>82</td><td>70</td><td>Small intestine NET</td><td>Nodes, mesentery, bone</td><td>7180</td><td>6642</td><td>2</td><td>4 h, 24 h, 72 h</td><td>4 h, 24 h, 72 h, 192 h</td></tr><tr><td>4</td><td>M</td><td>71</td><td>71</td><td>Pancreas NET</td><td>Liver</td><td>7054</td><td>7134</td><td>4</td><td>4 h, 24 h, 72 h, 192 h</td><td>4 h, 24 h, 192 h</td></tr><tr><td>5</td><td>F</td><td>63</td><td>56</td><td>Pancreas NET</td><td>Liver</td><td>7323</td><td>7071</td><td>4</td><td>4 h, 24 h, 72 h, 192 h</td><td>4 h, 24 h, 72 h</td></tr><tr><td>6</td><td>M</td><td>59</td><td>79</td><td>Small intestine NET</td><td>Mesentery, liver</td><td>7207</td><td>7188</td><td>4</td><td>4 h, 24 h, 72 h, 192 h</td><td>4 h, 24 h, 72 h, 192 h</td></tr><tr><td>7</td><td>F</td><td>82</td><td>57</td><td>Pancreas NET</td><td>Liver</td><td>7177</td><td>7239</td><td>4</td><td>4 h, 24 h, 72 h, 192 h</td><td>4 h, 24 h, 72 h, 192 h</td></tr><tr><td>8</td><td>M</td><td>61</td><td>74</td><td>Small intestine NET</td><td>Nodes, mesentery, liver, bone</td><td>7298</td><td>7210</td><td>4</td><td>4 h, 24 h, 72 h, 192 h</td><td>4 h, 24 h, 72 h, 192 h</td></tr><tr><td>9</td><td>M</td><td>70</td><td>88</td><td>Small intestine NET</td><td>Liver, bone</td><td>7222</td><td>7158</td><td>4</td><td>4 h, 24 h, 192 h</td><td>4 h, 24 h, 72 h, 192 h</td></tr><tr><td>10</td><td>M</td><td>78</td><td>89</td><td>Large intestine NET</td><td>Nodes, liver</td><td>7162</td><td>6620</td><td>2</td><td>4 h, 24 h, 72 h, 192 h</td><td>4 h, 24 h, 192 h</td></tr><tr><td>11</td><td>M</td><td>53</td><td>70</td><td>Small intestine NET</td><td>Liver, bone</td><td>7102</td><td>7260</td><td>3</td><td>4 h, 24 h, 72 h, 192 h</td><td>4 h, 24 h, 72 h, 192 h</td></tr><tr><td>12</td><td>M</td><td>74</td><td>70</td><td>Small intestine NET</td><td>Nodes, mesentery</td><td>7288</td><td>7578</td><td>3</td><td>4 h, 24 h, 72 h, 192 h</td><td>4 h, 24 h, 72 h, 192 h</td></tr></tbody></table>

NA not available, C cycle, NET neuroendocrine tumor, † patient dead
marrow absorbed dose was determined by delineating the trabecular section on the L2–L4 lumbar vertebrae, considering that this section represents 6.7% of the total red marrow [22]. For kidneys, the delineation encompassed the cortex and the medulla regions of the left and right kidney, as described by Sundlöv et al. [11]. For the kidney volume, the partial volume effect correction was considered negligible. For the trabecular section on the L2–L4 lumbar vertebrae, the obtained recovery coefficient was applied.

The residence times were calculated with the SPECT sensitivity factor, and then exported to OLINDA/EXM® V1.0 to calculate the AD/A admin values using the standard organ masses included in the software and the personalized organ masses. These were determined using the volume of each organ of interest defined on the CT images and the biological tissue density proposed by the Monte Carlo Gate database and by Vieira et al. [23].

The relative difference (in %) between the personalized and standard mass for organ i and patient j ($\Delta_{\text{mass i, j}}$) was defined as follows:

$$ \Delta_{\text{mass i, j}} = 100 \times \frac{\text{Personalized organ mass}_{ij} - \text{Standard organ mass}_{ij}}{\text{Standard organ mass}_{ij}} $$

where the personalized organ mass$_{ij}$ is the mass of organ i for patient j calculated from the CT images and the biological tissue density, and the standard organ mass$_{ij}$ is the mass of organ i for patient j included in the OLINDA/EXM®V1.0 software.

The AD/A admin values obtained with personalized and standard organ masses for each patient (calculated after cycle 1 and cycle 2) were compared with the paired Student's t test (n = 22). To study the intra-patient variability, the personalized absorbed doses after cycle 1 and 2 (n = 10) were compared using the paired Student's t test.

The difference (in %), between the organ i AD/A admin after cycle 1 and 2 and for the patient j ($\Delta_{\text{AD/Aadmin i, j}}$), was defined as follows:

$$ \Delta_{\text{AD/A admin i}} = 100 \times \frac{SD(\text{AD/A admin}_{1ij}, \text{AD/A admin}_{2ij})}{\text{Arithmetic mean}(\text{AD/A admin}_{1ij}, \text{AD/A admin}_{2ij})} $$

where SD is the standard deviation, and AD/A admin$_{1ij}$ and AD/A admin$_{2ij}$ are the organ i AD/A admin for patient j after cycle 1 and cycle 2, respectively.

$\bar{\Delta}_{\text{AD/A admin i}}$ (in %) represents the mean difference of the organ i AD/A admin values for all patients (j = 1 to j = 10) and was calculated as follows:

$$ \bar{\Delta}_{\text{AD/A admin i}} = \frac{1}{j} \sum_{j=1}^{j=10} \Delta_{\text{AD ij}} $$

# Results

## Preliminary imaging optimization/calibration with phantoms

The obtained planar sensitivity factor of 6.11 ± 0.01 counts/s/MBq did not vary significantly with the distance from the detectors (Fig. 1) and was then used to evaluate the optimized reconstruction parameters. The OSEM algorithm with six iterations and ten subsets and including all the corrections, associated with a Gaussian post-filter of 0.11 cm, gave the most accurate activity quantification with a relative difference of -5.2% from the expected concentration activity in the phantom (Table 2).

A SPECT/CT sensitivity factor of 5.67 ± 0.08 counts/s/MBq was obtained with these optimized reconstruction parameters.

The recovery coefficient obtained for the SPECT/CT acquisition of a 500 ml volume of $^{177}$Lu was 0.95. For smaller spheres (from 0.5 to 16 ml), the recovery coefficients were significantly lower, ranging from 0.43 to 0.78 for the phantom without background activity, and from 0.30 to 0.82 for the phantom with background activity (Fig. 2).

A recovery coefficient of 0.78 was obtained for the phantom with the central 80 ml bottle (L2–L4 red marrow model) flanked by two 200 ml bottles to simulate the kidneys.

## Patient dosimetry

### Patients

Dosimetry was performed after the first two treatment cycles, but for patients 1 and 2. Patient 1 did not have evaluable SPECT/CT data after cycle 1, and patient 2 died before the second cycle of $^{177}$Lu-[DOTA0, Tyr3]-octreotate. Patients received between two and four cycles of $^{177}$Lu-[DOTA0, Tyr3]-octreotate (Table 1). For six patients, SPECT/CT image acquisitions were not performed at all four time points after the first and second cycle (Table 1), mainly due to health reasons, technical issues, or logistic reasons.

### Dosimetry

Time activity curves and residence times (n = 12 patients) were computed for liver, kidneys, spleen, and red marrow (Fig. 3).

## Personalized organ mass versus standard organ mass

The mean personalized organ masses were 1962 g for liver, 470 g for kidneys, 296 g for spleen, and 891 g for red marrow. The maximum $\Delta_{\text{mass i, j}}$ between personalized and standard mass was 130.3% for liver (patient 2) and 307.2% for spleen (patient 11). The $\Delta_{\text{mass i, j}}$ ranged from 0.8 to 93.5% for kidneys, and from -55.1 to 11.9% for red marrow.
<figure><img src="image_5.png" /><figcaption>**Fig. 1** Planar sensitivity variation in function of the distance from the detectors of the gamma camera GE-Discovery NM/CT 670. For each distance (8 cm, 13 cm, 18 cm), the planar sensitivity was estimated from the geometric mean using the Xeleris® software. A circular region of interest that delineates the source distribution was used. The planar sensitivity factor remained constant (6.11 ± 0.01 counts/s/MBq) from 8 to 18 cm</figcaption></figure>

In agreement, the AD/A admin values for the four OARs varied when calculated using the personalized and standard organ mass, particularly for kidneys ($p = 1 \times 10^{-7}$), spleen ($p = 0.0069$), and red marrow ($p = 0.0027$) (Fig. 4). For example, the kidney AD/A admin values for patient 11 after cycle 2 were 0.91 mGy/MBq with the personalized and 1.50 mGy/MBq with the standard organ mass. In patient 12, the red marrow AD/A admin

values after cycle 2 reached a value of 0.09 mGy/MBq with the personalized organ mass and 0.06 mGy/MBq with the standard one. A large difference was observed for the spleen AD/A admin values, especially in patient 11 after cycle 1 (0.70 mGy/MBq versus 3.15 mGy/MBq).

**AD/A admin values after cycle 1 and 2** The AD/A admin values obtained using the personalized organ

<table><caption>**Table 2** Deviation between the activity concentration calculated from SPECT/CT images and the real activity concentration inside the bottle to simulate a kidney</caption><thead><tr><th>OSEM reconstruction corrections</th><th>Number of iterations (i) and subsets (ss)</th><th>Filter</th><th>Calculated activity concentration (MBq/ml)</th><th>Real activity concentration (MBq/ml)</th><th>Deviation (%)</th></tr></thead><tbody><tr><td rowspan="10">AC + RR + SC</td><td>4i*5ss</td><td rowspan="5">None</td><td>0.506</td><td rowspan="5">0.537</td><td>- 5.88</td></tr><tr><td>4i*10ss</td><td>0.495</td><td>- 7.88</td></tr><tr><td>6i*10ss</td><td>0.502</td><td>- 6.55</td></tr><tr><td>8i*10ss</td><td>0.503</td><td>- 6.33</td></tr><tr><td>16i*10ss</td><td>0.502</td><td>- 6.55</td></tr><tr><td>32i*10ss</td><td>0.507</td><td>- 5.66</td></tr><tr><td>4i*5ss</td><td rowspan="5">Gaussian 0.11 cm</td><td>0.507</td><td rowspan="5">0.537</td><td>- 5.66</td></tr><tr><td>4i*10ss</td><td>0.494</td><td>- 88.0</td></tr><tr><td>6i*10ss</td><td>0.509</td><td>- 5.22</td></tr><tr><td>8i*10ss</td><td>0.507</td><td>- 5.66</td></tr><tr><td>16i*10ss</td><td>0.505</td><td>- 66.0</td></tr></tbody></table>

AC attenuation correction, RR recovery resolution, SC scatter correction
<figure><img src="image_6.png" /><figcaption>**Fig. 2** Recovery coefficients in function of the volume using six iterations and ten subsets with a Gaussian post-filter of 0.11 cm, and attenuation, scatter and recovery resolution correction. For the objects ranging from 0.5 to 16 ml, the recovery coefficients ranged from 0.43 to 0.78 when using the phantom without background radioactivity (dash-dot line) and from 0.30 to 0.82 with the phantom with background radioactivity (dotted line)</figcaption></figure>

masses are presented in Table 3. The highest AD/A admin values were 2.26 mGy/MBq in liver (patient 7 after cycle 1), 0.91 mGy/MBq in kidneys (patient 11 after cycle 2), 0.88 mGy/MBq in spleen (patient 3 after cycle 1), and 0.09 mGy/MBq in red marrow (patient 12 after cycle 2).

Comparison of the AD/A admin values in ten patients (patients 3 to 12) after the first and second cycle of ¹⁷⁷Lu-[DOTA0, Tyr3]-octreotate (Fig. 4) highlighted some minor between-dose differences in all four OARs (p = 0.0563 for kidneys, p = 0.3414 for spleen, p = 0.5106 for liver, and p = 0.7220 for red marrow). The $\bar{\Delta}_{AD/A admin}$ was 14.4% (range 0–41.4%) for liver, 12.5% (range 2–41.4%) for kidneys, 14.4% (range 2.2–36.4%) for spleen, and 32.8% (range 0–80.9%) for red marrow. The highest between-dose differences were observed in liver, in patient 7 (2.26 and 1.22 mGy/MBq for cycle 1 and cycle 2; $\bar{\Delta}_{AD/A admin liver, P7}$ of 42.3%) and patient 11 (0.64 and 1.03 mGy/MBq; $\bar{\Delta}_{AD/A admin liver, P7}$ of 33%).

Similarly, when comparing the absorbed dose (in Gy) between cycles, intra-patient differences were observed especially in organs very close to or including tumor cells or metastases (Fig. 5). This was the case for liver in patient 7 (diffused liver metastases) and for red marrow in patients 3 and 11 (tumors in contact with L2–L4 and L2, L3, or L4 vertebral metastases) (data not shown). For organs without metastases or distant from the tumor, no significant between-cycle difference was observed (p = 0.0809, p = 0.4604, p = 0.5083, and p = 0.8720 for kidneys, spleen, liver, and red marrow respectively), except

for kidneys in patient 11 due to an obstruction of the double pigtail stent before the second injection.

**Mean absorbed dose per unit of administered activity to organs for all patients** The mean organ AD/A admin for all patients and the first two cycles were 0.54 ± 0.58 mGy/MBq for liver, 0.43 ± 0.13 mGy/MBq for kidneys, 0.61 ± 0.13 mGy/MBq for spleen, and 0.04 ± 0.02 mGy/MBq for red marrow (Fig. 6).

## Discussion

This study reports our experience with dosimetry in patients with NET after LUTATHERA® treatment, and the absorbed doses to the kidneys, liver, spleen, and red marrow, using SPECT/CT data and the commercial workstation Xeleris®.

Activity quantification from the SPECT images is a crucial step. Therefore, calibration of the gamma camera is a pre-requisite procedure in targeted radionuclide dosimetry [24, 25]. In our study, we reproduced the calibration method described in the MIRD pamphlet n°26 [18] using a large water cylinder containing a well-calibrated source of ¹⁷⁷Lu, a phantom often employed in clinical studies [15, 26, 27]. We used the SPECT sensitivity factor according to the study by Frey et al. [13], and to a recent multi-center evaluation on phantoms [28]. The obtained recovery coefficient results and reconstruction parameters were in agreement with those of recently published studies [29, 30]. For patient dosimetry, SPECT/CT image acquisition at different time points was
<figure><img src="image_7.png" /><figcaption>**Fig. 3** Dosimetry toolkit® workflow. **a** Example of segmentation of the organs at risk performed manually using the Dosimetry Toolkit® tools using SPECT or CT data. **b** Segmentation propagation, time activity curves and residence times for all the regions of interest at 4 h, 24 h, 72 h and 192 h after drug administration</figcaption></figure>

chosen to avoid the overlap with local high activity areas near the kidneys that could give an overestimation of the kidney activity, as described by Garkavij et al. [31].

In addition, to better evaluate the accuracy of the AD/ A admin values, especially for red marrow, we used an- other more suitable phantom to model the in vivo con- ditions in terms of organ volume and activity concentration. With this phantom, the recovery coeffi- cient obtained for the L2–L4 model of 80 ml was in the same order of magnitude as the value obtained for a 16 ml sphere without background activity, possibly due to the high concentration gradient between the L2–L4 model and kidney model.

To evaluate the quantification accuracy, other calibra-
tion methods with more complex phantoms have been de-
scribed. For example, Gnesin et al. [32] used an
anthropomorphic phantom (Kyoto liver/kidney phantom)
including kidneys and liver with spheres of 20, 30, and
40 mm in diameter to mimic lesions. Recently, Tran-Gia

et al. [33] presented a 3D-printed two-compartment kid-
ney phantom. A similar phantom for red marrow quantifi-
cation and dosimetry remains a real challenge. According
to Lassman et al. [16], the joint European project on MRT
dosimetry is working on a standardization of SPECT/CT
calibration procedures that will be available by the end of
2019. Given the diversity of dosimetry data published
using different methodologies, Gear et al. [34] recently
proposed a practical guide to express the accuracy of the
dosimetry results by determining the uncertainties at each
step of the dosimetry chain.

After this preliminary calibration step using phantoms
that was necessary to obtain accurate quantification of
the activity, we implemented the method in the clinic.
We observed considerable differences of AD/A admin
values for all the OARs, using the personalized and the
standard organ masses. The personalized organ masses
usually led to lower kidneys and spleen AD/A admin
values, as previously reported by Kupitz et al. [20], and
<figure><img src="image_8.png" /></figure>

higher red marrow AD/A admin values in most patients. Thereby a correction for the individual organ masses was applied.

Based on the report [5] by Eberlein et al. that summarizes data on the red marrow, kidneys, and tumor absorbed doses from various studies, the range of the kidney AD/A admin values in our patients was in concordance with the values reported by Sandström et al.

[27], Gupta et al. [35], and Heikkonen et al. [26]. The mean liver AD/A admin value was slightly higher than what previously published, possibly due to the presence of liver metastases in patient 7. Indeed, the mean and the range of liver AD/A admin values were in accordance with the literature after exclusion of this patient [35]. The mean red marrow AD/A admin was similar to published data. However, in our study, the dosimetry

<table><caption>**Table 3** Absorbed dose per administered activity (AD/A admin; in mGy/MBq) to the liver, kidneys, spleen, and red marrow after cycle 1 and cycle 2. The AD/A admin presented in the table were calculated with personalized organ masses</caption><thead><tr><th rowspan="3">Patients</th><th colspan="8">Absorbed dose per administered activity (mGy/MBq)</th></tr><tr><th colspan="2">Liver</th><th colspan="2">Kidneys</th><th colspan="2">Spleen</th><th colspan="2">Red marrow</th></tr><tr><th>Cycle 1</th><th>Cycle 2</th><th>Cycle 1</th><th>Cycle 2</th><th>Cycle 1</th><th>Cycle 2</th><th>Cycle 1</th><th>Cycle 2</th></tr></thead><tbody><tr><td>1</td><td>NA</td><td>0.29</td><td>NA</td><td>0.42</td><td>NA</td><td>0.89</td><td>NA</td><td>0.07</td></tr><tr><td>2</td><td>1.92</td><td>†</td><td>0.28</td><td>†</td><td>0.59</td><td>†</td><td>0.07</td><td>†</td></tr><tr><td>3</td><td>0.17</td><td>0.17</td><td>0.47</td><td>0.49</td><td>0.88</td><td>0.59</td><td>NA</td><td>NA</td></tr><tr><td>4</td><td>0.46</td><td>0.32</td><td>0.39</td><td>0.46</td><td>0.42</td><td>0.58</td><td>0.02</td><td>0.03</td></tr><tr><td>5</td><td>0.28</td><td>0.23</td><td>0.34</td><td>0.35</td><td>0.49</td><td>0.48</td><td>0.04</td><td>0.01</td></tr><tr><td>6</td><td>0.15</td><td>0.19</td><td>0.27</td><td>0.36</td><td>0.45</td><td>0.76</td><td>0.03</td><td>0.02</td></tr><tr><td>7</td><td>2.26</td><td>1.22</td><td>0.40</td><td>0.49</td><td>0.55</td><td>0.70</td><td>0.03</td><td>0.02</td></tr><tr><td>8</td><td>0.13</td><td>0.12</td><td>0.44</td><td>0.54</td><td>0.52</td><td>0.55</td><td>0.04</td><td>0.04</td></tr><tr><td>9</td><td>0.33</td><td>0.31</td><td>0.34</td><td>0.32</td><td>0.58</td><td>0.56</td><td>0.07</td><td>0.05</td></tr><tr><td>10</td><td>0.55</td><td>0.59</td><td>0.45</td><td>0.50</td><td>0.57</td><td>0.75</td><td>0.03</td><td>0.04</td></tr><tr><td>11</td><td>0.64</td><td>1.03</td><td>0.50</td><td>0.91</td><td>0.70</td><td>0.76</td><td>NA</td><td>NA</td></tr><tr><td>12</td><td>0.22</td><td>0.22</td><td>0.35</td><td>0.37</td><td>0.56</td><td>0.51</td><td>0.06</td><td>0.09</td></tr></tbody></table>

NA not available, † patient dead
<figure><img src="image_9.png" /><figcaption>**Fig. 5** Absorbed doses in Gy in liver, kidneys, spleen, and red marrow after cycle 1 and after cycle 2. As our dosimetry methodology was based on L2-L3-L4 red marrow delineation, it could not be used for patients 3 and 11 who had tumors in contact with L2, L3, L4 and L2, L3 or L4 vertebral metastases</figcaption></figure>

methodology was not based on blood sampling, but on L2–L4 red marrow delineation, and could not be applied to patient 3 and 11 who had tumors in contact with L2–L4 vertebrae and L2, L3, or L4 vertebral metastases. With our methodology, a high residence time for this volume of interest was extrapolated to the total red marrow, but in these two patients it was not representative. When excluding the data for these two patients from the analysis, the red marrow AD/A admin value was concomitant with previously published data. Ferrer et al. [36]

showed that although blood-based approaches are easier to implement, an imaging-based method better predicts the hematological toxicity. They reported that the approach based on L2–L4 vertebra imaging seemed to provide the best relation between absorbed dose to red marrow and platelet toxicity in radioimmunotherapy with ⁹⁰Y-epratuzumab tetraxetan. Blakkisrud et al. evaluated the correlation between absorbed dose and hematologic toxicity in patients treated with ¹⁷⁷Lu-lilotomab satetraxetan [37], and concluded that hematologic
<figure><img src="image_10.png" /><figcaption>**Fig. 6** Mean absorbed dose to organs per administered activity in mGy/MBq for all patients. The error bars show the range of the AD/A admin values for cycle 1 and 2</figcaption></figure>

toxicity can be predicted by calculating the red marrow absorbed dose from SPECT/CT images. Recently, Beykan et al [38] estimated that during PRRT with ¹⁷⁷Lu-DOTATATE, the blood-based absorbed dose to bone marrow was three times lower than the imaging-based absorbed dose, advocating the use of imaging-based dosimetry.

In our study, intra-patient differences of absorbed dose between cycle 1 and 2 were observed, especially for organs very close to, or including tumor cells or metastases. This could be explained by a lower uptake in the treated lesions between cycles. This suggests that the activity kinetic could be similar in organs containing or in contact with tumors and in the NET, and could explain the difference in absorbed dose between cycles.

For the kidneys, the reproducibility between absorbed dose after cycle 1 and 2 appeared at the limit of the significance. Indeed, intra-patient variability of the absorbed dose between cycles was previously reported by Sundlöv et al [11], possibly due to intra-patient pharmacokinetic variations between cycles. Particularly, variations of blood pressure and body hydration, and thus renal perfusion, could affect the residence time and the absorbed dose in kidneys.

Previous studies [35, 39] have shown that also in healthy tissues, the time-activity curve shape can differ considerably among patients and in the same patient between treatment cycles. The impact of kinetics and image timing on the absorbed dose to organs is crucial [40]. Consequently, our dosimetric results are clearly restricted to the use of the Dosimetry Toolkit® software

because the time activity curves for each segmented organ were fitted automatically by a mono-exponential function to calculate the residence times.

Although other analytic fit functions, such as bi- or tri-exponential functions, are known to better fit the time-activity curves [40, 41], these options were not implemented in the software. Moreover, the Dosimetry Toolkit® software has a limited set of segmentation tools, is not user-friendly, and not optimal. Another major drawback is that this software cannot save the images and volumes of interest in Dicom format from the beginning to the end of the process. For example, it is not possible to change organ delineation after the process. In addition, the Dosimetry Toolkit® software is associated with a specific gamma camera by a license system. The software only accepts images obtained with a gamma camera with a known license number. Consequently, imaging data obtained with other gamma cameras cannot be used. However, new solutions are commercially available with user-friendly tools for segmentation, including those commonly used for EBRT (such as Boolean operators, easy manual, semi-automatic, automatic segmentation options...) and for nuclear medicine. These software programs also propose a wide choice of interpolation methods to better fit the obtained time-activity data.

From a practical point of view, as reported in the recent review by Huizing et al [42], dosimetry in PRRT is not routinely performed due to implementation difficulties, the time-consuming, non-standard dosimetry approach, and the lack of evidence in literature. Our study
shows that routine patient dosimetry is achievable in a nuclear medicine department, if there is a strong medical willingness and awareness on how to perform patient dosimetry. As discussed by Flux et al. [43], a solid collaboration is needed between nuclear medicine physicians, medical physicists, and nuclear medicine technologists to implement dosimetry and provide the best patient care. Dosimetry of organs at risk and tumors brings real benefits for patients; however, health professionals must be aware that it is time-consuming using a gamma camera, and must anticipate the practical constraints before its routine implementation. In our clinical dosimetry protocol, each SPECT/CT acquisition by the gamma camera required about 25 min (i.e., an additional acquisition time of 1.25 h for the dosimetry of one patient after one treatment cycle). The registration, segmentation, and dosimetry steps also are time-consuming (around 1.5 h for one dosimetry). Despite the longer exam time, all imaging acquisitions were well received by our patients.

According to Hänscheid et al. [44], besides the additional burden for patients, this may also increase the complexity of data collection and analysis. Therefore, they proposed a simplified PRRT dosimetry method using a theoretical approach based on one single measurement of the activity retention at 96 h post-drug administration. Recently Sundlöv et al [45] evaluated in 22 patients (119 treatment cycles) a treatment planning schema that gives reliable dosimetric results with simplified acquisitions. They confirmed previous conclusions with the use of one SPECT/CT at 96 h after each treatment cycle. However, these results have to be confirmed with a large patient cohort.

## Conclusions

This personalized dosimetric study in patients with NET who underwent PPRT with ¹⁷⁷Lu-DOTATATE highlights the important inter-patient/intra-patient variations and consequently the need of individualized dosimetry. Additional studies are needed to estimate the absorbed doses to tumors and to perform dose-response correlations.

It also shows that implementing patient dosimetry in the clinical practice during PPRT with ¹⁷⁷Lu-DOTATA TE is feasible in a nuclear medicine department.

## Abbreviations

AD/A admin: Absorbed dose per administered activity; EBRT: External beam radiotherapy; MEGP: Medium-energy general purpose; MIRD: Medical internal radiation dose; NET: Neuroendocrine tumor; OAR: Organ at risk; OSEM: Ordered subset expectation maximization; PRRT: Peptide receptor radionuclide therapy; SPECT/CT: Single-photon emission computed tomography/computed tomography; WB: Whole body

## Acknowledgements

The authors wish to thank Dr. Hélène de Forges and Dr. Elisabetta Andermarcher for their help in writing and editing the manuscript.

## Funding

This study was supported by the Montpellier Cancer Institute.

## Availability of data and materials

Please contact the author for data requests.

## Authors' contributions

LS, EM-R, DT, and SC carried out the phantom measurements, image reconstruction, and data analysis. PE helped with phantom measurements. JPP, MB, and ED contributed to the study design. POK and ED were responsible for the patients and their treatment. LS managed the dosimetry software implementation and performed the patient dosimetry analyses with DT and SC. LS wrote the manuscript. All authors read and approved the final manuscript.

## Ethics approval and consent to participate

All procedures performed in this study involving human participants were in accordance with the ethical standards of the institutional research committee.

## Consent for publication

Consent has been obtained from participants to publish this work.

## Competing interests

The authors declare that they have no competing interests.

## Publisher's Note

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

## Author details

¹Nuclear Medicine Department, Montpellier Cancer Institute (ICM), University of Montpellier, 208 Avenue des Apothicaires, 34298 Montpellier Cedex5, France. ²Centre de Recherche en Cancérologie de Toulouse, Toulouse, France. ³INSERM, UMR 1037, Toulouse III Paul Sabatier University, Toulouse, France. ⁴University of Costa Rica, Physics School, CICANUM, San Jose, Costa Rica. ⁵Institut de Recherche en Cancérologie de Montpellier (IRCM), INSERM U1194, Université de Montpellier, Institut Régional du Cancer de Montpellier (ICM), Montpellier, France.

Received: 27 August 2018 Accepted: 20 November 2018

Published online: 29 November 2018

## References

1. Bodei L, Cremonesi M, Ferrari M, Pacifici M, Grana CM, Bartolomei M, et al. Long-term evaluation of renal toxicity after peptide receptor radionuclide therapy with 90Y-DOTATOC and 177Lu-DOTATATE: the role of associated risk factors. Eur J Nucl Med Mol Imaging. 2008;35(10):1928.

2. Strigari L, Konijnenberg M, Chiesa C, Bardies M, Du Y, Gleisner KS, et al. The evidence base for the use of internal dosimetry in the clinical practice of molecular radiotherapy. Eur J Nucl Med Mol Imaging. 2014;41(10):1976–88.

3. Ilan E, Sandstrom M, Wassberg C, Sundin A, Garske-Roman U, Eriksson B, et al. Dose response of pancreatic neuroendocrine tumors treated with peptide receptor radionuclide therapy using 177Lu-DOTATATE. J Nucl Med. 2015;56(2):177–82.

4. Bodei L, Kidd M, Paganelli G, Grana CM, Drozdov I, Cremonesi M, et al. Long-term tolerability of PRRT in 807 patients with neuroendocrine tumours: the value and limitations of clinical factors. Eur J Nucl Med Mol Imaging. 2015;42(1):5–19.

5. Eberlein U, Cremonesi M, Lassmann M. Individualized dosimetry for theranostics: necessary, nice to have, or counterproductive? J Nucl Med. 2017;58(Supplement 2):97S–103S.

6. Strosberg J, El-Haddad G, Wolin E, Hendifar A, Yao J, Chasen B, et al. Phase 3 trial of ¹⁷⁷Lu-Dotatate for midgut neuroendocrine tumors. N Engl J Med. 2017;376(2):125–35.

7. Cremonesi M, Ferrari M, Bodei L, Tosi G, Paganelli G. Dosimetry in peptide radionuclide receptor therapy: a review. J Nucl Med. 2006;47(9):1467–75.

8. Kwekkeboom DJ, Kam BL, van Essen M, Teunissen JJM, van Eijck CHJ, Valkema R, et al. Somatostatin receptor-based imaging and therapy of gastroenteropancreatic neuroendocrine tumors. Endocr Relat Cancer. 2010; 17(1):R53–73.
9. Wessels BW, Konijnenberg MW, Dale RG, Breitz HB, Cremonesi M, Meredith RF, et al. MIRD pamphlet no. 20: the effect of model assumptions on kidney dosimetry and response—implications for radionuclide therapy. J Nucl Med. 2008;49(11):1884-99.

10. Forrer F, Krenning EP, Kooij PP, Bernard BF, Konijnenberg M, Bakker WH, et al. Bone marrow dosimetry in peptide receptor radionuclide therapy with [177Lu-DOTA0, Tyr3] octreotate. Eur J Nucl Med Mol Imaging. 2009;36(7):1138-46.

11. Sundlöv A, Sjögreen-Gleisner K, Svensson J, Ljungberg M, Olsson T, Bernhardt P, et al. Individualised 177Lu-DOTATATE treatment of neuroendocrine tumours based on kidney dosimetry. Eur J Nucl Med Mol Imaging. 2017;44(9):1480-9.

12. Chiesa C, Sjogreen Gleisner K, Flux G, Gear J, Walrand S, Bacher K, et al. The conflict between treatment optimization and registration of radiopharmaceuticals with fixed activity posology in oncological nuclear medicine therapy. Eur J Nucl Med Mol Imaging. 2017;44(11):1783-6.

13. Frey EC, Humm JL, Ljungberg M. Accuracy and precision of radioactivity quantification in nuclear medicine images. Semin Nucl Med. 2012;42(3):208-18.

14. Marin G, Vanderlinden B, Karfis I, Guiot T, Wimana Z, Flamen P, et al. Accuracy and precision assessment for activity quantification in individualized dosimetry of 177Lu-DOTATATE therapy. EJNMMI Phys. 2017; 4(1):7 Available from: http://ejnmmiphys.springeropen.com/articles/10.1186/s40658-017-174-7.

15. Hippeläinen E, Tenhunen M, Mäenpää H, Sohlberg A. Quantitative accuracy of 177Lu SPECT reconstruction using different compensation methods: phantom and patient studies. EJNMMI Res. 2016;6(1):16 Available from: http://www.ejnmmires.com/content/6/1/16.

16. Lassmann M, Eberlein U. The relevance of dosimetry in precision medicine. J Nucl Med. 2018;59:1494-1499.

17. Dewaraja YK, Frey EC, Sgouros G, Brill AB, Roberson P, Zanzonico PB, et al. MIRD pamphlet no. 23: quantitative SPECT for patient-specific 3-dimensional dosimetry in internal radionuclide therapy. J Nucl Med. 2012;53(8):1310-25.

18. Ljungberg M, Celler A, Konijnenberg MW, Eckerman KF, Dewaraja YK, Sjogreen-Gleisner K. MIRD pamphlet no. 26: joint EANM/MIRD guidelines for quantitative 177Lu SPECT applied for dosimetry of radiopharmaceutical therapy. J Nucl Med. 2016;57(1):151-62.

19. GE Healthcare. Organ dose estimates for radio-isotope therapy treatment planning purposes. Dosimetry toolkit package. White Paper. 2011.

20. Kupitz D, Wetz C, Wissel H, Wedel F, Apostolova I, Wallbaum T, et al. Software-assisted dosimetry in peptide receptor radionuclide therapy with 177Lutetium-DOTATATE for various imaging scenarios. PLOS ONE. 2017; 12(11):e0187570.

21. Stabin MG, Sparks RB, Crowe E. OLINDA/EXM: the second-generation personal computer software for internal dose assessment in nuclear medicine. J Nucl Med. 2005;46(6):1023-7.

22. Shen S, Meredith RF, Duan J, Macey DJ, Khazaeli MB, Robert F, et al. Improved prediction of myelotoxicity using a patient-specific imaging dose estimate for non-marrow-targeting 90Y-antibody therapy. J Nucl Med. 2002; 43(9):1245-53.

23. Wieira JW, Lima FRA and Kramer, R. Calculation of red bone marrow equivalent dose to the MAX (Male Adult voXel) model for external exposures to photons. BONE. 2002;5500(1.92):2864.6.

24. Zhao W, Esquinas PL, Hou X, Uribe CF, Gonzalez M, Beauregard J-M, et al. Determination of gamma camera calibration factors for quantitation of therapeutic radioisotopes. EJNMMI Phys. 2018;5(1):8 Available from: https://ejnmmiphys.springeropen.com/articles/10.1186/s40658-018-0208-9.

25. Uribe CF, Esquinas PL, Tanguay J, Gonzalez M, Gaudin E, Beauregard J-M, et al. Accuracy of 177Lu activity quantification in SPECT imaging: a phantom study. EJNMMI Phys. 2017;4(1):2 Available from: http://ejnmmiphys. springeropen.com/articles/10.1186/s40658-016-0170-3.

26. Heikkonen J, Mäenpää H, Hippeläinen E, Reijonen V, Tenhunen M. Effect of calculation method on kidney dosimetry in 177 Lu-octreotate treatment. Acta Oncol. 2016;55(9-10):1069-76.

27. Sandstrom M, Garske-Roman U, Granberg D, Johansson S, Widstrom C, Eriksson B, et al. Individualized dosimetry of kidney and bone marrow in patients undergoing 177Lu-DOTA-Octreotate treatment. J Nucl Med. 2013; 54(1):33-41.

28. Zimmerman BE, Grošev D, Buvat I, Coca Pérez MA, Frey EC, Green A, et al. Multi-centre evaluation of accuracy and reproducibility of planar and SPECT image quantification: an IAEA phantom study. Z Für Med Phys. 2017;27(2):98-112.

29. Hippeläinen ET, Tenhunen MJ, Mäenpää HO, Heikkonen JJ, Sohlberg AO. Dosimetry software Hermes internal radiation dosimetry: from quantitative image reconstruction to voxel-level absorbed dose distribution. Nucl Med Commun. 2017;38(5):357-65.

30. Sanders JC, Kuwert T, Hornegger J, Ritt P. Quantitative SPECT/CT imaging of 177Lu with in vivo validation in patients undergoing peptide receptor radionuclide therapy. Mol Imaging Biol. 2015;17(4):585-93.

31. Garkavij M, Nickel M, Sjögreen-Gleisner K, Ljungberg M, Ohlsson T, Wingårdh K, et al. 177Lu-[DOTA0, Tyr3] octreotate therapy in patients with disseminated neuroendocrine tumors: Analysis of dosimetry with impact on future therapeutic strategy. Cancer. 2010;116(S4):1084-92.

32. Gnesin S. Quantitative Lu-177 SPECT/CT validation to assist theragnostic procedures. European journal of nuclear medicine and molecular imaging volume 44, supplement 2. Austria: EANM congress Vienna; 2017. https://doi.org/10.1007/s00259-017-3822-1.

33. Tran-Gia J, Lassmann M. Optimizing image quantification for 177 Lu SPECT/CT based on a 3D printed 2-compartment kidney phantom. J Nucl Med. 2018;59(4):616-24.

34. Gear JI, Cox MG, Gustafsson J, Gleisner KS, Murray I, Glatting G, et al. EANM practical guidance on uncertainty analysis for molecular radiotherapy absorbed dose calculations. Eur J Nucl Med Mol Imaging. 2018;45(13):2456-74 Available from: http://link.springer.com/10.1007/s00259-018-4136-7.

35. Gupta SK, Singla S, Thakral P, Bal C. Dosimetric analyses of kidneys, liver, spleen, pituitary gland, and neuroendocrine tumors of patients treated with 177Lu-DOTATATE. Clin Nucl Med. 2013;38(3):7.

36. Ferrer L, Kraeber-Bodéré F, Bodet-Millin C, Rousseau C, Gouill SL, Wegener WA, et al. Three methods assessing red marrow dosimetry in lymphoma patients treated with radioimmunotherapy. Cancer. 2010;116(S4):1093-100.

37. Blakkisrud J, Løndalen A, Dahle J, Turner S, Holte H, Kolstad A, et al. Red marrow-absorbed dose for non-Hodgkin lymphoma patients treated with 177 Lu-Lilotomab Satetraxetan, a novel anti-CD37 antibody-radionuclide conjugate. J Nucl Med. 2017;58(1):55-61.

38. Beykan, S et al. Patient-specific dosimetry of 177Lu-dotatate peptide receptor radionuclide therapy with high activities. Eur J Nucl Med Mol Imaging. 2018;45(Suppl 1):S33. OP-081 EANM congress Düsseldorf Germany.

39. Larsson M, Bernhardt P, Svensson JB, Wängberg B, Ahlman H, Forssell-Aronsson E. Estimation of absorbed dose to the kidneys in patients after treatment with 177 Lu-octreotate: comparison between methods based on planar scintigraphy. EJNMMI Res. 2012;2(1):49.

40. Guerriero F, Ferrari ME, Botta F, Fioroni F, Grassi E, Versari A, et al. Kidney dosimetry in 177 Lu and 90 Y peptide receptor radionuclide therapy: influence of image timing, time-activity integration method, and risk factors. Biomed Res Int. 2013;2013:1-12.

41. Kletting P, Schimmel S, Hänsccheid H, Luster M, Fernández M, Nosske D, et al. The NUKDOS software for treatment planning in molecular radiotherapy. Z Für Med Phys. 2015;25(3):264-74.

42. Huizing DMV, de Wit-van der Veen BJ, Verheij M, MPM S. Dosimetry methods and clinical applications in peptide receptor radionuclide therapy for neuroendocrine tumours: a literature review. EJNMMI Res. 2018;8(1):89 Available from: https://ejnmmires.springeropen.com/articles/10.1186/s13550-018-0443-z.

43. Flux GD, Sjogreen Gleisner K, Chiesa C, Lassmann M, Chouin N, Gear J, et al. From fixed activities to personalized treatments in radionuclide therapy: lost in translation? Eur J Nucl Med Mol Imaging. 2018;45(1):152-4.

44. Hänsccheid H, Lapa C, Buck AK, Lassmann M, Werner RA. Dose mapping after endoradiotherapy with 177 Lu-DOTATE-/TOC by one single measurement after four days. J Nucl Med. 2018;59(1):75-81.

45. Sundlöv A, Gustafsson J, Brolin G, Mortensen N, Hermann R, Bernhardt P, et al. Feasibility of simplifying renal dosimetry in 177Lu peptide receptor radionuclide therapy. EJNMMI Phys. 2018;5(1):12 Available from: https://ejnmmiphys.springeropen.com/articles/10.1186/s40658-018-0210-2.