Focal adhesion is associated with lithium response in bipolar disorder: evidence from a network-based multi-omics analysis

Lithium (Li) is one of the most effective drugs for treating bipolar disorder (BD), however, there is presently no way to predict response to guide treatment. The aim of this study is to identify functional genes and pathways that distinguish BD Li responders (LR) from BD Li non-responders (NR). An initial Pharmacogenomics of Bipolar Disorder study (PGBD) GWAS of lithium response did not provide any significant results. As a result, we then employed network-based integrative analysis of transcriptomic and genomic data. In transcriptomic study of iPSC-derived neurons, 41 significantly differentially expressed (DE) genes were identified in LR vs NR regardless of lithium exposure. In the PGBD, post-GWAS gene prioritization using the GWA-boosting (GWAB) approach identified 1119 candidate genes. Following DE-derived network propagation, there was a highly significant overlap of genes between the top 500- and top 2000-proximal gene networks and the GWAB gene list (Phypergeometric = 1.28E-09 and 4.10E-18, respectively). Functional enrichment analyses of the top 500 proximal network genes identified focal adhesion and the extracellular matrix (ECM) as the most significant functions. Our findings suggest that the difference between LR and NR was a much greater effect than that of lithium. The direct impact of dysregulation of focal adhesion on axon guidance and neuronal circuits could underpin mechanisms of response to lithium, as well as underlying BD. It also highlights the power of integrative multi-omics analysis of transcriptomic and genomic profiling to gain molecular insights into lithium response in BD.

A pilot pharmacogenetic study of calcium channel blocker treatment of bipolar mania

Common genetic variants located in calcium channel genes are important markers of genetic susceptibility for bipolar disorder (BD). Previous clinical trials with Calcium Channel Blocker (CCB) medication improved mood stability for some BD patients. We hypothesize that manic patients who carried calcium channel risk variants would differentially benefit from treatment with CCBs. In this pilot study, 50 BD patients (Chinese: 39; US: 11) who were hospitalized for manic episodes were given add-on CCB treatment. We determined genotypes for each patient. There was a significant decrease in the Young Mania Rating Scale (YMRS) after add-on medication treatment. Of note, two intronic variants of the Calcium Voltage-Gated Channel Subunit Alpha1 B (CACNA1B) were associated with treatment outcomes for manic patients: rs2739258 and rs2739260. BD rs2739258/rs2739260 AG-allele carriers had a better treatment response with add-on CCB than those carrying the AA or GG genotypes by survival analysis. Although these findings did not pass multiple testing correction, this study suggests that single-nucleotide polymorphisms (SNPs) residing in calcium channel genes could be predictors for response to add-on CCB treatment of bipolar mania patients, and that calcium channel genes may be involved in treatment responses for BD.

Characterization of childhood trauma, hippocampal mediation and Cannabis use in a large dataset of psychosis and non-psychosis individuals

BACKGROUND

Cannabis use (CA) and childhood trauma (CT) independently increase the risk of earlier psychosis onset; but their interaction in relation to psychosis risk and association with endocannabinoid-receptor rich brain regions, i.e. the hippocampus (HP), remains unclear. The objective was to determine whether lower age of psychosis onset (AgePsyOnset) is associated with CA and CT through mediation by the HP volumes, and genetic risk, as measured by schizophrenia polygene scores (SZ-PGRS).

METHODS

Cross-sectional, case-control, multicenter sample from 5 metropolitan US regions. Participants (n = 1185) included 397 controls not affected by psychosis (HC); 209 participants with bipolar disorder type-1; 279 with schizoaffective disorder; and 300 with schizophrenia (DSM IV-TR). CT was assessed using the Childhood Trauma Questionnaire (CTQ); CA was assessed by self-reports and trained clinical interviewers. Assessment included neuroimaging, symptomatology, cognition and calculation of the SZ polygenic risk score (SZ-PGRS).

RESULTS

In survival analysis, CT and CA exposure interact to be associated with lower AgePsyOnset. At high CT or CA, CT or CA are individually sufficient to affect AgePsyOnset. CT relation with AgePsyOnset is mediated in part by the HP in CA users before AgePsyOnset. CA before AgePsyOnset is associated with higher SZ-PGRS and correlated with younger age at CA usage.

DISCUSSION

CA and CT interact to increase risk when moderate; while severe CT and/or CA abuse/dependence are each sufficient to affect AgePsyOnset, indicating a ceiling effect. Probands with/out CA before AgePsyOnset differ on biological variables, suggesting divergent pathways to psychosis.

FUNDING

MH077945; MH096942; MH096913; MH077862; MH103368; MH096900; MH122759.

Single-neuron whole genome sequencing identifies increased somatic mutation burden in Alzheimer's disease related genes

Accumulation of somatic mutations in human neurons is associated with aging and neurodegeneration. To shed light on the somatic mutational burden in Alzheimer's disease (AD) neurons and get more insight into the role of somatic mutations in AD pathogenesis, we performed single-neuron whole genome sequencing to detect genome-wide somatic mutations (single nucleotide variants (SNVs) and Indels) in 96 single prefrontal cortex neurons from 8 AD patients and 8 elderly controls. We found that the mutational burden is ∼3000 somatic mutations per neuron genome in elderly subjects. AD patients have increased somatic mutation burden in AD-related annotation categories, including AD risk genes and differentially expressed genes in AD neurons. Mutational signature analysis showed somatic SNVs (sSNVs) primarily caused by aging and oxidative DNA damage processes but no significant difference was detected between AD and controls. Additionally, functional somatic mutations identified in AD patients showed significant enrichment in several AD-related pathways, including AD pathway, Notch-signaling pathway and Calcium-signaling pathway. These findings provide genetic insights into how somatic mutations may alter the function of single neurons and exert their potential roles in the pathogenesis of AD.

