- Research article
- Open Access
A searchable cross-platform gene expression database reveals connections between drug treatments and disease
BMC Genomics volume 13, Article number: 12 (2012)
Transcriptional data covering multiple platforms and species is collected and processed into a searchable platform independent expression database (SPIED). SPIED consists of over 100,000 expression fold profiles defined independently of control/treatment assignment and mapped to non-redundant gene lists. The database is thus searchable with query profiles defined over genes alone. The motivation behind SPIED is that transcriptional profiles can be quantitatively compared and ranked and thus serve as effective surrogates for comparing the underlying biological states across multiple experiments.
Drug perturbation, cancer and neurodegenerative disease derived transcriptional profiles are shown to be effective descriptors of the underlying biology as they return related drugs and pathologies from SPIED. In the case of Alzheimer's disease there is high transcriptional overlap with other neurodegenerative conditions and rodent models of neurodegeneration and nerve injury. Combining the query signature with correlating profiles allows for the definition of a tight neurodegeneration signature that successfully highlights many neuroprotective drugs in the Broad connectivity map.
Quantitative querying of expression data from across the totality of deposited experiments is an effective way of discovering connections between different biological systems and in particular that between drug action and biological disease state. Examples in cancer and neurodegenerative conditions validate the utility of SPIED.
There is a wealth of deposited gene expression data available for downloading and/or online interrogation. For example, the NCBI gene expression omnibus (GEO)  hosts over half a million single array chip expression profiles and the EBI hosts the ArrayExpress  database with a similar largely overlapping number of arrays. These data cannot be compared directly as they come from different array platforms covering many different species and a variety of normalisation schemes are used. In the overwhelming number of analyses expression profiles are compared within the given series and probed for the up or down regulation of single genes using volcano plot representations or other statistical filters . Alternatively, a larger set of responders can be scored against gene sets corresponding to pathways , interacting networks [5–7] or gene ontology classes [8, 9]. For large series it is possible to compile correlations of expression changes of individual gene pairs and groups of genes leading to a hierarchical clustering based network discovery and gene interaction prediction. To this end SOURCE  hosts gene expression profiles across a large collection of experimental series and profile correlations within a given series can be examined to predict genes with similar or related function. Many array analysis applications incorporate array derived network data that are valuable aids in characterising the expression profile data (Ingenuity Pathways Analysis (IPA) (Ingenuity® Systems, http://www.ingenuity.com); GeneGo (GeneGo, Inc., St. Joseph, MI)). However, these analyses do not allow for a direct quantitative comparison between separate expression studies and therefore a lot of the information contained in the experiment is effectively lost.
The idea that transcriptional change profiles can be directly compared to asses drug target specificity was demonstrated in yeast systems by Marton et al and later extended by Hughes et al. The connectivity map (CMAP) project sought to apply these ideas to generate a database of drug perturbagen transcriptional profiles that can be searched with transcriptional responder sets by third parties to match phenotype to drug treatment . In this methodology the expression change profile as a whole defines the biological perturbation and not a relatively small selection of down- or up-regulated genes. An important point here is that biological effects are not necessarily caused by the corresponding transcriptional changes. Rather, the underlying assumption is that correlations in transcriptional change profiles are reflected in similar biological responses. One powerful application of the CMAP is the matching of disease state to drug treatment. In simple terms, if a disease state is reflected in a well defined transcriptional response, then a drug that has the opposite effect on expression of these transcripts might be of therapeutic value. The fundamental assumption here is that there is a degree of overlap in the transcriptional changes induced by the same perturbagen in different cell contexts. In particular the CMAP consist of expression change data for human cancer cell lines and it is hoped that there is a degree of universality that will enable useful predictions to be made as to the action of the drugs in different cell types. Of course, the successful application of the CMAP should encourage rather than hinder the inclusion of other cell types more relevant to the type of biological system under investigation. At the present the CMAP consists of expression change fold profiles for 6,100 single treatments versus control pairs for a collection of 1,309 drug like perturbagens. Results are collected from treatments of four distinct types of human cancer cell lines. The CMAP database can be interrogated with expression change signatures consisting of lists of up and down regulated probe sets. Correlation both in the positive and negative sense are scored by means of a non-parametric Kolmogorov-Smirnov (KS) statistic . The remarkable observation was that signatures from published studies showed correlation with CMAP profiles for drugs known to act against the same targets. This has opened the way for the CMAP to be used as a drug discovery tool where it is probed with signatures encoding disease states.
If the CMAP methodology is accepted as a useful discovery tool then it is natural to look for ways of extending it to incorporate expression data from a wider set of experiments. There are obvious advantages to having this kind of database, for example it will open up a large number of different samples and treatment conditions for direct interrogation. This was the idea behind GEM-TREND , where 26,000 expression samples from various platforms and species were compiled into a searchable database. The search methodology mirrors that of CMAP in that the database consists of ranked lists of genes and it is interrogated with up and down regulated gene sets and query signatures are scored by a KS statistic with the significance based on reference to random gene set scores. One difference to the CMAP database is necessitated by the multiple origins of the expression profile data represented by multiple probe ID definitions. The problem of multiple probe IDs is solved by the GEM-TREND database having expression profiles mapped onto UniGene IDs. The database consists of experimental series where samples can be clearly assigned to treatment and control groups. Of course, this is not always the case and this limits the scope of the database.
In compiling the expression database SPIED we sought to loosen the restraints inherent in previous treatments and thereby open up a larger set of data for interrogation. In many expression series sets there is no clear control/treatment assignment or there could be multiple alternative reference profile definitions. To address this problem of generating fold change profiles without reference to a defined control, an effective fold (EF) has been introduced corresponding to the expression level relative to the experimental series average. In this way, data can be compiled automatically without the need for manual inspection. In cases where the experimental series consists of well defined multiple treatment and control samples the fold profiles are usually given by the ratio of the average treatment to average control values. In general this fold profile will have high positive correlation with the EF profiles from the treatment set and high negative correlations with the control set. In cases where there is no obvious way of separating samples into control and treatment sets, as with samples from multiple organ types or cell types, the EF representation can be viewed as a normalized expression value. In searching SPIED with a query profile one is not deriving any biological significance for non-correlating profiles as lack of correlation can be attributed to multiple factors such as bad experimental data or genuine lack of biological relevance. Rather significantly correlating or anti-correlating profiles are posited as having biological significance. The next objective was to reduce the expression profiles to non-redundant EF gene profiles by associating each gene with just one probe ID, so that the database can then be searched with gene set data alone. Here, for a given chip platform the distribution of each probe ID EF value across the totality of series was compiled and each gene was then assigned to the probe having the highest average fold magnitude. The gene names were unambiguously associated with the Entrez human gene list (http://www.ncbi.nlm.nih.gov), consisting of 24,764 genes and these were matched to probe IDs by inspection of the given platform annotation files. The final form of SPIED consists of individual files for each chip platform and these files are formatted starting with a gene list followed by the sample ID and corresponding EF profiles. This format lends itself to rapid searching in an analogous fashion to FASTA formatted sequence databases. In contrast to the KS query score scheme, which requires generating random reference gene list data, we adopted a simple regression scoring scheme with corresponding statistic. Searches can be performed on a standard desktop PC and take ~10 minutes per query. Although, the present database consisting of expression data for over 100,000 samples from five platforms covering three species (human, mouse and rat) is all from Affymetrix expression array chips, the methodology is truly platform independent and it is a straight forward matter to include data based on other array technologies. Other species and platform technologies will be added to SPIED in the future. For the present study Affymetrix was chosen because of the relatively large number of available samples. Further details are presented in the methods section below.
Drug treatment based profile SPIED queries
The CMAP contains expression change profiles as ranked array probe IDs for 6,100 individual treatments corresponding to 1,309 distinct drug-like compounds. Statistically filtered response profiles can be defined for 1,218 of the drugs as these have at least three instances in the database. The profiles can be mapped onto a non-redundant gene list by uniquely associating one probe ID to a given gene and dropping the other probe ID for this gene with less robust expression changes over the database. This is the same methodology underlying the SPIED database. We took the responder profiles for the 1,218 drugs and searched the SPIED for maximally (anti-)correlated expression change profiles. The objective is to see to what extent the CMAP transcriptional signatures correlate with transcriptional responses assimilated within our platform independent database of over 100,000 microarrays deposited by a very large number of groups to the public domain.
The CMAP is well populated with drugs that target the same or different steps in the PI3K-mTOR signalling cascade. In this context the results for LY-294002 (PI3K inhibitor) (61 instances), rapamycin (mTOR antagonist) (44 instances) and wortmannin (PI3K/mTOR antagonist) (18 instances) showed a high degree of overlap (rapamycin v LY-294002 r = 0.90 N = 3565, rapamycin v wotmannin r = 0.91 N = 1849, LY-294002 v wotmannin r = 0.88 N = 2217, where r is the Pearson correlation coefficient and N is the number of shared genes with significant fold values), see additional file 1 for the full fold change data. It is a straightforward matter to query the SPIED with these drug expression profiles. This is done by calculating the regression scores against the individual SPIED entries and retaining the top ~100 correlations, see Methods for details. For simplicity and uniformity of treatment, unless otherwise stated, we query SPIED with expression profiles containing 500 genes with the largest fold values passing the p < 0.05 significance threshold. It should be noted that results will be largely insensitive to the size of the query profile. The top SPIED correlate for all three drugs was the Pan-PI3K inhibitor GDC-0941 treated T47D breast cancer cells and the regression scores for the tree query signatures against all 6 samples in the series are shown in Figure 1A. The high degree of correlation is illustrated by regression plots for the three query profiles against the pooled GDC-0941 profile, see Figure 1B, C, D. All three inhibitor queries also pick out mTOR antagonist studies , but a more interesting correlation is with a glucocorticoid (dexamethasone) treatment of acute lymphoblastic leukaemia (ALL) cells , the rapamycin scores are shown in Figure 2A. The correlation increases with the length of drug treatment, being higher at 24 hours, Figure 2B, C. This result reveals another connection between mTOR antagonism and the corticosteroid mechanism as it has been shown that corticosteroid resistance in ALL can be overcome by mTOR antagonism . Chronic myeloid Leukaemia (CML) and some instances of ALL are the result of the ABL tyrosine kinase translocation and fusion to BCR, the BCR-ABL fusion event . This pathology has been targeted with rapamycin and our results support this approach based on the high degree of anti-correlation of the CMAP rapamycin profile with a transcriptional profile of BCR fusion construct transformed chord blood cells. The correlation scores are shown in Figure 3A. There is a clear anti-correlation of rapamycin profile with the BCR-ABL profiles pointing to a possible reversal of the phenotype, Figure 3B. Also, there is a high anti-correlation with the BCR-FGFR1 profile indicating a possible therapeutic role of rapamycin, Figure 3C.
In the original CMAP presentation  it was shown that meaningful results can be obtained from anti-correlating profiles. In particular the estrogen transcriptional response was shown to anti-correlate with the profiles of estrogen antagonists fulvestrant, tamoxifen and raloxifene. In this context it is of interest to note that high scoring SPIED hits for all three antagonists corresponded to anti-correlations with estrogen treatment samples. We have shown one example in Table 1 corresponding to a estrogen, tamoxifen and an extract from the cimicifuga plant .
For illustration purposes we have shown the common high correlating hits for three separate histone deacetylase (HDAC) inhibitor profiles in the CMAP series. These are vorinistat, trichostatin A and valporic acid. In Table 2 we have shown the regression scores for the multiple HDAC inhibitor study with a colorectal carcinoma cell line . The query results for all the above searches are given in additional file 2.
Next we consider profiles derived from disease states. For brevity we focus on two unrelated pathologies: cancer and neurodegeneration.
Querying SPIED with cancer derived profiles
The class of diseases with the most extensive repository of expression data is cancer and therefore a cancer disease profile search of SPIED will be an ideal testing ground for the methodology. The original CMAP disease application implicated mTOR inhibition as a target for imparting sensitivity to dexamethasone treatment resistant ALL . We searched the SPIED database with the dexamethasone resistant v sensitive profile to see if there are common features in published transcriptional studies. The query profile consisted of the 500 most highly regulated genes that passed the lowest significance test of p < 0.05, see additional file 1. As with the SPIED profiles the query profile also consists of a non-redundant gene list. Not surprisingly, the highest correlation scores came from the experiments from which the query profile was generated, see additional file 2 file. In addition, we found a high correlation to an independent later study of ALL sensitivity to corticosteroid (prednisolone) treatment . This study generated transcriptional profiles of ALL patient leukaemia cells with the objective of uncovering a gene signature that can predict the sensitivity to prednisolone treatment. Combining the 27 infant and non-infant corticosteroid sensitive samples and the 25 resistant samples we can define a statistically filtered sensitivity profile to make a direct comparison with the query profile and we find a high degree of correlation (r = 0.94), see additional file 2. When the high scoring sample belongs to a relatively large sample series and the phenotype is binary we can perform a non-parametric significance test to measure the extent of enrichment of the given phenotype for high or low correlation scores. For example in the last case there were 25 resistant and 27 sensitive samples. Ranking the samples according to their correlation with the resistant versus sensitive query profile we find 20 resistant samples in the top 25 and 22 sensitive samples in the bottom 27. This is highly significant and can be quantified with a simple Fisher exact test. Explicitly, the probability p of 20 or more resistant samples in the top 25 correlations is less than 9 × 10-7. The K-S significance score can be calculated by counting the number of times a random rearrangement of the samples gives a better enrichment, we find p < 3 × 10-6. The enrichment plot is given in Figure 4A. As expected the top scoring correlations were dominated by samples from blood derived cells, for simplicity we restricted our analysis to the top 100 most significantly correlating samples. However, two studies in unrelated tissue pathologies were highly correlated with the corticosteroid resistant profile. These were a comparison of lung epithelia with cancer in smokers  and a differential expression between healthy and cancerous pancreatic tissue . The smoking study consisted of non-diseased lung epithelia from 187 individual smokers 97 of whom were diagnosed with lung cancer. Ranking the samples according to query correlation score we find that in the top 97 there are 64 cancer cases and in the bottom 90 there are 57 non-cancer cases, with a significance score of p < 5 × 10-5. The K-S significance is p < 2 × 10-4. The enrichment for positive correlations with the corticosteroid resistance profile in the cancer cases is shown in Figure 4B. Interestingly, it has been shown that there is a down regulation of the glucocorticoid receptor (GR) in small cell lung cancer [25, 26] and reversing this promotes cancer cell apoptosis . The pancreatic cancer study sought to establish a transcriptional signature of tumour versus normal pancreatic tissue by laser capture of cancerous and normal tissue from the same pancreas. In total 39 sample pairs were published and we find a high positive correlation with the corticosteroid resistance profile, p < 2 × 10-6 and a K-S significance score of p < 3 × 10-7. The enrichment curve is shown in Figure 4C. In this context it has been reported that loss of GR expression has been seen in pancreatic carcinoma relative to normal tissue  and elevating GR expression has been shown to inhibit pancreatic tumour growth in a hamster model . The query results are given in the additional file 2.
The analysis of gene expression changes associated with neurodegenerative disease has been hampered by the difficulty of extracting high quality RNA from post-mortem tissue [30, 31]. One way of validating a disease associated gene expression profile is to show that it shares significant features with profiles derived from independent experiments on related pathologies. A positive result would validate the query profile and furthermore lead to a more robust core response profile based on multiple experiments. To this end we constructed three separate query profiles based on transcriptional profiles from the brains of patients with three degrees of severity of Alzheimer's disease (AD) , see additional file 1. The number of significant changes increases with severity of disease and we queried the SPIED with these three profiles, see additional file 2. Not surprisingly, the high scoring correlations are those from which the query profiles were derived. In addition to these the query returned correlations with other AD studies and various neurodegenerative diseases. The high scoring AD expression series was an extensive study of 161 samples from various brain regions of AD patients and age-matched controls . Ignoring brain regions for now, there are 87 AD samples and 74 controls. Ranking the samples according to correlation score against the severe AD query profile we find a very significant enrichment of positive correlations with AD samples (p < 10-8, based on the Fisher test as above). Pooling the samples from the different brain regions results in significant correlations for 5 out of the 6 brain regions, see Figure 5.
In addition to AD correlations we found high scoring correlations with samples derived from Huntington's disease (HD), Down's syndrome (DS), Parkinson's disease (PD) and bipolar disorder (BD) brains. In this sense the profile cannot be considered to distinguish AD pathology from other degenerative diseases. However, it is of interest to examine in greater detail these cross-disease similarities. In particular, the severe AD query had a high correlation with severe stage HD caudate nucleus (CN) samples. The HD study consisted of 404 samples split across two platforms (201 samples on GPL96 and 203 on GPL97) in three brain regions from control and HD individuals . The high correlation was with the GPL96 series. In terms of a binary Fisher analysis where brain region specificity is ignored, we get a small enrichment of p < 6 × 10-3. However, when the different brain regions are considered separately, we get significant regression scores in each region. The results are tabulated in Table 3.
The PD correlation was with a study of 94 samples from three different regions of diseased and normal brains . Pooling samples according to brain region we find that the severe AD profile had a high correlation with all three regions studied: superior frontal gyrus (SFG) r = 0.88; lateral substatia nigra (LSN) r = 0.77; medial substantia nigra (MSN) r = 0.82, see Table 3.
The chromosome 21 trisomy underlying DS leads to the development of many of the characteristics of AD pathology [36, 37]. Therefore, it is not surprising to find a high correlation in SPIED form a transcriptional profiling of DS brains. This study comprised 8 healthy and 7 DS individual brains . Combining the expression data into a thresholded fold change profile we find that there is a significant but small positive correlation with the severe AD profile, with r = 0.58. Interestingly, the correlation is higher with the moderate AD profile, with r = 0.68, see Table 3.
The first transcriptional profiling of BD brains pointed to the down regulation of synaptic and mitochondrial proteins in the orbital frontal cortex . This synaptic pathology picture of BD is further strengthened by our analysis of the AD profile correlates within SPIED. Pooling the 10 BD samples and 11 controls we find a high regression score with the severe AD profile, see Table 3. It is important to note that this correlation is with a subset of the BD signature as it consists of genes that are also altered in AD. However, it is outside the scope of the present paper to combine profiles into disease specific queries.
Not surprisingly the high correlations are dominated by experiments on human samples. Perhaps of greater interest to the biologist are animal models of neurodegeneration. There has indeed been a debate as to the relevance of animal models of neurodegeneration to drug discovery, as age-related neurodegenerative conditions are rare in nonhumans . In this context it is interesting to look at what correlations the AD query profile returns when we restrict the search to rodent platforms. The SPIED database contains samples from two murine and one rat platforms. Within the top 100 high scoring samples we have four separate studies directly relevant to neuropathology, see additional file 2. In particular, we find high scores with two separate spinal contusion models. The mouse experiments generated a post injury expression time series (GEO accession GSE5296, also http://pepr.cnmcresearch.org) and the AD profile correlation emerges at 72 hours post injury, see Table 4. The other spinal chord contusion study was in rats at 35 days post injury , see Table 4. In addition to these contusion models high scores were for a murine SOD1(G93A) mutant model of Amyotrophic lateral sclerosis (ALS) (GEO accession GSE18597) and a murine model of prion disease (GEO accession GSE23182). In the SOD1 (G93A) transcriptional profile series we found the correlation with AD emerging with older mice, with negligible correlation at the 28-70 day window and significant correlation with the 98-126 day late stage window profiles. This is consistent with the timescale of disease onset in the mouse model [42, 43]. Prion disease is modelled in mice through ME7 prion agent infection resulting in both a behavioural phenotype and synaptopathy . The transcriptional study corresponded to hippocampal profiles for ME7 v normal brain homogenate inoculated mice (GEO accession GSE23182). Pooling the treatment sets we get a good correlation with the AD profile, see Table 4. Thus it is clear that there is a core response profile shared across many neurodegenerative conditions and animal models of these conditions. Importantly, this core set is characterised by synaptic pathology and mitochondrial dysfunction, both of which are hypothesised to be causative of a number of neurodegenerative disease states.
It might be thought that we are getting further away from the specific pathology, in this case AD, and losing transcriptional information that could be of use in the hunt for a therapy. This is however not the case as can be seen when we search the CMAP with a profile composed of genes whose sense change is conserved across the rodent disease models. Combining the severe AD profile and the four rodent neurodegenerative disease model profiles we get a set of 24 genes whose sense change is conserved. This consists of 10 up regulated and 14 down regulated genes, which can be thought of as a binary signature for neuropathology, where +1 is assigned to up-regulated genes and -1 to down-regulated genes, see Table 5. The CMAP drugs with the highest anti-correlation with this signature are shown in Table 6. Remarkably, there are at least 9 neuroprotective agents in the top 22 hits. In particular, Galantamine, a plant alkaloid, is currently prescribed for early stage AD [45–48], it was originally studied for its acetylcholinesterase inhibitory activity, but it may also act on other targets . The flavones chrysin , apigenin [50, 51] and luteolin [50, 52, 53] have been reported to have neuroprotective activity. As have the two kinase inhibitors H-7 (PKA/PKC inhibitor) and GW-8510 (CDK inhibitor). The β-carboline plant alkaloid harmine has several neuronal actions. It acts to slow down the breakdown of monoamine neurotransmitters through inhibition of monoamine oxidase A [56, 57]. Also, it has been shown to specifically inhibit DYRK1A, an enzyme responsible for phosphorylation of tau and thereby may act to slow tau pathology in AD and DS [58, 59]. Nomifensine is a dopamine reuptake inhibitor originally prescribed as an anti-depressant  that has been shown to reverse dopaminergic neurotoxicity [61–63] and to have beneficial effects in Parkinson's disease . Carbachol is an acetylcholine receptor agonist, but with poor blood brain barrier penetration . The possible application of the other high scoring compounds remains to be determined.
Discussion and Conclusions
We have collected transcriptional data from diverse platform architectures corresponding to various species. By processing the data into effective fold profiles, with the expression levels factored by the average level over the experimental series and defined over a non-redundant gene list, we can directly compare transcriptional profiles from arbitrary sources. The fundamental principal underlying the utility of this approach is that biological effects can be compared through the corresponding transcriptional changes. This idea underlies the CMAP initiative for matching drug to phenotype by querying a database of drug induced transcriptional profiles with a profile defining the phenotype. We have extended this methodology to include potentially all available transcriptional data. In its current version SPIED contains transcriptional profiles for 106,101 arrays covering five platform architectures and three species. This can be easily extended to include other platforms and species. The results largely confirm the hypothesis that high scoring correlations correspond to similar biological processes. We have presented SPIED results for drug perturbagen induced profile queries and queries derived from disease states. For brevity we focussed on three sets of drug treatment profiles corresponding to mTOR/PI3K, estrogen and HDAC inhibitors. SPIED searches with these queries showed correlations with other drug treatments belonging to the same classes and in the case of the mTOR antagonist rapamycin we found high anti-correlations with the profile of a cancer inducing fusion transformation, suggesting a novel indication for rapamycin. Also, for brevity of exposition we focussed on two completely unrelated classes of pathology: cancer and neurodegeneration. In the case of leukaemia we show that a corticosteroid resistance signature derived from leukaemia cell cultures shows significant correlation with a lung cancer predisposition profile and a pancreatic cancer profile. Thereby implicating glucocorticoid resistance in these two pathologies. To illustrate the application of SPIED to neurodegenerative pathology we constructed a severe stage AD profile from a published study. Interestingly, querying SPIED resulted in high correlations with other neuropathological conditions indicating a common feature of synaptic loss and mitochondrial dysfunction. Restricting our searches to the rodent subset of SPIED returned expression profiles from animal models of neurodegeneration and neuronal injury. Combining the human and rodent signatures we obtained a core signature that we probed against CMAP for neuroprotective agents. Remarkably, we found at least 9 neuroprotective agents in the top 22 anti-correlating CMAP hits. These results motivate the extension of SPIED and the extension of the CMAP to include other cell types, for example a neuronal cell lineage will be more appropriate for generating drug profiles for neurological diseases. The correlation query scores maybe insensitive to a radical reduction in the number of probes and this should motivate the design of reduced and more cost effective arrays for more extensive data generation.
Compiling the data
Microarray sample files, GSM files, were downloaded form the NCBI GEO database. Individual GSM files were assigned to GSE series and log scaled values scaled to linear and low level responders dropped. EF profiles were then generated based on ratio of individual condition to the average across the series. Expression data from five Affymetrix GeneChip platforms corresponding to three species were collected. These were: all samples from two Human array platforms corresponding to Human Genome U133 Array Set HG-U133A GPL96 (27,337 samples) and U133 Plus 2.0 Array GPL570 (55,196 samples); all samples from the Mouse Genome 430 2.0 Array GPL1261 (21,219 samples) chip; all samples from two Rat chips corresponding to Rat Genome U34 Array GPL85 (2,827 samples) and Rat Genome 230 2.0 Array GPL1355 (7,476 samples). The database thus totals 106,101 samples. Of course, this can always be extended to include more platforms from the same species and/or other species.
The individual GSM sample file expression values were transformed into EF values corresponding to the expression relative to the series mean. Expression values that have been logarithmically transformed are transformed back to a linear scale and low expression values dropped, that is are set to zero and don't contribute to the fold profile. We found that the results were relatively insensitive to the cut-off value and we set this to be 10% of the average expression value. All sample expression profiles within a series were scaled to the same average. The fold vales are defined as, where s k is the expression level of the kthprobe set and is the average over the series. For the database to be searchable with cross platform response profiles and gene lists it has to be rewritten as a database of expression profiles over non-redundant gene lists. The EF profiles across the probe sets were therefore mapped onto expression profiles for a non-redundant gene list. In general each gene is represented by multiple probe sets. For each platform we generated the EF statistics for each probe set across the totality of samples. The probe set with the most robust response across the samples was chosen to represent the gene. Explicitly, the probe set with the highest root mean square deviation form zero was chosen to represent the given gene. The number of genes defined on each platform were as follows: GPL96 11,807, GPL570 15,983 genes, GPL1261 13,202 genes, GPL85 chip with 3,844 genes, GPL1355 chip with 6,341 genes. The database totals 106,101 samples and is searchable on a reasonably fast desktop PC in ~10 minutes per query.
Searching the database
The query profile is a statistically thresholded non-redundant list of genes and associated fold values. Statistical significance is assigned to a fold change based on a simple Student's t-test between multiple control and treatment sample expression values. This is compared to each profile in the database by means of a simple Pearson regression analysis, with a correlation coefficient r. The experiments are ranked according to the significance. The significance is measured by scaling the correlation to the normal by a Fisher transformation and measuring the number of standard deviations from the mean. The Fisher transformation is and the standard deviation is, where r is the Pearson correlation coefficient and N is the number of genes making up the correlation. The final ranking score is .
CMAP combined profiles
The CMAP contains ranked lists of probes for 6,100 separate perturbagen treatments of four different human cell lines, with the ranking based on response level relative to control. The treatments are various multiples of 1,306 different drug-like compounds. To generate responder sets that can be used to search SPIED we combined rankings for each separate compound treatment and converted these into pseudo-fold values with associated statistics. The pseudo-fold value is defined by, where r i is the rank of the ith gene and min/max are the minimal/maximal ranks. Remembering that the highest rank corresponds to the most up-regulated gene. The SPIED was searched with CMAP profiles corresponding to folds with a p < 0.05 threshold and with at least three replicates. This left 1,218 separate perturbagen probes. We sought to cluster the perturbagens based on predicted target and response profile similarity. The profiles are given in the additional file 1 file.
Availability of SPIED
The SPIED database and associated executables are available for download from ftp://ftp.hostedftp.com/~GWftpFILES/SPIED/. The download consists of the SPIED database together with executables for searching SPIED. Source code files to generate the database and perform query searches are provided together with the executables. Documentation on the database, the executables and source code files is also included.
Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Edgar R: NCBI GEO: mining tens of millions of expression profiles--database and tools update. Nucleic Acids Res. 2007, D760-765. 35 database
Parkinson H, Sarkans U, Kolesnikov N, Abeygunawardena N, Burdett T, Dylag M, Emam I, Farne A, Hastings E, Holloway E: ArrayExpress update--an archive of microarray and high-throughput sequencing-based functional genomics experiments. Nucleic Acids Res. 2011, D1002-1004. 39 database
Cui X, Churchill GA: Statistical tests for differential expression in cDNA microarray experiments. Genome Biol. 2003, 4 (4): 210-10.1186/gb-2003-4-4-210.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.
Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M: BioGRID: a general repository for interaction datasets. Nucleic Acids Res. 2006, D535-539. 34 database
Bader GD, Betel D, Hogue CW: BIND: the Biomolecular Interaction Network Database. Nucleic Acids Res. 2003, 31 (1): 248-250. 10.1093/nar/gkg056.
Xenarios I, Salwinski L, Duan XJ, Higney P, Kim SM, Eisenberg D: DIP, the Database of Interacting Proteins: a research tool for studying cellular networks of protein interactions. Nucleic Acids Res. 2002, 30 (1): 303-305. 10.1093/nar/30.1.303.
Zeeberg BR, Feng W, Wang G, Wang MD, Fojo AT, Sunshine M, Narasimhan S, Kane DW, Reinhold WC, Lababidi S: GoMiner: a resource for biological interpretation of genomic and proteomic data. Genome Biol. 2003, 4 (4): R28-10.1186/gb-2003-4-4-r28.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.
Diehn M, Sherlock G, Binkley G, Jin H, Matese JC, Hernandez-Boussard T, Rees CA, Cherry JM, Botstein D, Brown PO: SOURCE: a unified genomic resource of functional annotations, ontologies, and gene expression data. Nucleic Acids Res. 2003, 31 (1): 219-223. 10.1093/nar/gkg014.
Marton MJ, DeRisi JL, Bennett HA, Iyer VR, Meyer MR, Roberts CJ, Stoughton R, Burchard J, Slade D, Dai H: Drug target validation and identification of secondary drug target effects using DNA microarrays. Nat Med. 1998, 4 (11): 1293-1301. 10.1038/3282.
Hughes TR, Marton MJ, Jones AR, Roberts CJ, Stoughton R, Armour CD, Bennett HA, Coffey E, Dai H, He YD: Functional discovery via a compendium of expression profiles. Cell. 2000, 102 (1): 109-126. 10.1016/S0092-8674(00)00015-5.
Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, Lerner J, Brunet JP, Subramanian A, Ross KN: The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease. Science. 2006, 313 (5795): 1929-1935. 10.1126/science.1132939.
Hollander M, Wolfe D: Nonparametric Statistical Methods. 1973, New York: John Wiley & Sons, 2
Feng C, Araki M, Kunimoto R, Tamon A, Makiguchi H, Niijima S, Tsujimoto G, Okuno Y: GEM-TREND: a web tool for gene expression data mining toward relevant network discovery. BMC Genomics. 2009, 10: 411-10.1186/1471-2164-10-411.
Fernandez DR, Telarico T, Bonilla E, Li Q, Banerjee S, Middleton FA, Phillips PE, Crow MK, Oess S, Muller-Esterl W: Activation of mammalian target of rapamycin controls the loss of TCRzeta in lupus T cells through HRES-1/Rab4-regulated lysosomal degradation. J Immunol. 2009, 182 (4): 2063-2073. 10.4049/jimmunol.0803600.
Rainer J, Ploner C, Jesacher S, Ploner A, Eduardoff M, Mansha M, Wasim M, Panzer-Grumayer R, Trajanoski Z, Niederegger H: Glucocorticoid-regulated microRNAs and mirtrons in acute lymphoblastic leukemia. Leukemia. 2009, 23 (4): 746-752. 10.1038/leu.2008.370.
Wei G, Twomey D, Lamb J, Schlis K, Agarwal J, Stam RW, Opferman JT, Sallan SE, den Boer ML, Pieters R: Gene expression-based chemical genomics identifies rapamycin as a modulator of MCL1 and glucocorticoid resistance. Cancer Cell. 2006, 10 (4): 331-342. 10.1016/j.ccr.2006.09.006.
Chopra R, Pu QQ, Elefanty AG: Biology of BCR-ABL. Blood Rev. 1999, 13 (4): 211-229. 10.1054/blre.1999.0119.
Gaube F, Wolfl S, Pusch L, Kroll TC, Hamburger M: Gene expression profiling reveals effects of Cimicifuga racemosa (L.) NUTT. (black cohosh) on the estrogen receptor positive human breast cancer cell line MCF-7. BMC Pharmacol. 2007, 7: 11-
Liu Y, Salvador LA, Byeon S, Ying Y, Kwan JC, Law BK, Hong J, Luesch H: Anticolon cancer activity of largazole, a marine-derived tunable histone deacetylase inhibitor. J Pharmacol Exp Ther. 335 (2): 351-361.
Stam RW, Den Boer ML, Schneider P, de Boer J, Hagelstein J, Valsecchi MG, de Lorenzo P, Sallan SE, Brady HJ, Armstrong SA: Association of high-level MCL-1 expression with in vitro and in vivo prednisone resistance in MLL-rearranged infant acute lymphoblastic leukemia. Blood. 2010, 115 (5): 1018-1025. 10.1182/blood-2009-02-205963.
Spira A, Beane JE, Shah V, Steiling K, Liu G, Schembri F, Gilman S, Dumas YM, Calner P, Sebastiani P: Airway epithelial gene expression in the diagnostic evaluation of smokers with suspect lung cancer. Nat Med. 2007, 13 (3): 361-366. 10.1038/nm1556.
Badea L, Herlea V, Dima SO, Dumitrascu T, Popescu I: Combined gene expression analysis of whole-tissue and microdissected pancreatic ductal adenocarcinoma identifies genes specifically overexpressed in tumor epithelia. Hepatogastroenterology. 2008, 55 (88): 2016-2027.
Sommer P, Le Rouzic P, Gillingham H, Berry A, Kayahara M, Huynh T, White A, Ray DW: Glucocorticoid receptor overexpression exerts an antisurvival effect on human small cell lung cancer cells. Oncogene. 2007, 26 (50): 7111-7121. 10.1038/sj.onc.1210524.
Ray DW, Davis JR, White A, Clark AJ: Glucocorticoid receptor structure and function in glucocorticoid-resistant small cell lung carcinoma cells. Cancer Res. 1996, 56 (14): 3276-3280.
Sommer P, Cowen RL, Berry A, Cookson A, Telfer BA, Williams KJ, Stratford IJ, Kay P, White A, Ray DW: Glucocorticoid receptor over-expression promotes human small cell lung cancer apoptosis in vivo and thereby slows tumor growth. Endocr Relat Cancer. 17 (1): 203-213.
Norman J, Franz M, Schiro R, Nicosia S, Docs J, Fabri PJ, Gower WR: Functional glucocorticoid receptor modulates pancreatic carcinoma growth through an autocrine loop. J Surg Res. 1994, 57 (1): 33-38. 10.1006/jsre.1994.1105.
Evers BM, Thompson EB, Townsend CM, Lawrence JL, Johnson B, Srinivasan G, Thompson JC: Cortivazol increases glucocorticoid receptor expression and inhibits growth of hamster pancreatic cancer (H2T) in vivo. Pancreas. 1993, 8 (1): 7-14. 10.1097/00006676-199301000-00004.
Barton AJ, Pearson RC, Najlerahim A, Harrison PJ: Pre- and postmortem influences on brain RNA. J Neurochem. 1993, 61 (1): 1-11. 10.1111/j.1471-4159.1993.tb03532.x.
Webster MJ: Tissue preparation and banking. Prog Brain Res. 2006, 158: 3-14.
Blalock EM, Geddes JW, Chen KC, Porter NM, Markesbery WR, Landfield PW: Incipient Alzheimer's disease: microarray correlation analyses reveal major transcriptional and tumor suppressor responses. Proc Natl Acad Sci USA. 2004, 101 (7): 2173-2178. 10.1073/pnas.0308512100.
Liang WS, Dunckley T, Beach TG, Grover A, Mastroeni D, Walker DG, Caselli RJ, Kukull WA, McKeel D, Morris JC: Gene expression profiles in anatomically and functionally distinct regions of the normal aged human brain. Physiol Genomics. 2007, 28 (3): 311-322.
Hodges A, Strand AD, Aragaki AK, Kuhn A, Sengstag T, Hughes G, Elliston LA, Hartog C, Goldstein DR, Thu D: Regional and cellular gene expression changes in human Huntington's disease brain. Hum Mol Genet. 2006, 15 (6): 965-977. 10.1093/hmg/ddl013.
Moran LB, Duke DC, Deprez M, Dexter DT, Pearce RK, Graeber MB: Whole genome expression profiling of the medial and lateral substantia nigra in Parkinson's disease. Neurogenetics. 2006, 7 (1): 1-11. 10.1007/s10048-005-0020-2.
Tyrrell J, Cosgrave M, McCarron M, McPherson J, Calvert J, Kelly A, McLaughlin M, Gill M, Lawlor BA: Dementia in people with Down's syndrome. Int J Geriatr Psychiatry. 2001, 16 (12): 1168-1174. 10.1002/gps.502.
Zigman WB, Schupf N, Sersen E, Silverman W: Prevalence of dementia in adults with and without Down syndrome. Am J Ment Retard. 1996, 100 (4): 403-412.
Lockstone HE, Harris LW, Swatton JE, Wayland MT, Holland AJ, Bahn S: Gene expression profiling in the adult Down syndrome brain. Genomics. 2007, 90 (6): 647-660. 10.1016/j.ygeno.2007.08.005.
Ryan MM, Lockstone HE, Huffaker SJ, Wayland MT, Webster MJ, Bahn S: Gene expression analysis of bipolar disorder reveals downregulation of the ubiquitin cycle and alterations in synaptic genes. Mol Psychiatry. 2006, 11 (10): 965-978. 10.1038/sj.mp.4001875.
Jucker M: The benefits and limitations of animal models for translational research in neurodegenerative diseases. Nat Med. 16 (11): 1210-1214.
Aimone JB, Leasure JL, Perreau VM, Thallmair M: Spatial and temporal gene expression profiling of the contused rat spinal cord. Exp Neurol. 2004, 189 (2): 204-221. 10.1016/j.expneurol.2004.05.042.
Gurney ME, Pu H, Chiu AY, Dal Canto MC, Polchow CY, Alexander DD, Caliendo J, Hentati A, Kwon YW, Deng HX: Motor neuron degeneration in mice that express a human Cu, Zn superoxide dismutase mutation. Science. 1994, 264 (5166): 1772-1775. 10.1126/science.8209258.
Wong PC, Pardo CA, Borchelt DR, Lee MK, Copeland NG, Jenkins NA, Sisodia SS, Cleveland DW, Price DL: An adverse property of a familial ALS-linked SOD1 mutation causes motor neuron disease characterized by vacuolar degeneration of mitochondria. Neuron. 1995, 14 (6): 1105-1116. 10.1016/0896-6273(95)90259-7.
Cunningham C, Deacon R, Wells H, Boche D, Waters S, Diniz CP, Scott H, Rawlins JN, Perry VH: Synaptic changes characterize early behavioural signs in the ME7 model of murine prion disease. Eur J Neurosci. 2003, 17 (10): 2147-2155. 10.1046/j.1460-9568.2003.02662.x.
Tariot PN, Solomon PR, Morris JC, Kershaw P, Lilienfeld S, Ding C: A 5-month, randomized, placebo-controlled trial of galantamine in AD. The Galantamine USA-10 Study Group. Neurology. 2000, 54 (12): 2269-2276.
Wilcock GK, Lilienfeld S, Gaens E: Efficacy and safety of galantamine in patients with mild to moderate Alzheimer's disease: multicentre randomised controlled trial. Galantamine International-1 Study Group. Bmj. 2000, 321 (7274): 1445-1449. 10.1136/bmj.321.7274.1445.
Rockwood K, Mintzer J, Truyen L, Wessel T, Wilkinson D: Effects of a flexible galantamine dose in Alzheimer's disease: a randomised, controlled trial. J Neurol Neurosurg Psychiatry. 2001, 71 (5): 589-595. 10.1136/jnnp.71.5.589.
Prvulovic D, Hampel H, Pantel J: Galantamine for Alzheimer's disease. Expert Opin Drug Metab Toxicol. 6 (3): 345-354.
Ago Y, Koda K, Takuma K, Matsuda T: Pharmacological aspects of the acetylcholinesterase inhibitor galantamine. J Pharmacol Sci. 116 (1): 6-17.
Kang SS, Lee JY, Choi YK, Kim GS, Han BH: Neuroprotective effects of flavones on hydrogen peroxide-induced apoptosis in SH-SY5Y neuroblostoma cells. Bioorg Med Chem Lett. 2004, 14 (9): 2261-2264. 10.1016/j.bmcl.2004.02.003.
Liu R, Zhang T, Yang H, Lan X, Ying J, Du G: The flavonoid apigenin protects brain neurovascular coupling against amyloid-beta-induced toxicity in mice. J Alzheimers Dis. 24 (1): 85-100.
Chen HQ, Jin ZY, Wang XJ, Xu XM, Deng L, Zhao JW: Luteolin protects dopaminergic neurons from inflammation-induced injury through inhibition of microglial activation. Neurosci Lett. 2008, 448 (2): 175-179. 10.1016/j.neulet.2008.10.046.
Dragicevic N, Smith A, Lin X, Yuan F, Copes N, Delic V, Tan J, Cao C, Shytle RD, Bradshaw PC: Green Tea Epigallocatechin-3-Gallate (EGCG) and Other Flavonoids Reduce Alzheimer's Amyloid-Induced Mitochondrial Dysfunction. J Alzheimers Dis. 2011, 26 (3): 507-21.
Boniece IR, Wagner JA: Growth factors protect PC12 cells against ischemia by a mechanism that is independent of PKA, PKC, and protein synthesis. J Neurosci. 1993, 13 (10): 4220-4228.
Johnson K, Liu L, Majdzadeh N, Chavez C, Chin PC, Morrison B, Wang L, Park J, Chugh P, Chen HM: Inhibition of neuronal apoptosis by the cyclin-dependent kinase inhibitor GW8510: identification of 3' substituted indolones as a scaffold for the development of neuroprotective drugs. J Neurochem. 2005, 93 (3): 538-548. 10.1111/j.1471-4159.2004.03004.x.
Nelson DL, Herbet A, Petillot Y, Pichat L, Glowinski J, Hamon M: [3H]Harmaline as a specific ligand of MAO A--I. Properties of the active site of MAO A from rat and bovine brains. J Neurochem. 1979, 32 (6): 1817-1827. 10.1111/j.1471-4159.1979.tb02296.x.
Herraiz T, Gonzalez D, Ancin-Azpilicueta C, Aran VJ, Guillen H: beta-Carboline alkaloids in Peganum harmala and inhibition of human monoamine oxidase (MAO). Food Chem Toxicol. 48 (3): 839-845.
Frost D, Meechoovet B, Wang T, Gately S, Giorgetti M, Shcherbakova I, Dunckley T: beta-Carboline Compounds, Including Harmine, Inhibit DYRK1A and Tau Phosphorylation at Multiple Alzheimer's Disease-Related Sites. PLoS One. 6 (5): e19264-
Sitz JH, Baumgartel K, Hammerle B, Papadopoulos C, Hekerman P, Tejedor FJ, Becker W, Lutz B: The Down syndrome candidate dual-specificity tyrosine phosphorylation-regulated kinase 1A phosphorylates the neurodegeneration-related septin 4. Neuroscience. 2008, 157 (3): 596-605. 10.1016/j.neuroscience.2008.09.034.
Coccaro EF, Siever LJ: Second generation antidepressants: a comparative review. J Clin Pharmacol. 1985, 25 (4): 241-260.
Melamed E, Rosenthal J, Globus M, Cohen O, Uzzan A: Suppression of MPTP-induced dopaminergic neurotoxicity in mice by nomifensine and L-DOPA. Brain Res. 1985, 342 (2): 401-404. 10.1016/0006-8993(85)91146-1.
McKinley ET, Baranowski TC, Blavo DO, Cato C, Doan TN, Rubinstein AL: Neuroprotection of MPTP-induced toxicity in zebrafish dopaminergic neurons. Brain Res Mol Brain Res. 2005, 141 (2): 128-137.
Paris I, Martinez-Alvarado P, Perez-Pastene C, Vieira MN, Olea-Azar C, Raisman-Vozari R, Cardenas S, Graumann R, Caviedes P, Segura-Aguilar J: Monoamine transporter inhibitors and norepinephrine reduce dopamine-dependent iron toxicity in cells derived from the substantia nigra. J Neurochem. 2005, 92 (5): 1021-1032. 10.1111/j.1471-4159.2004.02931.x.
Bedard P, Parkes JD, Marsden CD: Nomifensine in Parkinson's disease. Br J Clin Pharmacol. 1977, 4 (Suppl 2): 187S-190S.
Hardman JG, Limbird LE, Gilman AG: Goodman & Gilman's The Pharmacological Basis of Therapeutics. 2001, New York: McGraw-Hill, 10
This work was supported by the Wolfson CARD. I would like to thank Owen Williams for help with the leukaemia analysis, Ramin Raouf and Patrick Doherty for general discussions on the database strategy.
GW is the sole author of the present study.
Electronic supplementary material
Additional file 1: Expression profiles derived from CMAP and published AD studies. Expression profile signatures derived from pooling CMAP compound treatment replicates, form corticosteroid resistance studies and from various AD stages referenced to age matched controls. These profiles are used to query SPIED. The data is given as excel spreadsheets. (XLS 2 MB)
Additional file 2: Top scoring correlations for various expression signature queries against SPIED. The results upon querying SPIED with expression profiles derived from the signatures in additional file 1. The data is given as excel spreadsheets. (XLS 178 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Williams, G. A searchable cross-platform gene expression database reveals connections between drug treatments and disease. BMC Genomics 13, 12 (2012). https://doi.org/10.1186/1471-2164-13-12
- Global gene expression
- connectivity map
- neurodegenerative disease
- Alzheimer's disease