Rgtsp: a generalized top scoring pairs package for class prediction.

SUMMARY: A top scoring pair (TSP) classifier consists of a pair of variables whose relative ordering can be used for accurately predicting the class label of a sample. This classification rule has the advantage of being easily interpretable and more robust against technical variations in data, as those due to different microarray platforms. Here we describe a parallel implementation of this classifier which significantly reduces the training time, and a number of extensions, including a multi-class approach, which has the potential of improving the classification performance. AVAILABILITY AND IMPLEMENTATION: Full C++ source code and R package Rgtsp are freely available from http://lausanne.isb-sib.ch/~vpopovic/research/. The implementation relies on existing OpenMP libraries

Cross-species analysis of genetically engineered mouse models of MAPK-driven colorectal cancer identifies hallmarks of the human disease

Effective treatment options for advanced colorectal cancer (CRC) are limited, survival rates are poor and this disease continues to be a leading cause of cancer-related deaths worldwide. Despite being a highly heterogeneous disease, a large subset of individuals with sporadic CRC typically harbor relatively few established ‘driver’ lesions. Here, we describe a collection of genetically engineered mouse models (GEMMs) of sporadic CRC that combine lesions frequently altered in human patients, including well-characterized tumor suppressors and activators of MAPK signaling. Primary tumors from these models were profiled, and individual GEMM tumors segregated into groups based on their genotypes. Unique allelic and genotypic expression signatures were generated from these GEMMs and applied to clinically annotated human CRC patient samples. We provide evidence that a Kras signature derived from these GEMMs is capable of distinguishing human tumors harboring KRAS mutation, and tracks with poor prognosis in two independent human patient cohorts. Furthermore, the analysis of a panel of human CRC cell lines suggests that high expression of the GEMM Kras signature correlates with sensitivity to targeted pathway inhibitors. Together, these findings implicate GEMMs as powerful preclinical tools with the capacity to recapitulate relevant human disease biology, and support the use of genetic signatures generated in these models to facilitate future drug discovery and validation efforts

DNA Topoisomerase I Gene Copy Number and mRNA Expression Assessed as Predictive Biomarkers for Adjuvant Irinotecan in Stage II/III Colon Cancer.

PURPOSE: Prospective-retrospective assessment of theTOP1gene copy number andTOP1mRNA expression as predictive biomarkers for adjuvant irinotecan in stage II/III colon cancer. EXPERIMENTAL DESIGN: Formalin-fixed, paraffin-embedded tissue microarrays were obtained from an adjuvant colon cancer trial (PETACC3) where patients were randomized to 5-fluorouracil/folinic acid with or without additional irinotecan.TOP1copy number status was analyzed by fluorescencein situhybridization (FISH) using aTOP1/CEN20 dual-probe combination.TOP1mRNA data were available from previous analyses. RESULTS: TOP1FISH and follow-up data were obtained from 534 patients.TOP1gain was identified in 27% using a single-probe enumeration strategy (≥4TOP1signals per cell) and in 31% when defined by aTOP1/CEN20 ratio ≥ 1.5. The effect of additional irinotecan was not dependent onTOP1FISH status.TOP1mRNA data were available from 580 patients with stage III disease. Benefit of irinotecan was restricted to patients characterized byTOP1mRNA expression ≥ third quartile (RFS: HRadjusted, 0.59;P= 0.09; OS: HRadjusted, 0.44;P= 0.03). The treatment byTOP1mRNA interaction was not statistically significant, but in exploratory multivariable fractional polynomial interaction analysis, increasingTOP1mRNA values appeared to be associated with increasing benefit of irinotecan. CONCLUSIONS: In contrast to theTOP1copy number, a trend was demonstrated for a predictive property ofTOP1mRNA expression. On the basis ofTOP1mRNA, it might be possible to identify a subgroup of patients where an irinotecan doublet is a clinically relevant option in the adjuvant setting of colon cancer.Clin Cancer Res; 22(7); 1621-31. ©2015 AACR

miR-345 in metastatic colorectal cancer: a non-invasive biomarker for clinical outcome in non-KRAS mutant patients treated with 3rd line cetuximab and irinotecan.

