Inference of Tumor Phylogenies from Genomic Assays on Heterogeneous Samples

Tumorigenesis can in principle result from many combinations of mutations, but only a few roughly equivalent sequences of mutations, or “progression pathways,” seem to account for most human tumors. Phylogenetics provides a promising way to identify common progression pathways and markers of those pathways. This approach, however, can be confounded by the high heterogeneity within and between tumors, which makes it difficult to identify conserved progression stages or organize them into robust progression pathways. To tackle this problem, we previously developed methods for inferring progression stages from heterogeneous tumor profiles through computational unmixing. In this paper, we develop a novel pipeline for building trees of tumor evolution from the unmixed tumor data. The pipeline implements a statistical approach for identifying robust progression markers from unmixed tumor data and calling those markers in inferred cell states. The result is a set of phylogenetic characters and their assignments in progression states to which we apply maximum parsimony phylogenetic inference to infer tumor progression pathways. We demonstrate the full pipeline on simulated and real comparative genomic hybridization (CGH) data, validating its effectiveness and making novel predictions of major progression pathways and ancestral cell states in breast cancers

The analysis of relapse-free survival curves: implications for evaluating intensive systemic adjuvant treatment regimens for breast cancer

Results of adjuvant dose intensification studies in patients with localised breast cancer have raised questions regarding the clinical usefulness of this treatment strategy. Here, we develop and fit a natural history model for the time to clinical tumour recurrence as a function of the number of involved lymph nodes, and derive plausible predictions of the effects of dose intensification under various conditions. The time to tumour recurrence is assumed to depend on the residual postoperative micrometastatic burden of tumour, the fractional reduction of residual tumour burden (RTB) by treatment, and the rate of regrowth of the RTB to a clinically detectable size. It is assumed that a proportion of micrometastatic tumours are unresponsive to adjuvant chemotherapy even at maximal dose intensity. Data fitted included the San Antonio Cancer Institute (SACI) database of untreated patients, and CALGB #9082, a study comparing a highly intensive and moderately intensity adjuvant regimen in patients with 10+ positive axillary nodes. The proportion of tumours unresponsive to maximally intensive adjuvant treatment is estimated to be 48% (29–67%). The estimated log kill for intermediate-dose therapy from CALGB #9082 was 6.5 logs, compared with 9 logs or greater for high-dose therapy. The model is consistent with a modest but nonnegligible advantage of dose intensification compared with standard therapies in patients with sensitive tumours who have 10+ positive axillary nodes, and suggests that much of this clinical benefit could be achieved using intermediate levels of treatment intensification. The model further suggests that, in patients with fewer than 10 involved axillary nodes, any advantage of treatment intensification over standard therapy would be much reduced, because in patients with smaller tumour burdens of sensitive tumour, a larger proportion of cures achievable with intensified therapy could be achieved as well with standard therapy

Applying unmixing to gene expression data for tumor phylogeny inference

<p>Abstract</p> <p>Background</p> <p>While in principle a seemingly infinite variety of combinations of mutations could result in tumor development, in practice it appears that most human cancers fall into a relatively small number of "sub-types," each characterized a roughly equivalent sequence of mutations by which it progresses in different patients. There is currently great interest in identifying the common sub-types and applying them to the development of diagnostics or therapeutics. Phylogenetic methods have shown great promise for inferring common patterns of tumor progression, but suffer from limits of the technologies available for assaying differences between and within tumors. One approach to tumor phylogenetics uses differences between single cells within tumors, gaining valuable information about intra-tumor heterogeneity but allowing only a few markers per cell. An alternative approach uses tissue-wide measures of whole tumors to provide a detailed picture of averaged tumor state but at the cost of losing information about intra-tumor heterogeneity.</p> <p>Results</p> <p>The present work applies "unmixing" methods, which separate complex data sets into combinations of simpler components, to attempt to gain advantages of both tissue-wide and single-cell approaches to cancer phylogenetics. We develop an unmixing method to infer recurring cell states from microarray measurements of tumor populations and use the inferred mixtures of states in individual tumors to identify possible evolutionary relationships among tumor cells. Validation on simulated data shows the method can accurately separate small numbers of cell states and infer phylogenetic relationships among them. Application to a lung cancer dataset shows that the method can identify cell states corresponding to common lung tumor types and suggest possible evolutionary relationships among them that show good correspondence with our current understanding of lung tumor development.</p> <p>Conclusions</p> <p>Unmixing methods provide a way to make use of both intra-tumor heterogeneity and large probe sets for tumor phylogeny inference, establishing a new avenue towards the construction of detailed, accurate portraits of common tumor sub-types and the mechanisms by which they develop. These reconstructions are likely to have future value in discovering and diagnosing novel cancer sub-types and in identifying targets for therapeutic development.</p

