Increased frequency of CHD1 deletions in prostate cancers of African American men is associated with rapid disease progression without inducing homologous recombination deficiency.

We analyzed genomic data derived from the prostate cancer of African and European American men in order to identify differences that may contribute to racial disparity of outcome and that could also define novel therapeutic strategies. In addition to analyzing patient derived next generation sequencing data, we performed FISH based confirmatory studies of Chromodomain helicase DNA-binding protein 1 (CHD1) loss on prostate cancer tissue microarrays. We created CRISPR edited, CHD1 deficient prostate cancer cell lines for genomic, drug sensitivity and functional homologous recombination (HR) activity analysis. We found that subclonal deletion of CHD1 is nearly three times as frequent in prostate tumors of African American men than in men of European ancestry and it associates with rapid disease progression. We further showed that CHD1 deletion is not associated with homologous recombination deficiency associated mutational signatures in prostate cancer. In prostate cancer cell line models CHD1 deletion did not induce HR deficiency as detected by RAD51 foci formation assay or mutational signatures, which was consistent with the moderate increase of olaparib sensitivity. CHD1 deficient prostate cancer cells, however, showed higher sensitivity to talazoparib. CHD1 loss may contribute to worse outcome of prostate cancer in African American men. A deeper understanding of the interaction between CHD1 loss and PARP inhibitor sensitivity will be needed to determine the optimal use of targeted agents such as talazoparib in the context of castration resistant prostate cancer.

A subset of lung cancer cases shows robust signs of homologous recombination deficiency associated genomic mutational signatures.

PARP inhibitors are approved for the treatment of solid tumor types that frequently harbor alterations in the key homologous recombination (HR) genes, BRCA1/2. Other tumor types, such as lung cancer, may also be HR deficient, but the frequency of such cases is less well characterized. Specific DNA aberration profiles (mutational signatures) are induced by homologous recombination deficiency (HRD) and their presence can be used to assess the presence or absence of HR deficiency in a given tumor biopsy even in the absence of an observed alteration of an HR gene. We derived various HRD-associated mutational signatures from whole-genome and whole-exome sequencing data in the lung adenocarcinoma and lung squamous carcinoma cases from TCGA, and in a patient of ours with stage IVA lung cancer with exceptionally good response to platinum-based therapy, and in lung cancer cell lines. We found that a subset of the investigated cases, both with and without biallelic loss of BRCA1 or BRCA2, showed robust signs of HR deficiency. The extreme platinum responder case also showed a robust HRD-associated genomic mutational profile. HRD-associated mutational signatures were also associated with PARP inhibitor sensitivity in lung cancer cell lines. Consequently, lung cancer cases with HRD, as identified by diagnostic mutational signatures, may benefit from PARP inhibitor therapy.

Strand Orientation Bias Detector to determine the probability of FFPE sequencing artifacts.

Formalin-fixed paraffin-embedded tissue, the most common tissue specimen stored in clinical practice, presents challenges in the analysis due to formalin-induced artifacts. Here, we present Strand Orientation Bias Detector (SOBDetector), a flexible computational platform compatible with all the common somatic SNV-calling pipelines, designed to assess the probability whether a given detected mutation is an artifact. The underlying predictor mechanism is based on the posterior distribution of a Bayesian logistic regression model trained on The Cancer Genome Atlas whole exomes. SOBDetector is a freely available cross-platform program, implemented in Java 1.8.

Publisher Correction: A clonal expression biomarker associates with lung cancer mortality.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Detection of Molecular Signatures of Homologous Recombination Deficiency in Prostate Cancer with or without BRCA1/2 Mutations.

