Predicting environmental chemical factors associated with disease-related gene expression data
© Patel and Butte. 2010
Received: 23 October 2009
Accepted: 6 May 2010
Published: 6 May 2010
Skip to main content
© Patel and Butte. 2010
Received: 23 October 2009
Accepted: 6 May 2010
Published: 6 May 2010
Many common diseases arise from an interaction between environmental and genetic factors. Our knowledge regarding environment and gene interactions is growing, but frameworks to build an association between gene-environment interactions and disease using preexisting, publicly available data has been lacking. Integrating freely-available environment-gene interaction and disease phenotype data would allow hypothesis generation for potential environmental associations to disease.
We integrated publicly available disease-specific gene expression microarray data and curated chemical-gene interaction data to systematically predict environmental chemicals associated with disease. We derived chemical-gene signatures for 1,338 chemical/environmental chemicals from the Comparative Toxicogenomics Database (CTD). We associated these chemical-gene signatures with differentially expressed genes from datasets found in the Gene Expression Omnibus (GEO) through an enrichment test.
We were able to verify our analytic method by accurately identifying chemicals applied to samples and cell lines. Furthermore, we were able to predict known and novel environmental associations with prostate, lung, and breast cancers, such as estradiol and bisphenol A.
We have developed a scalable and statistical method to identify possible environmental associations with disease using publicly available data and have validated some of the associations in the literature.
The etiology of many diseases results from interactions between environmental factors and biological factors . Our knowledge regarding interaction between environmental factors, such chemical exposure, and biological factors, such as genes and their products, is increasing with the advent of high-throughput measurement modalities. Building associations between environmental and genetic factors and disease is essential in understanding pathogenesis and creating hypotheses regarding disease etiology. However, it is currently difficult to ascertain multiple associations of chemicals to genes and disease without significant experimental investment or large-scale epidemiological study. Use of publicly-available environmental chemical factor and genomic data may facilitate the discovery of these associations.
We desired to use pre-existing datasets and knowledge-bases in order to derive hypotheses regarding chemical association to disease without upfront experimental design. Specifically, we asked what environmental chemicals could be associated with gene expression data of disease states such as cancer, and what analytic methods and data are required to query for such correlations. This study describes a method for answering these questions. We integrated publicly available data from gene expression studies of cancer and toxicology experiments to examine disease/environment associations. Central to our investigation was the Comparative Toxicogenomics Database (CTD) , which contains information about chemical/gene/protein interactions and chemical/gene/disease relationships, and the Gene Expression Omnibus (GEO) , the largest public gene expression data repository. Information in the CTD is curated from the peer-reviewed literature, while gene expression data in GEO is uploaded by submitters of manuscripts.
have tested a small number of chemicals on cells and measured responses on a genomic scale, or
used existing knowledge bases, such as Gene Ontology, to associate annotated pathways to environmental insult
The first method involves measuring physiological response on a gene expression microarray. This approach allows researchers to test chemical association on a genomic scale, but the breadth of discoveries is constrained by the number of chemicals tested against a cell line or model organism. These experiments are not intended for hypothesis generation across hundreds of potential chemical factors with multiple phenotypic states. Only a few chemicals can be tractably tested for association to gene activity [4, 5], or disease on cell lines , or on model organisms, including rat and mouse . In rare cases, this approach has reached the level of a hundred or thousand chemical compounds, such as the Connectivity Map, developed by Lamb, Golub, and colleagues , which attempts to associate drugs with gene expression changes. After measuring the genome-wide effect on gene expression after application of hundreds of drugs at various doses, drug signatures are calculated and are then queried with other datasets for which a potential therapeutic is desired. While this has proven to be an excellent system to find chemicals that essentially reverse the genome-wide effects seen in disease, the approach of measuring gene expression and calculating signatures across tens of thousands of environmental chemicals is not always feasible or scalable. Although other data-driven approaches have been described , few have given insight into external causes of disease.
A second approach has been to use knowledge bases, such as Gene Ontology  to aid in the interpretation of genomic results. For example, Gene Ontology analysis of a cancer experiment might elucidate a molecular mechanism related to an environmental chemical. Unfortunately, there is still a lack of methodology to derive hypotheses for environmental-genetic associations in disease pathogenesis, as Gene Ontology and general gene-set based approaches have limited information on environmental chemicals.
In contrast to the previous approaches, we claim that the integration of pre-existing data and knowledge bases can derive hypotheses regarding the association of chemicals to gene activity and disease from multiple datasets in a scalable manner. Gohlke et al have proposed an approach to predict environmental chemicals associated with phenotypes also using knowledge from the CTD . Their method utilizes the Genetic Association Database (GAD)  to associate phenotypes to genetic pathways and the CTD to link pathways to environmental factors. This method has proved its utility, allowing for production of hypotheses for chemicals associated with diseases categorized as metabolic or neuropsychiatric disorders. However, in its current configuration, their method is dependent on the GAD, which contains statically annotated phenotypes in relation to genes containing variants; such DNA changes are not likely to be reflective of molecular profiles of tissues being suspected for environmental influence. Unlike this method, our proposed approach is tissue- and data-driven in that the phenotype is determined by the individual measurements of gene expression in cells and tissues, allowing for the dynamic capture of phenotypes.
The approach we propose here is agnostic to experiment protocol, such as cell line or chemical agent tested, and provides for a less resource-intensive screening of chemicals to biologically validate. Our methodology essentially combines the best features of these current approaches. We start by compiling "chemical signatures" in a scalable way using the CTD. These chemical signatures capture known changes in gene expression secondary to hundreds of environmental chemicals. In a manner similar to how Gene Ontology categories are tested for over-representation, we then calculate the genes differentially expressed in disease-related experiments and determine which chemical signatures are significantly over-represented. We first verified the accuracy of our methodology by analyzing microarray data of samples with known chemical exposure. After these verification studies yielded positive results, we then applied the method to predict disease-chemical associations in breast, lung, and prostate cancer datasets. We validated some of these predictions with curated disease-chemical relations, warranting further study regarding pathogenesis and biological mechanism in context of environmental exposure. Our method appears to be a promising and scalable way to use existing datasets to predict environmental associations between genes and disease.
With the single gene, single chemical relationships, we created "chemical signatures", or gene sets associated with each chemical (Figure 1B). Gene sets were created from gene-expression relations spanning 249 species, but most relations came from H. sapiens, M. musculus, R. norvegicus, and D. rerio. We eliminated chemical-gene sets that had less than 5 genes in the set. This step yielded a total of 1,338 chemical-gene sets.
The CTD also contains curated data regarding the association of a diseases to chemicals. These associations are either shown in an experimental model physiological system or through epidemiological studies. We used these curated associations to validate our predicted factors associated to disease. There are 3,997 diseases-chemical associations in the CTD, consisting of 653 diseases (annotated by unique MeSH terms) and 1,515 chemicals (Figure 1C). The median, 70th, and 75th , and 80th percentile of the number of curated chemicals per disease is 2, 3, 4, and 5 respectively.
We used Significance Analysis of Microarrays (SAM) software to select differentially expressed genes from a microarray experiment . The FDR for SAM for all of our predictions were controlled up to a maximum of 5 to 7% in order to reduce false associations.
We mapped microarray annotations to other corresponding representative species, H. sapiens, M. musculus, and R. norvegicus using Homologene . In the CTD, gene identifiers were commonly associated with H. sapiens; however, some are mapped to specific organisms, such as M. musculus and R norvegicus. Most mappings in the CTD are among these 3 organisms. By mapping our expression annotation to these organisms, we ensured gene compatibility with a large portion of the CTD.
We checked for enrichment of differentially expressed genes among our 1,338 chemical-gene sets with the hypergeometric test. To account for multiple hypothesis testing, we computed the q-value, or false discovery rate for a given p-value, by using 100 random resamplings of genes from the microarray experiment and testing each of these random resamplings for enrichment against each of the 1,338 chemical-gene sets. This methodology is similar to the q-value estimation method described in "GoMiner", a gene ontology enrichment assessment tool . We assessed a positive prediction for those that had exceeded a certain p-value and q-value threshold in our list of 1,338 tested associations. All analyses were conducted using the R statistical environment .
For our verification phase, we surveyed publicly available data from the Gene Expression Omnibus (GEO) for experiments in which sets of samples exposed to chemicals were compared with controls. We found and used six datasets in the validation phase. Set 1 included GSE5145 (3 study samples and 3 controls) in which H. sapiens muscle cell samples were exposed to Vitamin D . Set 2 was GSE10082 (6 study samples and 5 controls) in which wild-type M. musculus were exposed to tetradibenzodioxin (TCDD) . Set 3 was GSE17624 in which H. sapiens Ishikawa cells (4 study samples and 4 controls) were exposed to high doses of bisphenol A (no reference). Set 4 was GSE2111 in which H. sapiens bronchial tissue (4 study samples and 4 controls) were exposed to zinc sulfate . The CTD had some chemical-gene relations based on this dataset; we removed these relations prior to computing the predictions for this dataset. Set 5 was GSE2889 in which M. musculus thymus tissues (2 study samples and 2 controls) were exposed to estradiol . Finally, set 6 was GSE11352 in which H. sapiens MCF-7 cell line was exposed to estradiol at 3 different time points . In all cases except for set 6, we treated SAM analysis as unpaired t-tests; for set 6, we used the time-course option in SAM. See Additional File 1 for the number of differentially expressed genes found for each dataset along with their median false discovery rate (Additional file 1, Supplementary Table S1).
We found previously measured cancer gene expression datasets to identify potential environmental associations with cancer. We used measurements from human prostate cancer from GSE6919 [23, 24], lung cancer from GSE10072 , and breast cancer from GSE6883 . We conducted all SAM analyses using an unpaired t-test between disease and control samples. Additional File 1 shows the number of differentially expressed genes measured for each dataset along with the level of FDR control (Additional file 1, Supplementary Table S2).
We deliberately chose cancer datasets that used a different population of controls rather than normal tissues from the same patients. The prostate cancer dataset (GSE6919) consisted of 65 prostate tissue cancer samples and 17 normal prostate tissue samples as controls.
The lung cancer dataset (GSE10072) consisted of two patient groups: non-smokers with cancer (historically and currently), and current smokers with cancer. We conducted the predictions on these groups separately. The cancer-non smoker group consisted of 16 samples and the cancer-smoker group had 24 samples. The control group consisted of 15 samples.
The breast cancer dataset (GSE6883) consisted of two distinct cancer sub-groups: non-tumorigenic and tumorigenic. As with the lung cancer data, we conducted our predictions on these groups separately. The non-tumorigenic group consisted of three samples and the tumorigenic group had six samples. The control group contained three samples.
We then validated our highly ranked factor predictions with disease-chemical knowledge from the CTD. In particular, we determined if the highly significant chemicals in our prediction list included those that had curated relationship with cancer in the CTD (disease-chemicak relation). This step was similar to measuring association to chemicals via enriched gene sets using the hypergeometric test as described above. We used curated factors associated with Prostatic Neoplasms (MeSH ID: D011471), Lung Neoplasms (D008175), and Breast Neoplasms (D001943), to validate our predictions generated with the prostate cancer, lung cancer, and breast cancer datasets respectively. Further, we assessed the validation by computing the actual number of false positives and true negatives. To compute this number, we assessed whether the prediction list was enriched for chemicals associated with any of the other diseases in the CTD at a higher significance level than the true disease; for this test, we chose diseases that had at least 5 chemical associations, a total of 141 diseases. As an example, to assess the false positive rate for the prostate cancer (MeSH ID: D011471) predictions, we determined the curated enrichment of our predictions for all 140 other disease-chemical sets and counted the number of diseases that had a lower p-value than that computed for D011471.
Chemical-gene sets derived from the CTD are but one representation of how a chemical might affect biological activity. Biological activity of chemicals may also be derived from high-throughput, in-vitro chemical screens such as those archived in PubChem [27, 28]. Specifically, the PubChem database provides a large number of phenotypic measurements (or "BioAssays") for many of the chemicals we predicted for cancer. In addition, PubChem provides tools to compare BioAssay measurements for different chemicals. Quantitative and standardized BioAssay measurements (normalized "scores") allow comparison of biological activities of chemicals and derivation of biological activity similarity between chemicals. For example, PubChem represents the biological activity of a compound through a vector of BioAssay scores and assembles a bioactivity similarity matrix between each pair of chemicals with this data.
We implemented a method to predict a list of environmental factors associated with differentially expressed genes (Figure 2). The method is centered on chemical-gene sets that are derived from single curated chemical-gene relationships in the CTD. We determine whether the differentially expressed genes are associated to a chemical by assessing if the expressed genes are enriched for a chemical-gene set, or contain more genes from the chemical-gene set than expected at random using the hypergeometric test. We applied this method in two phases, the first a verification phase in which we sought to rediscover known exposures applied to samples, and a query phase, in which we sought to find factors associated with cancer gene expression datasets. We refer to significant chemical-gene set associations to gene expression data as "associations" or "predictions" in the following.
We first applied our method to gene expression data from experiments in which samples were exposed to specific chemicals, reasoning that if our method could identify these known chemical exposures, we could use the method to predict chemicals that may have perturbed gene expression in unknown experimental or disease conditions. Our goal was to determine where a gene expression-altering chemical might lie in the range of significance rankings applied by the prediction method.
Chemical Prediction Results from the Verification Phase.
Actual Chemical Exposure (GEO accession)
Relevant Genes Expressed
Vitamin D3 on H. sapiens muscle cells (GSE5145)
1 × 10-23
VDR (25), CYP24A1 (14)
TCDD on M. musculus (GSE10082)
2 × 10-15
CYP1A1 (59), CYP1B1 (15), AHRR(6), CYP1A2 (14)
Bisphenol A on H. sapiens Ishikawa cells (GSE17624)
1 × 10-6
ESR1(31), ESR2(7), S100G (6)
Zinc sulfate on H. sapiens bronchial tissue (GSE2111)
3 × 10-3
SLC30A1 (3), MT1F(2), MT1G(2)
Estradiol on M. musculus thymus (GSE2889)
5 × 10-3
C3(6), LPL (4), CTSB (2)
Estradiol on H. sapiens MCF7 cell line (GSE11352)
5 × 10-3
ISG20 (2), MGP (2), SERPINA1 (2)
We were able to satisfactorily predict the exposures applied to the gene expression datasets. We ascertained a positive prediction if the exposure had a relatively high ranking (low p-value for enrichment) and if the q-value was lower than 0.1. For the datasets measuring expression after exposure to Vitamin D, calcitriol, a type of vitamin D, was ranked first in the list (p = 10-23, q = 0). Similarly, TCDD was predicted third in its respective list (p = 10-15, q = 0). The other exposures ranked within the top percentile, ranging from 15 to 19; the lower bound of p-values were between 10-6 and 0.01 and q-values less than 0.1. We reasoned that we could detect true associations between environmental chemicals and gene expression phenotypes provided they met these significance thresholds.
We applied our prediction methods to datasets measuring the gene expression for prostate, breast, and lung cancers. In particular, we computed predictions for prostate cancer from primary prostate tumor tissue, lung adenocarcinomas from lung tissue from non-smoking individuals, and non-tumorigenic breast cancer cells grown in mouse xenografts. Additional File 1 shows predictions for related data on tumorigenic breast cancer and smoker lung cancer samples (Additional file 1, Supplementary Tables S3 and S4). To validate and select specific predictions from our ranked list of 1,338 environmental chemicals, we measured how enriched top-ranking chemicals were for annotated disease-chemical citations in for diseases of interest ("Prostate Neoplasms", "Breast Neoplasms", and "Lung Neoplasms"). To call a positive chemical association or prediction to disease phenotype, we used p-value thresholds similar to what we observed during the verification phase (α ≤ 10-4, 0.001, 0.01) along with q-values as low as possible, specifically less than 0.1. For comparison, we also used the typical p-value threshold of 0.05.
Prediction of environmental chemicals associated with prostate cancer samples (GSE6919).
Relevant genes in set (number of references)
4 × 10-10
ESR2(37), PGR(34), MAPK1(14)
1 × 10-9
ESR2(6), IGF1(5), BCL2(4)
1 × 10-8
JUN(13), MAPK1(9), CCND1(8), FOS(6)
7 × 10-7
BCL2(23), MAPK1(14), TNF(10)
6 × 10-6
MT2A(14), MT1A(12), MT3(11), MT1(6)
3 × 10-5
6 × 10-4
ESR2(22), PGR (10), MAPK1 (5)
3 × 10-5
ESR2(8), FOS(8), HOXA10(4)
3 × 10-4
BCL2(3), ELF3(2), LDHA(2)
6 × 10-4
PGR(8), ESR2(7), IL4RA(2)
9 × 10-4
MT3(18), MT2A(13), MT1A(11)
Prediction of environmental chemicals associated with lung cancer samples (GSE10072).
Relevant genes in set (number of references)
1 × 10-6
4 × 10-4
CASP3(60), ABCB1(28), BAX(26), BCL2 (23)
8 × 10-6
4 × 10-4
JUN(13), NQ01(6), EGR1(6)
1 × 10-5
6 × 10-4
HBEGF(3), CDK7(1), CDKN1B (1), CDKN1C(1)
6 × 10-5
7 × 10-4
TGFB1(23), TIMP1(15), PCNA(6)
2 × 10-4
BIRC5(3), CDKN1B(2), MMP9(2)
3 × 10-4
ABCB1(4), ABCG2(4), KRT19(2)
4 × 10-4
IL6(2), MMP9(2), MMP12(2), PDGFB(2)
Prediction of environmental chemicals associated with breast cancer samples (GSE6883).
Relevant genes in set (number of references)
2 × 10-4
IL6(3), STC1(3), CEBPD(2)
6 × 10-4
CEBPD(1), APLP2(1), MLF1(1)
7 × 10-4
LPL(4), IL6(3), CEBPD(2)
3 × 10-3
CCDC50(1), BIRC3(1), DNAJB(1)
3 × 10-3
IL6(1), MARCKS(1), MXD1(1), MMP7(1)
4 × 10-4
IL6(3), MARCKS(1), PSMA5(1)
6 × 10-3
CEBPD(1), MLF1(1), DTL(1)
We have described a method of generating a list of chemical predictions associated with disease-annotated gene expression datasets and applied the method on gene expression data for several cancers. We have validated a subset of our predictions with evidence from the literature as described above (Tables 2, 3, 4).
We sought further evidence of the biological relevance of our predictions through internal comparison of their potential activity archived in PubChem. Specifically, we expected some degree of correlation between "similar" chemicals and their gene set significance to the cancer datasets. We opted to use PubChem BioActivity to assess chemical similarity, assuming this measure of phenotypic similarity would be representative of underlying biological pathways of action. We picked chemicals that were deemed significant for thresholds used above (p = 0.001, 0.001, 0.01, for the prostate, lung, and breast cancer datasets) for all of the cancer datasets. This resulted in a total of 130 chemicals, 66 of which had BioActivity data in PubChem. The BioActivity similarity for each of the 66 chemicals was computed through 790 BioAssay scores. Figure 5 shows the -log10 of significance for the highest ranked chemical predictions clustered by their BioActivity similarity.
We found some chemicals with similar biological activity profiles in PubChem had similar patterns of chemical-gene set association across the cancer datasets. For example, sodium arsenite, sodium arsenate, and doxorubicin have closely related biological profiles as well as high significance of chemical-gene set association for the prostate and lung cancer data (Figure 5, enclosed in orange box); however, we did not observe other biologically similar chemicals such as Tetradihydrobenzodioxin. On the other hand, we also observed correlation between the biological activity similarity and chemical-gene set association for hormone or steroidal chemicals such as ethinyl estradiol, estradiol, and diethylstilbestrol as well as progesterone and corticosterone (Figure 5, enclosed in purple boxes).
We have developed a knowledge- and data-driven method to predict chemical associations with gene expression datasets, using publicly available and previously disjoint datasets. To our knowledge, there are few methods that generate hypotheses regarding environmental associations with disease from gene expression data. Most current approaches in toxicology have focused on a small number of environmental influences on single or small groups of genes, while current approaches in toxicogenomics have been concentrated on measuring genome-wide responses for a few chemicals . Our prediction method enables the generation of hypotheses in a larger scalable manner using existing data, examining the potential role of hundreds of chemicals over thousands of genome-wide measurements and diseases.
As an example, we found predicted chemicals such as sodium arsenite in its association with prostate and lung cancers, estrogenic compounds such as bisphenol A and estradiol with prostate and breast cancers, and dimethylnitrosamine with lung cancer. Although each has curated knowledge behind the association in the CTD, mechanisms for the action are not well known and call for further study. So far, Benbrahim-Talaa et al have found hypomethylation patterns in the presence of arsenic in prostate cancer cells . Zanesi et al show a potential interaction role of FHIT gene and dimethylnitrosamine to produce lung cancers . Evidence of a complex mechanistic action of estrogens, such as estradiol, on breast cancer carcinogenesis has been established ; however the role of other estrogenic-like compounds have only recently been studied. For example, bisphenol A has been shown to invoke an aggressive response in cancer cell lines , possibly by affecting estrogen-dependent pathways . It is evident that more experimentation is required involving the measurements of exposure-affected proteins and genes and their activation state in cellular models and their relation to the chemical signatures.
An overlap of activity of the same genes induced by different chemicals would suggest a common physiological action by these chemicals. For example, the ESR2 and MAPK1 genes in the prostate cancer prediction, and the IL6 and CEBPD in the breast cancer predictions, were associated with several chemicals for each of the diseases. We also found an overlap between chemicals amongst different cancers. This result comes as a result of the correlation in the significant pathways shared by these cancers; however, it may also indicate a need to explore less significant associations in order to find unique and specific gene expression/chemical exposure relationships for a given disease. Furthermore, this result may also indicate a bias of gene and chemical relationships cataloged in the CTD. For example, it could be that genes specific to common cancer-related pathways are those that are well studied, such as BCL2 or ESR2.
Related to this, we have attempted to show how biological activity, as assayed in a high-throughput chemical screen in PubChem, can be correlated with chemical gene-set associations. Observing a correlation in both PubChem-derived bioactivity in addition to a chemical-gene set association from the CTD provides a way to identify shared modes of action among groups of similar or related chemicals. This data serves to both provide internal validation for list of predicted chemicals acting through similar pathways (such as those induced by estrogen) but also to prioritize hypotheses. For example, we did not find curated evidence in the CTD for association of the chemicals vinclozolin, tert-Butylhydroperoxide, and Carbon Tetrachloride to prostate or lung cancers; however, their similar bioactivity profiles (Figure 5, enclosed in blue box) and high chemical-gene set association calls for further review.
We do acknowledge some arbitrariness in our choice of methods and thresholds; most of these were chosen to show significance in our methodology without adding complexity. We could have chosen any of several alternative approaches to implementing our method; however, predictions made with the Gene Set Enrichment Analysis (GSEA)  method during the verification phase were not as sensitive (not shown). Another limitation in our first implementation is that in calculating the chemical signatures associating chemicals with gene sets, we ignored the specific degree of expression change (up or down) encoded in the CTD. We decided not to use this information due to the presence of contradictions (some references may point to an increase of exposure-induced gene expression while another reference might claim the opposite), and other preliminary work suggesting that filtering by the degree of change reduced sensitivity (data not shown). Because of these limitations, direction of association cannot be inferred. Further still, we acknowledge richer and more refined chemical signatures along with further integration with resources like PubChem will need to be built to make the most accurate predictions.
Another issue with querying the microarray data of any experiment is the lack of full sample information to stratify results; for example, different exposures may be associated with a subset of the samples. A related concern includes small sample sizes of some of the datasets used to evaluate the method. For example, the best predictive power was seen the largest dataset (prostate cancer, GSE6919), and the worst with one of the smallest, (breast cancer, GSE6883). Despite this heterogeneity and lack of power, we still arrived at noteworthy and literature-backed findings warranting further study. We also urge that more evaluation must occur with datasets that have a larger number of samples.
Most importantly, we stress that these types of association remain as predictions and hypotheses that need validation and verification. The method presented here is not a substitute for traditional toxicology or epidemiology. These studies are required to provide quantitative and population generalizable estimates of disease risk and dose-response relationships. However, as the space of potential environmental chemicals potentially causing biological effects is large, we suggest that this methodology would give investigators at least some clue where to start the search for environmental causal factors to study in these other modes. Furthermore, predicting a linkage between chemicals, genes, and clinically-relevant disease phenotypes using existing resources falls in agreement with the National Academies' vision of high-throughput efforts to decipher genetic pathways to toxicity .
We have described a novel and scalable method to associate changes in gene expression with environmental chemicals. While we successfully validated our methodology here and provide hypotheses regarding the potential association of chemicals in cancer development, these hypotheses would need to be carefully studied in controlled cellular experiments. Our method is limited by the lack of direction of association and effect size as typically ascertained in traditional toxicological and epidemiological studies; however, the vast number of chemicals that can be tested in silico is only limited by the amount of available data. This method is just one of potentially many tools that need to be built to predict environmental associations between genes and disease.
Comparative Toxicogenomics Database
Gene expression omnibus
Genetic Association Database
Gene Set Enrichment Analysis
Medical Subject Headings
Significance Analysis of Microarrays
False Discovery Rate
CJP was funded by National Library of Medicine (T15 LM 007033). AJB was funded by the Lucile Packard Foundation for Children's Health, the National Library of Medicine (R01 LM009719), the National Institute of General Medical Sciences (R01 GM079719), and the Howard Hughes Medical Institute. We thank Alex Skrenchuk and Boris Oskotsky from Stanford University for computer support and Rong Chen from Stanford University for critical review.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.