Reproducibility of measurements of potential doubling time of tumour cells in the multicentre National Cancer Institute protocol T92-0045

We compared the flow cytometric measurement and analysis of the potential doubling time (Tpot) between three centres involved in the National Cancer Institute (NCI) protocol T92-0045. The primary purpose was to understand and minimize the variation within the measurement. A total of 102 specimens were selected at random from patients entered into the trial. Samples were prepared, stained, run and analysed in each centre and a single set of data analysed by all three centres. Analysis of the disc data set revealed that the measurement of labelling index (LI) was robust and reproducible. The estimation of duration of S-phase (Ts) was subject to errors of profile interpretation, particularly DNA ploidy status, and analysis. The LI dominated the variation in Tpot such that the level of final agreement, after removal of outliers and ploidy agreement, reached correlation coefficients of 0.9. The sample data showed poor agreement within each of the components of the measurement. There was some improvement when ploidy was in agreement, but correlation coefficients failed to exceed values of 0.5 for Tpot. The data suggest that observer-associated analysis of Ts and tissue processing and tumour heterogeneity were the major causes of variability in the Tpot measurement. The first two aspects can be standardized and minimized, but heterogeneity will remain a problem with biopsy techniques. © 1999 Cancer Research Campaig

Progressive Telomere Dysfunction Causes Cytokinesis Failure and Leads to the Accumulation of Polyploid Cells

Most cancer cells accumulate genomic abnormalities at a remarkably rapid rate, as they are unable to maintain their chromosome structure and number. Excessively short telomeres, a known source of chromosome instability, are observed in early human-cancer lesions. Besides telomere dysfunction, it has been suggested that a transient phase of polyploidization, in most cases tetraploidization, has a causative role in cancer. Proliferation of tetraploids can gradually generate subtetraploid lineages of unstable cells that might fire the carcinogenic process by promoting further aneuploidy and genomic instability. Given the significance of telomere dysfunction and tetraploidy in the early stages of carcinogenesis, we investigated whether there is a connection between these two important promoters of chromosomal instability. We report that human mammary epithelial cells exhibiting progressive telomere dysfunction, in a pRb deficient and wild-type p53 background, fail to complete the cytoplasmatic cell division due to the persistence of chromatin bridges in the midzone. Flow cytometry together with fluorescence in situ hybridization demonstrated an accumulation of binucleated polyploid cells upon serial passaging cells. Restoration of telomere function through hTERT transduction, which lessens the formation of anaphase bridges by recapping the chromosome ends, rescued the polyploid phenotype. Live-cell imaging revealed that these polyploid cells emerged after abortive cytokinesis due to the persistence of anaphase bridges with large intervening chromatin in the cleavage plane. In agreement with a primary role of anaphase bridge intermediates in the polyploidization process, treatment of HMEC-hTERT cells with bleomycin, which produces chromatin bridges through illegimitate repair, resulted in tetraploid binucleated cells. Taken together, we demonstrate that human epithelial cells exhibiting physiological telomere dysfunction engender tetraploid cells through interference of anaphase bridges with the completion of cytokinesis. These observations shed light on the mechanisms operating during the initial stages of human carcinogenesis, as they provide a link between progressive telomere dysfunction and tetraploidy

CDH1 promoter hypermethylation and E-cadherin protein expression in infiltrating breast cancer

