- Software
- Open Access

# Instantaneous mutation rate in cancer initiation and progression

- Shuhao Sun
^{1}, - Fima Klebaner
^{1}, - Xinan Zhang
^{2}and - Tianhai Tian
^{1}Email author

**12 (Suppl 6)**:110

https://doi.org/10.1186/s12918-018-0629-z

© The Author(s) 2018

**Published:**22 November 2018

## Abstract

### Background

Cancer is one of the leading causes for the morbidity and mortality worldwide. Although substantial studies have been conducted theoretically and experimentally in recent years, it is still a challenge to explore the mechanisms of cancer initiation and progression. The investigation for these problems is very important for the diagnosis of cancer diseases and development of treatment schemes.

### Results

To accurately describe the process of cancer initiation, we propose a new concept of gene initial mutation rate based on our recently designed mathematical model using the non-constant mutation rate. Unlike the widely-used average gene mutation rate that depends on the number of mutations, the gene initial mutation rate can be used to describe the initiation process of a single patient. In addition, we propose the instantaneous tumour doubling time that is a continuous function of time based on the non-constant mutation rate. Our proposed concepts are supported by the clinic data of seven patients with advanced pancreatic cancer. The regression results suggest that, compared with the average mutation rate, the estimated initial mutation rate has a larger value of correlation coefficient with the patient survival time. We also provide the estimated tumour size of these seven patients over time.

### Conclusions

The proposed concepts can be used to describe the cancer initiation and progression for different patients more accurately. Since a quantitative understanding of cancer progression is important for clinical treatment, our proposed model and calculated results may provide insights into the development of treatment schemes and also have other clinic implications.

## Keywords

- Cancer
- Mathematical model
- Initial mutation rate
- Cancer cell doubling time

## Background

Pancreatic cancer is one of the most aggressive malignancies in humans, with a five-year relative survival rate of only 8% patients [1, 2]. This disease occurs when the damaged cells grow in an uncontrolled manner. The most common treatment options for pancreatic cancer patients include surgery, endoscopic treatment, chemotherapy and radiation therapy [3]. The design of a treatment plan is based on a number of factors, including the severity and spread of the tumour, as well as the patient’s health conditions and age [4, 5]. Due to the location of the pancreas which is deep within the abdomen, this type of cancer is difficult to diagnose and is often found at an advanced stage. The observed aggressive malignancy of this disease may be due to the factors such as the delay in diagnosis or early metastatic dissemination [6–8]. Therefore, it is important to conduct experimental and theoretical studies for the initiation of this disease and the dynamics of cancer progression, which may be helpful to diagnose this disease at an earlier stage [9].

It is widely recognised that cancer diseases are initiated from gene mutations that increase the fitness of cancer cells over that of the surrounding normal cells [10–12]. The recent advances in high-throughout technologies and systems biology approaches have provided huge amount of data showing the mutation heterogeneity in cancer cells [13–16]. A number of theoretical studies have been proposed to investigate the initiation and progression of cancer diseases [17–19]. Mathematical modelling has been used as a powerful tool to elucidate mechanisms in cancer initiation and progression [20–23]. These models have provided quantitative predictions that may be validated by experimental or clinical studies. For example, regarding the mechanisms of metastacis, a modelling study has suggested that more than 10 years might be needed to generate the first parental, non-metastatic founder cell from the first gene mutation [18]. In addition, more than five years may be needed for cancer cells to achieve the metastatic ability.

A key assumption in the mathematical models for cancer study is that the gene mutation rate is constant, which is defined as the number of mutations over a fixed time period. The value of mutation rate is obtained using the ratio of mutation number to the number of cells in the population [17, 18]. However, experimental studies suggests that the mutation rate of a cell may be dependent on the number of gene mutations inside the cell, though it is still in the debate regarding the effects of driver mutation and passenger mutation in gene mutation and cell growth [24, 25]. Growing evidence suggests that gene mutations can be deleterious to cancer cells and play an important role in both cancer progression and clinical treatment. In the process of tumorigenesis and the metastasis of pancreatic cancer, at least 3–7 driver mutations have been identified in the clinical and experimental studies. Experimental studies suggest that the telomeric shortening and mutations in Kras gene are among the earliest and most pervasive alterations [26–28].

