PTBP1 promotes hematopoietic stem cell maintenance and red blood cell development by ensuring sufficient availability of ribosomal constituents.

Ribosomopathies constitute a range of disorders associated with defective protein synthesis mainly affecting hematopoietic stem cells (HSCs) and erythroid development. Here, we demonstrate that deletion of poly-pyrimidine-tract-binding protein 1 (PTBP1) in the hematopoietic compartment leads to the development of a ribosomopathy-like condition. Specifically, loss of PTBP1 is associated with decreases in HSC self-renewal, erythroid differentiation, and protein synthesis. Consistent with its function as a splicing regulator, PTBP1 deficiency results in splicing defects in hundreds of genes, and we demonstrate that the up-regulation of a specific isoform of CDC42 partly mimics the protein-synthesis defect associated with loss of PTBP1. Furthermore, PTBP1 deficiency is associated with a marked defect in ribosome biogenesis and a selective reduction in the translation of mRNAs encoding ribosomal proteins. Collectively, this work identifies PTBP1 as a key integrator of ribosomal functions and highlights the broad functional repertoire of RNA-binding proteins.

Functional Genomic Identification of Predictors of Sensitivity and Mechanisms of Resistance to Multivalent Second-Generation TRAIL-R2 Agonists.

Multivalent second-generation TRAIL-R2 agonists are currently in late preclinical development and early clinical trials. Herein, we use a representative second-generation agent, MEDI3039, to address two major clinical challenges facing these agents: lack of predictive biomarkers to enable patient selection and emergence of resistance. Genome-wide CRISPR knockout screens were notable for the lack of resistance mechanisms beyond the canonical TRAIL-R2 pathway (caspase-8, FADD, BID) as well as p53 and BAX in TP53 wild-type models, whereas a CRISPR activatory screen identified cell death inhibitors MCL-1 and BCL-XL as mechanisms to suppress MEDI3039-induced cell death. High-throughput drug screening failed to identify genomic alterations associated with response to MEDI3039; however, transcriptomics analysis revealed striking association between MEDI3039 sensitivity and expression of core components of the extrinsic apoptotic pathway, most notably its main apoptotic effector caspase-8 in solid tumor cell lines. Further analyses of colorectal cell lines and patient-derived xenografts identified caspase-8 expression ratio to its endogenous regulator FLIP(L) as predictive of sensitivity to MEDI3039 in several major solid tumor types and a further subset indicated by caspase-8:MCL-1 ratio. Subsequent MEDI3039 combination screening of TRAIL-R2, caspase-8, FADD, and BID knockout models with 60 compounds with varying mechanisms of action identified two inhibitor of apoptosis proteins (IAP) that exhibited strong synergy with MEDI3039 that could reverse resistance only in BID-deleted models. In summary, we identify the ratios of caspase-8:FLIP(L) and caspase-8:MCL-1 as potential predictive biomarkers for second-generation TRAIL-R2 agonists and loss of key effectors such as FADD and caspase-8 as likely drivers of clinical resistance in solid tumors.

Kinome capture sequencing of high-grade serous ovarian carcinoma reveals novel mutations in the JAK3 gene.

High-grade serous ovarian carcinoma (HGSOC) remains the deadliest form of epithelial ovarian cancer and despite major efforts little improvement in overall survival has been achieved. Identification of recurring "driver" genetic lesions has the potential to enable design of novel therapies for cancer. Here, we report on a study to find such new therapeutic targets for HGSOC using exome-capture sequencing approach targeting all kinase genes in 127 patient samples. Consistent with previous reports, the most frequently mutated gene was TP53 (97% mutation frequency) followed by BRCA1 (10% mutation frequency). The average mutation frequency of the kinase genes mutated from our panel was 1.5%. Intriguingly, after BRCA1, JAK3 was the most frequently mutated gene (4% mutation frequency). We tested the transforming properties of JAK3 mutants using the Ba/F3 cell-based in vitro functional assay and identified a novel gain-of-function mutation in the kinase domain of JAK3 (p.T1022I). Importantly, p.T1022I JAK3 mutants displayed higher sensitivity to the JAK3-selective inhibitor Tofacitinib compared to controls. For independent validation, we re-sequenced the entire JAK3 coding sequence using tagged amplicon sequencing (TAm-Seq) in 463 HGSOCs resulting in an overall somatic mutation frequency of 1%. TAm-Seq screening of CDK12 in the same population revealed a 7% mutation frequency. Our data confirms that the frequency of mutations in kinase genes in HGSOC is low and provides accurate estimates for the frequency of JAK3 and CDK12 mutations in a large well characterized cohort. Although p.T1022I JAK3 mutations are rare, our functional validation shows that if detected they should be considered as potentially actionable for therapy. The observation of CDK12 mutations in 7% of HGSOC cases provides a strong rationale for routine somatic testing, although more functional and clinical characterization is required to understand which nonsynonymous mutations alterations are associated with homologous recombination deficiency.

