Skip to main content

Advertisement

Genetic and environmental pathways to complex diseases

BMC Systems Biology volume 3, Article number: 46 (2009) Cite this article

Article metrics

Abstract

Background

Pathogenesis of complex diseases involves the integration of genetic and environmental factors over time, making it particularly difficult to tease apart relationships between phenotype, genotype, and environmental factors using traditional experimental approaches.

Results

Using gene-centered databases, we have developed a network of complex diseases and environmental factors through the identification of key molecular pathways associated with both genetic and environmental contributions. Comparison with known chemical disease relationships and analysis of transcriptional regulation from gene expression datasets for several environmental factors and phenotypes clustered in a metabolic syndrome and neuropsychiatric subnetwork supports our network hypotheses. This analysis identifies natural and synthetic retinoids, antipsychotic medications, Omega 3 fatty acids, and pyrethroid pesticides as potential environmental modulators of metabolic syndrome phenotypes through PPAR and adipocytokine signaling and organophosphate pesticides as potential environmental modulators of neuropsychiatric phenotypes.

Conclusion

Identification of key regulatory pathways that integrate genetic and environmental modulators define disease associated targets that will allow for efficient screening of large numbers of environmental factors, screening that could set priorities for further research and guide public health decisions.

Background

Determining the extent to which environmental versus genetic factors are responsible for particular phenotypes is a central question in all of biological research. Elucidating associations between genotype and phenotype has been a central goal in human health research for some time, and has resulted in an impressive collection of research on genotype-phenotype relationships [1, 2]. While continued analysis of rare monogenic phenotypes is important for mechanistic discoveries [3], unraveling the interplay between genetic and environmental determinants of complex phenotypes will be critical for improving public health [4]. For example, gene-environment interactions have been shown to play a critical role in childhood leukemia and asthma [5–7]. However, much less is known about gene-environment interactions as they relate to the etiology of the common complex disease phenotypes such as unipolar depressive disorders, ischemic heart disease and cerebrovascular disease, all of which fall within the top six causes of the global burden of disease, and are projected to increase as the epidemiological transition continues in developing countries [8].

Network and bioinformatic methods have recently been applied to synthesize data on gene-disease relationships for those diseases that have a strong genetic component [9–11]. In addition, utilization of functional information to prioritize candidate driver genes in cancer has been advocated [12]. However, application of network theory to determine the interplay between genetics and environmental factors in complex diseases has been left unexplored. We hypothesize genetic and environmental factors involved in the progression of a particular complex phenotype are participants in the same underlying cellular processes. To test this hypothesis, we develop networks of complex diseases and environmental factors through linkage of human genetic association studies and mechanistic analyses of environmental factors, using evolutionarily conserved molecular pathways as the unifying system to define relationships. We further explore relationships identified by this method through comparison to known disease-chemical relationships and analysis of transcriptional regulation in gene expression datasets for metabolic syndrome phenotypes, neuropsychiatric phenotypes and several predicted environmental modulators.

Results and discussion

Clustering phenotypes by pathways

To identify common pathways between complex diseases, we annotated gene-phenotype relationships found in the Genetic Association Database (GAD) [1] (see Additional File 1), then analyzed these phenotype-associated gene lists using the Structurally Enhanced Pathway Enrichment Analysis (SEPEA) algorithm (See methods for summary and described in detail in [13]). This resulted in a clustergram of disease phenotypes based on the underlying pathways that are represented by the sum of polymorphic genes associated with a particular phenotype (See Methods) (Figure 1). Distinct clusters of phenotypes with similar broad clinical characteristics are evident such as cancers, cardiovascular and metabolic diseases, immune-related disorders, and neuropyschiatric disorders. Furthermore, the pathways that define these clusters are consistent with our current understanding of disease etiology. For example, the cancer cluster is defined by low p-values for Erbb, p53, and cell cycle pathways and the neuropsychiatric cluster is defined by low p-values for neuroactive ligand receptor interactions, calcium signaling, as well as tryptophan and tyrosine metabolism. Moreover, immune related pathways (e.g. Jak-STAT signaling, Toll-like receptor signaling, T cell receptor signaling etc.) contribute not only to classic autoimmune and infectious disease phenotypes, but also to a large proportion of the divergent phenotypes represented, such as cardiovascular and cerebrovascular disorders, kidney disease, as well as Alzheimer's disease and longevity, among others. Several unexpected results are also evident, such as the co-clustering of pregnancy loss and preeclampsia with immune phenotypes such as lupus erythematosus and Behcet's Disease and the co-clustering of asthma with Parkinson's disease.

Figure 1
figure1

Unsupervised Hierarchical Cluster of Phenotypes by Pathways. Genes associated with a particular phenotype were evaluated for enrichment in KEGG pathways using SEPEA. P-values for KEGG pathway enrichment were then clustered using Spearman rank correlation in Cluster and the graphic was prepared using TreeView [73], where color ranges linearly from blue (p = 1) to orange (p = 0). Phenotype-gene relationships were downloaded from the Genetic Association Database [1] in June 2007 and phenotypes were further grouped according to Additional file 1. Request the TreeView file of this cluster from Julia Gohlke gohlkej@niehs.nih.gov for more detailed exploration.

Full size image

Interesting relationships are observed through a comparison of pathways that are associated with preclinical phenotypes to those pathways that are significantly associated with outright disease. For example, when we look at common neuropsychiatric disorders, such as depression and anxiety disorder, we see that genes associated with these phenotypes are specifically associated with neuroactive ligand receptor interactions, calcium signaling, as well as tryptophan and tyrosine metabolism. However, we see that these pathways significantly associated with neuropsychiatric disorders are also associated with obesity, hypertension, and blood lipoprotein composition as well as substance abuse and smoking, all of which are significant risk factors for heart disease [14]. In contrast, genes associated with outright disease phenotypes (e.g. vascular disease, heart failure, myocardial infarction, and stroke), are significantly enriched in cardiovascular specific pathways such as the renin-angiotensin system and the VEGF signaling pathway, as well as immune related pathways, suggesting genetic susceptibility to outright heart failure can be distinguished from genetic susceptibility to risk factors for development of heart disease. Therefore, this phenotype-pathway cluster of genetic associations can delineate pathways that may be important at different points in the progression of complex chronic diseases.

An interaction network of phenotypes and environmental factors

Next, we sought to meld current knowledge of genetic susceptibility factors with environmental factors that contribute to a particular complex phenotype. To accomplish this, we identified enriched pathways based on compiled lists of environmental factor-gene/protein relationships described in the Comparative Toxicogenomics Database [15]. Networks between phenotypes (using genetic association studies as described above) and environmental factors were then developed where edges represent at least 2 significant pathways between a given phenotype-phenotype, phenotype-environmental factor, or environmental factor-environmental factor pair. In addition, 8 categories including neoplastic, cardiovascular, metabolic, immune, endocrine, neuropsychiatric, pulmonary, and hematologic were used to more broadly define phenotypes. These broad categorizations are important as many of the relationships found in the Comparative Toxicogenomics Database are derived from animal models.