Correction of depression-associated circadian rhythm abnormalities is associated with lithium response in bipolar disorder

BACKGROUND

Bipolar disorder (BD) is characterized by episodes of depression and mania and disrupted circadian rhythms. Lithium is an effective therapy for BD, but only 30%-40% of patients are fully responsive. Preclinical models show that lithium alters circadian rhythms. However, it is unknown if the circadian rhythm effects of lithium are essential to its therapeutic properties.

METHODS

In secondary analyses of a multi-center, prospective, trial of lithium for BD, we examined the relationship between circadian rhythms and therapeutic response to lithium. Using standardized instruments, we measured morningness, diurnal changes in mood, sleep, and energy (circadian rhythm disturbances) in a cross-sectional study of 386 BD subjects with varying lithium exposure histories. Next, we tracked symptoms of depression and mania prospectively over 12 weeks in a subset of 88 BD patients initiating treatment with lithium. Total, circadian, and affective mood symptoms were scored separately and analyzed.

RESULTS

Subjects with no prior lithium exposure had the most circadian disruption, while patients stable on lithium monotherapy had the least. Patients who were stable on lithium with another drug or unstable on lithium showed intermediate levels of disruption. Treatment with lithium for 12 weeks yielded significant reductions in total and affective depression symptoms. Lithium responders (Li-Rs) showed improvement in circadian symptoms of depression, but non-responders did not. There was no difference between Li-Rs and nonresponders in affective, circadian, or total symptoms of mania.

CONCLUSIONS

Exposure to lithium is associated with reduced circadian disruption. Lithium response at 12 weeks was selectively associated with the reduction of circadian depressive symptoms. We conclude that stabilization of circadian rhythms may be an important feature of lithium's therapeutic effects.

CLINICAL TRIALS REGISTRY

NCT0127253.

Inflammation subtypes in psychosis and their relationships with genetic risk for psychiatric and cardiometabolic disorders

Cardiometabolic disorders have known inflammatory implications, and peripheral measures of inflammation and cardiometabolic disorders are common in persons with psychotic disorders. Inflammatory signatures are also related to neurobiological and behavioral changes in psychosis. Relationships between systemic inflammation and cardiometabolic genetic risk in persons with psychosis have not been examined. Thirteen peripheral inflammatory markers and genome-wide genotyping were assessed in 122 participants (n ​= ​86 psychosis, n ​= ​36 healthy controls) of European ancestry. Cluster analyses of inflammatory markers classified higher and lower inflammation subgroups. Single-trait genetic risk scores (GRS) were constructed for each participant using previously reported GWAS summary statistics for the following traits: schizophrenia, bipolar disorder, major depressive disorder, coronary artery disease, type-2 diabetes, low-density lipoprotein, high-density lipoprotein, triglycerides, and waist-to-hip ratio. Genetic correlations across traits were quantified. Principal component (PC) analysis of the cardiometabolic GRSs generated six PC loadings used in regression models to examine associations with inflammation markers. Functional module discovery explored biological mechanisms of the inflammation association of cardiometabolic GRS genes. A subgroup of 38% persons with psychotic disorders was characterized with higher inflammation status. These higher inflammation individuals had lower BACS scores (p ​= ​0.038) compared to those with lower inflammation. The first PC of the cardiometabolic GRS matrix was related to higher inflammation status in persons with psychotic disorders (OR ​= ​2.037, p ​= ​0.001). Two of eight modules within the functional interaction network of cardiometabolic GRS genes were enriched for immune processes. Cardiometabolic genetic risk may predispose some individuals with psychosis to elevated inflammation which adversely impacts cognition associated with illness.

An opportunity for primary prevention research in psychotic disorders

An opportunity has opened for research into primary prevention of psychotic disorders, based on progress in endophenotypes, genetics, and genomics. Primary prevention requires reliable prediction of susceptibility before any symptoms are present. We studied a battery of measures where published data supports abnormalities of these measurements prior to appearance of initial psychosis symptoms. These neurobiological and behavioral measurements included cognition, eye movement tracking, Event Related Potentials, and polygenic risk scores. They generated an acceptably precise separation of healthy controls from outpatients with a psychotic disorder. METHODS: The Bipolar and Schizophrenia Network on Intermediate Phenotypes (B-SNIP) measured this battery in an ancestry-diverse series of consecutively recruited adult outpatients with a psychotic disorder and healthy controls. Participants include all genders, 16 to 50 years of age, 261 with psychotic disorders (Schizophrenia (SZ) 109, Bipolar with psychosis (BPP) 92, Schizoaffective disorder (SAD) 60), 110 healthy controls. Logistic Regression, and an extension of the Linear Mixed Model to include analysis of pairwise interactions between measures (Environmental kernel Relationship Matrices (ERM)) with multiple iterations, were performed to predict case-control status. Each regression analysis was validated with four-fold cross-validation. RESULTS AND CONCLUSIONS: Sensitivity, specificity, and Area Under the Curve of Receiver Operating Characteristic of 85%, 62%, and 86%, respectively, were obtained for both analytic methods. These prediction metrics demonstrate a promising diagnostic distinction based on premorbid risk variables. There were also statistically significant pairwise interactions between measures in the ERM model. The strong prediction metrics of both types of analytic model provide proof-of-principle for biologically-based laboratory tests as a first step toward primary prevention studies. Prospective studies of adolescents at elevated risk, vs. healthy adolescent controls, would be a next step toward development of primary prevention strategies.