ERG Controls B Cell Development by Promoting Igh V-to-DJ Recombination.

B cell development depends on the coordinated expression and cooperation of several transcription factors. Here we show that the transcription factor ETS-related gene (ERG) is crucial for normal B cell development and that its deletion results in a substantial loss of bone marrow B cell progenitors and peripheral B cells, as well as a skewing of splenic B cell populations. We find that ERG-deficient B lineage cells exhibit an early developmental block at the pre-B cell stage and proliferate less. The cells fail to express the immunoglobulin heavy chain due to inefficient V-to-DJ recombination, and cells that undergo recombination display a strong bias against incorporation of distal V gene segments. Furthermore, antisense transcription at PAX5-activated intergenic repeat (PAIR) elements, located in the distal region of the Igh locus, depends on ERG. These findings show that ERG serves as a critical regulator of B cell development by ensuring efficient and balanced V-to-DJ recombination.

Principal component analysis of binary genomics data.

Motivation: Genome-wide measurements of genetic and epigenetic alterations are generating more and more high-dimensional binary data. The special mathematical characteristics of binary data make the direct use of the classical principal component analysis (PCA) model to explore low-dimensional structures less obvious. Although there are several PCA alternatives for binary data in the psychometric, data analysis and machine learning literature, they are not well known to the bioinformatics community. Results: In this article, we introduce the motivation and rationale of some parametric and nonparametric versions of PCA specifically geared for binary data. Using both realistic simulations of binary data as well as mutation, CNA and methylation data of the Genomic Determinants of Sensitivity in Cancer 1000 (GDSC1000), the methods were explored for their performance with respect to finding the correct number of components, overfit, finding back the correct low-dimensional structure, variable importance, etc. The results show that if a low-dimensional structure exists in the data, that most of the methods can find it. When assuming a probabilistic generating process is underlying the data, we recommend to use the parametric logistic PCA model, while when such an assumption is not valid and the data are considered as given, the nonparametric Gifi model is recommended. Availability: The codes to reproduce the results in this article are available at the homepage of the Biosystems Data Analysis group (www.bdagroup.nl).

Assessment of PD-L1 expression across breast cancer molecular subtypes, in relation to mutation rate, BRCA1-like status, tumor-infiltrating immune cells and survival.

