cudaMap: a GPU accelerated program for gene expression connectivity mapping

BACKGROUND: Modern cancer research often involves large datasets and the use of sophisticated statistical techniques. Together these add a heavy computational load to the analysis, which is often coupled with issues surrounding data accessibility. Connectivity mapping is an advanced bioinformatic and computational technique dedicated to therapeutics discovery and drug re-purposing around differential gene expression analysis. On a normal desktop PC, it is common for the connectivity mapping task with a single gene signature to take > 2h to complete using sscMap, a popular Java application that runs on standard CPUs (Central Processing Units). Here, we describe new software, cudaMap, which has been implemented using CUDA C/C++ to harness the computational power of NVIDIA GPUs (Graphics Processing Units) to greatly reduce processing times for connectivity mapping. RESULTS: cudaMap can identify candidate therapeutics from the same signature in just over thirty seconds when using an NVIDIA Tesla C2050 GPU. Results from the analysis of multiple gene signatures, which would previously have taken several days, can now be obtained in as little as 10 minutes, greatly facilitating candidate therapeutics discovery with high throughput. We are able to demonstrate dramatic speed differentials between GPU assisted performance and CPU executions as the computational load increases for high accuracy evaluation of statistical significance. CONCLUSION: Emerging ‘omics’ technologies are constantly increasing the volume of data and information to be processed in all areas of biomedical research. Embracing the multicore functionality of GPUs represents a major avenue of local accelerated computing. cudaMap will make a strong contribution in the discovery of candidate therapeutics by enabling speedy execution of heavy duty connectivity mapping tasks, which are increasingly required in modern cancer research. cudaMap is open source and can be freely downloaded from http://purl.oclc.org/NET/cudaMap

Modelling the comet assay.

Abstract The single-cell gel electrophoresis technique or comet assay is widely regarded as a quick and reliable method of analysing DNA damage in individual cells. It has a proven track record from the fields of biomonitoring to nutritional studies. The assay operates by subjecting cells that are fixed in agarose to high salt and detergent lysis, thus removing all the cellular content except the DNA. By relaxing the DNA in an alkaline buffer, strands containing breaks are released from supercoiling. Upon electrophoresis, these strands are pulled out into the agarose, forming a tail which, when stained with a fluorescent dye, can be analysed by fluorescence microscopy. The intensity of this tail reflects the amount of DNA damage sustained. Despite being such an established and widely used assay, there are still many aspects of the comet assay which are not fully understood. The present review looks at how the comet assay is being used, and highlights some of its limitations. The protocol itself varies among laboratories, so results from similar studies may vary. Given such discrepancies, it would be attractive to break the assay into components to generate a mathematical model to investigate specific parameters. The technique The single-cell gel electrophoresis assay has been used for over 20 years to assess DNA damage. The assay is able to measure both single-and double-strand breaks with versatility and sensitivity The comet technique has been adapted over the years into different protocols such as the alkaline and neutral versions which were thought to discriminate between double-and single-strand breaks. This theory was revealed by Collins et al. [1] at the most recent International Comet Assay conference to be a misconception as alkaline conditions are not specific to detecting single-strand breaks. FISH (fluorescent in situ hybridization) techniques were also adapted for application to the comet where the DNA probe can be hybridized to a specific gene sequence or even to a whole chromosom

Immune-derived PD-L1 gene expression defines a subgroup of stage II/III colorectal cancer patients with favorable prognosis that may be harmed by adjuvant chemotherapy

Abstract 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 &amp;lt; 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-L1low vs. PD-L1high 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-L1high/MSI/immunehigh 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.</jats:p

Depletion of DNMT1 in differentiated human cells highlights key classes of sensitive genes and an interplay with polycomb repression

Additional file 3: Figure S2. Changes in methylation levels by genomic element. (A) Protein levels in knockdown lines by western blotting. As a control HCT116 colon cancer cells which are WT or have a homozygous mutation in DNMT1 (KO) are shown: the DNMT1-specific top band is indicated by the arrowhead at right. (B) Median levels of methylation are shown for each genomic element (listed at top). The positions of medians are also indicated at right (arrowheads). The differences between WT and KD medians were used to plot Fig. 1d. (C) Density distribution of methylation at the three main elements involved in gene regulation, shown by cell line. Demethylation seems most marked at gene bodies (Genes), indicated by increased density of probes at low methylation (β) values

Repurposing FDA approved drugs as radiosensitizers for treating hypoxic prostate cancer

