A pipeline for copy number profiling of single circulating tumour cells to assess intrapatient tumour heterogeneity.

Intrapatient tumour heterogeneity is likely a major determinant of clinical outcome in cancer patients. To assess heterogeneity in a minimally invasive manner, methods to perform single circulating tumour cell (CTC) genomics at high resolution are necessary. However, due to the rarity of CTCs, development of such methods is challenging. Here, we developed a modular single CTC analysis pipeline to assess intrapatient heterogeneity by copy number (CN) profiling. To optimize this pipeline, spike-in experiments using MCF-7 breast cancer cells were performed. The VyCAP puncher system was used to isolate single cells. The quality of whole genome amplification (WGA) products generated by REPLI-g and Ampli1™ methods, as well as the results from the Illumina Truseq and the Ampli1™ LowPass library preparation techniques, was compared. Moreover, a bioinformatic pipeline was designed to generate CN profiles from single CTCs. The optimal combination of Ampli1™ WGA and Illumina Truseq library preparation was successfully validated on patient-derived CTCs. In conclusion, we developed a novel modular pipeline to isolate single CTCs and subsequently generate detailed patient-derived CN profiles that allow assessment of intrapatient heterogeneity in future studies.

11-Ketotestosterone is the predominant active androgen in prostate cancer patients after castration.

BACKGROUNDContinued androgen receptor (AR) signaling constitutes a key target for treatment in metastatic castration-resistant prostate cancer (CRPC). Studies have identified 11-ketotestosterone (11KT) as a potent AR agonist, but it is unknown if 11KT is present at physiologically relevant concentrations in patients with CRPC to drive AR activation. The goal of this study was to investigate the circulating steroid metabolome including all active androgens in patients with CRPC.METHODSPatients with metastatic CRPC (n = 29) starting a new line of systemic therapy were included. Sequential plasma samples were obtained for measurement of circulating steroid concentrations by multisteroid profiling employing liquid chromatography-tandem mass spectrometry. Metastatic tumor biopsy samples were obtained at baseline and subjected to RNA sequencing.RESULTS11KT was the most abundant circulating active androgen in 97% of patients with CRPC (median 0.39 nmol/L, range: 0.03-2.39 nmol/L), constituting 60% (IQR 43%-79%) of the total active androgen (TA) pool. Treatment with glucocorticoids reduced 11KT by 84% (49%-89%) and testosterone by 68% (38%-79%). Circulating TA concentrations at baseline were associated with a distinct intratumor gene expression signature comprising AR-regulated genes.CONCLUSIONThe potent AR agonist 11KT is the predominant circulating active androgen in patients with CRPC and, therefore, one of the potential drivers of AR activation in CRPC. Assessment of androgen status should be extended to include 11KT, as current clinical approaches likely underestimate androgen abundance in patients with CRPC.TRIAL REGISTRATIONNetherlands Trial Register: NL5625 (NTR5732).FUNDINGDaniel den Hoed Foundation and Wellcome Trust (Investigator Award WT209492/Z/17/Z).

Detection of tumor-derived extracellular vesicles in plasma from patients with solid cancer.

BACKGROUND: Extracellular vesicles (EVs) are actively secreted by cells into body fluids and contain nucleic acids of the cells they originate from. The goal of this study was to detect circulating tumor-derived EVs (ctEVs) by mutant mRNA transcripts (EV-RNA) in plasma of patients with solid cancers and compare the occurrence of ctEVs with circulating tumor DNA (ctDNA) in cell-free DNA (cfDNA). METHODS: For this purpose, blood from 20 patients and 15 healthy blood donors (HBDs) was collected in different preservation tubes (EDTA, BCT, CellSave) and processed into plasma within 24 h from venipuncture. EVs were isolated with the ExoEasy protocol from this plasma and from conditioned medium of 6 cancer cell lines and characterized according to MISEV2018-guidelines. RNA from EVs was isolated with the ExoRNeasy protocol and evaluated for transcript expression levels of 96 genes by RT-qPCR and genotyped by digital PCR. RESULTS: Our workflow applied on cell lines revealed a high concordance between cellular mRNA and EV-RNA in expression levels as well as variant allele frequencies for PIK3CA, KRAS and BRAF. Plasma CD9-positive EV and GAPDH EV-RNA levels were significantly different between the preservation tubes. The workflow detected only ctEVs with mutant transcripts in plasma of patients with high amounts (> 20%) of circulating tumor DNA (ctDNA). Expression profiling showed that the EVs from patients resemble healthy donors more than tumor cell lines supporting that most EVs are derived from healthy tissue. CONCLUSIONS: We provide a workflow for ctEV detection by spin column-based generic isolation of EVs and PCR-based measurement of gene expression and mutant transcripts in EV-RNA derived from cancer patients' blood plasma. This workflow, however, detected tumor-specific mutations in blood less often in EV-RNA than in cfDNA.

