Determining breast cancer histological grade from RNA-sequencing data.

BACKGROUND: The histologic grade (HG) of breast cancer is an established prognostic factor. The grade is usually reported on a scale ranging from 1 to 3, where grade 3 tumours are the most aggressive. However, grade 2 is associated with an intermediate risk of recurrence, and carries limited information for clinical decision-making. Patients classified as grade 2 are at risk of both under- and over-treatment. METHODS: RNA-sequencing analysis was conducted in a cohort of 275 women diagnosed with invasive breast cancer. Multivariate prediction models were developed to classify tumours into high and low transcriptomic grade (TG) based on gene- and isoform-level expression data from RNA-sequencing. HG2 tumours were reclassified according to the prediction model and a recurrence-free survival analysis was performed by the multivariate Cox proportional hazards regression model to assess to what extent the TG model could be used to stratify patients. The prediction model was validated in N=487 breast cancer cases from the The Cancer Genome Atlas (TCGA) data set. Differentially expressed genes and isoforms associated with HGs were analysed using linear models. RESULTS: The classification of grade 1 and grade 3 tumours based on RNA-sequencing data achieved high accuracy (area under the receiver operating characteristic curve = 0.97). The association between recurrence-free survival rate and HGs was confirmed in the study population (hazard ratio of grade 3 versus 1 was 2.62 with 95 % confidence interval = 1.04-6.61). The TG model enabled us to reclassify grade 2 tumours as high TG and low TG gene or isoform grade. The risk of recurrence in the high TG group of grade 2 tumours was higher than in low TG group (hazard ratio = 2.43, 95 % confidence interval = 1.13-5.20). We found 8200 genes and 13,809 isoforms that were differentially expressed between HG1 and HG3 breast cancer tumours. CONCLUSIONS: Gene- and isoform-level expression data from RNA-sequencing could be utilised to differentiate HG1 and HG3 tumours with high accuracy. We identified a large number of novel genes and isoforms associated with HG. Grade 2 tumours could be reclassified as high and low TG, which has the potential to reduce over- and under-treatment if implemented clinically.

Majority of differentially expressed genes are down-regulated during malignant transformation in a four-stage model.

The transformation of normal cells to malignant, metastatic tumor cells is a multistep process caused by the sequential acquirement of genetic changes. To identify these changes, we compared the transcriptomes and levels and distribution of proteins in a four-stage cell model of isogenically matched normal, immortalized, transformed, and metastatic human cells, using deep transcriptome sequencing and immunofluorescence microscopy. The data show that ∼6% (n = 1,357) of the human protein-coding genes are differentially expressed across the stages in the model. Interestingly, the majority of these genes are down-regulated, linking malignant transformation to dedifferentiation. The up-regulated genes are mainly components that control cellular proliferation, whereas the down-regulated genes consist of proteins exposed on or secreted from the cell surface. As many of the identified gene products control basic cellular functions that are defective in cancers, the data provide candidates for follow-up studies to investigate their functional roles in tumor formation. When we further compared the expression levels of four of the identified proteins in clinical cancer cohorts, similar differences were observed between benign and cancer cells, as in the cell model. This shows that this comprehensive demonstration of the molecular changes underlying malignant transformation is a relevant model to study the process of tumor formation.

Exome sequencing of contralateral breast cancer identifies metastatic disease.

Women with contralateral breast cancer (CBC) have significantly worse prognosis compared to women with unilateral cancer. A possible explanation of the poor prognosis of patients with CBC is that in a subset of patients, the second cancer is not a new primary tumor but a metastasis of the first cancer that has potentially obtained aggressive characteristics through selection of treatment. Exome and whole-genome sequencing of solid tumors has previously been used to investigate the clonal relationship between primary tumors and metastases in several diseases. In order to assess the relationship between the first and the second cancer, we performed exome sequencing to identify somatic mutations in both first and second cancers, and compared paired normal tissue of 25 patients with metachronous CBC. For three patients, we identified shared somatic mutations indicating a common clonal origin thereby demonstrating that the second tumor is a metastasis of the first cancer, rather than a new primary cancer. Accordingly, these patients all developed distant metastasis within 3 years of the second diagnosis, compared with 7 out of 22 patients with non-shared somatic profiles. Genomic profiling of both tumors help the clinicians distinguish between true CBCs and subsequent metastases.

Sequencing of breast cancer stem cell populations indicates a dynamic conversion between differentiation states in vivo.

