A genetic network model of cellular responses to lithium treatment and cocaine abuse in bipolar disorder
BMC Systems Biology volume 4, Article number: 158 (2010)
Lithium is an effective treatment for Bipolar Disorder (BD) and significantly reduces suicide risk, though the molecular basis of lithium's effectiveness is not well understood. We seek to improve our understanding of this effectiveness by posing hypotheses based on new experimental data as well as published data, testing these hypotheses in silico, and posing new hypotheses for validation in future studies. We initially hypothesized a gene-by-environment interaction where lithium, acting as an environmental influence, impacts signal transduction pathways leading to differential expression of genes important in the etiology of BD mania.
Using microarray and rt-QPCR assays, we identified candidate genes that are differentially expressed with lithium treatment. We used a systems biology approach to identify interactions among these candidate genes and develop a network of genes that interact with the differentially expressed candidates. Notably, we also identified cocaine as having a potential influence on the network, consistent with the observed high rate of comorbidity for BD and cocaine abuse. The resulting network represents a novel hypothesis on how multiple genetic influences on bipolar disorder are impacted by both lithium treatment and cocaine use. Testing this network for association with BD and related phenotypes, we find that it is significantly over-represented for genes that participate in signal transduction, consistent with our hypothesized-gene-by environment interaction. In addition, it models related pharmacogenomic, psychiatric, and chemical dependence phenotypes.
We offer a network model of gene-by-environment interaction associated with lithium's effectiveness in treating BD mania, as well as the observed high rate of comorbidity of BD and cocaine abuse. We identified drug targets within this network that represent immediate candidates for therapeutic drug testing. Posing novel hypotheses for validation in future work, we prioritized SNPs near genes in the network based on functional annotation. We also developed a "concept signature" for the genes in the network and identified additional candidate genes that may influence the system because they are significantly associated with the signature.
Bipolar Disorder (BD) is characterized by severe mood swings, from deep depression to mania, and shows familial transmission patterns consistent with multiple genetic influences on susceptibility . It poses a significant cost to affected individuals, including ~15% rate of suicide, and to society as a whole, affecting 1% to 3% of the population. Lithium is effective in preventing mania in many BD patients, although the molecular basis of lithium's action is not well understood and not all BD patients respond to lithium treatment [2, 3]. Since lithium treatment exerts an environmental influence on cells, and differential gene expression is one important mechanism of cellular response to environmental influences, we hypothesized that lithium activates signal transduction pathways leading to differential expression of genes related to BD mania [3–5].
Microarray analysis provides an unbiased approach to identifying genes that are differentially expressed in cells with treatment, relative to untreated cells . BD is thought to be a brain disorder, however, accessing the most appropriate tissue (live human brain) for expression studies is unlikely, so a model system is necessary. Some researchers have used brain tissue from animal models [7, 8], although expression patterns may be different in animals than in humans. Alternately, studies have been conducted with postmortem human brain tissue [8, 9], although expression in postmortem tissue may differ from expression in living tissue. As a third option, researchers use peripheral blood cells from humans to extrapolate genetic variants associated with brain disease [10–12]. Since expression in brain may be different from expression in peripheral blood cells, this third approach necessitates a follow-up analysis to maximize the likelihood that differential expression seen in peripheral cells is consistent with differential expression in brain. In spite of this extra step in the analysis, we chose peripheral blood cells because they are readily available from live human participants. Equally, while these cells are available, each sample represents a finite resource, so transformation to form Lymphoblast Cell Lines (LCLs) produces a resource that can be used in follow-on studies. EBV-transformed lines exhibit chromosomal stability  while providing ease of handling and availability for repeated DNA preparations for follow-on studies . Further, they are consistent with samples derived from repositories such as the Coriell Institute  or the Rutgers University Cell and DNA Repository , so comparisons can be made based on samples from a range of phenotypes. Notably, transformation could influence gene expression. However, the virus is incorporated randomly into the genome, so this influence would not be likely to produce consistent changes of expression of any particular genes. As such, while imperfect, we believe that LCLs represent the most appropriate model for this work.
In this work (Figure 1), based on our initial hypothesis of lithium activating signal transduction pathways important in BD etiology, we used microarray analysis to look for differential expression in untreated LCLs, relative to the same LCLs treated with a therapeutic dose of lithium. We confirmed that these candidate genes are expressed in brain and confirmed differential expression with rt-QPCR. We then looked for interactions among these candidate genes and used a systems biology approach to reveal a genetic network consistent with the hypothesized response to lithium treatment. This lithium response network is consistent with the neurotransmitter theory of BD, models responders and non-responders, and shows statistically significant over-representation of genes annotated for signal transduction. Interestingly, we also find a relationship with comorbid cocaine abuse, consistent with high rates of substance abuse in BD patients . Based on this network, we nominate: a) known therapeutic targets for drug testing, b) SNPs for validation testing, and c) additional candidate genes that share the "concept signature" of the network genes.
Genes Differentially Expressed with Lithium Treatment
We used 10 cell lines derived from BD patients who provided informed consent for their samples to be used in genetic studies. The protocol was approved by the University of Michigan's, IRBMED, Prechter Bipolar Genetics Repository, HUM00010454. All ten cell lines were derived from adult patients, ages 27 to 60, who were diagnosed as BD1 based on DSM III criteria. Of the ten, 5 were male and 5 female, 6 were responsive to lithium (ranging from 10 to 120 months of treatment) and 4 were naive to lithium treatment. No phenotypic information on cocaine use was collected on these patients.
Whole blood was drawn from each patient and white blood cells were isolated. Samples were transformed with EBV  and frozen for storage. Cell lines were thawed and cultured under identical conditions: RPMI-1640 media supplemented with 12% FBS, 2 mM L-glutamine, and 1% ampicillin/streptomycin, in an incubator set at 37°C with 5% constant CO2. For comparison, each line was divided and the treatment samples received 1 mM LiCl in the media, over a period of 8 days. We isolated total RNA from the stocked cell pellets for each cell line, using TRizol reagent according to the supplier's protocol (Invitrogen, Carlsbad, CA). We treated RNA samples with the RNase-free DNase kit (Cat. no. 79254, Qiagen), and further cleaned samples with the RNeasy MinElute Cleanup kit (Qiagen), according to the manufacturer's handbook. Total RNA quality was assessed by loading 1 μg of RNA on 1% agarose gel, to check for potential visible DNA contamination via gel electrophoresis, and by checking for A260/A280 ratio in the 1.8-2.0 range, via the Agilent® 2100 Bioanalyzer.
For cRNA synthesis, we started with 250 ng of total RNA from each sample and synthesized cRNA via the Illumina® TotalPrep™ RNA Amplification Kit, following protocols provided by the supplier, then determined cRNA quantities via the ND-1000 spectrophotometer. We performed BeadChip microarray hybridization according to the protocol provided by Illumina, using 20 Sentrix HumanRef-8 v2 BeadChip microarrays (Illumina, CA), 10 samples treated and 10 paired controls, loading 750 ng of cRNA onto each BeadChip. We used the BeadArray Reader 5000X (Illumina) to scan post-hybridization images, and BeadStudio software version 3 (Illumina) to process the scanned raw image data. After background subtraction and quantile normalization, we exported text files of the intensities of each probe.
We used the IBMT (Intensity-Based Moderated T-statistic)  method to assess differential gene expression. This method is an extension of the eBayes function in the limma R package , using an empirical Bayes test that models the dependence of variance on absolute expression levels, estimated based on local regression. The posterior variance is a weighted average of the sample variance and an intensity level-specific background variance level, with the weights determined by the empirical Bayes model. IBMT has been shown to perform favorably in experiments with small sample size. All parameters used by IBMT are internally estimated, based on the expression data, so the only user input to the algorithm is the relevant expression data. We guarded against uncontrolled false positive rates common in high-throughput experimental data  by selecting only genes that showed both statistically significant (Bonferroni corrected p-value ≤ 0.05) and biologically relevant (fold change ≥ +/-20%) differential expression. One of these transcripts is a hypothetical gene (FLJ39653) and one is a model (LOC400986), so we removed them from further consideration. We assessed brain expression for each candidate gene via Unigene's EST Profile database , then returned to the LCLs and used rt-QPCR to confirm differential expression in 6 genes randomly selected from the set. Briefly, we used reverse-transcription reactions to convert total RNA to first-strand cDNA with the SuperScript Preamplification System according to its manual (Cat No. 18089-011, Invitrogen). Synthesized first-strand cDNA was diluted 25-fold and used to template TaqMan assays from Applied Biosystems, Inc (ABI, CA). We carried out the TaqMan assays in triplicate, according to ABI's recommendations, using an SDS7900 real-time quantitative PCR thermocycler (ABI, CA) and used the 2-ΔΔCt method to calculate fold change, with normalization to the Ct values of the internal reference gene TBP. We performed two-tailed, paired, t-tests and tested for differential expression consistent with expression changes seen in the microarray study.
Candidate Gene Analysis
In complex diseases, multiple genetic influences converge on a single phenotype, consistent with some interaction among the genes involved. We first used PDG-ACE (Prioritizing Disease Genes by Analysis of Common Elements) [21, 22] to look for statistically significant gene-gene interactions among the differentially expressed candidate genes. Briefly, we used PDG-ACE's Medical Subject Heading (MeSH)- derived controlled vocabulary of 2,531 keywords, at least 106 iterations in significance testing, repeated each test to confirm a satisfactory survey of the genome, and only accepted keywords over-represented at a Bonferroni corrected p-value ≤ 0.05. After identifying a significant interaction between FOS and NR4A2 with respect to the keyword "cocaine", we used NESARC (the National Epidemiological Survey of Alcohol and Related Conditions) , a population based survey of 43,093 subjects, to assess rates of comorbid BD and cocaine abuse. We used GRAIL (Gene Relationships Across Implicated Loci)  to assess the group of candidate genes as a whole, and identify statistically over-represented keywords associated with the set of candidate genes. To assess the potential roles of each of the differentially expressed candidate genes, we used the MiSearch  adaptive publications search tool, querying specifically for association of each of the candidates with bipolar disorder, lithium treatment, and cocaine abuse.
Network Hypothesis Generation
While the differentially expressed genes may individually or collectively influence lithium's therapeutic action, they also interact with genes that are not differentially expressed with lithium treatment. To understand how the differentially expressed candidates interact with this larger set of genes, we used a systems biology approach to build a network using GeneGo's MetaCore  database of gene-gene and gene-small molecule interactions (GeneGo Inc., St. Joseph, MI). Assuming that the genes most likely to influence lithium response are those most closely interacting with the differentially expressed genes, we set parameters for the shortest path algorithm and accepted only interactions that were manually curated, including both functional and binding interactions. To take advantage of known metabolic and signaling pathways data we included canonical pathways information, then built the smallest network that includes all of the differentially expressed genes. To assess the most likely influences of environmental lithium, we added lithium to the network, along with the nodes (genes) required to include lithium in the network. After identifying cocaine comorbidity as a significant factor in the network, we added environmental cocaine, along with the nodes required to include cocaine in the network.
Network Hypothesis Testing
Based on our initial hypothesis, we tested this network for over-representation of genes associated with signal transduction. GeneGo automatically provides an internal test for association of a given network with a range of annotations, based on manually curated data that is proprietary to GeneGo. In addition, ConceptGen  is a recently developed open access resource that allows the user to upload a set of genes, then look for over-representation of genes associated with a range of "concepts". Casting a broad net, we searched for MeSH concepts that may be significantly associated with our network. Finally, we used the Genetic Association Database (GAD)  via the DAVID interface  to test whether the network was over-represented for genes associated with related disease phenotypes. We searched for "GENETIC_ASSOCIATION_DB_DISEASE_CLASS" and accepted phenotypes over-represented at the FDR ≤ 5% level. Since we added lithium and cocaine, plus the nodes required to connect them to the network, and in consideration of the possibility that these additional nodes would bias the analysis towards lithium and/or cocaine related phenotypes, we did a second round of GAD hypothesis testing, excluding these nodes.
Additional Hypotheses for Follow-up Validation and Testing
Given the initial set of differentially expressed genes, as well as the larger network based on interaction with these candidate genes, we sought to generate additional hypotheses for follow-up work. We first prioritized Single Nucleotide Polymorphisms (SNPs) in and near the genes in the network using the GIN  approach, to systematically annotate the appropriate SNPs with functional data and allow future researchers to prioritize them for follow-up. GeneGo provides annotation of genes that are either therapeutic or experimental targets for known drugs, so we used this annotation to prioritize drugs that target genes in the network, based on their known functions. Finally, we used a novel resource called ConSig (Concept Signature) to hypothesize additional candidate genes, based on the characteristics of the network genes . ConSig creates a "signature" for a set of genes, based on annotation of concepts associated with the individual genes, then searches the genome for additional genes that are significantly associated with the signature. Since these genes have signatures similar to the set of genes in the network, they represent novel hypotheses on association with the phenotype.
Differentially expressed genes
RNA preprocessing produced high quality total RNA, with an A260/A280 ratio in the range of 1.8-2.0 as measured by an Agilent® 2100 Bioanalyzer, while maximally maintaining the integrity of the RNA. A total of 22,177 transcripts were detected on the Sentrix HumanRef-8 v2 BeadChip. Table 1 shows the 12 transcripts showing differential expression meeting a Bonferroni corrected p-value ≤ 0.05 and fold change ≥ +/-20% (excluding hypothetical gene FLJ39653 and model gene LOC400986). The Unigene database shows that all of these genes are expressed in the human brain, though FOS and NR4A2 are expressed at significantly higher levels than the rest. As seen in Table 1, all 6 of the rt-QPCR results show the same direction and approximate magnitude of effect as seen in the microarray study, while 5 of the 6 show p-value ≤ 0.05. The rt-QPCR data on DDR2 is somewhat noisy with an insignificant p-value, however it was kept in the set of differentially expressed genes because it is suggested that IBMT is superior to a t-test in managing variance in small samples . Under the null hypothesis of no differential expression, the probability of 6 out of 6 tests matching for direction, and 5 of 6 matching for significance (p-value ≤ 0.05) is ~2.8 E-8.
PDG-ACE identified a significant interaction between FOS and FOSB, both members of the FBJ murine osteosarcoma viral oncogene family, based on the keyword "forebrain". In addition, FOS and NR4A2 show significant interaction based on the keyword "cocaine". Probing this interaction between FOS and NR4A2, NESARC shows an Odds Ratio of 5.6 for cocaine abuse among patients with mania in bipolar disorder, relative to the general population. For the set of 12 differentially expressed genes, GRAIL identified common keywords "induced, complement, phosphatase, induction, activation, response, neurons, expression, dopaminergic, mapk, transcription, kinase, cycle, mitogen, after, midbrain, rats, dopamine, mice, and activated". MiSearch (Table 2) revealed that four of our differentially expressed genes were previously associated with lithium in the literature (FOS, FOSB, NR4A2, RGS2) and six were previously associated with cocaine (DUSP1, FOS, FOSB, HIG2, NR4A2, and RGS2).
Network hypothesis generation
We started building the network looking for direct interactions among the differentially expressed genes, but the resulting network did not include all of the genes. We expanded the network to include one node between the differentially expressed genes, and the resulting network connects the 12 genes. We then added lithium and cocaine, as well as the nodes required to link them to the network. The resulting network includes the 12 genes differentially expressed with lithium treatment, the genes most closely interacting with them, and the nodes required to include lithium and cocaine (Figure 2). This network has been organized to display major hubs including c-FOS (FOS), MKP-1, NR4A2 (NURR1), FOSB, and RGS2, as well as the multiple positive and negative feedback loops in the network. The 12 differentially expressed genes have colored circles surrounding them. This network represents our secondary hypothesis on the mechanism through which the environmental influences of lithium treatment (upper left corner) and cocaine use (upper right corner) influence BD, lithium response, and/or comorbid cocaine abuse. Additional details of this network are included in Additional file 1 - GeneGo network details.
Network hypothesis testing
Based on internal data, GeneGo shows annotation for "signal transduction" as being significantly over-represented among the network genes, relative to all genes, with a p-value of 3.5E-29 (Additional file 1 - GeneGo network details). Based on MeSH annotation, ConceptGen analysis also shows signal transduction to be the MeSH term most significantly over-represented, with a False Discovery Rate of ~4 E-55. Both of these results are consistent with our original hypothesis. Testing this network for related hypotheses via the DAVID interface for the Genetic Association Database, we found that the network is over represented for pharmacogenomic, psychiatric, and chemical dependency phenotypes (Table 3), at the FDR ≤ 5.0% significance level. For the network that excludes lithium and cocaine, pharmacogenomic and psychiatric phenotypes remain significant, though the chemical dependency phenotype is not significant. Notably, the network is also associated with disease phenotypes that are not normally considered psychiatric or substance use disorders (ageing, cancer, etc.).
Hypotheses for follow-up analysis and testing
Using the GIN algorithm, we nominate ~5,000 SNPs for validation testing, and prioritize them based on functional annotation (Additional file 2 - GIN details). ConSig nominates additional candidate genes based on the signature of the network genes (Additional file 3 - Candidate genes nominated by Concept Signature). As expected, a consistent theme among concepts significant in this signature is signal transduction, consistent with the original hypothesis. Finally, we nominate drugs for follow-up testing because they are immediately available and target genes in the network (Additional file 1 - GeneGo Network details).
Based on our initial hypothesis that lithium treatment poses an environmental influence on cells through activation of signal transduction pathways, we investigated differential gene expression in response to lithium treatment. Microarray analysis revealed 12 genes that were significantly differentially expressed in LCLs with lithium treatment. The resulting list of candidate genes may provide insight to the etiology of lithium's effectiveness in BD. Seeking to put these genes into context, we first looked for interactions (commonality) among them via PDG-ACE and GRAIL analyses. PDG-ACE identified significant commonality between FOS and FOSB, serving as a positive control consistent with their roles as members of a single gene family. Notably, the keyword "forebrain", common across the FOS/FOSB pair, is also consistent with impulsivity seen in BD patients where variations in neurotransmitter signaling within forebrain regions may influence impulsivity associated with mania [32, 33]. In addition, PDG-ACE revealed a significant interaction between FOS and NR4A2 based on the keyword "cocaine". Pursuing this result via the NESARC survey, we found that the odds ratio for cocaine abuse among "manic" BD patients is more than 5 times that of the general population . This is also consistent with a number of studies that have documented association between BD and cocaine abuse [34–42]. GRAIL results are consistent with the psychiatric implications of the candidate genes but do not suggest novel hypotheses.
Based on our MiSearch publications search, we find that both FOS and FOSB are well established as candidate genes for both lithium response and cocaine abuse. Relatively little research currently associates NR4A2 with BD, lithium response, or cocaine abuse. However, Xing, et al., showed a reduction of NR4A2 in the prefrontal cortex of patients with BD  and Buervenich, et al., found that NR4A2 mutations caused a 30-40% reduction of in-vitro transcriptional activity in one case of BD . In both cases, the data are consistent with lithium having a therapeutic effect in BD, by normalizing deficient NR4A2 levels. Contrary to this effect, chronic lithium treatment was shown to decrease NR4A2 expression in rat brain , though this effect was localized to the CA1 hippocampal subregion. Much more research has been published on potential effects of FOS on BD susceptibility and/or lithium response. Rao, et al.,  recently reported increased expression of FOS in postmortem brain tissue from BD patients, relative to controls. St. Andre, et al., reported induction of FOS in multiple brain regions with LiCl treatment , along with Spencer and Houpt , Hammamura, et al., , Swank , and Portillo . These results are consistent with the multiple genetic and environmental interactions influencing lithium etiology in BD. However, while lithium appears to influence FOS expression, the multiple positive and negative feedback loops evident in the network model likely make simple predictions unreliable.
PDG-ACE points us in the direction of lithium's influence on dopamine signaling via NR4A2 and FOS, consistent with the catecholamine theory of BD etiology. NR4A2 also regulates dopaminergic neuron development [52, 53]. Cocaine represses NR4A2 expression [54–57] as well as DAT expression . Lithium may increase NR4A2 expression [45–58] and reduces mania, so we speculate that lithium may function in BD by increasing expression of DAT and other neurotransmitter related genes though, as with FOS, it is likely that simple predictions are unreliable. Interestingly, both lithium and cocaine pose environmental influences on cells, consistent with the hypothesized antagonism of these two substances [59–61], as well as the potential for using lithium treatment in cocaine abuse [62–64]. Notably, the keyword "forebrain" characterizing the interaction between FOS and FOSB is also consistent with cocaine abuse. Zahm, et al., showed that FOS expression in basal forebrain was "recalibrated" with cocaine use .
The network developed provides a model of the multiple interacting genetic and environmental influences involved in lithium etiology, as well as the influence of cocaine on this system. It is strongly consistent with our initial hypothesis that signal transduction plays an important role in lithium etiology, and also models related pharmacogenomic, psychiatric, and chemical dependence phenotypes. Variation in any of the genes in this network could influence an individual's response to lithium treatment or susceptibility to substance abuse, explaining the approximate 70% rate of lithium response in BD patients, as well as high rates of comorbid substance use disorders. Since substance use poses an environmental influence on cells, signal transduction is implicated along with neurotransmitter signaling and metabolism, consistent with our previous work in BD comorbid with tobacco use disorder , and depression comorbid with alcohol use disorders . In addition, as we observed in these previous studies, this network is enriched for genes associated with phenotypes that are not normally considered psychiatric disorders (ageing, cancer, immune disorders, etc) (Table3). This result challenges our traditional view of psychiatric and substance use disorders as being distinct from medical disorders.
Clinical implications for BD
Our network model reveals multiple pathways for both lithium and cocaine to influence the genetic network, as well as both positive and negative feedback loops. This is consistent with BD's characterization as a complex disease, where multiple genetic and environmental influences interact in predisposition to disease and in modulating sensitivity to drug therapies. In addition, individuals with certain constellations of genetic variants in the network genes may be more likely to respond to lithium treatment and/or be more vulnerable to comorbid cocaine abuse. Of the genes in the network, 28 are known therapeutic drug targets (Additional file 1 - GeneGo Network details).
Given the network of interactions, the genes associated, and the drugs known to target these genes, the potential exists for novel applications of known drugs in BD and cocaine abuse treatment. We know that psychosis may be present in mood episodes of bipolar patients, including both mania and depression . In addition, there is evidence  linking the use of anti-psychotic medications to treatment efficacy in bipolar depression and bipolar mania. Lithium appears to have an interaction with genes whose expression is also believed to be altered by treatment with certain antipsychotic medications (i.e. perphenazine, thioridazine, asenapine, chlorpromazine, and clozaril) . Clozaril, known for its potential side effect of neutropenia, is a medication purported to be effective in treatment resistant schizophrenia illness . The onset of such neutropenia is frequently managed successfully with the stoppage of clozaril treatment and the initiation of lithium treatment. This is consistent with lithium's known effect on white blood cell (WBC) elevation and subsequent neutrophilia, believed to be induced via lithium's induction of Granulocyte Macrophage-Colony Stimulating Factor (GM-CSF).
Sex hormone changes are widely believed to be associated with mood changes, both pathological and non-pathological. Indeed, the dysphoria often present in relation to the female menstrual cycle is well known to clinicians. In addition, pathological changes in the post-partum period are often present and particularly insidious in their severity and presentation. It would thus be logical to believe that medications that can alter sex hormone levels would influence mood. Indeed clomiphene and diethylstilbestrol treatment have both been observed to be associated with potential side effects of depression [72, 73] and anxiety [73, 74]. Abortifacients such as mifepristone have been investigated for their potential effect in the treatment of neurocognitive functioning and mood stability in bipolar disorder . Prednisone treatment is believed to be associated with the side effects of development of both depression and manic psychosis . Notably, both measured cortisol levels and results of the dexamethasone suppression test (DST) have been shown in bipolar subjects to differ from healthy controls , but are somewhat consistent across either depressive or manic episodes in a particular bipolar patient.
Among the other drugs that may be important in the treatment of BD, scopolamine has received some attention as an investigation treatment in the relief of bipolar depressive episodes. Treatment with TNF-alpha (Tumor Necrosis Factor alpha) has the frequent side effect of inducing depression and there have been reports of mania . As such, anti-TNF-alpha medications would potentially reverse this side effect. In addition, two of the genes that link cocaine to the network, PKC-epsillon and PYK2, are therapeutic targets of KAI1455 and PF562271, respectively, and may offer insight into the comorbidity of BD with cocaine abuse.
Based on the genes in our network, GIN analysis nominates ~5,000 SNPs for follow-up analysis, prioritized by functional annotation (Additional file 2 - GIN details). In addition, depending on the threshold set by the user, ConSig identifies hundreds of candidate genes that are significantly associated with the signature of the network genes (Additional file 3 - Genes nominated by Concept Signature). These genes represent novel candidates for validation in follow-up studies.
In this work, we hypothesized that signal transduction is significant in the etiology of lithium response in BD. Consistent with this hypothesis, we identified 12 genes that are differentially expressed with lithium treatment, then explored their likely roles in BD. In the course of this exploration, we identified a significant influence on comorbid cocaine abuse in BD, consistent with the epidemiological evidence. Using a systems biology approach to place the 12 differentially expressed genes into context, we developed a novel network model of the multiple interacting genetic and environmental influences on BD, lithium response, and comorbid cocaine abuse. Consistent with our hypothesis, this network is significantly associated with signal transduction, as well as pharmacogenomic, psychiatric, and chemical dependence phenotypes. Of the genes in this network, 28 are therapeutic drug targets, making them immediate candidates for follow-up drug testing. We also pose new hypotheses in the form of SNPs prioritized by functional annotation and novel candidate genes that are significantly associated with the concept signature of our network candidates.
Prioritizing Disease Genes by Analysis of Common Elements
Gene Relationships Across Implicated Loci
Genetic Association Database
Escamilla MA, Zavala JM: Genetics of bipolar disorder. Dialogues Clin Neurosci. 2008, 10: 141-152.
Storosum JG, Wohlfarth T, Schene A, Elferink A, van Zwieten BJ, van den Brink W: Magnitude of effect of lithium in short-term efficacy studies of moderate to severe manic episode. Bipolar Disord. 2007, 9: 793-798.
Alda M, Grof P, Rouleau GA, Turecki G, Young LT: Investigating responders to lithium prophylaxis as a strategy for mapping susceptibility genes for bipolar disorder. Prog Neuropsychopharmacol Biol Psychiatry. 2005, 29: 1038-1045.
Rowe MK, Chuang DM: Lithium neuroprotection: molecular mechanisms and clinical implications. Expert Rev Mol Med. 2004, 6: 1-18.
Ruzickova M, Turecki G, Alda M: Pharmacogenetics and mood stabilization in bipolar disorder. Am J Med Genet C Semin Med Genet. 2003, 123C: 18-25.
Taib Z: Statistical analysis of oligonucleotide microarray data. C R Biol. 2004, 327: 175-180.
Niculescu AB, Kelsoe JR: Convergent functional genomics: application to bipolar disorder. Ann Med. 2001, 33: 263-271.
Zhou X, Barrett TB, Kelsoe JR: Promoter variant in the GRK3 gene associated with bipolar disorder alters gene expression. Biol Psychiatry. 2008, 64: 104-110.
Kim S, Choi KH, Baykiz AF, Gershenfeld HK: Suicide candidate genes associated with bipolar disorder and schizophrenia: an exploratory gene expression profiling analysis of post-mortem prefrontal cortex. BMC Genomics. 2007, 8: 413-
Huff MO, Li XP, Ginns E, El-Mallakh RS: Effect of ethacrynic acid on the sodium- and potassium-activated adenosine triphosphatase activity and expression in Old Order Amish bipolar individuals. J Affect Disord. 2010, 123: 303-307.
Sugawara H, Iwamoto K, Bundo M, Ishiwata M, Ueda J, Kakiuchi C, Ishigooka J, Kato T: Effect of mood stabilizers on gene expression in lymphoblastoid cells. J Neural Transm. 2010, 117: 155-164.
Washizuka S, Iwamoto K, Kakiuchi C, Bundo M, Kato T: Expression of mitochondrial complex I subunit gene NDUFV2 in the lymphoblastoid cells derived from patients with bipolar disorder and schizophrenia. Neurosci Res. 2009, 63: 199-204.
Neitzel H: A routine method for the establishment of permanent growing lymphoblastoid cell lines. Hum Genet. 1986, 73: 320-326.
Editor in the field. The Coriell Institute for Medical Research. Nat Genet. 2005, 37: 2-
Schulze TG, Buervenich S, Badner JA, Steele CJ, Detera-Wadleigh SD, Dick D, Foroud T, Cox NJ, MacKinnon DF, Potash JB, et al.: Loci on chromosomes 6q and 6p interact to increase susceptibility to bipolar affective disorder in the national institute of mental health genetics initiative pedigrees. Biol Psychiatry. 2004, 56: 18-23.
Brady KT, Lydiard RB: Bipolar affective disorder and substance abuse. J Clin Psychopharmacol. 1992, 12: 17S-22S.
Sartor MA, Tomlinson CR, Wesselkamper SC, Sivaganesan S, Leikauf GD, Medvedovic M: Intensity-based hierarchical Bayes method improves testing for differentially expressed genes in microarray experiments. BMC Bioinformatics. 2006, 7: 538-
Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-
Blalock EM, Chen KC, Stromberg AJ, Norris CM, Kadish I, Kraner SD, Porter NM, Landfield PW: Harnessing the power of gene microarrays for the study of brain aging and Alzheimer's disease: statistical reliability and functional correlation. Ageing Res Rev. 2005, 4: 481-512.
Zhuo D, Zhao WD, Wright FA, Yang HY, Wang JP, Sears R, Baer T, Kwon DH, Gordon D, Gibbs S, et al.: Assembly, annotation, and integration of UNIGENE clusters into the human genome draft. Genome Res. 2001, 11: 904-918.
Keller BJ, McEachin RC: Identifying hypothetical genetic influences on complex disease phenotypes. BMC Bioinformatics. 2009, 10 (Suppl 2): S13-
McEachin RC, Keller BJ, Zandi PP, Almani M, McInnis MG: Prioritizing Disease Genes by Analysis of Common Elements (PDG-ACE). AMIA Annu Symp Proc. 2007, 1068-
Conway KP, Compton W, Stinson FS, Grant BF: Lifetime comorbidity of DSM-IV mood and anxiety disorders and specific drug use disorders results from the National Epidemiologic Survey on Alcohol and Related Conditions. J Clin Psychiatry. 2006, 67: 247-257.
Raychaudhuri S, Plenge RM, Rossin EJ, Ng AC, Purcell SM, Sklar P, Scolnick EM, Xavier RJ, Altshuler D, Daly MJ: Identifying relationships among genomic disease regions: predicting genes at pathogenic SNP associations and rare deletions. PLoS Genet. 2009, 5: e1000534-
States DJ, Ade AS, Wright ZC, Bookvich AV, Athey BD: MiSearch adaptive pubMed search tool. Bioinformatics. 2009, 25: 974-976.
Dezso Z, Nikolsky Y, Nikolskaya T, Miller J, Cherba D, Webb C, Bugrim A: Identifying disease-specific genes based on their topological significance in protein networks. BMC Syst Biol. 2009, 3: 36-
Sartor MA, Mahavisno V, Keshamouni VG, Cavalcoli J, Wright Z, Karnovsky A, Kuick R, Jagadish HV, Mirel B, Weymouth T, et al.: ConceptGen: a gene set enrichment and gene set relation mapping tool. Bioinformatics. 2010, 26: 456-463.
Becker KG, Barnes KC, Bright TJ, Wang SA: The genetic association database. Nat Genet. 2004, 36: 431-432.
Huang da W, Sherman BT, Tan Q, Kir J, Liu D, Bryant D, Guo Y, Stephens R, Baseler MW, Lane HC, Lempicki RA: DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res. 2007, 35: W169-175.
Saccone SF, Saccone NL, Swan GE, Madden PA, Goate AM, Rice JP, Bierut LJ: Systematic biological prioritization after a genome-wide association study: an application to nicotine dependence. Bioinformatics. 2008, 24: 1805-1811.
Wang XS, Prensner JR, Chen G, Cao Q, Han B, Dhanasekaran SM, Ponnala R, Cao X, Varambally S, Thomas DG, et al.: An integrative approach to reveal driver gene fusions from paired-end sequencing data in cancer. Nat Biotechnol. 2009, 27: 1005-1011.
Najt P, Perez J, Sanches M, Peluso MA, Glahn D, Soares JC: Impulsivity and bipolar disorder. Eur Neuropsychopharmacol. 2007, 17: 313-320.
Hellemans KG, Nobrega JN, Olmstead MC: Early environmental experience alters baseline and ethanol-induced cognitive impulsivity relationship to forebrain 5-HT1A receptor binding. Behav Brain Res. 2005, 159: 207-220.
Albanese MJ, Clodfelter RC, Pardo TB, Ghaemi SN: Underdiagnosis of bipolar disorder in men with substance use disorder. J Psychiatr Pract. 2006, 12: 124-127.
Bauer MS, Altshuler L, Evans DR, Beresford T, Williford WO, Hauger R: Prevalence and distinct correlates of anxiety substance, and combined comorbidity in a multi-site public sector sample with bipolar disorder. J Affect Disord. 2005, 85: 301-315.
Cassidy F, Ahearn EP, Carroll BJ: Substance abuse in bipolar disorder. Bipolar Disord. 2001, 3: 181-188.
Chengappa KN, Levine J, Gershon S, Kupfer DJ: Lifetime prevalence of substance or alcohol abuse and dependence among subjects with bipolar I and II disorders in a voluntary registry. Bipolar Disord. 2000, 2: 191-195.
Fenn HH, Bauer MS, Altshuler L, Evans DR, Williford WO, Kilbourne AM, Beresford TP, Kirk G, Stedman M, Fiore L: Medical comorbidity and health-related quality of life in bipolar disorder across the adult age span. J Affect Disord. 2005, 86: 47-60.
McElroy SL, Altshuler LL, Suppes T, Keck PE, Frye MA, Denicoff KD, Nolen WA, Kupka RW, Leverich GS, Rochussen JR, et al.: Axis I psychiatric comorbidity and its relationship to historical illness variables in 288 patients with bipolar disorder. Am J Psychiatry. 2001, 158: 420-426.
Pini S, Dell'Osso L, Mastrocinque C, Marcacci G, Papasogli A, Vignoli S, Pallanti S, Cassano G: Axis I comorbidity in bipolar disorder with psychotic features. Br J Psychiatry. 1999, 175: 467-471.
Sonne SC, Brady KT, Morton WA: Substance abuse and bipolar affective disorder. J Nerv Ment Dis. 1994, 182: 349-352.
Ringen PA, Melle I, Birkenaes AB, Engh JA, Faerden A, Jonsdottir H, Nesvag R, Vaskinn A, Friis S, Larsen F, et al.: Illicit drug use in patients with psychotic disorders compared with that in the general population a cross-sectional study. Acta Psychiatr Scand. 2008, 117: 133-138.
Xing G, Zhang L, Russell S, Post R: Reduction of dopamine-related transcription factors Nurr1 and NGFI-B in the prefrontal cortex in schizophrenia and bipolar disorders. Schizophr Res. 2006, 84: 36-56.
Buervenich S, Carmine A, Arvidsson M, Xiang F, Zhang Z, Sydow O, Jonsson EG, Sedvall GC, Leonard S, Ross RG, et al.: NURR1 mutations in cases of schizophrenia and manic-depressive disorder. Am J Med Genet. 2000, 96: 808-813.
Al Banchaabouchi M, Pena de Ortiz S, Menendez R, Ren K, Maldonado-Vlaar CS: Chronic lithium decreases Nurr1 expression in the rat brain and impairs spatial discrimination. Pharmacol Biochem Behav. 2004, 79: 607-621.
Rao JS, Harry GJ, Rapoport SI, Kim HW: Increased excitotoxicity and neuroinflammatory markers in postmortem frontal cortex from bipolar disorder patients. Mol Psychiatry. 2010, 15 (4): 384-92.
St Andre J, Albanos K, Reilly S: C-fos expression in the rat brain following lithium chloride-induced illness. Brain Res. 2007, 1135: 122-128.
Spencer CM, Houpt TA: Dynamics of c-fos and ICER mRNA expression in rat forebrain following lithium chloride injection. Brain Res Mol Brain Res. 2001, 93: 113-126.
Hamamura T, Lee Y, Ohashi K, Fujiwara Y, Miki M, Suzuki H, Kuroda S: A low dose of lithium chloride selectively induces Fos protein in the central nucleus of the amygdala of rat brain. Prog Neuropsychopharmacol Biol Psychiatry. 2000, 24: 285-294.
Portillo F, Carrasco M, Vallo JJ: c-fos expression in the rat hypothalamic paraventricular nucleus induced by LiCl: descending projections to the dorsal vagal motor nucleus. Rev Esp Fisiol. 1997, 53: 361-366.
Bunney WE: The current status of research in the catecholamine theories of affective disorders. Psychopharmacol Commun. 1975, 1: 599-609.
Perrone-Capano C, Di Porzio U: Genetic and epigenetic control of midbrain dopaminergic neuron development. Int J Dev Biol. 2000, 44: 679-687.
Wallen A, Perlmann T: Transcriptional control of dopamine neuron development. Ann N Y Acad Sci. 2003, 991: 48-60.
Bannon MJ, Pruetz B, Barfield E, Schmidt CJ: Transcription factors specifying dopamine phenotype are decreased in cocaine users. Neuroreport. 2004, 15: 401-404.
Bannon MJ, Pruetz B, Manning-Bog AB, Whitty CJ, Michelhaugh SK, Sacchetti P, Granneman JG, Mash DC, Schmidt CJ: Decreased expression of the transcription factor NURR1 in dopamine neurons of cocaine abusers. Proc Natl Acad Sci USA. 2002, 99: 6382-6385.
Lemberger T, Parkitna JR, Chai M, Schutz G, Engblom D: CREB has a context-dependent role in activity-regulated transcription and maintains neuronal cholesterol homeostasis. Faseb J. 2008, 22: 2872-2879.
Leo D, di Porzio U, Racagni G, Riva MA, Fumagalli F, Perrone-Capano C: Chronic cocaine administration modulates the expression of transcription factors involved in midbrain dopaminergic neuron function. Exp Neurol. 2007, 203: 472-480.
Misiuta IE, Saporta S, Sanberg PR, Zigova T, Willing AE: Influence of retinoic acid and lithium on proliferation and dopaminergic potential of human NT2 cells. J Neurosci Res. 2006, 83: 668-679.
Flemenbaum A: Antagonism of behavioral effects of cocaine by lithium. Pharmacol Biochem Behav. 1977, 7: 83-85.
Mandell AJ, Knapp S: Current research in the indoleamine hypothesis of affective disorders. Psychopharmacol Commun. 1975, 1: 587-597.
Mandell AJ, Knapp S: A neurobiological model for the symmetrical prophylactic action of lithium in bipolar affective disorder. Pharmakopsychiatr Neuropsychopharmakol. 1976, 9: 116-126.
Gawin FH, Kleber HD: Cocaine abuse treatment. Open pilot trial with desipramine and lithium carbonate. Arch Gen Psychiatry. 1984, 41: 903-909.
Kleber H, Gawin F: Psychopharmacological trials in cocaine abuse treatment. Am J Drug Alcohol Abuse. 1986, 12: 235-246.
Nunes EV, McGrath PJ, Wager S, Quitkin FM: Lithium treatment for cocaine abusers with bipolar spectrum disorders. Am J Psychiatry. 1990, 147: 655-657.
Zahm DS, Becker ML, Freiman AJ, Strauch S, Degarmo B, Geisler S, Meredith GE, Marinelli M: Fos after single and repeated self-administration of cocaine and saline in the rat: emphasis on the Basal forebrain and recalibration of expression. Neuropsychopharmacology. 2010, 35: 445-463.
McEachin RC, Saccone NL, Saccone SF, Kleyman-Smith YD, Kar T, Kare RK, Ade AS, Sartor MA, Cavalcoli JD, McInnis MG: Modeling complex genetic and environmental influences on comorbid bipolar disorder with tobacco use disorder. BMC Med Genet. 2010, 11: 14-
McEachin RC, Keller BJ, Saunders EF, McInnis MG: Modeling gene-by-environment interaction in comorbid depression with alcohol use disorders via an integrated bioinformatics approach. BioData Min. 2008, 1: 2-
Pregelj P: Psychosis and depression - a neurobiological view. Psychiatr Danub. 2009, 21 (Suppl 1): 102-105.
Hammonds MD, Shim SS: Effects of 4-week treatment with lithium and olanzapine on levels of brain-derived neurotrophic factor, B-cell CLL/lymphoma 2 and phosphorylated cyclic adenosine monophosphate response element-binding protein in the sub-regions of the hippocampus. Basic Clin Pharmacol Toxicol. 2009, 105: 113-119.
Mattai A, Fung L, Bakalar J, Overman G, Tossell J, Miller R, Rapoport J, Gogtay N: Adjunctive use of lithium carbonate for the management of neutropenia in clozapine-treated children. Hum Psychopharmacol. 2009, 24: 584-589.
Focosi D, Azzara A, Kast RE, Carulli G, Petrini M: Lithium and hematology: established and proposed uses. J Leukoc Biol. 2009, 85: 20-28.
O'Reilly EJ, Mirzaei F, Forman MR, Ascherio A: Diethylstilbestrol exposure in utero and depression in women. Am J Epidemiol. 2010, 171: 876-882.
Saunders EJ: Physical and psychological problems associated with exposure to diethylstilbestrol (DES). Hosp Community Psychiatry. 1988, 39: 73-77.
Choi SH, Shapiro H, Robinson GE, Irvine J, Neuman J, Rosen B, Murphy J, Stewart D: Psychological side-effects of clomiphene citrate and human menopausal gonadotrophin. J Psychosom Obstet Gynaecol. 2005, 26: 93-100.
Gallagher P, Malik N, Newham J, Young AH, Ferrier IN, Mackin P: Antiglucocorticoid treatments for mood disorders. Cochrane Database Syst Rev. 2008, CD005168-
Brown ES: Effects of glucocorticoids on mood, memory, and the hippocampus. Treatment and preventive therapy. Ann N Y Acad Sci. 2009, 1179: 41-55.
Schule C, Baghai TC, Eser D, Hafner S, Born C, Herrmann S, Rupprecht R: The combined dexamethasone/CRH Test (DEX/CRH test) and prediction of acute treatment response in major depression. PLoS One. 2009, 4: e4324-
Drevets WC, Furey ML: Replication of scopolamine's antidepressant efficacy in major depressive disorder: a randomized, placebo-controlled clinical trial. Biol Psychiatry. 2010, 67: 432-438.
Clerici M, Arosio B, Mundo E, Cattaneo E, Pozzoli S, Dell'osso B, Vergani C, Trabattoni D, Altamura AC: Cytokine polymorphisms in the pathophysiology of mood disorders. CNS Spectr. 2009, 14: 419-425.
This work was supported by National Institutes of Health grant #U54 DA021519 and the Prechter Bipolar Research Fund at the University of Michigan. SFS was supported by K01 DA024722. HC is supported by grants from the NIH (MH064596), the University of Michigan Comprehensive Depression Center Rachel Upjohn Clinical Scholars Award (2007), and Stanley Medical Research Institute's research grant (2005).
Financial: MGM has served as a consultant or on the speaker's bureau for Lilly, Pfizer, Merck, Astra-Zeneca, and Janssen Pharmaceuticals. SFS is listed as an inventor on a patent (US 20070258898) covering the use of certain SNPs in determining the diagnosis, prognosis, and treatment of addiction.
RCM and HC conceived the study and drafted the manuscript, with significant input from JDC and MGM. HC conducted the microarray experiments. MAS analyzed the microarray data and performed the concept signature analyses. RCM and BJK performed the systems biology analyses. SFS performed the SNP prioritization analyses and ARP assessed the clinical implications of drugs that target the network candidate genes. JDC and MGM provided expertise on systems biology analysis and bipolar disorder, respectively. All authors read and approved the final manuscript.
Richard C McEachin, Haiming Chen contributed equally to this work.
Electronic supplementary material
Additional file 1:GeneGo network details. GeneGo_Li_Resp_lithium_cocaine_network_statistics.xls is a spreadsheet detailing input parameters, output, drug targets, and hypothesis testing of the network. (XLS 1 MB)
Additional file 2:GIN details. SNPs_Prioritized_via_GIN.xls is a spreadsheet that prioritizes SNPs in and near the network genes based on functional annotation. (XLS 3 MB)
Additional file 3:Candidate genes nominated by Concept Signature. ConSig_results.xls is a spreadsheet detailing Concept Signature prioritization for novel candidate genes, based on their similarity to genes in the network. (XLS 6 MB)
About this article
Cite this article
McEachin, R.C., Chen, H., Sartor, M.A. et al. A genetic network model of cellular responses to lithium treatment and cocaine abuse in bipolar disorder. BMC Syst Biol 4, 158 (2010). https://doi.org/10.1186/1752-0509-4-158