Tumor heterogeneity assessed by sequencing and fluorescence in situ hybridization (FISH) data.

MOTIVATION: Computational reconstruction of clonal evolution in cancers has become a crucial tool for understanding how tumors initiate and progress and how this process varies across patients. The field still struggles, however, with special challenges of applying phylogenetic methods to cancers, such as the prevalence and importance of copy number alteration (CNA) and structural variation events in tumor evolution, which are difficult to profile accurately by prevailing sequencing methods in such a way that subsequent reconstruction by phylogenetic inference algorithms is accurate. RESULTS: In this work, we develop computational methods to combine sequencing with multiplex interphase fluorescence in situ hybridization to exploit the complementary advantages of each technology in inferring accurate models of clonal CNA evolution accounting for both focal changes and aneuploidy at whole-genome scales. By integrating such information in an integer linear programming framework, we demonstrate on simulated data that incorporation of FISH data substantially improves accurate inference of focal CNA and ploidy changes in clonal evolution from deconvolving bulk sequence data. Analysis of real glioblastoma data for which FISH, bulk sequence and single cell sequence are all available confirms the power of FISH to enhance accurate reconstruction of clonal copy number evolution in conjunction with bulk and optionally single-cell sequence data. AVAILABILITY AND IMPLEMENTATION: Source code is available on Github at https://github.com/CMUSchwartzLab/FISH_deconvolution. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

High Levels of Chromosomal Copy Number Alterations and TP53 Mutations Correlate with Poor Outcome in Younger Breast Cancer Patients.

Prognosis in young patients with breast cancer is generally poor, yet considerable differences in clinical outcomes between individual patients exist. To understand the genetic basis of the disparate clinical courses, tumors were collected from 34 younger women, 17 with good and 17 with poor outcomes, as determined by disease-specific survival during a follow-up period of 17 years. The clinicopathologic parameters of the tumors were complemented with DNA image cytometry profiles, enumeration of copy numbers of eight breast cancer genes by multicolor fluorescence in situ hybridization, and targeted sequence analysis of 563 cancer genes. Both groups included diploid and aneuploid tumors. The degree of intratumor heterogeneity was significantly higher in aneuploid versus diploid cases, and so were gains of the oncogenes MYC and ZNF217. Significantly more copy number alterations were observed in the group with poor outcome. Almost all tumors in the group with long survival were classified as luminal A, whereas triple-negative tumors predominantly occurred in the short survival group. Mutations in PIK3CA were more common in the group with good outcome, whereas TP53 mutations were more frequent in patients with poor outcomes. This study shows that TP53 mutations and the extent of genomic imbalances are associated with poor outcome in younger breast cancer patients and thus emphasize the central role of genomic instability vis-a-vis tumor aggressiveness.

Single-cell genetic analysis of clonal dynamics in colorectal adenomas indicates CDX2 gain as a predictor of recurrence.

Colorectal adenomas are common precancerous lesions with the potential for malignant transformation to colorectal adenocarcinoma. Endoscopic polypectomy provides an opportunity for cancer prevention; however, recurrence rates are high. We collected formalin-fixed paraffin-embedded tissue of 15 primary adenomas with recurrence, 15 adenomas without recurrence, and 14 matched pair samples (primary adenoma and the corresponding recurrent adenoma). The samples were analysed by array-comparative genomic hybridisation (aCGH) and single-cell multiplex interphase fluorescence in situ hybridisation (miFISH) to understand clonal evolution, to examine the dynamics of copy number alterations (CNAs) and to identify molecular markers for recurrence prediction. The miFISH probe panel consisted of 14 colorectal carcinogenesis-relevant genes (COX2, PIK3CA, APC, CLIC1, EGFR, MYC, CCND1, CDX2, CDH1, TP53, HER2, SMAD7, SMAD4 and ZNF217), and a centromere probe (CEP10). The aCGH analysis confirmed the genetic landscape typical for colorectal tumorigenesis, that is, CNAs of chromosomes 7, 13q, 18 and 20q. Focal aberrations (≤10 Mbp) were mapped to chromosome bands 6p22.1-p21.33 (33.3%), 7q22.1 (31.4%) and 16q21 (29.4%). MiFISH detected gains of EGFR (23.6%), CDX2 (21.8%) and ZNF217 (18.2%). Most adenomas exhibited a major clone population which was accompanied by multiple smaller clone populations. Gains of CDX2 were exclusively seen in primary adenomas with recurrence (25%) compared to primary adenomas without recurrence (0%). Generation of phylogenetic trees for matched pair samples revealed four distinct patterns of clonal dynamics. In conclusion, adenoma development and recurrence are complex genetic processes driven by multiple CNAs whose evaluations by miFISH, with emphasis on CDX2, might serve as a predictor of recurrence.

