Abstract
Dosimetry of organs and tumors helps to assess risks and benefit of treatment with 177Lu-DOTATATE/DOTATOC. However, it is often not performed in clinical routine because of additional efforts, the complexity of data collection and analysis, and the additional burden for the patients. Aiming at a simplification of dosimetry, we analyzed the accuracy of a theoretically substantiated approximation, which allows the calculation of absorbed doses from a single measurement of the abdominal activity distribution. Methods: Activity kinetics were retrospectively assessed from planar images in 29 patients with neuroendocrine tumors (NETs; n = 21) or meningioma (n = 8) after the administration of 177Lu-DOTATATE (n = 22) or 177Lu-DOTATOC (n = 7). Mono- or biexponential functions were fitted to measured data in 54 kidneys, 25 livers, 27 spleens, and 30 NET lesions. It was evaluated for each fit function how well the integral over time was represented by an approximation calculated as the product of the time tl of a single measurement, the expected reading at time tl, and the factor 2/ln(2). Tissue-specific deviations of the approximation from the time integral were calculated for time points tl of 24, 48, 72, 96, 120, and 144 h. Results: The correlation between time integral and approximation improved with increasing time tl. Pearson r exceeded 0.95 for a tl of 96 h or more in all tissues. The lowest maximum errors were observed at a tl of 96 h, with deviations of the approximation from the time integral of median +5% (range, −9% to +17%) for kidneys, +6% (range, −7% to +12%) for livers, +8% (range, +2% to +20%) for spleens, and +6% (range, −11% to +16%) for NET lesions. Accuracy was reduced for measurements after 72 or 120 h. For measurements after 24, 48, and 144 h, the approximation led to large deviations for some of the patients, in particular unacceptable underestimates of the absorbed dose to the kidneys. Conclusion: A single quantitative measurement of the abdominal activity concentration by SPECT/CT 4 d after the administration of 177Lu-DOTATATE/DOTATOC provides a 3-dimensional dose map and can be used to estimate the doses actually absorbed in the treatment cycle with minor additional resources and effort.
Peptide receptor radionuclide therapy (PRRT) using 177Lu-DOTATATE/DOTATOC enables delivery of β-radiation to tumors with high expression of somatostatin receptors. Before the recent publication of the first data from a randomized, controlled, prospective study (NETTER-1) (1), nonrandomized observational studies had already demonstrated the superior overall benefit of PRRT (2–4). NETTER-1 has now confirmed these data, showing that treatment with 177Lu-DOTATATE resulted in markedly longer progression-free survival and a significantly higher response rate than high-dose octreotide alone among patients with advanced midgut neuroendocrine tumors (1). Rising numbers of PRRT are already foreseeable.
The currently routinely used medium-energy γ-/β-emitting nuclide 177Lu opens the window for posttherapeutic scintigraphy as well as pretreatment dosimetry (5). The assessment of absorbed doses to tumors and organs can be used to estimate the therapeutic effect to the target tissue and to reduce side effects to critical organs, such as bone marrow or kidneys (5–8).
Although there is evidence that individually tailored treatment accurately planned by dosimetry is advantageous for PRRT optimization (8–10), standard activities are usually administered in a series of treatment cycles. The 1-size-fits-all concept of administering a fixed activity of 7.4 GBq, as used in the NETTER-1 trial, is covered by the recent practical guidance of the European Association of Nuclear Medicine (5). The cumulated absorbed dose to the kidneys, which often are the organs at risk in endoradiotherapy of NET, should not exceed predefined limits in a series of PRRT cycles. Dose verification after each PRRT cycle is therefore desirable to ensure safety and to assess the tumor dose indicating treatment efficacy (11). However, repeated measurements of the activities retained in assessed tissues over several days and an elaborate data evaluation are necessary for an exact dosimetry (12,13). The procedures require several patient visits and result in additional burdens. Noncompliance with procedural instructions, the complexity of dosimetry data interpretation, and flaws of routinely key data recording as part of clinical processes are main factors that are related to low performance rates of PRRT dosimetry studies (14). Attempts to find methods to reduce the efforts for data collection and dosimetry are desirable and have already been performed (15,16).
Aiming at a simplified dosimetry, which is well founded by a theoretic approach, we analyzed the accuracy of an absorbed dose estimate based on a single measurement of the activity retention.
MATERIALS AND METHODS
Study Design
This study comprises a retrospective analysis of the activity kinetics in organs and tumors after 177Lu-DOTATATE/DOTATOC from routinely acquired patient data. Included patients gave written informed consent for the recording and anonymized analysis of their data. The local ethics committee waived the need for further approval.
Patient Measurements
Activity kinetics were analyzed in the kidneys (excluding renal pelvis), livers, spleens, and in up to 2 NET lesions in patients who underwent 177Lu-DOTATATE/DOTATOC. The first assessment for each patient with at least 4 planar whole-body scans, including measurements after 1–4 h, 1 d, 2 d, and 4 d or more after the activity administration, was analyzed. Data from 29 patients (17 women, 12 men) with histologically proven NET (n = 21) or meningioma (n = 8) who underwent PRRT at our institution between March 2010 and October 2016 were included. Patients underwent 177Lu-DOTATATE (n = 22) or 177Lu-DOTATOC (n = 7) for pretherapeutic dosimetry (n = 7; 164−225 MBq of 177Lu) or treatment (n = 22; 5.3–8.1 GBq of 177Lu). Only organ but not tumor kinetics were analyzed in the 8 patients with meningioma because the meningiomas showed significantly different activity kinetics as compared with NET lesions, in particular a shorter effective half-life.
Measurements were performed with dual head γ-cameras (DUET or Symbia; Siemens) equipped with medium-energy collimators. All whole-body scans in a series were acquired using the same camera and identical camera settings (matrix, 1,024 × 256; 20% window at 208 keV; identical measuring distance). Regions of interest were drawn including the organ or NET lesion, or part of the tissue in the case of overlapping accumulating tissues, under consideration and over an area with representative background. Identical regions were copied to each scan in the series, and background-corrected counting rates were extracted. Organs completely masked by uptake in other tissues were considered unmeasurable and were excluded from further analyses. The images were not corrected for attenuation and scatter because attenuation maps and triple-energy data were not available for most of the routinely acquired whole-body scans. Activities in the regions of interest were not quantified; only the temporal development of the counting rates was determined assuming identical attenuation, scatter, and sensitivity in all scans in a series.
Theory and Data Evaluation
A biexponential decay function was fitted to the counting rate kinetics of each measurable organ or tumor by ordinary least-squares regression, and the Pearson correlation coefficient r was calculated. Some tissues showed a monoexponential decay function, which was adopted if the correspondent fit resulted in an r greater than 0.98. The fit functions of 54 kidneys, 25 livers, 27 spleens, and 30 NET lesions were included in further analysis. Each fit function was normalized to 100% at its maximum to deduce the normalized uptake function 
, where rs denotes the tissue under consideration and t the time after the administration. Figure 1 shows, as an example, the functions 
of the kidneys in a patient with meningioma.
Normalized uptake in left (red; 
) and right (green; 
) kidneys of patient with meningioma after administration of 7 GBq of 177Lu-DOTATATE.
The total absorbed dose D(rS) to rS by self-irradiation is (17):
Eq. 1where 
is the activity time function and 
the dose rate per unit activity in rS. 
can be deduced from 
when the absolute activity 
at any given time tl is known from a quantitative assessment:
Eq. 2Equation 1 becomes:
Eq. 3with 
representing the time integral of 
:
Eq. 4The time integral of a function 
decaying monoexponentially with an effective half-life Teff is 
. 
is related to the uptake 
after time tl by:
Eq. 5Within a certain time interval after the activity administration, which is determined by the value of Teff, the factor 
in Equation 5 can be approximated by the factor 
(18,19) with less than 10% deviation for 0.75 < 
< 2.5 (Fig. 2). The time integral 
of a monoexponential decay function can therefore be estimated from only a single measurement at time tl by the approximation 
:
Eq. 6In PRRT, however, the fit functions 
are often biexponential (12,16,20). To investigate whether the approximation in Equation 6 can nonetheless be used to estimate absorbed doses in PRRT with an acceptable uncertainty from a single measurement, the expected 
at tl = 24, 48, 72, 96, 120, and 144 h were calculated for each uptake function and used to determine the error introduced by the approximation.
Deviation of ratio of 
and the factor 
from unity in time interval 
.
When Equations 3 and 6 are combined, the absorbed dose is approximated with the present approach by:
Eq. 7
RESULTS
Normalized Uptake Functions u(rS,t)
The correlation coefficient r exceeded 0.98 in monoexponential fits in 25 of 54 kidneys, 12 of 25 livers, 3 of 27 spleens, and 7 of 22 tumors. Biexponential functions 
were assumed for the remaining tissues. The biexponential functions generally showed a dominant component, with longer effective half-life and a short-lived component leading to some deviation from a monoexponential decay shortly after the administration. Compared with the time integral of the main component, the time integral of the short-lived component was a median of −5% (range, −17% to −1%) in kidneys, −2% (range, −5% to +7%) in livers, −3% (range, −14% to −1%) in spleens, and −2% (range, −10% to 0%) in tumors. The tissue-specific effective half-life of the main component in all evaluated tissues was a median of 51 h (range, 40–68 h with 1 outlier at 106 h) in kidneys, 67 h (range, 55–117 h) in livers, 68 h (range, 52–99 h) in spleens, and 77 h (range, 56–130 h) in tumors.
Error of Single Sample Estimates
Table 1 shows the observed empiric correlation between the approximation 
and the time integral 
. The Pearson correlation coefficient improved with increasing time tl. A strong linear relationship between 
and 
was observed only at late times tl, indicating that estimates based on early measurements are less reliable.
Correlation Between Approximation 
Deduced from Single Measurement After 24, 48, 72, 96, 120, and 144 Hours and Actual Time Integral 
Figure 3 shows the error introduced by the approximation in Equation 6 in dependence on the actual time integral 
for all evaluated tissues. The error is defined as the percentage deviation of the ratio 
/
from unity. Data points that originate from monoexponential decay functions are located on the gray lines (Fig. 3). Deviations of data points from the gray lines are due to the influence of a second short lived component. Quantiles of the distribution of deviations are listed in Table 2.
Percentage deviation of approximation 
from 
vs. actual time integral 
for single measurements after tl = 72, 96, and 120 h (top; scale, −30% to 25%) and tl = 24, 48, and 144 h (bottom; scale, −70% to 25%). Deviations are positive for overestimations and negative for underestimations of actually absorbed doses. Each point represents kidney (black), liver (green), spleen (yellow), or NET lesion (red). Gray line shows deviation expected for monoexponential decay functions.
Median, Range, and 0.1 and 0.9 Quantiles of Percentage Deviation of Approximation 
Deduced from Single Measurement after 24, 48, 72, 96, 120, and 144 Hours from Actual Time Integral 
The kidneys, which most often represent the dose-limiting organs and typically showed the shortest effective half-life, Teff, were best represented by approximations based on the measurements after 72 and 96 h. For the 96-h values, deviation was a median of +5%, with underestimates of more than 5% (−9% both) in only 2 of 54 kidneys (1 patient) and overestimates of more than 10% (up to 17%) in 6 kidneys (4 patients). The measurement at 120 h was superior for livers, spleens, and tumors, which typically showed longer effective half-lives.
The approximation in Equation 6 is inadequate for measurements after 24 h, which generally led to unacceptable underestimates of absorbed doses in all tissues. Large underestimates were also observed for measurements after 48 h in tissues with a long effective half-life and for 144 h in those with short half-life.
The approximation based on the 96-h measurement provided good agreement for all tissues. The median deviation for all organs and NET lesions was +6% (range, −11% to +20%), and 90% of deviations ranged between −4% and +15%.
DISCUSSION
The method presented in this work is a potential option to deduce an estimate of the absorbed doses in abdominal organs and NET lesions from a single late-activity measurement. It is not based on an empirically observed correlation, but rather provides an approximation that is well substantiated by theory when an adequate time of measurement is selected and the decay function shows only minor deviations from a monoexponential decay.
Time of Measurement
The error introduced by the approximation is limited to the deviations shown in Figure 2 for tissues with activities decaying monoexponentially. Figure 4 shows the time window for a reliable dose estimate in such tissues for different effective half-lives. Underestimates of the actually absorbed dose exceeding 10% occur only if the ratio tl/Teff is smaller than 0.75 or larger than 2.5. A measurement at tl = 96 h opens a window for reliable single-sample dosimetry for effective half-lives ranging from 38 to 128 h. In the present investigation, all tissues but 1 NET lesion (Teff = 130 h) met this requirement.
Ratio tl/Teff and deviation of ratio of 
and the factor 
from unity in dependence on time of measurement tl for monoexponential decay functions. Lines represent median effective half-lives of main component in kidneys (black bold), livers (green), spleens (yellow), and NET lesions (red bold) as well as shortest (black thin) and longest (red thin) measured half-life.
Especially absorbed doses in kidneys, which most often represent the dose-limiting organ in PRRT with 177Lu-DOTATATE/DOTATOC, must not be significantly underestimated by any dose approximation. Half-lives of the main decay component in kidneys were a median of 51 h and ranged from 40 to 68 h with 1 outlier at 106 h in a patient with only 1 kidney. These data are consistent to those reported by others (15,21). Given the actually observed half-lives, a measurement at tl = 72 h, which induces less than 10% error for 29 < Teff < 96 h, is also adequate for kidneys. However, the absorbed dose for the outlier at 106 h was underestimated by the 72-h uptake by 15%. Approximations based on earlier measurements led to an increased risk for higher deviations and a reduced correlation between dose estimates and actual doses. This finding is in agreement with the observation by Miederer et al. (16) that the correlation between the absorbed dose from PRRT with 177Lu and the kidney uptake measured in OctreoScan improved with increasing time of measurement.
Effective half-lives in livers, spleens, and NET lesions ranged from 52 to 130 h. Measurement should be performed after 96 or after 120 h (<10% error for 48 < Teff < 160 h) for these tissues.
Decay Functions
Most tissues evaluated in the present investigation showed a biexponential decay function. Consideration of potential deviations from a monoexponential decay is not uniform in PRRT dosimetry. Most often, any short-lived component is ignored (8,9,15,22–26), as is also the case with the single-sample approximation. A second decay component early after the activity administration can induce both over- and underestimates depending on its time integral in comparison to the time integral of the dominant long-lived component. Only few biexponential functions in our study showed biphasic clearance, with 2 components with positive amplitude. In most cases, a negative short-lived component was observed, indicating further accumulation early after the activity administration, as visible in the example shown in Figure 1. Such a negative component leads to an overestimation of the absorbed dose by the single-sample approximation. The contribution by the short-lived component in functions with biphasic clearance was generally too small to produce a considerable dose underestimate.
Uncertainties
Guerriero et al. (20) have analyzed the impacts of timing of measurements and the time–activity integration method. The authors found that these parameters strongly influenced the dosimetric calculations and that only a few measurements up to 4 d were often not sufficient for a uniquely defined fit of the activity time function. This was observable in the present study also; biexponential fitting often resulted in parameters with high uncertainty and covariance. The uncertainties can be reduced with later measurements. Heikkonen et al. (15), neglecting the effect of the short-lived component, acquired the last measurement after 168 h and reported a high correlation between simplified dosimetry from 2 measurements after 24 and 168 h and full dosimetry with an additional scan after 72 h. For highest accuracy, larger numbers of measurements especially during the first day and after a week or later after the activity administration would be desirable.
Kinetics were determined from counting rates observed in organ and background regions in planar images without correction for attenuation and scatter. The background and its change in time can be different in such regions, especially in the presence of overlapping tissues or variable intestinal activity, leading to a modification of the 
function and thus the time integral 
. Uncertainties would have been smaller for kinetics deduced from series of SPECT measurements.
It is, therefore, a shortcoming of the present study that the normalized uptake functions 
measured in the present study cannot be considered to be a gold standard reflecting the true kinetics with negligible error, although the accuracy of the approximation is not sensitive to the exact values of the fit parameters in 
. However, independent verification of the methodology in other patient groups is necessary before routine application.
Advantages of Single-Sample Dosimetry
The single-sample method enables dose estimates with even lower expense than earlier methods to simplify dosimetry. Our study indicates that it provides absorbed dose data with only minor additional errors compared with full dosimetry, thus enabling routinely performed dosimetry during each PRRT cycle with little additional effort.
In principle, activity data extracted from a single planar scan after 96 h can be used to approximate the absorbed dose. Most reports of absorbed doses from PRRT are based on planar imaging and conjugate-view activity quantification (13). However, accuracy of dosimetry from planar scans is limited, and large systematic errors are possible even after careful correction of background, scatter, and attenuation (13,22,25,27).
Any potential error from inadequate background or attenuation correction of planar images as well as unsuitable fitting and integration of 
due to insufficient data acquisition are excluded when the single-sample method is used together with high-quality SPECT/CT imaging with state-of-the-art quantitative reconstruction and partial-volume correction. Phantom measurements indicate that the activity concentration can be quantified by SPECT/CT with 177Lu with uncertainties of 10% or less (13,28).
It is an inherent benefit of the method that the absorbed dose estimate is largely independent of the effective half-life. While even skillful planar dosimetry with several measurements is usually unsuitable to identify differences in kinetics and absorbed doses in distinct areas of organs or tumors, an approximation based on a late quantitative SPECT/CT correctly reproduces the inhomogeneity of the distribution of absorbed doses. The counting rate information in the reconstructed SPECT images is proportional to the absorbed doses. Thus, apart from image distortion inherently present in SPECT, such as partial-volume effects, the data represent a 3-dimensional approximation of the dose distribution.
Dosimetric Procedure
Information on the best practice to produce reliable quantitative SPECT data for 177Lu were described by Ljungberg et al. (13). After patient imaging and quantitative SPECT data reconstruction, the counts in a volume of interest large enough to avoid major reconstruction errors were, with reasonable accuracy, proportional to the absorbed doses approximated by Equation 7. The product 
in Equation 7 is almost inversely proportional to the size of the tissue. It can be determined with the sphere model in OLINDA/EXM (29) that a constant factor of 
= 0.25 Gy·g/MBq/h (±2%) is valid for accumulating tissues with masses m between 10 g and 1 kg. The activity concentration in a volume of interest within such a tissue can be calculated from the counts CVOI and the mass mVOI in that volume, the sensitivity of the γ-camera ε (13), and the total acquisition duration tacq (frame duration times number of angular steps). The absorbed dose becomes:
Eq. 8The following example uses arbitrary values: a double-head γ-camera has a sensitivity of 10 cps/MBq per head (ε = 20 cps/MBq). SPECT is acquired with 60 angular steps (120 projections) and 20-s acquisition duration per angular step (tacq = 60 × 20 s = 1,200 s). According to Equation 8, the absorbed dose in a volume of interest with 40 g of healthy liver tissue with 84,000 counts in an acquisition at tl = 92 h after PRRT is 2 Gy.
CONCLUSION
The absorbed dose from PRRT with 177Lu to NET lesions and relevant abdominal organs can be deduced with reasonable accuracy from a single measurement 4 d after the activity administration by an approximation, which is substantiated by theory. All 177Lu-accumulating tissues showing monoexponential decay with effective half-lives between 38 and 128 h are well represented, with an error of less than 10%. Deviations from a monoexponential decay function, which most often are not considered in PRRT dosimetry due to a lack of data points, potentially induce additional errors.
The method to approximate absorbed doses from counting rates obtained by quantitative SPECT/CT imaging is easy to implement, provides a 3-dimensional dose map, and is a promising option to register cumulated doses actually absorbed in a series of PRRT treatment cycles with minor additional resources and effort.
DISCLOSURE
No potential conflict of interest relevant to this article was reported.
Footnotes
Published online Jun. 6, 2017.
- © 2018 by the Society of Nuclear Medicine and Molecular Imaging.
REFERENCES
- Received for publication February 7, 2017.
- Accepted for publication April 14, 2017.
Jump to section
Related Articles
Cited By...
- Impact of the Reference Multiple-Time-Point Dosimetry Protocol on the Validity of Single-Time-Point Dosimetry for [177Lu]Lu-PSMA-I&T Therapy
- Assessing Response to PSMA Radiopharmaceutical Therapies with Single SPECT Imaging at 24 Hours After Injection
- Safety and Efficacy of Para-Aminohippurate Coinfusion for Renal Protection During Peptide Receptor Radiotherapy in Patients with Neuroendocrine Tumors
- Single-Time-Point Renal Dosimetry Using Nonlinear Mixed-Effects Modeling and Population-Based Model Selection in [177Lu]Lu-PSMA-617 Therapy
- Impact of Single-Time-Point Estimates of 177Lu-PRRT Absorbed Doses on Patient Management: Validation of a Trained Multiple-Linear-Regression Model in 159 Patients and 477 Therapy Cycles
- Single-Time-Point Imaging for Dosimetry After [177Lu]Lu-DOTATATE: Accuracy of Existing Methods and Novel Data-Driven Models for Reducing Sensitivity to Time-Point Selection
- Toward Single-Time-Point Image-Based Dosimetry of 177Lu-PSMA-617 Therapy
- A Pipeline for Automated Voxel Dosimetry: Application in Patients with Multi-SPECT/CT Imaging After 177Lu-Peptide Receptor Radionuclide Therapy
- Dosimetry in Radiopharmaceutical Therapy
- Reply: Single-Time-Point Tumor Dosimetry Assuming Normal Distribution of Tumor Kinetics
- Single-Time-Point Tumor Dosimetry Assuming Normal Distribution of Tumor Kinetics
- Dosimetry in Clinical Radiopharmaceutical Therapy of Cancer: Practicality Versus Perfection in Current Practice
- Reimbursement Approaches for Radiopharmaceutical Dosimetry: Current Status and Future Opportunities
- Dosimetry for Radiopharmaceutical Therapy: Current Practices and Commercial Resources
- A Novel Time-Activity Information-Sharing Approach Using Nonlinear Mixed Models for Patient-Specific Dosimetry with Reduced Imaging Time Points: Application in SPECT/CT After 177Lu-DOTATATE
- Feasibility of Single-Time-Point Dosimetry for Radiopharmaceutical Therapies
- What You See Is Not What You Get: On the Accuracy of Voxel-Based Dosimetry in Molecular Radiotherapy
- Will SPECT/CT Cameras Soon Be Able to Display Absorbed Doses? Dosimetry from Single-Activity-Concentration Measurements
- Radiation Dosimetry in 177Lu-PSMA-617 Therapy Using a Single Posttreatment SPECT/CT Scan: A Novel Methodology to Generate Time- and Tissue-Specific Dose Factors
- Investigating PSMA-Targeted Radioligand Therapy Efficacy as a Function of Cellular PSMA Levels and Intratumoral PSMA Heterogeneity
- Side Effects of CXC-Chemokine Receptor 4-Directed Endoradiotherapy with Pentixather Before Hematopoietic Stem Cell Transplantation
- The Relevance of Dosimetry in Precision Medicine