BACKGROUND: The E-cadherin gene (CDH1) maps, at chromosome 16q22.1, a region often associated with loss of heterozygosity (LOH) in human breast cancer. LOH at this site is thought to lead to loss of function of this tumor suppressor gene and was correlated with decreased disease-free survival, poor prognosis, and metastasis. Differential CpG island methylation in the promoter region of the CDH1 gene might be an alternative way for the loss of expression and function of E-cadherin, leading to loss of tissue integrity, an essential step in tumor progression. METHODS: The aim of our study was to assess, by Methylation-Specific Polymerase Chain Reaction (MSP), the methylation pattern of the CDH1 gene and its possible correlation with the expression of E-cadherin and other standard immunohistochemical parameters (Her-2, ER, PgR, p53, and K-67) in a series of 79 primary breast cancers (71 infiltrating ductal, 5 infiltrating lobular, 1 metaplastic, 1 apocrine, and 1 papillary carcinoma). RESULTS: CDH1 hypermethylation was observed in 72% of the cases including 52/71 ductal, 4/5 lobular carcinomas and 1 apocrine carcinoma. Reduced levels of E-cadherin protein were observed in 85% of our samples. Although not statistically significant, the levels of E-cadherin expression tended to diminish with the CDH1 promoter region methylation. In the group of 71 ductal cancinomas, most of the cases of showing CDH1 hypermethylation also presented reduced levels of expression of ER and PgR proteins, and a possible association was observed between CDH1 methylation and ER expression (p = 0.0301, Fisher's exact test). However, this finding was not considered significant after Bonferroni correction of p-value. CONCLUSION: Our preliminary findings suggested that abnormal CDH1 methylation occurs in high frequencies in infiltrating breast cancers associated with a decrease in E-cadherin expression in a subgroup of cases characterized by loss of expression of other important genes to the mammary carcinogenesis process, probably due to the disruption of the mechanism of maintenance of DNA methylation in tumoral cells

Increasing genome instability in adrenocortical carcinoma progression with involvement of chromosomes 3, 9 and X at the adenoma stage

The investigation of chromosomal aberrations in adrenocortical tumours has been limited by the difficulties of applying classical cytogenetics to tumours with low levels of proliferation. We have therefore applied the technique of interphase cytogenetics to paraffin-embedded archival specimens of 14 adrenocortical adenomas and 13 carcinomas. Hybridizations were performed using centromere-specific probes to chromosomes 3, 4, 9, 17, 18 and X, which have been shown to be altered in other types of tumours. Chromosomal imbalance was defined on the basis of changes in both chromosome index (CI) and signal distribution (SD). Where only one of these was altered, this was classified as a tendency to gain or loss. On the basis of the analysis of optimal hybridizations, carcinomas showed gains in all chromosomes studied, five of nine showing gains in multiple chromosomes. Gains were most common in chromosomes 3, 9 and, in particular X, eight of 11 showing gain, and one a tendency to gain. Chromosomal gain was seen less commonly in adenomas, but again chromosomes 3, 9 and X were involved. Losses were infrequent, only one carcinoma showing loss of chromosome 18, and adenomas showing a tendency to loss of chromosomes 4 (two cases), 17 (one case) and 18 (two cases). Our data suggest that changes in chromosomes 3, 9 and X are early events in adrenocortical tumorigenesis, and that there is increasing chromosomal instability with tumour progression. © 1999 Cancer Research Campaig

Exploiting Clinical Trial Data Drastically Narrows the Window of Possible Solutions to the Problem of Clinical Adaptation of a Multiscale Cancer Model

The development of computational models for simulating tumor growth and response to treatment has gained significant momentum during the last few decades. At the dawn of the era of personalized medicine, providing insight into complex mechanisms involved in cancer and contributing to patient-specific therapy optimization constitute particularly inspiring pursuits. The in silico oncology community is facing the great challenge of effectively translating simulation models into clinical practice, which presupposes a thorough sensitivity analysis, adaptation and validation process based on real clinical data. In this paper, the behavior of a clinically-oriented, multiscale model of solid tumor response to chemotherapy is investigated, using the paradigm of nephroblastoma response to preoperative chemotherapy in the context of the SIOP/GPOH clinical trial. A sorting of the model's parameters according to the magnitude of their effect on the output has unveiled the relative importance of the corresponding biological mechanisms; major impact on the result of therapy is credited to the oxygenation and nutrient availability status of the tumor and the balance between the symmetric and asymmetric modes of stem cell division. The effect of a number of parameter combinations on the extent of chemotherapy-induced tumor shrinkage and on the tumor's growth rate are discussed. A real clinical case of nephroblastoma has served as a proof of principle study case, demonstrating the basics of an ongoing clinical adaptation and validation process. By using clinical data in conjunction with plausible values of model parameters, an excellent fit of the model to the available medical data of the selected nephroblastoma case has been achieved, in terms of both volume reduction and histological constitution of the tumor. In this context, the exploitation of multiscale clinical data drastically narrows the window of possible solutions to the clinical adaptation problem

HER2 expression as a potential marker for response to therapy targeted to the EGFR

Since human epidermal growth factor receptor 2 (HER2) is known to participate with the epidermal growth factor receptor (EGFR) in mitogenic signalling, we hypothesised that HER2 overexpression might indicate responsiveness to EGFR targeted therapies. MCF7 breast cancer cells transfected with the HER2 gene were subcloned to establish a set of genetically related cell lines expressing graded levels of HER2 by immunoblot analysis. The subcloned cell lines and parental MCF7 cells were characterised by their growth characteristics, and cell by cell patterns of EGFR, HER2 and HER3 expression as well as levels of phosphorylated mitogen-activated protein kinase (MAPK) and AKT by laser scanning cytometry (LSC). Growth inhibition assays were used to characterise response to EGFR targeted therapy, and to determine the relationship between therapeutic response and levels of tyrosine kinase expression. The levels of growth inhibition of AG1478 and of the AG1478-trastuzumab combinations were correlated with levels of HER2 expression among the different cell lines. Among EGFR, HER2 and HER3, HER2 overexpression was the best single predictive marker, but combinations of two markers provided additional predictive information