Aneuploidy, TP53 mutation, and amplification of MYC correlate with increased intratumor heterogeneity and poor prognosis of breast cancer patients.

The clinical course of breast cancer varies from one patient to another. Currently, the choice of therapy relies on clinical parameters and histological and molecular tumor features. Alas, these markers are informative in only a subset of patients. Therefore, additional predictors of disease outcome would be valuable for treatment stratification. Extensive studies showed that the degree of variation of the nuclear DNA content, i.e., aneuploidy, determines prognosis. Our aim was to further elucidate the molecular basis of aneuploidy. We analyzed five diploid and six aneuploid tumors with more than 20 years of follow-up. By performing FISH with a multiplexed panel of 10 probes to enumerate copy numbers in individual cells, and by sequencing 563 cancer-related genes, we analyzed how aneuploidy is linked to intratumor heterogeneity. In our cohort, none of the patients with diploid tumors died of breast cancer during follow-up in contrast to four of six patients with aneuploid tumors (mean survival 86.4 months). The FISH analysis showed markedly increased genomic instability and intratumor heterogeneity in aneuploid tumors. MYC gain was observed in only 20% of the diploid cancers, while all aneuploid cases showed a gain. The mutation burden was similar in diploid and aneuploid tumors, however, TP53 mutations were not observed in diploid tumors, but in all aneuploid tumors in our collective. We conclude that quantitative measurements of intratumor heterogeneity by multiplex FISH, detection of MYC amplification and TP53 mutation could augment prognostication in breast cancer patients.

The evolution of single cell-derived colorectal cancer cell lines is dominated by the continued selection of tumor-specific genomic imbalances, despite random chromosomal instability.

Intratumor heterogeneity is a major challenge in cancer treatment. To decipher patterns of chromosomal heterogeneity, we analyzed six colorectal cancer cell lines by multiplex interphase FISH (miFISH). The mismatch-repair-deficient cell lines DLD-1 and HCT116 had the most stable copy numbers, whereas aneuploid cell lines (HT-29, SW480, SW620 and H508) displayed a higher degree of instability. We subsequently assessed the clonal evolution of single cells in two colorectal carcinoma cell lines, SW480 and HT-29, which both have aneuploid karyotypes but different degrees of chromosomal instability. The clonal compositions of the single cell-derived daughter lines, as assessed by miFISH, differed for HT-29 and SW480. Daughters of HT-29 were stable, clonal, with little heterogeneity. Daughters of SW480 were more heterogeneous, with the single cell-derived daughter lines separating into two distinct populations with different ploidy (hyper-diploid and near-triploid), morphology, gene expression and tumorigenicity. To better understand the evolutionary trajectory for the two SW480 populations, we constructed phylogenetic trees which showed ongoing instability in the daughter lines. When analyzing the evolutionary development over time, most single cell-derived daughter lines maintained their major clonal pattern, with the exception of one daughter line that showed a switch involving a loss of APC. Our meticulous analysis of the clonal evolution and composition of these colorectal cancer models shows that all chromosomes are subject to segregation errors, however, specific net genomic imbalances are maintained. Karyotype evolution is driven by the necessity to arrive at and maintain a specific plateau of chromosomal copy numbers as the drivers of carcinogenesis.

Phylogenetic analysis of multiple FISH markers in oral tongue squamous cell carcinoma suggests that a diverse distribution of copy number changes is associated with poor prognosis.

