Meta-analysis of genome-wide association studies of anxiety disorders.

Anxiety disorders (ADs), namely generalized AD, panic disorder and phobias, are common, etiologically complex conditions with a partially genetic basis. Despite differing on diagnostic definitions based on clinical presentation, ADs likely represent various expressions of an underlying common diathesis of abnormal regulation of basic threat-response systems. We conducted genome-wide association analyses in nine samples of European ancestry from seven large, independent studies. To identify genetic variants contributing to genetic susceptibility shared across interview-generated DSM-based ADs, we applied two phenotypic approaches: (1) comparisons between categorical AD cases and supernormal controls, and (2) quantitative phenotypic factor scores (FS) derived from a multivariate analysis combining information across the clinical phenotypes. We used logistic and linear regression, respectively, to analyze the association between these phenotypes and genome-wide single nucleotide polymorphisms. Meta-analysis for each phenotype combined results across the nine samples for over 18 000 unrelated individuals. Each meta-analysis identified a different genome-wide significant region, with the following markers showing the strongest association: for case-control contrasts, rs1709393 located in an uncharacterized non-coding RNA locus on chromosomal band 3q12.3 (P=1.65 × 10(-8)); for FS, rs1067327 within CAMKMT encoding the calmodulin-lysine N-methyltransferase on chromosomal band 2p21 (P=2.86 × 10(-9)). Independent replication and further exploration of these findings are needed to more fully understand the role of these variants in risk and expression of ADs.

Behavioral metabolomics analysis identifies novel neurochemical signatures in methamphetamine sensitization.

Behavioral sensitization has been widely studied in animal models and is theorized to reflect neural modifications associated with human psychostimulant addiction. While the mesolimbic dopaminergic pathway is known to play a role, the neurochemical mechanisms underlying behavioral sensitization remain incompletely understood. In this study, we conducted the first metabolomics analysis to globally characterize neurochemical differences associated with behavioral sensitization. Methamphetamine (MA)-induced sensitization measures were generated by statistically modeling longitudinal activity data for eight inbred strains of mice. Subsequent to behavioral testing, nontargeted liquid and gas chromatography-mass spectrometry profiling was performed on 48 brain samples, yielding 301 metabolite levels per sample after quality control. Association testing between metabolite levels and three primary dimensions of behavioral sensitization (total distance, stereotypy and margin time) showed four robust, significant associations at a stringent metabolome-wide significance threshold (false discovery rate, FDR <0.05). Results implicated homocarnosine, a dipeptide of GABA and histidine, in total distance sensitization, GABA metabolite 4-guanidinobutanoate and pantothenate in stereotypy sensitization, and myo-inositol in margin time sensitization. Secondary analyses indicated that these associations were independent of concurrent MA levels and, with the exception of the myo-inositol association, suggest a mechanism whereby strain-based genetic variation produces specific baseline neurochemical differences that substantially influence the magnitude of MA-induced sensitization. These findings demonstrate the utility of mouse metabolomics for identifying novel biomarkers, and developing more comprehensive neurochemical models, of psychostimulant sensitization.

Genome-wide pharmacogenomic study of citalopram-induced side effects in STAR*D.

Affecting about 1 in 12 Americans annually, depression is a leading cause of the global disease burden. While a range of effective antidepressants are now available, failure and relapse rates remain substantial, with intolerable side effect burden the most commonly cited reason for discontinuation. Thus, understanding individual differences in susceptibility to antidepressant therapy side effects will be essential to optimize depression treatment. Here we perform genome-wide association studies (GWAS) to identify genetic variation influencing susceptibility to citalopram-induced side effects. The analysis sample consisted of 1762 depression patients, successfully genotyped for 421K single-nucleotide polymorphisms (SNPs), from the Sequenced Treatment Alternatives to Relieve Depression (STAR(*)D) study. Outcomes included five indicators of citalopram side effects: general side effect burden, overall tolerability, sexual side effects, dizziness and vision/hearing side effects. Two SNPs met our genome-wide significance criterion (q<0.1), ensuring that, on average, only 10% of significant findings are false discoveries. In total, 12 additional SNPs demonstrated suggestive associations (q<0.5). The top finding was rs17135437, an intronic SNP within EMID2, mediating the effects of citalopram on vision/hearing side effects (P=3.27 × 10(-8), q=0.026). The second genome-wide significant finding, representing a haplotype spanning ∼30 kb and eight genotyped SNPs in a gene desert on chromosome 13, was associated with general side effect burden (P=3.22 × 10(-7), q=0.096). Suggestive findings were also found for SNPs at LAMA1, AOX2P, EGFLAM, FHIT and RTP2. Although our findings require replication and functional validation, this study demonstrates the potential of GWAS to discover genes and pathways that potentially mediate adverse effects of antidepressant medications.

