Tumor-associated copy number changes in the circulation of patients with prostate cancer identified through whole-genome sequencing
- Ellen Heitzer†1,
- Peter Ulz†1,
- Jelena Belic1,
- Stefan Gutschi2,
- Franz Quehenberger3,
- Katja Fischereder2,
- Theresa Benezeder1,
- Martina Auer1,
- Carina Pischler1,
- Sebastian Mannweiler4,
- Martin Pichler5,
- Florian Eisner5,
- Martin Haeusler6,
- Sabine Riethdorf7,
- Klaus Pantel7,
- Hellmut Samonigg5,
- Gerald Hoefler4,
- Herbert Augustin2,
- Jochen B Geigl1Email author and
- Michael R Speicher1Email author
© Heitzer et al.; licensee BioMed Central Ltd. 2013
Received: 9 February 2013
Accepted: 5 April 2013
Published: 5 April 2013
Patients with prostate cancer may present with metastatic or recurrent disease despite initial curative treatment. The propensity of metastatic prostate cancer to spread to the bone has limited repeated sampling of tumor deposits. Hence, considerably less is understood about this lethal metastatic disease, as it is not commonly studied. Here we explored whole-genome sequencing of plasma DNA to scan the tumor genomes of these patients non-invasively.
We wanted to make whole-genome analysis from plasma DNA amenable to clinical routine applications and developed an approach based on a benchtop high-throughput platform, that is, Illuminas MiSeq instrument. We performed whole-genome sequencing from plasma at a shallow sequencing depth to establish a genome-wide copy number profile of the tumor at low costs within 2 days. In parallel, we sequenced a panel of 55 high-interest genes and 38 introns with frequent fusion breakpoints such as the TMPRSS2-ERG fusion with high coverage. After intensive testing of our approach with samples from 25 individuals without cancer we analyzed 13 plasma samples derived from five patients with castration resistant (CRPC) and four patients with castration sensitive prostate cancer (CSPC).
The genome-wide profiling in the plasma of our patients revealed multiple copy number aberrations including those previously reported in prostate tumors, such as losses in 8p and gains in 8q. High-level copy number gains in the AR locus were observed in patients with CRPC but not with CSPC disease. We identified the TMPRSS2-ERG rearrangement associated 3-Mbp deletion on chromosome 21 and found corresponding fusion plasma fragments in these cases. In an index case multiregional sequencing of the primary tumor identified different copy number changes in each sector, suggesting multifocal disease. Our plasma analyses of this index case, performed 13 years after resection of the primary tumor, revealed novel chromosomal rearrangements, which were stable in serial plasma analyses over a 9-month period, which is consistent with the presence of one metastatic clone.
The genomic landscape of prostate cancer can be established by non-invasive means from plasma DNA. Our approach provides specific genomic signatures within 2 days which may therefore serve as 'liquid biopsy'.
Prostate cancer is the most common malignancy in men. In Europe each year an estimated number of 2.6 million new cases is diagnosed . The wide application of PSA testing has resulted in a shift towards diagnosis at an early stage so that many patients do not need treatment or are cured by radical surgery . However, patients still present with metastatic or recurrent disease despite initial curative treatment . In these cases prostate-cancer progression can be inhibited by androgen-deprivation therapy (ADT) for up to several years. However, disease progression is invariably observed with tumor cells resuming proliferation despite continued treatment (termed castration-resistant prostate cancer or CRPC) . CRPC is a strikingly heterogeneous disease and the overall survival can be extremely variable . Scarcity of predictive and prognostic markers underlines the growing need for a better understanding of the molecular makeup of these lethal tumors.
However, acquiring tumor tissue from patients with metastatic prostate cancer often represents a challenge. Due to the propensity of metastatic prostate cancer to spread to bone biopsies can be technically challenging and limit repeated sampling of tumor deposits. As a consequence, considerably less is understood about the later acquired genetic alterations that emerge in the context of the selection pressure of an androgen-deprived milieu .
Consistent and frequent findings from recent genomic profiling studies in clinical metastatic prostate tumors include the TMPRSS2-ERG fusion in approximately 50%, 8p loss in approximately 30% to 50%, 8q gain in approximately 20% to 40% of cases, and the androgen receptor (AR) amplification in approximately 33% of CRPC cases [7–10]. Several whole-exome or whole-genome sequencing studies consistently reported low overall mutation rates even in heavily treated CRPCs [9–14].
The difficulties in acquiring tumor tissue can partly be addressed by elaborate procedures such as rapid autopsy programs to obtain high-quality metastatic tissue for analysis . However, this material can naturally only be used for research purposes, but not for biomarker detection for individualized treatment decisions. This makes blood-based assays crucially important to individualize management of prostate cancer . Profiling of blood offers several practical advantages, including the minimally invasive nature of sample acquisition, relative ease of standardization of sampling protocols, and the ability to obtain repeated samples over time. For example, the presence of circulating tumor cells (CTCs) in peripheral blood is a prognostic biomarker and a measure of therapeutic response in patients with prostate cancer [17–20]. Novel microfluidic devices enhance CTC capture [21–23] and allow to establish a non-invasive measure of intratumoral AR signaling before and after hormonal therapy . Furthermore, prospective studies have demonstrated that mRNA expression signatures from whole blood can be used to stratify patients with CRPC into high- and low-risk groups [25, 26].
Another option represents the analysis of plasma DNA . One approach is the identification of known alterations previously found in the resected tumors from the same patients in plasma DNA for monitoring purposes [28, 29]. Furthermore, recurrent mutations can be identified in plasma DNA in a subset of patients with cancer [30–32]. Given that chromosomal copy number changes occur frequently in human cancer, we developed an approach allowing the mapping of tumor-specific copy number changes from plasma DNA employing array-CGH . At the same time, massively parallel sequencing of plasma DNA from the maternal circulation is emerging to a clinical tool for the routine detection of fetal aneuploidy [34–36]. Using essentially the same approach, that is, next-generation sequencing from plasma, the detection of chromosomal alterations in the circulation of three patients with hepatocellular carcinoma and one patient with both breast and ovarian cancer  and from 10 patients with colorectal and breast cancer  was reported.
However, the costs of the aforementioned plasma sequencing studies necessary for detection of rearrangements were prohibitive for routine clinical implementation [37, 38]. In addition, these approaches are very time-consuming. Previously it had been shown that whole-genome sequencing with a shallow sequencing depth of about 0.1x is sufficient for a robust and reliable analysis of copy number changes from single cells . Hence, we developed a different whole-genome plasma sequencing approach employing a benchtop high-throughput sequencing instrument, that is, the Illumina MiSeq, which is based on the existing Solexa sequencing-by-synthesis chemistry, but has dramatically reduced run times compared to the Illumina HiSeq . Using this instrument we performed whole-genome sequencing from plasma DNA and measured copy number from sequence read depth. We refer to this approach as plasma-Seq. Furthermore, we enriched 1.3 Mbp consisting of exonic sequences of 55 high-interest cancer genes and 38 introns of genes, where fusion breakpoints have been described and subjected the DNA to next-generation sequencing at high coverage (approximately 50x). Here we present the implementation of our approach with 25 plasma samples from individuals without cancer and results obtained with whole genome sequencing of 13 plasma DNA samples derived from nine patients (five CRPC, four CSPC) with prostate cancer.
Patient eligibility criteria
This study was conducted among men with prostate cancer (Clinical data in Additional file 1, Table S1) who met the following criteria: histologically-proven, based on a biopsy, metastasized prostate cancer. We distinguished between CRPC and CSPC based on the guidelines on prostate cancer from the European Association of Urology , that is: 1, castrate serum levels of testosterone (testosterone <50 ng/dL or <1.7 nmol/L); 2, three consecutive rises of PSA, 1 week apart, resulting in two 50% increases over the nadir, with a PSA >2 ng/mL; 3, anti-androgen withdrawal for at least 4 weeks for flutamide and for at least 6 weeks for bicalutamide; 4, PSA progression, despite consecutive hormonal manipulations. Furthermore, we focused on patients who had ≥5 CTCs per 7.5 mL  and/or a biphasic plasma DNA size distribution as described previously by us .
The study was approved by the ethics committee of the Medical University of Graz (approval numbers 21-228 ex 09/10, prostate cancer, and 23-250 ex 10/11, prenatal plasma DNA analyses), conducted according to the Declaration of Helsinki, and written informed consent was obtained from all patients and healthy blood donors. Blood from prostate cancer patients and from male controls without malignant disease was obtained from the Department of Urology or the Division of Clinical Oncology, Department of Internal Medicine, at the Medical University of Graz. From prostate cancer patients we obtained a buccal swab in addition. Blood samples from pregnant females and from female controls without malignant disease were collected at the Department of Obstetrics and Gynecology, Medical University of Graz. The blood samples from the pregnant females were taken prior to an invasive prenatal diagnostic procedure.
Plasma DNA preparation
Plasma DNA was prepared using the QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany) as previously described . Samples selected for sequence library construction were analyzed by using the Bioanalyzer instrument (Agilent Technologies, Santa Clara, CA, USA) to observe the plasma DNA size distribution. In this study we included samples with a biphasic plasma DNA size distribution as previously described .
Enumeration of CTCs
We performed CTC enumeration using the automated and FDA approved CellSearch assay. Blood samples (7.5 mL each) were collected into CellSave tubes (Veridex, Raritan, NJ, USA). The Epithelial Cell Kit (Veridex) was applied for CTC enrichment and enumeration with the CellSearch system as described previously [42, 43].
Array-CGH was carried out using a genome-wide oligonucleotide microarray platform (Human genome CGH 60K microarray kit, Agilent Technologies, Santa Clara, CA, USA), following the manufacturer's instructions (protocol version 6.0) as described . Evaluation was done based on our previously published algorithm [33, 44, 45].
HT29 dilution series
Sensitivity of our plasma-Seq approach was determined using serial dilutions of DNA from HT29 cell line (50%, 20%, 15%, 10%, 5%, 1%, and 0%) in the background of normal DNA (Human Genomic DNA: Female; Promega, Fitchburg, WI, USA). Since quantification using absorption or fluorescence absorption is often not reliable we used quantitative PCR to determine the amount of amplifiable DNA and normalized the samples to a standard concentration using the Type-it CNV SYBR Green PCR Kits (Qiagen, Hilden, Germany). Dilution samples were then fragmented using the Covaris S220 System (Covaris, Woburn, MA, USA) to a maximum of 150-250 bp and 10 ng of each dilution were used for library preparation to simulate plasma DNA condition.
Shotgun libraries were prepared using the TruSeq DNA LT Sample preparation Kit (Illumina, San Diego, CA, USA) following the manufacturer´s instructions with three exceptions. First, due to limited amounts of plasma DNA samples we used 5-10 ng of input DNA. Second, we omitted the fragmentation step since the size distribution of the plasma DNA samples was analyzed on a Bioanalyzer High Sensitivity Chip (Agilent Technologies, Santa Clara, CA, USA) and all samples showed an enrichment of fragments in the range of 160 to 340 bp. Third, for selective amplification of the library fragments that have adapter molecules on both ends we used 20-25 PCR cycles. Four libraries were pooled equimolarily and sequenced on an Illumina MiSeq (Illumina, San Diego, CA, USA).
The MiSeq instrument was prepared following routine procedures. The run was initiated for 1x150 bases plus 1x25 bases of SBS sequencing, including on-board clustering and paired-end preparation, the sequencing of the respective barcode indices and analysis. On the completion of the run, data were base called and demultiplexed on the instrument (provided as Illumina FASTQ 1.8 files, Phred+33 encoding). FASTQ format files in Illumina 1.8 format were considered for downstream analysis.
Calculation of segments with identical log2 ratio values
We employed a previously published algorithm  to create a reference sequence. The pseudo-autosomal region (PAR) on the Y chromosome was masked and the mappability of each genomic position examined by creating virtual 150 bp reads for each position in the PAR-masked genome. Virtual sequences were mapped to the PAR-masked genome and mappable reads were extracted. Fifty thousand genomic windows were created (mean size, 56,344 bp) each having the same amount of mappable positions.
Low-coverage whole-genome sequencing reads were mapped to the PAR-masked genome and reads in different windows were counted and normalized by the total amount of reads. We further normalized read counts according to the GC content using LOWESS-statistics. In order to avoid position effects we normalized the sequencing data with GC-normalized read counts of plasma DNA of our healthy controls and calculated log2 ratios.
Calculation of z-scores for specific regions
Calculation of genome-wide z-scores
In order to establish a genome-wide z-score to detect aberrant genomic content in plasma, we divided the genome into equally-sized regions of 1 Mbp length and calculated z-scores therein.
Under the condition that all ratios were drawn from the same normal distribution, z-scores are distributed proportionally to Student's t-distribution with n-1 degrees of freedom. For controls, z-scores were calculated using cross-validation. In brief, z-score calculation of one control is based on means and standard deviation of the remaining controls. This prevents controls from serving as their own controls.
The variance of these cross-validated z-scores of controls is slightly higher than the variance of z-scores of tumor patients. Thus ROC performance is underestimated. This was confirmed in the simulation experiment described below.
Genome-wide z-scores were calculated from S-scores. Other methods of aggregation of z-score information, such as sums of absolute values or PA scores , performed poorer and were therefore not considered. Per window z-scores were clustered hierarchically by the hclust function of R using Manhattan distance that summed up the distance of each window.
In order to validate the diagnostic performance of the genome-wide z-score in silico, artificial cases and controls were simulated from mean and standard deviations of ratios from 10 healthy controls according to a normal distribution. Simulated tumor cases were obtained through multiplication of the mean by the empirical copy number ratio of 204 prostate cancer cases . Segmented DNA-copy-number data were obtained via the cBio Cancer Genomics Portal .
We performed ROC analyses of 500 simulated controls and 102 published prostate tumor data and their respective dilutions using the pROC R-package . The prostate tumor data were derived from a previously published dataset  and the 102 cases were selected based on their copy number profiles.
Gene-Breakpoint Panel: target enrichment of cancer genes, alignment and SNP-calling, SNP-calling results
We enriched 1.3 Mbp of seven plasma DNAs (four CRPC cases, CRPC1-3 and CRPC5; three CSPC cases, CSPC1-2 and CSPC4) including exonic sequences of 55 cancer genes and 38 introns of 18 genes, where fusion breakpoints have been described using Sure Select Custom DNA Kit (Agilent, Santa Clara, CA, USA) following the manufacturer's recommendations. Since we had very low amounts of input DNA we increased the number of cycles in the enrichment PCR to 20. Six libraries were pooled equimolarily and sequenced on an Illumina MiSeq (Illumina, San Diego, CA, USA).
We generated a mean of 7.78 million reads (range, 3.62-14.96 million), 150 bp paired-end reads on an Illumina MiSeq (Illumina, San Diego, CA, USA). Sequences were aligned using BWA  and duplicates were marked using picard . We subsequently performed realigning around known indels and applied the Unified Genotyper SNP-calling software provided by the GATK .
We further annotated resulting SNPs by employing annovar  and reduced the SNP call set by removing synonymous variants, variants in segmental duplications and variants listed in the 1000 Genome Project  and Exome sequencing (Project Exome Variant Server, NHLBI Exome Sequencing Project (ESP), Seattle, WA)  with allele frequency >0.01.
We set very stringent criteria to reduce false positives according to previously published values : a mutation had to be absent from the constitutional DNA sequencing and the sequencing depth for the particular nucleotide position had to be >20-fold. Furthermore, all putative mutations or breakpoint spanning regions were verified by Sanger sequencing.
Since plasma DNA is fragmented the read pair method is not suitable for identification of structural rearrangements  and therefore we performed split-read analysis of 150 bp reads. We used the first and the last 60 bp of each read (leaving a gap of 30 bp) and mapped these independently. We further analyzed discordantly mapped split-reads by focusing on targeted regions and filtering out split-reads mapping within repetitive regions and alignments having a low mapping quality (<25). Reads where discordantly mapped reads were found were aligned to the human genome using BLAT  to further specify putative breakpoints.
All sequencing raw data were deposited at the European Genome-phenome Archive (EGA) , which is hosted by the EBI, under accession numbers EGAS00001000451 (Plasma-Seq) and EGAS00001000453 (Gene-Breakpoint Panel).
Implementation of our approach
Previously, we demonstrated that tumor-specific, somatic chromosomal alterations can be detected from plasma of patients with cancer using array-CGH . In order to extend our method to a next-generation sequencing-based approach, that is, plasma-Seq, on a benchtop Illumina MiSeq instrument, we first analyzed plasma DNA from 10 men (M1 to M10) and nine women (F1 to F9) without malignant disease. On average we obtained 3.3 million reads per sample (range, 1.9-5.8 million; see Additional file 1, Table S2) and applied a number of filtering steps to remove sources of variation and to remove known GC bias effects [62–64] (for details see Material and Methods).
In addition, we calculated 'segmental z-scores' where the z-scores are not calculated for 1-Mbp windows but for chromosomal segments with identical copy number. In order to determine such segments we employed an algorithm for the assignment of segments with identical log2 ratios [39, 46] (Material and Methods) and calculated a z-score for each of these segments (hence, 'segmental z-scores'). As sequencing analyses of chromosome content in the maternal circulation are now frequently being used for detection of fetal aneuploidy [34, 36] and as our mean sequencing depth is lower compared to previous studies, we wanted to test whether our approach would be feasible for this application. To this end we obtained two plasma samples each of pregnancies with euploid and trisomy 21 fetuses and one each of pregnancies with trisomies of chromosomes 13 and 18, respectively. In the trisomy cases the respective chromosomes were identified as segments with elevated log2 ratios and accordingly also increased z-scores (Additional file 2).
Sensitivity and specificity of our approach
We wanted to gain insight into the sensitivity of our approach to detect tumor-derived sequences in a patient's plasma. To this end we calculated a genome-wide z-score for each sample (Material and Methods). The main purpose of the genome-wide z-score is to distinguish between aneuploid and euploid plasma samples. The genome-wide z-score from the plasma of male individuals ranged from -1.10 to 2.78 and for female individuals from -0.48 to 2.64. We performed receiver operating characteristic (ROC) analyses of simulated next-generation sequencing data from 102 published prostate cancer data and 500 simulated controls based on the data from our healthy individuals. Using the equivalent of one-quarter MiSeq run, these analyses suggested that using the genome-wide z-score tumor DNA concentrations at levels ≥10% could be detected in the circulation of patients with prostate cancers with a sensitivity of >80% and specificity of >80% (Figure 1c).
To test these estimates with actual data we fragmented DNA from the colorectal cancer cell line HT29 to sizes of approximately 150-250 bp to reflect the degree of fragmented DNA in plasma and performed serial dilution experiments with the fragmented DNA (that is, 50%, 20%, 15%, 10%, 5%, 1%, and 0%). We established the copy-number status of this cell line with undiluted, that is, 100%, DNA using both array-CGH and our next-generation sequencing approach (Additional file 3) and confirmed previously reported copy number changes [65, 66]. Calculating the genome-wide z-score for each dilution we noted its expected decrease with increasing dilution. Whereas the genome-wide z- score was 429.74 for undiluted HT29 DNA, it decreased to 7.75 for 1% (Additional file 1, Table S2). Furthermore, when we performed hierarchical cluster analysis the female controls were separated from the various HT29 dilutions, further confirming that our approach may indicate aneuploidy in the presence of 1% circulating tumor DNA (Figure 1d).
Plasma analysis from patients with cancer
Plasma-Seq and GB-panel of patients with prostate cancer
We then obtained 13 plasma samples from nine patients with metastatic prostate cancer (five with castration-resistant disease, CRPC1 to CRPC5, and four with castration-sensitive disease, CSPC1 to CSPC4. Furthermore, from each of patients CRPC1 and CSPC1 we obtained three samples at different time points (Clinical data in Additional file 1, Table S1). Applying plasma-Seq, we obtained on average 3.2 million reads (range, 1.1 (CSPC4) to 5.2 (CRPC5) million reads) for the plasma samples from patients with prostate cancer per sample (see Additional file 1, Table S2).
Applying the GB-panel, we achieved on average a coverage of ≥50x for 71.8% of target sequence (range, 67.3% (CSPC4) to 77.6% (CSPC2)) (see Additional file 1, Table S4). Using very stringent conditions (see Material & Methods) the GB-panel allowed us to identify 12 mutations in all seven patients for which the analyses were performed (that is, CRPC1-3, CRPC5, CSPC1-2, and CSPC4). Sanger sequencing confirmed the presence of five of these mutations in both plasma and the respective constitutional DNA, whereas seven mutations were only confirmed in plasma but not in the constitutional DNA. The latter mutations, which were observed in five patients (that is, CRPC2-3, CRPC5, CSPC2, CSPC4), are likely somatic mutations and occurred in genes previously implicated in prostate cancer tumorigenesis, such as TP53, BRCA1, BRCA2, and MLL3 (see Additional file 1, Table S4). We used these somatic mutations for ultra-deep sequencing with an average coverage of 362,016 (range, 307,592 to 485,467) to estimate the tumor fraction. Using these estimates the tumor fraction was lowest in CSPC4 with 30.75% and highest in CRPC5 with 54.49%.
Plasma-Seq from these patients exhibited a wide range of copy number aberrations indicative of malignant origin, including those that have been previously reported in prostate tumors. For example, the three CRPC patients (that is, CRPC2-3, CRPC5) had high-level gains in a region on chromosome × including the AR locus. Over-representation of 8q regions was observed in all five patients and loss of 8p regions in three patients (CRPC5, CSPC2, and CSPC4) (Figure 3c).
As control we performed array-CGH analyses of all plasma cases as described  in parallel (see Additional file 4). These array-CGH profiles had a great concordance with those obtained with plasma-Seq.
TMPRSS2-ERG fusion mapping
The fusion through deletion TMPRSS2-ERG rearrangement results in a well-defined 3-Mbp interstitial deletion on chromosome 21 [69, 70] and occurs in approximately 50% of prostate cancer cases . We tested whether our approach would allow distinguishing TMPRSS2-ERG-positive from TMPRSS2-ERG-negative prostate cancers.
Analyses of serial plasma samples
Plasma samples were taken at three different time points over a 9-month period (we refer to them in addition to CRPC1 as CRPC1_2 and CRPC1_3). At the time of our plasma collections the patient was castration resistant and had stable disease under ongoing ADT and chemotherapy. Plasma-Seq identified again multiple prostate cancer-associated chromosomal alterations, such as 8p loss, gain of 8q regions, the 3-Mbp TMPRSS2-ERG deletion on chromosome 21, and AR amplification (Figures 4 and 5). Thus, plasma-Seq identified multiple rearrangements, that is, the TMPRSS2-ERG deletion on chromosome 21, which had not been present in the primary tumor. Furthermore, plasma-Seq yielded remarkably similar results in our three analyses over the 9-month period (Figure 5), which is in agreement with the clinically stable disease and suggests the presence of one dominant clone releasing DNA into the circulation. This is consistent with the proposed monoclonal origin of metastatic prostate cancer . Hierarchical clustering confirmed the concordance between the three plasma-Seq copy number profiles and the tremendous differences to the various sectors of the primary tumor (Figure 5).
The high similarity of copy number changes at various time points is another confirmation of the high reliability and robustness of our approach.
Evaluation of copy number changes of prostate cancer genes
The evaluation of 1-Mbp or segmental z-scores each involves relatively large regions. We wanted to test whether z-scores can also be calculated for much smaller regions, that is specific genes, and calculated gene-specific z-scores (see Material & Methods).
For example, in prostate cancer, one of the most interesting regions is the AR-locus on chromosome Xq12, which is amplified in approximately 33% of patients with CRPC . As expected, none of the male healthy controls had an amplification of AR, whereas AR amplification was present in four of the five CRPC cases. In order to validate the plasma-Seq gene-specific copy number estimates with another approach we selected a subset of samples (CRPC1, CRPC2, CRPC5, CSPC1, CSPC1_2, and CSPC2) for validation of the AR copy-number status with qPCR. In fact, we observed a very close correlation between the plasma-Seq and the qPCR values (see Additional file 6). Interestingly, CRPC1 had only a duplication of the AR region and the AR copy number did not change over our observation period of 9 months, which was consistent with the clinically stable disease. One of the CSPC cases, CSPC4, had a slightly increased AR ratio (ratio, 1.46; z-score, 4.60). Whether such a value may indicate the beginning of ADT resistance remains presently unclear, as sufficient follow-up data were not available.
We also tested our approach for some other genes, which have frequently been implicated in prostate cancer. For example, evidence for cooperation between AR and NCOA2 amplifications on 8q13.3 in early prostate cancer was reported . However, alternatively it was suggested that tumors first acquire NCOA2 amplification along with broad amplifications on chromosome 8q . Our gene-specific z-score identified NCOA2 gene amplifications in five patients (CRPC1, CRPC5, CSPC1-3), thus, two CRPC and three CSPC cases, which may support the notion that NCOA2 amplifications may occur prior to AR amplification .
Loss of PTEN on 10q23.31 occurs in approximately 40% of prostate cancers [9, 74]. We observed PTEN loss in five patients (CRPC3-5, CSPC1, and CSPC3); that is, in three CRPC and two CSPC cases. The AKT-inactivating phosphatase PHLPP1 on 18q21.33 has recently been identified as a prostate tumor suppressor . We found that this gene was lost in four patients (CRPC1, CRPC3, CSPC1-2); that is, in two CRPC and two CSPC cases. Furthermore, it has recently been reported that the TMPRSS2-ERG fusion is associated with a deletion at chromosome 3p14 that includes the FOXP1 gene . In fact, we observed loss of this region in five of our patients (CRPC1-2, CRPC4, CSPC1, CSPC4) and four of these patients (CRPC1-2, CSPC1, CSPC4) did indeed have the TMPRSS2-ERG fusion, confirming the association between these two loci.
In summary, our results suggest that gene-specific information can be derived from plasma-Seq, which may facilitate the evaluation of pathways potentially comprised in prostate cancer.
This study represents the first whole-genome sequencing analysis from plasma DNA of patients with prostate cancer. Usually the identification of tumor genotypes that inform selection of targeted therapies is performed on the initial diagnostic specimen. However, these may not be readily available or in case of fine needle aspirates not sufficient for molecular analyses, as was the case for our patients who presented with metastatic disease. The only exception in our cohort was CRPC1, who had recurrent disease many years after initial operative treatment. We could demonstrate that the initial primary tumor specimen represented multifocal disease and none of the analyzed sectors was representative of the metastatic clone, which arose 13 years later. Thus, molecular analysis of plasma may provide a non-invasive approach for tumor cell genotyping, which can easily be repeated during the course of therapy.
Multiple lines of evidence support the copy number changes observed. First, the observation of known prostate cancer alterations in our dataset indicates successful performance of our assay. Second, our previously published array-based plasma method  was applied in parallel to confirm the copy number aberrations observed with plasma-Seq. Third, we identified the well characterized 3-Mbp interstitial 21q22.2-3 deletion spanning ERG and TMPRSS2 on chromosome 21 [69, 70] and confirmed its presence with our GB-panel and Sanger sequencing. Fourth, for two of our patients we were able to repeat our analysis at different time points. These repeated analyses revealed a high degree of similarity among samples from the same patient. The shared copy number aberrations were indicative of common lineage, which is consistent with the view that metastases in this disease are of monoclonal origin . Finally, implementation of our approach with 19 samples from individuals without cancer and five plasma samples from pregnant females with aneuploidy fetuses further confirmed the reliability and robustness of our approach.
Tests for sensitivity and specificity of our approach suggested that tumor DNA concentrations at levels ≥10% can be detected with a sensitivity of >80% and specificity of 80%. Furthermore, our simulations and HT29 dilution experiments suggested that the genome-wide z-score detects aneuploidy even at tumor DNA concentrations of only 1%. In general, the resolution of non-invasive tumor genome-wide scans from plasma is limited by the depth of the sequencing and the percentage of tumor fragments in the plasma. Therefore, previously published similar studies [37, 38] employed high-throughput sequencing platforms tailored chiefly toward large-scale applications. As a consequence, footprints, workflows, reagent costs, and run times are poorly matched to the needs of small laboratories and furthermore, the cost of the sequencing necessary for detection of rearrangements at this level is prohibitive for routine clinical implementation . In contrast, advantages of a benchtop high-throughput sequencing instrument include the speed of analyses and the reduced costs. A MiSeq run produces a throughput of 1.6 Gbp with a read length of 150 bp . As whole-genome sequencing with a 0.1x coverage was reported to yield robust and reliable copy-number measurements from single cells , we tested such a sequencing approach for our plasma analyses. Accordingly, we found that the characteristics of the MiSeq are sufficient for our plasma-Seq purposes. Especially attractive features of this strategy include the speed (library prep, approximately 24 h; sequencing of 150 bp single reads, approximately 12 h; identification of segments with identical log2 ratios, approximately 2 h; calculation of z-scores, 30 min) and the costs (approximately €300) with which the aneuploidy scoring by plasma-Seq can be performed. In contrast, completion of the GB-panel analysis, done at 50x coverage, will normally require at least 7 days (library prep, approximately 24 h; targeted enrichment, approximately 48-72 h; sequencing 150 bp paired end, approximately 26 h; evaluation and SNP calling, several hours) not including verifications of mutations by Sanger sequencing or estimation of the fractional load of tumor fragments by deep sequencing.
A disadvantage of low coverage whole-genome sequencing is that structural inter- and intrachromosomal rearrangements cannot be identified with high confidence. This is because plasma DNA fragments, whose paired-end reads map to different chromosomes or to the same chromosome but at large distances (several kb) apart, will likely not be detected in multiple reads. Another disadvantage is the reduced resolution for identification of mutations. However, several large scale whole-exome or whole-genome sequencing studies consistently reported low overall mutation rates even in heavily treated CRPCs ranging from 0.9 to 2.00 mutations per mega base [9, 11–14]. These studies confirmed that the most commonly mutated gene was AR, however no single gene other than AR had frequent mutations and even common, broadly mutated oncogenes such as PIK3CA, KRAS, and BRAF are not commonly mutated in prostate cancer . We addressed both issues, structural rearrangements and mutations, with a focused sequencing approach with higher coverage. Focused sequencing, such as our GB-panel, with tailored design and analytical prioritization strategies may represent an attractive alternative to large-scale whole-genome sequencing in terms of speed and costs. Such a focused approach is flexible and can easily be adapted if new, important genes or regions evolve from large-scale sequencing projects.
Another potential short-coming is that we do not know whether the changes observed in the plasma are related to the primary tumor or to any of the metastatic sites. In fact, it is currently unknown whether all tumor cells contribute to the plasma DNA equally and which factors influence the release of tumor DNA into the circulation. Further studies are needed to determine whether changes observed by plasma-Seq represent an average of the DNA alterations from all malignant sites or whether they show characteristic changes of the dominant tumor cell clone at the time of the blood collection.
At present we do not know how our plasma DNA signatures perform compared with other emerging candidate markers, for example, CTC analysis . However, our approach circumvents an inherent limitation of all published CTC-based studies, that is, it is not focused on EpCAM-positive CTCs. Furthermore, plasma isolation does not necessitate special equipment as usually required for CTC isolation [21–23]. As we already have plasma-Seq data from patients with colon and breast cancer our method may also be applicable to other tumor types.
Whether these blood copy-number signatures will be true game changers for the management of prostate cancer has to be further evaluated. Drug development for castration-resistant prostate cancer is an area of intensive research and several new agents are currently being tested in phase 3 clinical trials. Interrogation of the genomic signature may reveal whether those targeted therapies are effectively hitting their target in vivo, thus providing information that may be useful in guiding therapeutic decisions.
Our strategy may contribute to a better definition of the evolution towards a castration-resistant disease and could potentially aid in identifying patients more or less likely respond to AR-targeted therapies. The simplicity and the costs of our test are attractive and might ease the clinical translation. However, the extent to which these signatures contribute independent prognostic or predictive value beyond clinicopathological variables must be explored in more depth.
Androgen deprivation therapy
Basic Local Alignment Search Tool (BLAST)-like alignment tool
Comparative genomic hybridization
Castration-resistant prostate cancer
Castration-sensitive prostate cancer
Circulating tumor cell
European Bioinformatics Institute
European Genome-Phenome Archive
Exome Sequencing Project
Genome Analysis Toolkit
We thank Mag. Maria Langer-Winter for editing the manuscript. Funding was provided by the Austrian Science Fund (FWF) (grant#: P20338, P23284, and W 1226-B18, DKplus Metabolic and Cardiovascular Disease), the Oesterreichische Nationalbank (15093), and by the COMET Center ONCOTYROL.
- Bray F, Sankila R, Ferlay J, Parkin DM: Estimates of cancer incidence and mortality in Europe in 1995. Eur J Cancer. 2002, 38: 99-166.View ArticlePubMedGoogle Scholar
- Albertsen PC: Treatment of localized prostate cancer: when is active surveillance appropriate?. Nat Rev Clin Oncol. 2010, 7: 394-400.View ArticlePubMedGoogle Scholar
- Shariat SF, Kattan MW, Vickers AJ, Karakiewicz PI, Scardino PT: Critical review of prostate cancer predictive tools. Future Oncol. 2009, 5: 1555-1584.View ArticlePubMedPubMed CentralGoogle Scholar
- Scher HI, Sawyers CL: Biology of progressive, castration-resistant prostate cancer: directed therapies targeting the androgen-receptor signaling axis. J Clin Oncol. 2005, 23: 8253-8261.View ArticlePubMedGoogle Scholar
- Scher HI, Morris MJ, Basch E, Heller G: End points and outcomes in castration-resistant prostate cancer: from clinical trials to clinical practice. J Clin Oncology. 2011, 29: 3695-3704.View ArticleGoogle Scholar
- Friedlander TW, Roy R, Tomlins SA, Ngo VT, Kobayashi Y, Azameera A, Rubin MA, Pienta KJ, Chinnaiyan A, Ittmann MM, Ryan CJ, Paris PL: Common structural and epigenetic changes in the genome of castration-resistant prostate cancer. Cancer Res. 2012, 72: 616-625.View ArticlePubMedGoogle Scholar
- Holcomb IN, Young JM, Coleman IM, Salari K, Grove DI, Hsu L, True LD, Roudier MP, Morrissey CM, Higano CS, Nelson PS, Vessella RL, Trask BJ: Comparative analyses of chromosome alterations in soft-tissue metastases within and across patients with castration-resistant prostate cancer. Cancer Res. 2009, 69: 7793-7802.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu W, Laitinen S, Khan S, Vihinen M, Kowalski J, Yu G, Chen L, Ewing CM, Eisenberger MA, Carducci MA, Nelson WG, Yegnasubramanian S, Luo J, Wang Y, Xu J, Isaacs WB, Visakorpi T, Bova GS: Copy number analysis indicates monoclonal origin of lethal metastatic prostate cancer. Nat Med. 2009, 15: 559-565.View ArticlePubMedPubMed CentralGoogle Scholar
- Taylor BS, Schultz N, Hieronymus H, Gopalan A, Xiao Y, Carver BS, Arora VK, Kaushik P, Cerami E, Reva B, Antipin Y, Mitsiades N, Landers T, Dolgalev I, Major JE, Wilson M, Socci ND, Lash AE, Heguy A, Eastham JA, Scher HI, Reuter VE, Scardino PT, Sander C, Sawyers CL, Gerald WL: Integrative genomic profiling of human prostate cancer. Cancer Cell. 2010, 18: 11-22.View ArticlePubMedPubMed CentralGoogle Scholar
- Robbins CM, Tembe WA, Baker A, Sinari S, Moses TY, Beckstrom-Sternberg S, Beckstrom-Sternberg J, Barrett M, Long J, Chinnaiyan A, Lowey J, Suh E, Pearson JV, Craig DW, Agus DB, Pienta KJ, Carpten JD: Copy number and targeted mutational analysis reveals novel somatic events in metastatic prostate tumors. Genome Res. 2011, 21: 47-55.View ArticlePubMedPubMed CentralGoogle Scholar
- Barbieri CE, Baca SC, Lawrence MS, Demichelis F, Blattner M, Theurillat JP, White TA, Stojanov P, Van Allen E, Stransky N, Nickerson E, Chae SS, Boysen G, Auclair D, Onofrio RC, Park K, Kitabayashi N, MacDonald TY, Sheikh K, Vuong T, Guiducci C, Cibulskis K, Sivachenko A, Carter SL, Saksena G, Voet D, Hussain WM, Ramos AH, Winckler W, Redman MC, et al.: Exome sequencing identifies recurrent SPOP, FOXA1 and MED12 mutations in prostate cancer. Nat Genet. 2012, 44: 685-689.View ArticlePubMedPubMed CentralGoogle Scholar
- Berger MF, Lawrence MS, Demichelis F, Drier Y, Cibulskis K, Sivachenko AY, Sboner A, Esgueva R, Pflueger D, Sougnez C, Onofrio R, Carter SL, Park K, Habegger L, Ambrogio L, Fennell T, Parkin M, Saksena G, Voet D, Ramos AH, Pugh TJ, Wilkinson J, Fisher S, Winckler W, Mahan S, Ardlie K, Baldwin J, Simons JW, Kitabayashi N, MacDonald TY, et al.: The genomic complexity of primary human prostate cancer. Nature. 2011, 470: 214-220.View ArticlePubMedPubMed CentralGoogle Scholar
- Grasso CS, Wu YM, Robinson DR, Cao X, Dhanasekaran SM, Khan AP, Quist MJ, Jing X, Lonigro RJ, Brenner JC, Asangani IA, Ateeq B, Chun SY, Siddiqui J, Sam L, Anstett M, Mehra R, Prensner JR, Palanisamy N, Ryslik GA, Vandin F, Raphael BJ, Kunju LP, Rhodes DR, Pienta KJ, Chinnaiyan AM, Tomlins SA: The mutational landscape of lethal castration-resistant prostate cancer. Nature. 2012, 487: 239-243.View ArticlePubMedPubMed CentralGoogle Scholar
- Kumar A, White TA, MacKenzie AP, Clegg N, Lee C, Dumpit RF, Coleman I, Ng SB, Salipante SJ, Rieder MJ, Nickerson DA, Corey E, Lange PH, Morrissey C, Vessella RL, Nelson PS, Shendure J: Exome sequencing identifies a spectrum of mutation frequencies in advanced and lethal prostate cancers. Proc Natl Acad Sci USA. 2011, 108: 17087-17092.View ArticlePubMedPubMed CentralGoogle Scholar
- Rubin MA, Putzi M, Mucci N, Smith DC, Wojno K, Korenchuk S, Pienta KJ: Rapid ("warm") autopsy study for procurement of metastatic prostate cancer. Clin Cancer Res. 2000, 6: 1038-1045.PubMedGoogle Scholar
- Fiorentino M, Capizzi E, Loda M: Blood and tissue biomarkers in prostate cancer: state of the art. Urol Clin North Am. 2010, 37: 131-141. Table of Contents,View ArticlePubMedPubMed CentralGoogle Scholar
- Attard G, Swennenhuis JF, Olmos D, Reid AH, Vickers E, A'Hern R, Levink R, Coumans F, Moreira J, Riisnaes R, Oommen NB, Hawche G, Jameson C, Thompson E, Sipkema R, Carden CP, Parker C, Dearnaley D, Kaye SB, Cooper CS, Molina A, Cox ME, Terstappen LW, de Bono JS: Characterization of ERG, AR and PTEN gene status in circulating tumor cells from patients with castration-resistant prostate cancer. Cancer Res. 2009, 69: 2912-2918.View ArticlePubMedGoogle Scholar
- Danila DC, Fleisher M, Scher HI: Circulating tumor cells as biomarkers in prostate cancer. Clin Cancer Res. 2011, 17: 3903-3912.View ArticlePubMedPubMed CentralGoogle Scholar
- de Bono JS, Scher HI, Montgomery RB, Parker C, Miller MC, Tissing H, Doyle GV, Terstappen LW, Pienta KJ, Raghavan D: Circulating tumor cells predict survival benefit from treatment in metastatic castration-resistant prostate cancer. Clin Cancer Res. 2008, 14: 6302-6309.View ArticlePubMedGoogle Scholar
- Olmos D, Arkenau HT, Ang JE, Ledaki I, Attard G, Carden CP, Reid AH, A'Hern R, Fong PC, Oomen NB, Molife R, Dearnaley D, Parker C, Terstappen LW, de Bono JS: Circulating tumour cell (CTC) counts as intermediate end points in castration-resistant prostate cancer (CRPC): a single-centre experience. Ann Oncol. 2009, 20: 27-33.View ArticlePubMedGoogle Scholar
- Gleghorn JP, Pratt ED, Denning D, Liu H, Bander NH, Tagawa ST, Nanus DM, Giannakakou PA, Kirby BJ: Capture of circulating tumor cells from whole blood of prostate cancer patients using geometrically enhanced differential immunocapture (GEDI) and a prostate-specific antibody. Lab Chip. 2010, 10: 27-29.View ArticlePubMedPubMed CentralGoogle Scholar
- Stott SL, Hsu CH, Tsukrov DI, Yu M, Miyamoto DT, Waltman BA, Rothenberg SM, Shah AM, Smas ME, Korir GK, Floyd FP, Gilman AJ, Lord JB, Winokur D, Springer S, Irimia D, Nagrath S, Sequist LV, Lee RJ, Isselbacher KJ, Maheswaran S, Haber DA, Toner M: Isolation of circulating tumor cells using a microvortex-generating herringbone-chip. Proc Natl Acad Sci USA. 2010, 107: 18392-18397.View ArticlePubMedPubMed CentralGoogle Scholar
- Stott SL, Lee RJ, Nagrath S, Yu M, Miyamoto DT, Ulkus L, Inserra EJ, Ulman M, Springer S, Nakamura Z, Moore AL, Tsukrov DI, Kempner ME, Dahl DM, Wu CL, Iafrate AJ, Smith MR, Tompkins RG, Sequist LV, Toner M, Haber DA, Maheswaran S: Isolation and characterization of circulating tumor cells from patients with localized and metastatic prostate cancer. Sci Transl Med. 2010, 2: 25ra23-View ArticlePubMedPubMed CentralGoogle Scholar
- Miyamoto DT, Lee RJ, Stott SL, Ting DT, Wittner BS, Ulman M, Smas ME, Lord JB, Brannigan BW, Trautwein J, Bander NH, Wu CL, Sequist LV, Smith MR, Ramaswamy S, Toner M, Maheswaran S, Haber DA: Androgen receptor signaling in circulating tumor cells as a marker of hormonally responsive prostate cancer. Cancer Discov. 2012, 2: 995-1003.View ArticlePubMedPubMed CentralGoogle Scholar
- Olmos D, Brewer D, Clark J, Danila DC, Parker C, Attard G, Fleisher M, Reid AH, Castro E, Sandhu SK, Barwell L, Oommen NB, Carreira S, Drake CG, Jones R, Cooper CS, Scher HI, de Bono JS: Prognostic value of blood mRNA expression signatures in castration-resistant prostate cancer: a prospective, two-stage study. Lancet Oncol. 2012, 13: 1114-1124.View ArticlePubMedGoogle Scholar
- Ross RW, Galsky MD, Scher HI, Magidson J, Wassmann K, Lee GS, Katz L, Subudhi SK, Anand A, Fleisher M, Kantoff PW, Oh WK: A whole-blood RNA transcript-based prognostic model in men with castration-resistant prostate cancer: a prospective study. Lancet Oncol. 2012, 13: 1105-1113.View ArticlePubMedGoogle Scholar
- Schwarzenbach H, Hoon DS, Pantel K: Cell-free nucleic acids as biomarkers in cancer patients. Nat Rev Cancer. 2011, 11: 426-437.View ArticlePubMedGoogle Scholar
- Leary RJ, Kinde I, Diehl F, Schmidt K, Clouser C, Duncan C, Antipova A, Lee C, McKernan K, De La Vega FM, Kinzler KW, Vogelstein B, Diaz LA, Velculescu VE: Development of personalized tumor biomarkers using massively parallel sequencing. Sci Trans Med. 2010, 2: 20ra14-View ArticleGoogle Scholar
- McBride DJ, Orpana AK, Sotiriou C, Joensuu H, Stephens PJ, Mudie LJ, Hamalainen E, Stebbings LA, Andersson LC, Flanagan AM, Durbecq V, Ignatiadis M, Kallioniemi O, Heckman CA, Alitalo K, Edgren H, Futreal PA, Stratton MR, Campbell PJ: Use of cancer-specific genomic rearrangements to quantify disease burden in plasma from patients with solid tumors. Genes Chromosomes Cancer. 2010, 49: 1062-1069.View ArticlePubMedPubMed CentralGoogle Scholar
- Diehl F, Li M, Dressman D, He Y, Shen D, Szabo S, Diaz LA, Goodman SN, David KA, Juhl H, Kinzler KW, Vogelstein B: Detection and quantification of mutations in the plasma of patients with colorectal tumors. Proc Natl Acad Sci USA. 2005, 102: 16368-16373.View ArticlePubMedPubMed CentralGoogle Scholar
- Diehl F, Schmidt K, Choti MA, Romans K, Goodman S, Li M, Thornton K, Agrawal N, Sokoll L, Szabo SA, Kinzler KW, Vogelstein B, Diaz LA: Circulating mutant DNA to assess tumor dynamics. Nat Med. 2008, 14: 985-990.View ArticlePubMedPubMed CentralGoogle Scholar
- Forshew T, Murtaza M, Parkinson C, Gale D, Tsui DW, Kaper F, Dawson SJ, Piskorz AM, Jimenez-Linan M, Bentley D, Hadfield J, May AP, Caldas C, Brenton JD, Rosenfeld N: Noninvasive identification and monitoring of cancer mutations by targeted deep sequencing of plasma DNA. Sci Trans Med. 2012, 4: 136ra168-View ArticleGoogle Scholar
- Heitzer E, Auer M, Hoffmann EM, Pichler M, Gasch C, Ulz P, Lax S, Waldispuehl-Geigl J, Mauermann O, Mohan S, Pristauz G, Lackner C, Hofler G, Eisner F, Petru E, Sill H, Samonigg H, Pantel K, Riethdorf S, Bauernhofer T, Geigl JB, Speicher MR: Establishment of tumor-specific copy number alterations from plasma DNA of patients with cancer. Int J Cancer. 2013, 15: 346-356.View ArticleGoogle Scholar
- Chiu RW, Chan KC, Gao Y, Lau VY, Zheng W, Leung TY, Foo CH, Xie B, Tsui NB, Lun FM, Zee BC, Lau TK, Cantor CR, Lo YM: Noninvasive prenatal diagnosis of fetal chromosomal aneuploidy by massively parallel genomic sequencing of DNA in maternal plasma. Proc Natl Acad Sci USA. 2008, 105: 20458-20463.View ArticlePubMedPubMed CentralGoogle Scholar
- Chiu RW, Lo YM: Clinical applications of maternal plasma fetal DNA analysis: translating the fruits of 15 years of research. Clin Chem Lab Med. 2012, 51: 197-204.Google Scholar
- Fan HC, Blumenfeld YJ, Chitkara U, Hudgins L, Quake SR: Noninvasive diagnosis of fetal aneuploidy by shotgun sequencing DNA from maternal blood. Proc Natl Acad Sci USA. 2008, 105: 16266-16271.View ArticlePubMedPubMed CentralGoogle Scholar
- Chan KC, Jiang P, Zheng YW, Liao GJ, Sun H, Wong J, Siu SS, Chan WC, Chan SL, Chan AT, Lai PB, Chiu RW, Lo YM: Cancer genome scanning in plasma: detection of tumor-associated copy number aberrations, single-nucleotide variants, and tumoral heterogeneity by massively parallel sequencing. Clin Chem. 2013, 59: 211-224.View ArticlePubMedGoogle Scholar
- Leary RJ, Sausen M, Kinde I, Papadopoulos N, Carpten JD, Craig D, O'Shaughnessy J, Kinzler KW, Parmigiani G, Vogelstein B, Diaz LA, Velculescu VE: Detection of chromosomal alterations in the circulation of cancer patients with whole-genome sequencing. Sci Trans Med. 2012, 4: 162ra154-View ArticleGoogle Scholar
- Navin N, Kendall J, Troge J, Andrews P, Rodgers L, McIndoo J, Cook K, Stepansky A, Levy D, Esposito D, Muthuswamy L, Krasnitz A, McCombie WR, Hicks J, Wigler M: Tumour evolution inferred by single-cell sequencing. Nature. 2011, 472: 90-94.View ArticlePubMedPubMed CentralGoogle Scholar
- Loman NJ, Misra RV, Dallman TJ, Constantinidou C, Gharbia SE, Wain J, Pallen MJ: Performance comparison of benchtop high-throughput sequencing platforms. Nat Biotechnol. 2012, 30: 434-439.View ArticlePubMedGoogle Scholar
- Riethdorf S, Fritsche H, Muller V, Rau T, Schindlbeck C, Rack B, Janni W, Coith C, Beck K, Janicke F, Jackson S, Gornet T, Cristofanilli M, Pantel K: Detection of circulating tumor cells in peripheral blood of patients with metastatic breast cancer: a validation study of the CellSearch system. Clin Cancer Res. 2007, 13: 920-928.View ArticlePubMedGoogle Scholar
- Heitzer E, Auer M, Gasch C, Pichler M, Ulz P, Hoffmann EM, Lax S, Waldispuehl-Geigl J, Mauermann O, Lackner C, Hofler G, Eisner F, Sill H, Samonigg H, Pantel K, Riethdorf S, Bauernhofer T, Geigl JB, Speicher MR: Complex tumor genomes inferred from single circulating tumor cells by array-CGH and next-generation sequencing. Cancer Res. 2013, 73: 2965-2975.View ArticlePubMedGoogle Scholar
- Geigl JB, Obenauf AC, Waldispuehl-Geigl J, Hoffmann EM, Auer M, Hormann M, Fischer M, Trajanoski Z, Schenk MA, Baumbusch LO, Speicher MR: Identification of small gains and losses in single cells after whole genome amplification on tiling oligo arrays. Nucleic Acids Res. 2009, 37: e105-View ArticlePubMedPubMed CentralGoogle Scholar
- Geigl JB, Speicher MR: Single-cell isolation from cell suspensions and whole genome amplification from single cells to provide templates for CGH analysis. Nat Protoc. 2007, 2: 3173-3184.View ArticlePubMedGoogle Scholar
- Baslan T, Kendall J, Rodgers L, Cox H, Riggs M, Stepansky A, Troge J, Ravi K, Esposito D, Lakshmi B, Wigler M, Navin N, Hicks J: Genome-wide copy number analysis of single cells. Nat Protoc. 2012, 7: 1024-1041.View ArticlePubMedGoogle Scholar
- Olshen AB, Venkatraman ES, Lucito R, Wigler M: Circular binary segmentation for the analysis of array-based DNA copy number data. Biostatistics. 2004, 5: 557-572.View ArticlePubMedGoogle Scholar
- Hupe P, Stransky N, Thiery JP, Radvanyi F, Barillot E: Analysis of array CGH data: from signal ratio to gain and loss of DNA regions. Bioinformatics. 2004, 20: 3413-3422.View ArticlePubMedGoogle Scholar
- Lai W, Choudhary V, Park PJ: CGHweb: a tool for comparing DNA copy number segmentations from multiple algorithms. Bioinformatics. 2008, 24: 1014-1015.View ArticlePubMedPubMed CentralGoogle Scholar
- R: A language and environment for statistical computing.http://www.R-project.org
- Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, Muller M: pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011, 12: 77-View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760.View ArticlePubMedPubMed CentralGoogle Scholar
- McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA: The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20: 1297-1303.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang K, Li M, Hakonarson H: ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38: e164-View ArticlePubMedPubMed CentralGoogle Scholar
- Genomes Project C: A map of human genome variation from population-scale sequencing. Nature. 2010, 467: 1061-1073.View ArticleGoogle Scholar
- Le Scouarnec S, Gribble SM: Characterising chromosome rearrangements: recent technical advances in molecular cytogenetics. Heredity. 2012, 108: 75-85.View ArticlePubMedPubMed CentralGoogle Scholar
- Kent WJ: BLAT--the BLAST-like alignment tool. Genome Res. 2002, 12: 656-664.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen EZ, Chiu RW, Sun H, Akolekar R, Chan KC, Leung TY, Jiang P, Zheng YW, Lun FM, Chan LY, Jin Y, Go AT, Lau ET, To WW, Leung WC, Tang RY, Au-Yeung SK, Lam H, Kung YY, Zhang X, van Vugt JM, Minekawa R, Tang MH, Wang J, Oudejans CB, Lau TK, Nicolaides KH, Lo YM: Noninvasive prenatal diagnosis of fetal trisomy 18 and trisomy 13 by maternal plasma DNA sequencing. PloS One. 2011, 6: e21791-View ArticlePubMedPubMed CentralGoogle Scholar
- Fan HC, Quake SR: Sensitivity of noninvasive prenatal detection of fetal aneuploidy from maternal plasma using shotgun sequencing is limited only by counting statistics. PloS One. 2010, 5: e10439-View ArticlePubMedPubMed CentralGoogle Scholar
- Palomaki GE, Kloza EM, Lambert-Messerlian GM, Haddow JE, Neveux LM, Ehrich M, van den Boom D, Bombard AT, Deciu C, Grody WW, Nelson SF, Canick JA: DNA sequencing of maternal plasma to detect Down syndrome: an international clinical validation study. Genet Med. 2011, 13: 913-920.View ArticlePubMedGoogle Scholar
- Corzo C, Petzold M, Mayol X, Espinet B, Salido M, Serrano S, Real FX, Sole F: RxFISH karyotype and MYC amplification in the HT-29 colon adenocarcinoma cell line. Genes Chromosomes Cancer. 2003, 36: 425-426.View ArticlePubMedGoogle Scholar
- Kawai K, Viars C, Arden K, Tarin D, Urquidi V, Goodison S: Comprehensive karyotyping of the HT-29 colon adenocarcinoma cell line. Genes Chromosomes Cancer. 2002, 34: 1-8.View ArticlePubMedGoogle Scholar
- Tomlins SA, Laxman B, Dhanasekaran SM, Helgeson BE, Cao X, Morris DS, Menon A, Jing X, Cao Q, Han B, Yu J, Wang L, Montie JE, Rubin MA, Pienta KJ, Roulston D, Shah RB, Varambally S, Mehra R, Chinnaiyan AM: Distinct classes of chromosomal rearrangements create oncogenic ETS gene fusions in prostate cancer. Nature. 2007, 448: 595-599.View ArticlePubMedGoogle Scholar
- Tomlins SA, Rhodes DR, Perner S, Dhanasekaran SM, Mehra R, Sun XW, Varambally S, Cao X, Tchinda J, Kuefer R, Lee C, Montie JE, Shah RB, Pienta KJ, Rubin MA, Chinnaiyan AM: Recurrent fusion of TMPRSS2 and ETS transcription factor genes in prostate cancer. Science. 2005, 310: 644-648.View ArticlePubMedGoogle Scholar
- Tomlins SA, Bjartell A, Chinnaiyan AM, Jenster G, Nam RK, Rubin MA, Schalken JA: ETS gene fusions in prostate cancer: from discovery to daily clinical practice. Eur Urol. 2009, 56: 275-286.View ArticlePubMedGoogle Scholar
- Visakorpi T, Hyytinen E, Koivisto P, Tanner M, Keinanen R, Palmberg C, Palotie A, Tammela T, Isola J, Kallioniemi OP: In vivo amplification of the androgen receptor gene and progression of human prostate cancer. Nat Genet. 1995, 9: 401-406.View ArticlePubMedGoogle Scholar
- Lim J, Ghadessy FJ, Abdullah AA, Pinsky L, Trifiro M, Yong EL: Human androgen receptor mutation disrupts ternary interactions between ligand, receptor domains, and the coactivator TIF2 (transcription intermediary factor 2). Mol Endocrinol. 2000, 14: 1187-1197.View ArticlePubMedGoogle Scholar
- Pourmand G, Ziaee AA, Abedi AR, Mehrsai A, Alavi HA, Ahmadi A, Saadati HR: Role of PTEN gene in progression of prostate cancer. Urol J. 2007, 4: 95-100.PubMedGoogle Scholar
- Chen M, Pratt CP, Zeeman ME, Schultz N, Taylor BS, O'Neill A, Castillo-Martin M, Nowak DG, Naguib A, Grace DM, Murn J, Navin N, Atwal GS, Sander C, Gerald WL, Cordon-Cardo C, Newton AC, Carver BS, Trotman LC: Identification of PHLPP1 as a tumor suppressor reveals the role of feedback activation in PTEN-mutant prostate cancer progression. Cancer Cell. 2011, 20: 173-186.View ArticlePubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.