High-dimensional immune profiling identifies a biomarker to monitor dimethyl fumarate response in multiple sclerosis.

Dimethyl fumarate (DMF) is an immunomodulatory treatment for multiple sclerosis (MS). Despite its wide clinical use, the mechanisms underlying clinical response are not understood. This study aimed to reveal immune markers of therapeutic response to DMF treatment in MS. For this purpose, we prospectively collected peripheral blood mononuclear cells (PBMCs) from a highly characterized cohort of 44 individuals with MS before and at 12 and 48 wk of DMF treatment. Single cells were profiled using high-dimensional mass cytometry. To capture the heterogeneity of different immune subsets, we adopted a bioinformatic multipanel approach that allowed cell population-cluster assignment of more than 50 different parameters, including lineage and activation markers as well as chemokine receptors and cytokines. Data were further analyzed in a semiunbiased fashion implementing a supervised representation learning approach to capture subtle longitudinal immune changes characteristic for therapy response. With this approach, we identified a population of memory T helper cells expressing high levels of neuroinflammatory cytokines (granulocyte-macrophage colony-stimulating factor [GM-CSF], interferon γ [IFNγ]) as well as CXCR3, whose abundance correlated with treatment response. Using spectral flow cytometry, we confirmed these findings in a second cohort of patients. Serum neurofilament light-chain levels confirmed the correlation of this immune cell signature with axonal damage. The identified cell population is expanded in peripheral blood under natalizumab treatment, substantiating a specific role in treatment response. We propose that depletion of GM-CSF-, IFNγ-, and CXCR3-expressing T helper cells is the main mechanism of action of DMF and allows monitoring of treatment response.

Association of polygenic score and the involvement of cholinergic and glutamatergic pathways with lithium treatment response in patients with bipolar disorder.

Lithium is regarded as the first-line treatment for bipolar disorder (BD), a severe and disabling mental health disorder that affects about 1% of the population worldwide. Nevertheless, lithium is not consistently effective, with only 30% of patients showing a favorable response to treatment. To provide personalized treatment options for bipolar patients, it is essential to identify prediction biomarkers such as polygenic scores. In this study, we developed a polygenic score for lithium treatment response (Li+PGS) in patients with BD. To gain further insights into lithium's possible molecular mechanism of action, we performed a genome-wide gene-based analysis. Using polygenic score modeling, via methods incorporating Bayesian regression and continuous shrinkage priors, Li+PGS was developed in the International Consortium of Lithium Genetics cohort (ConLi+Gen: N = 2367) and replicated in the combined PsyCourse (N = 89) and BipoLife (N = 102) studies. The associations of Li+PGS and lithium treatment response - defined in a continuous ALDA scale and a categorical outcome (good response vs. poor response) were tested using regression models, each adjusted for the covariates: age, sex, and the first four genetic principal components. Statistical significance was determined at P < 0.05. Li+PGS was positively associated with lithium treatment response in the ConLi+Gen cohort, in both the categorical (P = 9.8 × 10-12, R2 = 1.9%) and continuous (P = 6.4 × 10-9, R2 = 2.6%) outcomes. Compared to bipolar patients in the 1st decile of the risk distribution, individuals in the 10th decile had 3.47-fold (95%CI: 2.22-5.47) higher odds of responding favorably to lithium. The results were replicated in the independent cohorts for the categorical treatment outcome (P = 3.9 × 10-4, R2 = 0.9%), but not for the continuous outcome (P = 0.13). Gene-based analyses revealed 36 candidate genes that are enriched in biological pathways controlled by glutamate and acetylcholine. Li+PGS may be useful in the development of pharmacogenomic testing strategies by enabling a classification of bipolar patients according to their response to treatment.

Relationships between neurotransmitter receptor densities and expression levels of their corresponding genes in the human hippocampus.