Oral tongue squamous cell carcinoma (OTSCC) is associated with poor prognosis. To improve prognostication, we analyzed four gene probes (TERC, CCND1, EGFR and TP53) and the centromere probe CEP4 as a marker of chromosomal instability, using fluorescence in situ hybridization (FISH) in single cells from the tumors of sixty-five OTSCC patients (Stage I, n = 15; Stage II, n = 30; Stage III, n = 7; Stage IV, n = 13). Unsupervised hierarchical clustering of the FISH data distinguished three clusters related to smoking status. Copy number increases of all five markers were found to be correlated to non-smoking habits, while smokers in this cohort had low-level copy number gains. Using the phylogenetic modeling software FISHtrees, we constructed models of tumor progression for each patient based on the four gene probes. Then, we derived test statistics on the models that are significant predictors of disease-free and overall survival, independent of tumor stage and smoking status in multivariate analysis. The patients whose tumors were modeled as progressing by a more diverse distribution of copy number changes across the four genes have poorer prognosis. This is consistent with the view that multiple genetic pathways need to become deregulated in order for cancer to progress.

Inferring models of multiscale copy number evolution for single-tumor phylogenetics.

MOTIVATION: Phylogenetic algorithms have begun to see widespread use in cancer research to reconstruct processes of evolution in tumor progression. Developing reliable phylogenies for tumor data requires quantitative models of cancer evolution that include the unusual genetic mechanisms by which tumors evolve, such as chromosome abnormalities, and allow for heterogeneity between tumor types and individual patients. Previous work on inferring phylogenies of single tumors by copy number evolution assumed models of uniform rates of genomic gain and loss across different genomic sites and scales, a substantial oversimplification necessitated by a lack of algorithms and quantitative parameters for fitting to more realistic tumor evolution models. RESULTS: We propose a framework for inferring models of tumor progression from single-cell gene copy number data, including variable rates for different gain and loss events. We propose a new algorithm for identification of most parsimonious combinations of single gene and single chromosome events. We extend it via dynamic programming to include genome duplications. We implement an expectation maximization (EM)-like method to estimate mutation-specific and tumor-specific event rates concurrently with tree reconstruction. Application of our algorithms to real cervical cancer data identifies key genomic events in disease progression consistent with prior literature. Classification experiments on cervical and tongue cancer datasets lead to improved prediction accuracy for the metastasis of primary cervical cancers and for tongue cancer survival. AVAILABILITY AND IMPLEMENTATION: Our software (FISHtrees) and two datasets are available at ftp://ftp.ncbi.nlm.nih.gov/pub/FISHtrees.

Single-cell genetic analysis reveals insights into clonal development of prostate cancers and indicates loss of PTEN as a marker of poor prognosis.

Gauging the risk of developing progressive disease is a major challenge in prostate cancer patient management. We used genetic markers to understand genomic alteration dynamics during disease progression. By using a novel, advanced, multicolor fluorescence in situ hybridization approach, we enumerated copy numbers of six genes previously identified by array comparative genomic hybridization to be involved in aggressive prostate cancer [TBL1XR1, CTTNBP2, MYC (alias c-myc), PTEN, MEN1, and PDGFB] in six nonrecurrent and seven recurrent radical prostatectomy cases. An ERG break-apart probe to detect TMPRSS2-ERG fusions was included. Subsequent hybridization of probe panels and cell relocation resulted in signal counts for all probes in each individual cell analyzed. Differences in the degree of chromosomal and genomic instability (ie, tumor heterogeneity) or the percentage of cells with TMPRSS2-ERG fusion between samples with or without progression were not observed. Tumors from patients that progressed had more chromosomal gains and losses, and showed a higher degree of selection for a predominant clonal pattern. PTEN loss was the most frequent aberration in progressers (57%), followed by TBL1XR1 gain (29%). MYC gain was observed in one progresser, which was the only lesion with an ERG gain, but no TMPRSS2-ERG fusion. According to our results, a probe set consisting of PTEN, MYC, and TBL1XR1 would detect progressers with 86% sensitivity and 100% specificity. This will be evaluated further in larger studies.

Algorithms to model single gene, single chromosome, and whole genome copy number changes jointly in tumor phylogenetics.