Impact of polygenic risk for coronary artery disease and cardiovascular medication burden on cognitive impairment in psychotic disorders

BACKGROUND

Cognitive impairment is a core deficit across psychotic disorders, the causes and therapeutics of which remain unclear. Epidemiological observations have suggested associations between cognitive dysfunction in psychotic disorders and cardiovascular risk factors, but an underlying etiology has not been established.

METHODS

Neuropsychological performance using the Brief Assessment of Cognition in Schizophrenia (BACS) was assessed in 616 individuals of European ancestry (403 psychosis, 213 controls). Polygenic risk scores for coronary artery disease (PRSCAD) were quantified for each participant across 13 p-value thresholds (PT 0.5-5e-8). Cardiovascular and psychotropic medications were categorized for association analyses. Each PRSCAD was examined in relation to the BACS and the optimized PT was confirmed with five-fold cross-validation and independent validation. Functional enrichment analyses were used to identify biological mechanisms linked to PRSCAD-cognition associations. Multiple regression analyses examined PRSCAD under the optimal PT and medication burden in relation to the BACS composite and subtest scores.

RESULTS

Higher PRSCAD was associated with lower BACS composite scores (p = 0.001) in the psychosis group, primarily driven by the Verbal Memory subtest (p < 0.001). Genes linked to multiple nervous system related processes and pathways were significantly enriched in PRSCAD. After controlling for PRSCAD, a greater number of cardiovascular medications was also correlated with worse BACS performance in patients with psychotic disorders (p = 0.029).

CONCLUSIONS

Higher PRSCAD and taking more cardiovascular medications were both significantly associated with cognitive impairment in psychosis. These findings indicate that cardiovascular factors may increase the risk for cognitive dysfunction and related functional outcomes among individuals with psychotic disorders.

Schizophrenia risk loci from xMHC region were associated with antipsychotic response in chronic schizophrenic patients with persistent positive symptom

We examined whether common variants from the extended major histocompatibility complex (xMHC) region contribute to the response to antipsychotic drugs (APDs) in patients with schizophrenia with persistent psychosis. Subjects participated in a prospective longitudinal study of the effect of APDs on psychopathology were temporally split into discovery (n = 88) and replication (n = 42) cohorts. The primary endpoint was a change in Brief Psychiatric Rating Scale at 6-week or 6-month after treatment. rs204991 (β = 3.917, p = 3.72 × 10^-6), the strongest signal associated with response at 6-week was located near C4A/C4B after a linear regression adjusted for covariates. xMHC SNP imputation disclosed much stronger signals (rs9268469, β = 5.140, p = 1.57 × 10^-7) and other weaker signals (p < 1 × 10^-5) spanning the entire xMHC region. All the variants were previously identified schizophrenia risk loci. Conditional fine-mapping revealed three subgroups of SNPs which were the eQTLs (p < 1 × 10^-7) for C4A, HLA-C, and BTN3A2 in disease-relevant tissue. Epistasis between HLA-C and C4A was observed (p = 0.019). Minor allele (G) carriers of rs204991, eQTL for C4A, having decreased risk for schizophrenia and lower imputed expression of C4A, had a better response to APDs. Some imputed HLA alleles associated with a decreased risk for schizophrenia had a positive association with improvement in psychotic symptoms. An independent cohort validated the association of change in psychosis with C4A. We provide evidence that genetic risk factors for schizophrenia from the xMHC region are associated with response to APDs and those variants significantly alter the imputed expression of C4A, HLA-C, and BTN3A2. The minor alleles predicting higher C4A level are associated with diminished improvement in psychotic symptoms after APD treatment.

Clinical predictors of non-response to lithium treatment in the Pharmacogenomics of Bipolar Disorder (PGBD) study

BACKGROUND

Lithium is regarded as a first-line treatment for bipolar disorder (BD), but partial response and non-response commonly occurs. There exists a need to identify lithium non-responders prior to initiating treatment. The Pharmacogenomics of Bipolar Disorder (PGBD) Study was designed to identify predictors of lithium response.

METHODS

The PGBD Study was an eleven site prospective trial of lithium treatment in bipolar I disorder. Subjects were stabilized on lithium monotherapy over 4 months and gradually discontinued from all other psychotropic medications. After ensuring a sustained clinical remission (defined by a score of ≤3 on the CGI for 4 weeks) had been achieved, subjects were followed for up to 2 years to monitor clinical response. Cox proportional hazard models were used to examine the relationship between clinical measures and time until failure to remit or relapse.

RESULTS

A total of 345 individuals were enrolled into the study and included in the analysis. Of these, 101 subjects failed to remit or relapsed, 88 achieved remission and continued to study completion, and 156 were terminated from the study for other reasons. Significant clinical predictors of treatment failure (p < 0.05) included baseline anxiety symptoms, functional impairments, negative life events and lifetime clinical features such as a history of migraine, suicidal ideation/attempts, and mixed episodes, as well as a chronic course of illness.

CONCLUSIONS

In this PGBD Study of lithium response, several clinical features were found to be associated with failure to respond to lithium. Future validation is needed to confirm these clinical predictors of treatment failure and their use clinically to distinguish who will do well on lithium before starting pharmacotherapy.

Genome-wide association study accounting for anticholinergic burden to examine cognitive dysfunction in psychotic disorders

