Particle swarm optimization artificial intelligence technique for gene signature discovery in transcriptomic cohorts.

The development of gene signatures is key for delivering personalized medicine, despite only a few signatures being available for use in the clinic for cancer patients. Gene signature discovery tends to revolve around identifying a single signature. However, it has been shown that various highly predictive signatures can be produced from the same dataset. This study assumes that the presentation of top ranked signatures will allow greater efforts in the selection of gene signatures for validation on external datasets and for their clinical translation. Particle swarm optimization (PSO) is an evolutionary algorithm often used as a search strategy and largely represented as binary PSO (BPSO) in this domain. BPSO, however, fails to produce succinct feature sets for complex optimization problems, thus affecting its overall runtime and optimization performance. Enhanced BPSO (EBPSO) was developed to overcome these shortcomings. Thus, this study will validate unique candidate gene signatures for different underlying biology from EBPSO on transcriptomics cohorts. EBPSO was consistently seen to be as accurate as BPSO with substantially smaller feature signatures and significantly faster runtimes. 100% accuracy was achieved in all but two of the selected data sets. Using clinical transcriptomics cohorts, EBPSO has demonstrated the ability to identify accurate, succinct, and significantly prognostic signatures that are unique from one another. This has been proposed as a promising alternative to overcome the issues regarding traditional single gene signature generation. Interpretation of key genes within the signatures provided biological insights into the associated functions that were well correlated to their cancer type.

Pathways, Processes, and Candidate Drugs Associated with a Hoxa Cluster-Dependency Model of Leukemia.

High expression of the HOXA cluster correlates with poor clinical outcome in acute myeloid leukemias, particularly those harboring rearrangements of the mixed-lineage-leukemia gene (MLLr). Whilst decreased HOXA expression acts as a readout for candidate experimental therapies, the necessity of the HOXA cluster for leukemia maintenance has not been fully explored. Primary leukemias were generated in hematopoietic stem/progenitor cells from Cre responsive transgenic mice for conditional deletion of the Hoxa locus. Hoxa deletion resulted in reduced proliferation and colony formation in which surviving leukemic cells retained at least one copy of the Hoxa cluster, indicating dependency. Comparative transcriptome analysis of Hoxa wild type and deleted leukemic cells identified a unique gene signature associated with key pathways including transcriptional mis-regulation in cancer, the Fanconi anemia pathway and cell cycle progression. Further bioinformatics analysis of the gene signature identified a number of candidate FDA-approved drugs for potential repurposing in high HOXA expressing cancers including MLLr leukemias. Together these findings support dependency for an MLLr leukemia on Hoxa expression and identified candidate drugs for further therapeutic evaluation.

NUQA: Estimating Cancer Spatial and Temporal Heterogeneity and Evolution through Alignment-Free Methods.

Longitudinal next-generation sequencing of cancer patient samples has enhanced our understanding of the evolution and progression of various cancers. As a result, and due to our increasing knowledge of heterogeneity, such sampling is becoming increasingly common in research and clinical trial sample collections. Traditionally, the evolutionary analysis of these cohorts involves the use of an aligner followed by subsequent stringent downstream analyses. However, this can lead to large levels of information loss due to the vast mutational landscape that characterizes tumor samples. Here, we propose an alignment-free approach for sequence comparison-a well-established approach in a range of biological applications including typical phylogenetic classification. Such methods could be used to compare information collated in raw sequence files to allow an unsupervised assessment of the evolutionary trajectory of patient genomic profiles. In order to highlight this utility in cancer research we have applied our alignment-free approach using a previously established metric, Jensen-Shannon divergence, and a metric novel to this area, Hellinger distance, to two longitudinal cancer patient cohorts in glioma and clear cell renal cell carcinoma using our software, NUQA. We hypothesize that this approach has the potential to reveal novel information about the heterogeneity and evolutionary trajectory of spatiotemporal tumor samples, potentially revealing early events in tumorigenesis and the origins of metastases and recurrences. Key words: alignment-free, Hellinger distance, exome-seq, evolution, phylogenetics, longitudinal.

ACE: A Workbench Using Evolutionary Genetic Algorithms for Analyzing Association in TCGA.

