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.

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.

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.

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.

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.

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.

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.

FP-TDI SNP scoring by manual and statistical procedures: a study of error rates and types.

For technologies that are commonly used in ordinary laboratories such as fluorescence-polarization detection with template-directed, dye-terminator incorporation (FP-TDI), SNP genotype scoring is usually done manually. Here we study rates of errors and missing genotypes obtained with this procedure. We also introduce three statistical genotype scoring methods to examine whether they form a viable alternative. Data consisted of eight SNPs typed in about 1400 individuals from 268 pedigrees. The statistical procedures performed better on several internal criteria, such as the number of Mendelian errors, and showed much higher agreement with discrepant genotypes re-scored by two raters. The best results were obtained with the statistical procedure that incorporated information about regularities in the error structure of the FP-TDI data. We estimated that there were about 1.6% more errors if genotypes were scored manually. About 0.6% of these errors could be explained by data manipulation errors, leaving 1% as the result of possible incorrect scoring. There were 3.3% more missing genotypes in the manual scoring due to errors in data manipulation (1.7%) and conservative scoring (1.6%).

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.

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.

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.

Will haplotype maps be useful for finding genes?

From its introduction into the literature, the idea of haplotype map-based linkage disequilibrium (LD) studies has been the subject of disputes. These queries involve the extent to which the haplotype blocks exist, the validity of fundamental concepts such as the recombination hotspot, and the application of this idea in the form of the HapMap project. In this article, we review the relevant literature to evaluate the potential importance of haplotype maps for psychiatric genetics. We first take a closer look at the nature of haplotype blocks and then address the impact of block definitions and methodological factors, such as single-nucleotide polymorphism density and sample size, on findings from haplotype block studies. After distinguishing between two types of haplotype map-based LD studies, we discuss the importance of the recombination hotspot and the nature of the disease mutations affecting complex traits. In the final section, we summarize our main conclusions and comment on the usefulness of haplotype maps for finding genes.

Using shared and unique parental views to study the etiology of 7-year-old twins' internalizing and externalizing problems.

In a sample of 1,940 Dutch 7-year-old twin pairs we studied the etiology of individual differences in Internalizing and Externalizing behavioral problems. For the majority of twins in the sample, both maternal and paternal ratings of behavioral problems were obtained from the Child Behavior Checklist. This made it possible to take into account processes underlying agreement and disagreement between maternal and paternal ratings. For both problem behaviors, a Psychometric model fitted the data better than a Rater Bias model, implying that parents, in addition to the behaviors they similarly observed, also assessed unique aspects of their children's behaviors. Relatively large genetic influences were found for Externalizing problems, explaining over 50% of the variance in both boys and girls. For internalizing problems, the heritability was over 30% in both sexes. Shared environmental factors were nearly as important as genetic influences in explaining the variation in behavioral problems. For both Externalizing and Internalizing problems, around 30% of the variance was accounted for by the shared environmental factors.

Identification of a high-risk haplotype for the dystrobrevin binding protein 1 (DTNBP1) gene in the Irish study of high-density schizophrenia families.

A recent report showed significant associations between several SNPs in a previously unknown EST cluster with schizophrenia. (1). The cluster was identified as the human dystrobrevin binding protein 1 gene (DTNBP1) by sequence database comparisons and homology with mouse DTNBP1. (2). However, the linkage disequilibrium (LD) among the SNPs in DTNBP1 as well as the pattern of significant SNP-schizophrenia association was complex. This raised several questions such as the number of susceptibility alleles that may be involved and the size of the region where the actual disease mutation(s) could be located. To address these questions, we performed different single-marker tests on the 12 previously studied and 2 new SNPs in DTNBP1 that were re-scored using an improved procedure, and performed a variety of haplotype analyses. The sample consisted of 268 Irish multiplex families selected for high density of schizophrenia. Results suggested a simple structure where the LD in the target region could be explained by 6 haplotypes that together accounted for 96% of haplotype diversity in the whole sample. From these six, a single high-risk haplotype was identified that showed a significant association with schizophrenia and explained the pattern of significant findings in the analyses with individual markers. This haplotype was 30 kb long, had a large effect, could be measured with two tag SNPs only, had a frequency of 6% in our sample, seemed to be of relatively recent origin in evolutionary terms, and was equally distributed over Ireland. Implications of these findings for follow-up and replication studies are discussed.

Co-occurrence of aggressive behavior and rule-breaking behavior at age 12: multi-rater analyses.

Aggressive Behavior (AGG) and Rule-Breaking Behavior (RB) are two of the eight CBCL syndromes. The phenotypic correlation between AGG and RB ranges from.48 to.76, and varies depending on the rater and the sex of the child. Prevalence of AGG and RB (i.e., T > or = 67) is in the range of 6%-7% in both boys and girls. Fifty percent to 60% of the children who are deviant on AGG are also deviant on RB and vice versa. Why so many children show problem behavior in the clinical range for both syndromes is unclear. This co-occurrence could be due to genetic factors influencing both traits, to environmental factors influencing both traits, or to both. The purpose of this study is to use a genetically informative sample to estimate genetic and environmental influences on AGG and RB and to investigate the etiology of the co-occurrence of both behaviors. We do this using multiple informants to take into account underlying sources of parental agreement and disagreement in ratings of their offspring. To this end, mother and father ratings of AGG and RB were collected by using the Child Behavior Checklist in a large sample of 12-year-old twins. Parental agreement is represented by an interparent correlation in the range of .53-.76, depending on phenotype (AGG or RB) and sex of the child. Genetic influences account for 79% and 69% of the individual differences in RB and AGG behavior (defined as AGG and RB on which both parents do agree) in boys. In girls 56% and 72% of the variance in RB and AGG are accounted for by genetic factors. Shared environmental influences are significant for RB in girls only, explaining 23% of the total variance. Eighty percent of the covariance between AGG and RB, similarly assessed by both parents, can be explained by genetic influences. So, co-occurrence in AGG and RB is mainly caused by a common set of genes. Parental disagreement seems to be a combination of so-called rater bias and of parental specific views.