More than 25 years of genetic studies of clozapine-induced agranulocytosis.

Clozapine is one of the most effective atypical antipsychotic drugs prescribed to patients with treatment-resistant schizophrenia. Approximately 1% of patients experience potential life-threatening adverse effects in the form of agranulocytosis, greatly hindering its applicability in clinical practice. The etiology of clozapine-induced agranulocytosis (CIA) remains unclear, but is thought to be a heritable trait. We reviewed the genetic studies of CIA published thus far. One recurrent finding from early candidate gene study to more recent genome-wide analysis is that of the involvement of human leukocyte antigen locus. We conclude that CIA is most likely a complex, polygenic trait, which may hamper efforts to the development of a genetic predictor test with clinical relevance. To decipher the genetic architecture of CIA, it is necessary to apply more rigorous standards of phenotyping and study much larger sample sizes.

Subcortical brain volume abnormalities in 2028 individuals with schizophrenia and 2540 healthy controls via the ENIGMA consortium.

The profile of brain structural abnormalities in schizophrenia is still not fully understood, despite decades of research using brain scans. To validate a prospective meta-analysis approach to analyzing multicenter neuroimaging data, we analyzed brain MRI scans from 2028 schizophrenia patients and 2540 healthy controls, assessed with standardized methods at 15 centers worldwide. We identified subcortical brain volumes that differentiated patients from controls, and ranked them according to their effect sizes. Compared with healthy controls, patients with schizophrenia had smaller hippocampus (Cohen's d=-0.46), amygdala (d=-0.31), thalamus (d=-0.31), accumbens (d=-0.25) and intracranial volumes (d=-0.12), as well as larger pallidum (d=0.21) and lateral ventricle volumes (d=0.37). Putamen and pallidum volume augmentations were positively associated with duration of illness and hippocampal deficits scaled with the proportion of unmedicated patients. Worldwide cooperative analyses of brain imaging data support a profile of subcortical abnormalities in schizophrenia, which is consistent with that based on traditional meta-analytic approaches. This first ENIGMA Schizophrenia Working Group study validates that collaborative data analyses can readily be used across brain phenotypes and disorders and encourages analysis and data sharing efforts to further our understanding of severe mental illness.

High educational performance is a distinctive feature of bipolar disorder: a study on cognition in bipolar disorder, schizophrenia patients, relatives and controls.

BACKGROUND: Schizophrenia is associated with lower intelligence and poor educational performance relative to the general population. This is, to a lesser degree, also found in first-degree relatives of schizophrenia patients. It is unclear whether bipolar disorder I (BD-I) patients and their relatives have similar lower intellectual and educational performance as that observed in schizophrenia. METHOD: This cross-sectional study investigated intelligence and educational performance in two outpatient samples [494 BD-I patients, 952 schizophrenia spectrum (SCZ) patients], 2231 relatives of BD-I and SCZ patients, 1104 healthy controls and 100 control siblings. Mixed-effects and regression models were used to compare groups on intelligence and educational performance. RESULTS: BD-I patients were more likely to have completed the highest level of education (odds ratio 1.88, 95% confidence interval 1.66-2.70) despite having a lower IQ compared to controls (ß = -9.09, S.E. = 1.27, p < 0.001). In contrast, SCZ patients showed both a lower IQ (ß = -15.31, S.E. = 0.86, p < 0.001) and lower educational levels compared to controls. Siblings of both patient groups had significantly lower IQ than control siblings, but did not differ on educational performance. IQ scores did not differ between BD-I parents and SCZ parents, but BD-I parents had completed higher educational levels. CONCLUSIONS: Although BD-I patients had a lower IQ than controls, they were more likely to have completed the highest level of education. This contrasts with SCZ patients, who showed both intellectual and educational deficits compared to healthy controls. Since relatives of BD-I patients did not demonstrate superior educational performance, our data suggest that high educational performance may be a distinctive feature of bipolar disorder patients.

Evaluating historical candidate genes for schizophrenia.

Prior to the genome-wide association era, candidate gene studies were a major approach in schizophrenia genetics. In this invited review, we consider the current status of 25 historical candidate genes for schizophrenia (for example, COMT, DISC1, DTNBP1 and NRG1). The initial study for 24 of these genes explicitly evaluated common variant hypotheses about schizophrenia. Our evaluation included a meta-analysis of the candidate gene literature, incorporation of the results of the largest genomic study yet published for schizophrenia, ratings from informed researchers who have published on these genes, and ratings from 24 schizophrenia geneticists. On the basis of current empirical evidence and mostly consensual assessments of informed opinion, it appears that the historical candidate gene literature did not yield clear insights into the genetic basis of schizophrenia. A likely reason why historical candidate gene studies did not achieve their primary aims is inadequate statistical power. However, the considerable efforts embodied in these early studies unquestionably set the stage for current successes in genomic approaches to schizophrenia.

Genome-wide association study of NMDA receptor coagonists in human cerebrospinal fluid and plasma.

The N-methyl-D-aspartate receptor (NMDAR) coagonists glycine, D-serine and L-proline play crucial roles in NMDAR-dependent neurotransmission and are associated with a range of neuropsychiatric disorders. We conducted the first genome-wide association study of concentrations of these coagonists and their enantiomers in plasma and cerebrospinal fluid (CSF) of human subjects from the general population (N=414). Genetic variants at chromosome 22q11.2, located in and near PRODH (proline dehydrogenase), were associated with L-proline in plasma (ß=0.29; P=6.38 × 10(-10)). The missense variant rs17279437 in the proline transporter SLC6A20 was associated with L-proline in CSF (ß=0.28; P=9.68 × 10(-9)). Suggestive evidence of association was found for the D-serine plasma-CSF ratio at the D-amino-acid oxidase (DAO) gene (ß=-0.28; P=9.08 × 10(-8)), whereas a variant in SRR (that encodes serine racemase and is associated with schizophrenia) constituted the most strongly associated locus for the L-serine to D-serine ratio in CSF. All these genes are highly expressed in rodent meninges and choroid plexus, anatomical regions relevant to CSF physiology. The enzymes and transporters they encode may be targeted to further construe the nature of NMDAR coagonist involvement in NMDAR gating. Furthermore, the highlighted genetic variants may be followed up in clinical populations, for example, schizophrenia and 22q11 deletion syndrome. Overall, this targeted metabolomics approach furthers the understanding of NMDAR coagonist concentration variability and sets the stage for non-targeted CSF metabolomics projects.

A genome-wide association study of rheumatoid arthritis without antibodies against citrullinated peptides.

INTRODUCTION: Rheumatoid arthritis (RA) patients can be classified based on presence or absence of anticitrullinated peptide antibodies (ACPA) in their serum. This heterogeneity among patients may reflect important biological differences underlying the disease process. To date, the majority of genetic studies have focused on the ACPA-positive group. Therefore, our goal was to analyse the genetic risk factors that contribute to ACPA-negative RA. METHODS: We performed a large-scale genome-wide association study (GWAS) in three Caucasian European cohorts comprising 1148 ACPA-negative RA patients and 6008 controls. All patients were screened using the Illumina Human Cyto-12 chip, and controls were genotyped using different genome-wide platforms. Population-independent analyses were carried out by means of logistic regression. Meta-analysis with previously published data was performed as follow-up for selected signals (reaching a total of 1922 ACPA-negative RA patients and 7087 controls). Imputation of classical HLA alleles, amino acid residues and single nucleotide polymorphisms was undertaken. RESULTS: The combined analysis of the studied cohorts resulted in identification of a peak of association in the HLA-region and several suggestive non-HLA associations. Meta-analysis with previous reports confirmed the association of the HLA region with this subset and an observed association in the CLYBL locus remained suggestive. The imputation and deep interrogation of the HLA region led to identification of a two amino acid model (HLA-B at position 9 and HLA-DRB1 at position 11) that accounted for the observed genome-wide associations in this region. CONCLUSIONS: Our study shed light on the influence of the HLA region in ACPA-negative RA and identified a suggestive risk locus for this condition.

Genome-wide study of association and interaction with maternal cytomegalovirus infection suggests new schizophrenia loci.

Genetic and environmental components as well as their interaction contribute to the risk of schizophrenia, making it highly relevant to include environmental factors in genetic studies of schizophrenia. This study comprises genome-wide association (GWA) and follow-up analyses of all individuals born in Denmark since 1981 and diagnosed with schizophrenia as well as controls from the same birth cohort. Furthermore, we present the first genome-wide interaction survey of single nucleotide polymorphisms (SNPs) and maternal cytomegalovirus (CMV) infection. The GWA analysis included 888 cases and 882 controls, and the follow-up investigation of the top GWA results was performed in independent Danish (1396 cases and 1803 controls) and German-Dutch (1169 cases, 3714 controls) samples. The SNPs most strongly associated in the single-marker analysis of the combined Danish samples were rs4757144 in ARNTL (P=3.78 × 10(-6)) and rs8057927 in CDH13 (P=1.39 × 10(-5)). Both genes have previously been linked to schizophrenia or other psychiatric disorders. The strongest associated SNP in the combined analysis, including Danish and German-Dutch samples, was rs12922317 in RUNDC2A (P=9.04 × 10(-7)). A region-based analysis summarizing independent signals in segments of 100 kb identified a new region-based genome-wide significant locus overlapping the gene ZEB1 (P=7.0 × 10(-7)). This signal was replicated in the follow-up analysis (P=2.3 × 10(-2)). Significant interaction with maternal CMV infection was found for rs7902091 (P(SNP × CMV)=7.3 × 10(-7)) in CTNNA3, a gene not previously implicated in schizophrenia, stressing the importance of including environmental factors in genetic studies.

Genome-wide association study of monoamine metabolite levels in human cerebrospinal fluid.

Studying genetic determinants of intermediate phenotypes is a powerful tool to increase our understanding of genotype-phenotype correlations. Metabolic traits pertinent to the central nervous system (CNS) constitute a potentially informative target for genetic studies of intermediate phenotypes as their genetic underpinnings may elucidate etiological mechanisms. We therefore conducted a genome-wide association study (GWAS) of monoamine metabolite (MM) levels in cerebrospinal fluid (CSF) of 414 human subjects from the general population. In a linear model correcting for covariates, we identified one locus associated with MMs at a genome-wide significant level (standardized ß=0.32, P=4.92 × 10(-8)), located 20 kb from SSTR1, a gene involved with brain signal transduction and glutamate receptor signaling. By subsequent whole-genome expression quantitative trait locus (eQTL) analysis, we provide evidence that this variant controls expression of PDE9A (ß=0.21; P unadjusted=5.6 × 10(-7); P corrected=0.014), a gene previously implicated in monoaminergic transmission, major depressive disorder and antidepressant response. A post hoc analysis of loci significantly associated with psychiatric disorders suggested that genetic variation at CSMD1, a schizophrenia susceptibility locus, plays a role in the ratio between dopamine and serotonin metabolites in CSF. The presented DNA and mRNA analyses yielded genome-wide and suggestive associations in biologically plausible genes, two of which encode proteins involved with glutamate receptor functionality. These findings will hopefully contribute to an exploration of the functional impact of the highlighted genes on monoaminergic transmission and neuropsychiatric phenotypes.

Genome-wide association study of lymphoblast cell viability after clozapine exposure.

Clozapine is an antipsychotic drug with proven efficacy in treatment-resistant schizophrenia but also known to induce potentially lethal agranulocytosis (CIA) in 1% of patients. Genetic factors are likely to play a role in the molecular basis of CIA. We explored an in vitro system to study the genetic susceptibility of CIA. Cell viability was measured in 90 lymphoblast cell lines exposed to a series of increasing concentrations of clozapine for 48 hr. Quantitative trait measures of cell viability as well as area under the survival curve were used in a linear mixed model for genome-wide association analyses. The estimated heritability of clozapine-induced cell viability reduction in these cell lines is h2=0.76. No genome-wide significant association was observed after correction for multiple testing. Two independent loci with nominal evidence of association were observed at 30× clinical clozapine concentration: rs2709505 (P=1.41×10(-8)) in an intron of MDFIC and rs10457252 (P=1.79×10(-8)) located in a gene desert at chromosome 6q21. We identified one locus (rs1293970) near PRG4 that was consistently associated for all separate concentration analyses at P<5×10(-5). PRG4 encodes hemangiopoietin, a growth stimulator for hematopoietic stem cells. No evidence was observed for involvement of the MHC region. Our results demonstrate that clozapine-induced viability reduction in lymphoblast cell lines is a heritable, polygenic trait. Thus, in vitro models of CIA might be a useful tool for future discovery of genetic risk factors, although larger sample sizes will be required to unambiguously identify these loci.

Schizophrenia genetic variants are not associated with intelligence.

BACKGROUND: Schizophrenia is associated with lower pre-morbid intelligence (IQ) in addition to (pre-morbid) cognitive decline. Both schizophrenia and IQ are highly heritable traits. Therefore, we hypothesized that genetic variants associated with schizophrenia, including copy number variants (CNVs) and a polygenic schizophrenia (risk) score (PSS), may influence intelligence. METHOD: IQ was estimated with the Wechsler Adult Intelligence Scale (WAIS). CNVs were determined from single nucleotide polymorphism (SNP) data using the QuantiSNP and PennCNV algorithms. For the PSS, odds ratios for genome-wide SNP data were calculated in a sample collected by the Psychiatric Genome-Wide Association Study (GWAS) Consortium (8690 schizophrenia patients and 11 831 controls). These were used to calculate individual PSSs in our independent sample of 350 schizophrenia patients and 322 healthy controls. RESULTS: Although significantly more genes were disrupted by deletions in schizophrenia patients compared to controls (p = 0.009), there was no effect of CNV measures on IQ. The PSS was associated with disease status (R 2 = 0.055, p = 2.1 × 10-7) and with IQ in the entire sample (R 2 = 0.018, p = 0.0008) but the effect on IQ disappeared after correction for disease status. CONCLUSIONS: Our data suggest that rare and common schizophrenia-associated variants do not explain the variation in IQ in healthy subjects or in schizophrenia patients. Thus, reductions in IQ in schizophrenia patients may be secondary to other processes related to schizophrenia risk.

A double hit implicates DIAPH3 as an autism risk gene.

Recent studies have shown that more than 10% of autism cases are caused by de novo structural genomic rearrangements. Given that some heritable copy number variants (CNVs) have been observed in patients as well as in healthy controls, to date little attention has been paid to the potential function of these non-de novo CNVs in causing autism. A normally intelligent patient with autism, with non-affected parents, was identified with a maternally inherited 10 Mb deletion at 13q21.2. Sequencing of the genes within the deletion identified a paternally inherited nonsynonymous amino-acid substitution at position 614 of diaphanous homolog 3 (DIAPH3) (proline to threonine; Pro614Thr). This variant, present in a highly conserved domain, was not found in 328 healthy subjects. Experiments showed a transient expression of Diaph3 in the developing murine cerebral cortex, indicating it has a function in brain development. Transfection of Pro614Thr in murine fibroblasts showed a significant reduction in the number of induced filopodia in comparison to the wild-type gene. DIAPH3 is involved in cell migration, axon guidance and neuritogenesis, and is suggested to function downstream of SHANK3. Our findings strongly suggest DIAPH3 as a novel autism susceptibility gene. Moreover, this report of a 'double-hit' compound heterozygote for a large, maternally inherited, genomic deletion and a paternally inherited rare missense mutation shows that not only de novo genomic variants in patients should be taken seriously in further study but that inherited CNVs may also provide valuable information.

Copy number variations of chromosome 16p13.1 region associated with schizophrenia.

Deletions and reciprocal duplications of the chromosome 16p13.1 region have recently been reported in several cases of autism and mental retardation (MR). As genomic copy number variants found in these two disorders may also associate with schizophrenia, we examined 4345 schizophrenia patients and 35,079 controls from 8 European populations for duplications and deletions at the 16p13.1 locus, using microarray data. We found a threefold excess of duplications and deletions in schizophrenia cases compared with controls, with duplications present in 0.30% of cases versus 0.09% of controls (P=0.007) and deletions in 0.12 % of cases and 0.04% of controls (P>0.05). The region can be divided into three intervals defined by flanking low copy repeats. Duplications spanning intervals I and II showed the most significant (P = 0.00010) association with schizophrenia. The age of onset in duplication and deletion carriers among cases ranged from 12 to 35 years, and the majority were males with a family history of psychiatric disorders. In a single Icelandic family, a duplication spanning intervals I and II was present in two cases of schizophrenia, and individual cases of alcoholism, attention deficit hyperactivity disorder and dyslexia. Candidate genes in the region include NTAN1 and NDE1. We conclude that duplications and perhaps also deletions of chromosome 16p13.1, previously reported to be associated with autism and MR, also confer risk of schizophrenia.

Expanding the range of ZNF804A variants conferring risk of psychosis.

A trio of genome-wide association studies recently reported sequence variants at three loci to be significantly associated with schizophrenia. No sequence polymorphism had been unequivocally (P<5 × 10(-8)) associated with schizophrenia earlier. However, one variant, rs1344706[T], had come very close. This polymorphism, located in an intron of ZNF804A, was reported to associate with schizophrenia with a P-value of 1.6 × 10(-7), and with psychosis (schizophrenia plus bipolar disorder) with a P-value of 1.0 × 10(-8). In this study, using 5164 schizophrenia cases and 20,709 controls, we replicated the association with schizophrenia (odds ratio OR = 1.08, P = 0.0029) and, by adding bipolar disorder patients, we also confirmed the association with psychosis (added N = 609, OR = 1.09, P = 0.00065). Furthermore, as it has been proposed that variants such as rs1344706[T]-common and with low relative risk-may also serve to identify regions harboring less common, higher-risk susceptibility alleles, we searched ZNF804A for large copy number variants (CNVs) in 4235 psychosis patients, 1173 patients with other psychiatric disorders and 39,481 controls. We identified two CNVs including at least part of ZNF804A in psychosis patients and no ZNF804A CNVs in controls (P = 0.013 for association with psychosis). In addition, we found a ZNF804A CNV in an anxiety patient (P = 0.0016 for association with the larger set of psychiatric disorders).

Genome arrays for the detection of copy number variations in idiopathic mental retardation, idiopathic generalized epilepsy and neuropsychiatric disorders: lessons for diagnostic workflow and research.

We review the contributions and limitations of genome-wide array-based identification of copy number variants (CNVs) in the clinical diagnostic evaluation of patients with mental retardation (MR) and other brain-related disorders. In unselected MR referrals a causative genomic gain or loss is detected in 14-18% of cases. Usually, such CNVs arise de novo, are not found in healthy subjects, and have a major impact on the phenotype by altering the dosage of multiple genes. This high diagnostic yield justifies array-based segmental aneuploidy screening as the initial genetic test in these patients. This also pertains to patients with autism (expected yield about 5-10% in nonsyndromic and 10-20% in syndromic patients) and schizophrenia (at least 5% yield). CNV studies in idiopathic generalized epilepsy, attention-deficit hyperactivity disorder, major depressive disorder and Tourette syndrome indicate that patients have, on average, a larger CNV burden as compared to controls. Collectively, the CNV studies suggest that a wide spectrum of disease-susceptibility variants exists, most of which are rare (<0.1%) and of variable and usually small effect. Notwithstanding, a rare CNV can have a major impact on the phenotype. Exome sequencing in MR and autism patients revealed de novo mutations in protein coding genes in 60 and 20% of cases, respectively. Therefore, it is likely that arrays will be supplanted by next-generation sequencing methods as the initial and perhaps ultimate diagnostic tool in patients with brain-related disorders, revealing both CNVs and mutations in a single test.

Cannabis use at a young age is associated with psychotic experiences.

BACKGROUND: Cannabis use is associated with psychosis and a range of subclinical psychiatric symptoms. The strength of this association depends on dosage and age at first use. The current study investigates whether level of cannabis exposure and starting age are associated with specific profiles of subclinical symptoms. METHOD: We collected cross-sectional data from a young adult population sample by administering an online version of the Community Assessment of Psychic Experiences (CAPE). Cannabis exposure was quantified as the amount of Euros spent on cannabis per week and the age of initial cannabis use. The primary outcome measure was the odds ratio (OR) to belong to the highest 10% of scores on the total CAPE and the positive-, negative- and depressive symptom dimensions. RESULTS: In 17 698 adolescents (mean age 21.6, s.d.=4.2 years), cannabis use at age 12 years or younger was strongly associated with a top 10% score on psychotic experiences [OR 3.1, 95% confidence interval (CI) 2.1-4.3] and to a lesser degree with negative symptoms (OR 1.7, 95% CI 1.1-2.5). The OR of heavy users (>€25/week) for negative symptoms was 3.4 (95% CI 2.9-4.1), for psychotic experiences 3.0 (95% CI 2.4-3.6), and for depressive symptoms 2.8 (95% CI 2.3-3.3). CONCLUSIONS: Early start of cannabis use is strongly associated with subclinical psychotic symptoms and to a lesser degree with negative symptoms, while smoking high amounts of cannabis is associated with increased levels of all three symptom dimensions: psychotic, negative and depressive. These results support the hypothesis that the impact of cannabis use is age specific.

Association between cannabis and psychiatric hospitalization.

OBJECTIVE: To investigate the relationship between cannabis use and mental health. METHOD: A cross-sectional analysis in a sample of 17 698 individuals with a mean age of 22 years (SD: 4.2). Participants provided information on the amount and initial age of cannabis use and history of psychiatric hospitalizations through a web-based questionnaire. To quantify Δ(9) -tetrahydrocannabinol exposure, we operationalized cannabis use as the amount of money spent on cannabis per week over the last month. The odds ratio of having a history of psychiatric hospitalizations was the primary outcome measure. RESULTS: We found a dose-response relationship between the amount of cannabis use and the odds for psychiatric hospitalization. Adjusted odds ratios for hospitalization increased with the amount of cannabis consumed from 1.6 (95% CI: 1.1-2.3) in incidental users to 6.2 (95% CI: 4.3-8.9) in heavy users (>€25/week). Our data suggested that concomitant drug use was an intermediate factor. Exposure to cannabis before the age of 12 years was found to carry a 4.8 (95% CI: 2.9-7.8) times increased odds for past psychiatric hospitalizations. CONCLUSION: We conclude that early and heavy uses of cannabis are each and independently associated with poor mental health in its users.

Recurrent inversion events at 17q21.31 microdeletion locus are linked to the MAPT H2 haplotype.

The chromosomal band 17q21.31, containing the microtubule-associated protein tau (MAPT) gene, is a hotspot for chromosomal rearrangements. It is known to contain a common inversion polymorphism of approximately 900 kb in populations with European ancestry. The inverted configuration is linked to a distinct MAPT haplotype, H2, which is relatively common in Europeans but nearly absent in Asian and African populations. Recent studies have demonstrated that the H2 haplotype is ancestral in hominoids, and under positive selection in Europeans. This haplotype is also linked to events leading to the 17q21.31 microdeletion syndrome, one of the most common causes of 'idiopathic' mental retardation in people of European descent. We performed direct analysis of the chromosome structure by fluorescence in situ hybridization and observed heterozygosity of the inversion status for the H2 chromosomes, but not for the H1 haplotype. Inversion heterozygosity was also observed in a mother homozygous for the H2 haplotype, who transmitted the chromosome with the deletion to a proband with 17q21.31 microdeletion syndrome. Our results highlight an allele-specific sensitivity to chromosome rearrangements and suggest that it is the heterozygosity of inversion status that predisposes to the 17q21.31 microdeletion syndrome.

"Forward genetics" as a method to maximize power and cost-efficiency in studies of human complex traits.

There is increasing interest in methods to disentangle the relationship between genotype and (endo)phenotypes in human complex traits. We present a population-based method of increasing the power and cost-efficiency of studies by selecting random individuals with a particular genotype and then assessing the accompanying quantitative phenotypes. Using statistical derivations, power- and cost graphs we show that such a "forward genetics" approach can lead to a marked reduction in sample size and costs. This approach is particularly apt for implementing in epidemiological studies for which DNA is already available but the phenotyping costs are high.

No association between anti-thyroidperoxidase antibodies and bipolar disorder: a study in the Dutch Bipolar Cohort and a meta-analysis.

BACKGROUND: Thyroid autoimmunity has been associated with bipolar disorder (BD). However, results from previous studies on the seroprevalence of anti-thyroid peroxidase antibodies (TPO-abs) in BD are inconsistent. OBJECTIVES: The aim of the present study is to investigate whether the seroprevalence and titer levels of TPO-abs are related to BD. METHOD: TPO-abs were measured in plasma samples of 760 patients with bipolar disorder, 261 first-degree relatives and 363 controls by enzyme-linked immunosorbent assay (ELISA). To address methodological limitations of previous studies, we assessed clinical characteristics with several (self-reported) questionnaires to investigate whether TPO-abs positivity is related to particular clinical subgroups of BD patients. We performed an additional meta-analysis of seroprevalences of TPO-abs in BD patients including data from present and previous studies. RESULTS: Seroprevalence or titer levels of TPO-abs did not significantly differ between patients with BD, their first-degree relatives, and controls. In BD patients, the prevalence of TPO-abs was unrelated to specific clinical factors, including lithium use. Our meta-analysis of twelve studies showed an overall odds ratio of 1.3 (CI 95 %: 0.7-2.3; p = 0.30), reaffirming the absence of an association of BD with TPO-abs. CONCLUSIONS: In the largest study of TPO-abs in BD to date, our findings indicate that TPO-abs are not associated with (the risk for) bipolar disorder.