We compared our predicted phenotype-environmental factor relationships to a set of 1084 manually curated direct chemical-disease relationships as reported in the Comparative Toxicogenomics Database [16]. The receiver-operator curve (ROC) is illustrated in Figure 2. This figure suggests the relative loss in specificity outweighs the gain in sensitivity at a SEPEA pathway enrichment p-value cutoff of approximately 0.003 for both specific and broad categorizations of environmental factor-phenotype relationships. At this p-value cutoff, 226 of the 10,793 predicted environmental factor-phenotype relationships are supported by manually curated evidence, demonstrating the majority of connections within our network define new hypotheses of environmental factor-phenotype relationships, yet this overlap is much higher than would be expected by chance (p < 10-16). When the diseases analyzed are collapsed into the eight broad disease categories, 48% (or 271 out of 567) of the manually curated relationships are captured in our analysis (p < 10-16). This suggests that our method is more sensitive in identifying known chemical-disease category relationships than in identifying known specific disease-chemical relationships. This result makes sense in light of the fact that environmental factor data is largely derived from animal models, where one would not predict strong concordance between phenotypes and specific human diseases. In addition, the pathways analyzed are not specific to tissue and/or life stage, suggesting a specific disease of a particular tissue or developmental stage will be hard to differentiate using this method. Based on the hypothesis that there are common pathways associated with both the genetic and environmental components of broad disease categories.

Figure 2
figure2

Receiver operating characteristic (ROC) curve. A graphical representation of the sensitivity versus (1-specificity) comparing environmental factor-phenotype predictions at different p-value cutoffs to a manually curated set of direct chemical-disease relationships from the Comparative Toxicogenomics Database [16] using either specific diseases or broad categorizations of diseases. The SEPEA pathway enrichment p-value cutoff of 0.003 is indicated with arrows for each analysis.

Full size image

A graphical representation of the predicted network is presented in Figure 3, where environmental factors with known physiological actions are colored coded based on MeSH annotation (see Additional File 2 for complete annotation of nodes). Because the above comparison determined the sensitivity and specificity of our environmental factor-phenotype predictions based on only those 1084 chemical-disease relationships manually curated within CTD, we wanted to determine if the phenotype-phenotype, environmental factor-environmental factor, and environmental factor-phenotype relationships predicted in this graph are supported by known broad categorizations of phenotypes and physiological actions of the environmental factors. Therefore, we computed the significance of the number of edges that are shared between nodes in a given category using the graph clustering coefficient [17]. Using this method, the clustering of the metabolic, immune, neoplastic, and neuropsychiatric phenotypes are considered significant (p ≤ 0.05). However, when the MeSH annotated environmental factors are added, only the immune and neoplastic categories are significant (p ≤ 0.05), suggesting the broad categorization used may not be suitable to describing endocrine and cardiovascular phenotypes, or the MeSH annotated physiological actions of many of the environmental factors in this network.

Figure 3
figure3

Interaction Network of Phenotypes and Environmental Factors. Phenotypes are represented as circular nodes and environmental factors as diamond shaped nodes. Edges represent sharing at least two significantly enriched pathways (p ≤ 0.003) using lists of genes associated with a particular phenotype or environmental factor, according to the phenotype-gene relationships in the Genetic Association Database [1] or the environmental factor-gene relationships found in the Comparative Toxicogenomics Database[15], respectively. MeSH IDs are used as environmental factor node labels. Environmental factors with pharmacological or toxicological action in the MeSH record are color coded based on broad phenotype categories according to annotation in Additional file 2. *Phenotypes which do not fit into a specific category or environmental factors with undetermined pharmacological or toxicological action are gray. Request the Cytoscape session file of this network from Julia Gohlke gohlkej@niehs.nih.gov for more detailed exploration. Annotation of the circled metabolic syndrome and neuropsychiatric subnetworks can be found in Additional file 3.

Full size image

An important application of this work is generating hypotheses of interacting environmental factors that may be important in the prevention, initiation, progression, or treatment of complex diseases based on the network relationships found between phenotypes and environmental factors. Therefore, the tight cluster of metabolic syndrome phenotypes and neuropsychiatric disorders identified in Figure 3 are examined in further detail through analysis of gene expression datasets.

Metabolic syndrome cluster

Significance in both PPAR signaling and adipocytokine signaling form the tight subnetwork of 93 environmental factors linked to several metabolic syndrome phenotypes such as serum lipoprotein and triglyceride levels, body mass index, insulin sensitivity, type II diabetes, and obesity (see Additional File 3). Consistent with our results, a recent network analysis of microarray datasets from diabetes patients suggests PPAR signaling is the key underlying pathway in the pathogenesis of Type II diabetes [18]. Thiazolidinediones, which are antidiabetic PPARγ agonists [19, 20], the PPARα agonist fenofibrate and the HMG-CoA reductase inhibitor atorvastatin, both of which are used in the treatment of hyperlipidemia [21, 22] are identified in this subnetwork. Furthermore, dopamine antagonists, which includes several antipsychotic medications known to cause weight gain [23] are identified in this cluster. Retinoids are also found in this cluster, which is particularly intriguing in light of novel research showing retinaldehyde represses diet-induced obesity [24]. In fact, the widely used antineoplastic synthetic retinoic acid receptor alpha agonist Ro 41–5253 has recently been shown to induce PPARγ activity [25]. In addition, the increasing body of evidence linking Body Mass Index, retinoids, and cancer risk was recently highlighted in the most comprehensive analysis to date on diet and cancer risk [26]. In addition to the pharmaceuticals identified, di-n-hexyl phthalate (DHP), a widely used plasticizer that has recently been shown to act as a PPAR agonist [27] in addition to previous findings that high levels of exposure cause reproductive toxicity in animal models [28]. The Omega 3 fatty acids present in fish oil are an important dietary environmental factor identified in this cluster [29].

To test the hypothesis that regulation via PPAR and adipocytokine signaling plays an important role in environmental and genetic factors influencing metabolic syndrome phenotypes, we analyzed gene expression datasets after exposure to several predicted environmental modulators, as well as gene expression datasets from Familial combined hyperlipidemia cases, obese versus lean Pima Indians and obese versus lean mice fed a controlled diet (Table 1) [30–44]. Lists of significantly up or down regulated genes were submitted to DiRE http://dire.dcode.org/, a transcription factor binding site (tfbs) enrichment optimization algorithm that identifies tfbs that are enriched in evolutionary conserved regions surrounding a given set of genes versus a randomly generated background set of genes [45]. Lists of the tfbs enriched in the evolutionarily conserved regions surrounding the significantly up or downregulated gene list for each dataset are compiled in Additional file 4. Across all of these independent datasets, binding sites for the three transcriptional regulators of PPAR and adipocytokine signaling, namely PPAR, NFkB, and STAT, are consistently enriched in the differentially expressed gene sets (p ≤ 0.005) (Figure 4A). Therefore, this alternative analysis supports our previous subnetwork predictions suggesting a variety of environmental factors as well as genetic contributions to metabolic syndrome phenotypes can be integrated at the level of PPAR and adipocytokine signaling pathways. When the enriched tfbs identified for these metabolic syndrome subnetwork datasets are compared to enriched tfbs identified in the neuropsychiatric datasets (described below), we see that PPAR, PU.1 and FREAC binding sites are significantly enriched in these metabolic syndrome datasets (p ≤ 0.05).

Figure 4
figure4

Gene expression regulation across microarray datasets. Enriched transcription factor binding sites (tfbs) were identified in evolutionarily conserved regions surrounding differentially up and downregulated genes from metabolic (A.) or neuropsychiatric (B.) microarray datasets identified in Table 1 (see Methods). Results for each microarray dataset are presented in Additional file 4. The mean frequency of identifying a particular tfbs enriched in a dataset was 13% (dotted line). Those enriched tfbs that are consistently identified across the metabolic (A.) or neuropsychiatric (B.) datasets are color coded red (p ≤ 0.005), whereas those tfbs that are specific to the metabolic datasets versus neuropsychiatric datasets and vice versa (p ≤ 0.05) are identified with an asterisk.

Full size image
Table 1 Global gene expression datasets utilized for validation of metabolic syndrome and neuropsychiatric subnetworks
Full size table