Genome-wide association study of major depressive disorder: new results, meta-analysis, and lessons learned.

Major depressive disorder (MDD) is a common complex disorder with a partly genetic etiology. We conducted a genome-wide association study of the MDD2000+ sample (2431 cases, 3673 screened controls and >1 M imputed single-nucleotide polymorphisms (SNPs)). No SNPs achieved genome-wide significance either in the MDD2000+ study, or in meta-analysis with two other studies totaling 5763 cases and 6901 controls. These results imply that common variants of intermediate or large effect do not have main effects in the genetic architecture of MDD. Suggestive but notable results were (a) gene-based tests suggesting roles for adenylate cyclase 3 (ADCY3, 2p23.3) and galanin (GAL, 11q13.3); published functional evidence relates both of these to MDD and serotonergic signaling; (b) support for the bipolar disorder risk variant SNP rs1006737 in CACNA1C (P=0.020, odds ratio=1.10); and (c) lack of support for rs2251219, a SNP identified in a meta-analysis of affective disorder studies (P=0.51). We estimate that sample sizes 1.8- to 2.4-fold greater are needed for association studies of MDD compared with those for schizophrenia to detect variants that explain the same proportion of total variance in liability. Larger study cohorts characterized for genetic and environmental risk factors accumulated prospectively are likely to be needed to dissect more fully the etiology of MDD.

Antipsychotic-induced vacuous chewing movements and extrapyramidal side effects are highly heritable in mice.

Pharmacogenomics is yet to fulfill its promise of manifestly altering clinical medicine. As one example, a predictive test for tardive dyskinesia (TD) (an adverse drug reaction consequent to antipsychotic exposure) could greatly improve the clinical treatment of schizophrenia but human studies are equivocal. A complementary approach is the mouse-then-human design in which a valid mouse model is used to identify susceptibility loci, which are subsequently tested in human samples. We used inbred mouse strains from the Mouse Phenome Project to estimate the heritability of haloperidol-induced activity and orofacial phenotypes. In all, 159 mice from 27 inbred strains were chronically treated with haloperidol (3 mg kg(-1) per day via subdermal slow-release pellets) and monitored for the development of vacuous chewing movements (VCMs; the mouse analog of TD) and other movement phenotypes derived from open-field activity and the inclined screen test. The test battery was assessed at 0, 30, 60, 90 and 120 days in relation to haloperidol exposure. As expected, haloperidol caused marked changes in VCMs, activity in the open field and extrapyramidal symptoms (EPS). Unexpectedly, factor analysis demonstrated that these measures were imprecise assessments of a latent construct rather than discrete constructs. The heritability of a composite phenotype was ∼0.9 after incorporation of the longitudinal nature of the design. Murine VCMs are a face valid animal model of antipsychotic-induced TD, and heritability estimates from this study support the feasibility of mapping of susceptibility loci for VCMs.

Genome-wide association study of antipsychotic-induced QTc interval prolongation.

QT prolongation is associated with increased risk of cardiac arrhythmias. Identifying the genetic variants that mediate antipsychotic-induced prolongation may help to minimize this risk, which might prevent the removal of efficacious drugs from the market. We performed candidate gene analysis and five drug-specific genome-wide association studies (GWASs) with 492K single-nucleotide polymorphisms to search for genetic variation mediating antipsychotic-induced QT prolongation in 738 schizophrenia patients from the Clinical Antipsychotic Trial of Intervention Effectiveness study. Our candidate gene study suggests the involvement of NOS1AP and NUBPL (P-values=1.45 × 10(-05) and 2.66 × 10(-13), respectively). Furthermore, our top GWAS hit achieving genome-wide significance, defined as a Q-value <0.10 (P-value=1.54 × 10(-7), Q-value=0.07), located in SLC22A23, mediated the effects of quetiapine on prolongation. SLC22A23 belongs to a family of organic ion transporters that shuttle a variety of compounds, including drugs, environmental toxins and endogenous metabolites, across the cell membrane. This gene is expressed in the heart and is integral in mouse heart development. The genes mediating antipsychotic-induced QT prolongation partially overlap with the genes affecting normal QT interval variation. However, some genes may also be unique for drug-induced prolongation. This study demonstrates the potential of GWAS to discover genes and pathways that mediate antipsychotic-induced QT prolongation.

Pharmacogenomic study of side-effects for antidepressant treatment options in STAR*D.