Identifying genetic contributors to cognitive impairments in psychosis-spectrum disorders can advance understanding of disease pathophysiology. Although CNS medications are known to affect cognitive performance, they are often not accounted for in genetic association studies. In this study, we performed a genome-wide association study (GWAS) of global cognitive performance, measured as composite z-scores from the Brief Assessment of Cognition in Schizophrenia (BACS), in persons with psychotic disorders and controls (N = 817; 682 cases and 135 controls) from the Bipolar-Schizophrenia Network on Intermediate Phenotypes (B-SNIP) study. Analyses accounting for anticholinergic exposures from both psychiatric and non-psychiatric medications revealed five significantly associated variants located at the chromosome 3p21.1 locus, with the top SNP rs1076425 in the inter-alpha-trypsin inhibitor heavy chain 1 (ITIH1) gene (P = 3.25×E-9). The inclusion of anticholinergic burden improved association models (P < 0.001) and the number of significant SNPs identified. The effect sizes and direction of effect of the top variants remained consistent when investigating findings within individuals receiving specific antipsychotic drugs and after accounting for antipsychotic dose. These associations were replicated in a separate study sample of untreated first-episode psychosis. The chromosome 3p21.1 locus was previously reported to have association with the risk for psychotic disorders and cognitive performance in healthy individuals. Our findings suggest that this region may be a psychosis risk locus that is associated with cognitive mechanisms. Our data highlight the general point that the inclusion of medication exposure information may improve the detection of gene-cognition associations in psychiatric genetic research.

Genome-wide association study of more than 40,000 bipolar disorder cases provides new insights into the underlying biology

Bipolar disorder is a heritable mental illness with complex etiology. We performed a genome-wide association study of 41,917 bipolar disorder cases and 371,549 controls of European ancestry, which identified 64 associated genomic loci. Bipolar disorder risk alleles were enriched in genes in synaptic signaling pathways and brain-expressed genes, particularly those with high specificity of expression in neurons of the prefrontal cortex and hippocampus. Significant signal enrichment was found in genes encoding targets of antipsychotics, calcium channel blockers, antiepileptics and anesthetics. Integrating expression quantitative trait locus data implicated 15 genes robustly linked to bipolar disorder via gene expression, encoding druggable targets such as HTR6, MCHR1, DCLK3 and FURIN. Analyses of bipolar disorder subtypes indicated high but imperfect genetic correlation between bipolar disorder type I and II and identified additional associated loci. Together, these results advance our understanding of the biological etiology of bipolar disorder, identify novel therapeutic leads and prioritize genes for functional follow-up studies.

GWAS significance thresholds for deep phenotyping studies can depend upon minor allele frequencies and sample size

An important issue affecting genome-wide association studies with deep phenotyping (multiple correlated phenotypes) is determining the suitable family-wise significance threshold. Straightforward family-wise correction (Bonferroni) of p < 0.05 for 4.3 million genotypes and 335 phenotypes would give a threshold of p < 3.46E-11. This would be too conservative because it assumes all tests are independent. The effective number of tests, both phenotypic and genotypic, must be adjusted for the correlations between them. Spectral decomposition of the phenotype matrix and LD-based correction of the number of tested SNPs are currently used to determine an effective number of tests. In this paper, we compare these calculated estimates with permutation-determined family-wise significance thresholds. Permutations are performed by shuffling individual IDs of the genotype vector for this dataset, to preserve correlation of phenotypes. Our results demonstrate that the permutation threshold is influenced by minor allele frequency (MAF) of the SNPs, and by the number of individuals tested. For the more common SNPs (MAF > 0.1), the permutation family-wise threshold was in close agreement with spectral decomposition methods. However, for less common SNPs (0.05 < MAF ≤ 0.1), the permutation threshold calculated over all SNPs was off by orders of magnitude. This applies to the number of individuals studied (here 777) but not to very much larger numbers. Based on these findings, we propose that the threshold to find a particular level of family-wise significance may need to be established using separate permutations of the actual data for several MAF bins.

Neuropsychiatric Symptoms Among Hispanics: Results of the Maracaibo Aging Study

BACKGROUND

Neuropsychiatric symptoms play an important role in diagnosing and clinical follow-up of cognitive impairment and dementia.

OBJECTIVE

We investigated the relationship between neuropsychiatric symptoms, cognitive impairment, and dementia in Hispanics.

METHODS

We included 529 participants (age ≥40 years) from the Maracaibo Aging Study with standardized neuropsychiatric assessments, including the Neuropsychiatric Inventory (NPI). Based on the Clinical Dementia Rating and the Mini-Mental State Examination scores, participants' cognitive status was categorized into normal cognition, mild/moderate, and severe cognitive impairment. Diagnosis of dementia was established in a consensus conference. Statistical analyses included multivariable logistic regression models and area under the curve (AUC).

RESULTS

The mean age of participants was 59.3 years, and 71.8%were women. The proportion of dementia was 6.8%. Disturbed sleep, anxiety, and depression were the most common neuropsychiatric symptoms in the study sample. In crude analyses, the proportions of hallucinations, aberrant motor behavior, agitation/aggression, apathy, delusions, irritability, eating disturbance, depression, and euphoria were differently distributed among cognitive status groups (p < 0.05). After accounting for confounders, aberrant motor behavior and agitation/aggression remained significantly associated with cognitive impairment and dementia (p < 0.05). The inclusion of the NPI domains significantly improved the AUC to discriminate severe cognitive impairment and dementia compared to a basic model that included sex, age, education, alcohol, obesity, serum glucose, total cholesterol, hypertension, and stroke.