We present methods to construct phylogenetic models of tumor progression at the cellular level that include copy number changes at the scale of single genes, entire chromosomes, and the whole genome. The methods are designed for data collected by fluorescence in situ hybridization (FISH), an experimental technique especially well suited to characterizing intratumor heterogeneity using counts of probes to genetic regions frequently gained or lost in tumor development. Here, we develop new provably optimal methods for computing an edit distance between the copy number states of two cells given evolution by copy number changes of single probes, all probes on a chromosome, or all probes in the genome. We then apply this theory to develop a practical heuristic algorithm, implemented in publicly available software, for inferring tumor phylogenies on data from potentially hundreds of single cells by this evolutionary model. We demonstrate and validate the methods on simulated data and published FISH data from cervical cancers and breast cancers. Our computational experiments show that the new model and algorithm lead to more parsimonious trees than prior methods for single-tumor phylogenetics and to improved performance on various classification tasks, such as distinguishing primary tumors from metastases obtained from the same patient population.

Phylogenetic analysis of multiprobe fluorescence in situ hybridization data from tumor cell populations.

MOTIVATION: Development and progression of solid tumors can be attributed to a process of mutations, which typically includes changes in the number of copies of genes or genomic regions. Although comparisons of cells within single tumors show extensive heterogeneity, recurring features of their evolutionary process may be discerned by comparing multiple regions or cells of a tumor. A useful source of data for studying likely progression of individual tumors is fluorescence in situ hybridization (FISH), which allows one to count copy numbers of several genes in hundreds of single cells. Novel algorithms for interpreting such data phylogenetically are needed, however, to reconstruct likely evolutionary trajectories from states of single cells and facilitate analysis of tumor evolution. RESULTS: In this article, we develop phylogenetic methods to infer likely models of tumor progression using FISH copy number data and apply them to a study of FISH data from two cancer types. Statistical analyses of topological characteristics of the tree-based model provide insights into likely tumor progression pathways consistent with the prior literature. Furthermore, tree statistics from the resulting phylogenies can be used as features for prediction methods. This results in improved accuracy, relative to unstructured gene copy number data, at predicting tumor state and future metastasis. AVAILABILITY: Source code for software that does FISH tree building (FISHtrees) and the data on cervical and breast cancer examined here are available at ftp://ftp.ncbi.nlm.nih.gov/pub/FISHtrees. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Single-cell genetic analysis of ductal carcinoma in situ and invasive breast cancer reveals enormous tumor heterogeneity yet conserved genomic imbalances and gain of MYC during progression.

Ductal carcinoma in situ (DCIS) is a precursor lesion of invasive ductal carcinoma (IDC) of the breast. To understand the dynamics of genomic alterations in this progression, we used four multicolor fluorescence in situ hybridization probe panels consisting of the oncogenes COX2, MYC, HER2, CCND1, and ZNF217 and the tumor suppressor genes DBC2, CDH1, and TP53 to visualize copy number changes in 13 cases of synchronous DCIS and IDC based on single-cell analyses. The DCIS had a lower degree of chromosomal instability than the IDC. Despite enormous intercellular heterogeneity in DCIS and IDC, we observed signal patterns consistent with a nonrandom distribution of genomic imbalances. CDH1 was most commonly lost, and gain of MYC emerged during progression from DCIS to IDC. Four of 13 DCISs showed identical clonal imbalances in the IDCs. Six cases revealed a switch, and in four of those, the IDC had acquired a gain of MYC. In one case, the major clone in the IDC was one of several clones in the DCIS, and in another case, the major clone in the DCIS became one of the two major clones in the IDC. Despite considerable chromosomal instability, in most cases the evolution from DCIS to IDC is determined by recurrent patterns of genomic imbalances, consistent with a biological continuum.

Multicentric dermatofibrosarcoma protuberans in patients with adenosine deaminase-deficient severe combined immune deficiency.

