Genetic Variants on 3q21 and in the Sp8 Transcription Factor Gene (SP8) as Susceptibility Loci for Psychotic Disorders: A Genetic Association Study

Abstract

Background

Recent genome-wide association studies (GWASs) investigating bipolar disorder (BD) have detected a number of susceptibility genes. These studies have also provided novel insight into shared genetic components between BD and schizophrenia (SCZ), two major psychotic disorders. To examine the replication of the risk variants for BD and the pleiotropic effect of the variants associated with BD, we conducted a genetic association study of single nucleotide polymorphisms (SNPs) that were selected based upon previous BD GWASs, which targeted psychotic disorders (BD and SCZ) in the Japanese population.

Methods

Forty-eight SNPs were selected based upon previous GWASs. A two-stage analysis was conducted using first-set screening (for all SNPs: BD = 1,012, SCZ = 1,032 and control = 993) and second-set replication samples (for significant SNPs in the screening analysis: BD = 821, SCZ = 1,808 and control = 2,149). We assessed allelic association between BD, SCZ, psychosis (BD+SCZ) and the SNPs selected for the analysis.

Results

Eight SNPs revealed nominal association signals for all comparisons (Puncorrected<0.05). Among these SNPs, the top two SNPs (associated with psychosis: Pcorrected = 0.048 and 0.037 for rs2251219 and rs2709722, respectively) were further assessed in the second-set samples, and we replicated the signals from the initial screening analysis (associated with psychosis: Pcorrected = 0.0070 and 0.033 for rs2251219 and rs2709722, respectively). The meta-analysis between the current and previous GWAS results showed that rs2251219 in Polybromo1 (PBRM1) was significant on genome-wide association level (P = 5×10−8) only for BD (P = 9.4×10−9) and psychosis (P = 2.0×10−10). Although the association of rs2709722 in Sp8 transcription factor (SP8) was suggestive in the Asian population (P = 2.1×10−7 for psychosis), this signal weakened when the samples size was increased by including data from a Caucasian population (P = 4.3×10−3).

Conclusions

We found 3p21.1 (including PBRM1, strong linkage disequilibrium made it difficult to pinpoint the risk genes) and SP8 as risk loci for BD, SCZ and psychosis. Further replication studies will be required for conclusive results.

About the Authors:

Kenji Kondo

Affiliation: Department of Psychiatry, Fujita Health University School of Medicine, Toyoake, Aichi, Japan

Masashi Ikeda

* E-mail: [email protected]

Affiliation: Department of Psychiatry, Fujita Health University School of Medicine, Toyoake, Aichi, Japan

Yusuke Kajio

Affiliation: Department of Psychiatry, Fujita Health University School of Medicine, Toyoake, Aichi, Japan

Takeo Saito

Affiliation: Department of Psychiatry, Fujita Health University School of Medicine, Toyoake, Aichi, Japan

Yoshimi Iwayama

Affiliation: Laboratory for Molecular Psychiatry, RIKEN Brain Science Institute, Wako, Saitama, Japan

Branko Aleksic

Affiliation: Department of Psychiatry, Nagoya University Graduate School of Medicine, Nagoya, Aichi, Japan

Kazuo Yamada

Affiliation: Laboratory for Molecular Psychiatry, RIKEN Brain Science Institute, Wako, Saitama, Japan

Tomoko Toyota

Affiliation: Laboratory for Molecular Psychiatry, RIKEN Brain Science Institute, Wako, Saitama, Japan

Eiji Hattori

Affiliation: Laboratory for Molecular Psychiatry, RIKEN Brain Science Institute, Wako, Saitama, Japan

Hiroshi Ujike

Affiliation: Department of Neuropsychiatry, Okayama University Graduate School of Medicine, Dentistry and Pharmaceutical Science, Okayama, Okayama, Japan

Toshiya Inada

Affiliation: Department of Psychiatry, Seiwa Hospital, Institute of Neuropsychiatry, Shinjuku, Tokyo, Japan

Hiroshi Kunugi

Affiliation: Department of Mental Disorder Research, National Institute of Neuroscience, National Center of Neurology and Psychiatry, Kodaira, Tokyo, Japan

Tadafumi Kato

Affiliation: Laboratory for Molecular Dynamics of Mental Disorders, RIKEN Brain Science Institute, Wako, Saitama, Japan

Takeo Yoshikawa