Neurotransmitter receptors are key molecules in signal transmission, their alterations are associated with brain dysfunction. Relationships between receptors and their corresponding genes are poorly understood, especially in humans. We combined in vitro receptor autoradiography and RNA sequencing to quantify, in the same tissue samples (7 subjects), the densities of 14 receptors and expression levels of their corresponding 43 genes in the Cornu Ammonis (CA) and dentate gyrus (DG) of human hippocampus. Significant differences in receptor densities between both structures were found only for metabotropic receptors, whereas significant differences in RNA expression levels mostly pertained ionotropic receptors. Receptor fingerprints of CA and DG differ in shapes but have similar sizes; the opposite holds true for their "RNA fingerprints", which represent the expression levels of multiple genes in a single area. In addition, the correlation coefficients between receptor densities and corresponding gene expression levels vary widely and the mean correlation strength was weak-to-moderate. Our results suggest that receptor densities are not only controlled by corresponding RNA expression levels, but also by multiple regionally specific post-translational factors.

Combined analysis of cytoarchitectonic, molecular and transcriptomic patterns reveal differences in brain organization across human functional brain systems.

Brain areas show specific cellular, molecular, and gene expression patterns that are linked to function, but their precise relationships are largely unknown. To unravel these structure-function relationships, a combined analysis of 53 neurotransmitter receptor genes, receptor densities of six transmitter systems and cytoarchitectonic data of the auditory, somatosensory, visual, motor systems was conducted. Besides covariation of areal gene expression with receptor density, the study reveals specific gene expression patterns in functional systems, which are most prominent for the inhibitory GABAA and excitatory glutamatergic NMDA receptors. Furthermore, gene expression-receptor relationships changed in a systematic manner according to information flow from primary to higher associative areas. The findings shed new light on the relationship of anatomical, functional, and molecular and transcriptomic principles of cortical segregation towards a more comprehensive understanding of human brain organization.

Effects of copy number variations on brain structure and risk for psychiatric illness: Large-scale studies from the ENIGMA working groups on CNVs.

The Enhancing NeuroImaging Genetics through Meta-Analysis copy number variant (ENIGMA-CNV) and 22q11.2 Deletion Syndrome Working Groups (22q-ENIGMA WGs) were created to gain insight into the involvement of genetic factors in human brain development and related cognitive, psychiatric and behavioral manifestations. To that end, the ENIGMA-CNV WG has collated CNV and magnetic resonance imaging (MRI) data from ~49,000 individuals across 38 global research sites, yielding one of the largest studies to date on the effects of CNVs on brain structures in the general population. The 22q-ENIGMA WG includes 12 international research centers that assessed over 533 individuals with a confirmed 22q11.2 deletion syndrome, 40 with 22q11.2 duplications, and 333 typically developing controls, creating the largest-ever 22q11.2 CNV neuroimaging data set. In this review, we outline the ENIGMA infrastructure and procedures for multi-site analysis of CNVs and MRI data. So far, ENIGMA has identified effects of the 22q11.2, 16p11.2 distal, 15q11.2, and 1q21.1 distal CNVs on subcortical and cortical brain structures. Each CNV is associated with differences in cognitive, neurodevelopmental and neuropsychiatric traits, with characteristic patterns of brain structural abnormalities. Evidence of gene-dosage effects on distinct brain regions also emerged, providing further insight into genotype-phenotype relationships. Taken together, these results offer a more comprehensive picture of molecular mechanisms involved in typical and atypical brain development. This "genotype-first" approach also contributes to our understanding of the etiopathogenesis of brain disorders. Finally, we outline future directions to better understand effects of CNVs on brain structure and behavior.

Using polygenic scores and clinical data for bipolar disorder patient stratification and lithium response prediction: machine learning approach.

BACKGROUND: Response to lithium in patients with bipolar disorder is associated with clinical and transdiagnostic genetic factors. The predictive combination of these variables might help clinicians better predict which patients will respond to lithium treatment. AIMS: To use a combination of transdiagnostic genetic and clinical factors to predict lithium response in patients with bipolar disorder. METHOD: This study utilised genetic and clinical data (n = 1034) collected as part of the International Consortium on Lithium Genetics (ConLi+Gen) project. Polygenic risk scores (PRS) were computed for schizophrenia and major depressive disorder, and then combined with clinical variables using a cross-validated machine-learning regression approach. Unimodal, multimodal and genetically stratified models were trained and validated using ridge, elastic net and random forest regression on 692 patients with bipolar disorder from ten study sites using leave-site-out cross-validation. All models were then tested on an independent test set of 342 patients. The best performing models were then tested in a classification framework. RESULTS: The best performing linear model explained 5.1% (P = 0.0001) of variance in lithium response and was composed of clinical variables, PRS variables and interaction terms between them. The best performing non-linear model used only clinical variables and explained 8.1% (P = 0.0001) of variance in lithium response. A priori genomic stratification improved non-linear model performance to 13.7% (P = 0.0001) and improved the binary classification of lithium response. This model stratified patients based on their meta-polygenic loadings for major depressive disorder and schizophrenia and was then trained using clinical data. CONCLUSIONS: Using PRS to first stratify patients genetically and then train machine-learning models with clinical predictors led to large improvements in lithium response prediction. When used with other PRS and biological markers in the future this approach may help inform which patients are most likely to respond to lithium treatment.