CONCLUSION

Neuropsychiatric symptoms are associated with severe cognitive impairment and dementia. The addition of NPI items to the global cognitive assessment might help early detection of dementia in primary care settings.

Biotyping in psychosis: using multiple computational approaches with one data set

Focusing on biomarker identification and using biomarkers individually or in clusters to define biological subgroups in psychiatry requires a re-orientation from behavioral phenomenology to quantifying brain features, requiring big data approaches for data integration. Much still needs to be accomplished, not only to refine but also to build support for the application and customization of such an analytical phenotypic approach. In this review, we present some of what Bipolar-Schizophrenia Network for Intermediate Phenotypes (B-SNIP) has learned so far to guide future applications of multivariate phenotyping and their analyses to understanding psychosis. This paper describes several B-SNIP projects that use phenotype data and big data computations to generate novel outcomes and glimpse what phenotypes contribute to disease understanding and, with aspiration, to treatment. The source of the phenotypes varies from genetic data, structural neuroanatomic localization, immune markers, brain physiology, and cognition. We aim to see guiding principles emerge and areas of commonality revealed. And, we will need to demonstrate not only data stability but also the usefulness of biomarker information for subgroup identification enhancing target identification and treatment development.

The association between lithium use and neurocognitive performance in patients with bipolar disorder

Lithium remains the gold standard for the treatment of bipolar disorder (BD); however, its use has declined over the years mainly due to the side effects and the subjective experience of cognitive numbness reported by patients. In the present study, we aim to methodically test the effects of lithium on neurocognitive functioning in the largest single cohort (n = 262) of BD patients reported to date by harnessing the power of a multi-site, ongoing clinical trial of lithium monotherapy. At the cross-sectional level, multivariate analysis of covariance (MANCOVA) was conducted to examine potential group differences across neurocognitive tests [California Verbal Learning Test (CVLT trials 1-5,CVLT delayed recall), Wechsler Digit Symbol, Trail-making Test parts A and B (TMT-A; TMT-B), and a global cognition index]. At the longitudinal level, on a subset of patients (n = 88) who achieved mood stabilization with lithium monotherapy, we explored the effect of lithium treatment across time on neurocognitive functioning. There were no differences at baseline between BD patients that were taking lithium compared with those that were not. At follow-up a significant neurocognitive improvement in the global cognitive index score [F = 31.69; p < 0.001], CVLT trials 1-5 [F = 29.81; p < 0.001], CVLT delayed recall [F = 15.27; p < 0.001], and TMT-B [F = 6.64, p = 0.012] was detected. The cross-sectional and longitudinal (on a subset of 88 patients) investigations suggest that lithium may be beneficial to neurocognitive functioning in patients with BD and that at the very least it does not seem to significantly impair cognition when used therapeutically.

NRXN1 is associated with enlargement of the temporal horns of the lateral ventricles in psychosis

Schizophrenia, Schizoaffective, and Bipolar disorders share behavioral and phenomenological traits, intermediate phenotypes, and some associated genetic loci with pleiotropic effects. Volumetric abnormalities in brain structures are among the intermediate phenotypes consistently reported associated with these disorders. In order to examine the genetic underpinnings of these structural brain modifications, we performed genome-wide association analyses (GWAS) on 60 quantitative structural brain MRI phenotypes in a sample of 777 subjects (483 cases and 294 controls pooled together). Genotyping was performed with the Illumina PsychChip microarray, followed by imputation to the 1000 genomes multiethnic reference panel. Enlargement of the Temporal Horns of Lateral Ventricles (THLV) is associated with an intronic SNP of the gene NRXN1 (rs12467877, P = 6.76E-10), which accounts for 4.5% of the variance in size. Enlarged THLV is associated with psychosis in this sample, and with reduction of the hippocampus and enlargement of the choroid plexus and caudate. Eight other suggestively significant associations (P < 5.5E-8) were identified with THLV and 5 other brain structures. Although rare deletions of NRXN1 have been previously associated with psychosis, this is the first report of a common SNP variant of NRXN1 associated with enlargement of the THLV in psychosis.

Genome-wide association study identifies 30 loci associated with bipolar disorder

Bipolar disorder is a highly heritable psychiatric disorder. We performed a genome-wide association study (GWAS) including 20,352 cases and 31,358 controls of European descent, with follow-up analysis of 822 variants with P < 1 × 10-4 in an additional 9,412 cases and 137,760 controls. Eight of the 19 variants that were genome-wide significant (P < 5 × 10-8) in the discovery GWAS were not genome-wide significant in the combined analysis, consistent with small effect sizes and limited power but also with genetic heterogeneity. In the combined analysis, 30 loci were genome-wide significant, including 20 newly identified loci. The significant loci contain genes encoding ion channels, neurotransmitter transporters and synaptic components. Pathway analysis revealed nine significantly enriched gene sets, including regulation of insulin secretion and endocannabinoid signaling. Bipolar I disorder is strongly genetically correlated with schizophrenia, driven by psychosis, whereas bipolar II disorder is more strongly correlated with major depressive disorder. These findings address key clinical questions and provide potential biological mechanisms for bipolar disorder.

Chronotype and cellular circadian rhythms predict the clinical response to lithium maintenance treatment in patients with bipolar disorder