Affiliation: Laboratory for Molecular Psychiatry, RIKEN Brain Science Institute, Wako, Saitama, Japan

Norio Ozaki

Affiliation: Department of Psychiatry, Nagoya University Graduate School of Medicine, Nagoya, Aichi, Japan

Nakao Iwata

Affiliation: Department of Psychiatry, Fujita Health University School of Medicine, Toyoake, Aichi, Japan

Introduction

Bipolar disorder (BD) is a severe mental condition, and the main symptom is associated with abnormal affective status (i.e., a patient's mood will swing from manic to depression or vice-versa). The prevalence of BD worldwide is greater than 1% [1], [2], but the precise molecular mechanism is largely unknown. Nevertheless, epidemiological surveys have suggested that genetic factors contribute substantially compared with the environmental factors, and heritability has been estimated at 80% [3].

The results of genetic association studies, particularly genome-wide association studies (GWASs), have identified an increasing number of risk genes for susceptibility to BD. The initial meta-analysis of the GWASs detected a single nucleotide polymorphism (SNP) in the α subunit of the L-type voltage-gated calcium channel (CACNA1C) with suggestive statistical evidence for association (P = 7×10−8) [4]. This gene remains one of the most promising genes for BD, as the most recent mega-analysis in the Psychiatric GWAS Consortium (PGC) Bipolar Working Group provided stronger evidence of association with the SNPs in CACNA1C (P = 1.5×10−8) [5]. This mega-analysis also detected another risk SNP in Drosophila pair-rule gene ten-m (ODZ4) with genome-wide significance (P = 4.4×10−8), proposing a novel candidate gene for BD [5]. Other studies based upon genome-wide screening methodology identified several possible candidate genes for BD, such as ankyrin 3 (ANK3) [4] and neurocan (NCAN) [6].

Interestingly, genetic association studies of schizophrenia (SCZ), which is another major psychotic disorder, have also revealed that the BD risk SNPs, such as CACNA1C and ANK3, are associated with SCZ [7]. This result indicates that there is a shared genetic component between BD and SCZ. Furthermore, several independent lines of evidence converge to support this hypothesis: first, polygenic component analysis, in which the cumulative number of liberal risk alleles for SCZ or BD were enriched in the patients with BD or SCZ [8], and second, epidemiological studies, in which relatives of a proband with BD also had elevated risk for SCZ in addition to BD [9]. Therefore, combining BD and SCZ samples into a single psychosis group [10] can provide increased statistical power, specifically for the detection of overlapping risk SNPs. Several studies have used this concept also enhance detection of risk SNPs for either BD or SCZ [11], [12].

The goal of the current study was to conduct replication and cross-phenotype analyses of the SNPs that were selected based upon BD GWAS findings with two psychotic disorders (BD and SCZ) in the Japanese population. We used a two-stage design (screening and replication samples) and a meta-analysis approach to obtain reliable results.

Materials and Methods

Ethics Statement

Written informed consent was obtained from each subject after the procedures had been fully explained. This study was performed in accordance with the World Medical Association's Declaration of Helsinki and approved by the ethics committees at Fujita Health University, RIKEN BSI and institutes participating in the Collaborative Study of Mood Disorder (COSMO) [13].

Subjects

We used two independent sample sets in this study. In the first-set screening analysis, we examined 1,012 patients with BD [51.8% female, age ± standard deviation (SD) = 50.7 ± 14.3 years, BD type I = 621, BD type II = 380, schizoaffective disorder (SA) = 7, BD type information not available = 4], 1,032 SCZ (48.3% female, mean ± SD = 46.8 ± 14.8 years) and 993 healthy controls (51.1% female, age ± SD = 49.7 ± 14.0 years).

For the two SNPs that showed a significant association in the screening analysis, we used an independent second-set of samples for replication analysis. This sample consisted of 821 patients with BD (54.6% female, age ± SD = 48.2±14.4, BD type I = 387, BD type II = 344, SA = 89, BD type information not available = 1), 1,808 patients with SCZ (45.1% female, age ± SD = 49.8±14.8 years) and 2,149 healthy controls (58.3% female, age ± SD = 42.3±14.2 years). A detailed description of our subjects, including a general characterization and psychiatric assessment, is described elsewhere [13].

SNP selection and quality control

We selected 48 SNPs from BD GWAS data published prior to September 2011 [4], [5], [6], [14], [15], [16], [17], [18], [19], [20]. Regarding SNP selection, we used the following inclusion criteria. The potential risk SNPs in autosomal chromosomes must have had a P-value less than 1×10−5 if the original GWAS was conducted using a Caucasian population. The P-value must have been less than 1×10−4 if the study was based upon Asian population or PGC [5]. The minor allele frequency (MAF) must not have been equal to zero based upon the HapMap JPT panel. We used a Sequenom iPLEX Gold System (Sequenom, San Diego, CA) genotyping platform. In the optimization step, two SNPs (rs10193871 and rs1012053) were excluded due to a primer design problem. Moreover, because a visual inspection of the clustering revealed that six SNPs did not yield acceptable genotyping calls, we designed new primers for their proxy SNPs (N = 8) based upon tight linkage disequilibrium (LD). However, at this stage, we could not obtain optimal clustering for three of these SNPs. In total, we analyzed 45 SNPs (Figure S1 and Table S1). The quality control (QC) was conducted based upon the following criteria: (1) the missing call rate per person (less than 10%); (2) the missing call rate per SNP (less than 5%); and (3) a Hardy-Weinberg Equilibrium (HWE) P>0.0001 threshold (Table S1).

Statistical analysis

We assessed the allelic association of the SNPs and the following three phenotypes: 1) BD (referred to as BD association); 2) SCZ (referred to as SCZ association); and 3) psychosis (BD+SCZ; referred to as psychosis association) (Figure S2 and S3).