PURPOSE: Prostate cancers with mutations in genes involved in homologous recombination (HR), most commonly BRCA2, respond favorably to PARP inhibition and platinum-based chemotherapy. We investigated whether other prostate tumors that do not harbor deleterious mutations in these particular genes can similarly be deficient in HR, likely rendering those sensitive to HR-directed therapies. EXPERIMENTAL DESIGN: Homologous recombination deficiency (HRD) levels can be estimated using various mutational signatures derived from next-generation sequencing data. We used this approach on whole-genome sequencing (WGS; n = 311) and whole-exome sequencing (WES) data (n = 498) of both primary and metastatic prostate adenocarcinomas to determine whether prostate cancer cases display clear signs of HRD in somatic tumor biopsies. RESULTS: Known BRCA-deficient samples showed all previously described HRD-associated mutational signatures in the WGS data. HRD-associated mutational signatures were also detected in a subset of patients who did not harbor germline or somatic mutations in BRCA1/2 or other HR-related genes. Similar results, albeit with lower sensitivity and accuracy, were also obtained from WES data. CONCLUSIONS: These findings may expand the number of cases likely to respond to PARP inhibitor treatment. On the basis of the HR-associated mutational signatures, 5% to 8% of localized prostate cancer cases may be good candidates for PARP-inhibitor treatment (including those with BRCA1/2 mutations).

A clonal expression biomarker associates with lung cancer mortality.

An aim of molecular biomarkers is to stratify patients with cancer into disease subtypes predictive of outcome, improving diagnostic precision beyond clinical descriptors such as tumor stage1. Transcriptomic intratumor heterogeneity (RNA-ITH) has been shown to confound existing expression-based biomarkers across multiple cancer types2-6. Here, we analyze multi-region whole-exome and RNA sequencing data for 156 tumor regions from 48 patients enrolled in the TRACERx study to explore and control for RNA-ITH in non-small cell lung cancer. We find that chromosomal instability is a major driver of RNA-ITH, and existing prognostic gene expression signatures are vulnerable to tumor sampling bias. To address this, we identify genes expressed homogeneously within individual tumors that encode expression modules of cancer cell proliferation and are often driven by DNA copy-number gains selected early in tumor evolution. Clonal transcriptomic biomarkers overcome tumor sampling bias, associate with survival independent of clinicopathological risk factors, and may provide a general strategy to refine biomarker design across cancer types.

The genomic imprint of cancer therapies helps timing the formation of metastases.

A retrospective determination of the time of metastasis formation is essential for a better understanding of the evolution of oligometastatic cancer. This study was based on the hypothesis that genomic alterations induced by cancer therapies could be used to determine the temporal order of the treatment and the formation of metastases. We analysed the whole genome sequence of a primary tumour sample and three metastatic sites derived from autopsy samples from a young never-smoker lung adenocarcinoma patient with an activating EGFR mutation. Mutation detection methods were refined to accurately detect and distinguish clonal and subclonal mutations. In comparison to a panel of samples from untreated smoker or never-smoker patients, we showed that the mutagenic effect of cisplatin treatment could be specifically detected from the base substitution mutations. Metastases that arose before or after chemotherapeutic treatment could be distinguished based on the allele frequency of cisplatin-induced dinucleotide mutations. In addition, genomic rearrangements and late amplification of the EGFR gene likely induced by afatinib treatment following the acquisition of a T790M gefitinib resistance mutation provided further evidence to tie the time of metastasis formation to treatment history. The established analysis pipeline for the detection of treatment-derived mutations allows the drawing of tumour evolutionary paths based on genomic data, showing that metastases may be seeded well before they become detectable by clinical imaging.

Migrating the SNP array-based homologous recombination deficiency measures to next generation sequencing data of breast cancer.

The first genomic scar-based homologous recombination deficiency (HRD) measures were produced using SNP arrays. As array-based technology has been largely replaced by next generation sequencing approaches, it has become important to develop algorithms that derive the same type of genomic scar scores from next generation sequencing (whole exome "WXS", whole genome "WGS") data. In order to perform this analysis, we introduce here the scarHRD R package and show that using this method the SNP array-based and next generation sequencing-based derivation of HRD scores show good correlation (Pearson correlation between 0.73 and 0.87 depending on the actual HRD measure) and that the NGS-based HRD scores distinguish similarly well between BRCA mutant and BRCA wild-type cases in a cohort of triple-negative breast cancer patients of the TCGA data set.

A comprehensive survey of the mutagenic impact of common cancer cytotoxics.