To better understand the expression pattern of programmed death-ligand 1 (PD-L1) expression in different breast cancer types, we characterized PD-L1 expression in tumor and tumor-infiltrating immune cells, in relation to mutation rate, BRCA1-like status and survival. We analyzed 410 primary treatment-naive breast tumors comprising 162 estrogen receptor-positive (ER+) and HER2-, 101 HER2+ and 147 triple-negative (TN) cancers. Pathologists quantified tumor-infiltrating lymphocytes (TILs) and PD-L1 expression in tumor cells and TILs using whole slides and tissue microarray. Mutation rate was assessed by DNA sequencing, BRCA1-like status using multiplex ligation-dependent probe amplification, and immune landscape by multiplex image analyses of CD4, CD68, CD8, FOXP3, cytokeratin, and PD-L1. Half of PD-L1 scores evaluated by tissue microarray were false negatives compared to whole slide evaluations. We observed at least 1% of PD-L1-positive (PD-L1+) cells in 53.1% of ER+HER2-, 73.3% of HER2+, and 84.4% of TN tumors. PD-L1 expression was higher in ductal compared to lobular carcinomas, also within ER+HER2- tumors (p = 0.04). High PD-L1+ TILs score (> 50%) was independently associated with better outcome in TN tumors (HR = 0.27; 95%CI = 0.10-0.69). Within TN tumors, PD-L1 and TIL scores showed a modest but significant positive association with the number of silent mutations, but no association with BRCA1-like status. Multiplex image analyses indicated that PD-L1 is expressed on multiple immune cells (CD68+ macrophages, CD4+, FOXP3+, and CD8+ T cells) in the breast tumor microenvironment, independent of the PD-L1 status of the tumor cells. We found no evidence that levels of PD-L1+ TILs in TN breast cancer are driven by high mutation rate or BRCA1-like status.

iTOP: inferring the topology of omics data.

Motivation: In biology, we are often faced with multiple datasets recorded on the same set of objects, such as multi-omics and phenotypic data of the same tumors. These datasets are typically not independent from each other. For example, methylation may influence gene expression, which may, in turn, influence drug response. Such relationships can strongly affect analyses performed on the data, as we have previously shown for the identification of biomarkers of drug response. Therefore, it is important to be able to chart the relationships between datasets. Results: We present iTOP, a methodology to infer a topology of relationships between datasets. We base this methodology on the RV coefficient, a measure of matrix correlation, which can be used to determine how much information is shared between two datasets. We extended the RV coefficient for partial matrix correlations, which allows the use of graph reconstruction algorithms, such as the PC algorithm, to infer the topologies. In addition, since multi-omics data often contain binary data (e.g. mutations), we also extended the RV coefficient for binary data. Applying iTOP to pharmacogenomics data, we found that gene expression acts as a mediator between most other datasets and drug response: only proteomics clearly shares information with drug response that is not present in gene expression. Based on this result, we used TANDEM, a method for drug response prediction, to identify which variables predictive of drug response were distinct to either gene expression or proteomics. Availability and implementation: An implementation of our methodology is available in the R package iTOP on CRAN. Additionally, an R Markdown document with code to reproduce all figures is provided as Supplementary Material. Supplementary information: Supplementary data are available at Bioinformatics online.

Loss of functional BAP1 augments sensitivity to TRAIL in cancer cells.

Malignant mesothelioma (MM) is poorly responsive to systemic cytotoxic chemotherapy and invariably fatal. Here we describe a screen of 94 drugs in 15 exome-sequenced MM lines and the discovery of a subset defined by loss of function of the nuclear deubiquitinase BRCA associated protein-1 (BAP1) that demonstrate heightened sensitivity to TRAIL (tumour necrosis factor-related apoptosis-inducing ligand). This association is observed across human early passage MM cultures, mouse xenografts and human tumour explants. We demonstrate that BAP1 deubiquitinase activity and its association with ASXL1 to form the Polycomb repressive deubiquitinase complex (PR-DUB) impacts TRAIL sensitivity implicating transcriptional modulation as an underlying mechanism. Death receptor agonists are well-tolerated anti-cancer agents demonstrating limited therapeutic benefit in trials without a targeting biomarker. We identify BAP1 loss-of-function mutations, which are frequent in MM, as a potential genomic stratification tool for TRAIL sensitivity with immediate and actionable therapeutic implications.

The BRCA1ness signature is associated significantly with response to PARP inhibitor treatment versus control in the I-SPY 2 randomized neoadjuvant setting.

