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

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 ((increment)5, (increment)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 (increment)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

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 (l et-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 m iR-30b expression with axon guidance and l et-7i expression with cell adhesion, migration, and motility

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 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 c

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

The 29.5 kb <i>APOBEC3B</i> Deletion Polymorphism Is Not Associated with Clinical Outcome of Breast Cancer

<div><p>Increased <i>APOBEC3B</i> 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 <i>APOBEC3B</i>, resulting in an <i>APOBEC3A-B</i> hybrid transcript, has been associated with an increased breast cancer risk and the hypermutator phenotype. Here we evaluated whether the <i>APOBEC3B</i> 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 <i>APOBEC3B</i> 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 <i>APOBEC3B</i> 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 <i>APOBEC3B</i> copy number values and the length of metastasis-free survival (MFS; hazard ratio (HR) = 1.00, 95% confidence interval (CI) = 0.90–1.11, <i>P</i> = 0.96). Subgroup analysis by ER status also did not reveal an association between <i>APOBEC3B</i> copy number values and the length of MFS. The predictive value of the <i>APOBEC3B</i> 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 <i>APOBEC3B</i> copy number values and the overall response rate (ORR) to either tamoxifen (odds ratio (OR) = 0.88, 95% CI = 0.69–1.13, <i>P</i> = 0.31) or chemotherapy (OR = 0.97, 95% CI = 0.71–1.33, <i>P</i> = 0.87) was found. Thus, in contrast to <i>APOBEC3B</i> mRNA levels, the <i>APOBEC3B</i> deletion polymorphism has neither a prognostic nor a predictive value for breast cancer patients. Although a correlation exists between <i>APOBEC3B</i> 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 <i>APOBEC3B</i> mRNA levels.</p></div

Univariate logistic regression analysis of the overall response rate in patients treated with first-line tamoxifen and in patients treated with first-line chemotherapy for recurrent disease.

<p>Univariate logistic regression analysis of the overall response rate in patients treated with first-line tamoxifen and in patients treated with first-line chemotherapy for recurrent disease.</p

Power as a function of the hazard ratio for <i>APOBEC3B</i> copy number in each of the analyzed subgroups.

<p>The dashed horizontal line crosses the curve of each subgroup at the minimal hazard ratio for which we had 80% power in our <i>APOBEC3B</i> copy number analyses.</p

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

BACKGROUND. Continued 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. METHODS. Patients 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. RESULTS. 11KT 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. CONCLUSION. The 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.</p