sscMap: An extensible Java application for connecting small-molecule drugs using gene-expression signatures

Background
Connectivity mapping is a process to recognize novel pharmacological and toxicological properties in small molecules by comparing their gene expression signatures with others in a database. A simple and robust method for connectivity mapping with increased specificity and sensitivity was recently developed, and its utility demonstrated using experimentally derived gene signatures.

Results
This paper introduces sscMap (statistically significant connections' map), a Java application designed to undertake connectivity mapping tasks using the recently published method. The software is bundled with a default collection of reference gene-expression profiles based on the publicly available dataset from the Broad Institute Connectivity Map 02, which includes data from over 7000 Affymetrix microarrays, for over 1000 small-molecule compounds, and 6100 treatment instances in 5 human cell lines. In addition, the application allows users to add their custom collections of reference profiles and is applicable to a wide range of other 'omics technologies.

Conclusion
The utility of sscMap is two fold. First, it serves to make statistically significant connections between a user-supplied gene signature and the 6100 core reference profiles based on the Broad Institute expanded dataset. Second, it allows users to apply the same improved method to custom-built reference profiles which can be added to the database for future referencing. The software can be freely downloaded from http://purl.oclc.org/NET/sscMap

Application of connectivity mapping in predictive toxicology based on gene-expression similarity.

Connectivity mapping is the process of establishing connections between different biological states using gene-expression profiles or signatures. There are a number of applications but in toxicology the most pertinent is for understanding mechanisms of toxicity. In its essence the process involves comparing a query gene signature generated as a result of exposure of a biological system to a chemical to those in a database that have been previously derived. In the ideal situation the query gene-expression signature is characteristic of the event and will be matched to similar events in the database. Key criteria are therefore the means of choosing the signature to be matched and the means by which the match is made. In this article we explore these concepts with examples applicable to toxicology.

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.

The molecular signature of the stroma response in prostate cancer-induced osteoblastic bone metastasis highlights expansion of hematopoietic and prostate epithelial stem cell niches

The reciprocal interaction between cancer cells and the tissue-specific stroma is critical for primary and metastatic tumor growth progression. Prostate cancer cells colonize preferentially bone (osteotropism), where they alter the physiological balance between osteoblast-mediated bone formation and osteoclast-mediated bone resorption, and elicit prevalently an osteoblastic response (osteoinduction). The molecular cues provided by osteoblasts for the survival and growth of bone metastatic prostate cancer cells are largely unknown. We exploited the sufficient divergence between human and mouse RNA sequences together with redefinition of highly species-specific gene arrays by computer-aided and experimental exclusion of cross-hybridizing oligonucleotide probes. This strategy allowed the dissection of the stroma (mouse) from the cancer cell (human) transcriptome in bone metastasis xenograft models of human osteoinductive prostate cancer cells (VCaP and C4-2B). As a result, we generated the osteoblastic bone metastasis-associated stroma transcriptome (OB-BMST). Subtraction of genes shared by inflammation, wound healing and desmoplastic responses, and by the tissue type-independent stroma responses to a variety of non-osteotropic and osteotropic primary cancers generated a curated gene signature ("Core" OB-BMST) putatively representing the bone marrow/bone-specific stroma response to prostate cancer-induced, osteoblastic bone metastasis. The expression pattern of three representative Core OB-BMST genes (PTN, EPHA3 and FSCN1) seems to confirm the bone specificity of this response. A robust induction of genes involved in osteogenesis and angiogenesis dominates both the OB-BMST and Core OB-BMST. This translates in an amplification of hematopoietic and, remarkably, prostate epithelial stem cell niche components that may function as a self-reinforcing bone metastatic niche providing a growth support specific for osteoinductive prostate cancer cells. The induction of this combinatorial stem cell niche is a novel mechanism that may also explain cancer cell osteotropism and local interference with hematopoiesis (myelophthisis). Accordingly, these stem cell niche components may represent innovative therapeutic targets and/or serum biomarkers in osteoblastic bone metastasis.

Natural Killer Cell-like signature observed in locally advanced rectal cancer after neoadjuvant chemoradiotherapy in a tumour regression grade dependent manner

Background: Around 10-15% of patients with locally advanced rectal cancer (LARC) undergo a pathologically complete response (TRG4) to neoadjuvant chemoradiotherapy; the rest of patients exhibit a spectrum of tumour regression (TRG1-3). Understanding therapy-related genomic alterations may help us to identify underlying biology or novel targets associated with response that could increase the efficacy of therapy in patients that do not benefit from the current standard of care.
Methods: 48 FFPE rectal cancer biopsies and matched resections were analysed using the WG-DASL HumanHT-12_v4 Beadchip array on the illumina iScan. Bioinformatic analysis was conducted in Partek genomics suite and R studio. Limma and glmnet packages were used to identify genes differentially expressed between tumour regression grades. Validation of microarray results will be carried out using IHC, RNAscope and RT-PCR.
Results: Immune response genes were observed from supervised analysis of the biopsies which may have predictive value. Differential gene expression from the resections as well as pre and post therapy analysis revealed induction of genes in a tumour regression dependent manner. Pathway mapping and Gene Ontology analysis of these genes suggested antigen processing and natural killer mediated cytotoxicity respectively. The natural killer-like gene signature was switched off in non-responders and on in the responders. IHC has confirmed the presence of Natural killer cells through CD56+ staining.
Conclusion: Identification of NK cell genes and CD56+ cells in patients responding to neoadjuvant chemoradiotherapy warrants further investigation into their association with tumour regression grade in LARC. NK cells are known to lyse malignant cells and determining whether their presence is a cause or consequence of response is crucial. Interrogation of the cytokines upregulated in our NK-like signature will help guide future in vitro models.

Development of a multi-gene PCR assay for the prediction of the response to hormone therapy in breast cancer

Deux tiers des cancers du sein expriment des récepteurs hormonaux ostrogéniques (tumeur ER-positive) et la croissance de ces tumeurs est stimulée par l’estrogène. Des traitements adjuvant avec des anti-estrogènes, tel que le Tamoxifen et les Inhibiteurs de l’Aromatase peuvent améliorer la survie des patientes atteinte de cancer du sein. Toutefois la thérapie hormonale n’est pas efficace dans toutes les tumeurs mammaires ER-positives. Les tumeurs peuvent présenter avec une résistance intrinsèque ou acquise au Tamoxifen. Présentement, c’est impossible de prédire quelle patiente va bénéficier ou non du Tamoxifen. Des études préliminaires du laboratoire de Dr. Mader, ont identifié le niveau d’expression de 20 gènes, qui peuvent prédire la réponse thérapeutique au Tamoxifen (survie sans récidive). Ces marqueurs, identifié en utilisant une analyse bioinformatique de bases de données publiques de profils d’expression des gènes, sont capables de discriminer quelles patientes vont mieux répondre au Tamoxifen. Le but principal de cette étude est de développer un outil de PCR qui peut évaluer le niveau d’expression de ces 20 gènes prédictif et de tester cette signature de 20 gènes dans une étude rétrospective, en utilisant des tumeurs de cancer du sein en bloc de paraffine, de patients avec une histoire médicale connue. Cet outil aurait donc un impact direct dans la pratique clinique. Des traitements futiles pourraient être éviter et l’indentification de tumeurs ER+ avec peu de chance de répondre à un traitement anti-estrogène amélioré. En conséquence, de la recherche plus appropriée pour les tumeurs résistantes au Tamoxifen, pourront se faire.

Evaluation of the psychrotrophic specific signatures for cspA gene and 16S rDNA on the phenotype of Bacillus cereus sensu strictu

This study characterised the psychrotrophic genotypes and phenotypes behaviour of 63 strains of Bacillus cereus sensu stricto isolated from dairy products. The presence of the cspA gene signature and the 16S rDNA mesophilic and/or psychrotrophic specific signatures was evaluated. Among the strains, 25 (39.7%) had the cspA gene signature, 38 (60.3%) had both mesophilic and psychrotrophic 16S rDNA signatures, 24 (38.1%) had only mesophilic and one exhibited only psychrotrophic. No strain grew at 7 °C. The results indicate that the presence of psychrotrophic signatures for cspA gene or the 16S rDNA did not ensure a psychrotrophic behaviour on a B. cereus phenotype. © 2013 Society of Dairy Technology.

Gene expression profiling in leiomyosarcomas and undifferentiated pleomorphic sarcomas: SRC as a new diagnostic marker

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Anti-HIV state but not apoptosis depends on IFN signature in CD4+ T cells

To gain insights into the molecular mechanisms underlying early host responses to HIV in the CD4(+) T cell target population, we examined gene expression in CD4(+) T cells isolated 24 h after ex vivo HIV infection of lymphocyte aggregate cultures derived from human tonsils. Gene profiling showed a distinct up-regulation of genes related to immune response and response to virus, notably of IFN-stimulated genes (ISGs), irrespective of the coreceptor tropism of the virus. This mostly IFN-alpha-dependent gene signature suggested the involvement of plasmacytoid dendritic cells, a principal component of the antiviral immune response. Indeed, depletion of plasmacytoid dendritic cells before HIV inoculation abrogated transcriptional up-regulation of several ISGs and resulted in increased levels of HIV replication. Treatment with a blocking anti-IFN-alphaR Ab yielded increased HIV replication; conversely, HIV replication was decreased in pDC-depleted cultures treated with IFN-alpha. Among up-regulated ISGs was also TRAIL, indicating a potential role of the IFN signature in apoptosis. However, a blocking anti-TRAIL Ab did not abrogate apoptosis of CD4(+) T cells in CXCR4-tropic HIV-infected cultures, suggesting the involvement of pathways other than TRAIL mediated. We conclude that acute HIV infection of lymphoid tissue results in up-regulation of ISGs in CD4(+) T cells, which induces an anti-HIV state but not apoptosis.

The molecular signature of the stroma response in prostate cancer-induced osteoblastic bone metastasis highlights expansion of hematopoietic and prostate epithelial stem cell niches

The reciprocal interaction between cancer cells and the tissue-specific stroma is critical for primary and metastatic tumor growth progression. Prostate cancer cells colonize preferentially bone (osteotropism), where they alter the physiological balance between osteoblast-mediated bone formation and osteoclast-mediated bone resorption, and elicit prevalently an osteoblastic response (osteoinduction). The molecular cues provided by osteoblasts for the survival and growth of bone metastatic prostate cancer cells are largely unknown. We exploited the sufficient divergence between human and mouse RNA sequences together with redefinition of highly species-specific gene arrays by computer-aided and experimental exclusion of cross-hybridizing oligonucleotide probes. This strategy allowed the dissection of the stroma (mouse) from the cancer cell (human) transcriptome in bone metastasis xenograft models of human osteoinductive prostate cancer cells (VCaP and C4-2B). As a result, we generated the osteoblastic bone metastasis-associated stroma transcriptome (OB-BMST). Subtraction of genes shared by inflammation, wound healing and desmoplastic responses, and by the tissue type-independent stroma responses to a variety of non-osteotropic and osteotropic primary cancers generated a curated gene signature ("Core" OB-BMST) putatively representing the bone marrow/bone-specific stroma response to prostate cancer-induced, osteoblastic bone metastasis. The expression pattern of three representative Core OB-BMST genes (PTN, EPHA3 and FSCN1) seems to confirm the bone specificity of this response. A robust induction of genes involved in osteogenesis and angiogenesis dominates both the OB-BMST and Core OB-BMST. This translates in an amplification of hematopoietic and, remarkably, prostate epithelial stem cell niche components that may function as a self-reinforcing bone metastatic niche providing a growth support specific for osteoinductive prostate cancer cells. The induction of this combinatorial stem cell niche is a novel mechanism that may also explain cancer cell osteotropism and local interference with hematopoiesis (myelophthisis). Accordingly, these stem cell niche components may represent innovative therapeutic targets and/or serum biomarkers in osteoblastic bone metastasis.

Homozygous deletion mapping in myeloma samples identifies genes and an expression signature relevant to pathogenesis and outcome.

PURPOSE: Myeloma is a clonal malignancy of plasma cells. Poor-prognosis risk is currently identified by clinical and cytogenetic features. However, these indicators do not capture all prognostic information. Gene expression analysis can be used to identify poor-prognosis patients and this can be improved by combination with information about DNA-level changes. EXPERIMENTAL DESIGN: Using single nucleotide polymorphism-based gene mapping in combination with global gene expression analysis, we have identified homozygous deletions in genes and networks that are relevant to myeloma pathogenesis and outcome. RESULTS: We identified 170 genes with homozygous deletions and corresponding loss of expression. Deletion within the "cell death" network was overrepresented and cases with these deletions had impaired overall survival. From further analysis of these events, we have generated an expression-based signature associated with shorter survival in 258 patients and confirmed this signature in data from two independent groups totaling 800 patients. We defined a gene expression signature of 97 cell death genes that reflects prognosis and confirmed this in two independent data sets. CONCLUSIONS: We developed a simple 6-gene expression signature from the 97-gene signature that can be used to identify poor-prognosis myeloma in the clinical environment. This signature could form the basis of future trials aimed at improving the outcome of poor-prognosis myeloma.

In vitro functional characterisation of IGF-I : VN-induced breast cancer progression

Members of the insulin-like growth factor (IGF) family have been shown to play critical roles in normal growth and development, as well as in tumour biology. The IGF system is complex and the biological effects of the IGFs are determined by their diverse interactions between many molecules, including their interactions with extracellular matrix (ECM) proteins. Recent studies have demonstrated that IGFs associate with the ECM protein vitronectin (VN) through IGF-binding proteins (IGFBP) and that this interaction modulates IGF-stimulated biological functions, namely cell migration and cell survival through the cooperative involvement of the type-I IGF receptor (IGF-1R) and VN-binding integrins. Since IGFs play important roles in the transformation and progression of breast cancer and VN has been found to be over-expressed at the leading edge of breast tumours, this project aimed to describe the effects of IGF-I:VN interactions on breast cell function. This was undertaken to dissect the molecular mechanisms underlying IGF-I:VN-induced responses and to design inhibitors to block the effects of such interactions. The studies described herein demonstrate that the increase in migration of MCF-7 breast cancer cells in response to the IGF-I:IGFBP-5:VN complex is accompanied by differential expression of genes known to be involved in migration, invasion and/or survival, including Tissue-factor (TF), Stratifin (SFN), Ephrin-B2, Sharp-2 and PAI-1. This „migration gene signature‟ was confirmed using real-time PCR analysis. Substitution of the native IGF-I within the IGF-I:IGFBP:VN complex with the IGF-I analogue, \[L24]\[A31]-IGF-I, which has a reduced affinity for the IGF-1R, failed to stimulate cell migration and interestingly, also failed to induce the differential gene expression. This supports the involvement of the IGF-1R in mediating these changes in gene expression. Furthermore, lentiviral shRNA-mediated stable knockdown of TF and SFN completely abrogated the increased cell migration induced by IGF-I:IGFBP:VN complexes in MCF-7 cells. Indeed, when these cells were grown in 3D Matrigel™ cultures a decrease in the overall size of the 3D spheroids in response to the IGF-I:IGFBP:VN complexes was observed compared to the parental MCF-7 cells. This suggests that TF and SFN have a role in complex-stimulated cell survival. Moreover, signalling studies performed on cells with the reduced expression of either TF or SFN had a decreased IGF-1R activation, suggesting the involvement of signalling pathways downstream of IGF-1R in TF- and/or SFN-mediated cell migration and cell survival. Taken together, these studies provide evidence for a common mechanism activated downstream of the IGF-1R that induces the expression of the „migration gene signature‟ in response to the IGF-I:IGFBP:VN complex that confers breast cancer cells the propensity to migrate and survive. Given the functional significance of the interdependence of ECM and growth factor (GF) interactions in stimulating processes key to breast cancer progression, this project aimed at developing strategies to prevent such growth factor:ECM interactions in an effort to inhibit the downstream functional effects. This may result in the reduction in the levels of ECM-bound IGF-I present in close proximity to the cells, thereby leading to a reduction in the stimulation of IGF-1R present on the cell surface. Indeed, the inhibition of IGF-I-mediated effects through the disruption of its association with ECM would not alter the physiological levels of IGF-I and potentially only exert effects in situations where abnormal over expression of ECM proteins are found; namely carcinomas and hyperproliferative diseases. In summary, this PhD project has identified novel, innovative and realistic strategies that can be used in vitro to inhibit the functions exerted by the IGF-I:IGFBP:VN multiprotein complexes critical for cancer progression, with a potential to be translated into in vivo investigations. Furthermore, TF and SFN were found to mediate IGF-I:IGFBP:VN-induced effects, thereby revealing their potential to be used as therapeutic targets or as predictive biomarkers for the efficacy of IGF-1R targeting therapies in breast cancer patients. In addition to its therapeutic and clinical scope, this PhD project has significantly contributed to the understanding of the role of the IGF system in breast tumour biology by providing valuable new information on the mechanistic events underpinning IGF-I:VN-mediated effects on breast cell functions. Furthermore, this is the first instance where favourable binding sites for IGF-II, IGFBP-3 and IGFBP-5 on VN have been identified. Taken together, this study has functionally characterised the interactions between IGF-I and VN and through innovative strategies has provided a platform for the development of novel therapies targeting these interactions and their downstream effects.

Relative abundance of full-length and truncated FOXP1 isoforms is associated with differential NFκB activity in Follicular Lymphoma

FOXP1 is a transcriptional repressor that has been proposed to repress the expression of some NFκB-responsive genes. Furthermore, truncated forms of FOXP1 have been associated with a subtype of Diffuse Large B-cell Lymphoma characterised by constitutive NFκB activity, indicating that they may inhibit this repression. We have shown that FL tumors have increased relative abundance of truncated FOXP1 isoforms and this is associated with increased expression of NFκB-associated genes. Our results provide strong evidence that relative FOXP1 isoform abundance is associated with NFκB activity in FL, and could potentially be used as a marker for this gene signature.

Estrogen receptor β activation impairs mitochondrial oxidative metabolism and affects malignant mesothelioma cell growth in vitro and in vivo

Estrogen receptor (ER)-β has been shown to possess a tumor suppressive effect, and is a potential target for cancer therapy. Using gene-expression meta-analysis of human malignant pleural mesothelioma, we identified an ESR2 (ERβ coding gene) signature. High ESR2 expression was strongly associated with low succinate dehydrogenase B (SDHB) (which encodes a mitochondrial respiratory chain complex II subunit) expression. We demonstrate that SDHB loss induced ESR2 expression, and that activated ERβ, by over-expression or by selective agonist stimulation, negatively affected oxidative phosphorylation compromising mitochondrial complex II and IV activity. This resulted in reduced mitochondrial ATP production, increased glycolysis dependence and impaired cell proliferation. The observed in vitro effects were phenocopied in vivo using a selective ERβ agonist in a mesothelioma mouse model. On the whole, our data highlight an unforeseen interaction between ERβ-mediated tumor suppression and energy metabolism that may be exploited to improve on the therapy for clinical management of malignant mesothelioma.

Epithelial mesenchymal transition traits in human breast cancer cell lines parallel the CD44HI/CD24lO/-stem cell phenotype in human breast cancer

We review here the recently emerging relationship between epithelial-mesenchymal transition (EMT) and breast cancer stem cells (BCSC), and provide analyses of published data on human breast cancer cell lines, supporting their utility as a model for the EMT/BCSC state. Genome-wide transcriptional profiling of these cell lines has confirmed the existence of a subgroup with mesenchymal tendencies and enhanced invasive properties ('Basal B'/Mesenchymal), distinct from subgroups with either predominantly luminal ('Luminal') or mixed basal/luminal ('Basal A') features (Neve et al. Cancer Cell, 2006). A literature-derived EMT gene signature has shown specific enrichment within the Basal B subgroup of cell lines, consistent with their over-expression of various EMT transcriptional drivers. Basal B cell lines are found to resemble BCSC, being CD44highCD24low. Moreover, gene products that distinguish Basal B from Basal A and Luminal cell lines (Basal B Discriminators) showed close concordance with those that define BCSC isolated from clinical material, as reported by Shipitsin et al. (Cancer Cell, 2007). CD24 mRNA levels varied across Basal B cell lines, correlating with other Basal B Discriminators. Many gene products correlating with CD24 status in Basal B cell lines were also differentially expressed in isolated BCSC. These findings confirm and extend the importance of the cellular product of the EMT with Basal B cell lines, and illustrate the value of analysing these cell lines for new leads that may improve breast cancer outcomes. Gene products specific to Basal B cell lines may serve as tools for the detection, quantification, and analysis of BCSC/EMT attributes.