Circulating Tumor Cell Enumeration and Characterization in Metastatic Castration-Resistant Prostate Cancer Patients Treated with Cabazitaxel.

(1) Background: Markers identifying which patients with metastatic, castration-resistant prostate cancer (mCRPC) will benefit from cabazitaxel therapy are currently lacking. Therefore, the aim of this study was to identify markers associated with outcome to cabazitaxel therapy based on counts and gene expression profiles of circulating tumor cells (CTCs). (2) Methods: From 120 mCRPC patients, CellSearch enriched CTCs were obtained at baseline and after 6 weeks of cabazitaxel therapy. Furthermore, 91 genes associated with prostate cancer were measured in mRNA of these CTCs. (3) Results: In 114 mCRPC patients with an evaluable CTC count, the CTC count was independently associated with poor progression-free survival (PFS) and overall survival (OS) in multivariable analysis with other commonly used variables associated with outcome in mCRPC (age, prostate specific antigen (PSA), alkaline phosphatase, lactate dehydrogenase (LDH), albumin, hemoglobin), together with alkaline phosphatase and hemoglobin. A five-gene expression profile was generated to predict for outcome to cabazitaxel therapy. However, even though this signature was associated with OS in univariate analysis, this was not the case in the multivariate analysis for OS nor for PFS. (4) Conclusion: The established five-gene expression profile in CTCs was not independently associated with PFS nor OS. However, along with alkaline phosphatase and hemoglobin, CTC-count is independently associated with PFS and OS in mCRPC patients who are treated with cabazitaxel.

The circular RNome of primary breast cancer.

Circular RNAs (circRNAs) are a class of RNAs that is under increasing scrutiny, although their functional roles are debated. We analyzed RNA-seq data of 348 primary breast cancers and developed a method to identify circRNAs that does not rely on unmapped reads or known splice junctions. We identified 95,843 circRNAs, of which 20,441 were found recurrently. Of the circRNAs that match exon boundaries of the same gene, 668 showed a poor or even negative (R < 0.2) correlation with the expression level of the linear gene. In silico analysis showed only a minority (8.5%) of circRNAs could be explained by known splicing events. Both these observations suggest that specific regulatory processes for circRNAs exist. We confirmed the presence of circRNAs of CNOT2, CREBBP, and RERE in an independent pool of primary breast cancers. We identified circRNA profiles associated with subgroups of breast cancers and with biological and clinical features, such as amount of tumor lymphocytic infiltrate and proliferation index. siRNA-mediated knockdown of circCNOT2 was shown to significantly reduce viability of the breast cancer cell lines MCF-7 and BT-474, further underlining the biological relevance of circRNAs. Furthermore, we found that circular, and not linear, CNOT2 levels are predictive for progression-free survival time to aromatase inhibitor (AI) therapy in advanced breast cancer patients, and found that circCNOT2 is detectable in cell-free RNA from plasma. We showed that circRNAs are abundantly present, show characteristics of being specifically regulated, are associated with clinical and biological properties, and thus are relevant in breast cancer.

High-throughput isolation of circulating tumor DNA: a comparison of automated platforms.

The emerging interest in circulating tumor DNA (ctDNA) analyses for clinical trials has necessitated the development of a high-throughput method for fast, reproducible, and efficient isolation of ctDNA. Currently, the majority of ctDNA studies use the manual QIAamp (QA) platform to isolate DNA from blood. The purpose of this study was to compare two competing automated DNA isolation platforms [Maxwell (MX) and QIAsymphony (QS)] to the current 'gold standard' QA to facilitate high-throughput processing of samples in prospective trials. We obtained blood samples from healthy blood donors and metastatic cancer patients for plasma isolation. Total cell-free DNA (cfDNA) quantity was assessed by TERT quantitative PCR. Recovery efficiency was investigated by quantitative PCR analysis of spiked-in synthetic plant DNA. In addition, a ß-actin fragmentation assay was performed to determine the amount of contamination by genomic DNA from lysed leukocytes. ctDNA quality was assessed by digital PCR for somatic variant detection. cfDNA quantity and recovery efficiency were lowest using the MX platform, whereas QA and QS showed a comparable performance. All platforms preferentially isolated small (136 bp) DNA fragments over large (420 and 2000 bp) DNA fragments. Detection of the number variant and wild-type molecules was most comparable between QA and QS. However, there was no significant difference in variant allele frequency comparing QS and MX to QA. In summary, we show that the QS platform has comparable performance to QA, the 'gold standard', and outperformed the MX platform depending on the readout used. We conclude that the QS can replace the more laborious QA platform, especially when high-throughput cfDNA isolation is needed.

