SATB2 organizes the 3D genome architecture of cognition in cortical neurons.

The DNA-binding protein SATB2 is genetically linked to human intelligence. We studied its influence on the three-dimensional (3D) epigenome by mapping chromatin interactions and accessibility in control versus SATB2-deficient cortical neurons. We find that SATB2 affects the chromatin looping between enhancers and promoters of neuronal-activity-regulated genes, thus influencing their expression. It also alters A/B compartments, topologically associating domains, and frequently interacting regions. Genes linked to SATB2-dependent 3D genome changes are implicated in highly specialized neuronal functions and contribute to cognitive ability and risk for neuropsychiatric and neurodevelopmental disorders. Non-coding DNA regions with a SATB2-dependent structure are enriched for common variants associated with educational attainment, intelligence, and schizophrenia. Our data establish SATB2 as a cell-type-specific 3D genome modulator, which operates both independently and in cooperation with CCCTC-binding factor (CTCF) to set up the chromatin landscape of pyramidal neurons for cognitive processes.

Direct targets of MEF2C are enriched for genes associated with schizophrenia and cognitive function and are involved in neuron development and mitochondrial function.

Myocyte Enhancer Factor 2C (MEF2C) is a transcription factor that plays a crucial role in neurogenesis and synapse development. Genetic studies have identified MEF2C as a gene that influences cognition and risk for neuropsychiatric disorders, including autism spectrum disorder (ASD) and schizophrenia (SCZ). Here, we investigated the involvement of MEF2C in these phenotypes using human-derived neural stem cells (NSCs) and glutamatergic induced neurons (iNs), which represented early and late neurodevelopmental stages. For these cellular models, MEF2C function had previously been disrupted, either by direct or indirect mutation, and gene expression assayed using RNA-seq. We integrated these RNA-seq data with MEF2C ChIP-seq data to identify dysregulated direct target genes of MEF2C in the NSCs and iNs models. Several MEF2C direct target gene-sets were enriched for SNP-based heritability for intelligence, educational attainment and SCZ, as well as being enriched for genes containing rare de novo mutations reported in ASD and/or developmental disorders. These gene-sets are enriched in both excitatory and inhibitory neurons in the prenatal and adult brain and are involved in a wide range of biological processes including neuron generation, differentiation and development, as well as mitochondrial function and energy production. We observed a trans expression quantitative trait locus (eQTL) effect of a single SNP at MEF2C (rs6893807, which is associated with IQ) on the expression of a target gene, BNIP3L. BNIP3L is a prioritized risk gene from the largest genome-wide association study of SCZ and has a function in mitophagy in mitochondria. Overall, our analysis reveals that either direct or indirect disruption of MEF2C dysregulates sets of genes that contain multiple alleles associated with SCZ risk and cognitive function and implicates neuron development and mitochondrial function in the etiology of these phenotypes.

Genes positively regulated by Mef2c in cortical neurons are enriched for common genetic variation associated with IQ and educational attainment.

The myocyte enhancer factor 2 C (MEF2C) gene encodes a transcription factor important for neurogenesis and synapse development and contains common variants associated with intelligence (IQ) and educational attainment (EA). Here, we took gene expression data from the mouse cortex of a Mef2c mouse model with a heterozygous DNA binding-deficient mutation of Mef2c (Mef2c-het) and combined these data with MEF2C ChIP-seq data from cortical neurons and single-cell data from the mouse brain. This enabled us to create a set of genes that were differentially regulated in Mef2c-het mice, represented direct target genes of MEF2C and had elevated in expression in cortical neurons. We found this gene-set to be enriched for genes containing common genetic variation associated with IQ and EA. Genes within this gene-set that were down-regulated, i.e. have reduced expression in Mef2c-het mice versus controls, were specifically significantly enriched for both EA and IQ associated genes. These down-regulated genes were enriched for functionality in the adenylyl cyclase signalling system, which is known to positively regulate synaptic transmission and has been linked to learning and memory. Within the adenylyl cyclase signalling system, three genes regulated by MEF2C, CRHR1, RGS6, and GABRG3, are associated at genome-wide significant levels with IQ and/or EA. Our results indicate that genetic variation in MEF2C and its direct target genes within cortical neurons contribute to variance in cognition within the general population, and the molecular mechanisms involved include the adenylyl cyclase signalling system's role in synaptic function.