BACKGROUND: Dermatofibrosarcoma protuberans (DFSP) is a rare malignant skin tumor associated with a characteristic chromosomal translocation (t[17;22][q22;q13]) resulting in the COL1A1-platelet-derived growth factor ß(PDGFB) fusion gene. This malignancy is rarely diagnosed in childhood. OBJECTIVE: We observed an unexpected high incidence of this DFSP in children affected with adenosine deaminase-deficient severe combined immunodeficiency (ADA-SCID) and set out to evaluate the association of these 2 clinical entities. METHODS: Twelve patients with ADA-SCID were evaluated with a complete dermatologic examination and skin biopsy when indicated. Conventional cytogenetic and molecular analyses (fluorescence in situ hybridization, RT-PCR, or both) were performed when possible. RESULTS: Eight patients were found to have DFSP. Six patients had multicentric involvement (4-15 lesions), primarily of the trunk and extremities. Most lesions presented as 2- to 15-mm, round atrophic plaques. Nodular lesions were present in 3 patients. In all cases CD34 expression was diffusely positive, and diagnosis was confirmed either by means of cytogenetic analysis, molecular testing, or both. The characteristic DFSP-associated translocation, t(17;22)(q22;q13), was identified in 6 patients; results of fluorescence in situ hybridization were positive for fusion of the COL1A1 and PDGFB loci in 7 patients; and RT-PCR showed the COL1A1-PDGFB fusion transcript in 6 patients. CONCLUSIONS: We describe a previously unrecognized association between ADA-SCID and DFSP with unique features, such as multicentricity and occurrence in early age. We hypothesize that the t(17;22)(q22;q13) translocation that results in dermal overexpression of PDGFB and favors the development of fibrotic tumors might arise because of the known DNA repair defect in patients with ADA-SCID. Although the natural course of DFSP in the setting of ADA-SCID is unknown, this observation should prompt regular screening for DFSP in patients with ADA-SCID.

Extrachromosomal Amplification of Human Papillomavirus Episomes Is a Mechanism of Cervical Carcinogenesis.

SIGNIFICANCE: Multimers of the HPV genome are generated in cervical tumors replicating as extrachromosomal episomes, which is associated with deletion and rearrangement of the HPV genome and provides a mechanism for oncogenesis without integration.

Horizontal transmission and retention of malignancy, as well as functional human genes, after spontaneous fusion of human glioblastoma and hamster host cells in vivo.

Cell fusion in vitro has been used to study cancer, gene mapping and regulation, and the production of antibodies via hybridomas. However, in-vivo heterosynkaryon formation by cell-cell fusion has received less attention. This investigation describes the spontaneous fusion of a human glioblastoma with normal hamster cells after xenogeneic transplantation, resulting in malignant cells that express both human and hamster genes and gene products, and retention of glioblastoma traits with an enhanced ability to metastasize. Three of 7 human genes found showed translation of their proteins during serial propagation in vivo or in vitro for years; namely, CD74, CXCR4 and PLAGL2, each implicated with malignancy or glioblastoma. This supports the thesis that genetic hybridization of cancer and normal cells can transmit malignancy and also, as first described herein, regulatory genes involved in the tumor's organotypic morphology. Evidence also is increasing that even cell-free human cancer DNA can induce malignancy and transfer genetic information to normal cells. Hence, we posit that the transfer of genetic information between tumor and stromal cells, whether by cell-cell fusion or other mechanisms, is implicated in the progression of malignancy, and may further define the crosstalk between cancer cells and their stromal neighbors.

Genome and transcriptome profiles of CD133-positive colorectal cancer cells.

Colorectal carcinomas (CRC) might be organized hierarchically and contain a subpopulation of tumorigenic, putative cancer stem cells that are CD133 positive. We studied the biological and genetic characteristics of such cells in CRC cell lines and primary tumors. Three CRC cell lines were sorted in CD133 positive and negative fractions. The respective genetic aberration profiles were studied using array comparative genomic hybridization (aCGH) and expression profiling. Tumorigenicity for each cellular population was tested by injection into nude mice. Additionally, we compared CD133+ and CD133- cells of 12 primary colorectal tumors using laser capture microdissection and aCGH. Three of five CRC cell lines displayed both CD133+ and CD133- cells, but tumorigenicity of these subfractions did not differ significantly and aCGH revealed essentially identical genomic imbalances. However, 96 genes were differentially expressed between the two populations. Array comparative genomic hybridization analysis after laser capture microdissection of CD133+ and CD133- areas in primary colorectal tumors revealed genetic differences in 7 of 12 cases. The use of cell lines for studying genomic alterations that define cancer stem cell characteristics, therefore, seems questionable. In contrast, CD133+ cells in primary cancer samples showed a unique genomic aberration profile. In conclusion, our data suggest that CD133 positivity defines a genetically distinct cellular compartment in primary CRC, which potentially includes tumor initiating cells.

