Proteomic analysis of brain metastatic lung adenocarcinoma reveals intertumoral heterogeneity and specific alterations associated with the timing of brain metastases.

BACKGROUND: Brain metastases are associated with considerable negative effects on patients' outcome in lung adenocarcinoma (LADC). Here, we investigated the proteomic landscape of primary LADCs and their corresponding brain metastases. MATERIALS AND METHODS: Proteomic profiling was conducted on 20 surgically resected primary and brain metastatic LADC samples via label-free shotgun proteomics. After sample processing, peptides were analyzed using an Ultimate 3000 pump coupled to a QExactive HF-X mass spectrometer. Raw data were searched using PD 2.4. Further data analyses were carried out using Perseus, RStudio and GraphPad Prism. Proteomic data were correlated with clinical and histopathological parameters and the timing of brain metastases. Mass spectrometry-based proteomic data are available via ProteomeXchange with identifier PXD027259. RESULTS: Out of the 6821 proteins identified and quantified, 1496 proteins were differentially expressed between primary LADCs and corresponding brain metastases. Pathways associated with the immune system, cell-cell/matrix interactions and migration were predominantly activated in the primary tumors, whereas pathways related to metabolism, translation or vesicle formation were overrepresented in the metastatic tumors. When comparing fast- versus slow-progressing patients, we found 454 and 298 differentially expressed proteins in the primary tumors and brain metastases, respectively. Metabolic reprogramming and ribosomal activity were prominently up-regulated in the fast-progressing patients (versus slow-progressing individuals), whereas expression of cell-cell interaction- and immune system-related pathways was reduced in these patients and in those with multiple brain metastases. CONCLUSIONS: This is the first comprehensive proteomic analysis of paired primary tumors and brain metastases of LADC patients. Our data suggest a malfunction of cellular attachment and an increase in ribosomal activity in LADC tissue, promoting brain metastasis. The current study provides insights into the biology of LADC brain metastases and, moreover, might contribute to the development of personalized follow-up strategies in LADC.

Breast cancer brain metastases show increased levels of genomic aberration-based homologous recombination deficiency scores relative to their corresponding primary tumors.

Background: Based on its mechanism of action, PARP inhibitor therapy is expected to benefit mainly tumor cases with homologous recombination deficiency (HRD). Therefore, identification of tumor types with increased HRD is important for the optimal use of this class of therapeutic agents. HRD levels can be estimated using various mutational signatures from next generation sequencing data and we used this approach to determine whether breast cancer brain metastases show altered levels of HRD scores relative to their corresponding primary tumor. Patients and methods: We used a previously published next generation sequencing dataset of 21 matched primary breast cancer/brain metastasis pairs to derive the various mutational signatures/HRD scores strongly associated with HRD. We also carried out the myChoice HRD analysis on an independent cohort of 17 breast cancer patients with matched primary/brain metastasis pairs. Results: All of the mutational signatures indicative of HRD showed a significant increase in the brain metastases relative to their matched primary tumor in the previously published whole exome sequencing dataset. In the independent validation cohort, the myChoice HRD assay showed an increased level in 87.5% of the brain metastases relative to the primary tumor, with 56% of brain metastases being HRD positive according to the myChoice criteria. Conclusions: The consistent observation that brain metastases of breast cancer tend to have higher HRD measures may raise the possibility that brain metastases may be more sensitive to PARP inhibitor treatment. This observation warrants further investigation to assess whether this increase is common to other metastatic sites as well, and whether clinical trials should adjust their strategy in the application of HRD measures for the prioritization of patients for PARP inhibitor therapy.

Overexpression of BLM promotes DNA damage and increased sensitivity to platinum salts in triple-negative breast and serous ovarian cancers.

Background: Platinum-based therapy is an effective treatment for a subset of triple-negative breast cancer and ovarian cancer patients. In order to increase response rate and decrease unnecessary use, robust biomarkers that predict response to therapy are needed. Patients and methods: We performed an integrated genomic approach combining differential analysis of gene expression and DNA copy number in sensitive compared with resistant triple-negative breast cancers in two independent neoadjuvant cisplatin-treated cohorts. Functional relevance of significant hits was investigated in vitro by overexpression, knockdown and targeted inhibitor treatment. Results: We identified two genes, the Bloom helicase (BLM) and Fanconi anemia complementation group I (FANCI), that have both increased DNA copy number and gene expression in the platinum-sensitive cases. Increased level of expression of these two genes was also associated with platinum but not with taxane response in ovarian cancer. As a functional validation, we found that overexpression of BLM promotes DNA damage and induces sensitivity to cisplatin but has no effect on paclitaxel sensitivity. Conclusions: A biomarker based on the expression levels of the BLM and FANCI genes is a potential predictor of platinum sensitivity in triple-negative breast cancer and ovarian cancer.