SATB2-LEMD2 interaction links nuclear shape plasticity to regulation of cognition-related genes.

SATB2 is a schizophrenia risk gene and is genetically associated with human intelligence. How it affects cognition at molecular level is currently unknown. Here, we show that interactions between SATB2, a chromosomal scaffolding protein, and the inner nuclear membrane protein LEMD2 orchestrate the response of pyramidal neurons to neuronal activation. Exposure to novel environment in vivo causes changes in nuclear shape of CA1 hippocampal neurons via a SATB2-dependent mechanism. The activity-driven plasticity of the nuclear envelope requires not only SATB2, but also its protein interactor LEMD2 and the ESCRT-III/VPS4 membrane-remodeling complex. Furthermore, LEMD2 depletion in cortical neurons, similar to SATB2 ablation, affects neuronal activity-dependent regulation of multiple rapid and delayed primary response genes. In human genetic data, LEMD2-regulated genes are enriched for de novo mutations reported in intellectual disability and schizophrenia and are, like SATB2-regulated genes, enriched for common variants associated with schizophrenia and cognitive function. Hence, interactions between SATB2 and the inner nuclear membrane protein LEMD2 influence gene expression programs in pyramidal neurons that are linked to cognitive ability and psychiatric disorder etiology.

Genetic and inflammatory effects on childhood trauma and cognitive functioning in patients with schizophrenia and healthy participants.

Recent studies have reported a negative association between exposure to childhood trauma, including physical neglect, and cognitive functioning in patients with schizophrenia. Childhood trauma has been found to influence immune functioning, which may contribute to the risk of schizophrenia and cognitive symptoms of the disorder. In this study, we aimed to test the hypothesis that physical neglect is associated with cognitive ability, and that this association is mediated by a combined latent measure of inflammatory response, and moderated by higher genetic risk for schizophrenia. The study included 279 Irish participants, comprising 102 patients and 177 healthy participants. Structural equation modelling was used to perform mediation and moderation analyses. Inflammatory response was measured via basal plasma levels of IL-6, TNF-α, and CRP, and cognitive performance was assessed across three domains: full-scale IQ, logical memory, and the emotion recognition task. Genetic variation for schizophrenia was estimated using a genome-wide polygenic score based on genome-wide association study summary statistics. The results showed that inflammatory response mediated the association between physical neglect and all measures of cognitive functioning, and explained considerably more variance than any of the inflammatory markers alone. Furthermore, genetic risk for schizophrenia was observed to moderate the direct pathway between physical neglect and measures of non-social cognitive functioning in both patient and healthy participants. However, genetic risk did not moderate the mediated pathway associated with inflammatory response. Therefore, we conclude that the mediating role of inflammatory response and the moderating role of higher genetic risk may independently influence the association between adverse early life experiences and cognitive function in patients and healthy participants.

Mendelian randomization analysis using GWAS and eQTL data to investigate the relationship between chronotype and neuropsychiatric disorders and their molecular basis.

Chronotype is a proxy sleep measure that has been associated with neuropsychiatric disorders. By investigating how chronotype influences risk for neuropsychiatric disorders and vice versa, we may identify modifiable risk factors for each phenotype. Here we used Mendelian randomization (MR), to explore causal effects by (1) studying the causal relationships between neuropsychiatric disorders and chronotype and (2) characterizing the genetic components of these phenotypes. Firstly, we investigated if a causal role exists between five neuropsychiatric disorders and chronotype using the largest genome-wide association studies (GWAS) available. Secondly, we integrated data from expression quantitative trait loci (eQTLs) to investigate the role of gene expression alterations on these phenotypes. Evening chronotype was causal for increased risk of schizophrenia and autism spectrum disorder and schizophrenia was causal for a tendency toward evening chronotype. We identified 12 eQTLs where gene expression changes in brain or blood were causal for one of the phenotypes, including two eQTLs for SNX19 in hippocampus and hypothalamus that were causal for schizophrenia. These findings provide important evidence for the complex, bidirectional relationship that exists between a sleep-based phenotype and neuropsychiatric disorders, and use gene expression data to identify causal roles for genes at associated loci.

Genes influenced by MEF2C contribute to neurodevelopmental disease via gene expression changes that affect multiple types of cortical excitatory neurons.