INTRODUCTION: The cancer stem cell model implies a hierarchical organization within breast tumors maintained by cancer stem-like cells (CSCs). Accordingly, CSCs are a subpopulation of cancer cells with capacity for self-renewal, differentiation and tumor initiation. These cells can be isolated through the phenotypic markers CD44+/CD24-, expression of ALDH1 and an ability to form nonadherent, multicellular spheres in vitro. However, controversies to describe the stem cell model exist; it is unclear whether the tumorigenicity of CSCs in vivo is solely a proxy for a certain genotype. Moreover, in vivo evidence is lacking to fully define the reversibility of CSC differentiation. METHODS: In order to answer these questions, we undertook exome sequencing of CSCs from 12 breast cancer patients, along with paired primary tumor samples. As suggested by stem classical cell biology, we assumed that the number of mutations in the CSC subpopulation should be lower and distinct compared to the differentiated tumor cells with higher proliferation. RESULTS: Our analysis revealed that the majority of somatic mutations are shared between CSCs and bulk primary tumor, with similar frequencies in the two. CONCLUSIONS: The data presented here exclude the possibility that CSCs are only a phenotypic consequence of certain somatic mutations, that is a distinct and non-reversible population of cells. In addition, our results imply that CSCs must be a population of cells that can dynamically switch from differentiated tumor cells, and vice versa. This finding increases our understanding of CSC function in tumor heterogeneity and the importance of identifying drugs to counter de-differentiation rather than targeting CSCs.

Antibody-based protein profiling of the human chromosome 21.

The Human Proteome Project has been proposed to create a knowledge-based resource based on a systematical mapping of all human proteins, chromosome by chromosome, in a gene-centric manner. With this background, we here describe the systematic analysis of chromosome 21 using an antibody-based approach for protein profiling using both confocal microscopy and immunohistochemistry, complemented with transcript profiling using next generation sequencing data. We also describe a new approach for protein isoform analysis using a combination of antibody-based probing and isoelectric focusing. The analysis has identified several genes on chromosome 21 with no previous evidence on the protein level, and the isoform analysis indicates that a large fraction of human proteins have multiple isoforms. A chromosome-wide matrix is presented with status for all chromosome 21 genes regarding subcellular localization, tissue distribution, and molecular characterization of the corresponding proteins. The path to generate a chromosome-specific resource, including integrated data from complementary assay platforms, such as mass spectrometry and gene tagging analysis, is discussed.

Contribution of antibody-based protein profiling to the human Chromosome-centric Proteome Project (C-HPP).

A gene-centric Human Proteome Project has been proposed to characterize the human protein-coding genes in a chromosome-centered manner to understand human biology and disease. Here, we report on the protein evidence for all genes predicted from the genome sequence based on manual annotation from literature (UniProt), antibody-based profiling in cells, tissues and organs and analysis of the transcript profiles using next generation sequencing in human cell lines of different origins. We estimate that there is good evidence for protein existence for 69% (n = 13985) of the human protein-coding genes, while 23% have only evidence on the RNA level and 7% still lack experimental evidence. Analysis of the expression patterns shows few tissue-specific proteins and approximately half of the genes expressed in all the analyzed cells. The status for each gene with regards to protein evidence is visualized in a chromosome-centric manner as part of a new version of the Human Protein Atlas ( www.proteinatlas.org ).

Prognostic implications of estrogen receptor pattern of both tumors in contralateral breast cancer.

Estrogen receptor (ER) status is important for breast cancer survival, it is however unclear how prognosis of contralateral breast cancer (CBC) is affected by ER-status of the two tumors. We conducted a large, population-based study of ER-status of both tumors in CBC patients and its influence on prognosis. The cohort consisted of all women diagnosed with CBC in Stockholm, Sweden during 1976-2005, with information on ER-status from medical records (N = 933). Prognosis was modeled as incidence rates of distant metastasis via Poisson regression. The proportion of CBCs with both cancers of the same ER-status was significantly larger than expected by chance. For synchronous (simultaneous) cancers the prognosis was significantly affected by the combined ER-status of both tumors (p = 0.01). Compared to unilateral breast cancer patients the incidence rate ratio (IRR) for patients with double ER-positive tumors was 1.25 (95 % CI: 0.88-1.76), for ER-discordant tumors 2.19 (95 % CI: 1.18-4.08) and for double ER-negative tumors 3.95 (95 % CI: 1.77-8.81). For metachronous (non-simultaneous) cancers, women with double ER-positive tumors had similarly bad prognosis (IRR = 2.95; 95 % CI: 2.39-3.64) as women with double ER-negative tumors (IRR = 2.88; 95 % CI: 1.83-4.52). Both shorter time span between first and second cancer and endocrine therapy for the first cancer further worsened prognosis of women with double ER-positive metachronous CBC. For synchronous CBC patients, ER-pattern of both tumors is an important prognosticator, while among metachronous CBC patients, double ER-positive tumors confer equally bad prognosis as double ER-negative cancers. Our results indicate that this might be due to endocrine therapy resistance.