Bipolar multiplex families have an increased burden of common risk variants for psychiatric disorders.

Multiplex families with a high prevalence of a psychiatric disorder are often examined to identify rare genetic variants with large effect sizes. In the present study, we analysed whether the risk for bipolar disorder (BD) in BD multiplex families is influenced by common genetic variants. Furthermore, we investigated whether this risk is conferred mainly by BD-specific risk variants or by variants also associated with the susceptibility to schizophrenia or major depression. In total, 395 individuals from 33 Andalusian BD multiplex families (166 BD, 78 major depressive disorder, 151 unaffected) as well as 438 subjects from an independent, BD case/control cohort (161 unrelated BD, 277 unrelated controls) were analysed. Polygenic risk scores (PRS) for BD, schizophrenia (SCZ), and major depression were calculated and compared between the cohorts. Both the familial BD cases and unaffected family members had higher PRS for all three psychiatric disorders than the independent controls, with BD and SCZ being significant after correction for multiple testing, suggesting a high baseline risk for several psychiatric disorders in the families. Moreover, familial BD cases showed significantly higher BD PRS than unaffected family members and unrelated BD cases. A plausible hypothesis is that, in multiplex families with a general increase in risk for psychiatric disease, BD development is attributable to a high burden of common variants that confer a specific risk for BD. The present analyses demonstrated that common genetic risk variants for psychiatric disorders are likely to contribute to the high incidence of affective psychiatric disorders in the multiplex families. However, the PRS explained only part of the observed phenotypic variance, and rare variants might have also contributed to disease development.

Association of polygenic score for major depression with response to lithium in patients with bipolar disorder.

Lithium is a first-line medication for bipolar disorder (BD), but only one in three patients respond optimally to the drug. Since evidence shows a strong clinical and genetic overlap between depression and bipolar disorder, we investigated whether a polygenic susceptibility to major depression is associated with response to lithium treatment in patients with BD. Weighted polygenic scores (PGSs) were computed for major depression (MD) at different GWAS p value thresholds using genetic data obtained from 2586 bipolar patients who received lithium treatment and took part in the Consortium on Lithium Genetics (ConLi+Gen) study. Summary statistics from genome-wide association studies in MD (135,458 cases and 344,901 controls) from the Psychiatric Genomics Consortium (PGC) were used for PGS weighting. Response to lithium treatment was defined by continuous scores and categorical outcome (responders versus non-responders) using measurements on the Alda scale. Associations between PGSs of MD and lithium treatment response were assessed using a linear and binary logistic regression modeling for the continuous and categorical outcomes, respectively. The analysis was performed for the entire cohort, and for European and Asian sub-samples. The PGSs for MD were significantly associated with lithium treatment response in multi-ethnic, European or Asian populations, at various p value thresholds. Bipolar patients with a low polygenic load for MD were more likely to respond well to lithium, compared to those patients with high polygenic load [lowest vs highest PGS quartiles, multi-ethnic sample: OR = 1.54 (95% CI: 1.18-2.01) and European sample: OR = 1.75 (95% CI: 1.30-2.36)]. While our analysis in the Asian sample found equivalent effect size in the same direction: OR = 1.71 (95% CI: 0.61-4.90), this was not statistically significant. Using PGS decile comparison, we found a similar trend of association between a high genetic loading for MD and lower response to lithium. Our findings underscore the genetic contribution to lithium response in BD and support the emerging concept of a lithium-responsive biotype in BD.

Efficacy and safety of exogenous beta-hydroxybutyrate for preventive treatment in episodic migraine: A single-centred, randomised, placebo-controlled, double-blind crossover trial.