Myocyte enhancer factor 2 C (MEF2C) is an important transcription factor during neurodevelopment. Mutation or deletion of MEF2C causes intellectual disability (ID), and common variants within MEF2C are associated with cognitive function and schizophrenia risk. We investigated if genes influenced by MEF2C during neurodevelopment are enriched for genes associated with neurodevelopmental phenotypes and if this can be leveraged to identify biological mechanisms and individual brain cell types affected. We used a set of 1055 genes that were differentially expressed in the adult mouse brain following early embryonic deletion of Mef2c in excitatory cortical neurons. Using genome-wide association studies data, we found these differentially expressed genes (DEGs) to be enriched for genes associated with schizophrenia, intelligence and educational attainment but not autism spectrum disorder (ASD). For this gene set, genes that overlap with target genes of the Fragile X mental retardation protein (FMRP) are a major driver of these enrichments. Using trios data, we found these DEGs to be enriched for genes containing de novo mutations reported in ASD and ID, but not schizophrenia. Using single-cell RNA sequencing data, we identified that a number of different excitatory glutamatergic neurons in the cortex were enriched for these DEGs including deep layer pyramidal cells and cells in the retrosplenial cortex, entorhinal cortex and subiculum, and these cell types are also enriched for FMRP target genes. The involvement of MEF2C and FMRP in synapse elimination suggests that disruption of this process in these cell types during neurodevelopment contributes to cognitive function and risk of neurodevelopmental disorders.

Altered gene regulation as a candidate mechanism by which ciliopathy gene SDCCAG8 contributes to schizophrenia and cognitive function.

Mutations in genes that encode centrosomal/ciliary proteins cause severe cognitive deficits, while common single-nucleotide polymorphisms in these genes are associated with schizophrenia (SZ) and cognition in genome-wide association studies. The role of these genes in neuropsychiatric disorders is unknown. The ciliopathy gene SDCCAG8 is associated with SZ and educational attainment (EA). Genome editing of SDCCAG8 caused defects in primary ciliogenesis and cilium-dependent cell signalling. Transcriptomic analysis of SDCCAG8-deficient cells identified differentially expressed genes that are enriched in neurodevelopmental processes such as generation of neurons and synapse organization. These processes are enriched for genes associated with SZ, human intelligence (IQ) and EA. Phenotypic analysis of SDCCAG8-deficent neuronal cells revealed impaired migration and neuronal differentiation. These data implicate ciliary signalling in the aetiology of SZ and cognitive dysfunction. We found that centrosomal/ciliary genes are enriched for association with IQ, suggesting altered gene regulation as a general model for neurodevelopmental impacts of centrosomal/ciliary genes.

Genes encoding SATB2-interacting proteins in adult cerebral cortex contribute to human cognitive ability.

During CNS development, the nuclear protein SATB2 is expressed in superficial cortical layers and determines projection neuron identity. In the adult CNS, SATB2 is expressed in pyramidal neurons of all cortical layers and is a regulator of synaptic plasticity and long-term memory. Common variation in SATB2 locus confers risk of schizophrenia, whereas rare, de novo structural and single nucleotide variants cause severe intellectual disability and absent or limited speech. To characterize differences in SATB2 molecular function in developing vs adult neocortex, we isolated SATB2 protein interactomes at the two ontogenetic stages and identified multiple novel SATB2 interactors. SATB2 interactomes are highly enriched for proteins that stabilize de novo chromatin loops. The comparison between the neonatal and adult SATB2 protein complexes indicates a developmental shift in SATB2 molecular function, from transcriptional repression towards organization of chromosomal superstructure. Accordingly, gene sets regulated by SATB2 in the neocortex of neonatal and adult mice show limited overlap. Genes encoding SATB2 protein interactors were grouped for gene set analysis of human GWAS data. Common variants associated with human cognitive ability are enriched within the genes encoding adult but not neonatal SATB2 interactors. Our data support a shift in the function of SATB2 in cortex over lifetime and indicate that regulation of spatial chromatin architecture by the SATB2 interactome contributes to cognitive function in the general population.

Genes regulated by SATB2 during neurodevelopment contribute to schizophrenia and educational attainment.

