Nonlinear Dynamics in Applied Sciences Systems: Advances and Perspectives
View this Special IssueResearch Article  Open Access
On Analytical Methods in Neuroblastoma Detection
Abstract
Nonlinear dynamics of cancer recurrence are known to be governed by several factors as initial tumour size, number of metastatic sites, or quantity of drug resistant cells. The precise extent and location of tumours are very important factors so quantitative and consistent methods of evaluation are needed to assess reponse to patient therapy. Wholebody metaiodobenzylguanidine (mIBG) scintigraphy is used as primary medical image modality to detect neuroblastoma tumours due to its high specificity and sensitivity. However, current oncological guidelines are based on qualitative observerdependent analysis. This fact makes it difficult to compare results of scintigraphies taken at different moments during therapy or at different institutions. In this paper, we review analytical methods used in neuroblastoma detection and propose an observerindependent method to quantitatively analyse a mIBG scintigraphy.
1. Introduction
Nonlinear dynamics of cancer recurrence are governed by a balance between several factors such as initial tumor size, mutation rates, and growth kinetics of drugsensitive and resistance cells [1]. In recent years, there has been an increasing interest for mathematical models in tumour growth, by taking into account the temporal evolution which can be controlled by slight changes in the sequence of tumour images. Based on these images, many theoretical models have been proposed [2–6].
Neuroblastoma is the third most common malignant solid tumor of childhood. It is metastatic or with high risk for recurrence in nearly 50% of cases [7]; therefore, successful planning of individual patient therapy requires precise delineation of the extent of the neuroblastoma.
Anatomic imaging methods, such as Computerized Tomography (CT) and Magnetic Resonance Imaging (MRI), are most useful for the evaluation of the primary tumor mass and nearby involved lymph nodes. Functional imaging radioactive tracers, such as mIBG, EDG, and MDP, are used to assess the extent of disease and to search for distant metastases [8].
In particular, neuroblastoma’s cells take up and metabolize norepinephrine, so its analogue metaiodobenzylguanidine (mIBG) is an ideal specific agent for imaging. mIBG is a gamma emitting radionuclide developed in the early 1980s to visualise tumours [9]. mIBG scintigraphy has been in clinical use for more than 20 years for diagnostic assessment. Current radiological guidelines for image acquisition with mIBG in [10] define that scanning is to be performed 24 hours after its injection. A single or multiple head gamma camera with a large field of view is necessary to acquire planar images. The patient should be placed in the supine position to perform a wholebody imaging. Gamma camera scans the whole patient body counting gamma photons detected at each position. The scan produces a planar image where each pixel’s intensity represents gamma photons detected at this position in patient’s body so that white peaks could be observed at patient’s mIBG avid tissues (e.g., tumour, metastases, etc.). Each wholebody scan includes two images, the front and the rear scan of the patient’s body.
To date, radiological guidelines in [10] define a qualitative interpretation criterion of image which has been defined in terms of visual differences between pathological uptake pattern against physiological (normal) pattern. Following these qualitative guidelines, physicians can perform a good evaluation of a given image in terms of clear positive findings or clear negative results. But they are unable to compare results of a given scintigraphy with other previous ones to evaluate the stage of disease when visual differences between scans are not so evident.
A major drawback to the diagnosis of successive scans is that they display variations in intensity, due to the difficulty of maintaining the same exact conditions of the patient and the gamma camera in every scanning session. This variation in intensity is caused by changes in the accumulation of radioisotopes during each examination, the image quality of the gamma camera, the volume of water intake, the amount of urine excreted, or the response to drug treatment of each patient.
Current oncological studies [11] are focused on finding some observerindependent method allowing of comparison of the results between scintigraphies that may help to stratify patient response to therapy. In this paper, we introduce a simple procedure for quantitative observerindependent measurement of uptake level in a single wholebody mIBG scintigraphy.
This paper is structured as follows. In Section 2, we detail the proposed observerindependent measurement method. Experimental results are provided in Section 3, and conclusions are drawn in Section 4.
2. Materials and Methods
mIBG scintigraphies were acquired in human studies by Nuclear Medicine Unit of Hospital La Fe following guidelines described in [10] after a slow intravenous injection of mIBG (dose: Ci/Kg) and using a scintillation camera set to a photopeak equal to the principal gamma photon emitted, that is, of 159 KeV. A minimum of 250,000 and a maximum of 500,000 photon counts were obtained for each front and rear view.
The group of mIBG scans from each patient were read in a blinded fashion by two trained clinicians without knowledge of the interpretation of the original report made by the nuclear medicine physicians with good correlation.
Each scintigraphic image is a representation of the radioactive distribution inside the body. Each pixel contains a discrete value, which is related to the number of gamma photons detected within the period of acquisition. These discrete values follow a statistical Poisson distribution due to the random nature of radioactive disintegration.
In other words, the number of discrete arrivals over a period of time is not constant but is modelled by Poisson’s equation as follows: where is the probability that the number of detected photons is equal to and is the mean value (which is unknown) of the distribution.
This formula implies in daily nuclear medicine practice that the number of photon counts has to be increased in order to reduce the effect of Poisson noise. This can be achieved by increasing the acquisition time, which increases the patient motion and decreases camera throughput, or increasing the amount of administered radioactive material, which will lead to serious patient’s health problems. Given these limitations, image processing techniques are needed to work without increasing photon counts.
First attempts to quantify uptake in regions of interest (ROI) in scintigraphies were made in 1973. In [12], an uptake measurement of in thyroid studies based on the number of gamma photons detected (counts) was proposed. The main idea was to calculate the ratio between counts in ROI divided by the total body counts.
Hearttomediastinal (H/M) uptake ratio derived from mIBG scintigraphy has been used as a quantitative measurement of uptake with success to monitor response to medical treatment in areas such as cardiac sympathetic imaging [13]. This measure of specific to nonspecific uptake is intrinsically very simple: the ratio of the counts per pixel in two regions of interest (ROIs), one drawn around the heart and the other drawn in the midline upper thorax. This simplicity has contributed to the widespread use of the H/M ratio for comparing different populations. However, there is great variation of H/M ratio among publications from various institutions.
We propose also a very simple approach based on the uptake ratio between neuroblastoma tumours and the rest of tissues to standarize the measurement procedure and to minimize variation in measurements between institutions.
Let us denote by the front (or rear) view of a wholebody mIBG scintigraphy of size pixels. Let be the intensity of pixel at coordinates in image that represents radionuclide uptake by tissue (photon counts) at these coordinates. Since we are interested in measuring the uptake inside the patient’s body, we will only take into account the set .
And let and be the mean and the standard deviation of intensity in , given, respectively, by (2) and (3) as
In order to be able to compare the uptake value at pixel with uptake value at the same patient’s body coordinates in another scintigraphy, we need both values to be transformed into the same normalized scale.
We will transform the intensity level at each pixel into its normalized uptake value with (4) as and we call it the normalized uptake value at coordinates referenced to . Notice that we are only interested in measuring those values that are greater than because they are the relevant uptake points. Uptake at values lower than has a value of .
Once every pixel in scintigraphy is transformed into its normalized uptake value based on its standard deviation from mean , we are able to take measurements that are comparable with the next scintigraphies of the same patient.
When measuring uptake at any body point, image is quantified by areas, but not by pixels. A single pixel may have a high scintigraphy value, but it is not taken into account because noise exists in the original image.
For noise reduction, we need to replace raw counts at region with processed counts having reduced noise. This can be achieved by smoothing counts over neighboring pixels. The smoothed area is a weighted average, in which raw counts at point are replaced by processed counts given by (5) where the weights are chosen to provide a smoothing scheme and are normalized by the sum rule in (6) as
Therefore, when measuring any image region of interest (e.g., liver, tumour, etc.), we will use weights equal to where is the number of pixels in the region. Then the uptake measurement, , of any region of interest is defined as the average of the positive normalized uptakes at ; that is,
3. Results and Discussion
We took measurements in a test set of wholebody mIBG scintigraphies of annoymous patients. This experimental set contains 23 wholebody mIBG scintigraphies belonging to 19 children under 14 labeled patients from P1 to P19 suffering from neuroblastoma taken with different gamma cameras.
We were interested in calculating the uptake present at the whole patient’s body frame which we call and is defined as follows:
At every wholebody scintigraphy in test set, physicians labeled each tumour region in each view (rear and front) and measured with (7) applied to tumour pixels of each . In case a tumour was visible in both views, it had two measurements, the referenced to of the front view and another measurement in the rear view.
We also calculated the wholebody background uptake mean in both views (rear and front) with (8) and .
A comparison of the measurement results of P4 patient’s mIBGs (only the rear view for simplicity) is shown in Figure 2 and Table 1 where we may appreciate that there exists a threshold of reference at that remains close to in all five P4 scintigraphies, mainly in and . Relation between and each (tumour) may be used as an evidence of tumour uptake evolution. Furthermore, values at both views are similar in all the P4 patient’s mIBGs.