BACKGROUND: Several studies propose that brain energy deficit might be partially involved in the pathophysiology of migraine. Previously, studies demonstrated that ketogenic diet causes a substantial reduction in migraine frequency. Since the ketogenic diet is restricting and its adherence is difficult, we proposed to supplement ketone bodies exogenously to provide a prophylactic effect in migraineurs. AIM: To evaluate the prophylactic effect of exogenous DL-beta-hydroxybutyrate supplementation in episodic migraineurs. METHODS: A double-blind, placebo-controlled, randomised crossover trial was conducted, involving 41 patients with episodic migraine. Patients were randomised 1:1 into placebo or beta-hydroxybutyrate group before entering the first treatment period. Each treatment period was 12 weeks long, followed by four weeks of washout phase and four weeks of run-in phase before entering into the corresponding second treatment period. The primary endpoint was the number of migraine days in the last four weeks of treatment, adjusted for baseline. RESULTS: We observed no clinically significant amelioration of migraine frequency or intensity under DL-beta-hydroxybutyrate treatment as compared to placebo regarding number of migraine days (mean difference [95% CI]: -1.1[-5.07, 2.85]), migraine intensity (0-10 VAS: 1.5[-0.8, 3.7]). CONCLUSION: The selected dose of supplemented exogenous DL-beta-hydroxybutyrate did not demonstrate efficacy in episodic migraineurs.ClinicalTrials.gov Identifier: NCT03132233.

Characterisation of age and polarity at onset in bipolar disorder.

BACKGROUND: Studying phenotypic and genetic characteristics of age at onset (AAO) and polarity at onset (PAO) in bipolar disorder can provide new insights into disease pathology and facilitate the development of screening tools. AIMS: To examine the genetic architecture of AAO and PAO and their association with bipolar disorder disease characteristics. METHOD: Genome-wide association studies (GWASs) and polygenic score (PGS) analyses of AAO (n = 12 977) and PAO (n = 6773) were conducted in patients with bipolar disorder from 34 cohorts and a replication sample (n = 2237). The association of onset with disease characteristics was investigated in two of these cohorts. RESULTS: Earlier AAO was associated with a higher probability of psychotic symptoms, suicidality, lower educational attainment, not living together and fewer episodes. Depressive onset correlated with suicidality and manic onset correlated with delusions and manic episodes. Systematic differences in AAO between cohorts and continents of origin were observed. This was also reflected in single-nucleotide variant-based heritability estimates, with higher heritabilities for stricter onset definitions. Increased PGS for autism spectrum disorder (ß = -0.34 years, s.e. = 0.08), major depression (ß = -0.34 years, s.e. = 0.08), schizophrenia (ß = -0.39 years, s.e. = 0.08), and educational attainment (ß = -0.31 years, s.e. = 0.08) were associated with an earlier AAO. The AAO GWAS identified one significant locus, but this finding did not replicate. Neither GWAS nor PGS analyses yielded significant associations with PAO. CONCLUSIONS: AAO and PAO are associated with indicators of bipolar disorder severity. Individuals with an earlier onset show an increased polygenic liability for a broad spectrum of psychiatric traits. Systematic differences in AAO across cohorts, continents and phenotype definitions introduce significant heterogeneity, affecting analyses.

Identification of ADHD risk genes in extended pedigrees by combining linkage analysis and whole-exome sequencing.

Attention-deficit/hyperactivity disorder (ADHD) is a common neurodevelopmental disorder with a complex genetic background, hampering identification of underlying genetic risk factors. We hypothesized that combining linkage analysis and whole-exome sequencing (WES) in multi-generation pedigrees with multiple affected individuals can point toward novel ADHD genes. Three families with multiple ADHD-affected members (Ntotal = 70) and apparent dominant inheritance pattern were included in this study. Genotyping was performed in 37 family members, and WES was additionally carried out in 10 of those. Linkage analysis was performed using multi-point analysis in Superlink Online SNP 1.1. From prioritized linkage regions with a LOD score ≥ 2, a total of 24 genes harboring rare variants were selected. Those genes were taken forward and were jointly analyzed in gene-set analyses of exome-chip data using the MAGMA software in an independent sample of patients with persistent ADHD and healthy controls (N = 9365). The gene-set including all 24 genes together, and particularly the gene-set from one of the three families (12 genes), were significantly associated with persistent ADHD in this sample. Among the latter, gene-wide analysis for the AAED1 gene reached significance. A rare variant (rs151326868) within AAED1 segregated with ADHD in one of the families. The analytic strategy followed here is an effective approach for identifying novel ADHD risk genes. Additionally, this study suggests that both rare and more frequent variants in multiple genes act together in contributing to ADHD risk, even in individual multi-case families.