In a previous study, we have proposed a model using the non-constant gene mutation rate [29]. However, this initial study did not provide the detailed method for determining the variations of gene mutations in different patients. In this work, based on the non-constant mutation rate, we further propose the concept of the gene initial mutation rate and instantaneous tumour double time. Using the data of seven patients, we calculate the initial mutation rate of each patient. Then we calculate the tumour size using the instantaneous tumour double time.

## Implementation and results

### Mathematical model

*N*is the number of mutations,

*μ*(

*t*) is the mutation rate at time

*t*,

*a*is the mutation rate of normal cells (namely the cells without any gene mutation) and

*b*is a constant. The solution of Eq. (1) with regarding to

*N*is

*c*is an arbitrary constant. Note that this solution is valid if

*b*≠0. Then the mutation rate can be written as

*b*≠0. Combined with the case of

*b*=0 (namely

*μ*=

*a*), the mutation rate is denoted as

Here mutation rate *μ*(0) is termed as the initial mutation rate which is the mutation rate for cells without any gene mutation. If the value of *b* is zero, the mutation rate is a constant which has been widely used in the literature.

*p*

_{j}(

*t*) be the fraction of cancer cells with

*j*mutations at time

*t*. In addition,

*t*

_{j}is the time point when the first cancer cell with exact

*j*mutations appears. Here we assume that

*t*

_{0}=0. When using the continuous mutation rate

*μ*(

*t*), the system is modelled by

*j*=1,...,

*k*−1. Figure 1 provides simulations of model (3) based on either a constant mutation rate (namely

*μ*(

*t*)=

*c*

*o*

*n*

*s*

*t*) or a rate function of time

*t*(2). Here we considered a cancer system with a maximal number of 8 mutations. Figure 1 suggest that the difference between simulations obtained by the two types of mutation rates is small if the number of mutations

*k*is small (see Fig. 1a). However, Fig. 1b shows that the difference may be large when the mutation number is large. In this simulation,

*a*=0.000001,

*b*=0.00003, and the initial condition is

*p*

_{0}(0)=1 and

*p*

_{i}(0)=0 (

*i*=1,…,8).

*t*

_{j−1},

*t*

_{j}] between the two consecutive mutations. Using this piece-wise mutation rate, we have the following model [29]

where *j*=1,...,*k*−1.

*t*

_{j}, at which the first

*j*-mutated cell occurs, satisfies

where *λ*_{j}=*λ*(*t*_{j}).

The clinic data may provide the gene mutation number *N* and the time point when the *N*-th mutation occurs only. That is why the average mutation rate has been widely used in studies. However, it is not obvious to determine the initial mutation rate and parameter *b* in the non-constant mutation rate model based on this limited information. The major contribution of this work is to derive the relationship between the initial mutation rate *μ*(0), parameter *b*, mutation number *N* and average mutation rate.

### Determination of non-constant gene mutation rate

*b*and gene initial mutation rate

*μ*(0) based on the gene mutation number

*N*and average mutation rate \(\overline {\mu }\). We first consider the average mutation rate, defined by

*μ*

_{0}is given by

*j*=1,...,

*N*, the analytical solution for

*b*may not exist. Thus our goal is to find an approximation of

*b*with good accuracy.

*j*=1,...,

*N*. Substituting the expression (5) of

*t*

_{N}into the above equation, we have that

*b*is very small, we have that

*b*with good accuracy, we assume that

*j*=1,2,…,

*N*−1, and

*λ*

_{N}=

*N*. Thus from Eq. (7), we again obtain Eq. (6), which implies that assumption (8) is reasonable. On the other hand, using the notation

###
**Theorem 1**

The relationship between the average mutation rate \(\overline {\mu }\), number of mutations *N* and parameter *b* is given by

###
*Proof*

The average mutation rate satisfies □

Thus this theorem is proved by using Eqs. (6) and (10).

Note that, when the values of \(\overline {\mu }_{N}\) and *N* are given, we can find the value of *b* by solving the nonlinear Eq. (9). We use MAPLE to solve this equation and obtain the value of *b*.

### Initial mutation rate

*b*for the seven patients are also given in Table 1. In addition, we calculate the initial mutation rate using the total number of mutations, parameter

*b*and average mutation rate, given by

Estimated average mutation rate, initial mutation rate and the value of parameter *b* based on the clinic data of seven patients from [18] (Survival from diagnosis: month)