Modern methods of acquiring molecular data have improved rapidly in recent years, making it easier for researchers to collect large volumes of information. However, this has increased the challenge of recognizing interesting patterns within the data. Atlas Correlation Explorer (ACE) is a user-friendly workbench for seeking associations between attributes in The Cancer Genome Atlas (TCGA) database. It allows any combination of clinical and genomic data streams to be searched using an evolutionary algorithm approach. To showcase ACE, we assessed which RNA sequencing transcripts were associated with estrogen receptor (ESR1) in the TCGA breast cancer cohort. The analysis revealed already well-established associations with XBP1 and FOXA1, but also identified a strong association with CT62, a potential immunotherapeutic target with few previous associations with breast cancer. In conclusion, ACE can produce results for very large searches in a short time and will serve as an increasingly useful tool for biomarker discovery in the big data era. SIGNIFICANCE: ACE uses an evolutionary algorithm approach to perform large searches for associations between any combinations of data in the TCGA database.

Impact of Variable RNA-Sequencing Depth on Gene Expression Signatures and Target Compound Robustness: Case Study Examining Brain Tumor (Glioma) Disease Progression.

PURPOSE: Gene expression profiling can uncover biologic mechanisms underlying disease and is important in drug development. RNA sequencing (RNA-seq) is routinely used to assess gene expression, but costs remain high. Sample multiplexing reduces RNAseq costs; however, multiplexed samples have lower cDNA sequencing depth, which can hinder accurate differential gene expression detection. The impact of sequencing depth alteration on RNA-seq-based downstream analyses such as gene expression connectivity mapping is not known, where this method is used to identify potential therapeutic compounds for repurposing. METHODS: In this study, published RNA-seq profiles from patients with brain tumor (glioma) were assembled into two disease progression gene signature contrasts for astrocytoma. Available treatments for glioma have limited effectiveness, rendering this a disease of poor clinical outcome. Gene signatures were subsampled to simulate sequencing alterations and analyzed in connectivity mapping to investigate target compound robustness. RESULTS: Data loss to gene signatures led to the loss, gain, and consistent identification of significant connections. The most accurate gene signature contrast with consistent patient gene expression profiles was more resilient to data loss and identified robust target compounds. Target compounds lost included candidate compounds of potential clinical utility in glioma (eg, suramin, dasatinib). Lost connections may have been linked to low-abundance genes in the gene signature that closely characterized the disease phenotype. Consistently identified connections may have been related to highly expressed abundant genes that were ever-present in gene signatures, despite data reductions. Potential noise surrounding findings included false-positive connections that were gained as a result of gene signature modification with data loss. CONCLUSION: Findings highlight the necessity for gene signature accuracy for connectivity mapping, which should improve the clinical utility of future target compound discoveries.

Natural killer-like signature observed post therapy in locally advanced rectal cancer is a determinant of pathological response and improved survival.

Around 12-15% of patients with locally advanced rectal cancer undergo a pathologically complete response (tumor regression grade 4) to long-course preoperative chemoradiotherapy; the remainder exhibit a spectrum of tumor regression (tumor regression grade 1-3). Understanding therapy-related transcriptional alterations may enable better prediction of response as measured by progression-free and overall survival, in addition to aiding the development of improved strategies based on the underlying biology of the disease. To this end, we performed high-throughput gene expression profiling in 40 pairs of formalin-fixed paraffin-embedded rectal cancer biopsies and matched resections following long-course preoperative chemoradiotherapy (discovery cohort). Differential gene expression analysis was performed contrasting tumor regression grades in resections. Enumeration of the tumor microenvironment cell population was undertaken using in silico analysis of the transcriptional data, and real-time PCR validation of NCR1 undertaken. Immunohistochemistry and survival analysis was used to measure CD56+ cell populations in an independent cohort (n=150). Gene expression traits observed following long-course preoperative chemoradiotherapy in the discovery cohort suggested an increased abundance of natural killer cells in tumors that displayed a clinical response to CRT in a tumor regression grade-dependent manner. CD56+ natural killer-cell populations were measured by immunohistochemistry and found to be significantly higher in tumor regression grade 3 patients compared with tumor regression grade 1-2 in the validation cohort. Furthermore, it was observed that patients positive for CD56 cells after therapy had a better overall survival (HR=0.282, 95% CI=0.109-0.729, χ2=7.854, P=0.005). In conclusion, we have identified a novel post-therapeutic natural killer-like transcription signature in patients responding to long-course preoperative chemoradiotherapy. Furthermore, patients with a higher abundance of CD56-positive natural killer cells post long-course preoperative chemoradiotherapy had better overall survival. Therefore, harnessing a natural killer-like response after therapy may improve outcomes for locally advanced rectal cancer patients. Finally, we hypothesize that future assessment of this natural killer-like response in on-treatment biopsy material may inform clinical decision-making for treatment duration.