Single-cell resolved ploidy and chromosomal aberrations in nonalcoholic steatohepatitis-(NASH) induced hepatocellular carcinoma and its precursor lesions.

Nonalcoholic steatohepatitis (NASH)-induced hepatocellular carcinoma (HCC) and its precursor, nonalcoholic fatty liver disease (NAFLD) are an unmet health issue due to widespread obesity. We assessed copy number changes of genes associated with hepatocarcinogenesis and oxidative pathways at a single-cell level. Eleven patients with NASH-HCC and 11 patients with NAFLD were included. Eight probes were analyzed using multiplex interphase fluorescence in situ hybridization (miFISH), single-cell imaging and phylogenetic tree modelling: Telomerase reverse transcriptase (TERT), C-Myc (MYC), hepatocyte growth factor receptor tyrosine kinase (MET), tumor protein 53 (TP53), cyclin D1 (CCND1), human epidermal growth factor receptor 2 (HER2), the fragile histidine triad gene (FHIT) and FRA16D oxidoreductase (WWOX). Each NASH-HCC tumor had up to 14 distinct clonal signal patterns indicating multiclonality, which correlated with high tumor grade. Changes frequently observed were TP53 losses, 45%; MYC gains, 36%; WWOX losses, 36%; and HER2 gains, 18%. Whole-genome duplications were frequent (82%) with aberrant tetraploid cells evolving from diploid ancestors. Non-tumorous NAFLD/NASH biopsies did not harbor clonal copy number changes. Fine mapping of NASH-HCC using single-cell multiplex FISH shows that branched tumor evolution involves genome duplication and that multiclonality increases with tumor grade. The loss of oxidoreductase WWOX and HER2 gains could be potentially associated with NASH-induced hepatocellular carcinoma.

Genomic instability and oncogene amplifications in colorectal adenomas predict recurrence and synchronous carcinoma.

Individual colorectal adenomas have different propensities to progress to invasive disease. In this study, we explored whether these differences could be explained by gene copy number alterations. We evaluated 18 adenomas of patients without synchronous or subsequent carcinoma (6.5 years follow-up), 23 adenomas of carcinoma patients, and 6 related carcinomas. All samples were measured for their DNA ploidy status. Centromere probes for chromosomes 17 and 18, as well as gene-specific probes for SMAD7, EGFR, NCOA3, TP53, MYC, and RAB20 were assessed by multicolor fluorescence in situ hybridization. An increased genomic instability index of CEP17, SMAD7, and EGFR, as well as TP53 deletions and MYC amplifications defined adenomas of patients with synchronous carcinoma (P<0.05). Diploid NCOA3 signal counts were associated with longer adenoma recurrence-free surveillance (P=0.042). In addition, NCOA3, MYC, EGFR, and RAB20 amplifications, as well as TP53 deletions correlated with increased DNA stem line values and/or aneuploidy in adenomas (P<0.05). Furthermore, aberrations of NCOA3, MYC, and RAB20 were associated with histopathologically defined high-risk adenomas (P<0.05). RAB20 amplifications were also correlated with high-grade dysplastic adenomas (P=0.002). We conclude that genomic instability in colorectal adenomas is reflected by EGFR, MYC, NCOA3, and RAB20 amplifications that do correlate with histomorphological features and are indicative for adenoma recurrence and the presence of synchronous carcinomas.

Joint Clustering of Single-Cell Sequencing and Fluorescence In Situ Hybridization Data for Reconstructing Clonal Heterogeneity in Cancers.