BACKGROUND: Patients with BRCA1-like tumors correlate with improved response to DNA double-strand break-inducing therapy. A gene expression-based classifier was developed to distinguish between BRCA1-like and non-BRCA1-like tumors. We hypothesized that these tumors may also be more sensitive to PARP inhibitors than standard treatments. METHODS: A diagnostic gene expression signature (BRCA1ness) was developed using a centroid model with 128 triple-negative breast cancer samples from the EU FP7 RATHER project. This BRCA1ness signature was then tested in HER2-negative patients (n = 116) from the I-SPY 2 TRIAL who received an oral PARP inhibitor veliparib in combination with carboplatin (V-C), or standard chemotherapy alone. We assessed the association between BRCA1ness and pathologic complete response in the V-C and control arms alone using Fisher's exact test, and the relative performance between arms (biomarker × treatment interaction, likelihood ratio p < 0.05) using a logistic model and adjusting for hormone receptor status (HR). RESULTS: We developed a gene expression signature to identify BRCA1-like status. In the I-SPY 2 neoadjuvant setting the BRCA1ness signature associated significantly with response to V-C (p = 0.03), but not in the control arm (p = 0.45). We identified a significant interaction between BRCA1ness and V-C (p = 0.023) after correcting for HR. CONCLUSIONS: A genomic-based BRCA1-like signature was successfully translated to an expression-based signature (BRC1Aness). In the I-SPY 2 neoadjuvant setting, we determined that the BRCA1ness signature is capable of predicting benefit of V-C added to standard chemotherapy compared to standard chemotherapy alone. TRIAL REGISTRATION: I-SPY 2 TRIAL beginning December 31, 2009: Neoadjuvant and Personalized Adaptive Novel Agents to Treat Breast Cancer (I-SPY 2), NCT01042379 .

Genomic Determinants of Protein Abundance Variation in Colorectal Cancer Cells.

Assessing the impact of genomic alterations on protein networks is fundamental in identifying the mechanisms that shape cancer heterogeneity. We have used isobaric labeling to characterize the proteomic landscapes of 50 colorectal cancer cell lines and to decipher the functional consequences of somatic genomic variants. The robust quantification of over 9,000 proteins and 11,000 phosphopeptides on average enabled the de novo construction of a functional protein correlation network, which ultimately exposed the collateral effects of mutations on protein complexes. CRISPR-cas9 deletion of key chromatin modifiers confirmed that the consequences of genomic alterations can propagate through protein interactions in a transcript-independent manner. Lastly, we leveraged the quantified proteome to perform unsupervised classification of the cell lines and to build predictive models of drug response in colorectal cancer. Overall, we provide a deep integrative view of the functional network and the molecular structure underlying the heterogeneity of colorectal cancer cells.

TANDEM: a two-stage approach to maximize interpretability of drug response models based on multiple molecular data types.