Intratumoural evolutionary landscape of high-risk prostate cancer: the PROGENY study of genomic and immune parameters.

BACKGROUND: Intratumoural heterogeneity (ITH) is well recognised in prostate cancer (PC), but its role in high-risk disease is uncertain. A prospective, single-arm, translational study using targeted multiregion prostate biopsies was carried out to study genomic and T-cell ITH in clinically high-risk PC aiming to identify drivers and potential therapeutic strategies. PATIENTS AND METHODS: Forty-nine men with elevated prostate-specific antigen and multiparametric-magnetic resonance imaging detected PC underwent image-guided multiregion transperineal biopsy. Seventy-nine tumour regions from 25 patients with PC underwent sequencing, analysis of mutations, copy number and neoepitopes combined with tumour infiltrating T-cell subset quantification. RESULTS: We demonstrated extensive somatic nucleotide variation and somatic copy number alteration heterogeneity in high-risk PC. Overall, the mutational burden was low (0.93/Megabase), but two patients had hypermutation, with loss of mismatch repair (MMR) proteins, MSH2 and MSH6. Somatic copy number alteration burden was higher in patients with metastatic hormone-naive PC (mHNPC) than in those with high-risk localised PC (hrlPC), independent of Gleason grade. Mutations were rarely ubiquitous and mutational frequencies were similar for mHNPC and hrlPC patients. Enrichment of focal 3q26.2 and 3q21.3, regions containing putative metastasis drivers, was seen in mHNPC patients. We found evidence of parallel evolution with three separate clones containing activating mutations of ß-catenin in a single patient. We demonstrated extensive intratumoural and intertumoural T-cell heterogeneity and high inflammatory infiltrate in the MMR-deficient (MMRD) patients and the patient with parallel evolution of ß-catenin. Analysis of all patients with activating Wnt/ß-catenin mutations demonstrated a low CD8+/FOXP3+ ratio, a potential surrogate marker of immune evasion. CONCLUSIONS: The PROGENY (PROstate cancer GENomic heterogeneitY) study provides a diagnostic platform suitable for studying tumour ITH. Genetic aberrations in clinically high-risk PC are associated with altered patterns of immune infiltrate in tumours. Activating mutations of Wnt/ß-catenin signalling pathway or MMRD could be considered as potential biomarkers for immunomodulation therapies. CLINICAL TRIALS.GOV IDENTIFIER: NCT02022371.

Loss of BRCA1 or BRCA2 markedly increases the rate of base substitution mutagenesis and has distinct effects on genomic deletions.

This corrects the article DOI: 10.1038/onc.2016.243.

Fast and accurate mutation detection in whole genome sequences of multiple isogenic samples with IsoMut.

BACKGROUND: Detection of somatic mutations is one of the main goals of next generation DNA sequencing. A wide range of experimental systems are available for the study of spontaneous or environmentally induced mutagenic processes. However, most of the routinely used mutation calling algorithms are not optimised for the simultaneous analysis of multiple samples, or for non-human experimental model systems with no reliable databases of common genetic variations. Most standard tools either require numerous in-house post filtering steps with scarce documentation or take an unpractically long time to run. To overcome these problems, we designed the streamlined IsoMut tool which can be readily adapted to experimental scenarios where the goal is the identification of experimentally induced mutations in multiple isogenic samples. METHODS: Using 30 isogenic samples, reliable cohorts of validated mutations were created for testing purposes. Optimal values of the filtering parameters of IsoMut were determined in a thorough and strict optimization procedure based on these test sets. RESULTS: We show that IsoMut, when tuned correctly, decreases the false positive rate compared to conventional tools in a 30 sample experimental setup; and detects not only single nucleotide variations, but short insertions and deletions as well. IsoMut can also be run more than a hundred times faster than the most precise state of art tool, due its straightforward and easily understandable filtering algorithm. CONCLUSIONS: IsoMut has already been successfully applied in multiple recent studies to find unique, treatment induced mutations in sets of isogenic samples with very low false positive rates. These types of studies provide an important contribution to determining the mutagenic effect of environmental agents or genetic defects, and IsoMut turned out to be an invaluable tool in the analysis of such data.

Loss of BRCA1 or BRCA2 markedly increases the rate of base substitution mutagenesis and has distinct effects on genomic deletions.