Aneuploidy and whole genome duplication (WGD) events are common features of cancers associated with poor outcomes, but the ways they influence trajectories of clonal evolution are poorly understood. Phylogenetic methods for reconstructing clonal evolution from genomic data have proven a powerful tool for understanding how clonal evolution occurs in the process of cancer progression, but extant methods so far have limited the ability to resolve tumor evolution via ploidy changes. This limitation exists in part because single-cell DNA-sequencing (scSeq), which has been crucial to developing detailed profiles of clonal evolution, has difficulty in resolving ploidy changes and WGD. Multiplex interphase fluorescence in situ hybridization (miFISH) provides a more unambiguous signal of single-cell ploidy changes but it is limited to profiling small numbers of single markers. Here, we develop a joint clustering method to combine these two data sources with the goal of better resolving ploidy changes in tumor evolution. We develop a probabilistic framework to maximize the probability of latent variables given the pre-clustered datasets, which we optimize via Markov chain Monte Carlo sampling combined with linear regression. We validate the method by using simulated data derived from a glioblastoma (GBM) case profiled by both scSeq and miFISH. We further apply the method to two GBM cases with scSeq and miFISH data by reconstructing a phylogenetic tree from the joint clustering results, demonstrating their synergistic value in understanding how focal copy number changes and WGD events can collectively contribute to tumor progression.

Single Cell Genetic Profiling of Tumors of Breast Cancer Patients Aged 50 Years and Older Reveals Enormous Intratumor Heterogeneity Independent of Individual Prognosis.

PURPOSE: Older breast cancer patients are underrepresented in cancer research even though the majority (81.4%) of women dying of breast cancer are 55 years and older. Here we study a common phenomenon observed in breast cancer which is a large inter- and intratumor heterogeneity; this poses a tremendous clinical challenge, for example with respect to treatment stratification. To further elucidate genomic instability and tumor heterogeneity in older patients, we analyzed the genetic aberration profiles of 39 breast cancer patients aged 50 years and older (median 67 years) with either short (median 2.4 years) or long survival (median 19 years). The analysis was based on copy number enumeration of eight breast cancer-associated genes using multiplex interphase fluorescence in situ hybridization (miFISH) of single cells, and by targeted next-generation sequencing of 563 cancer-related genes. RESULTS: We detected enormous inter- and intratumor heterogeneity, yet maintenance of common cancer gene mutations and breast cancer specific chromosomal gains and losses. The gain of COX2 was most common (72%), followed by MYC (69%); losses were most prevalent for CDH1 (74%) and TP53 (69%). The degree of intratumor heterogeneity did not correlate with disease outcome. Comparing the miFISH results of diploid with aneuploid tumor samples significant differences were found: aneuploid tumors showed significantly higher average signal numbers, copy number alterations (CNAs) and instability indices. Mutations in PIKC3A were mostly restricted to luminal A tumors. Furthermore, a significant co-occurrence of CNAs of DBC2/MYC, HER2/DBC2 and HER2/TP53 and mutual exclusivity of CNAs of HER2 and PIK3CA mutations and CNAs of CCND1 and PIK3CA mutations were revealed. CONCLUSION: Our results provide a comprehensive picture of genome instability profiles with a large variety of inter- and intratumor heterogeneity in breast cancer patients aged 50 years and older. In most cases, the distribution of chromosomal aneuploidies was consistent with previous results; however, striking exceptions, such as tumors driven by exclusive loss of chromosomes, were identified.

TP53 loss initiates chromosomal instability in fallopian tube epithelial cells.

High-grade serous ovarian cancer (HGSOC) originates in the fallopian tube epithelium and is characterized by ubiquitous TP53 mutation and extensive chromosomal instability (CIN). However, direct causes of CIN, such as mutations in DNA replication and mitosis genes, are rare in HGSOC. We therefore asked whether oncogenic mutations that are common in HGSOC can indirectly drive CIN in non-transformed human fallopian tube epithelial cells. To model homologous recombination deficient HGSOC, we sequentially mutated TP53 and BRCA1 then overexpressed MYC. Loss of p53 function alone was sufficient to drive the emergence of subclonal karyotype alterations. TP53 mutation also led to global gene expression changes, influencing modules involved in cell cycle commitment, DNA replication, G2/M checkpoint control and mitotic spindle function. Both transcriptional deregulation and karyotype diversity were exacerbated by loss of BRCA1 function, with whole-genome doubling events observed in independent p53/BRCA1-deficient lineages. Thus, our observations indicate that loss of the key tumour suppressor TP53 is sufficient to deregulate multiple cell cycle control networks and thereby initiate CIN in pre-malignant fallopian tube epithelial cells. This article has an associated First Person interview with the first author of the paper.