This article has Open Peer Review reports available.
Using gene expression signatures to identify novel treatment strategies in gulf war illness
© Craddock et al. 2015
Received: 26 February 2015
Accepted: 26 June 2015
Published: 9 July 2015
Gulf War Illness (GWI) is a complex multi-symptom disorder that affects up to one in three veterans of this 1991 conflict and for which no effective treatment has been found. Discovering novel treatment strategies for such a complex chronic illness is extremely expensive, carries a high probability of failure and a lengthy cycle time. Repurposing Food and Drug Administration approved drugs offers a cost-effective solution with a significantly abbreviated timeline.
Here, we explore drug re-purposing opportunities in GWI by combining systems biology and bioinformatics techniques with pharmacogenomic information to find overlapping elements in gene expression linking GWI to successfully treated diseases. Gene modules were defined based on cellular function and their activation estimated from the differential expression of each module’s constituent genes. These gene modules were then cross-referenced with drug atlas and pharmacogenomic databases to identify agents currently used successfully for treatment in other diseases. To explore the clinical use of these drugs in illnesses similar to GWI we compared gene expression patterns in modules that were significantly expressed in GWI with expression patterns in those same modules in other illnesses.
We found 19 functional modules with significantly altered gene expression patterns in GWI. Within these modules, 45 genes were documented drug targets. Illnesses with highly correlated gene expression patterns overlapping considerably with GWI were found in 18 of the disease conditions studied. Brain, muscular and autoimmune disorders composed the bulk of these.
Of the associated drugs, immunosuppressants currently used in treating rheumatoid arthritis, and hormone based therapies were identified as the best available candidates for treating GWI symptoms.
KeywordsGulf war illness Systems biology Bioinformatics Drug repurposing Pharmacogenomics Complex chronic illness
Novel drug discovery is a costly, high-risk and extremely time-consuming enterprise. Costing in the range of millions to billions of dollars [1–4], the process can take up to 15 years to complete with one of every two novel drug candidates failing in the later stages . As such, the repurposing of Food and Drug Administration (FDA) approved drugs offers a fast cost effective solution for discovering novel treatments. This is all the more true for conditions with no known treatments.
Gulf War Illness (GWI) is but one example of a complex chronic illness with no known cure and which requires long-term treatment and monitoring. This both prolongs patient suffering and increases the financial burden of the illness, for the individual, the family and society. GWI is a multi-symptom disorder exhibiting a complex constellation of symptoms that include fatigue, musculoskeletal pain, and cognitive dysfunction . Since returning from the first Gulf War over 20 years ago investments of nearly $1 billion have been made to Gulf War Veterans health , however there is still no effective treatment available for the nearly 250,000 veterans with GWI. Though the cause and illness mechanisms of GWI are largely unknown, a leading hypothesis points to the involvement of neuroinflammatory cascade possibly triggered by exposure to battlefield toxins and exacerbated by stress [7–9]. Thus, this illness has far-reaching consequences and the development of effective treatments promises to benefit not only this patient population but a host of others as well.
In an attempt to identify candidate treatment avenues for GWI involving currently approved pharmacological agents we used a combination of systems biology and bioinformatics techniques combined with pharmacogenomic information to compare gene expression patterns in GWI to known drug targets and expression patterns found in a set of human diseases. We identified a specific set of functional gene modules with altered expression in GWI pointing to the cellular processes affected in this disorder. Significant correlations between gene expression profiles in these gene modules were found between GWI and other disease conditions. Cross-referencing of disease and drug information in pharmacogenomic databases with druggable genes in these functional modules ultimately revealed novel drug repurposing avenues for the treatment of GWI.
As part of a larger ongoing study a subset of GWI male subjects (n = 17) and healthy but sedentary Gulf War era veterans (n = 22) were recruited from the Miami Veterans Administration Medical Centers, clinics and the local veteran community between April 2006 and May 2008. All subjects were comparable in age, body mass index (BMI), and ethnicity. Subjects were male and ranged in age between 30 and 55. Inclusion criteria was derived from Fukuda et al. , and consisted in identifying veterans deployed to the theater of operations between August 8, 1990 and July 31, 1991, with one or more symptoms present after 6 months from at least 2 of the following: fatigue; mood and cognitive complaints; and musculoskeletal complaints. Subjects were in good health prior to 1990, and had no current exclusionary diagnoses . Medications that could have impacted immune function were excluded. Use of the Fukuda definition in GWI is supported by Collins et al. . Additional details may be found in Broderick et al. .
All subjects signed an informed consent approved by the Institutional Review Board of the University of Miami. Ethics review and approval for data analysis was also obtained by the IRB of the University of Alberta.
Blood was drawn at rest at comparable times of day from each subject during the April 2006 to May 2008 period. At each blood draw three 8-mL tubes of blood were collected in CPT vacutainers (B-D- Biosciences, San Jose, CA). The peripheral blood mononuclear cells (PBMC) were isolated and stored in liquid nitrogen under conditions designed to maintain viability. Specifically, whole blood was added to Ficoll-Paque, centrifuged at 1000 g for 25 min. PBMC’s were isolated from the PBMC ring atop the Ficoll layer into a separate tube, centrifuged at 300 g for 10 min, then re-suspended in PBS. Cells were then counted using a Beckman Coulter viCell, and cryopreserved in freezing media (temperature lowered 1 oC per minute until -80 oC).
Total RNA was extracted using TRI Reagent (Molecular Research Center, Cincinnati, OH) following the manufacturer’s protocol. The quality and quantity of RNA was assessed using the Agilent Bioanalyzer 2100 RNA 6000 Nano Kit (Agilent Technologies, CA). From each sample, 300 ng of total RNA was converted into cDNA by reverse transcription using a T7-oligo(dT) primer and the Affymetrix 3′ IVT Express Kit (Affymetrix, Santa Clara, CA) according to standard manufacturer protocol. The generated cDNA was purified using the GeneChip Sample Cleaning Module (Affymetrix) and labeled cRNA was generated by in vitro transcription using the biotinylated nucleotide mix. This was then purified with the Cleaning Module and quantified using the Nanodrop ND-1000 spectrophotometer (NanoDrop Technologies, Inc., Wilmington, DE USA). In each preparation 11 μg cRNA was fragmented in Fragmentation Buffer (Affymetrix) in a final reaction volume of 25 μl.
Hybridization, washing, staining and scanning were done using Affymetrix GeneChip instruments (Hybridization Oven 640, Fluidics Station 450Dx, Scanner GCS3000Dx) and Affymetrix Human U133 2.0 arrays (Affymetrix) as per manufacturer’s standards. Microarray image files (.cel data) were generated using the Affymetrix GCOS software tool with default microarray analysis parameters to provide overall within chip normalization of the image intensity distribution. The quality parameters that were monitored besides cRNA total yield and cRNA A260/A280 ratio included: (i) background noise (Q value), (ii) percentage of present called probe sets, (iii) scaling factor, (iv) information about exogenous Bacillus subtilis control transcripts from the Affymetrix Poly-A control kit (lys, phe, thr, and dap), and (v) the ratio of intensities of 3′ probes to 5′ probes for a housekeeping gene (GAPDH).
To generate a broad comparison group the Gene Expression Omnibus (GEO) DataSets [14, 15] and data from Suthram et al.  were used to obtain a set of gene expression profiles describing a number of human disease conditions. We restricted our selection to include only those sets in which both disease and a corresponding healthy control group were measured in the same cell type or tissue in the same experimental conditions. Sets that included different exposure times, exposure concentrations or multiple cell/tissue types were each treated as a separate disease condition. For consistency, and to avoid complications arising due to cross platform comparisons, datasets were restricted to the Affymetrix Gene Chip Human Genome U133A, U133 Plus 2.0, U133A 2.0 and U95 Version 2 arrays, to align with our GWI gene expression data. All diseases affecting male subjects in the GEO database meeting these criteria were included in this study. Overall, this resulted in 101 human disease gene expression datasets (Additional file 1: Table S1). The disease profiles selected provide a comprehensive set of diseases and include various cancers, neurodegenerative diseases, autoimmune illnesses, chemical exposures, viral infections, neurological and neuromuscular disorders.
Data transformation and normalization
Gene expression data was Log2 transformed then normalized using a Z-score transformation for each microarray sample to allow for the direct comparison of values across various samples and diseases.
Gene functional modules
This is a modification of the module activity described in  were the signed T-score was used. When averaging over signed T-scores equal but opposite scores will tend to cancel, resulting in no module activity, when in fact there are genes significantly over or under expressed. Here the absolute value was used in the differential module activity to capture both up and down-expressed genes. A threshold of MA = 1.5, corresponding to a p-value of ~0.1, was chosen as a liberal cutoff value to identify differentially expressed gene modules. Only modules with MA values above threshold were taken to be differentially expressed in GWI and considered for further analysis, all others were discarded.
Individual gene expression analysis
The 202 genes in the 19 modules with MA values above threshold were individually compared between GWI and controls. Individual genes i were compared between GWI subjects and healthy controls using an unpaired t-test. To account for multiple comparisons false discovery rates (FDR) were then calculated for each comparison from these resulting p-values using the procedure introduced by Storey . Genes with FDR of less than or equal to 0.05 were taken to be significantly different in GWI compared to control. Fold change was calculated by calculating the ratio of average gene expression between GWI and health controls for non-Log2 transformed data. Positive values were taken to indicate the fold-increase while ratios that were less than 1 were inverted and given a negative value to denote the fold-decrease. An absolute fold change of 1.5 was taken as a cutoff .
Mining the pharmacogenomic knowledge base (PharmGKB) database
Genes from modules affected in GWI were screened against the PharmGKB database (8.1.2015)  to find gene-drug and gene-disease relationships supported with pharmacogenomics research reported in the literature.
Pathway-based functional annotation of gene modules
Functional annotation of gene modules with gene-drug relations was performed using the ConsensusPathDB [20–22] to provide biological pathway information for each gene set. Over-representation analysis  incorporating the Kyoto Encyclopedia of Genes and Genomes (KEGG) (73.0) , Netpath (1.1.2015) , the Small Molecule Pathway Database (SMPDB) (8.1.2015) , the Integrating Network Objects with Hierarchies (INOH) (1.1.2015) , Biocarta (2009_05_12) , Humancyc (18.5) , Signalink (8.1.2015) , Edinburgh human metabolic network (Ehmn) (1.1.2015) , Reactome (51) , PharmGKB (8.1.2015) , Wikipathways (9.1.2015)  and the Pathway Interaction Database (PID) (2014_02_14)  pathway sets was used to interpret the function of druggable gene modules. Here the significance of the observed overlap between the gene module and the members of known pathways, compared to random expectations, was calculated based on the hypergeometric distribution. A minimum overlap of 2 genes between the gene module and the pathway set at a p-value cutoff 0.01 was required. Specifically, the p-value is calculated as the probability of randomly finding k or more successes from the population in N total draws. Thus, small p-values indicate a greater overlap than expected by chance. Pathway sets containing the majority of the druggable genes, the highest number of module genes overall, and the lowest p-value were taken as the functional annotation of the module. Pathway annotation was performed only to provide biological pathway information for each gene module set. All subsequent analysis was performed on these gene module sets constructed on the bases of human protein-protein interactions  and not on the bases of known pathway membership.
Gene module alignment across illnesses
Significance of correlation
A random background distribution of disease correlations across each module was created to determine if module correlations between GWI and other diseases were significantly different from random chance. Gene to module assignments were first randomized, while preserving the total number of modules, the number of genes per module, and the number of modules to which each gene belongs as in Suthram et al. . Then disease and control datasets for the disease conditions were randomly shuffled before calculating partial correlations with intact GWI datasets. This process was performed 100 times for each module to generate the random background distributions. The cumulative distribution of the random background correlations per module was used to calculate the p-value of the true correlations in each module. FDR was calculated for each correlation from these resulting p-values using the procedure introduced by Storey . Correlation values with FDR of less than or equal to 0.05 were taken to be significantly different from random background. All others were removed.
Affected genes and modules
Drug targetable affected genes
Pathway representation in druggable gene modules
Of these 8 drug targetable GWI modules, 86.7 % of the genes in module 2 belonged to the Tumor Necrosis Factor (TNF) alpha signaling pathway as annotated in the Netpath database (1.1.2015)  (Additional file 3: Table S2). This pathway includes both drug targetable genes, TNF receptor-associated factor 1 (TRAF1) and the inhibitor of kappa light polypeptide gene enhancer in B-cells, kinase gamma (IKGBG). Overlapping with this module is module 14, which embodies much of the Immune System pathway as annotated in the Reactome database (51) . Overall, 73.3 % of the genes in module 14 are associated with this pathway including the caspase 9, apoptosis-related cysteine peptidase (CASP9) gene, and the CD40 molecule, TNF receptor superfamily member 5 (CD40) gene.
Modules 5, 8, 9, 17, and 18 overlap to form a cluster of targetable genes. Module 5 was annotated as the PID (2014_02_14)  pathway for the Validated Nuclear Estrogen Receptor Alpha Network. While only 21.4 % of the genes in module 5 are included in this pathway, two of the three drug targets identified are among these, namely the breast cancer 1, early onset (BRCA1), and estrogen receptor 1 (ESR1) genes. The third druggable target, the X-ray repair cross-complementing protein 1 (XRCC1) gene, however was not included in the Validated Nuclear Estrogen Receptor Alpha Network. The Netpath (1.1.2015)  Androgen Receptor (AR) pathway was representative of the genes present in module 8. A total of 62.4 % of the genes in this module were found in the AR pathway, including 6 of the 7 drug targetable genes: BRCA1, cyclin D1 (CCND1), epidermal growth factor receptor (EGFR), nuclear receptor coactivator 1 (NCOA1) catenin (cadherin-associated protein), beta 1, 88 kDa (CTNNB1), and signal transducer and activator of transcription 3 (acute-phase response factor) (STAT3). Only ESR1 was not affiliated with the AR pathway. The PID (2014_02_14)  pathway for Signaling Mediated by p38-alpha & p38 beta was assigned as the annotation for module 9. Both drug targetable genes, mitogen activated protein kinases 11 and 14, are included among the 42.9 % of module 9 genes found in this pathway. Genes in both modules 17 and 18 were identified as part of the PID (2014_02_14)  ATF-2 Transcription Factor Network. This pathway contains 57.1 % and 66.7 % of the genes in module 17 and 18, respectively, including all of the identified drug targets for these modules, namely MAPK 1, 8, 9 and 14, and ESR1.
Finally, module 13 was enriched in genes from the Reactome database (51)  Metabolism of Lipids & Lipoproteins pathway. Only 21.4 % of the genes in the module are included in this pathway, and the coiled-coil alpha-helical rod protein 1 (CCHCR1) drug targetable gene is not among them making this isolated pathway uninformative for the selection of treatment repurposing in GWI.
Similarity of GWI to known human diseases
Refinement of drug targets through illness similarity
GWI is a complex multi-symptom illness, which manifests with cognitive dysfunction, fatigue, and musculoskeletal pain . At this time it is unclear whether GWI represents several syndromes, or one syndrome with several subtypes. While stratification of the case defined GWI subjects into subtypes based on shared clusters of symptoms is an active area of research which may provide improved resolution on underlying gene functional involvement, this issue will only be resolved when objective markers are firmly established. The analysis presented here is a step towards this goal.
Cognitive difficulties, such as memory problems, constitute one of the most common unexplained impairments reported in GWI . Recent work examining brain changes in GWI support an underlying neurobiological underpinning to these problems. Results show that participants with GWI perform significantly slower and less accurately on working memory tasks than matched healthy veterans, and that this decrease correlates with lower levels of activity in prefrontal brain regions . Fatigue and pain, other commonly reported symptoms associated with GWI have been linked to alterations in the brain’s white matter in GWI subjects [36, 37]. Memory problems, as well as mood disorder, have also been linked to mild inflammation and degenerations in the hippocampus in a mouse model of GWI . Our findings, indicate that over a third (38.9 %) of the diseases aligning with greater than 50 % of the GWI affected modules are localized to the central nervous system.
Musculoskeletal pain is another frequently reported symptom of GWI . A higher proportion of veterans of the Persian Gulf War of 1991 reported symptoms of muscle and joint pain than a military comparison group . This chronic pain has been linked to abnormal central processing of sensory and painful stimuli in GWI subjects in brain areas normally serving sensory perception, and threat and arousal, as well as in thalamocortical circuits and cerebellum . GWI subjects suffering from chronic musculoskeletal pain also experience greater naturally occurring muscle pain during exercise compared to healthy veterans and become more sensitive to pain stimuli following acute exercise suggesting that acute exercise augments the central nervous system sensitivity to sensory information . In this work we find that more than 1 in 5 (22 %) of the diseases, which overlap in gene expression with GWI, are muscular disorders.
Evidence is mounting that there is also a significant immune component to GWI. Ongoing Th1-type immune activation, as measured by intracellular production of cytokines in peripheral blood, appears to be symptomatic in afflicted Gulf War Veterans when compared to healthy counterparts . More recently, this finding has been confirmed while also suggesting that this may occur in the more complex context of a mixed Th1:Th2 response . In addition to this, our group previously reported altered gene expression associated with NK cell function and decreased NK cell cytotoxicity in GWI subjects . Abnormal regulation of the immune system can result in autoimmune diseases , and allergies . Like GWI, these illnesses are characterized by abnormal resting levels of immune cells, cytokines and circulating hormones. Our results indicate that close to 1 in 6 (17 %) of the overlapping diseases are autoimmune disorders.
While gene expression in GWI was measured in PBMC, gene expression in some of the other human illnesses was obtained from brain tissue (e.g. hippocampus in FTLD, substantia nigra in Parkinson’s). Although these measures of pathway activation were expressed in different tissue, substantial alignment between gene expression in blood and brain biopsies has been shown in Parkinson’s, Alzheimer’s and Huntington’s [48–50] supporting the relevance of this profiling. When considered collectively, these findings support an illness model whereby GWI might be described primarily as a central nervous system disorder manifesting with musculoskeletal problems that is potentially fueled by a dysregulation of immunity. Central to this would be an ongoing neuroinflammatory process.
A common issue in pathway analysis is that genes shared between overlapping functional modules and pathways may cause an inflated or biased p-value of statistical significance for some annotations. In such situations highly influential genes that are shared across multiple pathways may bias the identification of some pathways unfairly. The current analysis is not immune to this problem. As such it is pertinent to examine the results in the context of the specific illness in question, and the individual genes highlighted by the analysis. Pathway annotation of the drug treatable GWI affected modules revealed two major pathway clusters. The first is dominated by the Netpath (1.1.2015)  annotated pathway for Tumor Necrosis Factor (TNF) alpha and the Reactome (51)  Immune System pathway. Overall, this cluster of pathways suggest that immune signaling, particularly TNF-α signaling, are prime targets for repurposed drug treatment in GWI. Projecting this onto known illnesses with similar pathway involvement pointed to commonalities with RA. The relation between GWI and RA in the Netpath (1.1.2015) : Tumor Necrosis Factor (TNF) alpha pathway identifies the gene encoding TRAF1 in conjunction with several TNF-α inhibitors. Suppressed activity of TNF receptor and apoptotic pathways , and higher responsiveness of TNF-α [43, 45] have previously been identified in GWI subjects by our group, consistent with the current findings, suggesting an autoimmune component to this illness. TNF-α blockers, such as infliximab, adalimumab and etanercept, are currently approved by the FDA for the treatment of chronic inflammatory diseases such as RA, Crohn’s disease and ankylosing spondylitis. While, such treatments may affect immune function, inducing autoantibodies leading to conditions such as drug-induced lupus, there is low risk of infliximab, etanercept, and adalimumab in association with drug-induced lupus despite the higher incidence of autoantibodies with its use . As chronic inflammation is a suggested component of GWI pathophysiology , these compounds may provide a promising novel treatment avenue for GWI, however further tests are required.
The second pathway cluster central to GWI that was identified here is comprised of the Netpath (1.1.2015)  Androgen Receptor (AR) pathway and PID (2014_02_14)  Validated Nuclear Estrogen Receptor Alpha Network as well as PID (2014_02_14)  pathways for the ATF-2 Transcription Factor Network, and Signaling Mediated by p38-alpha and p38-beta. This large cluster is targetable by the majority of identified drug classes including immunosuppressants, protein kinase inhibitors, hormone treatments, estrogen receptor antagonists, monoclonal antibodies, taxanes, anti-neoplastics and platinum compounds. The p-38 MAPKs participate in signaling cascades that control cellular responses to immune signals and stress, while ATF-2 is normally activated in response to signals that converge on the stress-activated p38 MAPKs. Both the PID (2014_02_14)  pathway Validated Nuclear Estrogen Receptor Alpha Network pathway and Netpath (1.1.2015) : Androgen Receptor (AR) pathway highlight sex steroid hormone pathways. As androgens are the precursor of all estrogens, and both estrogen  and androgen  receptor activity is mediated by p38 signaling, this pathway cluster suggests that sex steroid signaling, as mediated by p38, may also be prime targets for repurposed drug treatment in GWI.
Estrogen is also a potent, albeit complex, modulator of inflammation . The estrogen receptor 1 (ESR1) gene is identified in relation to RA in both the PID database (2014_02_14)  Validated Nuclear Estrogen Receptor Alpha Network, and ATF-2 Transcription Factor Network pathways. Leflunomide, an immunosuppressive disease-modifying anti-rheumatic drug whose efficacy is modulated by ESR1 gene polymorphisms , is the only drug identified with a direct use in active moderate to severe rheumatoid arthritis. The remainder of the drug candidates identified in association with ESR1 and RA are most commonly used in the treatment and prevention of estrogen receptor positive cancers. Cisplatin, a platinum-containing anti-cancer drug that binds DNA resulting in apoptosis, also has immunosuppressive effects and has been shown to be effective in managing RA in a recent single case study . Alkylating agents, drugs that modulate the immune system such as cyclophosphamide, are sometimes used to treat severe cases of rheumatoid arthritis . The remaining candidate drugs associated with RA and ESR1 directly modulate estrogen signaling. Tamoxifen and fulvestrant, are selective estrogen receptor modulators (SERMs), that act as antagonists in breast and agonists in other tissue. Tamoxifen has been shown to exert anti-inflammatory effects in human RA cell cultures [58, 59]. Additionally, other SERMs currently show promise in treating RA in mouse models of both genders [60, 61]. However, estrogen receptor antagonists have also been reported to induce RA like symptoms in female subjects  (female mice), . Likewise exemestane, an irreversible steroidal aromatase inactivator, has been suggested to induce or reveal RA in female subjects [64–66]. The discrepancy in effects may be attributable to the influence of gender on tamoxifen-induced biochemical changes [67, 68].
Though based on gene expression in peripheral blood, these results are nonetheless consistent with the involvement of a neuroinflammatory component in GWI . Estrogen plays a significant role in inflammation in general , and neuroinflammation specifically , mediating via astrocytes neuroprotective effects on clinical function, inflammation in the central nervous system, and axonal loss . Furthermore, estrogen receptor α moderates the cytokine- and chemokine-mediated neuroinflammatory response . SERMs, including tamoxifen and fulvestrant, appear to have an important anti-inflammatory role  and are thus a promising treatment avenue for GWI. However, as we have found GWI subjects under maximal exercise challenge exhibit increased activation of pathways involving NF-κB  the tendency of SERMs to increase activity of the nuclear NFκB complex  must be taken into consideration. Evidence also suggests a role for estrogen and estrogen receptors in Amyotrophic Lateral Sclerosis (ALS). The incidence and prevalence of ALS are greater in men than in women  and evidence supports the notion that endogenous female hormones have a neuroprotective effect on motor neurons [75, 76]. This is of import as Gulf War veterans have been shown to develop amyotrophic lateral sclerosis (ALS) at twice the rate of non-deployed era veterans [77, 78]. Our analysis is consistent with this showing a GA value of 49 % between ALS and GWI. Early phase II trials of tamoxifen in ALS have demonstrated preliminary efficacy [79, 80], and therefore may be of therapeutic benefit for ALS like symptoms associated with GWI. Additionally, estrogen receptor agonists and estrogen have been shown to attenuate TNF-α induced apoptosis in motor neurons [75, 81]. While SERMs and TNF-α inhibitors provide novel treatment avenues for GWI, further analysis must be performed.
Collectively this work supports a treatment strategy that would target both immune and sex hormone signaling in GWI. Interestingly this aligns with our recent studies of altered neuroendocrine-immune homeostasis in GWI using computational models of known immune, stress (hypothalamic-pituitary-adrenal) and sex (hypothalamic-pituitary-gonadal) hormone signaling [82, 83]. In both these studies we found the cross-talk between immune and sex steroid signaling could in principle support alternate neuroendocrine-immune regulatory modes that resemble chronic GWI-associated signatures. Specifically we found evidence in this analysis that GWI may be associated with dysregulation in genetic pathways involved in immune signaling, particularly TNF-α, signaling, and sex steroid signaling, as mediated by p38 MAPKs. When compared to a set of human disease conditions GWI most closely resembles brain disorders with musculoskeletal problems and dysregulation of immunity, consistent with previously observed reports. Overall, while specific drugs have been highlighted the general findings of our study suggest the use of TNF-α immunosuppressive agents and SERMs, either individually or in conjunction, in the treatment of GWI. Further analysis of these treatment avenues will ultimately reveal their efficacy.
This work was funded under U.S. Department of Defense (CDMRP program) grant W81XWH-10-1-0774 (G. Broderick, PI) and grant W81XWH-09-2-0071 (N. Klimas, PI), as well as the U.S. Department of Veterans Affairs (Merit Award, N. Klimas, PI). This research was conducted in collaboration with and using the resources of the University of Miami Center for Computational Science (CCS). The authors would like to extend special thanks to Dr. Joel Zysman, Director of High Performance Computing at the CCS, for his support and many helpful discussions.
- Paul SM, Mytelka DS, Dunwiddie CT, Persinger CC, Munos BH, Lindborg SR, et al. How to improve R&D productivity: the pharmaceutical industry’s grand challenge. Nat Rev Drug Discov. 2010;9:203–14.PubMedGoogle Scholar
- Adams C, Brantner V. Estimating the cost of new drug development: is it really 802 million dollars? Health Aff (Millwood). 2006;25(2):420–8.Google Scholar
- Adams CP, Brantner VV. Spending on new drug development. Health Econ. 2010;19(2):130–41.PubMedGoogle Scholar
- DiMasi JA, Hansen RW, Grabowski HG. The price of innovation: new estimates of drug development costs. J Health Econ. 2003;22(2):151–85.PubMedGoogle Scholar
- Wang Y. Extracting knowledge from failed development programmes. Pharm Med. 2012;26(2):91–6.Google Scholar
- RAC-GWVI (Research Advisory Committee on Gulf War Veterans’ Illnesses). Gulf War Illness and the Health of Gulf War Veterans. Washington, DC: U.S. Government Printing Office; 2008.Google Scholar
- Amourette C, Lamproglou I, Barbier L, Fauquette W, Zoppe A, Viret R, et al. Gulf War illness: effects of repeated stress and pyridostigmine treatment on blood–brain barrier permeability and cholinesterase activity in rat brain. Behav Brain Res. 2009;203(2):207–14.PubMedGoogle Scholar
- Lamproglou I, Barbier L, Diserbo M, Fauvelle F, Fauquette W, Amourette C. Repeated stress in combination with pyridostigmine: part I: long-term behavioural consequences. Behav Brain Res. 2009;197(2):301–10.PubMedGoogle Scholar
- Barbier L, Diserbo M, Lamproglou I, Amourette C, Peinnequin A, Fauquette W. Repeated stress in combination with pyridostigmine: part II: changes in cerebral gene expression. Behav Brain Res. 2009;197(2):292–300.PubMedGoogle Scholar
- Fukuda K, Nisenbaum R, Stewart G, Thompson WW, Robin L, Washko RM, et al. Chronic multisymptom illness affecting Air Force veterans of the Gulf War. J Am Med Assoc. 1998;280:981–98.Google Scholar
- Reeves WC, Lloyd A, Vernon SD, Klimas N, Jason LA, et al. Identification of ambiguities in the 1994 chronic fatigue syndrome research case definition and recommendations for resolution. BMC Health Serv Res. 2003;3:25.PubMedPubMed CentralGoogle Scholar
- Collins JF, Donta ST, Engel CC, Baseman JB, Dever LL, Taylor T, et al. The antibiotic treatment trial of Gulf War Veterans’ Illnesses: issues, design, screening, and baseline characteristics. Control Clin Trials. 2002;23:333–53.PubMedGoogle Scholar
- Broderick G, Ben-Hamo R, Vashishtha S, Efroni S, Nathanson L, Barnes Z, et al. Altered immune pathway activity under exercise challenge in Gulf War Illness: an exploratory analysis. Brain Behav Immun. 2013;28:159–69.PubMedGoogle Scholar
- Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets--update. Nucleic Acids Res. 2013;41(Database issue):D991–5.PubMedGoogle Scholar
- Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–10.PubMedPubMed CentralGoogle Scholar
- Suthram S, Dudley JT, Chiang AP, Chen R, Hastie TJ, Butte AJ. Network-based elucidation of human disease similarities reveals common functional modules enriched for pluripotent drug targets. PLoS Comput Biol. 2010;6(2):e1000662.PubMedPubMed CentralGoogle Scholar
- Storey JD. A direct approach to false discovery rates. J R Stat Soc B. 2002;64(3):479–98.Google Scholar
- Dalman MR, Deeter A, Nimishakavi G, Duan ZH. Fold change and p-value cutoffs significantly alter microarray interpretations. BMC Bioinformatics. 2012;13 Suppl 2:S11.PubMedPubMed CentralGoogle Scholar
- Whirl-Carrillo M, McDonagh EM, Hebert JM, Gong L, Sangkuhl K, Thorn CF, et al. Pharmacogenomics knowledge for personalized medicine. Clinic Pharmacol Therapeut. 2012;92(4):414–7.Google Scholar
- Kamburov A, Wierling C, Lehrach H, Herwig R. ConsensusPathDB–a database for integrating human interaction networks. Nucleic Acids Res. 2009;37 suppl 1:D623–8.PubMedGoogle Scholar
- Kamburov A, Pentchev K, Galicka H, Wierling C, Lehrach H, Herwig R. ConsensusPathDB: toward a more complete picture of cell biology. Nucleic Acids Res. 2011;39 suppl 1:D712–7.PubMedGoogle Scholar
- Kamburov A, Stelzl U, Lehrach H, Herwig R. The ConsensusPathDB interaction database: 2013 update. Nucleic Acids Res. 2013;41(D1):D793–800.PubMedGoogle Scholar
- Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40:D109–14.PubMedGoogle Scholar
- Kandasamy K, Mohan SS, Raju R, Keerthikumar S, Kumar GSS, Venugopal AK, et al. NetPath: a public resource of curated signal transduction pathways. Genome Biol. 2010;11:R3.PubMedPubMed CentralGoogle Scholar
- Frolkis A, Knox C, Lim E, Jewison T, Law V, Hau DD, et al. SMPDB: the small molecule pathway database. Nucleic Acids Res. 2010;38:D480–7.PubMedGoogle Scholar
- Yamamoto S, Sakai N, Nakamura H, Fukagawa H, Fukuda K, Takagi T. INOH: ontology-based highly structured database of signal transduction pathways. Database. 2011;2011:bar052.PubMedPubMed CentralGoogle Scholar
- Nishimura D. Biocarta. Biotech Software & Internet Report. 2001;2(3):117–20.Google Scholar
- Romero P, Wagg J, Green ML, Kaiser D, Krummenacker M, Karp PD. Computational prediction of human metabolic pathways from the complete human genome. Genome Biol. 2004;6:R2.PubMedPubMed CentralGoogle Scholar
- Fazekas D, Koltai M, Türei D, Módos D, Pálfy M, Dúl Z, et al. SignaLink 2.0 - a signaling pathway resource with multi-layered regulatory networks. BMC Syst Biol. 2013;7:7.PubMedPubMed CentralGoogle Scholar
- Hao T, Ma HW, Zhao XM, Goryanin I. Compartmentalization of the Edinburgh human metabolic network. BMC Bioinformatics. 2010;11:393.PubMedPubMed CentralGoogle Scholar
- Croft D, Mundo AF, Haw R, Milacic M, Weiser J, Wu G, et al. The Reactome pathway knowledgebase. Nucleic Acids Res. 2014;42(Database issue):D472–7.PubMedGoogle Scholar
- Kelder T, van Iersel MP, Hanspers K, Kutmon M, Conklin BR, Evelo CT, et al. WikiPathways: building research communities on biological pathways. Nucleic Acids Res. 2012;40:D1301–7.PubMedGoogle Scholar
- Schaefer CF, Anthony K, Krupa S, Buchoff J, Day M, Hannay T, et al. PID: the pathway interaction database. Nucleic Acids Res. 2009;37:D674–9.PubMedGoogle Scholar
- Institute of Medicine. Health Consequences of Service During the Persian Gulf War: Recommendations for Research and Information Systems. Washington DC: National Academies Press; 1996.Google Scholar
- Hubbard NA, Hutchison JL, Motes MA, Shokri-Kojori E, Bennett IJ, Brigante RM, et al. Central executive dysfunction and deferred prefrontal processing in veterans with gulf war illness. Clin Psychol Sci. 2014;2(3):319–27.PubMedPubMed CentralGoogle Scholar
- Rayhan RU, Stevens BW, Raksit MP, Ripple JA, Timbol CR, Adewuyi O, et al. Exercise challenge in Gulf War Illness reveals two subgroups with altered brain structure and function. PLoS ONE. 2013;8(6):e63903.PubMedPubMed CentralGoogle Scholar
- Rayhan RU, Stevens BW, Timbol CR, Adewuyi O, Walitt B, Van Meter JW, et al. Increased brain white matter axial diffusivity associated with fatigue, pain and hyperalgesia in Gulf War illness. PLoS ONE. 2013;8(3):e58493.PubMedPubMed CentralGoogle Scholar
- Parihar VK, Hattiangady B, Shuai B, Shetty AK. Mood and memory deficits in a model of gulf war illness are linked with reduced neurogenesis, partial neuron loss, and mild inflammation in the hippocampus. Neuropsychopharmacol. 2013;38(12):2348–62.Google Scholar
- Thomas HV, Stimpson NJ, Weightman A, Dunstan F, Lewis G. Pain in veterans of the Gulf War of 1991: a systematic review. BMC Musculoskel Dis. 2006;7(1):74.Google Scholar
- Gopinath K, Gandhi P, Goyal A, Jiang L, Fang Y, Ouyang L, et al. FMRI reveals abnormal central processing of sensory and pain stimuli in ill Gulf War veterans. Neurotoxicol. 2012;33(3):261–71.Google Scholar
- Cook DB, Stegner AJ, Ellingson LD. Exercise alters pain sensitivity in gulf war veterans with chronic musculoskeletal pain. J Pain. 2010;11(8):764–72.PubMedGoogle Scholar
- Skowera A, Cleare A, Blair D, Bevis L, Wessely SC, Peakman M. High levels of type 2 cytokine-producing cells in chronic fatigue syndrome. Clin Exp Immunol. 2004;135(2):294–302.PubMedPubMed CentralGoogle Scholar
- Broderick G, Kreitz A, Fuite J, Fletcher MA, Vernon SD, Klimas N. A pilot study of immune network remodeling under challenge in Gulf War Illness. Brain Behav Immun. 2011;25:302–13.PubMedGoogle Scholar
- Whistler T, Fletcher MA, Lonergan W, Zeng XR, Lin JM, Laperriere A, et al. Impaired immune function in Gulf War illness. BMC Med Genomics. 2009;2:12.PubMedPubMed CentralGoogle Scholar
- Warnatz K, Voll RE. Pathogenesis of autoimmunity in common variable immunodeficiency. Frontiers Immunol. 2012;3:210.Google Scholar
- Simpson CR, Anderson WJA, Helms PJ, Taylor MW, Watson L, Prescott GJ, et al. Coincidence of immune-mediated diseases driven by Th1 and Th2 subsets suggests a common aetiology. A population-based study using computerized General Practice data. Clin Exp Allergy. 2002;32(1):37–42.PubMedGoogle Scholar
- Chang C, Gershwin ME. Drugs and autoimmunity – a contemporary review and mechanistic approach. J Autoimmun. 2012;34:J266–75.Google Scholar
- Mutez E, Larvor L, Leprêtre F, Mouroux V, Hamalek D, Kerckaert JP, et al. Transcriptional profile of Parkinson blood mononuclear cells with LRRK2 mutation. Neurobiol Aging. 2011;32(10):1839–48.PubMedGoogle Scholar
- Maes OC, Xu S, Yu B, Chertkow HM, Wang E, Schipper HM. Transcriptional profiling of Alzheimer blood mononuclear cells by microarray. Neurobiol Aging. 2007;28(12):1795–809.PubMedGoogle Scholar
- Runne H, Kuhn A, Wild EJ, Pratyaksha W, Kristiansen M, Isaacs JD, et al. Analysis of potential transcriptomic biomarkers for Huntington’s disease in peripheral blood. Proc Natl Acad Sci U S A. 2007;104(36):14424–9.PubMedPubMed CentralGoogle Scholar
- Johnson G, Slater B, Leis L, Bach R. Elevated plasma biomarkers of chronic inflammation in Gulf War illness (591.3). FASEB J. 2014;28(1 Supplement):591–3.Google Scholar
- Bhatt S, Xiao Z, Meng Z, Katzenellenbogen BS. Phosphorylation by p38 mitogen-activated protein kinase promotes estrogen receptor α turnover and functional activity via the SCFSkp2 proteasomal complex. Mol Cell Biol. 2012;32(10):1928–43.PubMedPubMed CentralGoogle Scholar
- Gioeli D, Black BE, Gordon V, Spencer A, Kesler CT, Eblen ST, et al. Stress kinase signaling regulates androgen receptor phosphorylation, transcription, and localization. Mol Endocrinol. 2006;20(3):503–15.PubMedGoogle Scholar
- Straub RH. The complex role of estrogens in inflammation. Endocrine Rev. 2007;28(5):521–74.Google Scholar
- Dziedziejko V, Kurzawski M, Safranow K, Chlubek D, Pawlik A. The effect of ESR1 and ESR2 gene polymorphisms on the outcome of rheumatoid arthritis treatment with leflunomide. Pharmacogenom. 2011;12(1):41–7.Google Scholar
- Omar EK, Adimulam S. Does the key to treat Rheumatoid nodules lie with oncology?-Is Cisplatin an option? BMC Musculoskelet Disord. 2013;14 Suppl 1:A5.PubMed CentralGoogle Scholar
- Turesson C. Extra-articular rheumatoid arthritis. Curr Opin Rheumatol. 2013;25(3):360–6.PubMedGoogle Scholar
- Mitani M, Miura Y, Saura R, Kitagawa A, Fukuyama T, Hashiramoto A, et al. Estrogen specifically stimulates expression and production of osteoprotegerin from rheumatoid synovial fibroblasts. Int J Mol Med. 2005;15(5):827–32.PubMedGoogle Scholar
- Komi J, Möttönen M, Luukkainen R, Lassila O. Non-steroidal anti-oestrogens inhibit the differentiation of synovial macrophages into dendritic cells. Rheumatol. 2001;40(2):185–91.Google Scholar
- Jochems C, Islander U, Erlandsson M, Engdahl C, Lagerquist M, Gjertsson I, et al. Role of endogenous and exogenous female sex hormones in arthritis and osteoporosis development in B10. Q-ncf1*/* mice with collagen-induced chronic arthritis. BMC Musculoskelet Disord. 2010;11(1):284.PubMedPubMed CentralGoogle Scholar
- Jordan VC. Tamoxifen: catalyst for the change to targeted therapy. Eur J Cancer. 2008;44(1):30–8.PubMedPubMed CentralGoogle Scholar
- Jansson L, Holmdahl R. Enhancement of collagen-induced arthritis in female mice by estrogen receptor blockage. Arthritis Rheum. 2001;44(9):2168–75.PubMedGoogle Scholar
- Creamer P, Lim K, George E, Dieppe P. Acute inflammatory polyarthritis in association with tamoxifen. Br J Rheumatol. 1994;33(6):583–5.PubMedGoogle Scholar
- Bertolini E, Letho-Gyselinck H, Prati C, Wendling D. Rheumatoid arthritis and aromatase inhibitors. Joint Bone Spine. 2011;78(1):62–4.PubMedGoogle Scholar
- Bruzzese V, Hassan C, Zullo A, Zampa G. Rheumatoid arthritis: a complication of aromatase inhibitor therapy? Int J Immunopathol Pharmacol. 2010;24(4):1099–101.Google Scholar
- Morel B, Marotte H, Miossec P. Will steroidal aromatase inhibitors induce rheumatoid arthritis? Ann Rheumatic Dis. 2007;66(4):557–8.Google Scholar
- Hemieda FAKES. Influence of gender on tamoxifen-induced biochemical changes in serum of rats. Mol Cell Biochem. 2007;301(1-2):137–42.PubMedGoogle Scholar
- Dluzen DE, Mickley KR. Gender differences in modulatory effects of tamoxifen upon the nigrostriatal dopaminergic system. Pharmacol Biochem Behav. 2005;80(1):27–33.PubMedGoogle Scholar
- Osório J. Reproductive endocrinology: less estrogen, more neuroinflammation? Nat Rev Endocrinol. 2012;8:381.PubMedGoogle Scholar
- Spence RD, Hamby ME, Umeda E, Itoh N, Du S, Wisdom AJ, et al. Neuroprotection mediated through estrogen receptor-α in astrocytes. Proc Natl Acad Sci U S A. 2011;108(21):8867–72.PubMedPubMed CentralGoogle Scholar
- Brown CM, Mulcahey TA, Filipek NC, Wise PM. Production of proinflammatory cytokines and chemokines during neuroinflammation: novel roles for estrogen receptors α and β. Endocrinol. 2010;151:4916–25.Google Scholar
- Rubio N, Cerciat M, Unkila M, Garcia-Segura LM, Arevalo MA. An in vitro experimental model of neuroinflammation: the induction of interleukin-6 in murine astrocytes infected with Theiler’s murine encephalomyelitis virus, and its inhibition by oestrogenic receptor modulators. Immunol. 2011;133(3):360–9.Google Scholar
- Giraud SN, Caron CM, Pham-Dinh D, Kitabgi P, Nicot AB. Estradiol inhibits ongoing autoimmune neuroinflammation and NFκB-dependent CCL2 expression in reactive astrocytes. Proc Natl Acad Sci U S A. 2010;107(18):8416–21.PubMedPubMed CentralGoogle Scholar
- McCombe PA, Henderson RD. Effects of gender in amyotrophic lateral sclerosis. Gender Medicine. 2010;7(6):557–70.PubMedGoogle Scholar
- Cardona-Rossinyol A, Mir M, Caraballo-Miralles V, Lladó J, Olmos G. Neuroprotective effects of estradiol on motoneurons in a model of rat spinal cord embryonic explants. Cell Mol Neurobiol. 2013;33(3):421–32.PubMedGoogle Scholar
- de Jong S, Huisman M, Sutedja N, van der Kooi A, de Visser M, Schelhaas J, et al. Endogenous female reproductive hormones and the risk of amyotrophic lateral sclerosis. J Neurol. 2013;260(2):507–12.PubMedGoogle Scholar
- Haley RW. Excess incidence of ALS in young Gulf War veterans. Neurology. 2003;61(6):750–6.PubMedGoogle Scholar
- Horner RD, Kamins KG, Feussner JR, Grambow SC, Hoff-Lindquist J, Harati Y, et al. Occurrence of amyotrophic lateral sclerosis among Gulf War veterans. Neurology. 2003;61:742–9.PubMedGoogle Scholar
- Traynor BJ, Bruijn L, Conwit R, Beal F, O’Neill G, Fagan SC, et al. Neuroprotective agents for clinical trials in ALS: a systematic assessment. Neurology. 2006;67:20–7.PubMedGoogle Scholar
- Brooks B, Sanjak M, Roelke K. Phase 2B randomized dose-ranging clinical trial of tamoxifen, a selective estrogen receptor modulator [SERM], in amyotrophic lateral sclerosis (ALS): sensitivity analyses of discordance between survival and functional outcomes with long-term follow-up. Amyotroph Lateral Scler. 2005;6 Suppl 1:118.Google Scholar
- Das A, Smith JA, Gibson C, Varma AK, Ray SK, Banik NL. Estrogen receptor agonists and estrogen attenuate TNF-α-induced apoptosis in VSC4.1 motoneurons. J Endocrinol. 2011;208(2):171–82.PubMedGoogle Scholar
- Craddock TJA, Fritsch P, Rice Jr MA, del Rosario RM, Miller DB, Fletcher MA, et al. A role for homeostatic drive in the perpetuation of complex chronic illness: Gulf War illness and chronic fatigue syndrome. PLoS ONE. 2014;9(1):e84839.PubMedPubMed CentralGoogle Scholar
- Fritsch P, Craddock TJA, del Rosario RM, Rice Jr MA, Smylie AL, Folcik VA, et al. Succumbing to the laws of attraction: exploring the sometimes pathogenic versatility of discrete immune logic. Sys Biomed. 2014;1(3):179–91.Google Scholar
- Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;2009(25):1091–3.Google Scholar
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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.