Correction: Dose response of the 16p11.2 distal copy number variant on intracranial volume and basal ganglia.

Prior to and following the publication of this article the authors noted that the complete list of authors was not included in the main article and was only present in Supplementary Table 1. The author list in the original article has now been updated to include all authors, and Supplementary Table 1 has been removed. All other supplementary files have now been updated accordingly. Furthermore, in Table 1 of this Article, the replication cohort for the row Close relative in data set, n (%) was incorrect. All values have now been corrected to 0(0%). The publishers would like to apologise for this error and the inconvenience it may have caused.

Dose response of the 16p11.2 distal copy number variant on intracranial volume and basal ganglia.

Carriers of large recurrent copy number variants (CNVs) have a higher risk of developing neurodevelopmental disorders. The 16p11.2 distal CNV predisposes carriers to e.g., autism spectrum disorder and schizophrenia. We compared subcortical brain volumes of 12 16p11.2 distal deletion and 12 duplication carriers to 6882 non-carriers from the large-scale brain Magnetic Resonance Imaging collaboration, ENIGMA-CNV. After stringent CNV calling procedures, and standardized FreeSurfer image analysis, we found negative dose-response associations with copy number on intracranial volume and on regional caudate, pallidum and putamen volumes (ß = -0.71 to -1.37; P < 0.0005). In an independent sample, consistent results were obtained, with significant effects in the pallidum (ß = -0.95, P = 0.0042). The two data sets combined showed significant negative dose-response for the accumbens, caudate, pallidum, putamen and ICV (P = 0.0032, 8.9 × 10-6, 1.7 × 10-9, 3.5 × 10-12 and 1.0 × 10-4, respectively). Full scale IQ was lower in both deletion and duplication carriers compared to non-carriers. This is the first brain MRI study of the impact of the 16p11.2 distal CNV, and we demonstrate a specific effect on subcortical brain structures, suggesting a neuropathological pattern underlying the neurodevelopmental syndromes.

Pathway-Specific Genetic Risk for Alzheimer's Disease Differentiates Regional Patterns of Cortical Atrophy in Older Adults.

Brain aging is highly variable and represents a challenge to delimit aging from disease processes. Moreover, genetic factors may influence both aging and disease. Here we focused on this issue and investigated effects of multiple genetic loci previously identified to be associated with late-onset Alzheimer's disease (AD) on brain structure of older adults from a population sample. We calculated a genetic risk score (GRS) using genome-wide significant single-nucleotide polymorphisms from genome-wide association studies of AD and tested its effect on cortical thickness (CT). We observed a common pattern of cortical thinning (right inferior frontal, left posterior temporal, medial occipital cortex). To identify CT changes by specific biological processes, we subdivided the GRS effect according to AD-associated pathways and performed follow-up analyses. The common pattern from the main analysis was further differentiated by pathway-specific effects yielding a more bilateral pattern. Further findings were located in the superior parietal and mid/anterior cingulate regions representing 2 unique pathway-specific patterns. All patterns, except the superior parietal pattern, were influenced by apolipoprotein E. Our step-wise approach revealed atrophy patterns that partially resembled imaging findings in early stages of AD. Our study provides evidence that genetic burden for AD contributes to structural brain variability in normal aging.

Shared genetic risk between eating disorder- and substance-use-related phenotypes: Evidence from genome-wide association studies.