Defining the transcriptome and proteome in three functionally different human cell lines.

An essential question in human biology is how cells and tissues differ in gene and protein expression and how these differences delineate specific biological function. Here, we have performed a global analysis of both mRNA and protein levels based on sequence-based transcriptome analysis (RNA-seq), SILAC-based mass spectrometry analysis and antibody-based confocal microscopy. The study was performed in three functionally different human cell lines and based on the global analysis, we estimated the fractions of mRNA and protein that are cell specific or expressed at similar/different levels in the cell lines. A highly ubiquitous RNA expression was found with >60% of the gene products detected in all cells. The changes of mRNA and protein levels in the cell lines using SILAC and RNA ratios show high correlations, even though the genome-wide dynamic range is substantially higher for the proteins as compared with the transcripts. Large general differences in abundance for proteins from various functional classes are observed and, in general, the cell-type specific proteins are low abundant and highly enriched for cell-surface proteins. Thus, this study shows a path to characterize the transcriptome and proteome in human cells from different origins.

Analysis of transcript and protein overlap in a human osteosarcoma cell line.

BACKGROUND: An interesting field of research in genomics and proteomics is to compare the overlap between the transcriptome and the proteome. Recently, the tools to analyse gene and protein expression on a whole-genome scale have been improved, including the availability of the new generation sequencing instruments and high-throughput antibody-based methods to analyze the presence and localization of proteins. In this study, we used massive transcriptome sequencing (RNA-seq) to investigate the transcriptome of a human osteosarcoma cell line and compared the expression levels with in situ protein data obtained in-situ from antibody-based immunohistochemistry (IHC) and immunofluorescence microscopy (IF). RESULTS: A large-scale analysis based on 2749 genes was performed, corresponding to approximately 13% of the protein coding genes in the human genome. We found the presence of both RNA and proteins to a large fraction of the analyzed genes with 60% of the analyzed human genes detected by all three methods. Only 34 genes (1.2%) were not detected on the transcriptional or protein level with any method. Our data suggest that the majority of the human genes are expressed at detectable transcript or protein levels in this cell line. Since the reliability of antibodies depends on possible cross-reactivity, we compared the RNA and protein data using antibodies with different reliability scores based on various criteria, including Western blot analysis. Gene products detected in all three platforms generally have good antibody validation scores, while those detected only by antibodies, but not by RNA sequencing, generally consist of more low-scoring antibodies. CONCLUSION: This suggests that some antibodies are staining the cells in an unspecific manner, and that assessment of transcript presence by RNA-seq can provide guidance for validation of the corresponding antibodies.

Genome-wide profiling of populus small RNAs.

BACKGROUND: Short RNAs, and in particular microRNAs, are important regulators of gene expression both within defined regulatory pathways and at the epigenetic scale. We investigated the short RNA (sRNA) population (18-24 nt) of the transcriptome of green leaves from the sequenced Populus trichocarpa using a concatenation strategy in combination with 454 sequencing. RESULTS: The most abundant size class of sRNAs were 24 nt. Long Terminal Repeats were particularly associated with 24 nt sRNAs. Additionally, some repetitive elements were associated with 22 nt sRNAs. We identified an sRNA hot-spot on chromosome 19, overlapping a region containing both the proposed sex-determining locus and a major cluster of NBS-LRR genes. A number of phased siRNA loci were identified, a subset of which are predicted to target PPR and NBS-LRR disease resistance genes, classes of genes that have been significantly expanded in Populus. Additional loci enriched for sRNA production were identified and characterised. We identified 15 novel predicted microRNAs (miRNAs), including miRNA*sequences, and identified a novel locus that may encode a dual miRNA or a miRNA and short interfering RNAs (siRNAs). CONCLUSIONS: The short RNA population of P. trichocarpa is at least as complex as that of Arabidopsis thaliana. We provide a first genome-wide view of short RNA production for P. trichocarpa and identify new, non-conserved miRNAs.

