Analysis of inter-patient variations in tumour growth rate
- Esmaeil Mehrara^{1, 2}Email author and
- Eva Forssell-Aronsson^{1, 2}
https://doi.org/10.1186/1742-4682-11-21
© Mehrara and Forssell-Aronsson; licensee BioMed Central Ltd. 2014
Received: 6 April 2014
Accepted: 13 May 2014
Published: 20 May 2014
Abstract
Purpose
Inter-patient variations in tumour growth rate are usually interpreted as biological heterogeneity among patients due to, e.g., genetic variability. However, these variations might be a result of non-exponential, e.g. the Gompertzian, tumour growth kinetics. The aim was to study if the natural tumour growth deceleration, i.e. non-exponential growth, is a dominant factor in such variations.
Materials and methods
The correlation between specific growth rate (SGR) and the logarithm of tumour volume, Ln(V), was calculated for tumours in patients with meningioma, hepatocellular carcinoma, pancreatic carcinoma, primary lung cancer, post-chemotherapy regrowth of non-small cell lung cancer (NSCLC), and in nude mice transplanted with human midgut carcinoid GOT1, a tumour group which is biologically more homogeneous than patient groups.
Results
The correlation between SGR and Ln(V) was statistically significant for meningioma, post-chemotherapy regrowth of NSCLC, and the mouse model, but not for any other patient groups or subgroups, based on differentiation and clinical stage.
Conclusion
This method can be used to evaluate the homogeneity of tumour growth kinetics among patients. Homogeneity of post-chemotherapy regrowth pattern of NSCLC suggests that, in contrast to untreated tumours, the remaining resistant cells or stem cells (if exist) might have similar biological characteristics among these patients.
Keywords
Tumour Growth rate deceleration Cancer Modelling GompertzianIntroduction
Mathematical modelling of tumour growth can provide not only key insights into tumour biology but also tools for, e.g., optimization of screening programs, cancer patient prognosis[1], scheduling of chemotherapy[2], and assessment of tumour spread[3, 4]. For example, Norton et al. showed that, considering the mathematical growth model of breast cancer tumours, patients must be treated with condensed-dose chemotherapy, and a clinical trial showed the significant benefit for the patients treated with the new method compared with the patients treated with the standard treatment[5, 6].
where V_{0} and V are the volume of tumour at t = t_{0} and t, respectively.
where SGR_{0} is the initial growth rate of tumour at t = t_{0} and λ is the growth deceleration constant.
The above equation enabled us estimating the growth model, including Gompertzian growth deceleration constant in Equation B (λ) and formation times of metastases in individual patients.
Tumour response to a specific treatment varies largely among patients. Beside other biological factors, tumour growth kinetics is important for how a tumour responds to therapy[21, 22]. Tumour volume at base line and its growth rate have both been shown to be correlated with tumour response to therapy. Different factors are responsible for the varying tumour growth rates among patients, e.g., genetic factors, microenvironment in host tissue and growth deceleration as tumours grow. Furthermore, measurement uncertainties may also influence the reported variation in growth rate.
According to Equation C, a significant correlation between SGR and ln(V) in a group of tumours indicates that tumours probably follow the same growth curve and the difference in their growth rate is a result of difference in their volume. In other words, such a correlation is a measure of similarity of growth pattern among tumours, e.g., in a group of patients. Then, a unique curve can describe the growth of all tumours in the group, e.g., the curve that we derived for growth of liver metastases from a primary ileum carcinoid in a previous article[20].
The aim of this study was to assess the homogeneity of the growth kinetics of tumours of the same type in a population of patients. Data from different groups of patients, including meningioma, hepatocellular carcinoma, pancreatic carcinoma, and primary lung cancer, were analysed. The model was also applied to data from a mouse model bearing transplanted human midgut carcinoid GOT1 tumours, a tumour group which is biologically more homogeneous than patient groups.
Materials and methods
Calculation method
SGR was calculated for each pair of tumour volume measurements using Equation A. Correlation between SGR and ln(V) was calculated between all SGR values and their corresponding tumour volume, i.e., the geometric mean of the two volumes used for calculation of SGR.
Clinical data
Data from clinical studies were retrieved from the literature based on the availability of tumour volume estimates and corresponding measurement time intervals. Tumour volumes were estimated by delineating tumours in computed tomography (CT) and/or magnetic resonance imaging (MRI) slices and multiplying the measured area with slice thickness. Correlation between SGR and the logarithm of the volume of tumours was calculated for the following types of tumours.
Lung cancer
Data on the growth of non-small cell lung cancer (NSCLC) tumours in 18 patients were used[23]. Tumour growth was measured between the end of induction chemotherapy and the start of radiation therapy. The study showed that the regrowth of tumours after induction therapy (mean DT = 46 days, median DT = 29 days) was much faster than the untreated tumour growth rate found in the literature (mean DT range: 102–452 days). Of the potentially curable patients 41% became incurable in the waiting period between chemo- and radio-therapy. Tumour DT was shorter for smaller tumours compared to large tumour. Considering the fast regrowth of NSCLC tumours after induction chemotherapy, they recommended diminishing the time interval between chemo- and radiotherapy to as short as possible[23].
Another set of data was found for the growth rate and characteristics of small peripheral lung tumours as they appear on CT images. The tumour types included were rapidly growing (DT < 150 days) small cell lung cancer, adenocarcinoma, and squamous cell carcinoma tumours[24].
Pancreatic carcinoma
Data for untreated pancreatic carcinoma in nine patients who underwent serial examinations by helical computed tomography was used for calculations[25]. The mean DT of the nine primary lesions was 159 days (median 144 days). Survival time was significantly correlated with DT.
Hepatocellular carcinoma
The growth of hepatocellular carcinoma (HCC) was studied in 11 untreated patients where serial CT or MRI images were available (16 tumours in total)[26]. Calculated DT value range was 17.5-541 days and the mean DT was 127 days. DT was related to baseline volume as DT = 114 × (baseline volume)^{0.14}. This study showed that smaller tumours had shorter DT values, i.e. grew faster, than larger tumours and, therefore, may require shorter follow-up time to observe progression[26].
Correlation between the specific growth rate, SGR, and the logarithm of tumour volume in groups of patients diagnosed with the same type of tumour
Tumour type | Reference | n | R^{2} | p (1-tailed) | Median time interval (days) | Median time interval/DTe | Average SGR (%/d) | Relative uncertainty of SGR (%) | DTe (days) | Relative uncertainty of volume estimation (%) |
---|---|---|---|---|---|---|---|---|---|---|
Meningioma | Nakamura et al.[29] | 41 | 0.2424 | 0.0005 | 1230 | 0.71 | 0.04 | 105 | 1733 | 37 |
Hepatocellular carcinoma | Nakajima et al.[27] | 34 | 0.038 | NS (0.13) | 128 | 2.4 | 1.30 | 73 | 53 | 86 |
Saito et al.[28] | 21 | 0.0623 | NS (0.14) | 146 | 1.05 | 0.50 | 47 | 139 | 24 | |
Taouli et al.[26] | 16 | 0.0014 | NS (0.89) | NA | NA | 0.90 | 115 | 77 | NA | |
Pancreatic carcinoma | Furukawa et al.[25] | 9 | 0.0248 | NS (0.34) | 295 | 2.13 | 0.50 | 58 | 139 | 60 |
Primary Lung Cancer | Wang et al.[24] | 12 | 0.1619 | NS (0.19) | 365 | 4.21 | 0.80 | 36 | 87 | 74 |
El Sharouni et al.[23] | 18 | 0.2713 | 0.026 | 48 | 1.94 | 2.80 | 73 | 25 | 69 | |
REGROWTH |
Another data set on HCC included tumours less than 3 cm in diameter at first observation in 21 patients[28]. The natural progression of each lesion (DT) was observed by ultrasonography[28]. They concluded that data from cell kinetic parameters and histological grade are useful for estimating the natural growth rate of HCC. We determined SGR and tumour volumes based on data in Table 1 in that article[28].
Meningioma
The natural history of incidental meningiomas was studied in asymptomatic patients[29]. Average tumour volume was 9 cm^{3} and DT ranged from 1.27 to 143 years (mean, 21.6 y). DT was shorter for younger patients, but was not correlated with tumour size. We determined SGR and tumour volumes based on data from that article[29]. The relatively large number of data included enabled us dividing tumours into two groups and comparing SGR between large and small tumours.
Mouse model of human midgut carcinoid, GOT1
In order to test the validity of our method, the correlation between SGR and ln(V) was studied on tumours in the GOT1 nude mouse model (human midgut carcinoid tumours), which should biologically be more homogeneous compared with patients[30]. Data from 57 mice was gathered and the natural growth of tumours was followed for several weeks before any treatment.
Statistical analyses
To test the significance of correlations, the t-values for Pearson correlation were calculated and then converted to p-values using one tailed T.DIST function in Microsoft Excel. P-value < 0.05 was assumed to be statistically significant. All patient data used in this study were retrieved retrospectively from already published articles. This type of information is exempt from ethical approval. The animal experiments were approved by the Ethical Committee for Animal Research at University of Gothenburg.
Results
Correlation between the specific growth rate, SGR, and the logarithm of tumour volume in hepatocellular carcinoma patients
Tumour type | Group | n | R^{2} | p-value |
---|---|---|---|---|
Hepatocellular carcinoma, | ||||
Nakajima et al.[27] | WD | 19 | 0.001 | 0.5 |
MD | 9 | 0.063 | 0.3 | |
PD | 6 | 0.441 | 0.1 | |
CS I | 17 | 0.040 | 0.2 | |
CS II | 15 | 0.009 | 0.4 | |
WD & CS I | 8 | 0.030 | 0.3 | |
WD & CS II | 10 | 0.001 | 0.5 |
Correlation between the specific growth rate, SGR, and the logarithm of tumour volume in hepatocellular carcinoma patients
Tumour type | Group | n | R^{2} | p-value |
---|---|---|---|---|
Hepatocellular carcinoma, | ||||
Saito et al.[28] | WD | 15 | 0.088 | 0.1 |
MD | 6 | 0.007 | 0.4 |
The correlation between SGR of tumour and the logarithm of its volume was statistically significant for data from nude mice bearing GOT1 tumours (p < 1E-11) (Figure 1).
Discussion
Accurate quantification and analyses of the tumour growth rates is essential for understanding the biological variance of human cancers[31]. Differences in the observed growth rates of tumours of the same type in a population of patients can be due to: (a) measurement uncertainties, (b) growth deceleration with increasing tumour volume, or (c) other biological differences between tumours, e.g., their location and microenvironment. With regard to measurement uncertainties, we previously showed that SGR is the tumour growth rate measure that is least influenced by measurement errors, compared to other growth rate measures, e.g., tumour volume doubling time (DT)[7]. We also showed that SGR is least affected by variances due to biological factors[8]. In the present study, focused on the growth deceleration in tumours, we used the relation between SGR and the logarithm of tumour volume to assess the contribution of growth rate decline in the observed variances in tumour growth rate found in the selected clinical studies. It should be noted that the limited amount of data on natural tumour growth available makes this type of studies difficult.
Expected highly significant negative correlation between SGR and the logarithm of tumour volume in the mouse model (which should be biologically more homogeneous compared to tumours in groups of patients) showed that the presented method is a useful tool to assess if growth deceleration is an important factor influencing difference in tumour growth rates observed among patients.
A significant correlation between SGR and the logarithm of tumour volume in a group of tumours also indicates that the smaller tumours represent the growth of larger tumours when they were of small size and vice versa. This provides the possibility for further development of mathematical models for elaborating this correlation in more accurate efficacy assessment of new drugs or combination of treatments. However, lack of correlation between SGR and the logarithm of tumour volume indicates that biological factors other than growth deceleration are more important for explaining the differences in the tumour growth rate observed in a population of patients with the same tumour type. These tumours may grow exponentially with different growth rates, or according to the Gompertzian model, and the model constants, SGR_{0} and λ, are heterogeneously distributed among tumours[20].
In this study, the correlation between SGR and the logarithm of tumour volume was statistically significant for the growth of tumours in meningioma patients. Further analysis by dividing the material into small and large tumours also supported this result. A similar growth model observed for a group of patients with one tumour type corroborates that the response rate in this group might be a suitable measure to assess the efficacy of novel treatments. The correlation between tumour growth rate and the logarithm of its volume was, however, not statistically significant for the studied patient cohorts with hepatocellular carcinoma, pancreatic carcinoma, and primary lung cancer. This result was expected for the primary lung cancer cohort, since that included different types of lung cancer. A lack of correlation indicates that the growth of these types of tumours varies between patients with similar tumour type and, therefore, growth rate should be taken into account as an independent variable in efficacy assessment of treatment of patients with these types of tumours both in clinical trials and for individualized therapy planning.
However, in the present study we were in general not able to include factors such as histologic grade and differentiation in our analyses, due to the limited number of tumours/patients or lack of information for other types of tumours. The only exception was for hepatocellular carcinoma, where subgroups could be categorized according to their differentiation level and clinical stage. However the correlation between tumour growth rate and the logarithm of its volume was not statistically significant for any of the two patient cohorts studied. It would otherwise be possible that stratification of the tumours according to, e.g., histopathological information would have given better correlation within subgroups of each tumour type. The result of the animal study indicates such a situation.
Nevertheless, the significant correlation between tumour growth rate and the logarithm of its volume for regrowth of primary lung cancer tumours after chemotherapy was interesting, because this can be interpreted as homogeneity of these tumours among patients. Post-chemotherapy tumours consisting of resistant cancer cells or cancer stem cell clones (if exist) of non-small cell lung cancer might be more homogeneous in terms of cellular and histological characteristics and, therefore, respond more similarly to therapy. Such information is valuable for further treatments in, e.g., neoadjuvant therapy.
The main difference between meningioma and other types of tumours in the present study is that meningiomas are benign. It is possible that the reason for meningioma being the only tumour type that gave significant correlation between tumour size and tumour growth rate in the patient cohorts is that it is benign, with low proliferation rate and histologically more homogeneous than malignant tumours.
In the original study on meningioma, no correlation between tumour DT and its volume was found[29]. On the other hand, we found strong correlation between SGR and the logarithm of tumour volume, a correlation that, by definition, is the accurate measure of growth rate deceleration as tumour grows. The significant difference between growth rate of small and large tumours also supported our result. This is in line with our previous study where we showed that using DT for quantification of tumour growth rate can result in wrong conclusions[8]. This emphasizes again the need for developing new accurate mathematical tools to analyse clinical data.
In conclusion, the presented method, i.e. estimation of the Gompertzian growth deceleration constant using limited clinical data, can be used to evaluate the homogeneity of tumour growth pattern among patients. Tumour growth kinetics was largely heterogeneous among patients with the same type of tumour, except for the meningioma group. This implies that the tumour growth kinetics in each patient should be considered in efficacy assessment of new drugs and for optimization of treatment in individual patients using, e.g., the tumour response model in ref[21, 22]. Furthermore, homogeneity of post-chemotherapy regrowth pattern of non-small cell lung cancer tumours suggests that, in contrast to pre-treatment tumours, the remaining resistant cancer cells or cancer stem cells (if exist) might have similar biological characteristics among these patients, a factor that should be valuable to consider in, e.g., neoadjuvant therapy.
Declarations
Acknowledgements
This study was supported by grants from the Swedish Cancer Society, the Swedish Research Council, BioCARE, and a National Strategic Research Program at University of Gothenburg, the King Gustav V Jubilee Clinic Cancer Research Foundation, the Assar Gabrielsson Cancer Research Foundation, and Lions Cancerfond Väst, Sweden. The authors thank Prof. Ragnar Hultborn, Department of Oncology, University of Gothenburg, for valuable discussions.
Authors’ Affiliations
References
- Bassukas ID, Hofmockel G, Tsatalpas P, Eberle V, Maurer-Schultze B: Prognostic relevance of the intrinsic growth deceleration of the first passage xenografts of human renal cell carcinomas. Cancer. 1996, 78: 2170-2172. 10.1002/(SICI)1097-0142(19961115)78:10<2170::AID-CNCR19>3.0.CO;2-W.View ArticlePubMedGoogle Scholar
- Norton L: A Gompertzian model of human breast cancer growth. Cancer Res. 1988, 48: 7067-7071.PubMedGoogle Scholar
- Withers HR, Lee SP: Modeling growth kinetics and statistical distribution of oligometastases. Semin Radiat Oncol. 2006, 16: 111-119. 10.1016/j.semradonc.2005.12.006.View ArticlePubMedGoogle Scholar
- Iwata K, Kawasaki K, Shigesada N: A dynamical model for the growth and size distribution of multiple metastatic tumors. J Theor Biol. 2000, 203: 177-186. 10.1006/jtbi.2000.1075.View ArticlePubMedGoogle Scholar
- Piccart-Gebhart MJ: Mathematics and oncology: a match for life?. J Clin Oncol. 2003, 21: 1425-1428. 10.1200/JCO.2003.12.068.View ArticlePubMedGoogle Scholar
- Schmidt C: The Gompertzian view: Norton honored for role in establishing cancer treatment approach. J Natl Cancer Inst. 2004, 96: 1492-1493. 10.1093/jnci/96.20.1492.View ArticlePubMedGoogle Scholar
- Mehrara E, Forssell-Aronsson E, Ahlman H, Bernhardt P: Specific growth rate versus doubling time for quantitative characterization of tumor growth rate. Cancer Res. 2007, 67: 3970-3975. 10.1158/0008-5472.CAN-06-3822.View ArticlePubMedGoogle Scholar
- Mehrara E, Forssell-Aronsson E, Ahlman H, Bernhardt P: Quantitative analysis of tumor growth rate and changes in tumor marker level: specific growth rate versus doubling time. Acta Oncol. 2009, 48: 591-597. 10.1080/02841860802616736.View ArticlePubMedGoogle Scholar
- Afenya EK, Calderon CP: Diverse ideas on the growth kinetics of disseminated cancer cells. Bull Math Biol. 2000, 62: 527-542. 10.1006/bulm.1999.0165.View ArticlePubMedGoogle Scholar
- Bajzer Z: Gompertzian growth as a self-similar and allometric process. Growth Dev Aging. 1999, 63: 3-11.PubMedGoogle Scholar
- Hart D, Shochat E, Agur Z: The growth law of primary breast cancer as inferred from mammography screening trials data. Br J Cancer. 1998, 78: 382-387. 10.1038/bjc.1998.503.PubMed CentralView ArticlePubMedGoogle Scholar
- Wennerberg J, Willen R, Trope C: Changes in histology and cell kinetics during the growth course of xenografted squamous cell carcinoma. Arch Otolaryngol Head Neck Surg. 1988, 114: 781-787. 10.1001/archotol.1988.01860190085030.View ArticlePubMedGoogle Scholar
- Spratt JA, von Fournier D, Spratt JS, Weber EE: Decelerating growth and human breast cancer. Cancer. 1993, 71: 2013-2019. 10.1002/1097-0142(19930315)71:6<2013::AID-CNCR2820710615>3.0.CO;2-V.View ArticlePubMedGoogle Scholar
- Spratt JS, Meyer JS, Spratt JA: Rates of growth of human neoplasms: part II. J Surg Oncol. 1996, 61: 68-83. 10.1002/1096-9098(199601)61:1<68::AID-JSO2930610102>3.0.CO;2-E.View ArticlePubMedGoogle Scholar
- Pavelic ZP, Porter CW, Allen LM, Mihich E: Cell population kinetics of fast- and slow-growing transplantable tumors derived from spontaneous mammary tumors of the DBA/2 Ha-DD mouse. Cancer Res. 1978, 38: 1533-1538.PubMedGoogle Scholar
- Bassukas ID, Maurer-Schultze B: Mechanism of growth retardation of the adenocarcinoma EO 771. Radiat Environ Biophys. 1987, 26: 125-141. 10.1007/BF01211407.View ArticlePubMedGoogle Scholar
- DeWys WD: Studies correlating the growth rate of a tumor and its metastases and providing evidence for tumor-related systemic growth-retarding factors. Cancer Res. 1972, 32: 374-379.PubMedGoogle Scholar
- Prehn RT: The inhibition of tumor growth by tumor mass. Cancer Res. 1991, 51: 2-4.PubMedGoogle Scholar
- Araujo RP, McElwain DL: A history of the study of solid tumour growth: the contribution of mathematical modelling. Bull Math Biol. 2004, 66: 1039-1091. 10.1016/j.bulm.2003.11.002.View ArticlePubMedGoogle Scholar
- Mehrara E, Forssell-Aronsson E, Johanson V, Kolby L, Hultborn R, Bernhardt P: A new method to estimate parameters of the growth model for metastatic tumours. Theor Biol Med Model. 2013, 10: 31-10.1186/1742-4682-10-31.PubMed CentralView ArticlePubMedGoogle Scholar
- Mehrara E, Forssell-Aronsson E, Bernhardt P: Objective assessment of tumour response to therapy based on tumour growth kinetics. Brit J Cancer. 2011, 105: 682-686. 10.1038/bjc.2011.276.PubMed CentralView ArticlePubMedGoogle Scholar
- Mehrara E, Forssell-Aronsson E, Bernhardt P: Objective assessment of tumour response to therapy based on tumour growth kinetics (vol 105, pg 682, 2011). Brit J Cancer. 2011, 105: 1468-1468. 10.1038/bjc.2011.436.PubMed CentralView ArticleGoogle Scholar
- El Sharouni SY, Kal HB, Battermann JJ: Accelerated regrowth of non-small-cell lung tumours after induction chemotherapy. Br J Cancer. 2003, 89: 2184-2189. 10.1038/sj.bjc.6601418.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang JC, Sone S, Feng L, Yang ZG, Takashima S, Maruyama Y, Hasegawa M, Kawakami S, Honda T, Yamanda T: Rapidly growing small peripheral lung cancers detected by screening CT: correlation between radiological appearance and pathological features. Br J Radiol. 2000, 73: 930-937.View ArticlePubMedGoogle Scholar
- Furukawa H, Iwata R, Moriyama N: Growth rate of pancreatic adenocarcinoma: initial clinical experience. Pancreas. 2001, 22: 366-369. 10.1097/00006676-200105000-00005.View ArticlePubMedGoogle Scholar
- Taouli B, Goh JS, Lu Y, Qayyum A, Yeh BM, Merriman RB, Coakley FV: Growth rate of hepatocellular carcinoma: evaluation with serial computed tomography or magnetic resonance imaging. J Comput Assist Tomogr. 2005, 29: 425-429. 10.1097/01.rct.0000164036.85327.05.View ArticlePubMedGoogle Scholar
- Nakajima T, Moriguchi M, Mitsumoto Y, Katagishi T, Kimura H, Shintani H, Deguchi T, Okanoue T, Kagawa K, Ashihara T: Simple tumor profile chart based on cell kinetic parameters and histologic grade is useful for estimating the natural growth rate of hepatocellular carcinoma. Hum Pathol. 2002, 33: 92-99. 10.1053/hupa.2002.30194.View ArticlePubMedGoogle Scholar
- Saito Y, Matsuzaki Y, Doi M, Sugitani T, Chiba T, Abei M, Shoda J, Tanaka N: Multiple regression analysis for assessing the growth of small hepatocellular carcinoma: the MIB-1 labeling index is the most effective parameter. J Gastroenterol. 1998, 33: 229-235. 10.1007/s005350050075.View ArticlePubMedGoogle Scholar
- Nakamura M, Roser F, Michel J, Jacobs C, Samii M: The natural history of incidental meningiomas. Neurosurgery. 2003, 53: 62-70. 10.1227/01.NEU.0000068730.76856.58. discussion 70–61View ArticlePubMedGoogle Scholar
- Kolby L, Bernhardt P, Ahlman H, Wangberg B, Johanson V, Wigander A, Forssell-Aronsson E, Karlsson S, Ahren B, Stenman G, Nilsson O: A transplantable human carcinoid as model for somatostatin receptor-mediated and amine transporter-mediated radionuclide uptake. Am J Pathol. 2001, 158: 745-755. 10.1016/S0002-9440(10)64017-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Spratt JS, Meyer JS, Spratt JA: Rates of growth of human solid neoplasms: part I. J Surg Oncol. 1995, 60: 137-146. 10.1002/jso.2930600216.View ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.