Cancer-cell intrinsic gene expression signatures overcome intratumoural heterogeneity bias in colorectal cancer patient classification.

Stromal-derived intratumoural heterogeneity (ITH) has been shown to undermine molecular stratification of patients into appropriate prognostic/predictive subgroups. Here, using several clinically relevant colorectal cancer (CRC) gene expression signatures, we assessed the susceptibility of these signatures to the confounding effects of ITH using gene expression microarray data obtained from multiple tumour regions of a cohort of 24 patients, including central tumour, the tumour invasive front and lymph node metastasis. Sample clustering alongside correlative assessment revealed variation in the ability of each signature to cluster samples according to patient-of-origin rather than region-of-origin within the multi-region dataset. Signatures focused on cancer-cell intrinsic gene expression were found to produce more clinically useful, patient-centred classifiers, as exemplified by the CRC intrinsic signature (CRIS), which robustly clustered samples by patient-of-origin rather than region-of-origin. These findings highlight the potential of cancer-cell intrinsic signatures to reliably stratify CRC patients by minimising the confounding effects of stromal-derived ITH.

Embracing an integromic approach to tissue biomarker research in cancer: Perspectives and lessons learned.

Modern approaches to biomedical research and diagnostics targeted towards precision medicine are generating 'big data' across a range of high-throughput experimental and analytical platforms. Integrative analysis of this rich clinical, pathological, molecular and imaging data represents one of the greatest bottlenecks in biomarker discovery research in cancer and other diseases. Following on from the publication of our successful framework for multimodal data amalgamation and integrative analysis, Pathology Integromics in Cancer (PICan), this article will explore the essential elements of assembling an integromics framework from a more detailed perspective. PICan, built around a relational database storing curated multimodal data, is the research tool sitting at the heart of our interdisciplinary efforts to streamline biomarker discovery and validation. While recognizing that every institution has a unique set of priorities and challenges, we will use our experiences with PICan as a case study and starting point, rationalizing the design choices we made within the context of our local infrastructure and specific needs, but also highlighting alternative approaches that may better suit other programmes of research and discovery. Along the way, we stress that integromics is not just a set of tools, but rather a cohesive paradigm for how modern bioinformatics can be enhanced. Successful implementation of an integromics framework is a collaborative team effort that is built with an eye to the future and greatly accelerates the processes of biomarker discovery, validation and translation into clinical practice.

KRAS mutant colorectal cancer gene signatures identified angiotensin II receptor blockers as potential therapies.

Colorectal cancer (CRC) is a life-threatening disease with high prevalence and mortality worldwide. The KRAS oncogene is mutated in approximately 40% of CRCs. While antibody based EGFR inhibitors (cetuximab and panitumumab) represent a major treatment strategy for advanced KRAS wild type (KRAS-WT) CRCs, there still remains no effective therapeutic course for advanced KRAS mutant (KRAS-MT) CRC patients.In this study, we employed a novel and comprehensive approach of gene expression connectivity mapping (GECM) to identify candidate compounds to target KRAS-MT tumors. We first created a combined KRAS-MT gene signature with 248 ranked significant genes using 677 CRC clinical samples. A series of 248 sub-signatures was then created containing an increasing number of the top ranked genes. As an input to GECM analysis, each sub-signature was translated into a statistically significant therapeutic drugs list, which was finally combined to obtain a single list of significant drugs.We identify four antihypertensive angiotensin II receptor blockers (ARBs) within the top 30 significant drugs indicating that these drugs have a mechanism of action that can alter the KRAS-MT CRC oncogenic signaling. A hypergeometric test (p-value = 6.57 × 10-6) confirmed that ARBs are significantly enriched in our results. These findings support the hypothesis that ARB antihypertensive drugs may directly block KRAS signaling resulting in improvement in patient outcome or, through a reversion to a KRAS wild-type phenotype, improve the response to anti-EGFR treatment. Antihypertensive angiotensin II receptor blockers are therefore worth further investigation as potential therapeutic candidates in this difficult category of advanced colorectal cancers.