Automation of cDNA synthesis and labelling improves reproducibility.

BACKGROUND: Several technologies, such as in-depth sequencing and microarrays, enable large-scale interrogation of genomes and transcriptomes. In this study, we asses reproducibility and throughput by moving all laboratory procedures to a robotic workstation, capable of handling superparamagnetic beads. Here, we describe a fully automated procedure for cDNA synthesis and labelling for microarrays, where the purification steps prior to and after labelling are based on precipitation of DNA on carboxylic acid-coated paramagnetic beads. RESULTS: The fully automated procedure allows for samples arrayed on a microtiter plate to be processed in parallel without manual intervention and ensuring high reproducibility. We compare our results to a manual sample preparation procedure and, in addition, use a comprehensive reference dataset to show that the protocol described performs better than similar manual procedures. CONCLUSIONS: We demonstrate, in an automated gene expression microarray experiment, a reduced variance between replicates, resulting in an increase in the statistical power to detect differentially expressed genes, thus allowing smaller differences between samples to be identified. This protocol can with minor modifications be used to create cDNA libraries for other applications such as in-depth analysis using next-generation sequencing technologies.

Development and Validation of a Novel RNA Sequencing-Based Prognostic Score for Acute Myeloid Leukemia.

Background: Recent progress in sequencing technologies allows us to explore comprehensive genomic and transcriptomic information to improve the current European LeukemiaNet (ELN) system of acute myeloid leukemia (AML). Methods: We compared the prognostic value of traditional demographic and cytogenetic risk factors, genomic data in the form of somatic aberrations of 25 AML-relevant genes, and whole-transcriptome expression profiling (RNA sequencing) in 267 intensively treated AML patients (Clinseq-AML). Multivariable penalized Cox models (overall survival [OS]) were developed for each data modality (clinical, genomic, transcriptomic), together with an associated prognostic risk score. Results: Of the three data modalities, transcriptomic data provided the best prognostic value, with an integrated area under the curve (iAUC) of a time-dependent receiver operating characteristic (ROC) curve of 0.73. We developed a prognostic risk score (Clinseq-G) from transcriptomic data, which was validated in the independent The Cancer Genome Atlas AML cohort (RNA sequencing, n = 142, iAUC = 0.73, comparing the high-risk group with the low-risk group, hazard ratio [HR]OS = 2.42, 95% confidence interval [CI] = 1.51 to 3.88). Comparison between Clinseq-G and ELN score iAUC estimates indicated strong evidence in favor of the Clinseq-G model (Bayes factor = 26.78). The proposed model remained statistically significant in multivariable analysis including the ELN and other well-known risk factors (HRos = 2.34, 95% CI = 1.30 to 4.22). We further validated the Clinseq-G model in a second independent data set (n = 458, iAUC = 0.66, adjusted HROS = 2.02, 95% CI = 1.33 to 3.08; adjusted HREFS = 2.10, 95% CI = 1.42 to 3.12). Conclusions: Our results indicate that the Clinseq-G prediction model, based on transcriptomic data from RNA sequencing, outperforms traditional clinical parameters and previously reported models based on genomic biomarkers.

Expression levels of long non-coding RNAs are prognostic for AML outcome.

BACKGROUND: Long non-coding RNA (lncRNA) expression has been implicated in a range of molecular mechanisms that are central in cancer. However, lncRNA expression has not yet been comprehensively characterized in acute myeloid leukemia (AML). Here, we assess to what extent lncRNA expression is prognostic of AML patient overall survival (OS) and determine if there are indications of lncRNA-based molecular subtypes of AML. METHODS: We performed RNA sequencing of 274 intensively treated AML patients in a Swedish cohort and quantified lncRNA expression. Univariate and multivariate time-to-event analysis was applied to determine association between individual lncRNAs with OS. Unsupervised statistical learning was applied to ascertain if lncRNA-based molecular subtypes exist and are prognostic. RESULTS: Thirty-three individual lncRNAs were found to be associated with OS (adjusted p value < 0.05). We established four distinct molecular subtypes based on lncRNA expression using a consensus clustering approach. LncRNA-based subtypes were found to stratify patients into groups with prognostic information (p value < 0.05). Subsequently, lncRNA expression-based subtypes were validated in an independent patient cohort (TCGA-AML). LncRNA subtypes could not be directly explained by any of the recurrent cytogenetic or mutational aberrations, although associations with some of the established genetic and clinical factors were found, including mutations in NPM1, TP53, and FLT3. CONCLUSION: LncRNA expression-based four subtypes, discovered in this study, are reproducible and can effectively stratify AML patients. LncRNA expression profiling can provide valuable information for improved risk stratification of AML patients.