Bipolar disorder (BD) is a serious mood disorder associated with circadian rhythm abnormalities. Risk for BD is genetically encoded and overlaps with systems that maintain circadian rhythms. Lithium is an effective mood stabilizer treatment for BD, but only a minority of patients fully respond to monotherapy. Presently, we hypothesized that lithium-responsive BD patients (Li-R) would show characteristic differences in chronotype and cellular circadian rhythms compared to lithium non-responders (Li-NR). Selecting patients from a prospective, multi-center, clinical trial of lithium monotherapy, we examined morning vs. evening preference (chronotype) as a dimension of circadian rhythm function in 193 Li-R and Li-NR BD patients. From a subset of 59 patient donors, we measured circadian rhythms in skin fibroblasts longitudinally over 5 days using a bioluminescent reporter (Per2-luc). We then estimated circadian rhythm parameters (amplitude, period, phase) and the pharmacological effects of lithium on rhythms in cells from Li-R and Li-NR donors. Compared to Li-NRs, Li-Rs showed a difference in chronotype, with higher levels of morningness. Evening chronotype was associated with increased mood symptoms at baseline, including depression, mania, and insomnia. Cells from Li-Rs were more likely to exhibit a short circadian period, a linear relationship between period and phase, and period shortening effects of lithium. Common genetic variation in the IP₃ signaling pathway may account for some of the individual differences in the effects of lithium on cellular rhythms. We conclude that circadian rhythms may influence response to lithium in maintenance treatment of BD.

Efficient region-based test strategy uncovers genetic risk factors for functional outcome in bipolar disorder

Genome-wide association studies of case-control status have advanced the understanding of the genetic basis of psychiatric disorders. Further progress may be gained by increasing sample size but also by new analysis strategies that advance the exploitation of existing data, especially for clinically important quantitative phenotypes. The functionally-informed efficient region-based test strategy (FIERS) introduced herein uses prior knowledge on biological function and dependence of genotypes within a powerful statistical framework with improved sensitivity and specificity for detecting consistent genetic effects across studies. As proof of concept, FIERS was used for the first genome-wide single nucleotide polymorphism (SNP)-based investigation on bipolar disorder (BD) that focuses on an important aspect of disease course, the functional outcome. FIERS identified a significantly associated locus on chromosome 15 (hg38: chr15:48965004 - 49464789 bp) with consistent effect strength between two independent studies (GAIN/TGen: European Americans, BOMA: Germans; n = 1592 BD patients in total). Protective and risk haplotypes were found on the most strongly associated SNPs. They contain a CTCF binding site (rs586758); CTCF sites are known to regulate sets of genes within a chromatin domain. The rs586758 - rs2086256 - rs1904317 haplotype is located in the promoter flanking region of the COPS2 gene, close to microRNA4716, and the EID1, SHC4, DTWD1 genes as plausible biological candidates. While implication with BD is novel, COPS2, EID1, and SHC4 are known to be relevant for neuronal differentiation and function and DTWD1 for psychopharmacological side effects. The test strategy FIERS that enabled this discovery is equally applicable for tag SNPs and sequence data.

Detecting significant genotype-phenotype association rules in bipolar disorder: market research meets complex genetics

Background

Disentangling the etiology of common, complex diseases is a major challenge in genetic research. For bipolar disorder (BD), several genome-wide association studies (GWAS) have been performed. Similar to other complex disorders, major breakthroughs in explaining the high heritability of BD through GWAS have remained elusive. To overcome this dilemma, genetic research into BD, has embraced a variety of strategies such as the formation of large consortia to increase sample size and sequencing approaches. Here we advocate a complementary approach making use of already existing GWAS data: a novel data mining procedure to identify yet undetected genotype–phenotype relationships. We adapted association rule mining, a data mining technique traditionally used in retail market research, to identify frequent and characteristic genotype patterns showing strong associations to phenotype clusters. We applied this strategy to three independent GWAS datasets from 2835 phenotypically characterized patients with BD. In a discovery step, 20,882 candidate association rules were extracted.

Results

Two of these rules—one associated with eating disorder and the other with anxiety—remained significant in an independent dataset after robust correction for multiple testing. Both showed considerable effect sizes (odds ratio ~ 3.4 and 3.0, respectively) and support previously reported molecular biological findings.

Conclusion

Our approach detected novel specific genotype–phenotype relationships in BD that were missed by standard analyses like GWAS. While we developed and applied our method within the context of BD gene discovery, it may facilitate identifying highly specific genotype–phenotype relationships in subsets of genome-wide data sets of other complex phenotype with similar epidemiological properties and challenges to gene discovery efforts.

Electronic supplementary material

The online version of this article (10.1186/s40345-018-0132-x) contains supplementary material, which is available to authorized users.

Genetic analysis of deep phenotyping projects in common disorders

Several studies of complex psychotic disorders with large numbers of neurobiological phenotypes are currently under way, in living patients and controls, and on assemblies of brain specimens. Genetic analyses of such data typically present challenges, because of the choice of underlying hypotheses on genetic architecture of the studied disorders and phenotypes, large numbers of phenotypes, the appropriate multiple testing corrections, limited numbers of subjects, imputations required on missing phenotypes and genotypes, and the cross-disciplinary nature of the phenotype measures. Advances in genotype and phenotype imputation, and in genome-wide association (GWAS) methods, are useful in dealing with these challenges. As compared with the more traditional single-trait analyses, deep phenotyping with simultaneous genome-wide analyses serves as a discovery tool for previously unsuspected relationships of phenotypic traits with each other, and with specific molecular involvements.

Genome-wide association studies of smooth pursuit and antisaccade eye movements in psychotic disorders: findings from the B-SNIP study

