Abstract
Doubling time (DT) is widely used for quantification of tumor growth rate. DT is usually determined from two volume estimations with measurement time intervals comparable with or shorter than DT. Clinical data show that the frequency distribution of DT in patients is positively skewed, with some very long DT values compared with the average DT. Growth rate can also be quantified using specific growth rate (SGR; %/d), equal to ln2/DT. The aim of this work was to compare DT and SGR as growth rate variables. Growth rate calculations were computer simulated for a tumor with DT of 100 days, measurement time interval of 1 to 200 days, and volume estimation uncertainty of 5% to 20%. Growth rate variables were determined and compared for previously published clinical data. The study showed that DT is not a suitable variable for tumor growth rate because (a) for short measurement time intervals, or high volume uncertainties, mean DT can either overestimate or underestimate the average growth rate; (b) DT is not defined if the consecutively estimated volumes are equal; and (c) the asymmetrical frequency distribution of DT makes it unsuitable for common statistical testing. In contrast, mean SGR and its equivalent DT give the correct values for average growth rate, SGR is defined for all tumor volume changes, and it has a symmetrical frequency distribution. SGR is also more accurate to use when discussing, for example, growth fraction, cell loss rate, and growth rate heterogeneities within the tumor. SGR should thus be used, instead of DT, to quantify tumor growth rate. [Cancer Res 2007;67(8):3970–5]
Introduction
Tumor growth rate is widely used for prognostic purposes and can quantify therapeutic effects of different treatment modalities (1–9). Tumor cells duplicate during cell cycle, which theoretically causes an exponential growth. When a tumor becomes larger, its growth rate might decrease and change to nonexponential growth model (e.g., the Gompertzian model; refs. 10–12). To observe growth retardation during volume increase, the tumor must be followed for a long period with several volume measurements. However, in clinical studies, volume estimations of nontreated tumors are usually available only for short measurement time intervals and tumor growth may well be described by an exponential model.
Tumor growth rate is usually characterized by the tumor volume doubling time (DT). The term DT was introduced 50 years ago and a graphical method was proposed for its estimation (13). For an exponentially growing tumor, the growth rate is proportional to its volume (V):
or
where SGR and t are the specific growth rate and time, respectively. The solution of Eq. B gives the well-known exponential growth equation:
and
DT is the time period when V2 = 2V1, then
The right side of Eq. E is equal to the formula, which was introduced by Schwartz (14) to estimate the growth rate of tumors.
To estimate DT of a tumor, its volume should be measured at least on two different occasions (Eq. E). The average growth rate for several tumors in one patient, or for the same tumor type in many patients, may be estimated by the arithmetic mean value of DT, DTm (15, 16). Clinical studies have shown that the frequency distribution of tumor DT for the same tumor type is usually asymmetrical in relation to the mean value (i.e., there are tumors with very long DT values compared with the mean value; positively skewed; ref. 17). Therefore, DTm does not indicate the average growth rate and DT is not suitable for common statistical testing. Some researchers have tried to obtain a symmetrical distribution by assuming a log-normal distribution of DT (4–6, 8, 17, 18). The average growth rate is then estimated by DTlog, calculated as the antilog to the arithmetic mean of the logarithms of DTs (17, 19, 20). The logarithm of DT, log(DT), is also proposed to be more suitable for statistical testing (17). DTlog is mathematically equal to geometric mean DT, DTgm, which is also used to estimate the average growth rate (4, 8, 21). In this article, we propose the equivalent DT, DTe, which is the value of DT calculated from Eq. E for the mean SGR. To our knowledge, the reason for the asymmetry in the frequency distribution of DT and the correctness of mean DT as an average growth rate estimator were not studied.
Some degree of uncertainty is included in all tumor volume measurements, which depends on several factors related to measurement technique/investigator (21–25). The uncertainty of the volume measurement propagates to the estimated DT and SGR values. If the tumor volume at the first measurement is underestimated, or overestimated at the second measurement, then the tumor growth rate will be overestimated, and vice versa. If the estimated volume at the second measurement is smaller than the first measurement, then the DT and the SGR of tumor will be negative. Both DT and SGR can be used to quantify the growth rate of a tumor, but the relationship between these variables is not linear (Eq. E). If the growth rates of a set of tumors are quantified using DT and the frequency distribution of DT and the mean DT are determined, then the results will not necessarily be similar to the frequency distribution of SGR and the mean SGR of the same set of tumors.
The aim of this study was to compare the frequency distribution and the mean value of different growth rate variables in relation to the true growth rate of tumor, based on computer simulations and evaluation of previously published clinical data. The final goal was to propose a better way to quantify tumor growth rates.
Materials and Methods
Monte Carlo simulations. To compare different methods of expressing tumor growth rate, DT, log(DT), and SGR, the frequency distribution of these variables and variation of their means, DTm, DTlog, and DTe, were analyzed by computer simulations for an exponentially growing tumor (constant SGR). Computer simulations were done using a Monte Carlo code, written in visual basic 6.0 (Microsoft), with the following assumptions: (a) The estimated volume of the tumor at the first measurement, V1, was a normally distributed random variable with the mean value V1true and the SD σv1. V1true was supposed to represent one unit of volume. (b) The true volume of the tumor at the second measurement, V2true, was calculated using Eq. C. The measured volume, V2, was a normally distributed random variable with the mean value V2true and the SD σv2. (c) The relative uncertainties of V1 and V2 were equal: (σv1 / V1) = (σv2 / V2). (d) Because V1 and V2 cannot be negative, the symmetry of each distribution was kept undisturbed by symmetrical truncation in the positive range. (e) DTtrue was assumed to be 100 days (SGR = 0.7%/d), which is a typical value of the growth rate of tumors in a clinical setting (26). Simulations were done for different measurement time intervals varying from 1 to 200 days. Each step was 1 day and the relative uncertainty of the volume estimation was 5%, 10%, or 20%.
For each time interval, 105 simulations were done. In each simulation, V1 and V2 were generated and SGRi and DTi were estimated for the range of i indices 1 to 105. For each time interval, DTm, DTlog, mean SGR
and
The relative uncertainty of SGR
If both sides of the above equation are divided by SGR and SGR on the right side is replaced from Eq. E, then
Because DT is inversely proportional to SGR, the simulation will generate unstable results for SGR close to zero. Therefore, SGR values between −0.0000693 and +0.0000693, corresponding to DTi with absolute values >10,000 days, were excluded from the calculations.
Because the logarithm of negative and zero values are undefined, the following method was used in the calculation of DTlog in the presence of negative DTi values: the absolute value of the minimum possible DT (−10,000) plus one (i.e., 10,001) was added to all DT values, and the mean of their logarithms was calculated. Thereafter, DTlog was derived by subtracting 10,001 from the obtained mean value. To investigate how the exclusion of negative growth rate values can influence the average growth rate estimators, DTm, DTlog, and DTe, the simulation was then repeated excluding SGR values less than +0.0000693.
Evaluation of clinical data. DTm, DTlog, and DTe were calculated for several types of tumors using quantitative data from previously published clinical studies (Table 1; refs. 2, 3, 9, 21, 27–30). The results were then compared with computer simulations.
Results
Figure 1 shows the simulated frequency distributions of DT (Fig. 1A), log(DT) (Fig. 1B), and SGR (Fig. 1C) for different time intervals (1, 5, 10, 50, 100, and 200 days), when the relative uncertainty of the volume measurement was 10%. For a time interval of 200 days (2 DTtrue), all DT values were positive and the frequency distribution of DT was symmetrical and centered at DT of 100 days (Fig. 1A). When the time interval was 100 days (1 DTtrue), the frequency distribution of DT was positively skewed and the peak shifted toward lower DT values. When the time interval was 50 days (0.5 DTtrue), the peak shifted more toward lower DT values and negative DT values appeared in the data as a very small peak in the negative range. The peak in the negative range increased further with decreasing time interval. With a 1-day time interval, the two peaks were very close and symmetrical in relation to zero and appeared as a single peak centered at zero (Fig. 1A). Therefore, mean DT was close to zero for very short time intervals (Fig. 2A). Theoretically, when the time interval approaches zero, the position of two peaks asymptotically approaches zero with a height of infinity. If negative values of DT were excluded, the peaks on the negative side of the frequency distribution of DT disappeared. Variations in the frequency distribution of log(DT) (Fig. 1B) were comparable with that of DT. For the time intervals of 200 and 100 days, all DT values were positive and only one peak appeared in the frequency distribution of log(DT) centered at 4.6 (=log 100) for 200 days and slightly shifted to the left for 100 days. For shorter time intervals, where negative DT values appeared in data, the peak shifted more to lower values in relation to 4.6, when negative DT values were excluded (Fig. 1B). When negative DT values were included for 50-, 10-, 5-, and 1-day time intervals, the symmetry point was shifted to 9.21 (=log 10,001; see Materials and Methods) comparable with zero in the frequency distribution of DT (Fig. 1B,, inset). For a 1-day time interval, the two peaks looked like a single peak centered at 9.21. Therefore, also DTlog was close to zero for very short time intervals (Fig. 2). The frequency distribution of SGR was symmetrical for all time intervals studied. The mean SGR was equal to the true SGR (0.7%/d) and its uncertainty increased with decreased time interval. The expected uncertainty of SGR from Eq. I and the calculated uncertainty of SGR from the simulations were well correlated (R2 > 0.999).
The results of the computer simulations of DTm, DTlog, and DTe are shown in Fig. 2. When the time interval was very long compared with DTtrue, all DT estimators were equal to DTtrue of the tumor (Fig. 2A). When the time interval decreased, DTm overestimated DTtrue with a maximum deviation of ∼30%. For very short time intervals compared with DTtrue, DTm underestimated DTtrue and approached zero for time intervals down to a few days. DTlog showed a similar variation as DTm, with a maximum overestimation of ∼20% but a larger underestimation than DTm for short time intervals. Neglecting small fluctuations at very short time intervals, DTe was equal to DTtrue for all time intervals studied (Fig. 2A). When the negative growth rate values were excluded, DTm and DTlog followed a similar shape as when the negative values were included (i.e., overestimation up to a maximum and then decreasing with decreasing time interval; Fig. 2B). However, the ranges of deviation from DTtrue were different. When negative values were excluded, DTm was much higher, whereas DTlog was closer to the results when negative values were included. Furthermore, DTe decreased with decreasing time interval and was lower than DTlog, which in turn was lower than DTtrue when negative values were excluded. Then, DTe approached zero for very short time intervals down to a few days.
Figure 3 shows the values of DTm, DTlog, and DTe as a function of the relative uncertainty of SGR. It shows that, regardless of the uncertainty of the volume estimation, DTm, DTlog, and DTe gave similar results for similar uncertainties of SGR.
DTm, DTlog, and DTe values estimated from the previously published clinical data on several types of tumors are presented in Fig. 4. The measurement time intervals varied between 13 and 1,212 days. The estimated DTs from these articles were between −1,350 and 1,172 days. The only study containing negative growth rates was that of adenocarcinoma and squamous cell lung carcinoma (31). For all studies, including only positive growth rates, DTe was lower than DTlog, which was lower than DTm (Fig. 4A). On average, DTlog and DTm were 25% (range, 3–88%) and 76% (range, 6–317%) higher than DTe, respectively. If the negative growth rates were included, negative DTm was obtained, whereas DTe was still positive (data not shown).
The SGR values from clinical data are summarized in Fig. 4B. Because SGR and DT are reciprocally related (Eq. E), a higher SGR value in Fig. 4B corresponds to a shorter DTe in Fig. 4A and vice versa. Such trend was not always seen for DTm and DTlog values because they may overestimate or underestimate the true DT of tumors depending on volume measurement uncertainties and the time interval (Figs. 2 and 3).
Discussion
This study clearly shows the importance of selecting a correct quantification method for proper evaluation of tumor growth rate. There are at least three advantages of using SGR instead of DT:
(1) Reduced errors. Errors are encountered with DTm and DTlog, which may cause either overestimation or underestimation of DTtrue. DTm can correctly estimate the average growth rate of tumors when the frequency distribution of DT is symmetric, that is when the uncertainty of tumor volume estimation is relatively low (e.g., 10%), or the time interval is considerably longer than DTtrue (e.g., 2 DTtrue). With increasing volume uncertainty, or decreasing time interval, the frequency distribution of DT becomes positively skewed. DTm will then overestimate DTtrue due to this asymmetry and the size of the error depends on the uncertainty of volume estimations and the measurement time interval. The logarithmic transformation of DT will decrease this asymmetry but is not symmetrical either. With time intervals shorter than DTtrue, negative DT values appear and increase in number with decreasing time interval. Therefore, the overestimation of DTtrue by DTm is limited to a maximum of ∼30%. Further decrease of the time interval causes decrease of DTm toward zero. Negative growth rates would have been excluded if the tumor was not detected at the second investigation. If negative growth rates are excluded, none of the average growth rate estimators will indicate the true growth rate (Figs. 2B and 3B). Assuming that the volume of undetected tumor equals the detection limit, one may include negative values of DT for proper evaluation. However, DTm and DTlog may still give incorrect estimates of the true growth rate (Figs. 2 and 3). DTm and DTlog can thus overestimate or underestimate DTtrue, depending on the degree of uncertainty of the volume measurement and the time interval. DTe is the only estimator that can give the true average growth rate of tumors. All measurement results, including negative SGR values, should be included in the calculation of mean SGR (Fig. 2A). Otherwise, DTe will deviate from DTtrue, as DTm and DTlog do (Fig. 2B). It should be noted that the inclusion of negative values does not mean that such values must exist. Winer-Muram et al. (30) observed the erroneous estimation of the average growth rate with DTm and made a better estimation by the reciprocal of the average of reciprocals of DT values. The results of that method are comparable with those using SGR as growth rate variable.
(2) No problem when V2 ≈ V1. If the estimated volume at the second measurement is only slightly larger, or smaller, than the volume at the first measurement, then the calculated DT will result in a very high positive or negative value (e.g., −1,350 days for lung adenocarcinoma; ref. 30). Furthermore, DT is not defined if the estimated volumes are equal (V2 = V1 in Eq. E). No problem is encountered in the calculation of SGR and DTe when V2 = V1.
(3) Suitable frequency distribution and statistical testing. When DT is used for quantification of tumor growth rate at time intervals used clinically, with an uncertainty in the volume measurement, an asymmetrical frequency distribution of DT is obtained (Fig. 1). This fact was shown by the computer simulations and is similar to observed DT distributions from clinical studies (17). The frequency distribution of the logarithm of DT is not symmetrical either. Most of the common statistical tests are based on the assumption of normally distributed variables. Therefore, DT and log(DT) are not suitable variables for such statistical analyses. SGR is the growth rate estimator, which is least influenced by uncertainties of the measurement procedure. Its frequency distribution is symmetrical for an exponentially growing tumor (i.e., SGR is suitable for common statistical tests).
The difference between exponential and nonexponential growth models is that SGR is constant for exponential growth, whereas it is related to tumor volume for nonexponential growth. Therefore, to make statistical analysis, the SGR values can be normalized for a specific volume with access to the nonexponential growth variables of the tumor. Nevertheless, dependency or nondependency of SGR on tumor volume, or time, does not degrade the superiority of SGR over DT.
Several assumptions were made in the computer simulations done in this study. The relative uncertainty of volume was assumed to be equal for all tumors, realistic for some applications (25), but is probably not the case for all imaging modalities. However, the effect of volume uncertainty is modulated by the time interval (Eq. I). The variable, which generates the differences between different growth rate estimators, is the relative uncertainty of SGR (Fig. 3). If two different sets of volume uncertainties and time intervals generate the same relative uncertainty of SGR, the estimated values of DTm, DTlog, and DTe will be similar in both settings (Fig. 3). Therefore, different relative volume uncertainties at different volumes will only change the relative uncertainty of SGR without affecting the main results and conclusions of this study. The same discussion is valid for the assumption of DT equal to 100 days. According to Eq. I, the SGR uncertainty depends on the ratio between the measurement time interval and DT. Therefore, the results are also valid for tumors with different DTs but with equal ratios between time interval and DT. Furthermore, the exclusion of SGR values between −0.0000693 and +0.0000693 (i.e., truncation) seems realistic because it practically excludes zero values of SGR, which may be excluded by investigators as nongrowing tumors (i.e., DT = ± ∞; ref. 30). When the simulation was repeated without truncation, DTm generated extremely large positive or negative values, whereas there were no problems to calculate SGR and DTe without truncation.
Differences in growth rates of tumors is mainly a result of different growth fractions (GF; fraction of tumor volume that consists of proliferating cells) and cell loss rates (CLR). Thus, the duration of cell cycle does not play a major role in the varying kinetics of tumor growth (26, 31). Such a growth pattern can quantitatively be described using SGR. If the absolute value of CLR is added to the SGR of tumor volume, the SGR of the entire tumor in the absence of cell loss is obtained. The SGR of the proliferating cells can then be obtained by dividing the result by GF. In general, the volumetric SGR of a tumor with polyclonal cell population (a heterogeneous SGR distribution within the tumor) is the mean of SGR values weighted by the fraction of each cell component, including stromal and tumor cell populations (i.e., for a tumor with n cell components: V = V1+V2…+Vn, then V·SGRV = V1·SGR1+V2·SGR2…+Vn·SGRn). Similar to DTm within a population, calculation of volumetric growth rate of tumor as mean DT of its components will result in erroneous estimations of the growth rates.
The SGR values for the clinical data ranged between 0.5%/d and 3%/d. Our recalculations of data from these studies clearly show that the DTe values always were the lowest ones for each study and tumor type. These results are comparable with simulation results in Fig. 3B. This means that when using routine methods today, DT might be overestimated; hence, the tumor might in reality grow faster than estimated when using DTm and DTlog.
In conclusion, to correctly quantify the growth rate of tumors, the variable SGR must be used. The average growth rate must also be estimated by the mean SGR or DTe. The uncertainty of SGR can be reduced with increasing measurement time interval or decreasing volume uncertainties. In addition, SGR is a suitable variable for common statistical testing based on the assumption of normally distributed variables.
Acknowledgments
Grant support: Swedish Cancer Society grant 4956 and the King Gustav V Jubilee Clinic Cancer Research Foundation (Göteborg, Sweden).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
We thank Prof. Ragnar Hultborn (Department of Oncology, Göteborg University, Göteborg, Sweden) for valuable discussion of the manuscript.