Eating disorders and substance use disorders frequently co-occur. Twin studies reveal shared genetic variance between liabilities to eating disorders and substance use, with the strongest associations between symptoms of bulimia nervosa and problem alcohol use (genetic correlation [rg ], twin-based = 0.23-0.53). We estimated the genetic correlation between eating disorder and substance use and disorder phenotypes using data from genome-wide association studies (GWAS). Four eating disorder phenotypes (anorexia nervosa [AN], AN with binge eating, AN without binge eating, and a bulimia nervosa factor score), and eight substance-use-related phenotypes (drinks per week, alcohol use disorder [AUD], smoking initiation, current smoking, cigarettes per day, nicotine dependence, cannabis initiation, and cannabis use disorder) from eight studies were included. Significant genetic correlations were adjusted for variants associated with major depressive disorder and schizophrenia. Total study sample sizes per phenotype ranged from ~2400 to ~537 000 individuals. We used linkage disequilibrium score regression to calculate single nucleotide polymorphism-based genetic correlations between eating disorder- and substance-use-related phenotypes. Significant positive genetic associations emerged between AUD and AN (rg = 0.18; false discovery rate q = 0.0006), cannabis initiation and AN (rg = 0.23; q < 0.0001), and cannabis initiation and AN with binge eating (rg = 0.27; q = 0.0016). Conversely, significant negative genetic correlations were observed between three nondiagnostic smoking phenotypes (smoking initiation, current smoking, and cigarettes per day) and AN without binge eating (rgs = -0.19 to -0.23; qs < 0.04). The genetic correlation between AUD and AN was no longer significant after co-varying for major depressive disorder loci. The patterns of association between eating disorder- and substance-use-related phenotypes highlights the potentially complex and substance-specific relationships among these behaviors.

Mast cells increase vascular permeability by heparin-initiated bradykinin formation in vivo.

Activated mast cells trigger edema in allergic and inflammatory disease. We report a paracrine mechanism by which mast cell-released heparin increases vascular permeability in vivo. Heparin activated the protease factor XII, which initiates bradykinin formation in plasma. Targeting factor XII or kinin B2 receptors abolished heparin-triggered leukocyte-endothelium adhesion and interfered with a mast cell-driven drop in blood pressure in rodents. Intravital laser scanning microscopy and tracer measurements showed heparin-driven fluid extravasation in mouse skin microvessels. Ablation of factor XII or kinin B2 receptors abolished heparin-induced skin edema and protected mice from allergen-activated mast cell-driven leakage. In contrast, heparin and activated mast cells induced excessive edema in mice deficient in the major inhibitor of factor XII, C1 esterase inhibitor. Allergen exposure triggered edema attacks in hereditary angioedema patients, lacking C1 esterase inhibitor. The data indicate that heparin-initiated bradykinin formation plays a fundamental role in mast cell-mediated diseases.

Investigating polygenic burden in age at disease onset in bipolar disorder: Findings from an international multicentric study.

OBJECTIVES: Bipolar disorder (BD) with early disease onset is associated with an unfavorable clinical outcome and constitutes a clinically and biologically homogenous subgroup within the heterogeneous BD spectrum. Previous studies have found an accumulation of early age at onset (AAO) in BD families and have therefore hypothesized that there is a larger genetic contribution to the early-onset cases than to late onset BD. To investigate the genetic background of this subphenotype, we evaluated whether an increased polygenic burden of BD- and schizophrenia (SCZ)-associated risk variants is associated with an earlier AAO in BD patients. METHODS: A total of 1995 BD type 1 patients from the Consortium of Lithium Genetics (ConLiGen), PsyCourse and Bonn-Mannheim samples were genotyped and their BD and SCZ polygenic risk scores (PRSs) were calculated using the summary statistics of the Psychiatric Genomics Consortium as a training data set. AAO was either separated into onset groups of clinical interest (childhood and adolescence [≤18 years] vs adulthood [>18 years]) or considered as a continuous measure. The associations between BD- and SCZ-PRSs and AAO were evaluated with regression models. RESULTS: BD- and SCZ-PRSs were not significantly associated with age at disease onset. Results remained the same when analyses were stratified by site of recruitment. CONCLUSIONS: The current study is the largest conducted so far to investigate the association between the cumulative BD and SCZ polygenic risk and AAO in BD patients. The reported negative results suggest that such a polygenic influence, if there is any, is not large, and highlight the importance of conducting further, larger scale studies to obtain more information on the genetic architecture of this clinically relevant phenotype.

Leveraging Genomic Annotations and Pleiotropic Enrichment for Improved Replication Rates in Schizophrenia GWAS.