BACKGROUND: Understanding individual differences in susceptibility to antidepressant therapy side-effects is essential to optimize the treatment of depression. METHOD: We performed genome-wide association studies (GWAS) to search for genetic variation affecting the susceptibility to side-effects. The analysis sample consisted of 1439 depression patients, successfully genotyped for 421K single nucleotide polymorphisms (SNPs), from the Sequenced Treatment Alternatives to Relieve Depression (STAR*D) study. Outcomes included four indicators of side-effects: general side-effect burden, sexual side-effects, dizziness and vision/hearing-related side-effects. Our criterion for genome-wide significance was a prespecified threshold ensuring that, on average, only 10% of the significant findings are false discoveries. RESULTS: Thirty-four SNPs satisfied this criterion. The top finding indicated that 10 SNPs in SACM1L mediated the effects of bupropion on sexual side-effects (p = 4.98 × 10(-7), q = 0.023). Suggestive findings were also found for SNPs in MAGI2, DTWD1, WDFY4 and CHL1. CONCLUSIONS: Although our findings require replication and functional validation, this study demonstrates the potential of GWAS to discover genes and pathways that could mediate adverse effects of antidepressant medication.

Genomewide pharmacogenomic study of metabolic side effects to antipsychotic drugs.

Understanding individual differences in the susceptibility to metabolic side effects as a response to antipsychotic therapy is essential to optimize the treatment of schizophrenia. Here, we perform genomewide association studies (GWAS) to search for genetic variation affecting the susceptibility to metabolic side effects. The analysis sample consisted of 738 schizophrenia patients, successfully genotyped for 492K single nucleotide polymorphisms (SNPs), from the genomic subsample of the Clinical Antipsychotic Trial of Intervention Effectiveness study. Outcomes included 12 indicators of metabolic side effects, quantifying antipsychotic-induced change in weight, blood lipids, glucose and hemoglobin A1c, blood pressure and heart rate. Our criterion for genomewide significance was a pre-specified threshold that ensures, on average, only 10% of the significant findings are false discoveries. A total of 21 SNPs satisfied this criterion. The top finding indicated that a SNP in Meis homeobox 2 (MEIS2) mediated the effects of risperidone on hip circumference (q=0.004). The same SNP was also found to mediate risperidone's effect on waist circumference (q=0.055). Genomewide significant finding were also found for SNPs in PRKAR2B, GPR98, FHOD3, RNF144A, ASTN2, SOX5 and ATF7IP2, as well as in several intergenic markers. PRKAR2B and MEIS2 both have previous research indicating metabolic involvement, and PRKAR2B has previously been shown to mediate antipsychotic response. Although our findings require replication and functional validation, this study shows the potential of GWAS to discover genes and pathways that potentially mediate adverse effects of antipsychotic medication.

Genome-wide pharmacogenomic analysis of response to treatment with antipsychotics.

Schizophrenia is an often devastating neuropsychiatric illness. Understanding the genetic variation affecting response to antipsychotics is important to develop novel diagnostic tests to match individual schizophrenia patients to the most effective and safe medication. In this study, we use a genome-wide approach to detect genetic variation underlying individual differences in response to treatment with the antipsychotics olanzapine, quetiapine, risperidone, ziprasidone and perphenazine. Our sample consisted of 738 subjects with DSM-IV schizophrenia who took part in the Clinical Antipsychotic Trials of Intervention Effectiveness. Subjects were genotyped using the Affymetrix 500 K genotyping platform plus a custom 164 K chip to improve genome-wide coverage. Treatment outcome was measured using the Positive and Negative Syndrome Scale. Our criterion for genome-wide significance was a prespecified threshold that ensures that, on an average, only 10% of the significant findings are false discoveries. The top statistical result reached significance at our prespecified threshold and involved a single-nucleotide polymorphism (SNP) in an intergenic region on chromosome 4p15. In addition, SNPs in Ankyrin Repeat and Sterile Alpha Motif Domain-Containing Protein 1B (ANKS1B) and in the Contactin-Associated Protein-Like 5 gene (CNTNAP5), which mediated the effects of olanzapine and risperidone on Negative symptoms, were very close to our threshold for declaring significance. The most significant SNP in CNTNAP5 is nonsynonymous, giving rise to an amino-acid substitution. In addition to highlighting our top results, we provide all P-values for download as a resource for investigators with the requisite samples to carry out replication. This study demonstrates the potential of genome-wide association studies to discover novel genes that mediate the effects of antipsychotics, which could eventually help to tailor drug treatment to schizophrenic patients.

Association study of SNAP25 and schizophrenia in Irish family and case-control samples.