INTRODUCTION: MicroRNAs (miRNAs) have important regulatory functions in cellular processes and have shown promising potential as prognostic markers for disease outcome in patients with cancer. The aim of the present study was to find miRNA expression profiles in whole blood that were prognostic for overall survival (OS) in patients with metastatic colorectal cancer (mCRC) treated with cetuximab and irinotecan. METHODS: From 138 patients with mCRC in 3rd line therapy with cetuximab and irinotecan in a prospective phase II study, 738 pretreatment miRNAs were isolated and profiled from whole blood using the TaqMan MicroRNA Array v2.0. Mutation status of KRAS, BRAF, and PI3KCA was known. RESULTS: After Bonferroni adjustment, 6 miRNAs: (miR-345, miR-143, miR-34a*, miR-628-5p, miR-886-3p and miR-324-3p), were found associated with short OS. miR-345 was the strongest prognostic miRNA, significant in the full cohort and in the non-KRAS mutant population. miR-345, as a continuous variable in the full cohort, resulted in a hazard ratio (HR) of 2.38 per IQR (CI 95%: 1.8-3.1, P-value = 2.86e-07, Bonferroni adjusted, univariable analysis) and a HR = 1.75 per IQR (CI 95%: 1.24-2.48, P-Wald = 1.45e-03) in the multivariable analysis adjusted for gender, age, KRAS, PI3KCA and performance status. miR-345 was prognostic in progression-free survival (PFS) with a HR = 1.63 per IQR (CI 95%: 1.25-2.114, P-Wald = 2.92e-4) in the multivariable analysis. In addition, high miR-345 expression was associated with lack of response to treatment with cetuximab and irinotecan. CONCLUSION: We identified miR-345 in whole blood as a potential biomarker for clinical outcome. MiR-345 was a single prognostic biomarker for both OS and PFS in all patients and also in the non-KRAS mutant population

Context-dependent interpretation of the prognostic value of BRAF and KRAS mutations in colorectal cancer.

The mutation status of the BRAF and KRAS genes has been proposed as prognostic biomarker in colorectal cancer. Of them, only the BRAF V600E mutation has been validated independently as prognostic for overall survival and survival after relapse, while the prognostic value of KRAS mutation is still unclear. We investigated the prognostic value of BRAF and KRAS mutations in various contexts defined by stratifications of the patient population. We retrospectively analyzed a cohort of patients with stage II and III colorectal cancer from the PETACC-3 clinical trial (N = 1,423), by assessing the prognostic value of the BRAF and KRAS mutations in subpopulations defined by all possible combinations of the following clinico-pathological variables: T stage, N stage, tumor site, tumor grade and microsatellite instability status. In each such subpopulation, the prognostic value was assessed by log rank test for three endpoints: overall survival, relapse-free survival, and survival after relapse. The significance level was set to 0.01 for Bonferroni-adjusted p-values, and a second threshold for a trend towards statistical significance was set at 0.05 for unadjusted p-values. The significance of the interactions was tested by Wald test, with significance level of 0.05. In stage II-III colorectal cancer, BRAF mutation was confirmed a marker of poor survival only in subpopulations involving microsatellite stable and left-sided tumors, with higher effects than in the whole population. There was no evidence for prognostic value in microsatellite instable or right-sided tumor groups. We found that BRAF was also prognostic for relapse-free survival in some subpopulations. We found no evidence that KRAS mutations had prognostic value, although a trend was observed in some stratifications. We also show evidence of heterogeneity in survival of patients with BRAF V600E mutation. The BRAF mutation represents an additional risk factor only in some subpopulations of colorectal cancers, in others having limited prognostic value. However, in the subpopulations where it is prognostic, it represents a marker of much higher risk than previously considered. KRAS mutation status does not seem to represent a strong prognostic variable

Gene expression patterns unveil a new level of molecular heterogeneity in colorectal cancer.