Patients | Survival from diagnosis | Clone time diagnosis | Mutations | Mutation rate | Initial rate | b |
---|---|---|---|---|---|---|

Pa01C | 6 | 9.8 | 49 | 0.0192 | 0.01917 | 0.00001 |

Pa02C | 8 | 9.4 | 35 | 0.019 | 0.0189 | 0.000018 |

Pa03C | 1 | 2.4 | 28 | 0.0223 | 0.0201 | 0.000030 |

Pa04C | 7 | 7.9 | 34 | 0.0188 | 0.0185 | 0.000019 |

Pa05C | 10 | 4.3 | 28 | 0.0194 | 0.0189 | 0.000029 |

Pa07C | 3 | 3.1 | 50 | 0.0198 | 0.0198 | 0.00001 |

Pa08C | 15 | 10.6 | 35 | 0.0193 | 0.0190 | 0.000018 |

*S*] is the survival time of patients from diagnosis. We use the least square regression method to find the values of coefficients of

*a*and

*b*which are −0.0019 and 0.0194, respectively. The values in Fig. 2 clearly show the negative correlation between the initial mutation rate and survival time of patients.

*c*and

*d*which are −0.0019 and 0.0211, respectively. Although Fig. 2 suggests that the average mutation rate also has a consistent negative correlation with the survival time, the scaled mean-square error, defined by

*μ*

_{i}is the estimated mutation rate in Table 1 for patient

*i*and \(\mu _{i}^{*}\) is the predicted mutation rate by using either the regression (12) or (13). The error for the predicted initial mutation rate and that for average mutation rate is 0.1007 and 0.1197, respectively. This result suggests that the initial mutation rate provides a better indicator for patient survivor time than the average mutation rate. We have also find the correlation regression relationship between the survival time and the age of patient, as well as the relationship between the survival time and gene mutation number. Numerical results suggest that these variables are not as good as the gene mutation rate for the indicator of survival time (Results not shown).

### Instantaneous tumour double time

Pancreatic cancer has an extremely poor prognosis. Factors that appear to be important in predicting long-term survival following resection include clear surgical margins, small tumour size (2 cm), negative lymph nodes, and reduced perioperative morbidity [7, 30, 31]. It is well known that metastasis accounts for 90% of cancer deaths [32]. The challenge is whether we can detect the tumour during the stage T1 (namely the time between tumour initiation and the birth of the cell giving rise to the parental clone), or even after stage T1 but before seeding of metastases. Advanced imaging methods, as well as other test methods to detect cancer-specific proteins, transcripts, or genes, tumor markers, offer hopes for such non-invasive early detection.

To address this issue, we study the time required for cancer cells successfully to leave the primary tumour based on our non-constant mutation rate. We first refine the notion of tumour doubling time (DT) and propose the concept of instantaneous tumour double time, which is defined as the average tumour double time over a very short period of time. Although the concept of tumour doubling time is widely used for the quantification of tumour growth rate, it involves a number of factors for the growth of cancer cells, such as the tumour type, growth stage, presence of symptoms, and the patient’s lifestyle. Thus it is not easy to determine and analyse the instantaneous tumour doubling time accurately. Based on the research results in [5], Yachida et al. considered the doubling time curve as a piecewise linear function [18]. It assumes that the cell division time is 2.3 days for pancreatic cancer when the tumour size is less than 1 mm, and after that (namely the tumour size is larger than 1 mm) the cell division time is 56 days [18]. The key factor behind this assumption is angiogenesis. Here we consider a continuous curve which fits the Amikura-Yachida curve with good accuracy.

*t*,

where *μ*_{0} and *b* are the initial mutation number and mutation parameter of a particular patient, \(\overline {\mu }_{0}\) and \(\overline {b}\) are the average parameters based on the corresponding values of all patients, and *a* and *K* are parameters to match the data in [18]. Thus it is clear the proposed method can make more accurate prediction if the number of patients is larger.

*a*and

*K*, we remove the initial mutation rate and parameter

*b*, the average value of the doubling time is given by

*a*=40.74 and

*K*=62.04. Figure 3a gives the Amikura-Yachida curve and our approximated curve (14). When

*t*<53 days, the difference between these two curves is small. However, when

*t*>53 days, this approximated curve provides a more reasonable prediction regarding the doubling time. Using the estimated values of