Gene length corrected trimmed mean of M-values (GeTMM) processing of RNA-seq data performs similarly in intersample analyses while improving intrasample comparisons.

BACKGROUND: Current normalization methods for RNA-sequencing data allow either for intersample comparison to identify differentially expressed (DE) genes or for intrasample comparison for the discovery and validation of gene signatures. Most studies on optimization of normalization methods typically use simulated data to validate methodologies. We describe a new method, GeTMM, which allows for both inter- and intrasample analyses with the same normalized data set. We used actual (i.e. not simulated) RNA-seq data from 263 colon cancers (no biological replicates) and used the same read count data to compare GeTMM with the most commonly used normalization methods (i.e. TMM (used by edgeR), RLE (used by DESeq2) and TPM) with respect to distributions, effect of RNA quality, subtype-classification, recurrence score, recall of DE genes and correlation to RT-qPCR data. RESULTS: We observed a clear benefit for GeTMM and TPM with regard to intrasample comparison while GeTMM performed similar to TMM and RLE normalized data in intersample comparisons. Regarding DE genes, recall was found comparable among the normalization methods, while GeTMM showed the lowest number of false-positive DE genes. Remarkably, we observed limited detrimental effects in samples with low RNA quality. CONCLUSIONS: We show that GeTMM outperforms established methods with regard to intrasample comparison while performing equivalent with regard to intersample normalization using the same normalized data. These combined properties enhance the general usefulness of RNA-seq but also the comparability to the many array-based gene expression data in the public domain.

Confirmation of a metastasis-specific microRNA signature in primary colon cancer.

The identification of patients with high-risk stage II colon cancer who may benefit from adjuvant therapy may allow the clinical approach to be tailored for these patients based on an understanding of tumour biology. MicroRNAs have been proposed as markers of the prognosis or treatment response in colorectal cancer. Recently, a 2-microRNA signature (let-7i and miR-10b) was proposed to identify colorectal cancer patients at risk of developing distant metastasis. We assessed the prognostic value of this signature and additional candidate microRNAs in an independent, clinically well-defined, prospectively collected cohort of primary colon cancer patients including stage I-II colon cancer without and stage III colon cancer with adjuvant treatment. The 2-microRNA signature specifically predicted hepatic recurrence in the stage I-II group, but not the overall ability to develop distant metastasis. The addition of miR-30b to the 2-microRNA signature allowed the prediction of both distant metastasis and hepatic recurrence in patients with stage I-II colon cancer who did not receive adjuvant chemotherapy. Available gene expression data allowed us to associate miR-30b expression with axon guidance and let-7i expression with cell adhesion, migration, and motility.

An In-Depth Evaluation of the Validity and Logistics Surrounding the Testing of AR-V7 mRNA Expression in Circulating Tumor Cells.

Recent reports have emphasized the clinical relevance of detecting AR-V7 in circulating tumor cells (CTCs). Our aim was to set up a validated multicenter pipeline to measure AR-V7 by quantitative RT-PCR (RT-qPCR) in RNA isolated from CellSearch-enriched CTCs to provide an AR-V7-positive or AR-V7-negative score in a clinically acceptable time range. CellSearch-enirched CTCs from patients with metastatic castration-resistant prostate cancer were characterized by RT-qPCR. After optimization, it was prospectively tested whether it was possible to report the AR-V7 status within 11 days (PRELUDE study). In the range of the RNA equivalent of 0.2 to 12 VCaP cells, the CV for AR-V7 was 9% (n = 37). The limit of detection was 0.3, and the limit of quantitation was 3 cells in the final RT-qPCR. No differences were observed between AR-V7 data generated by five technicians or in two different laboratories. For the 45 patients in PRELUDE, 13 patients were ineligible, 22 patients were AR-V7 negative, and 10 were AR-V7 positive. The median time to inform the physician of the test result was 7 days (range, 2 to 11 days). This assay can establish the AR-V7 status in CTCs from patients with metastatic castration-resistant prostate cancer. Furthermore, it was possible to provide an AR-V7 outcome within 11 days, indicating that it may be used to choose between an anti-androgen receptor or taxane-based cabazitaxel treatment.