Loss-of-function mutations in the BRCA1 and BRCA2 genes increase the risk of cancer. Owing to their function in homologous recombination repair, much research has focused on the unstable genomic phenotype of BRCA1/2 mutant cells manifest mainly as large-scale rearrangements. We used whole-genome sequencing of multiple isogenic chicken DT40 cell clones to precisely determine the consequences of BRCA1/2 loss on all types of genomic mutagenesis. Spontaneous base substitution mutation rates increased sevenfold upon the disruption of either BRCA1 or BRCA2, and the arising mutation spectra showed strong and specific correlation with a mutation signature associated with BRCA1/2 mutant tumours. To model endogenous alkylating damage, we determined the mutation spectrum caused by methyl methanesulfonate (MMS), and showed that MMS also induces more base substitution mutations in BRCA1/2-deficient cells. Spontaneously arising and MMS-induced insertion/deletion mutations and large rearrangements were also more common in BRCA1/2 mutant cells compared with the wild-type control. A difference in the short deletion phenotypes of BRCA1 and BRCA2 suggested distinct roles for the two proteins in the processing of DNA lesions, as BRCA2 mutants contained more short deletions, with a wider size distribution, which frequently showed microhomology near the breakpoints resembling repair by non-homologous end joining. An increased and prolonged gamma-H2AX signal in MMS-treated BRCA1/2 cells suggested an aberrant processing of stalled replication forks as the cause of increased mutagenesis. The high rate of base substitution mutagenesis demonstrated by our experiments is likely to significantly contribute to the oncogenic effect of the inactivation of BRCA1 or BRCA2.

SETD2 loss-of-function promotes renal cancer branched evolution through replication stress and impaired DNA repair.

Defining mechanisms that generate intratumour heterogeneity and branched evolution may inspire novel therapeutic approaches to limit tumour diversity and adaptation. SETD2 (Su(var), Enhancer of zeste, Trithorax-domain containing 2) trimethylates histone-3 lysine-36 (H3K36me3) at sites of active transcription and is mutated in diverse tumour types, including clear cell renal carcinomas (ccRCCs). Distinct SETD2 mutations have been identified in spatially separated regions in ccRCC, indicative of intratumour heterogeneity. In this study, we have addressed the consequences of SETD2 loss-of-function through an integrated bioinformatics and functional genomics approach. We find that bi-allelic SETD2 aberrations are not associated with microsatellite instability in ccRCC. SETD2 depletion in ccRCC cells revealed aberrant and reduced nucleosome compaction and chromatin association of the key replication proteins minichromosome maintenance complex component (MCM7) and DNA polymerase δ hindering replication fork progression, and failure to load lens epithelium-derived growth factor and the Rad51 homologous recombination repair factor at DNA breaks. Consistent with these data, we observe chromosomal breakpoint locations are biased away from H3K36me3 sites in SETD2 wild-type ccRCCs relative to tumours with bi-allelic SETD2 aberrations and that H3K36me3-negative ccRCCs display elevated DNA damage in vivo. These data suggest a role for SETD2 in maintaining genome integrity through nucleosome stabilization, suppression of replication stress and the coordination of DNA repair.

Glioblastoma adaptation traced through decline of an IDH1 clonal driver and macro-evolution of a double-minute chromosome.

BACKGROUND: Glioblastoma (GBM) is the most common malignant brain cancer occurring in adults, and is associated with dismal outcome and few therapeutic options. GBM has been shown to predominantly disrupt three core pathways through somatic aberrations, rendering it ideal for precision medicine approaches. METHODS: We describe a 35-year-old female patient with recurrent GBM following surgical removal of the primary tumour, adjuvant treatment with temozolomide and a 3-year disease-free period. Rapid whole-genome sequencing (WGS) of three separate tumour regions at recurrence was carried out and interpreted relative to WGS of two regions of the primary tumour. RESULTS: We found extensive mutational and copy-number heterogeneity within the primary tumour. We identified a TP53 mutation and two focal amplifications involving PDGFRA, KIT and CDK4, on chromosomes 4 and 12. A clonal IDH1 R132H mutation in the primary, a known GBM driver event, was detectable at only very low frequency in the recurrent tumour. After sub-clonal diversification, evidence was found for a whole-genome doubling event and a translocation between the amplified regions of PDGFRA, KIT and CDK4, encoded within a double-minute chromosome also incorporating miR26a-2. The WGS analysis uncovered progressive evolution of the double-minute chromosome converging on the KIT/PDGFRA/PI3K/mTOR axis, superseding the IDH1 mutation in dominance in a mutually exclusive manner at recurrence, consequently the patient was treated with imatinib. Despite rapid sequencing and cancer genome-guided therapy against amplified oncogenes, the disease progressed, and the patient died shortly after. CONCLUSION: This case sheds light on the dynamic evolution of a GBM tumour, defining the origins of the lethal sub-clone, the macro-evolutionary genomic events dominating the disease at recurrence and the loss of a clonal driver. Even in the era of rapid WGS analysis, cases such as this illustrate the significant hurdles for precision medicine success.