Other tfbs beyond PPAR and adipocytokine signaling regulators that are highly enriched across these datasets offer hypotheses for future experimental research in the transcriptional regulation of metabolic syndrome phenotypes. For example, EBOX sites for basic helix-loop-helix transcription factor and PU.1, an ETS like tf, are important in cell fate programs in hematopoesis, particularly in the monocyte/macrophage lineage [46, 47]. This is intriguing in light of numerous studies showcasing the importance of macrophages not only in cardiovascular disease, but in the development of obesity as well[48, 49] and their connectivity to PPAR signaling [50, 51]. ZIC1 is a zinc finger transcription factor known to be important during early developmental programs [52], while preliminary genetic association work suggests ATF/CREB tfs may also play a role in obesity[53]. Finally, FREAC sites bind several forkhead members (FOXF2, FOXC1, FOXD1, AND FOXL1), which have been shown to be important in the regulation of gut-associated lymphoid organ development and regulation of intestinal glucose uptake in mice [54, 55].

Neuropsychiatric cluster

Our results suggest data from genetic association studies for several neuropsychiatric diseases (autism, schizophrenia, depression, bipolar disorder, attention deficit hyperactivity disorder, anxiety disorder, obsessive compulsive disorder, and Huntington's disease) converge on tyrosine metabolism and neuroactive ligand receptor interactions, forming a tight cluster of these phenotypes linked by significance in these two pathways. In fact, genes that code receptors and metabolic enzymes of the dopamine and serotonin signaling systems form the basis of this result. In contrast to the metabolic syndrome cluster, very few environmental factors (11) are found in this tight cluster and include the opiate pentazocine, the muscarinic receptor agonist pilocarpine, and the GABA modulator pentobarbital (Additional File 3). In addition, the acetylcholinesterase inhibiting organophosphates, well known for their use as pesticides, are identified in this cluster.

We analyzed gene expression datasets from case versus control studies for several of the phenotypes, as well as gene expression datasets generated from fetal astrocytes or rat forebrain after exposure to the organophosphate pesticide chlorpyrifos (Table 1). Following the method described for analysis of the gene expression datasets for the metabolic syndrome cluster, lists of significantly up or down regulated genes were submitted to DiRE [45]. Lists of enriched tfbs in regions surrounding the significantly up or downregulated gene lists for each dataset are available in Additional file 4. Across all of these datasets, enrichment for EBOX and MEF2 binding sites are found most consistently in the differentially expressed genes for the neuropsychiatric cluster datasets (Figure 4B). Consistent with the result, several studies suggest coordinated action of the EBOX binding proneural bHLH transcriptional activators and Mef2c in the differentiation of neuronal subtypes in the developing mammalian forebrain [56–59]. When the enriched tfbs identified for these neuropsychiatric subnetwork datasets are compared to enriched tfbs identified in the metabolic syndrome datasets described above, we see that only chicken ovalbumin upstream promoter transcription factor (COUP) binding sites are significantly enriched in neuropsychiatric datasets (p ≤ 0.05). COUP-TFs are members of the steroid receptor superfamily in which dopamine is thought to be a physiological activator [60].

Consideration of bias associated with genetic association studies

One potential source of bias is the likelihood of false positive associations represented in the GAD database. For example, a large multi-center study could not validate several previously reported genetic risk factors for acute coronary syndrome [61]. In addition, publication of false positive events could lead to more extensive publication bias as these results are followed up for related phenotypes. To address this potential for bias, we have re-evaluated 3 phenotypes for which extensive meta-analyses of genome-wide association (GWA) findings exist [62–64]. The list of genes from these recent meta-analyses of GWA studies for Alzheimer's Disease, Parkinson's Disease, and Schizophrenia represent a subset of those KEGG represented genes found within GAD, as only one novel gene for each phenotype was identified by the GWA meta-analysis across these phenotypes (Additional file 5). Subsequently, the SEPEA algorithm was re-run using the genes associated with each phenotype based on these meta-analyses of GWA results. In general, the predicted enriched pathways were consistent across results obtained with the GAD and GWA generated lists, however there were some notable differences (Additional file 5). For example, for both the GAD and GWA Alzheimer's Disease gene lists, the Renin-Angiotensin System pathway ranked the highest, however, using the results from the meta-analysis of GWA studies suggests folate related pathways may be important whereas the GAD results suggest tyrosine metabolism is altered. In fact, tyrosine metabolism ranks high for all three phenotypes using the GAD generated gene list, whereas Parkinson's Disease is the only phenotype in which the GWA studies confirm this result. As more GWA results become available for other phenotypes, this potential limitation of the current analysis can be more fully evaluated.

Pathways to disease

Ultimately, a particular phenotype is produced by the integration of outputs from a multitude of molecular pathways within an organism. Therefore, we explored the higher order structure of pathway networks by overlaying our analysis onto the network structure of interconnected KEGG pathways (Figure 5). This analysis allows us to simultaneously visualize the key pathways to complex disease progression from the genetic standpoint by adjusting node size to reflect the number of human phenotypes associated with a particular pathway based on the sum of disease associated genetic polymorphisms, as well as from the environmental standpoint, by adjusting the color of the pathway node to reflect the number of environmental factors associated with a particular pathway.

Figure 5
figure5

Pathway Interaction Network. Nodes represent KEGG metabolic and signaling pathways and are connected based on the KEGG database [72]. Node size is reflective of the number of phenotypes associated with the particular pathway based on application of SEPEA (p ≤ 0.003) to gene lists annotated from The Genetic Association Database[1]. Intensity of node color is reflective of the number of environmental factors associated with the particular pathway based on enrichment of gene lists annotated from The Comparative Toxicogenomics Database[15].

Full size image

Looking at the intersection of the top 15 pathways most often enriched in genetic association studies and environmental factor research (Table 2), suggests metabolism of xenobiotics by cytochrome P450, retinol metabolism, Jak-STAT signaling, Toll-like receptor signaling, and adipocytokine signaling may be five critical pathways important to disease progression from both a genetic and environmental standpoint. From our analysis of phenotypes illustrated in Figure 1, we see that metabolism of xenobiotics by cytochrome P450 is significantly enriched in genetic association datasets for several phenotypes including cancers, cardiovascular disease, and immune related disorders. Adipocytokine signaling defines the cardiovascular and metabolic syndrome phenotypes, many of which have reached epidemic levels over the last 30 years [65], suggesting environmental components are critical in the etiology of these phenotypes. Retinol metabolism is significantly enriched in genetic polymorphism lists for hormonally regulated cancers such as breast, endometrial, testicular, prostate and thyroid, as well as pregnancy complications, reproductive dysfunction, and cardiovascular and endocrine disorders. This group of phenotypes is particularly interesting in light of the latest time trend statistics from the National Cancer Institute and Centers for Disease Control. As a whole, cancer incidence rates have been declining over the last decade, with the exception of 5 sites (thyroid, liver, kidney, skin, and testis). Thyroid cancer has by far the largest increase in incidence over the last decade, with an annual percent change of 5.3 between 1994 to 2004 [66]. In addition, pregnancy complications and endocrine disorders account for 5 of the 6 primary diagnoses with the greatest percent increase in ambulatory care visits over the last decade [67]. These time trends suggest environmental components are critical in rising incidence of endocrine related phenotypes, as the timeframe hardly crosses a generation, highlighting the importance of continued research on exposure routes and health effects of potential endocrine disrupters found in our environment [68, 69].

Table 2 Top pathways enriched using genetic association research or environmental factor research.
Full size table

Finally, we note the centrality of PPAR and adipocytokine signaling in the pathway network as the primary linkage between metabolism and cellular signaling pathways (Figure 5). As mentioned previously, these two pathways define the metabolic syndrome cluster in Figure 3. Several genetic and environmental factors are associated with each of these pathways, suggesting genetic and environmental modulators are critical to the role of these pathways in human disease progression, such as in metabolic syndrome phenotypes and cardiovascular disease.

Conclusion

According to systems theory, although individual genes or environmental factors may be a critical component in the pathogenesis of a particular complex disease, it is ultimately the modulation of underlying pathways that the particular gene/environmental factor is a part of that determines the resultant phenotype. Here we have integrated gene centered knowledge from epidemiological and mechanistic environmental research in an attempt to discover the interplay between genetic and environmental mediators of phenotype at the pathway level. In addition, we have provided a higher order structure of pathway interconnectivity to build hypotheses of disease progression based on clusters of pathways defining phenotypes.

The methods and findings presented here open the door to a number of new hypotheses that can be explored regarding the genetic and environmental factors governing human disease. The results suggest retinol metabolism, Jak-STAT signaling, Toll-like receptor signaling, and adipocytokine signaling are key pathways that should be prioritized targets for high-throughput screening currently being implemented to improve toxicity testing [70, 71]. For example, analysis of the metabolic syndrome subnetwork highlights the need for further epidemiological and mechanistic analyses of several compounds for their potential modulation of metabolic syndrome phenotypes, including plastic derivatives, synthetic and natural retinoids, pyrethrins and antipsychotic medications. In addition, the role of endocrine pathways in numerous phenotypes for which rates have increased over the last 30 years indicates a continued need to evaluate in greater detail the role of endocrine disruption in cancer, pregnancy and reproductive complications, and metabolic syndrome phenotypes. The multifactorial nature of complex diseases necessitates using knowledge-based, systems-driven evaluations, like the one presented here, for uncovering promising hypotheses for future research aimed at improving public health.

Methods

Characterization of Phenotype-Gene Relationships

The Genetic Association Database is an NIH supported gene-centered public repository of human association studies examining a wide range of human phenotypes, including non-mendelian common diseases, and is one of the largest databases of human disease associated polymorphisms currently available. All gene-phenotype relationships (N = 28,341) in the Genetic Association Database were downloaded (June 8, 2007). Phenotypes were further annotated to collapse synonyms, as well as group similar phenotypes into categories (see mapping used in Additional file 1). Only those phenotypes with at least 3 unique genes associated with it were analyzed further, resulting in 10,089 unique phenotype-gene relationships used in subsequent analyses.

Characterization of Environment-Gene Relationships

The Comparative Toxicogenomics Database is an NIH supported public database that provides curated interactions between environmental factors and genes or proteins. Using either a MeSH concept or descriptor as the environmental factor identifier, all unique environmental factor-gene/protein relationships as of June 2007 (N = 47,025) were evaluated to define relationships between environmental exposures and human genes.

Annotation of MeSH concepts or descriptors was performed using the 2007–2008 MeSH browser http://www.nlm.nih.gov/mesh/MBrowser.html to identify any known biological actions of the environmental factors within the MeSH record. All environmental factors described in the Comparative Toxicogenomics Database fall within the Chemicals and Drugs [D] heading. To identify the most biologically relevant categorization, priority for annotation was set as follows: Noxae [D27.888.569], Physiological Effects of Drugs [D27.505.696], Therapeutic Uses [D27.505.954], Molecular Mechanisms of Action [D27.505.519]. If no information was available within these categories, then annotation by substance structure using all other trees under Chemicals and Drugs was implemented to annotate the given environmental factor.

Evaluation of Gene-Pathway relationships

All sets of genes associated with a particular phenotype or environmental factor were analyzed for over-representation in specific molecular pathways found in the KEGG database [72] using Structurally Enhanced Pathway Enrichment Analysis (SEPEA), a novel pathway enrichment algorithm that incorporates relationships between nodes within a pathway using specific scoring rules described in detail elsewhere [13]. Briefly, the heavy ends scoring rule gives more importance to genes at the beginning (e.g. receptors) or end (e.g. transcription factors) of a pathway and the distance scoring rule gives more importance to those pathways where the perturbed genes (for a given condition) are close relative to each other in the pathway network. In this application, we use the SEPEA_NT3 method (see [13] for a detailed description). Broadly, the null hypothesis states that the distribution of the number of perturbed genes for a given condition in a specific pathway is not different from the distribution of a random set of genes chosen from all the genes involved in the KEGG pathways analyzed, in the context of the rules described above. Here we are incorporating a heavy ends scoring rule using a power function (0.5δ, with δ being the distance from a terminal node in the pathway), instead of a linear function as described in [13] Equation 8. Utilization of this power function emphasizes the underlying hypothesis of the biological importance of this rule in the final significance evaluation. We found that this emphasis resulted in a network with more clear separation of phenotypes and chemicals. In this analysis, those KEGG pathways developed based on a particular phenotype (disease pathways) were eliminated based on the potential redundancy of information found in the Genetic Association Database. Furthermore, only those KEGG pathways that had at least 3 human genes associated with them were analyzed.

Pathway-Phenotype-Cluster

To determine relationships between human phenotypes based on polymorphisms associated with those phenotypes, phenotype p-values for KEGG pathway enrichment were clustered using Spearman Rank correlation with average linkage using Cluster version 2.11 and viewed using TreeView [73] downloaded July 2007 from http://rana.lbl.gov/EisenSoftware.htm

Phenotype-Environmental factor Network

Each phenotype-phenotype, phenotype-environmental factor, or environmental factor-environmental factor pair with at least two common significant pathways was assigned an edge. Network connectivity between phenotypes and environmental factors were determined for a range of SEPEA pathway enrichment p-value cut-offs and the sensitivity and specificity of the results as compared to manually curated, direct chemical-disease relationships found in the Comparative Toxicogenomics Database (CTD) database (downloaded in September 2008)[16]. This dataset contains direct chemical-disease relationships reported in the literature in human and animal model studies. We reduced this CTD database to those diseases found in GAD, which resulted in 1084 CTD relationships. We further collapsed these into the 8 broad disease categories. This analysis was used to establish the optimal pathway enrichment significance p-value cutoff of less than 0.00321 (which corresponds to a FDR ≤ 0.32 as computed using the Benjamini-Hochberg method [74]). This FDR is comparable to other pathway enrichment algorithm cut-offs (e.g. 0.25 for GSEA [75]), and is considered acceptable if one is primarily interested in hypothesis generation. The significance of these results was evaluated using the binomial cumulative distribution function where the probability (pk) of observing at least 226 (or 271) significant chemical-disease (or chemical-disease category) relationships by random chance was determined using the data for 4559 chemicals from CTD and for 204 diseases in GAD. Graphical representation of the network was determined using the edge weighted spring embedded algorithm in Cytoscape 2.5.0 downloaded Aug. 2007 [76] with the following parameters: spring strength = 5.0, spring rest length = 10.0, rest length of a disconnected spring = 1500, and strength of a disconnected spring = 0.06. Only those environmental factors with at least 2 genes associated with them were included in the final representation.

Statistical Evaluation of Network

Each environmental factor or phenotype was labeled with one of 9 categories (Cardiovascular, Neurologic/Psychiatric, Neoplastic, Metabolic/Gastrointestinal, Immunologic, Hematologic, Endocrine/Reproductive, Pulmonary, Other) based on their classification in knowledge of the phenotype or Mesh categorizations for environmental factors as described above (See Additional file 2). The graph clustering coefficient method described in [17] was used to statistically evaluate the network generated. Briefly, p-values for the observed graph clustering coefficient for a given disease category (i) which has n(i) nodes associated with it in the network are based on choosing n(i) nodes randomly from the large network and computing the clustering coefficient for this random subset of nodes (based on 1000 permutations).

Identification of enriched transcription factor binding sites in differentially expressed genes from microarray datasets

Up and downregulated gene lists from the microarray data as described in Table 1 was accessed from the publication associated with the datasets [30–36, 77–84], or via downloading from GEO in Feb. or Nov. 2008. In the latter case, differentially expressed genes were identified using mattest and mafdr (Matlab 7.4.0.287 (R2007a)) with a fold change cutoff of 1.5 and a q value cutoff of 0.10. Lists of up or downregulated genes for each dataset were then submitted to DiRE http://dire.dcode.org/, a transcription factor binding site (tfbs) enrichment optimization algorithm that identifies tfbs that are enriched in evolutionary conserved regions surrounding a given set of genes versus a background set of genes [45]. Based on analyses using tissue specificity gene expression datasets, the developers of DiRE show that an importance score cutoff of 0.10 is reasonable for achieving good specificity and precision [45], therefore we used this cutoff to identify those tfbs that are significantly enriched in sequence surrounding up or downregulated genes from each dataset (Additional File 4). Similar TRANSFAC binding sites [85] were then collapsed to avoid double counting tfs with similar binding sites (Additional File 6). Significance of consistently identifying a particular tf class across the metabolic or psychiatric datasets was tested using a binomial cumulative distribution, where the probability of observing at least x number of significant tfbs for a given tf across significance lists generated from the 12 (metabolic) or 9 (neuropsychiatric) datasets with a probability parameter equal to the mean frequency of occurrence of any tf that satisfied the importance score cutoff of 0.1 among all the datasets (0.137) was computed. To test the specificity of our findings as they relate to the metabolic syndrome datasets versus neuropsychiatric datasets analyzed, we compared the probability of finding a specific tfbs at least the observed number of times using the metabolic syndrome datasets versus the neuropsychiatric datasets and vice versa using a hypergeometric distribution.

Generation of Pathway Interconnectivity Network

Connectivity between pathways was downloaded from KEGG (Aug. 16, 2007). Pathway network layout was generated using the force directed algorithm in Cytoscape 2.5.0 [76], where node size reflects the number of phenotypes associated with each pathway and node color gradation reflects the number of environmental factors associated with each pathway using a p-value cutoff of 0.003.

References

  1. 1.

    Becker KG, Barnes KC, Bright TJ, Wang SA: The Genetic Association Database. Nature Genetics. 2004, 36 (5): 431-432. 10.1038/ng0504-431

  2. 2.

    McKusick VA: Mendelian inheritance in man and its online version, OMIM. American Journal of Human Genetics. 2007, 80 (4): 588-604. 10.1086/514346

  3. 3.

    Antonarakis SE, Beckmann JS: Opinion – Mendelian disorders deserve more attention. Nature Reviews Genetics. 2006, 7 (4): 277-282. 10.1038/nrg1826

  4. 4.

    Gwinn M, Khoury MJ: Genomics and public health in the United States: Signposts on the translation highway. Community Genetics. 2006, 9 (1): 21-26. 10.1159/000090689

  5. 5.

    Kleeberger SR, Peden D: Gene-environment interactions in asthma and other respiratory diseases. Annual Review of Medicine. 2005, 56: 383-400. 10.1146/annurev.med.56.062904.144908

  6. 6.

    Infante-Rivard C, Mathonnet G, Sinnett D: Risk of childhood leukemia associated with diagnostic irradiation and polymorphisms in DNA repair genes. Environmental Health Perspectives. 2000, 108 (6): 495-498. 10.2307/3454609

  7. 7.

    Infante-Rivard C, Labuda D, Krajinovic M, Sinnett D: Risk of childhood leukemia associated with exposure to pesticides and with gene polymorphisms. Epidemiology. 1999, 10 (5): 481-487. 10.1097/00001648-199909000-00004

  8. 8.

    Mathers CD, Loncar D: Projections of global mortality and burden of disease from 2002 to 2030. Plos Medicine. 2006, 3 (11):

  9. 9.

    Butte AJ, Kohane IS: Creation and implications of a phenome-genome network. Nature Biotechnology. 2006, 24 (1): 55-62. 10.1038/nbt1150

  10. 10.

    Lage K, Karlberg EO, Storling ZM, Olason PI, Pedersen AG, Rigina O, Hinsby AM, Tumer Z, Pociot F, Tommerup N, et al.: A human phenome-interactome network of protein complexes implicated in genetic disorders. Nature Biotechnology. 2007, 25 (3): 309-316. 10.1038/nbt1295

  11. 11.

    Goh KI, Cusick ME, Valle D, Childs B, Vidal M, Barabasi AL: The human disease network. Proceedings of the National Academy of Sciences of the United States of America. 2007, 104 (21): 8685-8690. 10.1073/pnas.0701361104

  12. 12.

    Parmigiani G, Boca S, Lin J, Kinzler KW, Velculescu V, Vogelstein B: Design and analysis issues in genome-wide somatic mutation studies of cancer. Genomics. 2009, 93 (1): 17-21. 10.1016/j.ygeno.2008.07.005

  13. 13.

    Thomas R, Gohlke JM, Stopper GF, Parham FM, Portier CJ: Choosing the right path: enhancement of biologically-relevant sets of genes or proteins using pathway structure. Genome Biology. 2009, 10: R44- 10.1186/gb-2009-10-4-r44

  14. 14.

    Ritchie SA, Connell JMC: The link between abdominal obesity, metabolic syndrome and cardiovascular disease. Nutrition Metabolism and Cardiovascular Diseases. 2007, 17 (4): 319-326. 10.1016/j.numecd.2006.07.005.

  15. 15.

    Mattingly CJ, Rosenstein MC, Davis AP, Colby GT, Forrest JN, Boyer JL: The Comparative Toxicogenomics Database: A cross-species resource for building chemical-gene interaction networks. Toxicological Sciences. 2006, 92 (2): 587-595. 10.1093/toxsci/kfl008

  16. 16.

    Davis AP, Murphy CG, Saraceni-Richards CA, Rosenstein MC, Wiegers TC, Mattingly CJ: Comparative Toxicogenomics Database: a knowledgebase and discovery tool for chemical-gene-disease networks. Nucleic Acids Res. 2008, D786-92. Epub 2008 Sep 9, 37 Database

  17. 17.

    Watts DJ, Strogatz SH: Collective dynamics of 'small-world' networks. Nature. 1998, 393 (6684): 440-442. 10.1038/30918

  18. 18.

    Liu M, Liberzon A, Kong SW, Lai WR, Park PJ, Kohane IS, Kasif S: Network-Based Analysis of Affected Biological Processes in Type 2 Diabetes Models. PLoS Genetics. 2007, 3 (6): e96- 10.1371/journal.pgen.0030096

  19. 19.

    Sorbera LA, Leeson PA, Martin L, Castaner J: Farglitazar – Antidiabetic, PPAR gamma agonist. Drugs of the Future. 2001, 26 (4): 354-363. 10.1358/dof.2001.026.04.617323.

  20. 20.

    Oberfield JL, Collins JL, Holmes CP, Goreham DM, Cooper JP, Cobb JE, Lenhard JM, Hull-Ryde EA, Mohr CP, Blanchard SG, et al.: A peroxisome proliferator-activated receptor gamma ligand inhibits adipocyte differentiation. Proceedings of the National Academy of Sciences of the United States of America. 1999, 96 (11): 6102-6106. 10.1073/pnas.96.11.6102

  21. 21.

    Guay DRP: Update on fenofibrate. Cardiovascular Drug Reviews. 2002, 20 (4): 281-302.

  22. 22.

    Bakker-Arkema RG, Davidson MH, Goldstein RJ, Davignon J, Isaacsohn JL, Weiss SR, Keilson LM, Brown WV, Miller VT, Shurzinske LJ, et al.: Efficacy and safety of a new HMG-CoA reductase inhibitor, atorvastatin, in patients with hypertriglyceridemia. Jama. 1996, 275 (2): 128-133. 10.1001/jama.275.2.128

  23. 23.

    Baptista T, Kin N, Beaulieu S, de Baptista EA: Obesity and related metabolic abnormalities during antipsychotic drug administration: Mechanisms, management and research perspectives. Pharmacopsychiatry. 2002, 35 (6): 205-219. 10.1055/s-2002-36391

  24. 24.

    Ziouzenkova O, Orasanu G, Sharlach M, Akiyama TE, Berger JP, Viereck J, Hamilton JA, Tang GW, Dolnikowski GG, Vogel S, et al.: Retinaldehyde represses adipogenesis and diet-induced obesity. Nature Medicine. 2007, 13 (6): 695-702. 10.1038/nm1587

  25. 25.

    Schupp M, Curtin JC, Kim RJ, Billin AN, Lazar MA: A widely used retinoic acid receptor antagonist induces peroxisome proliferator-activated receptor-gamma activity. Molecular Pharmacology. 2007, 71 (5): 1251-1257. 10.1124/mol.106.033662

  26. 26.

    , : Food, Nutrition, Physical Activity, and the Prevention of Cancer: a Global Perspective. 2007, Washington, DC: World Cancer Research Fund/American Institute for Cancer Research

  27. 27.

    Boberg J, Metzdorff S, Wortziger R, Axelstad M, Brokken L, Vinggaard AM, Dalgaard M, Nellemann C: Impact of diisobutyl phthalate and other PPAR agonists on steroidogenesis and plasma insulin and leptin levels in fetal rats. Toxicology. 2008, 250 (2–3): 75-81. 10.1016/j.tox.2008.05.020

  28. 28.

    NTP-CERHR Monograph on the Potential Human Reproductive and Developmental Effects of Di-n-Hexyl Phthalate (DnHP). Ntp Cerhr Mon. 2003, i-III90. 7

  29. 29.

    Neschen S, Morino K, Rossbacher JC, Pongratz RL, Cline GW, Sono S, Gillum M, Shulman GI: Fish oil regulates adiponectin secretion by a peroxisome proliferator-activated receptor-gamma-dependent mechanism in mice. Diabetes. 2006, 55 (4): 924-928. 10.2337/diabetes.55.04.06.db05-0985

  30. 30.

    Lee YH, Nair S, Rousseau E, Allison DB, Page GP, Tataranni PA, Bogardus C, Permana PA: Microarray profiling of isolated abdominal subcutaneous adipocytes from obese vs non-obese Pima Indians: increased expression of inflammation-related genes. Diabetologia. 2005, 48 (9): 1776-1783. 10.1007/s00125-005-1867-3

  31. 31.

    Koza RA, Nikonova L, Hogan J, Rim JS, Mendoza T, Faulk C, Skaf J, Kozak LP: Changes in gene expression foreshadow diet-induced obesity in genetically identical mice. Plos Genetics. 2006, 2 (5): 769-780. 10.1371/journal.pgen.0020081.

  32. 32.

    Llaverias G, Pou J, Ros E, Zambon D, Cofan M, Sanchez A, Vazquez-Carrera M, Sanchez RM, Laguna JC, Alegret M: Monocyte gene-expression profile in men with familial combined hyperlipidemia and its modification by atorvastatin treatment. Pharmacogenomics. 2008, 9 (8): 1035-1054. 10.2217/14622416.9.8.1035

  33. 33.

    Fielden MR, Brennan R, Gollub J: A gene expression biomarker provides early prediction and mechanistic assessment of hepatic tumor induction by nongenotoxic chemicals. Toxicological Sciences. 2007, 99 (1): 90-100. 10.1093/toxsci/kfm156

  34. 34.

    Wang Y, Yao RS, Maciag A, Grubbs CJ, Lubet RA, You M: Organ-specific expression profiles of rat mammary gland, liver, and lung tissues treated with targretin, 9-cis retinoic acid, and 4-hydroxyphenylretinamide. Molecular Cancer Therapeutics. 2006, 5 (4): 1060-1072. 10.1158/1535-7163.MCT-05-0322

  35. 35.

    McClintick JN, Crabb DW, Tian HJ, Pinaire J, Smith JR, Jerome RE, Edenberg HJ: Global effects of vitamin A deficiency on gene expression in rat liver: evidence for hypoandrogenism. Journal of Nutritional Biochemistry. 2006, 17 (5): 345-355. 10.1016/j.jnutbio.2005.08.006

  36. 36.

    Bordoni A, Astolfi A, Morandi L, Pession A, Danesi F, Di Nunzio M, Franzoni M, Biagi P, Pession A: N-3 PUFAs modulate global gene expression profile in cultured rat cardiomyocytes. Implications in cardiac hypertrophy and heart failure. Febs Letters. 2007, 581 (5): 923-929. 10.1016/j.febslet.2007.01.070

  37. 37.

    Hsiao A, Worrall DS, Olefsky JM, Subramaniam S: Variance-modeled posterior inference of microarray data: detecting gene-expression changes in 3T3-L1 adipocytes. Bioinformatics. 2004, 20 (17): 3108-3127. 10.1093/bioinformatics/bth371

  38. 38.

    Mense SM, Sengupta A, Lan C, Zhou M, Bentsman G, Volsky DJ, Whyatt RM, Perera FP, Zhang L: The common insecticides cyfluthrin and chlorpyrifos alter the expression of a subset of genes with diverse functions in primary human astrocytes. Toxicol Sci. 2006, 93 (1): 125-135. 10.1093/toxsci/kfl046

  39. 39.

    Iwamoto K, Kakiuchi C, Bundo M, Ikeda K, Kato T: Molecular characterization of bipolar disorder by comparing gene expression profiles of postmortem brains of major mental disorders. Mol Psychiatry. 2004, 9 (4): 406-416. 10.1038/sj.mp.4001437

  40. 40.

    Saetre P, Emilsson L, Axelsson E, Kreuger J, Lindholm E, Jazin E: Inflammation-related genes up-regulated in schizophrenia brains. BMC Psychiatry. 2007, 7: 46- 10.1186/1471-244X-7-46

  41. 41.

    Hovatta I, Tennant RS, Helton R, Marr RA, Singer O, Redwine JM, Ellison JA, Schadt EE, Verma IM, Lockhart DJ, et al.: Glyoxalase 1 and glutathione reductase 1 regulate anxiety in mice. Nature. 2005, 438 (7068): 662-666. 10.1038/nature04250

  42. 42.

    Nishimura Y, Martin CL, Vazquez-Lopez A, Spence SJ, Alvarez-Retuerto AI, Sigman M, Steindler C, Pellegrini S, Schanen NC, Warren ST, et al.: Genome-wide expression profiling of lymphoblastoid cell lines distinguishes different forms of autism and reveals shared pathways. Hum Mol Genet. 2007, 16 (14): 1682-1698. 10.1093/hmg/ddm116

  43. 43.

    Gregg JP, Lit L, Baron CA, Hertz-Picciotto I, Walker W, Davis RA, Croen LA, Ozonoff S, Hansen R, Pessah IN, et al.: Gene expression changes in children with autism. Genomics. 2008, 91 (1): 22-29. 10.1016/j.ygeno.2007.09.003

  44. 44.

    Stapleton AR, Chan VT: Subtoxic chlorpyrifos treatment resulted in differential expression of genes implicated in neurological functions and development. Arch Toxicol. 2008, 83 (4): 319-33. Epub 2008 Jul 31. 10.1007/s00204-008-0346-2

  45. 45.

    Pennacchio LA, Loots GG, Nobrega MA, Ovcharenko I: Predicting tissue-specific enhancers in the human genome. Genome Research. 2007, 17 (2): 201-211. 10.1101/gr.5972507

  46. 46.

    O'Neil J, Look AT: Mechanisms of transcription factor deregulation in lymphoid cell transformation. Oncogene. 2007, 26 (47): 6838-6849. 10.1038/sj.onc.1210766

  47. 47.

    Rosenbauer F, Tenen DG: Transcription factors in myeloid development: balancing differentiation with transformation. Nature Reviews Immunology. 2007, 7 (2): 105-117. 10.1038/nri2024

  48. 48.

    Cancello R, Clement K: Is obesity an inflammatory illness? Role of low-grade inflammation and macrophage infiltration in human white adipose tissue. Bjog-an International Journal of Obstetrics and Gynaecology. 2006, 113 (10): 1141-1147. 10.1111/j.1471-0528.2006.01004.x.

  49. 49.

    Bastard JP, Maachi M, Lagathu C, Kim MJ, Caron M, Vidal H, Capeau J, Feve B: Recent advances in the relationship between obesity, inflammation, and insulin resistance. European Cytokine Network. 2006, 17 (1): 4-12.

  50. 50.

    Li AC, Palinski W: Peroxisome proliferator-activated receptors: How their effects on macrophages can lead to the development of a new drug therapy against atherosclerosis. Annual Review of Pharmacology and Toxicology. 2006, 46: 1-39. 10.1146/annurev.pharmtox.46.120604.141247.

  51. 51.

    Lehrke M, Lazar MA: The many faces of PPAR gamma. Cell. 2005, 123 (6): 993-999. 10.1016/j.cell.2005.11.026

  52. 52.

    Merzdorf CS: Emerging roles for zic genes in early development. Developmental Dynamics. 2007, 236 (4): 922-940. 10.1002/dvdy.21098

  53. 53.

    Rousset S, Gonzalez-Barroso MD, Gelly C, Pecqueur C, Bouillaud F, Ricquier D, Cassard-Doulcier AM: A new polymorphic site located in the human UCP1 gene controls the in vitro binding of CREB-like factor. International Journal of Obesity. 2002, 26 (5): 735-738. 10.1038/sj.ijo.0801973

  54. 54.

    Katz JP, Perreault N, Goldstein BG, Chao HH, Ferraris RP, Kaestner KH: Foxl1 null mice have abnormal intestinal epithelia, postnatal growth retardation, and defective intestinal glucose uptake. Am J Physiol Gastrointest Liver Physiol. 2004, 287 (4): G856-864. 10.1152/ajpgi.00136.2004

  55. 55.

    Fukuda K, Yoshida H, Sato T, Furumoto TA, Mizutani-Koseki Y, Suzuki Y, Saito Y, Takemori T, Kimura M, Sato H, et al.: Mesenchymal expression of Foxl1, a winged helix transcriptional factor, regulates generation and maintenance of gut-associated lymphoid organs. Dev Biol. 2003, 255 (2): 278-289. 10.1016/S0012-1606(02)00088-X

  56. 56.

    Li H, Radford JC, Ragusa MJ, Shea KL, McKercher SR, Zaremba JD, Soussou W, Nie Z, Kang YJ, Nakanishi N, et al.: Transcription factor MEF2C influences neural stem/progenitor cell differentiation and maturation in vivo. Proc Natl Acad Sci USA. 2008, 105 (27): 9397-9402. 10.1073/pnas.0802876105

  57. 57.

    Skerjanc IS, Wilton S: Myocyte enhancer factor 2C upregulates MASH-1 expression and induces neurogenesis in P19 cells. FEBS Lett. 2000, 472 (1): 53-56. 10.1016/S0014-5793(00)01438-1

  58. 58.

    Mattar P, Britz O, Johannes C, Nieto M, Ma L, Rebeyka A, Klenin N, Polleux F, Guillemot F, Schuurmans C: A screen for downstream effectors of Neurogenin2 in the embryonic neocortex. Developmental Biology. 2004, 273 (2): 373-389. 10.1016/j.ydbio.2004.06.013

  59. 59.

    Gohlke JM, Armant O, Parham FM, Smith MV, Zimmer C, Castro DS, Nguyen L, Parker JS, Gradwohl G, Portier CJ, et al.: Characterization of the proneural gene regulatory network during mouse telencephalon development. Bmc Biology. 2008, 6: 15- 10.1186/1741-7007-6-15

  60. 60.

    Power RF, Lydon JP, Conneely OM, O'Malley BW: Dopamine activation of an orphan of the steroid receptor superfamily. Science. 1991, 252 (5012): 1546-1548. 10.1126/science.2047861

  61. 61.

    Morgan TM, Krumholz HM, Lifton RP, Spertus JA: Nonvalidation of reported genetic risk factors for acute coronary syndrome in a large-scale replication study. Jama. 2007, 297 (14): 1551-1561. 10.1001/jama.297.14.1551

  62. 62.

    Allen NC, Bagade S, McQueen MB, Ioannidis JP, Kavvoura FK, Khoury MJ, Tanzi RE, Bertram L: Systematic meta-analyses and field synopsis of genetic association studies in schizophrenia: the SzGene database. Nat Genet. 2008, 40 (7): 827-834. 10.1038/ng.171

  63. 63.

    Bertram L, McQueen MB, Mullin K, Blacker D, Tanzi RE: Systematic meta-analyses of Alzheimer disease genetic association studies: the AlzGene database. Nat Genet. 2007, 39 (1): 17-23. 10.1038/ng1934

  64. 64.

    Bagade S, Allen NC, Tanzi RE, Bertram L: The PDGene Database. Alzheimer Research Forum. 2008,http://www.pdgene.org/dbindex.asp

  65. 65.

    Gillespie KM, Bain SC, Barnett AH, Bingley PJ, Christie MR, Gill GV, Gale EAM: The rising incidence of childhood type 1 diabetes and reduced contribution of high-risk HLA haplotypes. Lancet. 2004, 364 (9446): 1699-1700. 10.1016/S0140-6736(04)17357-1

  66. 66.

    Ries LAG, Melbert D, Krapcho M, Mariotto A, Miller BA, Feuer EJ, Clegg L, Horner MJ, Howlader N, Eisner MP, et al.: SEER Cancer Statistics Revew, 1975–2004. 2007, Bethesda, MD: National Cancer Institute

  67. 67.

    Schappert SM, Burt CW: Ambulatory care visits to physician offices, hospital outpatient departments, and emergency departments: United States, 2001–02. Vital Health Stat 13. 2006, 1-66. 159

  68. 68.

    Aragon A, Martinez E, Sanderson J, Frausto S, Wolff C, Savage DD: Effects of prenatal ethanol exposure on G-protein coupled receptor function. Alcoholism-Clinical and Experimental Research. 2004, 28 (5): 43A-43A.

  69. 69.

    Daston GP, Cook JC, Kavlock RJ: Uncertainties for endocrine disrupters: Our view on progress. Toxicological Sciences. 2003, 74 (2): 245-252. 10.1093/toxsci/kfg105

  70. 70.

    , : Toxicity Testing in the 21st Century, A vision and a strategy. Edited by: Council NR. 2007, 196-Washington, DC: National Academies Press

  71. 71.

    , : A National Toxicology Program for the 21st Century:A roadmap to achieve the NTP vision. 2004, Research Triangle Park, NC: National Toxicology Program/National Institute of Environmental Health Sciences

  72. 72.

    Kanehisa M, Goto S, Hattori M, Aoki-Kinoshita KF, Itoh M, Kawashima S, Katayama T, Araki M, Hirakawa M: From genomics to chemical genomics: new developments in KEGG. Nucleic Acids Research. 2006, 34: D354-D357. 10.1093/nar/gkj102

  73. 73.

    Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proceedings of the National Academy of Sciences of the United States of America. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863

  74. 74.

    Benjamini Y, Hochberg Y: Controlling the False Discovery Rate – a Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B-Methodological. 1995, 57 (1): 289-300.

  75. 75.

    Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al.: Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102

  76. 76.

    Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Research. 2003, 13 (11): 2498-2504. 10.1101/gr.1239303

  77. 77.

    Stapleton AR, Chan VT: Subtoxic chlorpyrifos treatment resulted in differential expression of genes implicated in neurological functions and development. Archives of Toxicology. 2009, 83 (4): 319-333. 10.1007/s00204-008-0346-2

  78. 78.

    Gregg JP, Lit L, Baron CA, Hertz-Picciotto I, Walker W, Davis RA, Croen LA, Ozonoff S, Hansen R, Pessah IN, et al.: Gene expression changes in children with autism. Genomics. 2008, 91 (1): 22-29. 10.1016/j.ygeno.2007.09.003

  79. 79.

    Saetre P, Emilsson L, Axelsson E, Kreuger J, Lindholm E, Jazin E: Inflammation-related genes up-regulated in schizophrenia brains. Bmc Psychiatry. 2007, 7: 46- 10.1186/1471-244X-7-46

  80. 80.

    Nishimura Y, Martin CL, Lopez AV, Spence SJ, Alvarez-Retuerto AI, Sigman M, Steindler C, Pellegrini S, Schanen NC, Warren ST, et al.: Genome-wide expression profiling of lymphoblastoid cell lines distinguishes different forms of autism and reveals shared pathways. Human Molecular Genetics. 2007, 16 (14): 1682-1698. 10.1093/hmg/ddm116

  81. 81.

    Mense SM, Sengupta A, Lan CG, Zhou M, Bentsman G, Volsky DJ, Whyatt RM, Perera FP, Zhang L: The common insecticides cyfluthrin and chlorpyrifos alter the expression of a subset of genes with diverse functions in primary human astrocytes. Toxicological Sciences. 2006, 93 (1): 125-135. 10.1093/toxsci/kfl046

  82. 82.

    Hovatta I, Tennant RS, Helton R, Marr RA, Singer O, Redwine JM, Ellison JA, Schadt EE, Verma IM, Lockhart DJ, et al.: Glyoxalase 1 and glutathione reductase 1 regulate anxiety in mice. Nature. 2005, 438 (7068): 662-666. 10.1038/nature04250

  83. 83.

    Hsiao A, Worrall DS, Olefsky JM, Subramaniam S: Variance-modeled posterior inference of microarray data: detecting gene-expression changes in 3T3-L1 adipocytes. Bioinformatics. 2004, 20 (17): 3108-3127. 10.1093/bioinformatics/bth371

  84. 84.

    Iwamoto K, Kakiuchi C, Bundo M, Ikeda K, Kato T: Molecular characterization of bipolar disorder by comparing gene expression profiles of postmortem brains of major mental disorders. Molecular Psychiatry. 2004, 9 (4): 406-416. 10.1038/sj.mp.4001437

  85. 85.

    Matys V, Kel-Margoulis OV, Fricke E, Liebich I, Land S, Barre-Dirrie A, Reuter I, Chekmenev D, Krull M, Hornischer K: TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes. Nucleic Acids Res. 2006, D108-110. 34 Database

Download references

Acknowledgements

This research was supported in part by the Intramural Research Program of the NIH, National Institute of Environmental Health Sciences (J.G., R.T., C.P.) and National Institute of Aging (Y.Z., K.B.) as well as NIEHS extramural funding under grant 5R01ES014065-02 (M.R., A.D., C.M. and C.J.M.).

Author information

Author notes

    Affiliations

    1. Environmental Systems Biology Group, Laboratory of Molecular Toxicology, National Institute of Environmental Health Sciences, Research Triangle Park, NC, 27709, USA

      • Julia M Gohlke
      • , Reuben Thomas
      •  & Christopher J Portier

    2. Gene Expression and Genomics Unit, National Institute on Aging, National Institutes of Health, Baltimore, MD, 21224, USA

      • Yonqing Zhang
      •  & Kevin G Becker

    3. Department of Bioinformatics, Mount Desert Island Biological Laboratory, Old Bar Harbor Road, Salisbury Cove, ME, 04672, USA

      • Michael C Rosenstein
      • , Allan P Davis
      • , Cynthia Murphy
      •  & Carolyn J Mattingly

    Authors

    1. Search for Julia M Gohlke in:

    2. Search for Reuben Thomas in:

    3. Search for Yonqing Zhang in:

    4. Search for Michael C Rosenstein in:

    5. Search for Allan P Davis in:

    6. Search for Cynthia Murphy in:

    7. Search for Kevin G Becker in:

    8. Search for Carolyn J Mattingly in:

    9. Search for Christopher J Portier in:

    Corresponding author

    Correspondence to Christopher J Portier.

    Additional information

    Authors' contributions

    JG and RT designed and implemented the research with important suggestions from CP. YZ and KB provided interpretation and management of GAD database and MR, AD, CM and CJM provided interpretation, data management, and annotation of CTD database. All authors have read and approved the final manuscript.

    Julia M Gohlke, Reuben Thomas contributed equally to this work.

    Electronic supplementary material

    Additional file 1: Phenotype annotation of Genetic Association Database. The spreadsheet provides original phenotype names from the Genetic Association Database mapped to the annotated phenotype names used in the present analysis. (XLS 356 KB)

    Additional file 2: Annotation of nodes illustrated in Figure 3. The spreadsheet provides full chemical name annotation of the MeSH ID labels of nodes shown in Figure 3 as well as mapping to the 9 broad categories shown as colors in Figure 3. (XLS 163 KB)

    Additional file 3: Annotation of metabolic syndrome and neuropsychiatric subnetworks. The spreadsheet provides phenotypes and full chemical name annotation of the MeSH ID labels of nodes found within the 2 encircled subnetworks diagrammed in Figure 3. (XLS 29 KB)

    Additional file 4: DiRE results for gene expression datasets described in Table 1. Full output from the DiRE program showing enriched tfbs for each of the microarray datasets represented in Figure 4. (XLS 52 KB)

    Additional file 5: Comparison of results for Alzheimer's disease, Parkinson's disease and schizophrenia based on GAD or Genome Wide Association studies. Full sets of genes found to be associated with Alzheimer's disease, Parkinson's disease, or schizophrenia based either on the GAD databse or a database only representing results from Genome Wide Association studies. (XLS 34 KB)

    Additional file 6: Collapsed transcription factor binding site (tfbs) annotation based on similarity in matrices found in TRANSFAC. A file containing original TRANSFAC matrix names collapsed to annotated matrix names used in the current analysis. (XLS 42 KB)

    Authors’ original submitted files for images

    Rights and permissions

    Reprints and Permissions

    About this article

    Cite this article

    Gohlke, J.M., Thomas, R., Zhang, Y. et al. Genetic and environmental pathways to complex diseases. BMC Syst Biol 3, 46 (2009) doi:10.1186/1752-0509-3-46

    Download citation

    Keywords

    • Gene Expression Dataset
    • PPAR Signaling
    • Comparative Toxicogenomics Database
    • Tyrosine Metabolism
    • Familial Combine Hyperlipidemia

    BMC Systems Biology

    ISSN: 1752-0509

    Contact us