SNAP25 occurs on chromosome 20p12.2, which has been linked to schizophrenia in some samples, and recently linked to latent classes of psychotic illness in our sample. SNAP25 is crucial to synaptic functioning, may be involved in axonal growth and dendritic sprouting, and its expression may be decreased in schizophrenia. We genotyped 18 haplotype-tagging SNPs in SNAP25 in a sample of 270 Irish high-density families. Single marker and haplotype analyses were performed in FBAT and PDT. We adjusted for multiple testing by computing q values. Association was followed up in an independent sample of 657 cases and 411 controls. We tested for allelic effects on the clinical phenotype by using the method of sequential addition and 5 factor-derived scores of the OPCRIT. Nine of 18 SNPs had P values <0.05 in either FBAT or PDT for one or more definitions of illness. Several two-marker haplotypes were also associated. Subjects inheriting the risk alleles of the most significantly associated two-marker haplotype were likely to have higher levels of hallucinations and delusions. The most significantly associated marker, rs6039820, was observed to perturb 12 transcription-factor binding sites in in silico analyses. An attempt to replicate association findings in the case-control sample resulted in no SNPs being significantly associated. We observed robust association in both single marker and haplotype-based analyses between SNAP25 and schizophrenia in an Irish family sample. Although we failed to replicate this in an independent sample, this gene should be further tested in other samples.

Association study between the serotonin 1A receptor (HTR1A) gene and neuroticism, major depression, and anxiety disorders.

The serotonin neurotransmitter system in general, and the serotonin 1A receptor in particular, has been broadly implicated in the pathophysiology of mood and anxiety disorders, although the results of genetic association studies have been mixed. In this study, we examined the serotonin 1A receptor gene, HTR1A, for its association with shared genetic risk across a range of anxiety and depression-related phenotypes. Using multivariate structural equation modeling, we selected twin pairs from the population-based Virginia Adult Twin Study of Psychiatric and Substance Use Disorders scoring at the extremes of a latent genetic risk factor that underlies susceptibility to neuroticism, major depression, and several anxiety disorders. One member from each selected pair was entered into a 2-stage, case-control association study for the HTR1A gene. In the resulting sample of 589 cases and 539 controls, four SNPs spanning the HTR1A locus, including the C(-1019)G functional promoter polymorphism (rs6295), were screened in stage 1, the positive results of which were tested for replication in stage 2. While one marker met threshold significance criteria in stage 1, this association was not replicated in stage 2. Post-hoc analyses did not reveal association to any of the specific psychiatric phenotypes. Our data suggests that the HTR1A gene may not play a major role in the genetic susceptibility underlying depressive and anxiety-related phenotypes.

A whole genome association study of neuroticism using DNA pooling.

We describe a multistage approach to identify single nucleotide polymorphisms (SNPs) associated with neuroticism, a personality trait that shares genetic determinants with major depression and anxiety disorders. Whole genome association with 452 574 SNPs was performed on DNA pools from approximately 2000 individuals selected on extremes of neuroticism scores from a cohort of 88 142 people from southwest England. The most significant SNPs were then genotyped on independent samples to replicate findings. We were able to replicate association of one SNP within the PDE4D gene in a second sample collected by our laboratory and in a family-based test in an independent sample; however, the SNP was not significantly associated with neuroticism in two other independent samples. We also observed an enrichment of low P-values in known regions of copy number variations. Simulation indicates that our study had approximately 80% power to identify neuroticism loci in the genome with odds ratio (OR)>2, and approximately 50% power to identify small effects (OR=1.5). Since we failed to find any loci accounting for more than 1% of the variance, the heritability of neuroticism probably arises from many loci each explaining much less than 1%. Our findings argue the need for much larger samples than anticipated in genetic association studies and that the biological basis of emotional disorders is extremely complex.

Catechol-O-methyltransferase and the clinical features of psychosis.

A functional polymorphism (Val-158-Met) at the Catechol-O-methyltransferase (COMT) locus has been identified as a potential etiological factor in schizophrenia. Yet the association has not been convincingly replicated across independent samples. We hypothesized that phenotypic heterogeneity might be diluting the COMT effect. To clarify the putative association, we performed an exploratory analysis to test for association between COMT and five psychosis symptom scales. These were derived through factor analysis of the Operational Criteria Checklist for Psychiatric Illness. Our sample was the Irish Study of High Density Schizophrenia Families, a large collection consisting of 268 multiplex families. This sample has previously shown a small but significant effect of the COMT Val allele in conferring risk for schizophrenia. We tested for preferential transmission of COMT alleles from parent to affected offspring (n = 749) for each of the five factor-derived scales (negative symptoms, delusions, hallucinations, mania, and depression). Significant overtransmission of the Val allele was found for mania (P < 0.05) and depression (P = 0.01) scales. Examination of odds ratios (ORs) revealed a heterogeneous effect of COMT, whereby it had no effect on Negative Symptoms, but largest impact on Depression (OR = 1.4). These results suggest a modest affective vulnerability conferred by this allele in psychosis, but will require replication.