Immune-Derived PD-L1 Gene Expression Defines a Subgroup of Stage II/III Colorectal Cancer Patients with Favorable Prognosis Who May Be Harmed by Adjuvant Chemotherapy.

A recent phase II study of patients with metastatic colorectal carcinoma showed that mismatch repair gene status was predictive of clinical response to PD-1-targeting immune checkpoint blockade. Further examination revealed strong correlation between PD-L1 protein expression and microsatellite instability (MSI) in stage IV colorectal carcinoma, suggesting that the amount of PD-L1 protein expression could identify late-stage patients who might benefit from immunotherapy. To assess whether the clinical associations between PD-L1 gene expression and MSI identified in metastatic colorectal carcinoma are also present in stage II/III colorectal carcinoma, we used in silico analysis to elucidate the cell types expressing the PD-L1 gene. We found a statistically significant association of PD-L1 gene expression with MSI in early-stage colorectal carcinoma (P < 0.001) and show that, unlike in non-colorectal carcinoma tumors, PD-L1 is derived predominantly from the immune infiltrate. We demonstrate that PD-L1 gene expression has positive prognostic value in the adjuvant disease setting (PD-L1(low) vs. PD-L1(high) HR = 9.09; CI, 2.11-39.10). PD-L1 gene expression had predictive value, as patients with high PD-L1 expression appear to be harmed by standard-of-care treatment (HR = 4.95; CI, 1.10-22.35). Building on the promising results from the metastatic colorectal carcinoma PD-1-targeting trial, we provide compelling evidence that patients with PD-L1(high)/MSI/immune(high) stage II/III colorectal carcinoma should not receive standard chemotherapy. This conclusion supports the rationale to clinically evaluate this patient subgroup for PD-1 blockade treatment. Cancer Immunol Res; 4(7); 582-91. ©2016 AACR.

QUADrATiC: scalable gene expression connectivity mapping for repurposing FDA-approved therapeutics.