Abstract Background The presence of hypoxia is a poor prognostic factor in prostate cancer and the hypoxic tumor microenvironment promotes radioresistance. There is potential for drug radiotherapy combinations to improve the therapeutic ratio. We aimed to investigate whether hypoxia-associated genes could be used to identify FDA approved drugs for repurposing for the treatment of hypoxic prostate cancer. Methods Hypoxia associated genes were identified and used in the connectivity mapping software QUADrATIC to identify FDA approved drugs as candidates for repurposing. Drugs identified were tested in vitro in prostate cancer cell lines (DU145, PC3, LNCAP). Cytotoxicity was investigated using the sulforhodamine B assay and radiosensitization using a clonogenic assay in normoxia and hypoxia. Results Menadione and gemcitabine had similar cytotoxicity in normoxia and hypoxia in all three cell lines. In DU145 cells, the radiation sensitizer enhancement ratio (SER) of menadione was 1.02 in normoxia and 1.15 in hypoxia. The SER of gemcitabine was 1.27 in normoxia and 1.09 in hypoxia. No radiosensitization was seen in PC3 cells. Conclusion Connectivity mapping can identify FDA approved drugs for potential repurposing that are linked to a radiobiologically relevant phenotype. Gemcitabine and menadione could be further investigated as potential radiosensitizers in prostate cancer

Transcriptional upregulation of c-MET is associated with invasion and tumor budding in colorectal cancer

c-MET and its ligand HGF are frequently overexpressed in colorectal cancer (CRC) and increased c-MET levels are found in CRC liver metastases. This study investigated the role of the HGF/c-MET axis in regulating migration/invasion in CRC, using pre-clinical models and clinical samples. Pre-clinically, we found marked upregulation of c-MET at both protein and mRNA levels in several invasive CRC cells. Down-regulation of c-MET using RNAi suppressed migration/invasion of parental and invasive CRC cells. Stimulation of CRC cells with rh-HGF or co-culture with HGF-expressing colonic myofibroblasts, resulted in significant increases in their migratory/invasive capacity. Importantly, HGF-induced c-MET activation promoted rapid downregulation of c-MET protein levels, while the MET transcript remained unaltered. Using RNA in situ hybridization (RNA ISH), we further showed that MET mRNA, but not protein levels, were significantly upregulated in tumor budding foci at the invasive front of a cohort of stage III CRC tumors (p < 0.001). Taken together, we show for the first time that transcriptional upregulation of MET is a key molecular event associated with CRC invasion and tumor budding. This data also indicates that RNA ISH, but not immunohistochemistry, provides a robust methodology to assess MET levels as a potential driving force of CRC tumor invasion and metastasis

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

Development and Validation of a 28-gene Hypoxia-related Prognostic Signature for Localized Prostate Cancer.

BACKGROUND: Hypoxia is associated with a poor prognosis in prostate cancer. This work aimed to derive and validate a hypoxia-related mRNA signature for localized prostate cancer. METHOD: Hypoxia genes were identified in vitro via RNA-sequencing and combined with in vivo gene co-expression analysis to generate a signature. The signature was independently validated in eleven prostate cancer cohorts and a bladder cancer phase III randomized trial of radiotherapy alone or with carbogen and nicotinamide (CON). RESULTS: A 28-gene signature was derived. Patients with high signature scores had poorer biochemical recurrence free survivals in six of eight independent cohorts of prostatectomy-treated patients (Log rank test P \u3c .05), with borderline significances achieved in the other two (P \u3c .1). The signature also predicted biochemical recurrence in patients receiving post-prostatectomy radiotherapy (n = 130, P = .007) or definitive radiotherapy alone (n = 248, P = .035). Lastly, the signature predicted metastasis events in a pooled cohort (n = 631, P = .002). Prognostic significance remained after adjusting for clinic-pathological factors and commercially available prognostic signatures. The signature predicted benefit from hypoxia-modifying therapy in bladder cancer patients (intervention-by-signature interaction test P = .0026), where carbogen and nicotinamide was associated with improved survival only in hypoxic tumours. CONCLUSION: A 28-gene hypoxia signature has strong and independent prognostic value for prostate cancer patients

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

This work was supported by a very generous grant from the Sean Crummey Memorial Fund. The staff and infrastructure provided by the N. Ireland Biobank and the Belfast Experimental Cancer Medicine Centre allowed this research to take place. These are supported by the Research and Development Division of the N. Ireland Public Health Agency and Cancer Research UK. We would also like to express our gratitude to the staff at the Grampian Biorepository for providing the Tissue microarray materials for validation purposes.Peer reviewedPostprin