Association between glutamic acid decarboxylase genes and anxiety disorders, major depression, and neuroticism.

Abnormalities in the gamma-aminobutyric acid (GABA) neurotransmitter system have been noted in subjects with mood and anxiety disorders. Glutamic acid decarboxylase (GAD) enzymes synthesize GABA from glutamate, and, thus, are reasonable candidate susceptibility genes for these conditions. In this study, we examined the GAD1 and GAD2 genes for their association with genetic risk across a range of internalizing disorders. We used multivariate structural equation modeling to identify common genetic risk factors for major depression, generalized anxiety disorder, panic disorder, agoraphobia, social phobia and neuroticism (N) in a sample of 9270 adult subjects from the population-based Virginia Adult Twin Study of Psychiatric and Substance Use Disorders. One member from each twin pair for whom DNA was available was selected as a case or control based on scoring at the extremes of the genetic factor extracted from the analysis. The resulting sample of 589 cases and 539 controls was entered into a two-stage association study in which candidate loci were screened in stage 1, the positive results of which were tested for replication in stage 2. Several of the six single-nucleotide polymorphisms tested in the GAD1 region demonstrated significant association in both stages, and a combined analysis in all 1128 subjects indicated that they formed a common high-risk haplotype that was significantly over-represented in cases (P=0.003) with effect size OR=1.23. Out of 14 GAD2 markers screened in stage 1, only one met the threshold criteria for follow-up in stage 2. This marker, plus three others that formed significant haplotype combinations in stage 1, did not replicate their association with the phenotype in stage 2. Subject to confirmation in an independent sample, our study suggests that variations in the GAD1 gene may contribute to individual differences in N and impact susceptibility across a range of anxiety disorders and major depression.

Genomewide linkage study in the Irish affected sib pair study of alcohol dependence: evidence for a susceptibility region for symptoms of alcohol dependence on chromosome 4.

Alcoholism is a relatively common, chronic, disabling and often treatment-resistant disorder. Evidence from twin and adoption studies indicates a substantial genetic influence, with heritability estimates of 50-60%. We conducted a genome scan in the Irish Affected Sib Pair Study of Alcohol Dependence (IASPSAD). Most probands were ascertained through alcoholism treatment settings and were severely affected. Probands, affected siblings and parents were evaluated by structured interview. A 4 cM genome scan was conducted using 474 families of which most (96%) were comprised by affected sib pairs. Nonparametric and quantitative linkage analyses were conducted using DSM-IV alcohol dependence (AD) and number of DSM-IV AD symptoms (ADSX). Quantitative results indicate strong linkage for number of AD criteria to a broad region of chromosome 4, ranging from 4q22 to 4q32 (peak multipoint LOD=4.59, P=2.1 x 10(-6), at D4S1611). Follow-up analyses suggest that the linkage may be due to variation in the symptoms of tolerance and out of control drinking. There was evidence of weak linkage (LODs of 1.0-2.0) to several other regions, including 1q44, 13q31, and 22q11 for AD along with 2q37, 9q21, 9q34 and 18p11 for ADSX. The location of the chromosome 4 peak is consistent with results from prior linkage studies and includes the alcohol dehydrogenase gene cluster. The results of this study suggest the importance of genetic variation in chromosome 4 in the etiology and severity of alcoholism in Caucasian populations.

Genetic and environmental mechanisms underlying stability and change in problem behaviors at ages 3, 7, 10, and 12.

Maternal ratings on internalizing (INT) and externalizing (EXT) behaviors were collected in a large, population-based longitudinal sample. The numbers of participating twin pairs at ages 3, 7, 10, and 12 were 5,602, 5,115, 2,956, and 1,481, respectively. Stability in both behaviors was accounted for by genetic and shared environmental influences. The genetic contribution to stability (INT: 43%; EXT: 60%) resulted from the fact that a subset of genes expressed at an earlier age was still active at the next time point. A common set of shared environmental factors operated at all ages (INT: 47%; EXT: 34%). The modest contribution of nonshared environmental factors (INT: 10%; EXT: 6%) could not be captured by a simple model. Significant age-specific influences were found for all components, indicating that genetic and environmental factors also contributed to changes in problem behavior.