BACKGROUND: Genomic mutations caused by cytotoxic agents used in cancer chemotherapy may cause secondary malignancies as well as contribute to the evolution of treatment-resistant tumour cells. The stable diploid genome of the chicken DT40 lymphoblast cell line, an established DNA repair model system, is well suited to accurately assay genomic mutations. RESULTS: We use whole genome sequencing of multiple DT40 clones to determine the mutagenic effect of eight common cytotoxics used for the treatment of millions of patients worldwide. We determine the spontaneous mutagenesis rate at 2.3 × 10(-10) per base per cell division and find that cisplatin, cyclophosphamide and etoposide induce extra base substitutions with distinct spectra. After four cycles of exposure, cisplatin induces 0.8 mutations per Mb, equivalent to the median mutational burden in common leukaemias. Cisplatin-induced mutations, including short insertions and deletions, are mainly located at sites of putative intrastrand crosslinks. We find two of the newly defined cisplatin-specific mutation types as causes of the reversion of BRCA2 mutations in emerging cisplatin-resistant tumours or cell clones. Gemcitabine, 5-fluorouracil, hydroxyurea, doxorubicin and paclitaxel have no measurable mutagenic effect. The cisplatin-induced mutation spectrum shows good correlation with cancer mutation signatures attributed to smoking and other sources of guanine-directed base damage. CONCLUSION: This study provides support for the use of cell line mutagenesis assays to validate or predict the mutagenic effect of environmental and iatrogenic exposures. Our results suggest genetic reversion due to cisplatin-induced mutations as a distinct mechanism for developing resistance.

A robust prognostic gene expression signature for early stage lung adenocarcinoma.

BACKGROUND: Stage I lung adenocarcinoma is usually not treated with adjuvant chemotherapy; however, around half of these patients do not survive 5 years. Therefore, a reliable prognostic biomarker for early stage patients would be critical to identify those most likely to benefit from early additional treatments. Several studies have searched for gene expression prognostic biomarkers for lung adenocarcinoma, but these have not yielded a widely accepted prognosticator. RESULTS: We analyzed gene expression from seven published lung adenocarcinoma cohorts for which we included only stage I and II patients who were not given adjuvant therapy. Seven genes consistently obtained statistical significance in Cox regression for overall survival. The combined signature has a weighted mean hazard ratio of 3.2 in all cohorts and 3.0 (C.I. 1.3-7.4, p < 0.01) in an independent validation cohort and is strongly correlated with previously published signatures of chromosomal instability and cell cycle progression. CONCLUSIONS: The new prognostic signature, if validated prospectively, may enable better stratification and treatment of early stage lung cancer patients.

TumorTracer: a method to identify the tissue of origin from the somatic mutations of a tumor specimen.

BACKGROUND: A substantial proportion of cancer cases present with a metastatic tumor and require further testing to determine the primary site; many of these are never fully diagnosed and remain cancer of unknown primary origin (CUP). It has been previously demonstrated that the somatic point mutations detected in a tumor can be used to identify its site of origin with limited accuracy. We hypothesized that higher accuracy could be achieved by a classification algorithm based on the following feature sets: 1) the number of nonsynonymous point mutations in a set of 232 specific cancer-associated genes, 2) frequencies of the 96 classes of single-nucleotide substitution determined by the flanking bases, and 3) copy number profiles, if available. METHODS: We used publicly available somatic mutation data from the COSMIC database to train random forest classifiers to distinguish among those tissues of origin for which sufficient data was available. We selected feature sets using cross-validation and then derived two final classifiers (with or without copy number profiles) using 80 % of the available tumors. We evaluated the accuracy using the remaining 20 %. For further validation, we assessed accuracy of the without-copy-number classifier on three independent data sets: 1669 newly available public tumors of various types, a cohort of 91 breast metastases, and a set of 24 specimens from 9 lung cancer patients subjected to multiregion sequencing. RESULTS: The cross-validation accuracy was highest when all three types of information were used. On the left-out COSMIC data not used for training, we achieved a classification accuracy of 85 % across 6 primary sites (with copy numbers), and 69 % across 10 primary sites (without copy numbers). Importantly, a derived confidence score could distinguish tumors that could be identified with 95 % accuracy (32 %/75 % of tumors with/without copy numbers) from those that were less certain. Accuracy in the independent data sets was 46 %, 53 % and 89 % respectively, similar to the accuracy expected from the training data. CONCLUSIONS: Identification of primary site from point mutation and/or copy number data may be accurate enough to aid clinical diagnosis of cancers of unknown primary origin.