SATB2 is associated with schizophrenia and is an important transcription factor regulating neocortical organization and circuitry. Rare mutations in SATB2 cause a syndrome that includes developmental delay, and mouse studies identify an important role for SATB2 in learning and memory. Interacting partners BCL11B and GATAD2A are also schizophrenia risk genes indicating that other genes interacting with or are regulated by SATB2 are making a contribution to schizophrenia and cognition. We used data from Satb2 mouse models to generate three gene-sets that contain genes either functionally related to SATB2 or targeted by SATB2 at different stages of development. Each was tested for enrichment using the largest available genome-wide association studies (GWAS) datasets for schizophrenia and educational attainment (EA) and enrichment analysis was also performed for schizophrenia and other neurodevelopmental disorders using data from rare variant sequencing studies. These SATB2 gene-sets were enriched for genes containing common variants associated with schizophrenia and EA, and were enriched for genes containing rare variants reported in studies of schizophrenia, autism and intellectual disability. In the developing cortex, genes targeted by SATB2 based on ChIP-seq data, and functionally affected when SATB2 is not expressed based on differential expression analysis using RNA-seq data, show strong enrichment for genes associated with EA. For genes expressed in the hippocampus or at the synapse, those targeted by SATB2 are more strongly enriched for genes associated EA than gene-sets not targeted by SATB2. This study demonstrates that single gene findings from GWAS can provide important insights to pathobiological processes. In this case we find evidence that genes influenced by SATB2 and involved in synaptic transmission, axon guidance and formation of the corpus callosum are contributing to schizophrenia and cognition.

Childhood trauma, parental bonding, and social cognition in patients with schizophrenia and healthy adults.

OBJECTIVE: This study investigated associations between childhood trauma, parental bonding, and social cognition (i.e., Theory of Mind and emotion recognition) in patients with schizophrenia and healthy adults. METHODS: Using cross-sectional data, we examined the recollections of childhood trauma experiences and social cognitive abilities in 74 patients with schizophrenia and 116 healthy adults. RESULTS: Patients had significantly higher scores compared with healthy participants on childhood trauma, and lower scores on parental bonding and social cognitive measures. Physical neglect was found to be the strongest predictor of emotion recognition impairments in both groups. Optimal parental bonding attenuated the impact of childhood trauma on emotion recognition. CONCLUSION: The present study provides evidence of an association between physical neglect and emotion recognition in patients with schizophrenia and healthy individuals and shows that both childhood trauma and parental bonding may influence social cognitive development. Psychosocial interventions should be developed to prevent and mitigate the long-term effects of childhood adversities.

Cognitive Genomics: Recent Advances and Current Challenges.

PURPOSE OF REVIEW: We review recent progress in uncovering the complex genetic architecture of cognition, arising primarily from genome-wide association studies (GWAS). We explore the genetic correlations between cognitive performance and neuropsychiatric disorders, the genetic and environmental factors associated with age-related cognitive decline, and speculate about the future role of genomics in the understanding of cognitive processes. RECENT FINDINGS: Improvements in genomic methods, and the increasing availability of large datasets via consortia cooperation, have led to a greater understanding of the role played by common and rare variants in the genomics of cognition, the highly polygenic basis of cognitive function and dysfunction, and the multiple biological processes involved. Recent research has aided in our understanding of the complex biological nature of genomics of cognition. Further development of data banks and techniques to analyze this data hold significant promise for understanding cognitive ability, and for treating cognitively related disability.

Genes regulated by BCL11B during T-cell development are enriched for de novo mutations found in schizophrenia patients.

While abnormal neurodevelopment contributes to schizophrenia (SCZ) risk, there is also evidence to support a role for immune dysfunction in SCZ. BCL11B, associated with SCZ in genome-wide association study (GWAS), is a transcription factor that regulates the differentiation and development of cells in the central nervous and immune systems. Here, we use functional genomics data from studies of BCL11B to investigate the contribution of neuronal and immune processes to SCZ pathophysiology. We identified the gene targets of BCL11B in brain striatal cells (n = 223 genes), double negative 4 (DN4) developing T cells (n = 114 genes) and double positive (DP) developing T cells (n = 518 genes) using an integrated analysis of RNA-seq and ChIP-seq data. No gene-set was enriched for genes containing common variants associated with SCZ but the DP gene-set was enriched for genes containing missense de novo mutations (DNMs; p = .001) using data from 3,447 SCZ trios. Post hoc analysis revealed the enrichment to be stronger for DP genes negatively regulated by BCL11B. Biological processes enriched for genes negatively regulated by BCL11B in DP gene-set included immune system development and cytokine signaling. These analyses, leveraging a GWAS-identified SCZ risk gene and data on gene expression and transcription factor binding, indicate that DNMs in immune pathways contribute to SCZ risk.