Estrogen receptor mutations and splice variants determined in liquid biopsies from metastatic breast cancer patients.

Mutations and splice variants in the estrogen receptor (ER) gene, ESR1, may yield endocrine resistance in metastatic breast cancer (MBC) patients. These putative endocrine resistance markers are likely to emerge during treatment, and therefore, its detection in liquid biopsies, such as circulating tumor cells (CTCs) and cell-free DNA (cfDNA), is of great interest. This research aimed to determine whether ESR1 mutations and splice variants occur more frequently in CTCs of MBC patients progressing on endocrine treatment. In addition, the presence of ESR1 mutations was evaluated in matched cfDNA and compared to CTCs. CellSearch-enriched CTC fractions (≥5/7.5 mL) of two MBC cohorts were evaluated, namely (a) patients starting first-line endocrine therapy (n = 43, baseline cohort) and (b) patients progressing on any line of endocrine therapy (n = 40, progressing cohort). ESR1 hotspot mutations (D538G and Y537S/N/C) were evaluated in CTC-enriched DNA using digital PCR and compared with matched cfDNA (n = 18 baseline cohort; n = 26 progressing cohort). Expression of ESR1 full-length and 4 of its splice variants (∆5, ∆7, 36 kDa, and 46 kDa) was evaluated in CTC-enriched mRNA. It was observed that in the CTCs, the ESR1 mutations were not enriched in the progressing cohort (8%), when compared with the baseline cohort (5%) (P = 0.66). In the cfDNA, however, ESR1 mutations were more prevalent in the progressing cohort (42%) than in the baseline cohort (11%) (P = 0.04). Three of the same mutations were observed in both CTCs and cfDNA, 1 mutation in CTCs only, and 11 in cfDNA only. Only the ∆5 ESR1 splice variant was CTC-specific expressed, but was not enriched in the progressing cohort. In conclusion, sensitivity for detecting ESR1 mutations in CTC-enriched fractions was lower than for cfDNA. ESR1 mutations detected in cfDNA, rarely present at the start of first-line endocrine therapy, were enriched at progression, strongly suggesting a role in conferring endocrine resistance in MBC.

The 29.5 kb APOBEC3B Deletion Polymorphism Is Not Associated with Clinical Outcome of Breast Cancer.

Increased APOBEC3B mRNA levels are associated with a hypermutator phenotype and poor prognosis in ER-positive breast cancer patients. In addition, a 29.5 kb deletion polymorphism of APOBEC3B, resulting in an APOBEC3A-B hybrid transcript, has been associated with an increased breast cancer risk and the hypermutator phenotype. Here we evaluated whether the APOBEC3B deletion polymorphism also associates with clinical outcome of breast cancer. Copy number analysis was performed by quantitative PCR (qPCR) in primary tumors of 1,756 Dutch breast cancer patients. The APOBEC3B deletion was found in 187 patients of whom 16 carried a two-copy deletion and 171 carried a one-copy deletion. The prognostic value of the APOBEC3B deletion for the natural course of the disease was evaluated among 1,076 lymph-node negative (LNN) patients who did not receive adjuvant systemic treatment. No association was found between APOBEC3B copy number values and the length of metastasis-free survival (MFS; hazard ratio (HR) = 1.00, 95% confidence interval (CI) = 0.90-1.11, P = 0.96). Subgroup analysis by ER status also did not reveal an association between APOBEC3B copy number values and the length of MFS. The predictive value of the APOBEC3B deletion was assessed among 329 ER-positive breast cancer patients who received tamoxifen as the first-line therapy for recurrent disease and 226 breast cancer patients who received first-line chemotherapy for recurrent disease. No association between APOBEC3B copy number values and the overall response rate (ORR) to either tamoxifen (odds ratio (OR) = 0.88, 95% CI = 0.69-1.13, P = 0.31) or chemotherapy (OR = 0.97, 95% CI = 0.71-1.33, P = 0.87) was found. Thus, in contrast to APOBEC3B mRNA levels, the APOBEC3B deletion polymorphism has neither a prognostic nor a predictive value for breast cancer patients. Although a correlation exists between APOBEC3B copy number and mRNA expression, it is relatively weak. This suggests that other mechanisms exist that may affect and therefore determine the prognostic value of APOBEC3B mRNA levels.