Sequenza: allele-specific copy number and mutation profiles from tumor sequencing data.

BACKGROUND: Exome or whole-genome deep sequencing of tumor DNA along with paired normal DNA can potentially provide a detailed picture of the somatic mutations that characterize the tumor. However, analysis of such sequence data can be complicated by the presence of normal cells in the tumor specimen, by intratumor heterogeneity, and by the sheer size of the raw data. In particular, determination of copy number variations from exome sequencing data alone has proven difficult; thus, single nucleotide polymorphism (SNP) arrays have often been used for this task. Recently, algorithms to estimate absolute, but not allele-specific, copy number profiles from tumor sequencing data have been described. MATERIALS AND METHODS: We developed Sequenza, a software package that uses paired tumor-normal DNA sequencing data to estimate tumor cellularity and ploidy, and to calculate allele-specific copy number profiles and mutation profiles. We applied Sequenza, as well as two previously published algorithms, to exome sequence data from 30 tumors from The Cancer Genome Atlas. We assessed the performance of these algorithms by comparing their results with those generated using matched SNP arrays and processed by the allele-specific copy number analysis of tumors (ASCAT) algorithm. RESULTS: Comparison between Sequenza/exome and SNP/ASCAT revealed strong correlation in cellularity (Pearson's r = 0.90) and ploidy estimates (r = 0.42, or r = 0.94 after manual inspecting alternative solutions). This performance was noticeably superior to previously published algorithms. In addition, in artificial data simulating normal-tumor admixtures, Sequenza detected the correct ploidy in samples with tumor content as low as 30%. CONCLUSIONS: The agreement between Sequenza and SNP array-based copy number profiles suggests that exome sequencing alone is sufficient not only for identifying small scale mutations but also for estimating cellularity and inferring DNA copy number aberrations.

The comparison of thrombocytosis and platelet-lymphocyte ratio as potential prognostic markers in colorectal cancer.

The aim of the present study was to analyse the preoperative platelet count and the platelet-lymphocyte ratio (PLR) in patients with colorectal cancer (CRC) of different stages and with hepatic metastasis of CRC (mCRC) and to compare these factors as potential prognostic markers. Clinicopathological data of 10 years were collected retrospectively from 336 patients with CRC and 118 patients with mCRC. Both in the CRC and the mCRC group overall survival (OS) was significantly worse in patients who had elevated platelet count (hazard ratio [HR] = 2.2, p < 0.001 and HR = 2.9, p = 0.018, respectively). Multivariate analysis indicated that elevated platelet count was an independent prognostic factor of CRC (HR = 1.7, p = 0.035) and mCRC (HR = 3.1, p = 0.017). Disease-free survival (DFS) was significantly worse in patients with elevated platelet count in the CRC group (HR = 2.0, p = 0.011). In the multivariate analysis the PLR was not a prognostic factor in either of the two cohorts (HR = 0.92, p < 0.001 and HR = 0.89, p = 0.789, respectively). The platelet count is a valuable prognostic marker for the survival in patients both with CRC and mCRC while the PLR is not prognostic in either group.

Down-regulation of the desmosomal cadherin desmocollin 3 in human breast cancer.

In previous studies using cDNA microarray analysis, we have identified an expressed sequence tag which is consistently down-regulated in six human breast tumor cell lines. In the current study, we have determined this tag to be part of the mRNA sequence of human desmocollin 3, a member of the cadherin superfamily of proteins and an integral component of desmosomes. Desmosomes are sites of adhesion between adjacent cells in layers of epithelia, as well as in some non-epithelial tissues, and play an important role in the maintenance of tissue structure. Northern analysis, quantitative real-time polymerase chain reaction assay and Western blot analysis showed that desmocollin 3 is present in normal and immortalized human mammary epithelial cells, but consistently exhibits a significant, and often complete, down-regulation in breast cancer cell lines and primary breast tumors, both at the mRNA and protein levels.

Mutual information analysis as a tool to assess the role of aneuploidy in the generation of cancer-associated differential gene expression patterns.

Most human tumors are characterized by: (1) an aberrant set of chromosomes, a state termed aneuploidy; (2) an aberrant gene expression pattern; and (3) an aberrant phenotype of uncontrolled growth. One of the goals of cancer research is to establish causative relationships between these three important characteristics. In this paper we were searching for evidence that aneuploidy is a major cause of differential gene expression. We describe how mutual information analysis of cancer-associated gene expression patterns could be exploited to answer this question. In addition to providing general guidelines, we have applied the proposed analysis to a recently published breast cancer-associated gene expression matrix. The results derived from this particular data set provided preliminary evidence that mutual information analysis may become a useful tool to investigate the link between differential gene expression and aneuploidy.