*a*and

*K*together with the gene mutation rate of each patient, Fig. 3b gives the estimated doubling time of the seven patients.

### Tomour growth curves in metastasis

Estimated tumour size for seven patients. Tumour-size (L) (M) (H) is based on the lower bound, average time and higher bound of the clone time

Patients | Pa01c | Pa02c | Pa03c | Pa04c | Pa05x | Pa07c | Pa08c |
---|---|---|---|---|---|---|---|

Clone time(y) | 4.05-5.94 | 3.3-4.84 | 3.9-5.72 | 3.45-5.06 | 3.3-4.8 | 3.7-5.4 | 3.1-4.6 |

Tumor-size(L) | 4.55 | 4.12 | 9.08 | 3.52 | 4.63 | 5.95 | 4.20 |

Tumor-size(M) | 4.73 | 4.46 | 10.15 | 3.80 | 5.15 | 6.25 | 4.58 |

Tumor-size(H) | 4.82 | 4.71 | 11.09 | 4.01 | 5.59 | 6.42 | 4.86 |

Time for 2cm | 0.22 | 0.28 | 0.17 | 0.325 | 0.275 | 0.19 | 0.29 |

Numerical results in Table 2 show that the lowest value of tumour size for these seven patients all are greater than 2cm. Our calculation presents a theoretic support to the empirical tumour size (<2 cm) hypothesis in [5]. In addition, Japan Pancreatic Cancer Registry [4] reported 11,317 patients with carcinoma of the pancreas during the past decade and 3743 patients underwent pancreatectomy. The 5-yr survival rate of all patients undergoing pancreatectomy, including those with malignant islet cell tumour and cystadenocarcinoma, was 16.6%. The 5-yr survival rate of patients with carcinoma of the pancreas of 2 cm or less, excluding malignant islet cell tumour and cystadenocarcinoma, was 36.2% after pancreatectomy. This registry shows that survival rate is on the increase in patients with carcinoma of the pancreas after pancreatectomy, especially in resettable cases with tumour of 2 cm or less. Further improvement in survival rate is expected in carcinoma of pancreas. Our quantitative analysis provides a theoretic support to the above conclusion.

## Discussion and conclusion

In this work, we have proposed new concepts for the initial mutation rate and instantaneous tumour doubling time based on our recently designed mathematical model using non-constant mutation rates. These concepts are aimed at replacing the widely used definitions of the average mutation rate and average doubling time. The proposed initial mutation rate is independent of the mutation number which determines the average mutation rate. Our regression results have suggested that, compared with the average mutation rate, the estimated initial mutation rate has a larger value of correlation coefficient with the patient survival time. In addition, our instantaneous tumour doubling time is a continuous function of time and considers the effect of initial mutation rate and non-constant mutation rate. Thus, compared with the existing model in which the doubling time is a piece-wise linear function, our continuous model may be able to provide more reasonable estimate of tumour growth process. Since a quantitative understanding of cancer progression is important for clinical treatment, our proposed model and calculated results may provide insights into the dynamics of cancer metastasis and hence have clinic implications.

Since the gene mutation rate is usually very small, gene mutations are normally observed in experiments over a long time period. The estimated values of the mutation rate constant in literature all are based on the number of gene mutations over the given time span. Thus, the average mutation rate strongly relies on the observed number of gene mutations. The calculated mutation rate is highly stochastic due to the very small value of mutation number. Our proposed continuous model for the mutation rate considered not only the total mutation number but also the average dynamics of gene mutation which is realized by the model of cancer cell progression. The major contribution of this work is to derive an analytic expression for the instantaneous mutation rate based on the total number of mutation.

## Availability and requirements

**Project name**: Not applicable.

**Project home page**: Not applicable.

**Operating system(s)**: Not applicable.

**Programming language**: MATLAB.

**Other requirements**: Not applicable.

**License**: Not applicable.

**Any restrictions to use by non-academics**: Not applicable.

## Declarations

### Funding

This work is supported by the National Natural Science Foundation of China (No11371161) and the self-determined research funds of CCNU from the colleges basic research and operation of MOE (Grant No.CCNU16JCZX10). The publication costs were funded by Monash University.

### Availability of data and materials

The dataset was derived from the published literature that can be found in the corresponding references.

### About this supplement