Eye movement deviations, particularly deficits of initial sensorimotor processing and sustained pursuit maintenance, and antisaccade inhibition errors, are established intermediate phenotypes for psychotic disorders. We here studied eye movement measures of 849 participants from the Bipolar-Schizophrenia Network on Intermediate Phenotypes (B-SNIP) study (schizophrenia N=230, schizoaffective disorder N=155, psychotic bipolar disorder N=206 and healthy controls N=258) as quantitative phenotypes in relation to genetic data, while controlling for genetically derived ancestry measures, age and sex. A mixed-modeling genome-wide association studies approach was used including ~4.4 million genotypes (PsychChip and 1000 Genomes imputation). Across participants, sensorimotor processing at pursuit initiation was significantly associated with a single nucleotide polymorphism in IPO8 (12p11.21, P=8 × 10^-11), whereas suggestive associations with sustained pursuit maintenance were identified with SNPs in SH3GL2 (9p22.2, P=3 × 10^-8). In participants of predominantly African ancestry, sensorimotor processing was also significantly associated with SNPs in PCDH12 (5q31.3, P=1.6 × 10^-10), and suggestive associations were observed with NRSN1 (6p22.3, P=5.4 × 10^-8) and LMO7 (13q22.2, P=7.3x10^-8), whereas antisaccade error rate was significantly associated with a non-coding region at chromosome 7 (P=6.5 × 10^-9). Exploratory pathway analyses revealed associations with nervous system development and function for 40 top genes with sensorimotor processing and pursuit maintenance (P=4.9 × 10^-2-9.8 × 10^-4). Our findings suggest novel patterns of genetic variation relevant for brain systems subserving eye movement control known to be impaired in psychotic disorders. They include genes involved in nuclear trafficking and gene silencing (IPO8), fast axonal guidance and synaptic specificity (PCDH12), transduction of nerve signals (NRSN1), retinal degeneration (LMO7), synaptic glutamate release (SH3GL2), and broader nervous system development and function.

Toxoplasma Modulates Signature Pathways of Human Epilepsy, Neurodegeneration & Cancer

One third of humans are infected lifelong with the brain-dwelling, protozoan parasite, Toxoplasma gondii. Approximately fifteen million of these have congenital toxoplasmosis. Although neurobehavioral disease is associated with seropositivity, causality is unproven. To better understand what this parasite does to human brains, we performed a comprehensive systems analysis of the infected brain: We identified susceptibility genes for congenital toxoplasmosis in our cohort of infected humans and found these genes are expressed in human brain. Transcriptomic and quantitative proteomic analyses of infected human, primary, neuronal stem and monocytic cells revealed effects on neurodevelopment and plasticity in neural, immune, and endocrine networks. These findings were supported by identification of protein and miRNA biomarkers in sera of ill children reflecting brain damage and T. gondii infection. These data were deconvoluted using three systems biology approaches: "Orbital-deconvolution" elucidated upstream, regulatory pathways interconnecting human susceptibility genes, biomarkers, proteomes, and transcriptomes. "Cluster-deconvolution" revealed visual protein-protein interaction clusters involved in processes affecting brain functions and circuitry, including lipid metabolism, leukocyte migration and olfaction. Finally, "disease-deconvolution" identified associations between the parasite-brain interactions and epilepsy, movement disorders, Alzheimer's disease, and cancer. This "reconstruction-deconvolution" logic provides templates of progenitor cells' potentiating effects, and components affecting human brain parasitism and diseases.

Genome-wide association study of borderline personality disorder reveals genetic overlap with bipolar disorder, major depression and schizophrenia

Borderline personality disorder (BOR) is determined by environmental and genetic factors, and characterized by affective instability and impulsivity, diagnostic symptoms also observed in manic phases of bipolar disorder (BIP). Up to 20% of BIP patients show comorbidity with BOR. This report describes the first case-control genome-wide association study (GWAS) of BOR, performed in one of the largest BOR patient samples worldwide. The focus of our analysis was (i) to detect genes and gene sets involved in BOR and (ii) to investigate the genetic overlap with BIP. As there is considerable genetic overlap between BIP, major depression (MDD) and schizophrenia (SCZ) and a high comorbidity of BOR and MDD, we also analyzed the genetic overlap of BOR with SCZ and MDD. GWAS, gene-based tests and gene-set analyses were performed in 998 BOR patients and 1545 controls. Linkage disequilibrium score regression was used to detect the genetic overlap between BOR and these disorders. Single marker analysis revealed no significant association after correction for multiple testing. Gene-based analysis yielded two significant genes: DPYD (P=4.42 × 10-7) and PKP4 (P=8.67 × 10-7); and gene-set analysis yielded a significant finding for exocytosis (GO:0006887, PFDR=0.019; FDR, false discovery rate). Prior studies have implicated DPYD, PKP4 and exocytosis in BIP and SCZ. The most notable finding of the present study was the genetic overlap of BOR with BIP (rg=0.28 [P=2.99 × 10-3]), SCZ (rg=0.34 [P=4.37 × 10-5]) and MDD (rg=0.57 [P=1.04 × 10-3]). We believe our study is the first to demonstrate that BOR overlaps with BIP, MDD and SCZ on the genetic level. Whether this is confined to transdiagnostic clinical symptoms should be examined in future studies.

Ethical and public policy challenges for pharmacogenomics