MOTIVATION: Clinical response to anti-cancer drugs varies between patients. A large portion of this variation can be explained by differences in molecular features, such as mutation status, copy number alterations, methylation and gene expression profiles. We show that the classic approach for combining these molecular features (Elastic Net regression on all molecular features simultaneously) results in models that are almost exclusively based on gene expression. The gene expression features selected by the classic approach are difficult to interpret as they often represent poorly studied combinations of genes, activated by aberrations in upstream signaling pathways. RESULTS: To utilize all data types in a more balanced way, we developed TANDEM, a two-stage approach in which the first stage explains response using upstream features (mutations, copy number, methylation and cancer type) and the second stage explains the remainder using downstream features (gene expression). Applying TANDEM to 934 cell lines profiled across 265 drugs (GDSC1000), we show that the resulting models are more interpretable, while retaining the same predictive performance as the classic approach. Using the more balanced contributions per data type as determined with TANDEM, we find that response to MAPK pathway inhibitors is largely predicted by mutation data, while predicting response to DNA damaging agents requires gene expression data, in particular SLFN11 expression. AVAILABILITY AND IMPLEMENTATION: TANDEM is available as an R package on CRAN (for more information, see http://ccb.nki.nl/software/tandem). CONTACT: m.michaut@nki.nl or l.wessels@nki.nl SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

OncoScape: Exploring the cancer aberration landscape by genomic data fusion.

Although large-scale efforts for molecular profiling of cancer samples provide multiple data types for many samples, most approaches for finding candidate cancer genes rely on somatic mutations and DNA copy number only. We present a new method, OncoScape, which exploits five complementary data types across 11 cancer types to identify new candidate cancer genes. We find many rarely mutated genes that are strongly affected by other aberrations. We retrieve the majority of known cancer genes but also new candidates such as STK31 and MSRA with very high confidence. Several genes show a dual oncogene- and tumor suppressor-like behavior depending on the tumor type. Most notably, the well-known tumor suppressor RB1 shows strong oncogene-like signal in colon cancer. We applied OncoScape to cell lines representing ten cancer types, providing the most comprehensive comparison of aberrations in cell lines and tumor samples to date. This revealed that glioblastoma, breast and colon cancer show strong similarity between cell lines and tumors, while head and neck squamous cell carcinoma and bladder cancer, exhibit very little similarity between cell lines and tumors. To facilitate exploration of the cancer aberration landscape, we created a web portal enabling interactive analysis of OncoScape results (http://ccb.nki.nl/software/oncoscape).

Integration of genomic, transcriptomic and proteomic data identifies two biologically distinct subtypes of invasive lobular breast cancer.

Invasive lobular carcinoma (ILC) is the second most frequently occurring histological breast cancer subtype after invasive ductal carcinoma (IDC), accounting for around 10% of all breast cancers. The molecular processes that drive the development of ILC are still largely unknown. We have performed a comprehensive genomic, transcriptomic and proteomic analysis of a large ILC patient cohort and present here an integrated molecular portrait of ILC. Mutations in CDH1 and in the PI3K pathway are the most frequent molecular alterations in ILC. We identified two main subtypes of ILCs: (i) an immune related subtype with mRNA up-regulation of PD-L1, PD-1 and CTLA-4 and greater sensitivity to DNA-damaging agents in representative cell line models; (ii) a hormone related subtype, associated with Epithelial to Mesenchymal Transition (EMT), and gain of chromosomes 1q and 8q and loss of chromosome 11q. Using the somatic mutation rate and eIF4B protein level, we identified three groups with different clinical outcomes, including a group with extremely good prognosis. We provide a comprehensive overview of the molecular alterations driving ILC and have explored links with therapy response. This molecular characterization may help to tailor treatment of ILC through the application of specific targeted, chemo- and/or immune-therapies.

ERBB2 Mutations Characterize a Subgroup of Muscle-invasive Bladder Cancers with Excellent Response to Neoadjuvant Chemotherapy.

UNLABELLED: A pathologic complete response to neoadjuvant chemotherapy (NAC) containing platinum is a strong prognostic determinant for patients with muscle-invasive bladder cancer (MIBC). Despite comprehensive molecular characterization of bladder cancer, associations of molecular alterations with treatment response are still largely unknown. We selected pathologic complete responders (ypT0N0; n=38) and nonresponders (higher than ypT2; n=33) from a cohort of high-grade MIBC patients treated with NAC. DNA was isolated from prechemotherapy tumor tissue and used for next-generation sequencing of 178 cancer-associated genes (discovery cohort) or targeted sequencing (validation cohort). We found that 9 of 38 complete responders had erb-b2 receptor tyrosine kinase 2 (ERBB2) missense mutations, whereas none of 33 nonresponders had ERBB2 mutations (p=0.003). ERBB2 missense mutations in complete responders were mostly confirmed activating mutations. ERCC2 missense mutations, recently found associated with response to NAC, were more common in complete responders; however, this association did not reach statistical significance in our cohort. We conclude that ERBB2 missense mutations characterize a subgroup of MIBC patients with an excellent response to NAC. PATIENT SUMMARY: In this report we looked for genetic alterations that can predict the response to neoadjuvant chemotherapy (NAC) in bladder cancer. We found that mutations in the gene ERBB2 are exclusively present in patients responding to NAC.

Next generation sequencing of triple negative breast cancer to find predictors for chemotherapy response.

INTRODUCTION: In triple negative breast cancers (TNBC) the initial response to chemotherapy is often favorable, but relapse and chemotherapy resistance frequently occur in advanced disease. Hence there is an urgent need for targeted treatments in this breast cancer subtype. In the current study we deep sequenced DNA of tumors prior to chemotherapy to search for predictors of response or resistance. METHODS: Next generation sequencing (NGS) was performed for 1,977 genes involved in tumorigenesis. DNA from 56 pre-treatment TNBC-biopsies was sequenced, as well as matched normal DNA. Following their tumor biopsy, patients started neoadjuvant chemotherapy with doxorubicin and cyclophosphamide. We studied associations between genetic alterations and three clinical variables: chemotherapy response, relapse-free survival and BRCA proficiency. RESULTS: The mutations observed were diverse and few recurrent mutations were detected. Most mutations were in TP53, TTN, and PIK3CA (55 %, 14 %, and 9 %, respectively). The mutation rates were similar between responders and non-responders (average mutation rate 9 vs 8 mutations). No recurrent mutations were associated with chemotherapy response or relapse. Interestingly, PIK3CA mutations were exclusively observed in patients proficient for BRCA1. Samples with a relapse had a higher copy number alteration rate, and amplifications of TTK and TP53BP2 were associated with a poor chemotherapy response. CONCLUSIONS: In this homogenous cohort of TNBCs few recurrent mutations were found. However, PIK3CA mutations were associated with BRCA proficiency, which can have clinical consequences in the near future.

BRCA1-like signature in triple negative breast cancer: Molecular and clinical characterization reveals subgroups with therapeutic potential.

Triple negative (TN) breast cancers make up some 15% of all breast cancers. Approximately 10-15% are mutant for the tumor suppressor, BRCA1. BRCA1 is required for homologous recombination-mediated DNA repair and deficiency results in genomic instability. BRCA1-mutated tumors have a specific pattern of genomic copy number aberrations that can be used to classify tumors as BRCA1-like or non-BRCA1-like. BRCA1 mutation, promoter methylation, BRCA1-like status and genome-wide expression data was determined for 112 TN breast cancer samples with long-term follow-up. Mutation status for 21 known DNA repair genes and PIK3CA was assessed. Gene expression and mutation frequency in BRCA1-like and non-BRCA1-like tumors were compared. Multivariate survival analysis was performed using the Cox proportional hazards model. BRCA1 germline mutation was identified in 10% of patients and 15% of tumors were BRCA1 promoter methylated. Fifty-five percent of tumors classified as BRCA1-like. The functions of genes significantly up-regulated in BRCA1-like tumors included cell cycle and DNA recombination and repair. TP53 was found to be frequently mutated in BRCA1-like (P < 0.05), while PIK3CA was frequently mutated in non-BRCA1-like tumors (P < 0.05). A significant association with worse prognosis was evident for patients with BRCA1-like tumors (adjusted HR = 3.32, 95% CI = 1.30-8.48, P = 0.01). TN tumors can be further divided into two major subgroups, BRCA1-like and non-BRCA1-like with different mutation and expression patterns and prognoses. Based on these molecular patterns, subgroups may be more sensitive to specific targeted agents such as PI3K or PARP inhibitors.

Parallel in vivo and in vitro melanoma RNAi dropout screens reveal synthetic lethality between hypoxia and DNA damage response inhibition.

To identify factors preferentially necessary for driving tumor expansion, we performed parallel in vitro and in vivo negative-selection short hairpin RNA (shRNA) screens. Melanoma cells harboring shRNAs targeting several DNA damage response (DDR) kinases had a greater selective disadvantage in vivo than in vitro, indicating an essential contribution of these factors during tumor expansion. In growing tumors, DDR kinases were activated following hypoxia. Correspondingly, depletion or pharmacologic inhibition of DDR kinases was toxic to melanoma cells, including those that were resistant to BRAF inhibitor, and this could be enhanced by angiogenesis blockade. These results reveal that hypoxia sensitizes melanomas to targeted inhibition of the DDR and illustrate the utility of in vivo shRNA dropout screens for the identification of pharmacologically tractable targets.

Distinct types of disorder in the human proteome: functional implications for alternative splicing.

Intrinsically disordered regions have been associated with various cellular processes and are implicated in several human diseases, but their exact roles remain unclear. We previously defined two classes of conserved disordered regions in budding yeast, referred to as "flexible" and "constrained" conserved disorder. In flexible disorder, the property of disorder has been positionally conserved during evolution, whereas in constrained disorder, both the amino acid sequence and the property of disorder have been conserved. Here, we show that flexible and constrained disorder are widespread in the human proteome, and are particularly common in proteins with regulatory functions. Both classes of disordered sequences are highly enriched in regions of proteins that undergo tissue-specific (TS) alternative splicing (AS), but not in regions of proteins that undergo general (i.e., not tissue-regulated) AS. Flexible disorder is more highly enriched in TS alternative exons, whereas constrained disorder is more highly enriched in exons that flank TS alternative exons. These latter regions are also significantly more enriched in potential phosphosites and other short linear motifs associated with cell signaling. We further show that cancer driver mutations are significantly enriched in regions of proteins associated with TS and general AS. Collectively, our results point to distinct roles for TS alternative exons and flanking exons in the dynamic regulation of protein interaction networks in response to signaling activity, and they further suggest that alternatively spliced regions of proteins are often functionally altered by mutations responsible for cancer.

A robust genomic signature for the detection of colorectal cancer patients with microsatellite instability phenotype and high mutation frequency.

Microsatellite instability (MSI) occurs in 10-20% of colorectal tumours and is associated with good prognosis. Here we describe the development and validation of a genomic signature that identifies colorectal cancer patients with MSI caused by DNA mismatch repair deficiency with high accuracy. Microsatellite status for 276 stage II and III colorectal tumours has been determined. Full-genome expression data was used to identify genes that correlate with MSI status. A subset of these samples (n = 73) had sequencing data for 615 genes available. An MSI gene signature of 64 genes was developed and validated in two independent validation sets: the first consisting of frozen samples from 132 stage II patients; and the second consisting of FFPE samples from the PETACC-3 trial (n = 625). The 64-gene MSI signature identified MSI patients in the first validation set with a sensitivity of 90.3% and an overall accuracy of 84.8%, with an AUC of 0.942 (95% CI, 0.888-0.975). In the second validation, the signature also showed excellent performance, with a sensitivity 94.3% and an overall accuracy of 90.6%, with an AUC of 0.965 (95% CI, 0.943-0.988). Besides correct identification of MSI patients, the gene signature identified a group of MSI-like patients that were MSS by standard assessment but MSI by signature assessment. The MSI-signature could be linked to a deficient MMR phenotype, as both MSI and MSI-like patients showed a high mutation frequency (8.2% and 6.4% of 615 genes assayed, respectively) as compared to patients classified as MSS (1.6% mutation frequency). The MSI signature showed prognostic power in stage II patients (n = 215) with a hazard ratio of 0.252 (p = 0.0145). Patients with an MSI-like phenotype had also an improved survival when compared to MSS patients. The MSI signature was translated to a diagnostic microarray and technically and clinically validated in FFPE and frozen samples.

Mapping the Hsp90 genetic interaction network in Candida albicans reveals environmental contingency and rewired circuitry.

The molecular chaperone Hsp90 regulates the folding of diverse signal transducers in all eukaryotes, profoundly affecting cellular circuitry. In fungi, Hsp90 influences development, drug resistance, and evolution. Hsp90 interacts with -10% of the proteome in the model yeast Saccharomyces cerevisiae, while only two interactions have been identified in Candida albicans, the leading fungal pathogen of humans. Utilizing a chemical genomic approach, we mapped the C. albicans Hsp90 interaction network under diverse stress conditions. The chaperone network is environmentally contingent, and most of the 226 genetic interactors are important for growth only under specific conditions, suggesting that they operate downstream of Hsp90, as with the MAPK Hog1. Few interactors are important for growth in many environments, and these are poised to operate upstream of Hsp90, as with the protein kinase CK2 and the transcription factor Ahr1. We establish environmental contingency in the first chaperone network of a fungal pathogen, novel effectors upstream and downstream of Hsp90, and network rewiring over evolutionary time.