Most of the genetic architecture of schizophrenia (SCZ) has not yet been identified. Here, we apply a novel statistical algorithm called Covariate-Modulated Mixture Modeling (CM3), which incorporates auxiliary information (heterozygosity, total linkage disequilibrium, genomic annotations, pleiotropy) for each single nucleotide polymorphism (SNP) to enable more accurate estimation of replication probabilities, conditional on the observed test statistic ("z-score") of the SNP. We use a multiple logistic regression on z-scores to combine information from auxiliary information to derive a "relative enrichment score" for each SNP. For each stratum of these relative enrichment scores, we obtain nonparametric estimates of posterior expected test statistics and replication probabilities as a function of discovery z-scores, using a resampling-based approach that repeatedly and randomly partitions meta-analysis sub-studies into training and replication samples. We fit a scale mixture of two Gaussians model to each stratum, obtaining parameter estimates that minimize the sum of squared differences of the scale-mixture model with the stratified nonparametric estimates. We apply this approach to the recent genome-wide association study (GWAS) of SCZ (n = 82,315), obtaining a good fit between the model-based and observed effect sizes and replication probabilities. We observed that SNPs with low enrichment scores replicate with a lower probability than SNPs with high enrichment scores even when both they are genome-wide significant (p < 5x10-8). There were 693 and 219 independent loci with model-based replication rates ≥80% and ≥90%, respectively. Compared to analyses not incorporating relative enrichment scores, CM3 increased out-of-sample yield for SNPs that replicate at a given rate. This demonstrates that replication probabilities can be more accurately estimated using prior enrichment information with CM3.

No Reliable Association between Runs of Homozygosity and Schizophrenia in a Well-Powered Replication Study.

It is well known that inbreeding increases the risk of recessive monogenic diseases, but it is less certain whether it contributes to the etiology of complex diseases such as schizophrenia. One way to estimate the effects of inbreeding is to examine the association between disease diagnosis and genome-wide autozygosity estimated using runs of homozygosity (ROH) in genome-wide single nucleotide polymorphism arrays. Using data for schizophrenia from the Psychiatric Genomics Consortium (n = 21,868), Keller et al. (2012) estimated that the odds of developing schizophrenia increased by approximately 17% for every additional percent of the genome that is autozygous (ß = 16.1, CI(ß) = [6.93, 25.7], Z = 3.44, p = 0.0006). Here we describe replication results from 22 independent schizophrenia case-control datasets from the Psychiatric Genomics Consortium (n = 39,830). Using the same ROH calling thresholds and procedures as Keller et al. (2012), we were unable to replicate the significant association between ROH burden and schizophrenia in the independent PGC phase II data, although the effect was in the predicted direction, and the combined (original + replication) dataset yielded an attenuated but significant relationship between Froh and schizophrenia (ß = 4.86,CI(ß) = [0.90,8.83],Z = 2.40,p = 0.02). Since Keller et al. (2012), several studies reported inconsistent association of ROH burden with complex traits, particularly in case-control data. These conflicting results might suggest that the effects of autozygosity are confounded by various factors, such as socioeconomic status, education, urbanicity, and religiosity, which may be associated with both real inbreeding and the outcome measures of interest.

Genome-wide association study of 40,000 individuals identifies two novel loci associated with bipolar disorder.

Bipolar disorder (BD) is a genetically complex mental illness characterized by severe oscillations of mood and behaviour. Genome-wide association studies (GWAS) have identified several risk loci that together account for a small portion of the heritability. To identify additional risk loci, we performed a two-stage meta-analysis of >9 million genetic variants in 9,784 bipolar disorder patients and 30,471 controls, the largest GWAS of BD to date. In this study, to increase power we used ∼2,000 lithium-treated cases with a long-term diagnosis of BD from the Consortium on Lithium Genetics, excess controls, and analytic methods optimized for markers on the X-chromosome. In addition to four known loci, results revealed genome-wide significant associations at two novel loci: an intergenic region on 9p21.3 (rs12553324, P = 5.87 × 10 - 9; odds ratio (OR) = 1.12) and markers within ERBB2 (rs2517959, P = 4.53 × 10 - 9; OR = 1.13). No significant X-chromosome associations were detected and X-linked markers explained very little BD heritability. The results add to a growing list of common autosomal variants involved in BD and illustrate the power of comparing well-characterized cases to an excess of controls in GWAS.