- Open Access
Perforin-dependent direct cytotoxicity in natural killer cells induces considerable knockdown of spontaneous lung metastases and computer modelling-proven tumor cell dormancy in a HT29 human colon cancer xenograft mouse model
Molecular Cancervolume 13, Article number: 244 (2014)
For long, natural killer (NK) cells have been suspected to play a critical role in suppressing the development of spontaneous metastases in cancer patients. Despite a wide range of studies it remains unclear so far to what extent primary tumor growth together with formation of distant metastases and NK cell activity influence each other.
To precisely investigate the role of NK cells with a perforin-deficiency in cancer growth and metastasis formation, human HT29 colon cancer cells were subcutaneously grafted into pore forming protein and recombination activating gene 2 double knock out (pfp/rag2) mice and in recombination activating gene 2 only knock out (rag2) mice both with black six background. Both mice lack B and T cell functions due to the absence of rag2.
Primary tumors developed in 16/16 in pfp/rag2 and 20/20 rag2 mice. At sacrifice primary tumor weight did not differ significantly. However, tumors grew faster in pfp/rag2 mice (50 days) than in pfp/rag2 mice (70 days). Circulating tumor cells (CTC) in murine blood were nearly three times higher in pfp/rag2 (68 cells/ml) than in rag2 mice (24 cells/ml). Lung metastases occurred frequently in pfp/rag2 mice (13/16) and infrequently in rag2 mice (5/20). The mean number of metastases was 789 in pfp/rag2 mice compared to 210 in rag2 mice. Lung metastases in pfp/rag2 mice consisted of 10–100 tumor cells while those in rag2 mice were generally disseminated tumor cells (DTCs).
Computer modelling showed that perforin-dependent killing of NK cells decelerates the growth of the primary tumour and kills 80% of CTCs. Furthermore, perforin-mediated cytotoxicity hampers the proliferation of the malignant cells in host tissue forcing them to stay dormant for at least 30 days.
The results exactly quantified the effect of perforin-dependent direct cytotoxicity of NK cells on HT29 on primary tumor growth, number of CTCs in the blood and the number of metastases. The largest effects were seen in the number of mice developing spontaneous lung metastases and the mean number of lung metastases. Hence, perforin-mediated cytotoxicity used for direct killing by NK cells is more important than indirect killing by secretion of death-inducing ligands by NK cells.
In the developed world cancer represents the second most common cause of death and more than 90% of these deaths can be attributed to its metastatic spread[1–3]. Therefore control of the metastatic process is the key determinant to an improved survival of cancer patients. The control of cancer cells in the body is complex but one cell population, natural killer (NK) cells, represents the prime position in the defense against cancer cells as it can kill malignant cells without the necessity of prior exposure[4, 5].
NK cells constitute a heterogeneous population of large granular lymphocytes that comprise approximately 10-15% of peripheral blood mononuclear cells (PBMC) in humans. Belonging to the innate immune system NK cells are considered as lymphocytes capable to discover and destroy tumor and virus-infected cells. The name natural killer cells has its origins in their ability to kill virus infected cells spontaneously without the necessity for priming as generally required in the specific immune system[7–9].
Their phenotype is characterized by the surface expression of CD56 (NCAM) and/or CD16 as well as the lack of expression of CD3[8, 10]. In comparison to T cells NK cells do not require the presence of antibodies and surface antigens stringently for their activation. Furthermore they have the unique skill to recognize and respond to the absence of MHC class I antigens on the target cell which is an immunoescape mechanism frequently used by tumor and virus-infected cells to avoid CD8 T-Cell mediated lysis[11–13]. Once activated NK cells have a multitude of mechanisms at their disposal to attack neoplastic cells. One of this mechanisms is the direct release of cytolytic granules containing perforin, granzymes and granulysin by exocytosis to kill target cancer cells[14, 15].
Because of their anti-cancer cell activity a wide range of studies postulates relationships between decreased NK-cell activity and advanced tumor growth as well as the formation of distant metastases but it remains unclear so far to what extent primary tumor growth together with formation of distant metastases and NK cell activity influence each other. Numerous studies indicate that metastatic spread is associated with decreased NK cell activity[16–19] and clinical observations suggest NK-cell activity to be a strong predictor for overall and progression-free survival in carcinoma patients[20, 21]. These clinical findings may be given support by a xenograft mouse model where it could be shown that in small cell lung cancer (SCLC) metastatic spread increased drastically when perforin-dependent killing was disabled in NK cells. However, the two mouse strains used in this experiment had disparate genetic backgrounds. Since the genetic background can significantly impact tumor development[23, 24] the crucial role of NK cells on metastatic spread has not been able to clearly demonstrate as yet.
In order to study the influence of perforin-dependent direct cytotoxicity of NK cells on primary tumor growth and the number of metastatic deposits formed, we used two different mouse strains, namely rag2 and pfp/rag2 mice on the same genetic background. Both mouse strains lack the recombination activating gene (rag2) and therefore share the inability to produce functioning T- and B-cell receptors by V(D)J-recombination. The special characteristic of the pfp/rag2 mice is the additional absence of the pore forming protein (pfp) leading to a suppression of the cytotoxicity of NK cells. The crucial advantage of this xenograft mouse model in comparison to other models investigating the mechanisms of metastatic spread is the fact that these two mouse strains have a completely identical genetic background and are only different in the pfp knock out[25–28]. This direct approach is superior to the use of an antibody against NK cells as the antigen recognized on NK cells by the anti-NK cell antibody is also expressed on tumor cells in general and especially also expressed on HT29 cells. Because of this antibody cross-reactivity we used this knock out model.
Furthermore, the results of the animal experiments will be analyzed via computer simulations. Computer models are valuable tools to quantitatively examine the observed results and to offer various explanations for different scenarios of the metastatic process. Thus, they allow additional valuable perspectives in details of the metastatic process hitherto unexpected.
Engraftment of primary tumors
Overall tumors developed in 100% (n = 16) of the pfp/rag2 mice and in 100% (n = 20) of the rag2 mice. The duration of the experiment ranged from 45 to 106 days in the rag2 mice while it ranged from 35 to 88 days in the pfp/rag2 mice. The mean growth period in the rag2 mice was 69.4 days compared to 49.9 days in the pfp/rag2 mice (Figure 1A, B). The tumor weight range of the rag2 mice varied between 0.08 g and 2.61 g (mean 1.16 g) whereas in the pfp/rag2 mice the weight range varied between 0.20 g and 2.74 g (mean 1.23 g) (Figure 1C). Variances were not statistically significantly different (p = 0.7353).
Circulating tumor cells in blood
The number of circulating tumor cells detected in murine blood by using qRT-PCR varied between 0 and 124 cells per ml (mean 24 cells per ml) in the rag2 mice compared to the pfp/rag2 mice where numbers varied between 1 and 368 cells per ml (mean 68 cells per ml) (Figure 1D). Variances were statistically significantly different (p = 0.0004).
Number of spontaneous lung metastasis
25% of the rag2 mice engrafted with HT29 cells (n = 20) developed spontaneous lung metastases while 81% of pfp/rag2 mice equally engrafted with HT29 cells (n = 16) developed spontaneous lung metastasis detected in HE stained sections (Figure 1E). Within the five rag2 mice that had developed spontaneous lung metastasis the number of metastasis varied between 132 and 264 (n = 5, mean 210). The number of spontaneous lung metastasis in 13 pfp/rag2 mice where metastatic spread of the primary tumors had taken place varied between 67 and 2990 (n = 13, mean 789) (Figure 1F). In total 1048 spontaneous lung metastases were detected in rag2 mice while 10251 spontaneous lung metastases were detectable in pfp/rag2 mice (Figure 1G). Variances were statistically significantly different (p = 0.0001).
The metastasis detected in the two different mouse strains did not only differ in their number but also in their morphology. Typical metastasis detected in rag2 mice consisted of only 1–10 cells at the maximum whereas disseminated tumor cells (DTC) were the most common presentation of metastasis in the rag2 mice. All cells were clearly extravasated. In pfp/rag2 mice typical metastases consisted of approximately 10–100 tumor cells and only sparse DTCs occurred (see Figure 2A, B)
In rag2 mice 2,9% of the tumor cells were mitotic compared to 3.9% of the tumor cells in pfp/rag2 mice (Figure 1H). Variances were statistically significantly different (p = 0.0107).
Computer modelling of cancer spread
In addition to the sole statistical analysis of the data the experimental data was further analyzed utilizing a computer model of the process of the metastatic progression[30, 31]. This computer model allows simulating various different scenarios that can be compared with the experimental data providing a quantitative analysis of the processes underlying metastatic progression in order to generate new hypotheses.
Using the a gompertzian growth function (see equation (1) in the Methods section) and the observed mean values for the tumor weight and the duration of the experiment (1.23 g and 49.9 days for pfp/rag2; 1.16 g and 69.4 days for rag2 mice) a growth rate constant of 0.0462 day-1 was determined for pfp/rag2 mice and 0.0326 day-1 for rag2 mice.
The spread of metastases was modelled via a colonization rate (see equation (2) in the Methods section). Based on the number of metastases in the lung (788 after 49.9 days in pfp/rag2) a colonization coefficient of 5.2e-5 (cell*day)-1 was computed for pfp/rag2. Since the same cell type was used in both mouse types, the primary tumor should exhibit the same spreading behavior in the pfp/rag2 and the rag2 mouse stain. Therefore, the determined colonization coefficient for pfp/rag2 mice was also applied for rag2 mice. To comprise the lesser amount of metastases in the rag2 mouse stain (209 after 69.4 days), the mortality of cycling tumor cells was increased. Different values were tested in the computer model, but only a mortality rate of 80% resulted in the observed number of metastases for rag2 mice, implicating that 80% of those cells that would have been able to found a metastasis in pfp/rag2 mice are killed by NK cells in rag2 mice.
At distant sites NK cells hamper the establishment of DTCs into a metastasis (see Figure 3). Metastases in pfp/rag2 mice have approximately 10–100 cells while in rag2 mice mostly DTCs were observed. The simulation results where no sort of dormancy was applied (black bars in Figure 3) do not fit with the observed morphology of the metastases in pfp/rag2 mice and rag2 mice.
Disseminated tumor cells in rag2 mice remain dormant for at least 30 days, before they are able to successfully proliferate (grey bars in Figure 3B). Applying shorter dormancy time spans, e.g. 21 days (white bars in Figure 3B), several metastases comprising more than 10 cells occur. Simulations with a dormancy of 30 days resulted in mostly disseminated tumor cells and only few multicellular metastases. The different growth rates applied on the metastases (1/3, 1/2 and the same growth rate as the primary tumor) had nearly no effect on the simulation results, as can be seen in Figure 3B.
In pfp/rag2 mice neither the simulation results with no dormancy (black bars in Figure 3A) nor the simulation results with a dormancy of 21 and 30 days (white and grey bars in Figure 3A) fit with the observed morphology. Therefore, a late dormancy was introduced: After the metastasis reaches a random size between 10 and 100 cells it stops growing for a certain time span. An explanation could be that the metastasis has to initiate neoangiogenesis to ensure nutrition and oxygen supply for further growth. The simulation results with an applied late dormancy of 30 days and a metastases growth rate the same as the primary tumors (cross-striped bars in the rightmost graph in Figure 3A), fit quite well with the observed morphology, indicating that the metastases in pfp/rag2 grow at the same rate as the primary tumor.
Summarizing, according to results of the computer simulations, NK cells decelerate the growth of the primary tumor and kill 80% of those cycling tumor cells which could have otherwise established a new metastasis. Furthermore, they hamper the proliferation of the malignant cells in distant tissue forcing them to stay dormant for at least 30 days. Without NK cells it is much easier for malignant cells to establish themselves in the distant tissue and start proliferating. However, they also undergo a dormant phase of at least 30 days after they reach a size between 10 and 100 cells.
The major aim of this study was to investigate to which extent perforin-dependent direct killing of natural killer cells influences both the growth of a primary tumor and the formation of metastases in vivo by using a HT29 xenograft mouse model. Experimental findings that natural killer cells might play a crucial role in the development of primary tumors as well as in spontaneous metastases formation have been described in earlier studies although their exact mode of action has remained elusive so far[9, 32–34]. In an earlier xenograft model using scid and pfp/rag2 mice it has been shown that the number of spontaneous SCLC metastasis increased when perforin-mediated direct cytotoxicity of NK cells was disabled. However the genetic backgrounds of the used mouse strains were not identical as scid mouse in this study had a Balb/cBy background whereas the background of the pfp/rag2 mouse strain was C57BL/6 J.
In order to exclude this strain difference as a possible reason for the difference in metastasis formation and to examine whether natural killer cells and perforin-dependent direct cytotoxicity in particular do play a crucial role in the growth of a primary tumor and its spread to distinct organs we used pfp/rag2 and rag2 mice of the same genetic background to investigate the importance of perforin-dependent killing by NK cells upon primary tumor growth and formation of distant metastasis in vivo.
The average survival of the pfp/rag2 mice was shortened to 49.9 days compared to 69.4 days in the rag2 mice as the mice had reached the termination criteria. This improved survival rate is validated by clinical data showing that colorectal cancer patients with extensive NK cell infiltration had a 5-year cumulative survival rate of 92% compared to only 53.7% of patients with little or moderate NK cell infiltration of their tumors. The mean weight of the primary tumors was nearly identical in both mouse strains with differences in tumor weight being not statistically significant. The primary tumors of pfp/rag2 mice grew to a similar weight as the primary tumors of rag2 mice in only 72% of rag2 mice’s average lifetime. Thus the diminished direct NK cell cytotoxicity apparently enables primary tumors to grow considerably faster even if NK cells were only scarce in colorectal cancer tissue from early stages on as demonstrated in clinical studies. Feulgen immunohistostaining also revealed that primary tumors of pfp/rag2 mice additionally showed statistically significantly more mitotic cells than primary tumors of rag2 mice. Hence direct killing of NK cells considerably influences primary tumor growth in immunodeficient mice.
Furthermore the 2.8 times higher number of circulating tumor cells in murine blood of pfp/rag2 mice compared to rag2 mice indicates that perforin-dependent killing of NK cells is very important for the elimination of circulating tumor cells in murine blood in order to prevent the formation of spontaneous metastases in distant organs. This is in line with the data generated by a computer modelling of cancer spread which revealed that 80% of the circulating tumor cells in the blood of rag2 mice had been eliminated by NK cells before extravasation was possible. These observations are corroborated by clinical findings showing that SCLC patients with high NK cell numbers in the blood had a relatively better prognosis with less distant metastasis than those patients with fewer NK cell numbers.
Unsurprisingly direct cytotoxicity of NK cells does have a major effect on the formation of spontaneous lung metastases as only 25% (n = 20) of the rag2 mice exhibited metastatic spread into the lung at all and the presence of NK cells completely prevented the formation of spontaneous lung metastasis in 75% of the rag2 mice. In contrast 81% (n = 16) of pfp/rag2 mice presented spontaneous lung metastases and only 3 animals, 19%, did not develop lung metastases. The absolute number of spontaneous metastases detected in pfp/rag2 mice is 9.78 times higher than the number of spontaneous metastases of rag2 mice. The almost ten times higher numbers of spontaneous lung metastasis found in pfp/rag2 mice in comparison to 3 times higher numbers of circulating tumor cells in the bloodstream of pfp/rag2 mice indicate that perforin dependent direct cytotoxicity of NK cells reduces the number of malignant cells significantly in both the bloodstream and the connective tissue stroma in the distant organ.
In addition to the number of spontaneous metastases the morphology of the spontaneous lung metastases also differed considerably: metastases in rag2 mice were only 1–10 tumor cells in line with single cell metastases (DTCs) being the most frequent detected metastatic deposits while metastases in pfp/rag2 mice generally consisted of 10–100 tumor cells on average presumably showing that metastasis cell killing also takes place in the metastatic organ. It can be excluded that the observed DTCs in rag2 originate from the initial injection of one million tumor cells, since only five from 20 animals developed DTCs in the lung. Otherwise, all mice should have displayed DTCs in the lung since all mice were injected equally by the same person.
The results of the computer simulations lead to the assumption that direct perforin-mediated killing of NK cells hampers the establishment of metastases in the lung, forcing the malignant cells to stay dormant for at least 30 days. In contrast, when perforin-mediated killing of NK cells is suppressed tumor cells proliferate much easier in the distant tissue. However, after the metastases reaches a size between 10 and 100 cells they also seem to undergo a dormant phase of at least 30 days. This result is unexpected and we do not have an evidence based interpretation why this dormancy occurs, however, this result shows that computer analyses can identify critical points in the growth and development of metastases previously overlooked by less sophisticated statistical analyses.
Considering the fact that the rag2 mice lived 30% longer than the pfp/rag2 mice to develop spontaneous metastasis the difference between the discovered number of metastasis in pfp/rag2 and rag2 mice becomes even more impressive.
Our findings are very similar to the observations made by Sodeur et al. where most notably the total number of metastases of all cell lines used in their experiment were considerably higher in pfp/rag2 mice than in scid mice which leaves hardly any doubt that perforin-dependent direct cytotoxicity of NK cells plays a crucial role in the anti-cancer control of our immune system indicating that the presence of NK cells with the abilitiy of direct perforin-dependent killing is more important in determining the number of metastases than the mouse strain background.
There is clinical uncertainty to what extent tumor growth and NK-cell-activity influence each other. Our results show that perforin-dependent direct cytotoxicity of NK-cells is a major component for combatting tumors and the down regulation of perforin-dependent direct cytotoxicity of NK-cells may lead to amplified primary tumor growth, facilitated earlier dissemination of tumor cells in the bloodstream and the formation of spontaneous distant metastases in particular. Hence the lack of NK-cells enables tumors to spread more easily and expansively.
The enormous differences in the formation of spontaneous metastasis and the dissemination in the bloodstream show that perforin-dependent killing is the major factor of NK cell cytotoxicity in this model as perforin-dependent killing reduced the number of spontaneous metastases by 90% in total.
Altogether our results can confirm the hypothesis that NK cells and in particular their perforin-dependent direct cytotoxicity play a very important role in preventing metastasis formation as they are able to reduce the number of cycling tumor cells in the blood and prevent the formation of distant metastasis to a high degree. Furthermore the mathematical treatment of the data generated in our mouse model has allowed the discovery that NK cells are able to restrict metastatic proliferation by inducing dormancy of malignant cells which is necessary to further explore and may yield an improved treatment of cancer patients in the future.
Animals and experiments
Sixteen eight week old pfp/rag2 mice and twenty eight week old rag2 mice were used for this experiment. The mice were C57BL6 (pfp/rag-2) mice obtained from Taconic, Hudson, NY (# 001177-MM; B6.129S6-Pfp(tm1)Cirk-RAG2(tm1)Fwa, N12). Both strains of mice were kept under pathogen-free conditions in individually ventilated cages (IVC-Rack, Techniplast Germany) and were fed with sterile standard food (ssniff, Soest, Germany) and water ad libitum. The animals were killed when the tumors started to ulcerate or when the tumor weight exceeded 20% of the original mouse body weight at the beginning of the experiment.
The experiment was supervised by the institutional animal welfare officer and approved by the local licensing authority (Behörde für Soziales, Gesundheit, Familie, Verbraucherschutz; Amt für Gesundheit und Verbraucherschutz, Hamburg, Germany, project No. G 09/58).
The human colon adenocarcinoma cell line HT29 was purchased from the European Cell Culture Collection (Porton Down, Wiltshire, UK). Cells were grown in RPMI 1690 – L-Glutamine (GIBCO, Invitrogen Corp., Grand Island, NY) medium supplemented with 10% fetal calf serum, 1% penicillin and streptomycin and cultured in a humidified atmosphere of 37°C and 5% carbon dioxide.
Each mouse was injected subcutaneously with one million viable tumor cells suspended in 200 μl culture medium RPMI 1690 + L-Glutamine (GIBCO, Invitrogen Corp., Grand Island, NY) between the scapulae.
After sacrificing the animals primary tumors were excised and fixed in 4% buffered formaldehyde for 24 h and rinsed with phosphate buffered saline. The tissues were then dehydrated in a series of graded ethanol and embedded in paraffin wax. Five μm thick sections were cut and stained with haematoxylin and eosin (H.E.).
In order to achieve a random distribution of the lung, the lungs were fixed en bloc and were sectioned after fixation into one mm thick lung slices. The slices were placed in warm agar and pressed down with a glass piston. After hardening of the agar the lung slices were processed to paraffin wax.
The agar blocks containing the lung slices were sectioned into five μm thick sections. The total number of sections of each lung was noted. In addition to every 10th section, two series of serial sections (n = 30) out of the middle of the paraffin wax block were preserved for further immunohistological evaluation. Ten of the 10th sections containing the most lung pieces of each wax block were H.E. stained. Metastases were counted in each of the ten stained sections under a microscope (Zeiss, Axioplan2). The number of metastases of each mouse was calculated by using the following term (Mean number of metastasis * total number of sections – 20%), according to an earlier established formula.
DNA extraction from murine blood and cultured human tumor cells
Blood samples were withdrawn from each mouse by puncturing the heart after deep general anesthesia with CO2. Approximately 1 ml of blood could be extracted from each mouse and a total of 200 μl mouse blood was prepared for DNA extraction using the High Pure PCR Template Preparation Kit. 1 × 106 HT29 tumor cells were isolated by using the same kit in order to establish a standard curve for calibration of the DNA.
The extracted DNA was resuspended in 200 μl elution buffer and a sequential dilution series was prepared. To enable comparability the standard solutions that were previously prepared as a bulk preparation were diluted into murine DNA from untreated pfp/rag2 and rag2 mice.
DNA extraction and real-time PCR for detection of circulating tumor cells
DNA extraction of murine blood was performed using the QIAamp DNA Blood Mini Kit and for DNA isolation of cell culture cells the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany) was used, according to manufacturers′ instructions. Real-time polymerase chain reaction (PCR) and melting curve analyses were performed in glass capillaries with the Light Cycler 2.0 System. For the real-time PCR, the LightCycler Fast Start DNA MasterPLUS SYBR-Green I Kit (Roche Diagnostics GmbH, Mannheim, Germany) was used. Two μl of DNA solution was used as a template for the PCR reaction and incubated in a total reaction volume of 10 μl, containing 1x SYBR-Green I Master mix including Taq DNA polymerase, Taq PCR buffer, a dNTP mixture and 1 mmol/l MgCl2-, 10 pmol specific Alu primers. Forward Alu primer (TGG CTC ACG CCT GTA ATC CCA) and reverse Alu primer (GCC ACT ACG CCC GGC TAA TTT) were synthesized by MWG Biotech AG (Ebersberg, Germany). The PCR conditions were initially 10 min. 95°C, followed by 50 cycles of 5 s 95°C, 5 s 67°C and 20 s 72°C (measurement of fluorescence). Melting curve analysis (0 s 95°C, 12 s 65°C and 0 s 95°C) was performed directly after each PCR run. To quantify circulating tumor cells a standard curve with 10 fold dilution of extracted DNA from 1×106 cell culture cells HT29 was established. Control probes were isolated from mouse blood without inoculated tumor cells. Human tumor cells were quantified by real-time polymerase chain reaction (PCR) using established primers specific for the human alu sequences as previous described. For each sample, analyses were performed in triplicates and performed as independent experiments at least twice.
The assessment of cell proliferation in the tumor population was made by a Feulgen stain of primary tumors performed as previously described.
Quantitative methods and statistics
The percentage of cells showing mitotic figures in ten different areas of the tumor, delineated by an eyepiece graticule (310 μm2), was determined by counting a minimum of 500 cells from each animal. The areas of measurement were standardized: one corner of the eyepiece graticule was positioned at the tumor-host interface with an objective lens of magnification 10 and counting of the mitotic figures was carried out at the same site using an objective lens of magnification 400. Only cells that were in easily recognizable meta- and anaphases were counted as mitotic.
All values are presented as mean values. Statistically significant differences between both samples were calculated by a Mann–Whitney U test. Graph Pad Prism 6.0 (Intuitive Software for Science, San Diego, CA, USA) was used for statistical calculations. Differences were considered significant at p <0.05.
Simulating the cancer spread
Computer simulations of cancer spread were performed in order to identify parts of the metastatic progress most influenced by NK cells. A previously developed computer model uses a discrete event simulation approach to analyze the metastatic progression[30, 31].
The main components of the computer model are so called compartments. They describe all parts that contain malignant cells such as the primary tumor, blood and distant metastases. The primary tumor and metastases are modelled as “continuous” compartments utilizing mathematical functions to describe the growth and spreading behavior of the tumor.
The growth of the primary tumor and metastases is modelled by a Gompertzian growth function that describes a sigmoid course:
The function x(t) provides the number of cells in the tumor at the time t. The parameter a represents the growth rate constant while b represents the size of the tumor at its saturated level. The parameter t 0 allows to comprise a start size of the tumor, e.g. if the primary tumor starts as a single cell t 0 will be 0. If the tumor starts as a cluster of cells, due to the injection of tumor cells into the mouse, t 0 can be parameterized to display the size of the cluster. If a start size is given, t 0 is automatically computed via an inverse function by the simulation software.
In this work it was assumed that 104 cells of the injected one million tumor cells survived in the mice to form the primary tumor. Simulations with 103 and 105 cells were equally performed, but since the results do not differ significantly only results for 104 cells are shown.
For the maximum tumor size b a value of 4.5 g was assumed. This value was estimated based on the experimental data. The value of the primary tumor growth rate constant a was computed using the determined mean values for the primary tumor weight and the duration of the experiment. The values are presented in the Results section.
The spread of metastases is described by the colonization rate β(x):
where x is the number of cells in the tumor, m is the colonization constant and α is the fractal dimension of blood vessels infiltrating the tumor which describes how well the tumor is supplied with blood. This value was assumed to be 0.663 which describes a superficial vascularity of the primary tumor. This seems plausible since the primary tumor grows very fast. The colonization constant m was derived from the experimental data and is presented in the Results section.
To save computation time the colonization rate β(x) models only those malignant cells that eventually are able to establish a new metastasis. The so described cells can be killed by NK cells in the model. All disseminated tumor cells that die due to other factors than NK cells, either in the blood stream or in distant tissue are not displayed by β(x).
The blood is modelled as a “discrete” compartment where the behavior of cells is described employing events only. An event describes what happens in a compartment at a specific time, e.g. intravasation, apoptosis or extravasation. The probability with which different types of events occur can be parameterized for each discrete compartment.
Starting with the primary tumor that grows according to the function x(t), the first intravasation event is created conforming to the colonization rate β(x). After processing the intravasation event a new event is created which describes what happens next to the cell in the blood compartment, e.g. get killed by NK cells or extravasate and create a new metastasis. The probability of NK cell induced cell death in the blood compartment was determined based on the experimental data, as described in the Results section. Furthermore, the next intravasation event for the primary tumor is created conforming to the colonization rate β(x).
New metastases also grow according to the function x(t). Three different growth rate constants were considered for the lung metastases: 1/3, 1/2 and the same growth rate constant as the primary tumor. Metastases from metastases were neglected in the simulations, since the metastases are too small to be able to spawn metastases by themselves in the short duration of the experiment.
To determine why the metastases in pfp/rag2 and rag2 display the observed difference in size a potential dormancy of 21 and 30 days with a standard deviation of 7 days of the lung metastases was simulated. For modelling details on dormancy see next section.
A snapshot of the simulated system, containing the actual time, the size of the primary tumor, the number of metastases, the number of cells in all metastases and a size histogram of all metastases, is saved at a parameterizable interval. After the simulation covered a determined time span it will stop. Each configuration is simulated 100 times. Afterwards the mean and the standard deviation are computed.
Simulating dormancy and late dormancy
To simulate dormancy an extended version of the growth function stated in equation (1) was introduced to the computer model:
The parameter t d and t ld allow comprising dormancy and late dormancy[48–50] into the simulation. They display the duration of the dormancy phases. The duration can be parameterized with mean and standard deviation.
When a new metastasis is created which undergoes a dormancy phase, its status is set to “dormant” and the exact duration of the dormancy t d is computed based on the stated mean and standard deviation. As long as the metastasis is in dormant state the simulation software will return a value of 1, whenever the size of the metastasis is enquired. As soon as the computed duration of the dormancy elapsed, the status of the metastasis is reset. The metastasis will now start growing conforming to the growth function stated in equation (3). The parameter t d represents the offset between the creation of the metastasis and the time point when it starts to grow.
When a new metastasis is created which undergoes a late dormancy, the first step is to compute the size at which the metastasis passes into the late dormancy phase. The size is computed based on a parameterizable mean and standard deviation. The time point when the metastasis reaches the computed size is computed via an inverse of the growth function. Until this time point the metastasis will grow unrestricted conforming to the growth function stated in equation (3) with a value of 0 for the parameter t ld . Similar to normal dormancy its status is then changed into dormant state and the duration of the late dormancy phase t ld is computed based on the stated mean and standard deviation. The simulation software will return the computed size whenever the size of the metastasis is enquired until the late dormancy phase elapsed. After the status of the metastasis is reset it will continue growing conforming to the growth function equation (3), but with an updated value for t ld to include the accrued offset.
- NK cells:
Natural killer cells
Pore forming protein and recombination activating gene 2 double knock out
Recombination activating gene 2 only knock out
Circulating tumor cells
Disseminated tumor cells.
Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D: Global cancer statistics. CA: A Cancer J Clin. 2011, 61: 69-90. 10.3322/caac.20107.
Hanahan D, Weinberg RA: Hallmarks of cancer: the next generation. Cell. 2011, 144: 646-674. 10.1016/j.cell.2011.02.013
Watson AJM, Collins PD: Colon cancer: a civilization disorder. Dig Dis. 2011, 29: 222-228. 10.1159/000323926
Trinchieri G: Biology of natural killer cells. Adv Immunol. 1989, 47: 187-376.
Biron CA: Activation and function of natural killer cell responses during viral infections. Curr Opin Immunol. 1997, 9: 24-34. 10.1016/S0952-7915(97)80155-0
Arina A, Murillo O, Dubrot J, Azpilikueta A, Alfaro C, Pérez-Gracia JL, Bendandi M, Palencia B, Hervás-Stubbs S, Melero I: Cellular liaisons of natural killer lymphocytes in immunology and immunotherapy of cancer. Expert Opin Biol Ther. 2007, 7: 599-615. 10.1517/14712522.214.171.1249
Bryceson YT, Long EO: Line of attack: NK cell specificity and integration of signals. Curr Opin Immunol. 2008, 20: 344-352. 10.1016/j.coi.2008.03.005
Farag SS, Caligiuri MA: Human natural killer cell development and biology. Blood Rev. 2006, 20: 123-137. 10.1016/j.blre.2005.10.001
Sanchez-Correa B, Morgado S, Gayoso I, Bergua JM, Casado JG, Arcos MJ, Bengochea ML, Duran E, Solana R, Tarazona R: Human NK cells in acute myeloid leukaemia patients: analysis of NK cell-activating receptors and their ligands. Cancer Immunol Immunother. 2011, 60: 1195-1205. 10.1007/s00262-011-1050-2
Seebach JD, Waneck GL: Natural killer cells in xenotransplantation. Xenotransplantation. 1997, 4: 201-211. 10.1111/j.1399-3089.1997.tb00184.x.
Ljunggren HG, Kärre K: In search of the “missing self”: MHC molecules and NK cell recognition. Immunol Today. 1990, 11: 237-244.
Vivier E, Nunès JA, Vély F: Natural killer cell signaling pathways. Science. 2004, 306: 1517-1519. 10.1126/science.1103478
Vivier E, Raulet DH, Moretta A, Caligiuri MA, Zitvogel L, Lanier LL, Yokoyama WM, Ugolini S: Innate or adaptive immunity? The example of natural killer cells. Science. 2011, 331: 44-49. 10.1126/science.1198687
Cullen SP, Martin SJ: Mechanisms of granule-dependent killing. Cell Death Differ. 2008, 15: 251-262. 10.1038/sj.cdd.4402244
Pardo J, Aguilo JI, Anel A, Martin P, Joeckel L, Borner C, Wallich R, Müllbacher A, Froelich CJ, Simon MM: The biology of cytotoxic cell granule exocytosis pathway: granzymes have evolved to induce cell death and inflammation. Microbes Infect. 2009, 11: 452-459. 10.1016/j.micinf.2009.02.004
Nüssler NC, Stange BJ, Petzold M, Nussler AK, Glanemann M, Guckelberger O: Reduced NK-cell activity in patients with metastatic colon cancer. EXCLI J. 2007, 6: 1-9.
Saijo N, Ozaki A, Beppu Y, Takahashi K, Fujita J, Sasaki Y, Nomori H, Kimata M, Shimizu E, Hoshi A: Analysis of metastatic spread and growth of tumor cells in mice with depressed natural killer activity by anti-asialo GM1 antibody or anticancer agents. J Cancer Res Clin Oncol. 1984, 107: 157-163. 10.1007/BF01032600
Gorelik E, Wiltrout RH, Okumura K, Habu S, Herberman RB: Role of NK cells in the control of metastatic spread and growth of tumor cells in mice. Int J Cancer. 1982, 30: 107-112. 10.1002/ijc.2910300118
Ben-Eliyahu S, Yirmiya R, Liebeskind JC, Taylor AN, Gale RP: Stress increases metastatic spread of a mammary tumor in rats: evidence for mediation by the immune system. Brain Behav Immun. 1991, 5: 193-205. 10.1016/0889-1591(91)90016-4
Liljefors M, Nilsson B, Hjelm Skog A-L, Ragnhammar P, Mellstedt H, Frödin J-E: Natural killer (NK) cell function is a strong prognostic factor in colorectal carcinoma patients treated with the monoclonal antibody 17-1A. Int J Cancer. 2003, 105: 717-723. 10.1002/ijc.11139
Ménard C, Blay J-Y, Borg C, Michiels S, Ghiringhelli F, Robert C, Nonn C, Chaput N, Taïeb J, Delahaye NF, Flament C, Emile J-F, Le Cesne A, Zitvogel L: Natural killer cell IFN-gamma levels predict long-term survival with imatinib mesylate therapy in gastrointestinal stromal tumor-bearing patients. Cancer Res. 2009, 69: 3563-3569. 10.1158/0008-5472.CAN-08-3807
Sodeur S, Ullrich S, Gustke H, Zangemeister-Wittke U, Schumacher U: Increased numbers of spontaneous SCLC metastasis in absence of NK cells after subcutaneous inoculation of different SCLC cell lines into pfp/rag2 double knock out mice. Cancer Lett. 2009, 282: 146-151. 10.1016/j.canlet.2009.03.005
Eccles SA: Basic principles for the study of metastasis using animal models. Metastasis Research Protocols: Volume II. Edited by: Brooks SA, Schumacher U. 2001, 161-171. [Methods in Molecular Medicine™, vol. 58], Totowa, New Jersey: Humana Press,
Puccini J, Dorstyn L, Kumar S: Genetic background and tumour susceptibility in mouse models. Cell Death Differ. 2013, 20: 964- 10.1038/cdd.2013.35
Gellert M: V(D)J recombination: RAG proteins, repair factors, and regulation. Annu Rev Biochem. 2002, 71: 101-132. 10.1146/annurev.biochem.71.090501.150203
Walsh CM, Matloubian M, Liu CC, Ueda R, Kurahara CG, Christensen JL, Huang MT, Young JD, Ahmed R, Clark WR: Immune function in mice lacking the perforin gene. Proc Natl Acad Sci U S A. 1994, 91: 10854-10858. 10.1073/pnas.91.23.10854
Shinkai Y, Rathbun G, Lam KP, Oltz EM, Stewart V, Mendelsohn M, Charron J, Datta M, Young F, Stall AM: RAG-2-deficient mice lack mature lymphocytes owing to inability to initiate V(D)J rearrangement. Cell. 1992, 68: 855-867. 10.1016/0092-8674(92)90029-C
Lowin B, Peitsch MC, Tschopp J: Perforin and granzymes: crucial effector molecules in cytolytic T lymphocyte and natural killer cell-mediated cytotoxicity. Curr Top Microbiol Immunol. 1995, 198: 1-24.
Magnusson KE, Gustafsson M, Holmgren K, Johansson B: Small intestinal differentiation in human colon carcinoma HT29 cells has distinct effects on the lateral diffusion of lipids (ganglioside GM1) and proteins (HLA class 1, HLA class 2, and neoplastic epithelial antigens) in the apical cell membrane. J Cell Physiol. 1990, 143: 381-390. 10.1002/jcp.1041430224
Wedemann G, Bethge A, Haustein V, Schumacher U: Computer simulation of the metastatic progression. Metastasis Research Protocols. 2014, 107-116. [Methods in Molecular Biology, vol. 1070], New York: Springer, 2,
Bethge A, Schumacher U, Wree A, Wedemann G: Are metastases from metastases clinical relevant? Computer modelling of cancer spread in a case of hepatocellular carcinoma. PLoS One. 2012, 7: e35689- 10.1371/journal.pone.0035689
Caligiuri MA: Human natural killer cells. Blood. 2008, 112: 461-469. 10.1182/blood-2007-09-077438
Smyth MJ, Hayakawa Y, Takeda K, Yagita H: New aspects of natural-killer-cell surveillance and therapy of cancer. Nat Rev Cancer. 2002, 2: 850-861. 10.1038/nrc928
Wu J, Lanier LL: Natural killer cells and cancer. Adv Cancer Res. 2003, 90: 127-156.
Coca S, Perez-Piqueras J, Martinez D, Colmenarejo A, Saez MA, Vallejo C, Martos JA, Moreno M: The prognostic significance of intratumoral natural killer cells in patients with colorectal carcinoma. Cancer. 1997, 79: 2320-2328. 10.1002/(SICI)1097-0142(19970615)79:12<2320::AID-CNCR5>3.0.CO;2-P
Halama N, Braun M, Kahlert C, Spille A, Quack C, Rahbari N, Koch M, Weitz J, Kloor M, Zoernig I, Schirmacher P, Brand K, Grabe N, Falk CS: Natural killer cells are scarce in colorectal carcinoma tissue despite high levels of chemokines and cytokines. Clin Cancer Res. 2011, 17: 678-689. 10.1158/1078-0432.CCR-10-2173
Nakamura H, Kawasaki N, Hagiwara M, Saito M, Konaka C, Kato H: Cellular immunologic parameters related to age, gender, and stage in lung cancer patients. Lung Cancer. 2000, 28: 139-145. 10.1016/S0169-5002(99)00133-6
Kajitani K, Tanaka Y, Arihiro K, Kataoka T, Ohdan H: Mechanistic analysis of the antitumor efficacy of human natural killer cells against breast cancer cells. Breast Cancer Res Treat. 2012, 134: 139-155. 10.1007/s10549-011-1944-x
Sjöström-Mattson J, Von Boguslawski K, Bengtsson N-O, Mjaaland I, Salmenkivi K, Blomqvist C: The expression of p53, bcl-2, bax, fas and fasL in the primary tumour and lymph node metastases of breast cancer. Acta Oncol. 2009, 48: 1137-1143. 10.3109/02841860902988688
Lahiry L, Saha B, Chakraborty J, Adhikary A, Mohanty S, Hossain DMS, Banerjee S, Das K, Sa G, Das T: Theaflavins target Fas/caspase-8 and Akt/pBad pathways to induce apoptosis in p53-mutated human breast cancer cells. Carcinogenesis. 2010, 31: 259-268. 10.1093/carcin/bgp240
Trapani JA, Smyth MJ: Functional significance of the perforin/granzyme cell death pathway. Nat Rev Immunol. 2002, 2: 735-747. 10.1038/nri911
Mitchell MJ, Wayne E, Rana K, Schaffer CB, King MR: TRAIL-coated leukocytes that kill cancer cells in the circulation. Proc Natl Acad Sci U S A. 2014, 111: 930-935. 10.1073/pnas.1316312111
Workman P, Aboagye EO, Balkwill F, Balmain A, Bruder G, Chaplin DJ, Double JA, Everitt J, Farningham DAH, Glennie MJ, Kelland LR, Robinson V, Stratford IJ, Tozer GM, Watson S, Wedge SR, Eccles SA, Navaratnam V, Ryder S: Guidelines for the welfare and use of animals in cancer research. Br J Cancer. 2010, 102: 1555-1577. 10.1038/sj.bjc.6605642
Jojovic M, Schumacher U: Quantitative assessment of spontaneous lung metastases of human HT29 colon cancer cells transplanted into SCID mice. Cancer Lett. 2000, 152: 151-156. 10.1016/S0304-3835(99)00443-7
Nehmann N, Wicklein D, Schumacher U, Müller R: Comparison of two techniques for the screening of human tumor cells in mouse blood: quantitative real-time polymerase chain reaction (qRT-PCR) versus laser scanning cytometry (LSC). Acta Histochem. 2010, 112: 489-496. 10.1016/j.acthis.2009.05.004
Mulisch M, Welsch U: Romeis Mikroskopische Technik. 2010, Heidelberg: Spektrum Akademischer Verlag GmbH,
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
Aguirre-Ghiso JA: Models, mechanisms and clinical evidence for cancer dormancy. Nat Rev Cancer. 2007, 7: 834-846. 10.1038/nrc2256
Chiarella P, Bruzzo J, Meiss RP, Ruggiero RA: Concomitant tumor resistance. Cancer Lett. 2012, 324: 133-141. 10.1016/j.canlet.2012.05.021
Wikman H, Vessella R, Pan℡ K: Cancer micrometastasis and tumour dormancy. APMIS. 2008, 116: 754-770. 10.1111/j.1600-0463.2008.01033.x
We thank Christine Knies for her help with the PCR and histology. AB would like to thank Andreas Wree for supervision and support. The work of TB was supported by a fellowship from the Studienstiftung des Deutschen Volkes. The work of AB was supported by the European Social Fund by a scholarship of the State Graduate Funding Program (Landesgraduiertenförderung) of the University of Rostock. The funders had no involvement in study design, collection, analysis and interpretation of data, decision to publish or preparation of the manuscript.
The authors declare that they have no competing interests.
TB collected, analyzed and interpreted the data and drafted the manuscript. NN and US conceived and designed the experiments. They also revised the manuscript critically for important intellectual content. US also participated in the design and analysis of the computer simulations. AB enhanced the computer model, designed and carried out the computer simulations, analyzed the simulation data and drafted the manuscript. GW developed the computer model, participated in the design and analysis of the computer simulations and revised the manuscript critically for important intellectual content. All authors read and approved the final manuscript.
About this article
- Natural killer cell
- Perforin-dependent direct cytotoxicity
- Spontaneous lung metastases xenograft
- HT29 colon cancer cell line
- Computer modelling