A comparison between multiple variables is a major concern to be addressed in a genetic study in which multiple SNPs and phenotypes are analyzed. However, thus far, no gold standard has been established. Therefore, we used a two-stage analysis and stringent cut-off level for the type I error rate in the first-set screening sample. LD between SNPs selected for analysis was calculated by SNPSpD program [21], [22] to establish an effective number of independent variables (N = 36.06). We used an adjusted statistical significance level (P<0.00138) based upon this number of independent variables. The associated SNPs from the first-set screening samples were followed-up and genotyped to replicate the association in the independent second sample set. In these analyses (first-set and second-set analyses), a one-tailed analysis was applied under a unidirectional hypothesis that risk alleles identified in the original studies were associated with risk in our dataset. We assumed this association because most of the original studies that we referred to used a larger number of subjects than those in our screening datasets [11].

A meta-analysis was conducted by combining the screening, the replication and/or the original datasets. It is worth noting that if the original dataset was involved in the PGC mega-analysis of BD, we used PGC results for the meta-analysis. A fixed model (if the I2 heterogeneity index was less than 50) or random effect model (if the I2 heterogeneity index was greater than 50) was applied in each analysis. All of the statistical procedures were calculated using PLINK version 1.07 [23].

Results

Replication analysis of the BD GWAS SNPs with BD, SCZ and psychosis in the Japanese population

After the QC calculations, 42 SNPs and 2,759 samples in the first-set screening samples were eligible for the association analysis (916 patients with BD, 946 patients with SCZ and 897 healthy controls). Table 1 lists the results, which indicated a nominal association signal (uncorrected P<0.05). Complete results are presented in Table S2. It is of note that all of the SNPs that had a nominal association with the BD sample (BD association), were also associated with SCZ (SCZ association), thus the P-values for psychosis (psychosis association) are more significant (Table 1).

[Figure omitted. See PDF.]

Table 1. Association analysis in the first-set screening samples (SNPs with nominal significant association).

https://doi.org/10.1371/journal.pone.0070964.t001

In the analysis of BD association and psychosis association, the most significant association maps to the Sp8 transcription factor (SP8) locus (rs2709736: uncorrected P = 0.0055 for BD, and rs2709722, uncorrected P = 0.0010 for psychosis), which is the same direction as the original Taiwanese population-based study [16]. However, in the analysis of SCZ association, the strongest association signal maps to chromosome 3 (52–53 Mb) and rs2251219 (uncorrected P = 0.0018) in the polybromo 1 gene (PBRM1).

Only two SNPs (rs2709722 in SP8 and rs2251219 in PBRM1) within the psychosis set remained significant after the multiple comparison correction (corrected P = 0.037 and 0.048 for rs2709722 and rs2251219, respectively). Thus, we performed a replication analysis of these two SNPs using an independent second set of samples with BD and SCZ. The associations of both SNPs were replicated, indicating a significant association with psychosis (corrected P = 0.033 and 0.0070 for rs2709722 and rs2251219, respectively). However, the SP8 SNP (rs2709722) and PBRM1 SNP (rs2251219) revealed only a nominal association level with BD or SCZ (0.01<corrected P<0.1, Table 1).