In Figure 3 and Table 2, evolution of values in scintigraphies and is consecutive in time of patients P4, P10, P11, and P16. In all four patients, of scintigraphy is close to of previous scintigraphy , specially in patients P10 and P16. may be used as a threshold of reference to compare results.

Comparison between all patients’ results is shown in Figure 4. We show in red squares of rear and front views of the scintigraphies of each patient and in blue diamonds (informed tumours) and . In the whole test set, remains close to , in both rear and front views, that will allow us to compare results between different patients.
4. Conclusions
A simple observerindependent method for measuring the uptake level of mIBG in any body region of children has been presented. This technique enables us to measure the relative amount of uptake with a quantitative and objective unit based on standard deviation of the uptake area. The lesion grey level is used as an uptake measurement evidence. Any region of interest could be measured in any other patient and compared with the estimated normal grey level. The larger the standard deviation value of grey level is at a given region, the greater the possibility of tumour existence is.
Preliminary results in test set show that the relationship between the background uptake mean and the tumours uptake mean may be used as evidence of a patient’s tumours uptake evolution due to the fact that background uptake mean tends to be constant in every analysed mIBG scintigraphy of the same patient at both rear and front views.
Moreover, the results of the wholebody background uptake mean obtained at rear and front views are similar in the whole test set, even though different gamma cameras were used on different patients and years. This indicates that the method may be useful to compare results between different patients and even between different medical institutions.
These quantitative measurements add more descriptive information of tumours that may be included in models of nonlinear dynamics of cancer recurrence to study correlation between uptake and tumour relapse.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors would like to acknowledge the important contribution of doctors A. Rivas and A. Cañete of the Hospital La Fe in Valencia (Spain) for the development of this work.
References
 J. Foo and K. Leder, “Dynamics of cancer recurrence,” Annals of Applied Probability, vol. 23, no. 4, pp. 1437–1468, 2013. View at: Google Scholar
 N. Bellomo and M. Delitala, “From the mathematical kinetic, and stochastic game theory to modelling mutations, onset, progression and immune competition of cancer cells,” Physics of Life Reviews, vol. 5, no. 4, pp. 183–206, 2008. View at: Publisher Site  Google Scholar
 C. Bianca and M. Delitala, “On the modelling of genetic mutations and immune system competition,” Computers & Mathematics with Applications, vol. 61, no. 9, pp. 2362–2375, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. Bianca and M. Pennisi, “Immune systems modelling by topdown and bottomup approaches,” International Mathematical Forum, vol. 7, no. 3, pp. 109–128, 2012. View at: Google Scholar  MathSciNet
 C. Cattani and A. Ciancio, “Separable transition density in the hybrid model for tumorimmune system competition,” Computational and Mathematical Methods in Medicine, vol. 2012, Article ID 610124, 6 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 C. Cattani, A. Ciancio, and A. d'Onofrio, “Metamodeling the learninghiding competition between tumours and the immune system: a kinematic approach,” Mathematical and Computer Modelling, vol. 52, no. 12, pp. 62–69, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. P. Mueller, E. Coppenrath, and T. P. Uger, “Nuclear medicine and multimodality imaging of pediatric neuroblastoma,” Pediatric Radiology, vol. 43, no. 4, pp. 418–427, 2013. View at: Google Scholar
 S. E. Sharp, M. T. Parisi, M. J. Gelfand, G. A. Yanik, and B. L. Shulkin, “Functionalmetabolic imaging of neuroblastoma,” Quarterly Journal of Nuclear Medicine and Molecular Imaging, vol. 57, pp. 6–20, 2013. View at: Google Scholar
 M. Nakajo, B. Shapiro, and J. Copp, “The normal and abnormal distribution of the adrenomedullary imaging agent m[i131]iodobenzylguanidine (i131 mibg) in man: evaluation by scintigraphy,” Journal of Nuclear Medicine, vol. 24, no. 8, pp. 672–682, 1983. View at: Google Scholar
 E. Bombardieri, F. Giammarile, C. Aktolun et al., “131i/123iMetaiodobenzylguanidine (mibg) scintigraphy: procedure guidelines for tumour imaging,” European Journal of Nuclear Medicine and Molecular Imaging, vol. 37, no. 12, pp. 2436–2446, 2010. View at: Publisher Site  Google Scholar
 K. K. Matthay, B. Shulkin, R. Ladenstein et al., “Criteria for evaluation of disease extent by 123imetaiodobenzylguanidine scans in neuroblastoma: a report for the International neuroblastoma risk group (inrg) task force,” British Journal of Cancer, vol. 102, no. 9, pp. 1319–1326, 2010. View at: Publisher Site  Google Scholar
 M. N. Maisey, T. K. Natarajan, P. J. Hurley, and H. N. Wagner Jr., “Validation of a rapid computerized method of measuring ^{99}mTc pertechnetate uptake for routine assessment of thyroid structure and function,” Journal of Clinical Endocrinology and Metabolism, vol. 36, no. 2, pp. 317–322, 1973. View at: Publisher Site  Google Scholar
 W. Chen, Q. Cao, and V. Dilsizian, “Variation of hearttomediastinal ratio in^{123}ImiBG cardiac sympathetic imaging: its affecting factors and potential corrections,” Current Cardiology Reports, vol. 13, no. 2, pp. 132–137, 2011. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 R. MartínezDíaz et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.