|Home | About | Journals | Submit | Contact Us | Français|
Anxiety disorders are prevalent and disabling yet understudied from a genetic standpoint, compared with other major psychiatric disorders such as bipolar disorder and schizophrenia. The fact that they are more common, diverse and perceived as embedded in normal life may explain this relative oversight. In addition, as for other psychiatric disorders, there are technical challenges related to the identification and validation of candidate genes and peripheral biomarkers. Human studies, particularly genetic ones, are susceptible to the issue of being underpowered, because of genetic heterogeneity, the effect of variable environmental exposure on gene expression, and difficulty of accrual of large, well phenotyped cohorts. Animal model gene expression studies, in a genetically homogeneous and experimentally tractable setting, can avoid artifacts and provide sensitivity of detection. Subsequent translational integration of the animal model datasets with human genetic and gene expression datasets can ensure cross-validatory power and specificity for illness. We have used a pharmacogenomic mouse model (involving treatments with an anxiogenic drug—yohimbine, and an anti-anxiety drug—diazepam) as a discovery engine for identification of anxiety candidate genes as well as potential blood biomarkers. Gene expression changes in key brain regions for anxiety (prefrontal cortex, amygdala and hippocampus) and blood were analyzed using a convergent functional genomics (CFG) approach, which integrates our new data with published human and animal model data, as a translational strategy of cross-matching and prioritizing findings. Our work identifies top candidate genes (such as FOS, GABBR1, NR4A2, DRD1, ADORA2A, QKI, RGS2, PTGDS, HSPA1B, DYNLL2, CCKBR and DBP), brain–blood biomarkers (such as FOS, QKI and HSPA1B), pathways (such as cAMP signaling) and mechanisms for anxiety disorders—notably signal transduction and reactivity to environment, with a prominent role for the hippocampus. Overall, this work complements our previous similar work (on bipolar mood disorders and schizophrenia) conducted over the last decade. It concludes our programmatic first pass mapping of the genomic landscape of the triad of major psychiatric disorder domains using CFG, and permitted us to uncover the significant genetic overlap between anxiety and these other major psychiatric disorders, notably the under-appreciated overlap with schizophrenia. PDE10A, TAC1 and other genes uncovered by our work provide a molecular basis for the frequently observed clinical co-morbidity and interdependence between anxiety and other major psychiatric disorders, and suggest schizo-anxiety as a possible new nosological domain.
‘Worry is a thin stream of fear trickling through the mind. If encouraged, it cuts a channel into which all other thoughts are drained.'
–Arthur Somers Roche
Anxiety disorders are prevalent and disabling. Approximately 30 million people are affected with anxiety disorders in United States1, 2 and the 12-month prevalence rate is estimated to be 18.1%.3 Anxiety disorders, under DSM classification, include generalized anxiety disorder (GAD), panic disorder, specific phobias, post-traumatic stress disorder (PTSD) and obsessive-compulsive disorder (OCD). They can be grouped into those without an obvious external trigger (GAD, panic disorder), those with an obvious external trigger (PTSD, phobias) and those that are more of a mixed picture, like OCD. Anxiety disorders are often co-morbid with other psychiatric disorders such as depression, bipolar disorder, schizophrenia and substance abuse.4, 5 Phenomenologicaly, anxiety disorders seem to have in common an increased reactivity to the environment, driven by uncertainty and fear of perceived threats.6 Stress is a common trigger and/or exacerbator.
Despite their prevalence and clinical impact, anxiety disorders are understudied from a genetic standpoint, compared with other major psychiatric disorders. Twin, adoption and familial studies have suggested a role for heritability in anxiety disorders.7, 8 Human genetic linkage studies have identified some susceptibility loci.9, 10, 11 Genetic association studies have identified polymorphisms in genes such as corticotropin-releasing hormone (CRH),12 glutamate transporter (SLC1A1),13 adenosine A2a receptor (ADORA2A),14 regulator of G-protein signaling 2 (RGS2),15, 16 delta-aminolevulinate dehydratase (ALAD),17 dynein light chain 2 (DYNLL2)17 and others as possibly involved in anxiety disorders, but with limited reproducibility. There are few published human postmortem brain gene expression studies to date on anxiety and related disorders.18, 19
To overcome this suboptimal state of affairs, we employed a comprehensive convergent functional genomics (CFG)20, 21, 22, 23 approach as a way of identifying and prioritizing candidate genes and blood biomarkers for anxiety disorders, as we did in our previous work on bipolar disorder,24, 25, 26, 27, 28 schizophrenia29, 30 and alcoholism.31 As a first step, we used drug effects on gene expression in mice in key brain regions for anxiety (prefrontal cortex (PFC), amygdala (AMY) and hippocampus (HIP)),32 as well as blood (BLD), as a way to tag genes that may have pathophysiological relevance. We then cross-matched and integrated that gene-level data with multiple other lines of evidence (genetic and gene expression) from human studies and other animal model studies (Figure 1).
For our mouse brain and blood gene expression studies, we used an agonist drug, which induces symptoms of anxiety (yohimbine),33, 34, 35 and a gold standard antagonist drug, which is used to treat anxiety disorders (diazepam)36, 37 (Figure 1). From the range of doses of the drugs that had been reported in the literature to have our desired behavioral effects, we chose doses at the low end of the range, to minimize potential supraphysiological dosing artifacts and side-effects. We also employed a behavioral readout to make sure the drugs were absorbed and doing what they were supposed to do (Figure 2).
Changes in gene expression in response to each of the two drugs, yohimbine and diazepam, would be of interest in and of themselves, in terms of candidate gene generation and CFG. However, not all genes that show changes in expression in response to either of the drugs are necessarily germane to the pathophysiology of anxiety and related disorders. It is likely that some of the gene expression changes have to do with other effects of the drugs, particularly their individual side-effects. We reasoned, first, that genes that change in expression in response to both drugs are more likely to be involved in the core pathophysiology we are modeling, and are higher probability candidate genes. Second, co-treatment with the two drugs, one an anxiogenic, and the other one an anxiolytic, could arguably show interference effects, and some of the genes that would be changed by single drug treatment would be ‘nipped in the bud' and show no changes in expression in response to co-treatment. Those genes would also be deemed higher probability candidate genes than the genes that still change during co-treatment.
As external cross-validators, for each gene changed in expression in our pharmacogenomics studies, we used six independent lines of evidence in our CFG analyses (Figure 1d). First, we assessed if there was any published genetic evidence—human genetic evidence of association with anxiety, or at least if it mapped to a linkage locus that had been implicated in anxiety disorders. We also looked at mouse transgenic or quantitative trait loci (QTL) studies relevant to anxiety. Second, we assessed if there was any published gene expression evidence in brain or blood in anxiety disorders, from human studies and, more broadly, from other animal models of anxiety.38 These external lines of evidence suffer from the obvious drawback of being constrained by what has been published so far, limiting novelty, and to the inherent biases and limitations of those particular lines of work.
According to Bayesian theory, an optimal estimate results from combining previous information with new evidence. Although we cannot exclude that some of the candidate genes we have identified are false positives because of potential biological or technical limitations of the methodology and approach we employed, the higher the number of independent lines of evidence (i.e. the higher the CFG score), the lower the likelihood of that being the case. The CFG scoring is arguably a reasonable compromise between specificity and sensitivity, between focus and broadness.
Our approach identifies and prioritizes an extensive series of candidate genes, some of which have already been reported using various related treatments or paradigms, as well as many others which are novel. Moreover, the coalescence of the candidate genes into pathways and mechanisms is of particular importance and opens new directions. Finally, we compared our results with our previous similar work in bipolar disorder,25, 26 schizophrenia29 and alcoholism,31 and were able to analyze the significant genetic overlap between anxiety and these other disorders, providing a molecular basis for the frequently observed clinical co-morbidity.
All experiments were performed with male C57/BL6 mice, 8–12 weeks of age, obtained from Jackson Laboratories (Bar Harbor, ME, USA), and acclimated for at least 2 weeks in our animal facility (Indiana University School of Medicine LARC) on reverse light cycle (1000 to 2200 hours) before any experimental manipulation. All experiments were conducted at the same time of day—between 1400 and 1600 hours. Mice were treated by intraperitoneal injection with single-dose of yohimbine (1mgkg–1), diazepam (0.3mgkg–1), a combination of yohimbine and diazepam (1 and 0.3mgkg–1), or control (vehicle) solution only. The control solution, which was also used to dissolve the drugs, consisted of 0.325% Tween 80 in 0.9% phosphate-buffered saline and alcohol (EtOH) at a final concentration of 10μlml–1 EtOH.
A SMART II Video Tracker system (San Diego Instruments, San Diego, CA, USA) was used to track movement of mice under normal light immediately after drug administration. After injection, mice were placed in the lower right-hand corner of one of four adjacent, 41_41_34-cm3 enclosures. Mice had no physical contact with other mice during testing. Each enclosure has nine pre-defined areas, that is, center area, corner areas and wall areas. After an initial 15min of adaptation, measures of locomotor activity were obtained from the second half (15min) of the total 30-min time recorded immediately after injection of the drugs, with a focus on behavior in the open field center area.
Three independent de novo biological experiments, performed at different times, were used for gene expression studies. Each experiment consisted of three mice per treatment condition, for a total of nine mice per condition across the three experiments (Figure 1b). Brain and blood from the same de novo experiment were used for microarray studies.
Twenty-four hours after drug administration, mice were sacrificed by cervical dislocation. The brains of the mice were harvested, stereotactically sliced and hand micro-dissected using Paxinos mouse anatomical atlas coordinates, to isolate anatomical regions of interest—PFC, AMY and HIP.25, 27, 29 Tissues were flash frozen in liquid nitrogen and stored at −80°C until future processing for RNA extraction and gene expression analyses. Approximately 1ml of blood/mouse was collected in PAXgene blood RNAcollection tubes (PreAnalytix, Qiagen, San Jose, CA, USA). The PAXgene tubes were stored at 4°C overnight, and then at −80°C until future processing for RNA extraction.
We used Mouse Genome 430 2.0 arrays (Affymetrix, Santa Clara, CA, USA). The GeneChip Mouse Genome 430 2.0 Array contain over 45000 probe sets that analyze the expression level of over 39000 transcripts and variants from over 34000 well-characterized mouse genes. Microarrays used in each independent experiment were derived from the same manufacturing lot.
For each brain region (PFC, AMY and HIP) and blood, equal amounts of total RNA extracted from tissue samples were pooled within each biological experiment (three mice per treatment group), and then used for labeling and microarray assays.
Standard techniques were used to obtain total RNA (22 gauge syringe homogenization in RLT buffer) and to purify the RNA (RNeasy mini kit, Qiagen) from micro-dissected mouse brain regions. For the whole mouse blood RNA extraction, PAXgene blood RNA extraction kit (PreAnalytiX, a Qiagen/BD Biosciences, San Jose, CA, USA) was used, followed by GLOBINclearTM–Mouse/Rat (Ambion/Applied Biosystems, Austin, TX, USA) to remove the globin mRNA. All the methods and procedures were carried out as per the manufacturer's instructions. The quality of the total RNA was confirmed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). The quantity and quality of total RNA was also independently assessed by 260nm ultraviolet absorption and by 260/280 ratios, respectively. Starting material of total RNA labeling reactions was kept consistent within each independent microarray experiment.
Standard Affymetrix protocols were used to reverse transcribe the messenger RNA and generate biotinlylated complementary RNA (http://www.affymetrix.com/support). The amount of complementary RNA used to prepare the hybridization cocktail was kept constant within each experiment. Samples were hybridized at 45°C for 17h under constant rotation. Arrays were washed and stained using the Affymetrix Fluidics Station 400 and scanned using the Affymetrix Model 3000 Scanner controlled by GCOS software. All sample labeling, hybridization, staining and scanning procedures were carried out as per the manufacturer's recommendations.
All arrays were scaled to a target intensity of 1000 using Affymetrix MASv 5.0 array analysis software. Quality control measures including 3′/5′ ratios for glyceraldehyde 3-phosphate dehydrogenase and β-actin, scaling factors, background and Q values were within acceptable limits.
Data analysis was performed using Affymetrix Microarray Suite 5.0 software (MAS v5.0). Default settings were used to define transcripts as present (P), marginal (M) or absent (A). A comparison analysis was performed for each drug treatment, using its corresponding saline vehicle treatment as the baseline. ‘Signal', ‘detection', ‘signal log ratio', ‘change' and ‘change P-value,' were obtained from this analysis. An empirical P-value threshold for change of P<0.00025 was used. Only transcripts that were called present in at least one of the two samples (saline vehicle or drug) intra-experiment, and that were reproducibly changed in the same direction in at least two out of three independent experiments, were analyzed further.
The identities of transcripts were established using NetAFFX (Affymetrix). Probe-sets that did not have a known gene were labeled ‘EST' and their accession numbers kept as identifiers.
We have established in our laboratory (Laboratory of Neurophenomics, Indiana University School of Medicine, www.neurophenomics.info) manually curated databases of all the human gene expression (postmortem brain, blood), human genetic (association, linkage) and animal model gene expression studies published to date on psychiatric disorders.21 Only the findings deemed significant in the primary publication, by the study investigators, using their particular experimental design and thresholds, are included in our databases. These constantly updated large databases have been used in our CFG cross-validation (Figure 1).
To designate convergence for a particular gene, the gene had to have published evidence of association or linkage for anxiety disorders, including PTSD, OCD, panic disorder and phobias. For linkage, the location of each gene was obtained through GeneCards (http://www.genecards.org), and the sex averaged cM location of the start of the gene was then obtained through http://compgen.rutgers.edu/old/map-interpolator/. For convergence, per our previously published criteria,25 the start of the gene had to map within 10cM of the location of a marker linked to the disorder.
Information about genes was obtained and imported in our databases searching the primary literature with PubMed (http://ncbi.nlm.nih.gov/PubMed), using various combinations of keywords (gene name, anxiety, stress, phobia, panic, PTSD, OCD, human, brain, postmortem, blood, lymphocytes, fibroblasts). Convergence was deemed to occur for a gene if there were published human postmortem brain data (or, rarely, blood and other tissue data) showing changes in expression of that gene in tissue from patients with anxiety and related disorders.
To search for mouse genetic evidence—QTL or transgenic—for our candidate genes, we utilized the MGI_3.54—Mouse Genome Informatics (http://www.informatics.jax.org). (Jackson Laboratory) and used the search ‘Genes and Markers' form to find QTL or transgenic for Mammalian Phenotype Ontology category ‘abnormal emotion/affect behavior', which includes the following sub-categories: abnormal fear/anxiety-related behavior, abnormal response to novelty and aggression-related behavior. To designate convergence for a particular gene, the gene had to map within 10 cM of a QTL marker for the abnormal behavior, or a transgenic mouse of the gene itself displayed that behavior.
For animal model brain and blood gene expression evidence, we have used in addition to our own data, published reports from the literature, curated in our databases.
Only genes reproducibly changed in expression in the same mouse tissue (PFC, AMY, HIP and blood), in the same direction, in two out of three independent experiments, were analyzed further. The three internal lines of evidence (pharmacological treatments—changed in yohimbine, changed in diazepam, no change in co-treatment) were scored with 1 point each. The six external cross-validating lines of evidence (three animal models, three human) were: animal model genetic data, animal model brain gene expression data, animal model blood gene expression data, human genetic data, human brain gene expression data and human blood gene expression data (Figure 1d). The lines of evidence received a maximum of 1 point each (for animal model genetic data, 0.5 points if it was QTL, 1 point if it was transgenic; for human genetic data, 0.5 points if it was linkage, 1 point if it was association). Thus the maximum possible CFG score for each gene was 3+6=9.
The more lines of evidence, that is, the more times a gene shows up as a positive finding across independent studies, platforms, methodologies and species, the higher its CFG score (Figure 1d). This is very similar conceptually to a Google PageRank algorithm, in which the more links to a page, the higher it comes up on the search prioritization list.23 Human and animal model, genetic and gene expression, data sets were integrated and tabulated. It has not escaped our attention that other ways of weighing the scores of line of evidence may give slightly different results in terms of prioritization, if not in terms of the list of genes per se. Nevertheless, this simple scoring system, where the different independent lines of evidence are weighted equally, and more of the lines of evidence are related to gene expression rather than genetics, arguably provides a good separation and prioritization of genes and blood biomarkers that are changed in expression and disease relevant, our stated focus.
Ingenuity 8.5 (Ingenuity Systems, Redwood City, CA, USA) was used to analyze the biological roles, including top canonical pathways, of the candidate genes resulting from our work (Table 5, Supplementary Table S2), as well as employed to identify genes in our data sets that are the target of existing drugs (Supplementary Table S4). GeneGo (Thompson Reuters) was used to analyze the disease categories of the genes identified (Table 7, Supplementary Table S3).
Our pharmacogenomics animal model displays a behavioral readout consistent with the drugs having an impact and their intended effects—anxiogenic for yohimbine, anxiolytic for diazepam and mitigation of effects for co-treatment (Figure 2).
We have a relatively large number of genes changed in expression in the mouse tissues examined (three brain regions and blood) (Table 1).
To start with, we have grouped the mouse model gene expression changes into categories I–IV, as described in Figure 1c and Table 1. We reasoned that genes that are category I genes, which are changed in expression by both the agonist and antagonist, as well as not changed (‘nipped in the bud') by co-treatment, are more likely to be involved in the core biology of anxiety disorders rather than be pleiotropic effects/side-effects of the drugs we used. Of note, the HIP and the blood have a relatively greater proportion of category I genes than the other brain regions (Table 1), suggesting an important role in anxiety disorders for the HIP, and a possible peripheral effect/biomarker readout for the blood.
For CFG scoring, each internal pharmacological line of evidence (changed in expression by yohimbine, changed by diazepam and not changed by co-treatment) was scored separately, along with each of the six external lines of evidence (three from animal model studies, and three from human studies), resulting in a maximum possible CFG score of 9 (Figure 1). Genes that have a CFG score of 4 or above, i.e. they have at least one full external line of evidence in addition to the maximal possible score of 3 from the internal evidence, were prioritized and shown in Table 2 and Figure 3. The average CFG score for the top candidate genes (Table 2) was again highest for HIP (4.4), followed by AMY (4.37), PFC (4.34) and blood (4.26). The relative role of HIP in anxiety disorders may thus be more important than previously appreciated, consistent with recent work in the field.32, 39
Our analysis identified and prioritized a number of top candidate genes (Figure 3 and Table 2), some well-known for involvement in anxiety, some less well known, such as FOS, GABBR1, NR4A2, DRD1, ADORA2A, QKI, RGS2, PTGDS, DYNLL2 and CCKBR. FOS (FBJ murine osteosarcoma viral oncogene homolog) is an oncogene as well as an immediate early response gene. It is a transcription factor involved in cellular reactivity to external signals. In our studies, it is also a top brain–blood biomarker for anxiety, concordantly changed in the AMY, HIP and blood. Interestingly, there is previous evidence of increase in expression of FOS in blood from PTSD patients.40 GABBR1 (gamma-aminobutyric acid (GABA) B receptor, 1) has a critical role in the fine-tuning of inhibitory synaptic transmission mediated by GABA. Our work provides evidence for its involvement in the AMY in anxiety (Table 2). Like FOS, it is also changed (decreased) in expression in blood from PTSD patients.40 GABBR1 has previous evidence suggestive for genetic association with OCD41 and with schizophrenia.42 NR4A2 (nuclear receptor subfamily 4, group A, member 2) is a steroid receptor family member, as well as immediate early response gene. It is a transcription factor involved in cellular reactivity to external signals, with a role in dopaminergic neuron development. Our work provides evidence for its involvement in the HIP in anxiety (Table 2). NR4A2 has previous evidence suggestive for genetic mutations43 and brain expression changes44 in schizophrenia and bipolar disorder. DRD1 (dopamine receptor 1), for which our work provides evidence for its involvement in the PFC and HIP in anxiety (Table 2), has previous evidence suggestive for genetic association in panic disorder.45 ADORA2A (adenosine A2a receptor), is a receptor for adenosine. The activity of this receptor is mediated by G proteins which activate adenylyl cyclase. Our work provides evidence for the involvement of ADORA2A in the AMY in anxiety (Table 2). There is previous evidence suggestive for genetic association in panic disorder.46, 47 Notably, with the exception of FOS, all the above discussed top candidate genes for anxiety have also been previously identified by our CFG work as being among top candidate genes for schizophrenia29 (Table 6, Supplementary Figure S1). QKI (quaking homolog, KH domain RNA binding), a RNA-binding protein, has a central role in myelination. In our studies, it is also a top brain–blood biomarker for anxiety, concordantly changed in the AMY and blood. Interestingly, there is previous evidence of increase in expression of QKI in blood from humans subjected to chronic stress.48 Finally, among our top candidate genes are RGS2, DYNLL2, PTGDS and CCKBR, all of which have previous human genetic association evidence for anxiety disorders and thus serve as a de facto positive control for our pharmacogenomic approach. Of note, PTGDS and CCKBR are also top candidate genes for schizophrenia in our previous work29 (Table 6, Supplementary Figure S1).
In addition, we have looked at what genes were changed in expression in all three brain regions studied (Table 3), on the premise they are more likely to be involved in the core biology of anxiety. Notably, EGR2 (early growth response 2) and SGK1 (serum/glucocorticoid-regulated kinase 1), which are involved in cellular reactivity to external signals and stress, have high CFG scores (i.e. multiple converging lines of evidence) for involvement in anxiety disorders.
Genes that are changed in expression in one of the key brain regions studied and in blood are candidate blood biomarkers.22 We used a narrow interpretation of what can constitute a candidate blood biomarker (Table 4), i.e. the change in gene expression in brain and blood has to be co-directional, inside the same drug treatment arm. FOS, QKI, HSPA1B and DNAJB1 are the top candidate biomarkers under this definition. There is more overlap between brain and blood if co-directionality of expression is not a criterion (Supplementary Table S1), as different tissues (and brain regions) can show different directions of gene expression changes. Moreover, there may be an even more significant overlap between brain and blood at a biological pathway level (Table 5 and Figure 4), where the same top pathways, if not necessarily the same genes, show alterations. Notably the glucocorticoid receptor signaling pathway and the CCR5 signaling pathway are altered in anxiety in both AMY and blood. In the end, panels of biomarkers and pathways need to be clinically validated, i.e. show predictive ability for anxiety state or response to treatment in independent human studies.
First, we carried out biological pathway analyses on all the genes that were changed in expression in our pharmacogenomic model, without any CFG prioritization (Supplementary Table S2). This may give a view of pathways involved in anxiety in the brain, but probably includes other pleiotropic effects of the drugs used.
Next, we carried out pathway analyses on the top candidate genes prioritized by CFG (CFG score of 4.0 and above) (Table 5). The resulting pathways are likely more specific to the core illness phenomenology, and less pleiotropic. Among these top biological pathways altered in anxiety, cAMP is changed in common in all three brain regions studied (Figure 4). cAMP signaling is fundamental to cellular reactivity to external signals. Previous evidence has been suggestive of a role for cAMP signaling pathways in anxiety disorders,49, 50, 51 but our work is the first to identify it as a core mechanism for anxiety across different brain regions.
We also identified biological pathways involved in anxiety specific to the different brain regions we studied. In the PFC, after cAMP signaling, the top pathway is Huntington's disease signaling. This pathway is also a top pathway altered in the blood in our analyses, suggesting its potential as a biomarker repository (Figure 4). In the AMY, the top pathway is aldosterone signaling. Previous work in animal models has suggested a role for the mineralocorticoid pathway in anxiety and stress response.52 Glucocorticoid receptor signaling and CCR5 signaling are other top pathways in the AMY, as well as in blood (Figure 4). In the HIP, after cAMP signaling, the top pathway is corticotropin-releasing hormone signaling. This pathway is well established in anxiety and stress response,53, 54 and serves as a reassuring positive control for our own work and analyses.
We have used a comprehensive, CFG approach for identifying high probability candidate genes, pathways and mechanisms for anxiety and related disorder, by the integration in a Bayesian fashion of multiple independent converging lines of evidence. This mapping of the genomic landscape of anxiety disorders completes our triad of first-pass mapping efforts of major psychiatric disorders domains—bipolar disorder,25, 26, 55 schizophrenia,29 and now anxiety disorders.
Our convergent approach emphasizes gene expression evidence more than genetic evidence, i.e. more of the scored lines of evidence come from gene expression studies than from genetic studies (Figure 1). Gene expression studies are arguably a better way to understand biology than genetic studies. After all, gene expression is the result of integration of the effects of many genetic polymorphisms, epigenetic changes and environmental effects, whereas genetics looks too early in this chain of events, and in a narrow fashion. Biologically important genes can thus be identified and studied in action at a gene expression level, whereas at a genetic level the complexity and heterogeneity of genetic polymorphisms precludes easy identification and gives no indication of their actual biological activity. The advantage of gene expression studies over genetic studies, including sequencing, may be magnified by evolutionary considerations of increased genetic heterogeneity in highly biologically active and environmentally reactive genes, such as brain and immune system genes, as a way of permitting adaptation to the environment.23 Moreover, as per our earlier formulation that ‘genes that change together (may) act together',24 the co-expression data sets we have generated in various brain regions offer testable hypotheses for transcriptional co-regulation, and for epistatic interactions among the corresponding loci.56
An acute treatment model like the one we are using is not necessarily inductive to assessing the long-term changes associated with anxiety, such as functional and structural changes apparent on imaging. Although we have no direct way of knowing if some of the genes we captured with our screen are involved or not in setting in motion such long-term changes, it is to be noted that some of these gene changes have also been reported in genetic studies of anxiety and anxiety-related disorders. Moreover, we have candidate genes in our data set with roles in brain infrastructure, including myelination (Table 2). More chronic treatments should, nevertheless, be pursued to verify and expand the findings presented in this paper.
Different combinations of anxiogenic and anxiolytic agents could be used in a comprehensive functional pharmacogenomic approach such as the one we have described. They could conceivably lead to different results, which would be of interest and welcome, since it is unlikely we are capturing with our model the full spectrum of gene expression changes involved in anxiety. However, if those drug combinations indeed mimic and modulate the same core phenomenology, the Venn diagrams of the overlap between different drug treatments will be of high interest in terms of identifying the key molecular players involved in the effects, as opposed to those involved in the (very different) side-effects of the individual drugs.
It is to be noted that our experimental approach for detecting gene expression changes relies on a single methodology, Affymetrix GeneChip oligonucleotide microarrays. It is possible that some of the gene expression changes detected from a single biological experiment, with a one-time assay with this technology, are biological or technical artifacts. With that in mind, we have designed our experiments to minimize the likelihood of having false positives, even at the expense of having false negatives. Working with an isogenic mouse strain affords us an ideal control baseline of saline vehicle injected animals for our drug-injected animals. We performed three independent de novo biological experiments, at different times, with different batches of mice (Figure 1b). We have pooled material from three mice in each experiment, and carried out microarray studies. The pooling process introduces a built in averaging of signal. We used a Venn diagram approach and only considered the genes that were reproducibly changed in the same direction in at least two out of three independent experiments. This overall design is geared to factor out both biological and technical variability. It is to be noted that the concordance between reproducible microarray experiments using the latest generations of oligonucleotide microarrays and other methodologies such as quantitative PCR, with their own attendant technical limitations, is estimated to be over 90%.57 Moreover, our CFG approach, as described above, is predicated on the existence of multiple internal and external cross-validators for each gene that is reproducibly changed in expression (Figure 1). These cross-validators are derived from independent gene expression or genetic experiments.
The results presented in this paper have a series of direct implications. First, in terms of pharmacotherapy and drug development, some of the candidate genes in our data set encode for proteins that are modulated by existing pharmacological agents (Supplementary Table S4), which may suggest future avenues for rational polypharmacy using currently available agents. Notably, existing drugs approved for other indications, such as dopaminergic agents, ion channel blockers, baclofen, nitrates, lipid modulators and disulfiram (Antabuse) are potential augmentation options for existing first-line anxiolytics and merit careful exploration as such. Some of the top anxiety candidate genes (FOS, PTGDS, HOMER1, NR4A2, GSK3B and LPL) are also modulated by the omega-3 fatty acid DHA in recent animal model studies carried by us (Le-Niculescu et al., Transl Psychiatry (2011) 1, e4; doi:10.1038/tp.2011.1), providing a potential non-pharmacological alternative for treatment. Our data sets of the effects of yohimbine and diazepam on gene expression in different key brain regions (Table 2) may be used as a source of new targets for drug development. The candidate biomarkers identified by us may, upon future validation, aid with drug development, monitoring response to treatment and early clinical intervention. Heterogeneity is possible, indeed likely, in individual human subjects—a fertile direction for future studies.26, 28, 30 Targeting key pathways identified by us (Figure 4) may provide broader options than targeting individual genes, for both drug development and peripheral blood readouts.
Second, despite using lines of evidence for our CFG approach that have to do only with anxiety disorders, the list of genes identified has a notable overlap with other psychiatric disorders, and with medical disorders (Tables 5 and 7, Supplementary Tables S2 and S3). This is a topic of major interest and debate in the field.58, 59 We demonstrate an overlap between top candidate genes for anxiety and candidate genes for schizophrenia and bipolar disorder, as well as alcoholism previously identified by us through CFG (Table 6 and Supplementary Figure S1), thus providing a possible molecular basis for the frequently observed clinical co-morbidity and interdependence between anxiety and those other major psychiatric disorders. Notably, PDE10A and TAC1 are at the overlap of all three major psychiatric domains, and may be of major interest for drug development.60, 61, 62 Among our top candidate genes for anxiety are DBP and RORB, circadian clock genes previously identified by us as candidate genes for bipolar disorder27, 55, 63 (Table 6 and Supplementary Figure S1). In addition to mood symptoms, we had previously demonstrated that DBP knock-out mice exhibit increased reactivity to stress, as well as increased alcohol consumption.27 NPAS4, another circadian gene in our anxiety dataset, is changed in expression in all three brain regions studied (Table 3). NPAS4 is a transcription factor that acts as a heterodimer partner for ARNTL, another top candidate gene for bipolar disorder identified by our previous work.25, 26, 64 The involvement of circadian genes in anxiety may underlie anxiety effects on sleep, diurnal variations in anxiety (for example, higher at night), and cycling in levels of anxiety symptoms in some patients—similar too, driven by or driving mood symptoms (cycloanxiety vs cyclothymia).6 Another top candidate gene at the overlap of bipolar disorder and anxiety is PENK (preproenkephalin). Our work provides evidence for the involvement of PENK in the PFC and HIP in anxiety (Table 2). Endogenous opiates may signal that the environment is favorable, improving mood and decreasing anxiety. As such, exogenous opiate drugs may be effective for treatment, but highly addictive. Unexpectedly, there is a major overlap between schizophrenia and anxiety, both at a top candidate genes level (Supplementary Figure S1 and Table 6) and at a pathway analyses level (Table 7 and Supplementary Table S3). Clinically, while there are some reports of co-morbidity between schizophrenia and anxiety,65, 66, 67 it is an area that has possibly been under-appreciated and understudied. Based on our work and the body of evidence in the field, we propose that a new diagnostic category of ‘schizoanxiety disorder' may have heuristic value and pragmatic clinical utility, similar to schizoaffective disorder.
Third, the mechanistic understanding and model for anxiety that emerges out of the candidate gene identified and the analyses of biological pathways involved points to signal transduction and reactivity to signals from the external environment and internal milieu (Figure 5). Notably, pathways involved in cellular stress and heat shock response (involving HSPA1B, HSPA8, HSPA4, HSPA13) seem to have been recruited by evolution for higher whole-body and mental functions68 such as anxiety. The cybernetic-like simplicity of the model should not overshadow the important fact that it is the result of the empirical coalescence of data in a non-hypothesis driven, discovery type approach. The implications for understanding the pathophysiology and treatment of anxiety and related disorders are profound. One needs to correct cellular, brain and whole organism reactivity to the environment in the treatment of these disorders. It is a place where psychopharmacology, management of medical problems, cognitive–behavioral therapy and social integration can and should go hand in hand.
In conclusion, we propose that our comprehensive CFG approach is a useful starting point in helping unravel the complex genetic code and neurobiology of anxiety and related disorders, and generates a series of leads for both future research and clinical practice.
We would like to acknowledge our debt of gratitude for the efforts and results of the many other groups, cited in our paper, who have conducted empirical studies (animal model and human, genetic and gene expression) in anxiety disorders. Without their arduous and careful work, a convergent approach such as ours would not be possible. This work was supported by an NIH Directors' New Innovator Award (1DP2OD007363) and a VA Merit Award (1I01CX000139-01) to ABN. Microarray studies were carried out in the Center for Medical Genomics at Indiana University School of Medicine.
The authors declare no conflict of interest.
Supplementary Information accompanies the paper on the Translational Psychiatry website (http://www.nature.com/tp)