This article has been published as part of *BMC Systems Biology Volume 12 Supplement 6, 2018: Selected articles from the IEEE BIBM International Conference on Bioinformatics & Biomedicine (BIBM) 2017: systems biology.* The full contents of the supplement are available online at https://bmcsystbiol.biomedcentral.com/articles/supplements/volume-12-supplement-6.

### Authors’ contributions

Conceived and designed the experiments: TT. Performed the experiments: SS, TT. Analyzed the data: SS. Wrote the paper: SS, FK, XZ, TT. Revising the draft: SS, FK, XZ, TT. All authors read and approved the final manuscript.

### Ethics approval and consent to participate

Not applicable.

### Consent for publication

Not applicable.

### Competing interests

The authors declare that they have no competing interests.

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

## Authors’ Affiliations

## References

- Kamisawa T, Wood LD, Itoi T, Takaori K. Pancreatic cancer. Lancet. 2016; 388(10039):73–85.View ArticleGoogle Scholar
- Jemal A, Siegel R, Xu J, Ward E. Cancer statistics, 2010. CA Cancer J Clin. 2010; 60(5):277–300.View ArticleGoogle Scholar
- Deplanque G, Demartines N. Pancreatic cancer: are more chemotherapy and surgery needed?. Lancet. 2017; 389(10073):985–6.View ArticleGoogle Scholar
- Committee of Pancreatic Cancer Registry (Chairman: Saitoh Y). Report of pancreatic cancer registry, 1981-1990. Jpn Pancreas Soc Kobe. 1991 (Japanese).Google Scholar
- Amikura K, Kobari M, Matsuno S. The time of occurrence of liver metastasis in carcinoma of the pancreas. Int J Pancreatol. 1995; 17(2):139–46.PubMedGoogle Scholar
- Haeno H, Gonen M, Davis MB, Herman JM, Iacobuzio-Donahue CA, Michor F. Computational modeling of pancreatic cancer reveals kinetics of metastasis suggesting optimum treatment strategies. Cell. 2012; 148(1-2):362–75.View ArticleGoogle Scholar
- Richter A, Niedergethmann M, Sturm JW, Lorenz D, Post S, Trede M. Long-term results of partial pancreaticoduodenectomy for ductal adenocarcinoma of the pancreatic head: 25-year experience. World J Surg. 2003; 27(3):324–9.View ArticleGoogle Scholar
- Gundara JS, Gill AJ, Samra JS. Efficacy of primary tumour volume as a predictor of survival compared with size alone in pancreatic ductal adenocarcinoma. Oncol Lett. 2015; 10(2):744–8.View ArticleGoogle Scholar
- Siegel RL, Miller KD, Jemal A. Cancer statistics, 2016. CA Cancer J Clin. 2016; 66(1):7–30.View ArticleGoogle Scholar
- Vogelstein B, Kinzler KW. Cancer genes and the pathways they control. Nat Med. 2004; 10(8):789.View ArticleGoogle Scholar
- Greenman C, Stephens P, Smith R, Dalgliesh GL, Hunter C, Bignell G, Davies H, Teague J, Butler A, Stevens C, et al. Patterns of somatic mutation in human cancer genomes. Nature. 2007; 446(7132):153.View ArticleGoogle Scholar
- Tian T, Olson S, Whitacre JM, Harding A. The origins of cancer robustness and evolvability. Integr Biol. 2011; 3(1):17–30.View ArticleGoogle Scholar
- Parsons DW, Jones S, Zhang X, Lin JC-H, Leary RJ, Angenendt P, Mankoo P, Carter H, Siu I-M, Gallia GL, et al. An integrated genomic analysis of human glioblastoma multiforme. Science. 2008; 321(5897):1807–12.View ArticleGoogle Scholar
- Sjöblom T, Jones S, Wood LD, Parsons DW, Lin J, Barber TD, Mandelker D, Leary RJ, Ptak J, Silliman N, et al. The consensus coding sequences of human breast and colorectal cancers. science. 2006; 314(5797):268–74.View ArticleGoogle Scholar
- Leary RJ, Lin JC, Cummins J, Boca S, Wood LD, Parsons DW, Jones S, Sjöblom T, Park B-H, Parsons R, et al. Integrated analysis of homozygous deletions, focal amplifications, and sequence alterations in breast and colorectal cancers. Proc Natl Acad Sci. 2008; 105(42):16224–9.View ArticleGoogle Scholar
- Bozic I, Antal T, Ohtsuki H, Carter H, Kim D, Chen S, Karchin R, Kinzler KW, Vogelstein B, Nowak MA. Accumulation of driver and passenger mutations during tumor progression. Proc Natl Acad Sci. 2010; 107(43):18545–50.View ArticleGoogle Scholar
- Jones S, Chen W-d, Parmigiani G, Diehl F, Beerenwinkel N, Antal T, Traulsen A, Nowak MA, Siegel C, Velculescu VE, et al. Comparative lesion sequencing provides insights into tumor evolution. Proc Natl Acad Sci. 2008; 105(11):4283–8.View ArticleGoogle Scholar
- Yachida S, Jones S, Bozic I, Antal T, Leary R, Fu B, Kamiyama M, Hruban RH, Eshleman JR, Nowak MA, et al. Distant metastasis occurs late during the genetic evolution of pancreatic cancer. Nature. 2010; 467(7319):1114.View ArticleGoogle Scholar
- Sun S, Klebaner F, Tian T. A new model of time scheme for progression of colorectal cancer. BMC Syst Biol. 2014; 8(3):2.View ArticleGoogle Scholar
- Altrock PM, Liu LL, Michor F. The mathematics of cancer: integrating quantitative models. Nat Rev Cancer. 2015; 15(12):730.View ArticleGoogle Scholar
- Louzoun Y, Xue C, Lesinski GB, Friedman A. A mathematical model for pancreatic cancer growth and treatments. J Theor Biol. 2014; 351:74–82.View ArticleGoogle Scholar
- Li L, Tian T, Zhang X. The impact of radiation on the development of lung cancer. J Theor Biol. 2017; 428:147–52.View ArticleGoogle Scholar
- Li L, Tian T, Zhang X. Mutation mechanisms of human breast cancer. J Comput Biol. 2018; 25(4):396–404.View ArticleGoogle Scholar
- Loeb LA. Human cancers express mutator phenotypes: origin, consequences and targeting. Nat Rev Cancer. 2011; 11(6):450.View ArticleGoogle Scholar
- Duesberg P, Stindl R, Hehlmann R. Explaining the high mutation rates of cancer cells to drug and multidrug resistance by chromosome reassortments that are catalyzed by aneuploidy. Proc Natl Acad Sci. 2000; 97(26):14295–300.View ArticleGoogle Scholar
- Hruban RH, Goggins M, Parsons J, Kern SE. Progression model for pancreatic cancer. Clin Cancer Res. 2000; 6(8):2969–72.PubMedGoogle Scholar
- Iacobuzio-Donahue CA. Genetic evolution of pancreatic cancer: lessons learnt from the pancreatic cancer genome sequencing project. Gut. 2012; 61(7):1085–94.View ArticleGoogle Scholar
- van Heek NT, Meeker AK, Kern SE, Yeo CJ, Lillemoe KD, Cameron JL, Offerhaus GJA, Hicks JL, Wilentz RE, Goggins MG, et al. Telomere shortening is nearly universal in pancreatic intraepithelial neoplasia. Am J Path. 2002; 161(5):1541–7.View ArticleGoogle Scholar
- Sun S, Klebaner F, Tian T. Mathematical model for pancreatic cancer progression using non-constant gene mutation rate. In: 2017 IEEE International Conference on Bioinformatics and Biomedicine (BIBM). Piscataway: IEEE: 2017. p. 138–141.Google Scholar
- Paik KY, Choi SH, Heo JS, Choi DW. Analysis of liver metastasis after resection for pancreatic ductal adenocarcinoma. World J Gastrointest Oncol. 2012; 4(5):109.View ArticleGoogle Scholar
- Campbell PJ, Yachida S, Mudie LJ, Stephens PJ, Pleasance ED, Stebbings LA, Morsberger LA, Latimer C, McLaren S, Lin M-L, et al. The patterns and dynamics of genomic instability in metastatic pancreatic cancer. Nature. 2010; 467(7319):1109.View ArticleGoogle Scholar
- Nguyen DX, Massagué J. Genetic determinants of cancer metastasis. Nat Rev Genet. 2007; 8(5):341.View ArticleGoogle Scholar