Effects of complement gene-set polygenic risk score on brain volume and cortical measures in patients with psychotic disorders and healthy controls.

Multiple genome-wide association studies of schizophrenia have reported associations between genetic variants within the MHC region and disease risk, an association that has been partially accounted for by alleles of the complement component 4 (C4) gene. Following on previous findings of association between both C4 and other complement-related variants and memory function, we tested the hypothesis that polygenic scores calculated based on identified schizophrenia risk alleles within the "complement" system would be broadly associated with memory function and associated brain structure. We tested this using a polygenic risk score (PRS) calculated for complement genes, but excluding C4 variants. Higher complement-based PRS scores were observed to be associated with lower memory scores for the sample as a whole (N = 620, F change = 8.25; p = .004). A significant association between higher PRS and lower hippocampal volume was also observed (N = 216, R2 change = 0.016, p = .015). However, after correcting for further testing of association with the more general indices of cortical thickness, surface area or total brain volume, none of which were associated with complement, the association with hippocampal volume became non-significant. A post-hoc analysis of hippocampal subfields suggested an association between complement PRS and several hippocampal subfields, findings that appeared to be particularly driven by the patient sample. In conclusion, our study yielded suggestive evidence of association between complement-based schizophrenia PRS and variation in memory function and hippocampal volume.

The effect of breed and diet type on the global transcriptome of hepatic tissue in beef cattle divergent for feed efficiency.

BACKGROUND: Feed efficiency is an important economic and environmental trait in beef production, which can be measured in terms of residual feed intake (RFI). Cattle selected for low-RFI (feed efficient) have similar production levels but decreased feed intake, while also emitting less methane. RFI is difficult and expensive to measure and is not widely adopted in beef production systems. However, development of DNA-based biomarkers for RFI may facilitate its adoption in genomic-assisted breeding programmes. Cattle have been shown to re-rank in terms of RFI across diets and age, while also RFI varies by breed. Therefore, we used RNA-Seq technology to investigate the hepatic transcriptome of RFI-divergent Charolais (CH) and Holstein-Friesian (HF) steers across three dietary phases to identify genes and biological pathways associated with RFI regardless of diet or breed. RESULTS: Residual feed intake was measured during a high-concentrate phase, a zero-grazed grass phase and a final high-concentrate phase. In total, 322 and 33 differentially expressed genes (DEGs) were identified across all diets for CH and HF steers, respectively. Three genes, GADD45G, HP and MID1IP1, were differentially expressed in CH when both the high-concentrate zero-grazed grass diet were offered. Two canonical pathways were enriched across all diets for CH steers. These canonical pathways were related to immune function. CONCLUSIONS: The absence of common differentially expressed genes across all dietary phases and breeds in this study supports previous reports of the re-ranking of animals in terms of RFI when offered differing diets over their lifetime. However, we have identified biological processes such as the immune response and lipid metabolism as potentially associated with RFI divergence emphasising the previously reported roles of these biological processes with respect to RFI.

Population-based identity-by-descent mapping combined with exome sequencing to detect rare risk variants for schizophrenia.

Genome-wide association studies (GWASs) are highly effective at identifying common risk variants for schizophrenia. Rare risk variants are also important contributors to schizophrenia etiology but, with the exception of large copy number variants, are difficult to detect with GWAS. Exome and genome sequencing, which have accelerated the study of rare variants, are expensive so alternative methods are needed to aid detection of rare variants. Here we re-analyze an Irish schizophrenia GWAS dataset (n = 3,473) by performing identity-by-descent (IBD) mapping followed by exome sequencing of individuals identified as sharing risk haplotypes to search for rare risk variants in coding regions. We identified 45 rare haplotypes (>1 cM) that were significantly more common in cases than controls. By exome sequencing 105 haplotype carriers, we investigated these haplotypes for functional coding variants that could be tested for association in independent GWAS samples. We identified one rare missense variant in PCNT but did not find statistical support for an association with schizophrenia in a replication analysis. However, IBD mapping can prioritize both individual samples and genomic regions for follow-up analysis but genome rather than exome sequencing may be more effective at detecting risk variants on rare haplotypes.

Genome-wide comparative analysis of atopic dermatitis and psoriasis gives insight into opposing genetic mechanisms.