Somatic Cell Fusions Reveal Extensive Heterogeneity in Basal-like Breast Cancer.

Basal-like and luminal breast tumors have distinct clinical behavior and molecular profiles, yet the underlying mechanisms are poorly defined. To interrogate processes that determine these distinct phenotypes and their inheritance pattern, we generated somatic cell fusions and performed integrated genetic and epigenetic (DNA methylation and chromatin) profiling. We found that the basal-like trait is generally dominant and is largely defined by epigenetic repression of luminal transcription factors. Definition of super-enhancers highlighted a core program common in luminal cells but a high degree of heterogeneity in basal-like breast cancers that correlates with clinical outcome. We also found that protein extracts of basal-like cells are sufficient to induce a luminal-to-basal phenotypic switch, implying a trigger of basal-like autoregulatory circuits. We determined that KDM6A might be required for luminal-basal fusions, and we identified EN1, TBX18, and TCF4 as candidate transcriptional regulators of the luminal-to-basal switch. Our findings highlight the remarkable epigenetic plasticity of breast cancer cells.

Pan-cancer analysis of genomic scar signatures associated with homologous recombination deficiency suggests novel indications for existing cancer drugs.

BACKGROUND: Ovarian and triple-negative breast cancers with BRCA1 or BRCA2 loss are highly sensitive to treatment with PARP inhibitors and platinum-based cytotoxic agents and show an accumulation of genomic scars in the form of gross DNA copy number aberrations. Cancers without BRCA1 or BRCA2 loss but with accumulation of similar genomic scars also show increased sensitivity to platinum-based chemotherapy. Therefore, reliable biomarkers to identify DNA repair-deficient cancers prior to treatment may be useful for directing patients to platinum chemotherapy and possibly PARP inhibitors. Recently, three SNP array-based signatures of chromosomal instability were published that each quantitate a distinct type of genomic scar considered likely to be caused by improper DNA repair. They measure telomeric allelic imbalance (named NtAI), large scale transition (named LST), and loss of heterozygosity (named HRD-LOH), and it is suggested that these signatures may act as biomarkers for the state of DNA repair deficiency in a given cancer. RESULTS: We explored the pan-cancer distribution of scores of the three signatures utilizing a panel of 5371 tumors representing 15 cancer types from The Cancer Genome Atlas, and found a good correlation between scores of the three signatures (Spearman's ρ 0.73-0.87). In addition we found that cancer types ordinarily receiving platinum as standard of care have higher median scores of all three signatures. Interestingly, we also found that smaller subpopulations of high-scoring tumors exist in most cancer types, including those for which platinum chemotherapy is not standard therapy. CONCLUSIONS: Within several cancer types that are not ordinarily treated with platinum chemotherapy, we identified tumors with high levels of the three genomic biomarkers. These tumors represent identifiable subtypes of patients which may be strong candidates for clinical trials with PARP inhibitors or platinum-based chemotherapeutic regimens.

Relationship of postoperative thrombocytosis and survival of patients with colorectal cancer.

INTRODUCTION: Thrombocytosis accompanying solid tumors and predicting the prognosis of malignant tumors has been the subject of intensive research lately. Reports so far have evaluated the role of preoperative platelet count. In our present study we looked at the effect of tumor removal on platelet count and the predictive power of postoperative thrombocytosis on the survival of patients with colorectal cancer (CRC). METHODS: We retrospectively evaluated the clinical and histopathological data of 336 patients operated due to CRC between 2001 and 2011. Thrombocytosis was defined as a platelet count exceeding 400 × 10(3)/µL. Preoperative platelet count was compared with the value measured 1 month postoperatively. RESULTS: The platelet count significantly decreased after the removal of the primary tumor (paired Wilcoxon test p < 0.001). In univariate analysis preoperative thrombocytosis was a significant marker of overall survival (OS) with HR 2.2, p < 0.001 while the postoperative thrombocytosis was nearly significant with HR = 1.59, p = 0.087. In multivariate setting, when corrected for location, stage, tumor size and controlling for gender and age (> 65 years vs. ≤ 65 years), both pre- and postoperative thrombocytosis were significant independent prognostic markers with HR 1.80, p = 0.20 and HR = 1.98, p = 0.018, respectively. DISCUSSION AND CONCLUSION: Although the pathomechanism of thrombocytosis related to solid tumors is not known the decrease of platelet count after the removal of the primary tumor raises the possibility that the tumor may play an active role in the development of thrombocytosis. Furthermore, the observation of postoperative thrombocytosis with significant worse outcome underlines the predictive power of elevated platelet count.

Thrombocytosis of Liver Metastasis from Colorectal Cancer as Predictive Factor.

There is increasing evidence that thrombocytosis is associated with tumor invasion and metastasis formation. It was shown in several solid tumor types that thrombocytosis prognosticates cancer progression. The aim of this study was to evaluate preoperative thrombocytosis as a potential prognostic biomarker in isolated metastases, in patients with liver metastasis of colorectal cancer (mCRC). Clinicopathological data of 166 patients with mCRC who had surgical resection between 2001 and 2011 were collected retrospectively. All primary tumors have been already resected. The platelet count was evaluated based on the standard preoperative blood profile. The patients were followed-up on average for 28 months. Overall survival (OS) of patients with thrombocytosis was significantly worse both in univariate (HR = 3.00, p = 0.03) and in multivariate analysis (HR = 4.68, p = 0.056) when adjusted for gender, age, tumor size and surgical margin. Thrombocytosis was also a good prognosticator of disease-free survival (DFS) with HR = 2.7, p = 0.018 and nearly significant in multivariate setting (HR = 2.26, p = 0.073). The platelet count is a valuable prognostic marker for the survival in patients with mCRC.

Thrombocytosis portends adverse prognostic significance in patients with stage II colorectal carcinoma.

Thrombocytosis portends adverse prognostic significance in many types of cancers including ovarian and lung carcinoma. In this study, we determined the prevalence and prognostic significance of thrombocytosis (defined as platelet count in excess of 400 K/µl) in patients with colorectal cancer. We performed a retrospective analysis of 310 consecutive patients diagnosed at our institution between 2004 and 2013. The patients (48.7% male and 51.3% female) had a mean age of 69.9 years (+/- 12.7 years) at diagnosis. Thrombocytosis was found in a total of 25 patients, with a higher incidence in those with stage III and IV disease (14.4% of patients). Although the mean platelet count increased with the depth of tumor invasion (pT), its values remained within normal limits in the whole patient cohort. No patient with stage I cancer (n=57) had elevated platelet count at diagnosis. By contrast, five of the 78 patients (6.4%) with stage II cancer showed thrombocytosis, and four of these patients showed early recurrence and/or metastatic disease, resulting in shortened survival (they died within one year after surgery). The incidence of thrombocytosis increased to 12.2% and 20.6%, respectively, in patients with stage III and IV disease. The overall survival rate of patients with thrombocytosis was lower than those without thrombocytosis in the stage II and III disease groups, but this difference disappeared in patients with stage IV cancer who did poorly regardless of their platelet count. We concluded that thrombocytosis at diagnosis indicates adverse clinical outcome in colorectal cancer patients with stage II or III disease. This observation is especially intriguing in stage II patients because the clinical management of these patients is controversial. If our data are confirmed in larger studies, stage II colon cancer patients with thrombocytosis should be upstaged and treated as stage III/IV disease patients.

The genome of the chicken DT40 bursal lymphoma cell line.

The chicken DT40 cell line is a widely used model system in the study of multiple cellular processes due to the efficiency of homologous gene targeting. The cell line was derived from a bursal lymphoma induced by avian leukosis virus infection. In this study we characterized the genome of the cell line using whole genome shotgun sequencing and single nucleotide polymorphism array hybridization. The results indicate that wild-type DT40 has a relatively normal karyotype, except for whole chromosome copy number gains, and no karyotype variability within stocks. In a comparison to two domestic chicken genomes and the Gallus gallus reference genome, we found no unique mutational processes shaping the DT40 genome except for a mild increase in insertion and deletion events, particularly deletions at tandem repeats. We mapped coding sequence mutations that are unique to the DT40 genome; mutations inactivating the PIK3R1 and ATRX genes likely contributed to the oncogenic transformation. In addition to a known avian leukosis virus integration in the MYC gene, we detected further integration sites that are likely to de-regulate gene expression. The new findings support the hypothesis that DT40 is a typical transformed cell line with a relatively intact genome; therefore, it is well-suited to the role of a model system for DNA repair and related processes. The sequence data generated by this study, including a searchable de novo genome assembly and annotated lists of mutated genes, will support future research using this cell line.

Lactate dehydrogenase as a selection criterion for ipilimumab treatment in metastatic melanoma.

INTRODUCTION: Ipilimumab, a cytotoxic T lymphocyte-associated antigen-4 blocking antibody, has improved overall survival (OS) in metastatic melanoma in phase III trials. However, about 80 % of patients fail to respond, and no predictive markers for benefit from therapy have been identified. We analysed a 'real world' population of patients treated with ipilimumab to identify markers for treatment benefit. METHODS: Patients with advanced cutaneous melanoma were treated in the Netherlands (NL) and the United Kingdom (UK) with ipilimumab at 3 mg/kg. Baseline characteristics and peripheral blood parameters were assessed, and patients were monitored for the occurrence of adverse events and outcomes. RESULTS: A total of 166 patients were treated in the Netherlands. Best overall response and disease control rates were 17 and 35 %, respectively. Median follow-up was 17.9 months, with a median progression-free survival of 2.9 months. Median OS was 7.5 months, and OS at 1 year was 37.8 % and at 2 years was 22.9 %. In a multivariate model, baseline serum lactate dehydrogenase (LDH) was demonstrated to be the strongest predictive factor for OS. These findings were validated in an independent cohort of 64 patients from the UK. CONCLUSION: In both the NL and UK cohorts, long-term benefit of ipilimumab treatment was unlikely for patients with baseline serum LDH greater than twice the upper limit of normal. In the absence of prospective data, clinicians treating melanoma may wish to consider the data presented here to guide patient selection for ipilimumab therapy.

[Evaluation of thrombocytosis as predictive factor in colorectal cancer]. / A thrombocytosis mint prediktív faktor értékelése colorectalis tumorokban.

UNLABELLED: INTRODUCTION/AIM OF THE STUDY: Preoperative thrombocytosis proved to be a negative prognostic factor in several solid tumor. However, there is still debate in the literature regarding colorectal cancer. The aim of our study was to examine whether thrombocytosis is an independent risk factor for metastasis development and predictor of survival in colorectal cancer. MATERIALS AND METHODS: Clinicopathological data of 336 patients with colorectal cancer (CRC) and 118 patients with liver metastasis of colorectal cancer (mCRC) who had operation between 2001 and 2011 were collected retrospectively. Thrombocytosis was defined as 400 G/L < platelet count. Disease-free survival (DFS) and overall survival (OS) were determined with Kaplan-Meier method supported by log-rank test. RESULTS: Both in the CRC and the mCRC group OS was significantly shorter in patients who had elevated platelet count (HR = 2.2, p < 0.001 and HR = 2.9, p = 0.018, respectively). Multivariate analysis confirmed that elevated platelet count was an independent prognostic factor of both CRC (HR = 1.7, p = 0.035) and mCRC (HR = 3.1, p = 0.017). DFS was significantly shorter in patients with elevated platelet count in the CRC group (HR = 2.0, p = 0.011). DISCUSSION: The platelet count is a valuable and cheap prognostic marker for the prediction of survival in patients both with CRC and mCRC.

Biasogram: visualization of confounding technical bias in gene expression data.

Gene expression profiles of clinical cohorts can be used to identify genes that are correlated with a clinical variable of interest such as patient outcome or response to a particular drug. However, expression measurements are susceptible to technical bias caused by variation in extraneous factors such as RNA quality and array hybridization conditions. If such technical bias is correlated with the clinical variable of interest, the likelihood of identifying false positive genes is increased. Here we describe a method to visualize an expression matrix as a projection of all genes onto a plane defined by a clinical variable and a technical nuisance variable. The resulting plot indicates the extent to which each gene is correlated with the clinical variable or the technical variable. We demonstrate this method by applying it to three clinical trial microarray data sets, one of which identified genes that may have been driven by a confounding technical variable. This approach can be used as a quality control step to identify data sets that are likely to yield false positive results.