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

A simple and robust method for connecting small-molecule drugs using gene-expression signatures

Background

Interaction of a drug or chemical with a biological system can result in a gene-expression profile or signature characteristic of the event. Using a suitably robust algorithm these signatures can potentially be used to connect molecules with similar pharmacological or toxicological properties by gene expression profile. Lamb et al first proposed the Connectivity Map [Lamb et al (2006), Science 313, 1929–1935] to make successful connections among small molecules, genes, and diseases using genomic signatures.

Results

Here we have built on the principles of the Connectivity Map to present a simpler and more robust method for the construction of reference gene-expression profiles and for the connection scoring scheme, which importantly allows the valuation of statistical significance of all the connections observed. We tested the new method with two randomly generated gene signatures and three experimentally derived gene signatures (for HDAC inhibitors, estrogens, and immunosuppressive drugs, respectively). Our testing with this method indicates that it achieves a higher level of specificity and sensitivity and so advances the original method.

Conclusion

The method presented here not only offers more principled statistical procedures for testing connections, but more importantly it provides effective safeguard against false connections at the same time achieving increased sensitivity. With its robust performance, the method has potential use in the drug development pipeline for the early recognition of pharmacological and toxicological properties in chemicals and new drug candidates, and also more broadly in other 'omics sciences.

Association of gene expression with sequential proliferation, differentiation and tumor formation in murine skin

Differential gene expression in two established initiation and promotion skin carcinogenesis models during promotion and tumor formation was determined by microarray technology with the purpose of distinguishing the genes more associated with neoplastic transformation from those linked with proliferation and differentiation. The first model utilized dimethylbenz[a]anthracene initiation and 12-O-tetradecanoylphorbol 13-acetate (TPA) promotion in the FVB/N mouse, and the second TPA promotion of the Tg.Ac mouse, which is endogenously initiated by virtue of an activated Ha-ras transgene. Comparison of gene expression profiles across the two models identified genes whose altered expression was associated with papilloma formation rather than TPA-induced proliferation and differentiation. DMBA suppressed TPA-induced differentiation which allowed identification of those genes associated more specifically with differentiation rather than proliferation. EASE (Expression Analysis Systemic Explorer) indicated a correlation between muscle-associated genes and skin differentiation, whereas genes involved with protein biosynthesis were strongly correlated with proliferation. For verification the altered expression of selected genes were confirmed by RT-PCR; Carbonic anhydrase 2, Thioredoxin 1 and Glutathione S-transferase omega 1 associated with papilloma formation and Enolase 3, Cystatin 6 and Filaggrin associated with TPA-induced proliferation and differentiation. In situ analysis located the papillomas Glutathione S-transferase omega 1 expression to the proliferating areas of the papillomas. Thus we have identified profiles of differential gene expression associated with the tumorigenesis and promotion stages for skin carcinogenesis in the mouse.

A statistical framework for the design of microarray experiments and effective detection of differential gene expression

Motivation: Microarray experiments generate a high data volume. However, often due to financial or experimental considerations, e.g. lack of sample, there is little or no replication of the experiments or hybridizations. These factors combined with the intrinsic variability associated with the measurement of gene expression can result in an unsatisfactory detection rate of differential gene expression (DGE). Our motivation was to provide an easy to use measure of the success rate of DGE detection that could find routine use in the design of microarray experiments or in post-experiment assessment.

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.

Advanced glycation end products cause increased CCN family and extracellular matrix gene expression in the diabetic rodent retina

Aims/hypothesis: Referred to as CCN, the family of growth factors consisting of cystein-rich protein 61 (CYR61, also known as CCN1), connective tissue growth factor (CTGF, also known as CCN2), nephroblastoma overexpressed gene (NOV, also known as CCN3) and WNT1-inducible signalling pathway proteins 1, 2 and 3 (WISP1, -2 and -3; also known as CCN4, -5 and -6) affects cellular growth, differentiation, adhesion and locomotion in wound repair, fibrotic disorders, inflammation and angiogenesis. AGEs formed in the diabetic milieu affect the same processes, leading to diabetic complications including diabetic retinopathy. We hypothesised that pathological effects of AGEs in the diabetic retina are a consequence of AGE-induced alterations in CCN family expression.

Materials and methods: CCN gene expression levels were studied at the mRNA and protein level in retinas of control and diabetic rats using real-time quantitative PCR, western blotting and immunohistochemistry at 6 and 12 weeks of streptozotocin-induced diabetes in the presence or absence of aminoguanidine, an AGE inhibitor. In addition, C57BL/6 mice were repeatedly injected with exogenously formed AGE to establish whether AGE modulate retinal CCN growth factors in vivo.

Results: After 6 weeks of diabetes, Cyr61 expression levels were increased more than threefold. At 12 weeks of diabetes, Ctgf expression levels were increased twofold. Treatment with aminoguanidine inhibited Cyr61 and Ctgf expression in diabetic rats, with reductions of 31 and 36%, respectively, compared with untreated animals. Western blotting showed a twofold increase in CTGF production, which was prevented by aminoguanidine treatment. In mice infused with exogenous AGE, Cyr61 expression increased fourfold and Ctgf expression increased twofold in the retina.

Conclusions/interpolation: CTGF and CYR61 are downstream effectors of AGE in the diabetic retina, implicating them as possible targets for future intervention strategies against the development of diabetic retinopathy.

Immune activation in Retinal Aging: A Gene Expression Study.

PURPOSE
To investigate changes in gene expression during aging of the retina in the mouse.

METHODS
Total RNA was extracted from the neuroretina of young (3-month-old) and old (20-month-old) mice and processed for microarray analysis. Age-related, differentially expressed genes were assessed by the empiric Bayes shrinkagemoderated t-statistics method. Statistical significance was based on dual criteria of a ratio of change in gene expression >2 and a P < 0.01. Differential expression in 11 selected genes was further verified by real-time PCR. Functional pathways involved in retinal ageing were analyzed by an online software package (DAVID-2008) in differentially expressed gene lists. Age-related changes in differential expression in the identified retinal molecular pathways were further confirmed by immunohistochemical staining of retinal flat mounts and retinal cryosections.

RESULTS
With ageing of the retina, 298 genes were upregulated and 137 genes were downregulated. Functional annotation showed that genes linked to immune responses (Ir genes) and to tissue stress/injury responses (TS/I genes) were most likely to be modified by ageing. The Ir genes affected included those regulating leukocyte activation, chemotaxis, endocytosis, complement activation, phagocytosis, and myeloid cell differentiation, most of which were upregulated, with only a few downregulated. Increased microglial and complement activation in the aging retina was further confirmed by confocal microscopy of retinal tissues. The most strongly upregulated gene was the calcitonin receptor (Calcr; >40-fold in old versus young mice).

CONCLUSIONS
The results suggest that retinal ageing is accompanied by activation of gene sets, which are involved in local inflammatory responses. A modified form of low-grade chronic inflammation (para-inflammation) characterizes these aging changes and involves mainly the innate immune system. The marked upregulation of Calcr in ageing mice most likely reflects this chronic inflammatory/stress response, since calcitonin is a known systemic biomarker of inflammation/sepsis. © Association for Research in Vision and Ophthalmology.