Atopic dermatitis and psoriasis are the two most common immune-mediated inflammatory disorders affecting the skin. Genome-wide studies demonstrate a high degree of genetic overlap, but these diseases have mutually exclusive clinical phenotypes and opposing immune mechanisms. Despite their prevalence, atopic dermatitis and psoriasis very rarely co-occur within one individual. By utilizing genome-wide association study and ImmunoChip data from >19,000 individuals and methodologies developed from meta-analysis, we have identified opposing risk alleles at shared loci as well as independent disease-specific loci within the epidermal differentiation complex (chromosome 1q21.3), the Th2 locus control region (chromosome 5q31.1), and the major histocompatibility complex (chromosome 6p21-22). We further identified previously unreported pleiotropic alleles with opposing effects on atopic dermatitis and psoriasis risk in PRKRA and ANXA6/TNIP1. In contrast, there was no evidence for shared loci with effects operating in the same direction on both diseases. Our results show that atopic dermatitis and psoriasis have distinct genetic mechanisms with opposing effects in shared pathways influencing epidermal differentiation and immune response. The statistical analysis methods developed in the conduct of this study have produced additional insight from previously published data sets. The approach is likely to be applicable to the investigation of the genetic basis of other complex traits with overlapping and distinct clinical features.

Effects of MiR-137 genetic risk score on brain volume and cortical measures in patients with schizophrenia and controls.

Multiple genome-wide association studies of schizophrenia have implicated genetic variants within the gene encoding microRNA-137. As risk variants within or regulated by MIR137 have been implicated in memory performance, we investigated the additive effects of schizophrenia-associated risk variants in genes empirically regulated by MIR137 on brain regions associated with memory function. A polygenic risk score (PRS) was calculated (at a p = 0.05 threshold), using this empirically regulated MIR137 gene set, to investigate associations between this PRS and structural brain measures. These measures included total brain volume, cortical thickness, cortical surface area, and hippocampal volume, in a sample of 216 individuals consisting of healthy participants (n = 171) and patients with psychosis (n = 45). We did not observe a significant association between MIR137 PRS and these cortical thickness, surface area or hippocampal volume measures linked to memory function; a significant association between increasing PRS and decreasing total brain volume, independent of diagnosis status (R2 = 0.008, Beta = -0.09, p = 0.029), was observed. This did not survive correction for multiple testing. In conclusion, our study yielded only suggestive evidence that risk variants interacting with MIR137 impacts on cortical structure.

Duplications of the neuropeptide receptor gene VIPR2 confer significant risk for schizophrenia.

Rare copy number variants (CNVs) have a prominent role in the aetiology of schizophrenia and other neuropsychiatric disorders. Substantial risk for schizophrenia is conferred by large (>500-kilobase) CNVs at several loci, including microdeletions at 1q21.1 (ref. 2), 3q29 (ref. 3), 15q13.3 (ref. 2) and 22q11.2 (ref. 4) and microduplication at 16p11.2 (ref. 5). However, these CNVs collectively account for a small fraction (2-4%) of cases, and the relevant genes and neurobiological mechanisms are not well understood. Here we performed a large two-stage genome-wide scan of rare CNVs and report the significant association of copy number gains at chromosome 7q36.3 with schizophrenia. Microduplications with variable breakpoints occurred within a 362-kilobase region and were detected in 29 of 8,290 (0.35%) patients versus 2 of 7,431 (0.03%) controls in the combined sample. All duplications overlapped or were located within 89 kilobases upstream of the vasoactive intestinal peptide receptor gene VIPR2. VIPR2 transcription and cyclic-AMP signalling were significantly increased in cultured lymphocytes from patients with microduplications of 7q36.3. These findings implicate altered vasoactive intestinal peptide signalling in the pathogenesis of schizophrenia and indicate the VPAC2 receptor as a potential target for the development of new antipsychotic drugs.

Identification of loci associated with schizophrenia by genome-wide association and follow-up.

We carried out a genome-wide association study of schizophrenia (479 cases, 2,937 controls) and tested loci with P < 10(-5) in up to 16,726 additional subjects. Of 12 loci followed up, 3 had strong independent support (P < 5 x 10(-4)), and the overall pattern of replication was unlikely to occur by chance (P = 9 x 10(-8)). Meta-analysis provided strongest evidence for association around ZNF804A (P = 1.61 x 10(-7)) and this strengthened when the affected phenotype included bipolar disorder (P = 9.96 x 10(-9)).