Circulating tumor cells (CTC) released into blood from primary cancers and metastases reflect the current status of tumor genotypes, which are prone to changes. Here, we conducted the first comprehensive genomic profiling of CTCs using array-comparative genomic hybridization (CGH) and next-generation sequencing. We used the U.S. Food and Drug Administration-cleared CellSearch system, which detected CTCs in 21 of 37 patients (range, 1-202/7.5 mL sample) with stage IV colorectal carcinoma. In total, we were able to isolate 37 intact CTCs from six patients and identified in those multiple colorectal cancer-associated copy number changes, many of which were also present in the respective primary tumor. We then used massive parallel sequencing of a panel of 68 colorectal cancer-associated genes to compare the mutation spectrum in the primary tumors, metastases, and the corresponding CTCs from two of these patients. Mutations in known driver genes [e.g., adenomatous polyposis coli (APC), KRAS, or PIK3CA] found in the primary tumor and metastasis were also detected in corresponding CTCs. However, we also observed mutations exclusively in CTCs. To address whether these mutations were derived from a small subclone in the primary tumor or represented new variants of metastatic cells, we conducted additional deep sequencing of the primary tumor and metastasis and applied a customized statistical algorithm for analysis. We found that most mutations initially found only in CTCs were also present at subclonal level in the primary tumors and metastases from the same patient. This study paves the way to use CTCs as a liquid biopsy in patients with cancer, providing more effective options to monitor tumor genomes that are prone to change during progression, treatment, and relapse. Cancer Res; 73(10); 2965-75. Ó2013 AACR.

With the increasing number of available predictive biomarkers, clinical management of cancer is becoming increasingly reliant on the accurate serial monitoring of tumor genotypes. We tested whether tumor-specific copy number changes can be inferred from the peripheral blood of patients with cancer. To this end, we determined the plasma DNA size distribution and the fraction of mutated plasma DNA fragments with deep sequencing and an ultrasensitive mutation-detection method, i.e., the Beads, Emulsion, Amplification, and Magnetics (BEAMing) assay. When analyzing the plasma DNA of 32 patients with Stage IV colorectal carcinoma, we found that a subset of the patients (34.4%) had a biphasic size distribution of plasma DNA fragments that was associated with increased circulating tumor cell numbers and elevated concentration of mutated plasma DNA fragments. In these cases, we were able to establish genome-wide tumor-specific copy number alterations directly from plasma DNA. Thus, we could analyze the current copy number status of the tumor genome, which was in some cases many years after diagnosis of the primary tumor. An unexpected finding was that not all patients with progressive metastatic disease appear to release tumor DNA into the circulation in measurable quantities. When we analyzed plasma DNA from 35 patients with metastatic breast cancer, we made similar observations suggesting that our approach may be applicable to a variety of tumor entities. This is the first description of such a biphasic distribution in a surprisingly high proportion of cancer patients which may have important implications for tumor diagnosis and monitoring.

Clinical DNA is often available in limited quantities requiring whole-genome amplification for subsequent genome-wide assessment of copy-number variation (CNV) by array-CGH. In pre-implantation diagnosis and analysis of micrometastases, even merely single cells are available for analysis. However, procedures allowing high-resolution analyses of CNVs from single cells well below resolution limits of conventional cytogenetics are lacking. Here, we applied amplification products of single cells and of cell pools (5 or 10 cells) from patients with developmental delay, cancer cell lines and polar bodies to various oligo tiling array platforms with a median probe spacing as high as 65 bp. Our high-resolution analyses reveal that the low amounts of template DNA do not result in a completely unbiased whole genome amplification but that stochastic amplification artifacts, which become more obvious on array platforms with tiling path resolution, cause significant noise. We implemented a new evaluation algorithm specifically for the identification of small gains and losses in such very noisy ratio profiles. Our data suggest that when assessed with sufficiently sensitive methods high-resolution oligo-arrays allow a reliable identification of CNVs as small as 500 kb in cell pools (5 or 10 cells), and of 2.6–3.0 Mb in single cells.