Meta-analysis of the significant SNPs detected in the first-set screening samples combining the results from our current study and the original study

In the meta-analysis, we combined our two datasets (the first-set screening samples and second-set replication samples) for the two SNPs (rs2709722 and rs2251219) to assess the association for only the Japanese population. We obtained stronger evidence of association in all of the sample sets (Table 2). Specifically, results from the psychosis sample had the most significant association (P = 8.0×10−5 for rs2251219 and P = 4.0×10−4 for rs2709722).

[Figure omitted. See PDF.]

Table 2. meta analysis of the two SNPs detected in the first-set screening analysis.

https://doi.org/10.1371/journal.pone.0070964.t002

We then combined the results from the original study [16] and/or PGC [5] datasets. For rs2251219, the original study by McMahon et al. [17] reported that the association was included in the PGC [5]; thus, we only combined the PGC results. The original study by Lee et al. [16] showed significance for rs2709722 based upon a dominant model. To conduct this meta-analysis, we recalculated the allele-wise P-value based upon the allele frequency information. For rs2251219, we detected an association with a genome-wide significance level only in the BD (P = 9.4×10−9), SCZ (P = 4.3×10−10) and psychosis sets (P = 2.0×10−10). Stronger evidence for rs2709722 was obtained in the meta-analysis (the current study and Lee's results [16]) (BD only: P = 5.1×10−7, SCZ only: P = 1.6×10−6, psychosis: P = 2.1×10−7), although this signal weakened the sample size was increased by including data from a Caucasian population (results from the current, Lee et al. and PGC studies: BD only: P = 0.0070 SCZ only: P = 0.0074 and psychosis: P = 0.0043: Table 2).

Discussion

In this study, we conducted a two-stage association analysis of the promising risk SNPs based upon BD GWASs with BD, SCZ and psychosis samples from a Japanese population. Two SNPs were detected with significant associations in the all of the phenotypes from the first-set screening (if uncorrected for multiple comparison) and second-set replication samples, indicating that these SNPs may play a role as a common risk factor for both BD and SCZ. Furthermore, we detected an association on a genome-wide significance level within the PBRM1 locus (rs2251219) by combining the results from the recent mega-analysis, which used the largest sample size thus far [5].

The SNP rs2251219, which maps to the PBRM1 locus, was originally reported in a meta-analysis of BD and major depressive disorders in a Caucasian population (P = 1.7×10−9) [17]. This SNP was also reported in the PGC BD as possessing a suggestive level of association (P = 5.5×10−7) [5]. Our meta-analysis supports this finding regarding BD because we detected a genome-wide significance (P = 9.4×10−9), even when the BD set was analyzed alone. It is also of note that the results of this meta-analysis merging our SCZ/psychosis sets and PGC BD showed genome-wide significance (SCZ and PGC BD: P = 4.3×10−10, psychosis and PGC BD: P = 2.0×10−10).

Interestingly, rs2251219 is in LD with another SNP (rs1042779, base position  = 52.8 Mb) in inter-alpha-trypsin inhibitor heavy chain 1 (ITIH1), which was examined in our screening sample (D' = 0.96 and r2 = 0.84 in an Asian population in the 1000 genome database as a reference panel) and revealed a nominal association (uncorrected P<0.05 for all phenotypes). This region (3p21.1) is one of the most attractive loci as a candidate region for risk for psychotic disorders. A recent study by Hamshere et al. [11] reported that SNP (rs2239547) in inter-alpha-trypsin inhibitor heavy chain 3 and 4 (ITIH3-4) was significantly associated with BD and/or SCZ at a genome-wide significance level (D' = 0.95 and r2 = 0.58 with rs2251219, Asian population in the 1000-genome database as a reference panel). The most recent study by the PGC Cross-Disorder Group (CDG) also reported that SNP (rs2535629) in ITIH3-4 was associated with five psychiatric disorders, and the strongest association was observed in the BD set [12]. These SNPs located in and around PBRM1, ITIH1 and ITIH3-4 are in strong LD. Thus, all of the SNPs in this LD block are promising candidates for genetic risk for BD, SCZ or psychosis. This LD structure in turn indicates that it is difficult to narrow down the true susceptibility variants within this region (Figure S4).

A Taiwanese BD GWAS found a suggestive association signal that maps to the SP8 locus (P = 4.8×10−7 in dominant model) [16]. Our result supports the association of this gene with BD, specifically for Asian populations. The meta-analysis of our results and Lee's results indicated a stronger association signal that was weakened when the sample size was increased by including data from a Caucasian population (Table 2). Considering the results of the heterogeneity test (in which the I2 score significantly increased by combining PGC BD dataset), the SP8 gene may play a role as a population-specific risk gene in individuals of Asian ancestry. SP8 is a SP transcription factor and functions in neural development by interacting with Wnt/beta-catenin and fibroblast growth factor (FGF) signaling [24]. Because there are no studies that have examined the association between SP8 and BD/SCZ, further research is needed to better understand the relationship between SP8 and psychiatric disorders. Special attention is needed regarding the population-specific effect that rs2709722 might have.

SNPs located in other promising candidate genes were not significantly associated with BD, SCZ or psychosis in the Japanese population. Although a trans-population effect is a likely explanation, our sample size may not have sufficient power to observe the association compared with the estimated effect size (odds ratio (OR) of ∼1.2). The analysis of the power of our study [25] indicated that our sample had 25% power for BD/SCZ and 40% for psychosis to detect significance (assuming an OR of 1.2) of risk with 25% MAF (average MAF in our examined SNPs in the control subjects) under an additive model (type I error rate = 0.00138). Therefore, a larger sample size in future work is essential.

In conclusion, we found two loci, PBRM1 (and neighboring genes) and SP8, that were replicated in psychotic disorders in a Japanese population. Specifically, a SNP within PBRM1 revealed genome-wide significance in the meta-analysis, suggesting promising candidate genes for BD. SP8, which was not significant on a genome-wide level, is still a candidate for a population-specific risk factor for BD.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

SNP selection strategy. JPT: HapMap Japanese Tokyo sample HWE: Hardy-Weinberg Equilibrium.

https://doi.org/10.1371/journal.pone.0070964.s001

(TIF)

Figure S2.

Sample numbers in the first-set screening samples. BD: Bipolar disorder SCZ: Schizophrenia.

https://doi.org/10.1371/journal.pone.0070964.s002

(TIF)

Figure S3.

Sample numbers in the second-set replication samples. BD: Bipolar disorder SCZ: Schizophrenia.

https://doi.org/10.1371/journal.pone.0070964.s003

(TIF)

Figure S4.

Linkage disequilibrium structure around 3q21 in the Asian population. Data are from HapMap JPT and CHB. The LD measure is based upon D'.

https://doi.org/10.1371/journal.pone.0070964.s004

(TIF)

Table S1.

SNPs selected based upon previous BD GWAS.

https://doi.org/10.1371/journal.pone.0070964.s005

(XLSX)

Table S2.

Association analysis for all SNPs.

https://doi.org/10.1371/journal.pone.0070964.s006

(XLSX)

Acknowledgments

We thank Ms. M. Miyata, Ms. M. Aizawa, Ms. Y. Umekage and Ms. M. Uchida for their technical support.

Author Contributions

Conceived and designed the experiments: KK MI NI. Performed the experiments: KK MI YK TS YI BA KY TT EH. Analyzed the data: KK MI. Contributed reagents/materials/analysis tools: TI HK TK HU TY NO. Wrote the paper: KK MI BA NI.

Citation: Kondo K, Ikeda M, Kajio Y, Saito T, Iwayama Y, Aleksic B, et al. (2013) Genetic Variants on 3q21 and in the Sp8 Transcription Factor Gene (SP8) as Susceptibility Loci for Psychotic Disorders: A Genetic Association Study. PLoS ONE 8(8): e70964. https://doi.org/10.1371/journal.pone.0070964

References

1. Merikangas KR, Jin R, He JP, Kessler RC, Lee S, et al. (2011) Prevalence and correlates of bipolar spectrum disorder in the world mental health survey initiative. Arch Gen Psychiatry 68: 241–251.

2. Craddock N, Sklar P (2009) Genetics of bipolar disorder: successful start to a long journey. Trends Genet 25: 99–105.

3. Smoller JW, Finn CT (2003) Family, twin, and adoption studies of bipolar disorder. Am J Med Genet C Semin Med Genet 123C: 48–58.

4. Ferreira MA, O'Donovan MC, Meng YA, Jones IR, Ruderfer DM, et al. (2008) Collaborative genome-wide association analysis supports a role for ANK3 and CACNA1C in bipolar disorder. Nat Genet 40: 1056–1058.

5. Psychiatric GWAS Consortium Bipolar Disorder Working Group (2011) Large-scale genome-wide association analysis of bipolar disorder identifies a new susceptibility locus near ODZ4. Nat Genet 43: 977–983.

6. Cichon S, Muhleisen TW, Degenhardt FA, Mattheisen M, Miro X, et al. (2011) Genome-wide association study identifies genetic variation in neurocan as a susceptibility factor for bipolar disorder. Am J Hum Genet 88: 372–381.

7. Williams HJ, Craddock N, Russo G, Hamshere ML, Moskvina V, et al. (2011) Most genome-wide significant susceptibility loci for schizophrenia and bipolar disorder reported to date cross-traditional diagnostic boundaries. Hum Mol Genet 20: 387–391.

8. Purcell SM, Wray NR, Stone JL, Visscher PM, O'Donovan MC, et al. (2009) Common polygenic variation contributes to risk of schizophrenia and bipolar disorder. Nature 460: 748–752.

9. Lichtenstein P, Yip BH, Bjork C, Pawitan Y, Cannon TD, et al. (2009) Common genetic determinants of schizophrenia and bipolar disorder in Swedish families: a population-based study. Lancet 373: 234–239.

10. Craddock N, Owen MJ (2010) The Kraepelinian dichotomy – going, going... but still not gone. Br J Psychiatry 196: 92–95.

11. Hamshere ML, Walters JT, Smith R, Richards AL, Green E, et al.. (2013) Genome-wide significant associations in schizophrenia to ITIH3/4, CACNA1C and SDCCAG8, and extensive replication of associations reported by the Schizophrenia PGC. Mol Psychiatry.

12. Smoller JW, Craddock N, Kendler K, Lee PH, Neale BM, et al. (2013) Identification of risk loci with shared effects on five major psychiatric disorders: a genome-wide analysis. Lancet 381: 1371–1379.

13. Matsunaga S, Ikeda M, Kishi T, Fukuo Y, Aleksic B, et al. (2012) An evaluation of polymorphisms in casein kinase 1 delta and epsilon genes in major psychiatric disorders. Neurosci Lett 529: 66–69.

14. Wellcome_Trust_Case_Control_Consortium (2007) Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 447: 661–678.

15. Baum AE, Akula N, Cabanero M, Cardona I, Corona W, et al. (2008) A genome-wide association study implicates diacylglycerol kinase eta (DGKH) and several other genes in the etiology of bipolar disorder. Mol Psychiatry 13: 197–207.

16. Lee MT, Chen CH, Lee CS, Chen CC, Chong MY, et al. (2011) Genome-wide association study of bipolar I disorder in the Han Chinese population. Mol Psychiatry 16: 548–556.

17. McMahon FJ, Akula N, Schulze TG, Muglia P, Tozzi F, et al. (2010) Meta-analysis of genome-wide association data identifies a risk locus for major mood disorders on 3p21.1. Nat Genet 42: 128–131.

18. Scott LJ, Muglia P, Kong XQ, Guan W, Flickinger M, et al. (2009) Genome-wide association and meta-analysis of bipolar disorder in individuals of European ancestry. Proc Natl Acad Sci U S A 106: 7501–7506.

19. Sklar P, Smoller JW, Fan J, Ferreira MA, Perlis RH, et al. (2008) Whole-genome association study of bipolar disorder. Mol Psychiatry 13: 558–569.

20. Smith EN, Bloss CS, Badner JA, Barrett T, Belmonte PL, et al. (2009) Genome-wide association study of bipolar disorder in European American and African American individuals. Mol Psychiatry 14: 755–763.

21. Nyholt DR (2004) A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other. Am J Hum Genet 74: 765–769.

22. Li J, Ji L (2005) Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix. Heredity (Edinb) 95: 221–227.

23. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.

24. Sahara S, Kawakami Y, Izpisua Belmonte JC, O'Leary DD (2007) Sp8 exhibits reciprocal induction with Fgf8 but has an opposing effect on anterior-posterior cortical area patterning. Neural Dev 2: 10.

25. Purcell S, Cherny SS, Sham PC (2003) Genetic Power Calculator: design of linkage and association genetic mapping studies of complex traits. Bioinformatics 19: 149–150.