The recognition that colorectal cancer (CRC) is a heterogeneous disease in terms of clinical behaviour and response to therapy translates into an urgent need for robust molecular disease subclassifiers that can explain this heterogeneity beyond current parameters (MSI, KRAS, BRAF). Attempts to fill this gap are emerging. The Cancer Genome Atlas (TGCA) reported two main CRC groups, based on the incidence and spectrum of mutated genes, and another paper reported an EMT expression signature defined subgroup. We performed a prior free analysis of CRC heterogeneity on 1113 CRC gene expression profiles and confronted our findings to established molecular determinants and clinical, histopathological and survival data. Unsupervised clustering based on gene modules allowed us to distinguish at least five different gene expression CRC subtypes, which we call surface crypt-like, lower crypt-like, CIMP-H-like, mesenchymal and mixed. A gene set enrichment analysis combined with literature search of gene module members identified distinct biological motifs in different subtypes. The subtypes, which were not derived based on outcome, nonetheless showed differences in prognosis. Known gene copy number variations and mutations in key cancer-associated genes differed between subtypes, but the subtypes provided molecular information beyond that contained in these variables. Morphological features significantly differed between subtypes. The objective existence of the subtypes and their clinical and molecular characteristics were validated in an independent set of 720 CRC expression profiles. Our subtypes provide a novel perspective on the heterogeneity of CRC. The proposed subtypes should be further explored retrospectively on existing clinical trial datasets and, when sufficiently robust, be prospectively assessed for clinical relevance in terms of prognosis and treatment response predictive capacity. Original microarray data were uploaded to the ArrayExpress database (http://www.ebi.ac.uk/arrayexpress/) under Accession Nos E-MTAB-990 and E-MTAB-1026. © 2013 Swiss Institute of Bioinformatics. Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland

Modulation of Angiogenic and Inflammatory Response in Glioblastoma by Hypoxia

Glioblastoma are rapidly proliferating brain tumors in which hypoxia is readily recognizable, as indicated by focal or extensive necrosis and vascular proliferation, two independent diagnostic criteria for glioblastoma. Gene expression profiling of glioblastoma revealed a gene expression signature associated with hypoxia-regulated genes. The correlated gene set emerging from unsupervised analysis comprised known hypoxia-inducible genes involved in angiogenesis and inflammation such as VEGF and BIRC3, respectively. The relationship between hypoxia-modulated angiogenic genes and inflammatory genes was associated with outcome in our cohort of glioblastoma patients treated within prospective clinical trials of combined chemoradiotherapy. The hypoxia regulation of several new genes comprised in this cluster including ZNF395, TNFAIP3, and TREM1 was experimentally confirmed in glioma cell lines and primary monocytes exposed to hypoxia in vitro. Interestingly, the cluster seems to characterize differential response of tumor cells, stromal cells and the macrophage/microglia compartment to hypoxic conditions. Most genes classically associated with the inflammatory compartment are part of the NF-kappaB signaling pathway including TNFAIP3 and BIRC3 that have been shown to be involved in resistance to chemotherapy

Chromosome 7 gain and DNA hypermethylation at the HOXA10 locus are associated with expression of a stem cell related HOX-signature in glioblastoma.

BACKGROUND: HOX genes are a family of developmental genes that are expressed neither in the developing forebrain nor in the normal brain. Aberrant expression of a HOX-gene dominated stem-cell signature in glioblastoma has been linked with increased resistance to chemo-radiotherapy and sustained proliferation of glioma initiating cells. Here we describe the epigenetic and genetic alterations and their interactions associated with the expression of this signature in glioblastoma. RESULTS: We observe prominent hypermethylation of the HOXA locus 7p15.2 in glioblastoma in contrast to non-tumoral brain. Hypermethylation is associated with a gain of chromosome 7, a hallmark of glioblastoma, and may compensate for tumor-driven enhanced gene dosage as a rescue mechanism by preventing undue gene expression. We identify the CpG island of the HOXA10 alternative promoter that appears to escape hypermethylation in the HOX-high glioblastoma. An additive effect of gene copy gain at 7p15.2 and DNA methylation at key regulatory CpGs in HOXA10 is significantly associated with HOX-signature expression. Additionally, we show concordance between methylation status and presence of active or inactive chromatin marks in glioblastoma-derived spheres that are HOX-high or HOX-low, respectively. CONCLUSIONS: Based on these findings, we propose co-evolution and interaction between gene copy gain, associated with a gain of chromosome 7, and additional epigenetic alterations as key mechanisms triggering a coordinated, but inappropriate, HOX transcriptional program in glioblastoma

FOXC2 controls adult lymphatic endothelial specialization, function, and gut lymphatic barrier preventing multiorgan failure.

The mechanisms maintaining adult lymphatic vascular specialization throughout life and their role in coordinating inter-organ communication to sustain homeostasis remain elusive. We report that inactivation of the mechanosensitive transcription factor Foxc2 in adult lymphatic endothelium leads to a stepwise intestine-to-lung systemic failure. Foxc2 loss compromised the gut epithelial barrier, promoted dysbiosis and bacterial translocation to peripheral lymph nodes, and increased circulating levels of purine metabolites and angiopoietin-2. Commensal microbiota depletion dampened systemic pro-inflammatory cytokine levels, corrected intestinal lymphatic dysfunction, and improved survival. Foxc2 loss skewed the specialization of lymphatic endothelial subsets, leading to populations with mixed, pro-fibrotic identities and to emergence of lymph node-like endothelial cells. Our study uncovers a cross-talk between lymphatic vascular function and commensal microbiota, provides single-cell atlas of lymphatic endothelial subtypes, and reveals organ-specific and systemic effects of dysfunctional lymphatics. These effects potentially contribute to the pathogenesis of diseases, such as inflammatory bowel disease, cancer, or lymphedema

CC+: a relational database of coiled-coil structures

We introduce the CC+ Database, a detailed, searchable repository of coiled-coil assignments, which is freely available at http://coiledcoils.chm.bris.ac.uk/ccplus. Coiled coils were identified using the program SOCKET, which locates coiled coils based on knobs-into-holes packing of side chains between α-helices. A method for determining the overall sequence identity of coiled-coil sequences was introduced to reduce statistical bias inherent in coiled-coil data sets. There are two points of entry into the CC+ Database: the ‘Periodic Table of Coiled-coil Structures’, which presents a graphical path through coiled-coil space based on manually validated data, and the ‘Dynamic Interface’, which allows queries of the database at different levels of complexity and detail. The latter entry level, which is the focus of this article, enables the efficient and rapid compilation of subsets of coiled-coil structures. These can be created and interrogated with increasingly sophisticated pull-down, keyword and sequence-based searches to return detailed structural and sequence information. Also provided are means for outputting the retrieved coiled-coil data in various formats, including PyMOL and RasMol scripts, and Position-Specific Scoring Matrices (or amino-acid profiles), which may be used, for example, in protein-structure prediction