It is timely to consider the ethical and social questions raised by progress in pharmacogenomics, based on the current importance of pharmacogenomics for avoidance of predictable side effects of drugs, and for correct choice of medications in certain cancers. It has been proposed that the entire population be genotyped for drug-metabolizing enzyme polymorphisms, as a measure that would prevent many untoward and dangerous drug reactions. Pharmacologic treatment targeting based on genomics of disease can be expected to increase greatly in the coming years. Policy and ethical issues exist on consent for large-scale genomic pharmacogenomic data collection, public vs corporate ownership of genomic research results, testing efficacy and safety of drugs used for rare genomic indications, and accessibility of treatments based on costly research that is applicable to relatively few patients. In major psychiatric disorders and intellectual deficiency, rare and de novo deletion or duplication of chromosomal segments (copy number variation), in the aggregate, are common causes of increased risk. This implies that the policy problems of pharmacogenomics will be particularly important for the psychiatric disorders.

New ethical issues for genetic counseling in common mental disorders

OBJECTIVE

Recent genetic findings of high-impact genetic variants in bipolar disorder, schizophrenia, and autism spectrum disorder (ASD) must lead to profound changes in genetic and family counseling. The authors present risk calculations, discuss the ethical implications of these findings, and outline the changes now required in the risk counseling process.

METHOD

The authors use data from recent mega-analyses and reviews of common and rare risk variants in bipolar disorder, schizophrenia, and ASD to calculate risks of illness based on genetic marker tests. They then consider new ethical issues in mental disorders presented by these risks, including within-family conflicts over genetic testing; effects of genetic discoveries on stigma, abortion, preimplantation procedures, and population screening for susceptibility; and genetic tests as a factor in marital choice.

RESULTS

New structural mutations (de novo copy number variants [CNVs], which are chromosomal microdeletions and micro-duplications) are present in 4%27% of patients with bipolar disorder, schizophrenia, or ASD and can occur almost anywhere in the genome. For a person with a de novo CNV, the absolute risk of bipolar disorder, schizophrenia, or ASD is 14%, much higher than the population risk. Rare CNVs have also been identified that are generally not new mutations but constitute very high-effect risk factors, ranging up to 82%.

CONCLUSIONS

A substantial minority of patients with bipolar disorder, schizophrenia, and ASD have high-impact detectable genetic events. This greatly changes psychiatric genetic counseling for these patients and families. A psychotherapeutic approach may be needed as a routine part of risk counseling, particularly for resolution of ethical issues and for within-family stigma and conflicts over genetic test results.

After GWAS: searching for genetic risk for schizophrenia and bipolar disorder

Ten years ago it was widely expected that the genetic basis of common disease would be resolved by genome-wide association studies (GWAS), large-scale studies in which the entire genome is covered by genetic markers. However, the bulk of heritable variance remains unexplained. The authors consider several alternative research strategies. For instance, whereas it has been hypothesized that a common disease is associated primarily with common genetic variants, it is now plausible that multiple rare variants each have a potent effect on disease risk and that they could accumulate to become a substantial component of common disease risk. This idea has become more appealing since the discovery that copy number variants (CNVs) are a substantial source of human mutations and are associated with multiple common diseases. CNVs are structural genomic variants consisting of microinsertions, microdeletions, and transpositions in the human genome. It has been argued that numerous rare CNVs are plausible causes of a substantial proportion of common disease, and rare CNVs have been found to be potent risk factors in schizophrenia and autism. Another approach is to "parse the genome," i.e., reanalyze subsets of current GWAS data, since the noise inherent in genome-wide approaches may be hiding valid associations. Lastly, technological advances and declining costs may allow large-scale genome-wide sequencing that would comprehensively identify all genetic variants. Study groups even larger than the 10,000 subjects in current meta-analyses would be required, but the outcomes may lead to resolution of our current dilemma in common diseases: Where is the missing heritability?

Accuracy of CNV Detection from GWAS Data

Several computer programs are available for detecting copy number variants (CNVs) using genome-wide SNP arrays. We evaluated the performance of four CNV detection software suites—Birdsuite, Partek, HelixTree, and PennCNV-Affy—in the identification of both rare and common CNVs. Each program's performance was assessed in two ways. The first was its recovery rate, i.e., its ability to call 893 CNVs previously identified in eight HapMap samples by paired-end sequencing of whole-genome fosmid clones, and 51,440 CNVs identified by array Comparative Genome Hybridization (aCGH) followed by validation procedures, in 90 HapMap CEU samples. The second evaluation was program performance calling rare and common CNVs in the Bipolar Genome Study (BiGS) data set (1001 bipolar cases and 1033 controls, all of European ancestry) as measured by the Affymetrix SNP 6.0 array. Accuracy in calling rare CNVs was assessed by positive predictive value, based on the proportion of rare CNVs validated by quantitative real-time PCR (qPCR), while accuracy in calling common CNVs was assessed by false positive/false negative rates based on qPCR validation results from a subset of common CNVs. Birdsuite recovered the highest percentages of known HapMap CNVs containing >20 markers in two reference CNV datasets. The recovery rate increased with decreased CNV frequency. In the tested rare CNV data, Birdsuite and Partek had higher positive predictive values than the other software suites. In a test of three common CNVs in the BiGS dataset, Birdsuite's call was 98.8% consistent with qPCR quantification in one CNV region, but the other two regions showed an unacceptable degree of accuracy. We found relatively poor consistency between the two “gold standards,” the sequence data of Kidd et al., and aCGH data of Conrad et al. Algorithms for calling CNVs especially common ones need substantial improvement, and a “gold standard” for detection of CNVs remains to be established.