BACKGROUND: Gene expression connectivity mapping has proven to be a powerful and flexible tool for research. Its application has been shown in a broad range of research topics, most commonly as a means of identifying potential small molecule compounds, which may be further investigated as candidates for repurposing to treat diseases. The public release of voluminous data from the Library of Integrated Cellular Signatures (LINCS) programme further enhanced the utilities and potentials of gene expression connectivity mapping in biomedicine. RESULTS: We describe QUADrATiC ( http://go.qub.ac.uk/QUADrATiC ), a user-friendly tool for the exploration of gene expression connectivity on the subset of the LINCS data set corresponding to FDA-approved small molecule compounds. It enables the identification of compounds for repurposing therapeutic potentials. The software is designed to cope with the increased volume of data over existing tools, by taking advantage of multicore computing architectures to provide a scalable solution, which may be installed and operated on a range of computers, from laptops to servers. This scalability is provided by the use of the modern concurrent programming paradigm provided by the Akka framework. The QUADrATiC Graphical User Interface (GUI) has been developed using advanced Javascript frameworks, providing novel visualization capabilities for further analysis of connections. There is also a web services interface, allowing integration with other programs or scripts. CONCLUSIONS: QUADrATiC has been shown to provide an improvement over existing connectivity map software, in terms of scope (based on the LINCS data set), applicability (using FDA-approved compounds), usability and speed. It offers potential to biological researchers to analyze transcriptional data and generate potential therapeutics for focussed study in the lab. QUADrATiC represents a step change in the process of investigating gene expression connectivity and provides more biologically-relevant results than previous alternative solutions.

Stratified analysis reveals chemokine-like factor (CKLF) as a potential prognostic marker in the MSI-immune consensus molecular subtype CMS1 of colorectal cancer.

The Colorectal Cancer (CRC) Subtyping Consortium (CRCSC) recently published four consensus molecular subtypes (CMS's) representing the underlying biology in CRC. The Microsatellite Instable (MSI) immune group, CMS1, has a favorable prognosis in early stage disease, but paradoxically has the worst prognosis following relapse, suggesting the presence of factors enabling neoplastic cells to circumvent this immune response. To identify the genes influencing subsequent poor prognosis in CMS1, we analyzed this subtype, centered on risk of relapse. In a cohort of early stage colon cancer (n=460), we examined, in silico, changes in gene expression within the CMS1 subtype and demonstrated for the first time the favorable prognostic value of chemokine-like factor (CKLF) gene expression in the adjuvant disease setting [HR=0.18, CI=0.04-0.89]. In addition, using transcription profiles originating from cell sorted CRC tumors, we delineated the source of CKLF transcription within the colorectal tumor microenvironment to the leukocyte component of these tumors. Further to this, we confirmed that CKLF gene expression is confined to distinct immune subsets in whole blood samples and primary cell lines, highlighting CKLF as a potential immune cell-derived factor promoting tumor immune-surveillance of nascent neoplastic cells, particularly in CMS1 tumors. Building on the recently reported CRCSC data, we provide compelling evidence that leukocyte-infiltrate derived CKLF expression is a candidate biomarker of favorable prognosis, specifically in MSI-immune stage II/III disease.

Challenging the Cancer Molecular Stratification Dogma: Intratumoral Heterogeneity Undermines Consensus Molecular Subtypes and Potential Diagnostic Value in Colorectal Cancer.

PURPOSE: A number of independent gene expression profiling studies have identified transcriptional subtypes in colorectal cancer with potential diagnostic utility, culminating in publication of a colorectal cancer Consensus Molecular Subtype classification. The worst prognostic subtype has been defined by genes associated with stem-like biology. Recently, it has been shown that the majority of genes associated with this poor prognostic group are stromal derived. We investigated the potential for tumor misclassification into multiple diagnostic subgroups based on tumoral region sampled. EXPERIMENTAL DESIGN: We performed multiregion tissue RNA extraction/transcriptomic analysis using colorectal-specific arrays on invasive front, central tumor, and lymph node regions selected from tissue samples from 25 colorectal cancer patients. RESULTS: We identified a consensus 30-gene list, which represents the intratumoral heterogeneity within a cohort of primary colorectal cancer tumors. Using a series of online datasets, we showed that this gene list displays prognostic potential HR = 2.914 (confidence interval 0.9286-9.162) in stage II/III colorectal cancer patients, but in addition, we demonstrated that these genes are stromal derived, challenging the assumption that poor prognosis tumors with stem-like biology have undergone a widespread epithelial-mesenchymal transition. Most importantly, we showed that patients can be simultaneously classified into multiple diagnostically relevant subgroups based purely on the tumoral region analyzed. CONCLUSIONS: Gene expression profiles derived from the nonmalignant stromal region can influence assignment of colorectal cancer transcriptional subtypes, questioning the current molecular classification dogma and highlighting the need to consider pathology sampling region and degree of stromal infiltration when employing transcription-based classifiers to underpin clinical decision making in colorectal cancer. Clin Cancer Res; 22(16); 4095-104. ©2016 AACRSee related commentary by Morris and Kopetz, p. 3989.

The Ste20-like kinase SLK is required for cell cycle progression through G2.

We have previously shown that the Ste20-like kinase SLK is a microtubule-associated protein that can regulate actin reorganization during cell adhesion and spreading (Wagner, S., Flood, T. A., O'Reilly, P., Hume, K., and Sabourin, L. A. (2002) J. Biol. Chem. 277, 37685-37692). Because of its association with the microtubule network, we investigated whether SLK plays a role in cell cycle progression, a process that requires microtubule dynamics during mitosis. Consistent with microtubule association in exponentially growing cells, our results showed that SLK co-localizes with the mitotic spindle in cells undergoing mitosis. Expression of a kinase-inactive mutant or SLK small interfering RNAs inhibited cell proliferation and resulted in an accumulation of quiescent cells stimulated to re-enter the cell cycle in the G2 phase. Cultures expressing the mutant SLK displayed a normal pattern of cyclin D, E, and B expression but failed to down-regulate cyclin A levels, suggesting that they cannot proceed through M phase. In addition, these cultures displayed low levels of both phospho-H3 and active p34/cdc2 kinase. Overexpression of active SLK resulted in ectopic spindle assembly and the induction of cell cycle re-entry of Xenopus oocytes, suggesting that SLK is required for progression through G2 upstream of H1 kinase activation.