Publisher Correction: Evidence that breast cancer risk at the 2q35 locus is mediated through IGFBP5 regulation.

This corrects the article DOI: 10.1038/ncomms5999.

Molecular Differences between Screen-Detected and Interval Breast Cancers Are Largely Explained by PAM50 Subtypes.

Purpose: Interval breast cancer is of clinical interest, as it exhibits an aggressive phenotype and evades detection by screening mammography. A comprehensive picture of somatic changes that drive tumors to become symptomatic in the screening interval can improve understanding of the biology underlying these aggressive tumors.Experimental Design: Initiated in April 2013, Clinical Sequencing of Cancer in Sweden (Clinseq) is a scientific and clinical platform for the genomic profiling of cancer. The breast cancer pilot study consisted of women diagnosed with breast cancer between 2001 and 2012 in the Stockholm/Gotland regions. A subset of 307 breast tumors was successfully sequenced, of which 113 were screen-detected and 60 were interval cancers. We applied targeted deep sequencing of cancer-related genes; low-pass, whole-genome sequencing; and RNA sequencing technology to characterize somatic differences in the genomic and transcriptomic architecture by interval cancer status. Mammographic density and PAM50 molecular subtypes were considered.Results: In the univariate analyses, TP53, PPP1R3A, and KMT2B were significantly more frequently mutated in interval cancers than in screen-detected cancers. Acquired somatic copy number aberrations with a frequency difference of at least 15% between the two groups included gains in 17q23-q25.3 and losses in 16q24.2. Gene expression analysis identified 447 significantly differentially expressed genes, of which 120 were replicated in an independent microarray dataset. After adjusting for PAM50, most differences were no longer significant.Conclusions: Molecular differences by interval cancer status were observed, but they were largely explained by PAM50 subtypes. This work offers new insights into the biological differences between the two tumor groups. Clin Cancer Res; 23(10); 2584-92. ©2016 AACR.

Study design requirements for RNA sequencing-based breast cancer diagnostics.

Sequencing-based molecular characterization of tumors provides information required for individualized cancer treatment. There are well-defined molecular subtypes of breast cancer that provide improved prognostication compared to routine biomarkers. However, molecular subtyping is not yet implemented in routine breast cancer care. Clinical translation is dependent on subtype prediction models providing high sensitivity and specificity. In this study we evaluate sample size and RNA-sequencing read requirements for breast cancer subtyping to facilitate rational design of translational studies. We applied subsampling to ascertain the effect of training sample size and the number of RNA sequencing reads on classification accuracy of molecular subtype and routine biomarker prediction models (unsupervised and supervised). Subtype classification accuracy improved with increasing sample size up to N = 750 (accuracy = 0.93), although with a modest improvement beyond N = 350 (accuracy = 0.92). Prediction of routine biomarkers achieved accuracy of 0.94 (ER) and 0.92 (Her2) at N = 200. Subtype classification improved with RNA-sequencing library size up to 5 million reads. Development of molecular subtyping models for cancer diagnostics requires well-designed studies. Sample size and the number of RNA sequencing reads directly influence accuracy of molecular subtyping. Results in this study provide key information for rational design of translational studies aiming to bring sequencing-based diagnostics to the clinic.

Sequencing-based breast cancer diagnostics as an alternative to routine biomarkers.

Sequencing-based breast cancer diagnostics have the potential to replace routine biomarkers and provide molecular characterization that enable personalized precision medicine. Here we investigate the concordance between sequencing-based and routine diagnostic biomarkers and to what extent tumor sequencing contributes clinically actionable information. We applied DNA- and RNA-sequencing to characterize tumors from 307 breast cancer patients with replication in up to 739 patients. We developed models to predict status of routine biomarkers (ER, HER2,Ki-67, histological grade) from sequencing data. Non-routine biomarkers, including mutations in BRCA1, BRCA2 and ERBB2(HER2), and additional clinically actionable somatic alterations were also investigated. Concordance with routine diagnostic biomarkers was high for ER status (AUC = 0.95;AUC(replication) = 0.97) and HER2 status (AUC = 0.97;AUC(replication) = 0.92). The transcriptomic grade model enabled classification of histological grade 1 and histological grade 3 tumors with high accuracy (AUC = 0.98;AUC(replication) = 0.94). Clinically actionable mutations in BRCA1, BRCA2 and ERBB2(HER2) were detected in 5.5% of patients, while 53% had genomic alterations matching ongoing or concluded breast cancer studies. Sequencing-based molecular profiling can be applied as an alternative to histopathology to determine ER and HER2 status, in addition to providing improved tumor grading and clinically actionable mutations and molecular subtypes. Our results suggest that sequencing-based breast cancer diagnostics in a near future can replace routine biomarkers.

Gene regulatory mechanisms underpinning prostate cancer susceptibility.

Molecular characterization of genome-wide association study (GWAS) loci can uncover key genes and biological mechanisms underpinning complex traits and diseases. Here we present deep, high-throughput characterization of gene regulatory mechanisms underlying prostate cancer risk loci. Our methodology integrates data from 295 prostate cancer chromatin immunoprecipitation and sequencing experiments with genotype and gene expression data from 602 prostate tumor samples. The analysis identifies new gene regulatory mechanisms affected by risk locus SNPs, including widespread disruption of ternary androgen receptor (AR)-FOXA1 and AR-HOXB13 complexes and competitive binding mechanisms. We identify 57 expression quantitative trait loci at 35 risk loci, which we validate through analysis of allele-specific expression. We further validate predicted regulatory SNPs and target genes in prostate cancer cell line models. Finally, our integrated analysis can be accessed through an interactive visualization tool. This analysis elucidates how genome sequence variation affects disease predisposition via gene regulatory mechanisms and identifies relevant genes for downstream biomarker and drug development.

Reprogramming Tumor-Associated Macrophages by Antibody Targeting Inhibits Cancer Progression and Metastasis.

Tumors are composed of multiple cell types besides the tumor cells themselves, including innate immune cells such as macrophages. Tumor-associated macrophages (TAMs) are a heterogeneous population of myeloid cells present in the tumor microenvironment (TME). Here, they contribute to immunosuppression, enabling the establishment and persistence of solid tumors as well as metastatic dissemination. We have found that the pattern recognition scavenger receptor MARCO defines a subtype of suppressive TAMs and is linked to clinical outcome. An anti-MARCO monoclonal antibody was developed, which induces anti-tumor activity in breast and colon carcinoma, as well as in melanoma models through reprogramming TAM populations to a pro-inflammatory phenotype and increasing tumor immunogenicity. This anti-tumor activity is dependent on the inhibitory Fc-receptor, FcγRIIB, and also enhances the efficacy of checkpoint therapy. These results demonstrate that immunotherapies using antibodies designed to modify myeloid cells of the TME represent a promising mode of cancer treatment.

Evaluation of exome sequencing to estimate tumor burden in plasma.

Accurate estimation of systemic tumor load from the blood of cancer patients has enormous potential. One avenue is to measure the presence of cell-free circulating tumor DNA in plasma. Various approaches have been investigated, predominantly covering hotspot mutations or customized, patient-specific assays. Therefore, we investigated the utility of using exome sequencing to monitor circulating tumor DNA levels through the detection of single nucleotide variants in plasma. Two technologies, claiming to offer efficient library preparation from nanogram levels of DNA, were evaluated. This allowed us to estimate the proportion of starting molecules measurable by sequence capture (<5%). As cell-free DNA is highly fragmented, we designed and provide software for efficient identification of PCR duplicates in single-end libraries with a varying size distribution. On average, this improved sequence coverage by 38% in comparison to standard tools. By exploiting the redundant information in PCR-duplicates the background noise was reduced to ∼1/35,000. By applying our optimized analysis pipeline to a simulation analysis, we determined the current sensitivity limit to ∼1/2400, starting with 30 ng of cell-free DNA. Subsequently, circulating tumor DNA levels were assessed in seven breast- and one prostate cancer patient. One patient carried detectable levels of circulating tumor DNA, as verified by break-point specific PCR. These results demonstrate exome sequencing on cell-free DNA to be a powerful tool for disease monitoring of metastatic cancers. To enable a broad implementation in the diagnostic settings, the efficiency limitations of sequence capture and the inherent noise levels of the Illumina sequencing technology must be further improved.