Lithium response in bipolar disorder is associated with focal adhesion and PI3K-Akt networks: a multi-omics replication study

Lithium is the gold standard treatment for bipolar disorder (BD). However, its mechanism of action is incompletely understood, and prediction of treatment outcomes is limited. In our previous multi-omics study of the Pharmacogenomics of Bipolar Disorder (PGBD) sample combining transcriptomic and genomic data, we found that focal adhesion, the extracellular matrix (ECM), and PI3K-Akt signaling networks were associated with response to lithium. In this study, we replicated the results of our previous study using network propagation methods in a genome-wide association study of an independent sample of 2039 patients from the International Consortium on Lithium Genetics (ConLiGen) study. We identified functional enrichment in focal adhesion and PI3K-Akt pathways, but we did not find an association with the ECM pathway. Our results suggest that deficits in the neuronal growth cone and PI3K-Akt signaling, but not in ECM proteins, may influence response to lithium in BD.

A resource of induced pluripotent stem cell (iPSC) lines including clinical, genomic, and cellular data from genetically isolated families with mood and psychotic disorders

Genome-wide (GWAS) and copy number variant (CNV) association studies have reproducibly identified numerous risk alleles associated with bipolar disorder (BD), major depressive disorder (MDD), and schizophrenia (SCZ), but biological characterization of these alleles lags gene discovery, owing to the inaccessibility of live human brain cells and inadequate animal models for human psychiatric conditions. Human-derived induced pluripotent stem cells (iPSCs) provide a renewable cellular reagent that can be differentiated into living, disease-relevant cells and 3D brain organoids carrying the full complement of genetic variants present in the donor germline. Experimental studies of iPSC-derived cells allow functional characterization of risk alleles, establishment of causal relationships between genes and neurobiology, and screening for novel therapeutics. Here we report the creation and availability of an iPSC resource comprising clinical, genomic, and cellular data obtained from genetically isolated families with BD and related conditions. Results from the first 324 study participants, 61 of whom have validated pluripotent clones, show enrichment of rare single nucleotide variants and CNVs overlapping many known risk genes and pathogenic CNVs. This growing iPSC resource is available to scientists pursuing functional genomic studies of BD and related conditions.

Exploring the genetics of lithium response in bipolar disorders

Background

Lithium (Li) remains the treatment of choice for bipolar disorders (BP). Its mood-stabilizing effects help reduce the long-term burden of mania, depression and suicide risk in patients with BP. It also has been shown to have beneficial effects on disease-associated conditions, including sleep and cardiovascular disorders. However, the individual responses to Li treatment vary within and between diagnostic subtypes of BP (e.g. BP-I and BP-II) according to the clinical presentation. Moreover, long-term Li treatment has been linked to adverse side-effects that are a cause of concern and non-adherence, including the risk of developing chronic medical conditions such as thyroid and renal disease. In recent years, studies by the Consortium on Lithium Genetics (ConLiGen) have uncovered a number of genetic factors that contribute to the variability in Li treatment response in patients with BP. Here, we leveraged the ConLiGen cohort (N=2,064) to investigate the genetic basis of Li effects in BP. For this, we studied how Li response and linked genes associate with the psychiatric symptoms and polygenic load for medical comorbidities, placing particular emphasis on identifying differences between BP-I and BP-II.

Results

We found that clinical response to Li treatment, measured with the Alda scale, was associated with a diminished burden of mania, depression, substance and alcohol abuse, psychosis and suicidal ideation in patients with BP-I and, in patients with BP-II, of depression only. Our genetic analyses showed that a stronger clinical response to Li was modestly related to lower polygenic load for diabetes and hypertension in BP-I but not BP-II. Moreover, our results suggested that a number of genes that have been previously linked to Li response variability in BP differentially relate to the psychiatric symptomatology, particularly to the numbers of manic and depressive episodes, and to the polygenic load for comorbid conditions, including diabetes, hypertension and hypothyroidism.

Conclusions

Taken together, our findings suggest that the effects of Li on symptomatology and comorbidity in BP are partially modulated by common genetic factors, with differential effects between BP-I and BP-II.

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.

Lithium Response in Bipolar Disorder is Associated with Focal Adhesion and PI3K-Akt Networks: A Multi-omics Replication Study

Lithium is the gold standard treatment for bipolar disorder (BD). However, its mechanism of action is incompletely understood, and prediction of treatment outcomes is limited. In our previous multi-omics study of the Pharmacogenomics of Bipolar Disorder (PGBD) sample combining transcriptomic and genomic data, we found that focal adhesion, the extracellular matrix (ECM), and PI3K-Akt signaling networks were associated with response to lithium. In this study, we replicated the results of our previous study using network propagation methods in a genome-wide association study of an independent sample of 2,039 patients from the International Consortium on Lithium Genetics (ConLiGen) study. We identified functional enrichment in focal adhesion and PI3K-Akt pathways, but we did not find an association with the ECM pathway. Our results suggest that deficits in the neuronal growth cone and PI3K-Akt signaling, but not in ECM proteins, may influence response to lithium in BD.

Exome-wide association study of treatment-resistant depression suggests novel treatment targets

Treatment-resistant depression (TRD) is a severe form of major depressive disorder (MDD) with substantial public health impact and poor treatment outcome. Treatment outcome in MDD is significantly heritable, but genome-wide association studies have failed to identify replicable common marker alleles, suggesting a potential role for uncommon variants. Here we investigated the hypothesis that uncommon, putatively functional genetic variants are associated with TRD. Whole-exome sequencing data was obtained from 182 TRD cases and 2021 psychiatrically healthy controls. After quality control, the remaining 149 TRD cases and 1976 controls were analyzed with tests designed to detect excess burdens of uncommon variants. At the gene level, 5 genes, ZNF248, PRKRA, PYHIN1, SLC7A8, and STK19 each carried exome-wide significant excess burdens of variants in TRD cases (q < 0.05). Analysis of 41 pre-selected gene sets suggested an excess of uncommon, functional variants among genes involved in lithium response. Among the genes identified in previous TRD studies, ZDHHC3 was also significant in this sample after multiple test correction. ZNF248 and STK19 are involved in transcriptional regulation, PHYIN1 and PRKRA are involved in immune response, SLC7A8 is associated with thyroid hormone transporter activity, and ZDHHC3 regulates synaptic clustering of GABA and glutamate receptors. These results implicate uncommon, functional alleles in TRD and suggest promising novel targets for future research.

Immunogenetics of lithium response and psychiatric phenotypes in patients with bipolar disorder

The link between bipolar disorder (BP) and immune dysfunction remains controversial. While epidemiological studies have long suggested an association, recent research has found only limited evidence of such a relationship. To clarify this, we investigated the contributions of immune-relevant genetic factors to the response to lithium (Li) treatment and the clinical presentation of BP. First, we assessed the association of a large collection of immune-related genes (4,925) with Li response, defined by the Retrospective Assessment of the Lithium Response Phenotype Scale (Alda scale), and clinical characteristics in patients with BP from the International Consortium on Lithium Genetics (ConLi+Gen, N = 2,374). Second, we calculated here previously published polygenic scores (PGSs) for immune-related traits and evaluated their associations with Li response and clinical features. We found several genes associated with Li response at p < 1×10- 4 values, including HAS3, CNTNAP5 and NFIB. Network and functional enrichment analyses uncovered an overrepresentation of pathways involved in cell adhesion and intercellular communication, which appear to converge on the well-known Li-induced inhibition of GSK-3β. We also found various genes associated with BP's age-at-onset, number of mood episodes, and presence of psychosis, substance abuse and/or suicidal ideation at the exploratory threshold. These included RTN4, XKR4, NRXN1, NRG1/3 and GRK5. Additionally, PGS analyses suggested serum FAS, ECP, TRANCE and cytokine ligands, amongst others, might represent potential circulating biomarkers of Li response and clinical presentation. Taken together, our results support the notion of a relatively weak association between immunity and clinically relevant features of BP at the genetic level.

Aneuploidy effects on human gene expression across three cell types

Aneuploidy syndromes impact multiple organ systems but understanding of tissue-specific aneuploidy effects remains limited-especially for the comparison between peripheral tissues and relatively inaccessible tissues like brain. Here, we address this gap in knowledge by studying the transcriptomic effects of chromosome X, Y, and 21 aneuploidies in lymphoblastoid cell lines, fibroblasts and iPSC-derived neuronal cells (LCLs, FCL, and iNs, respectively). We root our analyses in sex chromosome aneuploidies, which offer a uniquely wide karyotype range for dosage effect analysis. We first harness a large LCL RNA-seq dataset from 197 individuals with one of 6 sex chromosome dosages (SCDs: XX, XXX, XY, XXY, XYY, and XXYY) to i) validate theoretical models of SCD sensitivity and ii) define an expanded set of 41 genes that show obligate dosage sensitivity to SCD and are all in cis (i.e., reside on the X or Y chromosome). We then use multiple complementary analyses to show that cis effects of SCD in LCLs are preserved in both FCLs (n = 32) and iNs (n = 24), whereas trans effects (i.e., those on autosomal gene expression) are mostly not preserved. Analysis of additional datasets confirms that the greater cross-cell type reproducibility of cis vs. trans effects is also seen in trisomy 21 cell lines. These findings i) expand our understanding of X, Y, and 21 chromosome dosage effects on human gene expression and ii) suggest that LCLs may provide a good model system for understanding cis effects of aneuploidy in harder-to-access cell types.

Cellular Diversity in Human Subgenual Anterior Cingulate and Dorsolateral Prefrontal Cortex by Single-Nucleus RNA-Sequencing

Regional cellular heterogeneity is a fundamental feature of the human neocortex; however, details of this heterogeneity are still undefined. We used single-nucleus RNA-sequencing to examine cell-specific transcriptional features in the dorsolateral PFC (DLPFC) and the subgenual anterior cingulate cortex (sgACC), regions implicated in major psychiatric disorders. Droplet-based nuclei-capture and library preparation were performed on replicate samples from 8 male donors without history of psychiatric or neurologic disorder. Unsupervised clustering identified major neural cell classes. Subsequent iterative clustering of neurons further revealed 20 excitatory and 22 inhibitory subclasses. Inhibitory cells were consistently more abundant in the sgACC and excitatory neuron subclusters exhibited considerable variability across brain regions. Excitatory cell subclasses also exhibited greater within-class transcriptional differences between the two regions. We used these molecular definitions to determine which cell classes might be enriched in loci carrying a genetic signal in genome-wide association studies or for differentially expressed genes in mental illness. We found that the heritable signals of psychiatric disorders were enriched in neurons and that, while the gene expression changes detected in bulk-RNA-sequencing studies were dominated by glial cells, some alterations could be identified in specific classes of excitatory and inhibitory neurons. Intriguingly, only two excitatory cell classes exhibited concomitant region-specific enrichment for both genome-wide association study loci and transcriptional dysregulation. In sum, by detailing the molecular and cellular diversity of the DLPFC and sgACC, we were able to generate hypotheses on regional and cell-specific dysfunctions that may contribute to the development of mental illness.SIGNIFICANCE STATEMENT Dysfunction of the subgenual anterior cingulate cortex has been implicated in mood disorders, particularly major depressive disorder, and the dorsolateral PFC, a subsection of the PFC involved in executive functioning, has been implicated in schizophrenia. Understanding the cellular composition of these regions is critical to elucidating the neurobiology underlying psychiatric and neurologic disorders. We studied cell type diversity of the subgenual anterior cingulate cortex and dorsolateral PFC of humans with no neuropsychiatric illness using a clustering analysis of single-nuclei RNA-sequencing data. Defining the transcriptomic profile of cellular subpopulations in these cortical regions is a first step to demystifying the cellular and molecular pathways involved in psychiatric disorders.

Antipsychotic drug use complicates assessment of gene expression changes associated with schizophrenia

Recent postmortem transcriptomic studies of schizophrenia (SCZ) have shown hundreds of differentially expressed genes. However, the extent to which these gene expression changes reflect antipsychotic drug (APD) exposure remains uncertain. We compared differential gene expression in the prefrontal cortex of SCZ patients who tested positive for APDs at the time of death with SCZ patients who did not. APD exposure was associated with numerous changes in the brain transcriptome, especially among SCZ patients on atypical APDs. Brain transcriptome data from macaques chronically treated with APDs showed that APDs affect the expression of many functionally relevant genes, some of which show expression changes in the same directions as those observed in SCZ. Co-expression modules enriched for synaptic function showed convergent patterns between SCZ and some of the APD effects, while those associated with inflammation and glucose metabolism exhibited predominantly divergent patterns between SCZ and APD effects. In contrast, major cell-type shifts inferred in SCZ were primarily unaffected by APD use. These results show that APDs may confound SCZ-associated gene expression changes in postmortem brain tissue. Disentangling these effects will help identify causal genes and improve our neurobiological understanding of SCZ.

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 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=2,367) 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<����������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������.

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.

Association of Attention-Deficit/Hyperactivity Disorder and Depression Polygenic Scores with Lithium Response: A Consortium for Lithium Genetics Study

Response to lithium varies widely between individuals with bipolar disorder (BD). Polygenic risk scores (PRSs) can uncover pharmacogenomics effects and may help predict drug response. Patients (N = 2,510) with BD were assessed for long-term lithium response in the Consortium on Lithium Genetics using the Retrospective Criteria of Long-Term Treatment Response in Research Subjects with Bipolar Disorder score. PRSs for attention-deficit/hyperactivity disorder (ADHD), major depressive disorder (MDD), and schizophrenia (SCZ) were computed using lassosum and in a model including all three PRSs and other covariates, and the PRS of ADHD (β = -0.14; 95% confidence interval [CI]: -0.24 to -0.03; p value = 0.010) and MDD (β = -0.16; 95% CI: -0.27 to -0.04; p value = 0.005) predicted worse quantitative lithium response. A higher SCZ PRS was associated with higher rates of medication nonadherence (OR = 1.61; 95% CI: 1.34-1.93; p value = 2e-7). This study indicates that genetic risk for ADHD and depression may influence lithium treatment response. Interestingly, a higher SCZ PRS was associated with poor adherence, which can negatively impact treatment response. Incorporating genetic risk of ADHD, depression, and SCZ in combination with clinical risk may lead to better clinical care for patients with BD.

HLA-DRB1 and HLA-DQB1 genetic diversity modulates response to lithium in bipolar affective disorders

Bipolar affective disorder (BD) is a severe psychiatric illness, for which lithium (Li) is the gold standard for acute and maintenance therapies. The therapeutic response to Li in BD is heterogeneous and reliable biomarkers allowing patients stratification are still needed. A GWAS performed by the International Consortium on Lithium Genetics (ConLiGen) has recently identified genetic markers associated with treatment responses to Li in the human leukocyte antigens (HLA) region. To better understand the molecular mechanisms underlying this association, we have genetically imputed the classical alleles of the HLA region in the European patients of the ConLiGen cohort. We found our best signal for amino-acid variants belonging to the HLA-DRB1*11:01 classical allele, associated with a better response to Li (p < 1 × 10-3; FDR < 0.09 in the recessive model). Alanine or Leucine at position 74 of the HLA-DRB1 heavy chain was associated with a good response while Arginine or Glutamic acid with a poor response. As these variants have been implicated in common inflammatory/autoimmune processes, our findings strongly suggest that HLA-mediated low inflammatory background may contribute to the efficient response to Li in BD patients, while an inflammatory status overriding Li anti-inflammatory properties would favor a weak response.

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.

Prediction of lithium response using genomic data

Predicting lithium response prior to treatment could both expedite therapy and avoid exposure to side effects. Since lithium responsiveness may be heritable, its predictability based on genomic data is of interest. We thus evaluate the degree to which lithium response can be predicted with a machine learning (ML) approach using genomic data. Using the largest existing genomic dataset in the lithium response literature (n = 2210 across 14 international sites; 29% responders), we evaluated the degree to which lithium response could be predicted based on 47,465 genotyped single nucleotide polymorphisms using a supervised ML approach. Under appropriate cross-validation procedures, lithium response could be predicted to above-chance levels in two constituent sites (Halifax, Cohen’s kappa 0.15, 95% confidence interval, CI [0.07, 0.24]; and Würzburg, kappa 0.2 [0.1, 0.3]). Variants with shared importance in these models showed over-representation of postsynaptic membrane related genes. Lithium response was not predictable in the pooled dataset (kappa 0.02 [− 0.01, 0.04]), although non-trivial performance was achieved within a restricted dataset including only those patients followed prospectively (kappa 0.09 [0.04, 0.14]). Genomic classification of lithium response remains a promising but difficult task. Classification performance could potentially be improved by further harmonization of data collection procedures.

Polygenic risk for anxiety influences anxiety comorbidity and suicidal behavior in bipolar disorder

Bipolar disorder is often comorbid with anxiety, which is itself associated with poorer clinical outcomes, including suicide. A better etiologic understanding of this comorbidity could inform diagnosis and treatment. The present study aims to test whether comorbid anxiety in bipolar disorder reflects shared genetic risk factors. We also sought to assess the contribution of genetic risk for anxiety to suicide attempts in bipolar disorder. Polygenic risk scores (PRS) were calculated from published genome-wide association studies of samples of controls and cases with anxiety (n = 83,566) or bipolar disorder (n = 51,710), then scored in independent target samples (total n = 3369) of individuals with bipolar disorder who reported or denied lifetime anxiety disorders or suicidal attempts in research interviews. Participants were recruited from clinical and nonclinical settings and genotyped for common genetic variants. The results show that polygenic risk for anxiety was associated with comorbid anxiety disorders and suicide attempts in bipolar disorder, while polygenic risk for bipolar disorder was not associated with any of these variables. Our findings point out that comorbid anxiety disorders in bipolar disorder reflect a dual burden of bipolar and anxiety-related genes; the latter may also contribute to suicide attempts. Clinical care that recognizes and addresses this dual burden may help improve outcomes in people living with comorbid bipolar and anxiety disorders.

Genomics of esophageal squamous cell carcinoma in African Americans

e16552

Background: Esophageal cancer (EC) is one of the most lethal cancers with 90% mortality rate. Uneven Worldwide geographical distribution and varying etiologic factors contribute to the heterogenic nature of EC worldwide. In the US, 75% esophageal squamous cell carcinoma (ESCC) is diagnosed in African Americans (AA). However, despite comprehensive and multi-geographical efforts to understand the genomic basis of ESCC, the knowledge of molecular and genetic mechanisms leading to genomic subtyping of AA has been limited due to the underrepresentation of AAs in these studies. We hypothesized that the differences in mutational events in AA ESCC might, in part, cause the disparity in outcome and the mutational events may determine treatment options. Methods: Whole exome sequencing (WES) of matched- tumor and normal-cell DNA from late-stage ESCC of ten AA patients were performed by using Agilent SureSelect XT Human All Exon V6+UTR. Somatic variants per tumor samples called by two or more methods (Mutect2, VarScan2, and Strelka2) were combined and used for downstream population filtering. Somatic copy number changes were determined by using CNVKit in SevenBridges Genomics. Results: Somatic variants called by two or more algorithms filtered and sorted for rare mutations (equal to or less than 1% African originated population frequency) demonstrated an average of 23 nonsynonymous mutations per MB in high mutation rate AA ESCC. The remaining four samples had two or fewer nonsynonymous mutations per MB. TP53 was one of the most frequently mutated gene with 50% of mutation frequency. ARID2 mutated in 30%, EP300 in 20%, and RB1 in 10% of our cohort. Multiple amplified and deleted regions ranging from 94 to 46 were observed in seven samples in contrast to three samples that were silent in terms of copy number variations ranging from 29 to 9. Significant CNV were mostly seen in proliferation, cell cycle and checkpoint genes, squamous cell homeostasis, epithelial to mesenchymal transition, invasion and metastasis, receptor tyrosine kinase and signaling pathways, chromatin remodeling, detoxification, RNA/DNA editing and angiogenesis. Conclusions: WES of ten AA ESCC samples demonstrated a higher mutation rate in a group of patients with many passenger mutations and complex and recurrent copy number changes that affect oncogenic driver genes, which might suggest enhanced subtyping of ESCC in AA.

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.

Exploratory genome-wide association analysis of response to ketamine and a polygenic analysis of response to scopolamine in depression

Growing evidence suggests that the glutamatergic modulator ketamine has rapid antidepressant effects in treatment-resistant depressed subjects. The anticholinergic agent scopolamine has also shown promise as a rapid-acting antidepressant. This study applied genome-wide markers to investigate the role of genetic variants in predicting acute antidepressant response to both agents. The ketamine-treated sample included 157 unrelated European subjects with major depressive disorder (MDD) or bipolar disorder (BD). The scopolamine-treated sample comprised 37 unrelated European subjects diagnosed with either MDD or BD who had a current Major Depressive Episode (MDE), and had failed at least two adequate treatment trials for depression. Change in Montgomery-Asberg Depression Rating Scale (MADRS) or the 17-item Hamilton Depression Rating Scale (HAM-D) scale scores at day 1 (24 h post-treatment) was considered the primary outcome. Here, we conduct pilot genome-wide association study (GWAS) analyses to identify potential markers of ketamine response and dissociative side effects. Polygenic risk score analysis of SNPs ranked by the strength of their association with ketamine response was then calculated in order to assess whether common genetic markers from the ketamine study could predict response to scopolamine. Findings require replication in larger samples in light of low power of analyses of these small samples. Neverthless, these data provide a promising illustration of our future potential to identify genetic variants underlying rapid treatment response in mood disorders and may ultimately guide individual patient treatment selection in the future.

Basic Molecular Genetics Concepts and Tools

Progress in genetics in the last decades has been spurred largely by the Human Genome Project. Successful sequencing of the human genome and genomes of organisms across the evolutionary scale has been achieved. Scientific discoveries and innovations steered the refinement and annotation of sequences, cataloging of variation including SNPs and CNVs, development of fast, precise, high throughput methods and improved bioinformatics and statistical tools. These advances boosted GWAS on large samples identifying multiple susceptibility loci in diverse complex diseases and traits. Exome, whole genome and RNA sequencing has uncovered intricate elements of coding, noncoding and regulatory regions and the richness of the transcriptome. Evolution in genome editing has heightened the ability to alter any nucleotide and discern its biological effect. Exploiting these innovations in studies involving cells derived from patients through the iPSC technology offers manifold possibilities for functional studies that could enhance disease modeling and yield new therapeutic targets.

Analysis of the Influence of microRNAs in Lithium Response in Bipolar Disorder

Bipolar disorder (BD) is a common, highly heritable neuropsychiatric disease characterized by recurrent episodes of mania and depression. Lithium is the best-established long-term treatment for BD, even though individual response is highly variable. Evidence suggests that some of this variability has a genetic basis. This is supported by the largest genome-wide association study (GWAS) of lithium response to date conducted by the International Consortium on Lithium Genetics (ConLiGen). Recently, we performed the first genome-wide analysis of the involvement of miRNAs in BD and identified nine BD-associated miRNAs. However, it is unknown whether these miRNAs are also associated with lithium response in BD. In the present study, we therefore tested whether common variants at these nine candidate miRNAs contribute to the variance in lithium response in BD. Furthermore, we systematically analyzed whether any other miRNA in the genome is implicated in the response to lithium. For this purpose, we performed gene-based tests for all known miRNA coding genes in the ConLiGen GWAS dataset (n = 2,563 patients) using a set-based testing approach adapted from the versatile gene-based test for GWAS (VEGAS2). In the candidate approach, miR-499a showed a nominally significant association with lithium response, providing some evidence for involvement in both development and treatment of BD. In the genome-wide miRNA analysis, 71 miRNAs showed nominally significant associations with the dichotomous phenotype and 106 with the continuous trait for treatment response. A total of 15 miRNAs revealed nominal significance in both phenotypes with miR-633 showing the strongest association with the continuous trait (p = 9.80E-04) and miR-607 with the dichotomous phenotype (p = 5.79E-04). No association between miRNAs and treatment response to lithium in BD in either of the tested conditions withstood multiple testing correction. Given the limited power of our study, the investigation of miRNAs in larger GWAS samples of BD and lithium response is warranted.

Association of Polygenic Score for Schizophrenia and HLA Antigen and Inflammation Genes With Response to Lithium in Bipolar Affective Disorder: A Genome-Wide Association Study

IMPORTANCE

Lithium is a first-line mood stabilizer for the treatment of bipolar affective disorder (BPAD). However, the efficacy of lithium varies widely, with a nonresponse rate of up to 30%. Biological response markers are lacking. Genetic factors are thought to mediate treatment response to lithium, and there is a previously reported genetic overlap between BPAD and schizophrenia (SCZ).

OBJECTIVES

To test whether a polygenic score for SCZ is associated with treatment response to lithium in BPAD and to explore the potential molecular underpinnings of this association.

DESIGN, SETTING, AND PARTICIPANTS

A total of 2586 patients with BPAD who had undergone lithium treatment were genotyped and assessed for long-term response to treatment between 2008 and 2013. Weighted SCZ polygenic scores were computed at different P value thresholds using summary statistics from an international multicenter genome-wide association study (GWAS) of 36 989 individuals with SCZ and genotype data from patients with BPAD from the Consortium on Lithium Genetics. For functional exploration, a cross-trait meta-GWAS and pathway analysis was performed, combining GWAS summary statistics on SCZ and response to treatment with lithium. Data analysis was performed from September 2016 to February 2017.

MAIN OUTCOMES AND MEASURES

Treatment response to lithium was defined on both the categorical and continuous scales using the Retrospective Criteria of Long-Term Treatment Response in Research Subjects with Bipolar Disorder score. The effect measures include odds ratios and the proportion of variance explained.

RESULTS

Of the 2586 patients in the study (mean [SD] age, 47.2 [13.9] years), 1478 were women and 1108 were men. The polygenic score for SCZ was inversely associated with lithium treatment response in the categorical outcome, at a threshold P < 5 × 10-2. Patients with BPAD who had a low polygenic load for SCZ responded better to lithium, with odds ratios for lithium response ranging from 3.46 (95% CI, 1.42-8.41) at the first decile to 2.03 (95% CI, 0.86-4.81) at the ninth decile, compared with the patients in the 10th decile of SCZ risk. In the cross-trait meta-GWAS, 15 genetic loci that may have overlapping effects on lithium treatment response and susceptibility to SCZ were identified. Functional pathway and network analysis of these loci point to the HLA antigen complex and inflammatory cytokines.

CONCLUSIONS AND RELEVANCE

This study provides evidence for a negative association between high genetic loading for SCZ and poor response to lithium in patients with BPAD. These results suggest the potential for translational research aimed at personalized prescribing of lithium.

The antidepressant efficacy of subanesthetic-dose ketamine does not correlate with baseline subcortical volumes in a replication sample with major depressive disorder

BACKGROUND

This study sought to reproduce, in a larger sample, previous findings of a correlation between smaller raw 3-Tesla (3T) hippocampal volumes and improved antidepressant efficacy of ketamine in individuals with major depressive disorder (MDD). A secondary analysis stratified subjects according to functional BDNF rs6265 (val66met) genotype.

METHODS

Unmedicated subjects with treatment-resistant MDD ( n=55) underwent baseline structural 3T MRI. Data processing was conducted with FSL/FIRST and Freesurfer software. The amygdala, hippocampus, and thalamus were selected a priori for analysis. All subjects received a single 0.5mg/kg × 40-minute ketamine infusion. Pearson correlations were performed with subcortical volumes and percent change in MADRS score (from baseline to 230 minutes, 1 day, and 1 week post-infusion).

RESULTS

Raw and corrected subcortical volumes did not correlate with antidepressant response at any timepoint. In val/val subjects ( n=23), corrected left and right thalamic volume positively correlated with antidepressant response to ketamine at 230 minutes post-infusion but did not reach statistical significance. In met carriers ( n=14), corrected left and right thalamic volume negatively correlated with antidepressant response to ketamine.

CONCLUSION

Baseline subcortical volumes implicated in MDD did not correlate with ketamine's antidepressant efficacy. Baseline thalamic volume and BDNF genotype may be a combinatorial rapid antidepressant response biomarker.

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.

Genetic variants associated with response to lithium treatment in bipolar disorder: a genome-wide association study

BACKGROUND

Lithium is a first-line treatment in bipolar disorder, but individual response is variable. Previous studies have suggested that lithium response is a heritable trait. However, no genetic markers of treatment response have been reproducibly identified.

METHODS

Here, we report the results of a genome-wide association study of lithium response in 2563 patients collected by 22 participating sites from the International Consortium on Lithium Genetics (ConLiGen). Data from common single nucleotide polymorphisms (SNPs) were tested for association with categorical and continuous ratings of lithium response. Lithium response was measured using a well established scale (Alda scale). Genotyped SNPs were used to generate data at more than 6 million sites, using standard genomic imputation methods. Traits were regressed against genotype dosage. Results were combined across two batches by meta-analysis.

FINDINGS

A single locus of four linked SNPs on chromosome 21 met genome-wide significance criteria for association with lithium response (rs79663003, p=1·37 × 10(-8); rs78015114, p=1·31 × 10(-8); rs74795342, p=3·31 × 10(-9); and rs75222709, p=3·50 × 10(-9)). In an independent, prospective study of 73 patients treated with lithium monotherapy for a period of up to 2 years, carriers of the response-associated alleles had a significantly lower rate of relapse than carriers of the alternate alleles (p=0·03268, hazard ratio 3·8, 95% CI 1·1-13·0).

INTERPRETATION

The response-associated region contains two genes for long, non-coding RNAs (lncRNAs), AL157359.3 and AL157359.4. LncRNAs are increasingly appreciated as important regulators of gene expression, particularly in the CNS. Confirmed biomarkers of lithium response would constitute an important step forward in the clinical management of bipolar disorder. Further studies are needed to establish the biological context and potential clinical utility of these findings.

FUNDING

Deutsche Forschungsgemeinschaft, National Institute of Mental Health Intramural Research Program.

An Integrative Genomic Study Implicates the Postsynaptic Density in the Pathogenesis of Bipolar Disorder

Genome-wide association studies (GWAS) have identified several common variants associated with bipolar disorder (BD), but the biological meaning of these findings remains unclear. Integrative genomics-the integration of GWAS signals with gene expression data-may illuminate genes and gene networks that have key roles in the pathogenesis of BD. We applied weighted gene co-expression network analysis (WGCNA), which exploits patterns of co-expression among genes, to brain transcriptome data obtained by sequencing of poly-A RNA derived from postmortem dorsolateral prefrontal cortex from people with BD, along with age- and sex-matched controls. WGCNA identified 33 gene modules. Many of the modules corresponded closely to those previously reported in human cortex. Three modules were associated with BD, enriched for genes differentially expressed in BD, and also enriched for signals in prior GWAS of BD. Functional analysis of genes within these modules revealed significant enrichment of several functionally related sets of genes, especially those involved in the postsynaptic density (PSD). These results provide convergent support for the hypothesis that dysregulation of genes involved in the PSD is a key factor in the pathogenesis of BD. If replicated in larger samples, these findings could point toward new therapeutic targets for BD.

RNA-sequencing of the brain transcriptome implicates dysregulation of neuroplasticity, circadian rhythms and GTPase binding in bipolar disorder

RNA-sequencing (RNA-seq) is a powerful technique to investigate the complexity of gene expression in the human brain. We used RNA-seq to survey the brain transcriptome in high-quality postmortem dorsolateral prefrontal cortex from 11 individuals diagnosed with bipolar disorder (BD) and from 11 age- and gender-matched controls. Deep sequencing was performed, with over 350 million reads per specimen. At a false discovery rate of <5%, we detected five differentially expressed (DE) genes and 12 DE transcripts, most of which have not been previously implicated in BD. Among these, Prominin 1/CD133 and ATP-binding cassette-sub-family G-member2 (ABCG2) have important roles in neuroplasticity. We also show for the first time differential expression of long noncoding RNAs (lncRNAs) in BD. DE transcripts include those of serine/arginine-rich splicing factor 5 (SRSF5) and regulatory factor X4 (RFX4), which along with lncRNAs have a role in mammalian circadian rhythms. The DE genes were significantly enriched for several Gene Ontology categories. Of these, genes involved with GTPase binding were also enriched for BD-associated SNPs from previous genome-wide association studies, suggesting that differential expression of these genes is not simply a consequence of BD or its treatment. Many of these findings were replicated by microarray in an independent sample of 60 cases and controls. These results highlight common pathways for inherited and non-inherited influences on disease risk that may constitute good targets for novel therapies.

Molecular genetic overlap in bipolar disorder, schizophrenia, and major depressive disorder

OBJECTIVES

Genome-wide association studies (GWAS) in complex phenotypes, including psychiatric disorders, have yielded many replicated findings, yet individual markers account for only a small fraction of the inherited differences in risk. We tested the performance of polygenic models in discriminating between cases and healthy controls and among cases with distinct psychiatric diagnoses.

METHODS

GWAS results in bipolar disorder (BD), major depressive disorder (MDD), schizophrenia (SZ), and Parkinson's disease (PD) were used to assign weights to individual alleles, based on odds ratios. These weights were used to calculate allele scores for individual cases and controls in independent samples, summing across many single nucleotide polymorphisms (SNPs). How well allele scores discriminated between cases and controls and between cases with different disorders was tested by logistic regression.

RESULTS

Large sets of SNPs were needed to achieve even modest discrimination between cases and controls. The most informative SNPs were overlapping in BD, SZ, and MDD, with correlated effect sizes. Little or no overlap was seen between allele scores for psychiatric disorders and those for PD.

CONCLUSIONS

BD, SZ, and MDD all share a similar polygenic component, but the polygenic models tested lack discriminative accuracy and are unlikely to be useful for clinical diagnosis.

Race, genetic ancestry and response to antidepressant treatment for major depression

The Sequenced Treatment Alternatives to Relieve Depression (STAR*D) Study revealed poorer antidepressant treatment response among black compared with white participants. This racial disparity persisted even after socioeconomic and baseline clinical factors were taken into account. Some studies have suggested genetic contributions to this disparity, but none have attempted to disentangle race and genetic ancestry. Here we used genome-wide single-nucleotide polymorphism (SNP) data to examine independent contributions of race and genetic ancestry to citalopram response. Secondary data analyses included 1877 STAR*D participants who completed an average of 10 weeks of citalopram treatment and provided DNA samples. Participants reported their race as White (n=1464), black (n=299) or other/mixed (n=114). Genetic ancestry was estimated by multidimensional scaling (MDS) analyses of about 500 000 SNPs. Ancestry proportions were estimated by STRUCTURE. Structural equation modeling was used to examine the direct and indirect effects of observed and latent predictors of response, defined as change in the Quick Inventory of Depressive Symptomatology (QIDS) score from baseline to exit. Socioeconomic and baseline clinical factors, race, and anxiety significantly predicted response, as previously reported. However, direct effects of race disappeared in all models that included genetic ancestry. Genetic African ancestry predicted lower treatment response in all models. Although socioeconomic and baseline clinical factors drive racial differences in antidepressant response, genetic ancestry, rather than self-reported race, explains a significant fraction of the residual differences. Larger samples would be needed to identify the specific genetic mechanisms that may be involved, but these findings underscore the importance of including more African-American patients in drug trials.

Association study of 83 candidate genes for bipolar disorder in chromosome 6q selected using an evidence-based prioritization algorithm

BACKGROUND

Prior genome-scans of bipolar disorder have revealed chromosome 6q22 as a promising candidate region. However, linkage disequilibrium (LD) mapping studies have yet to identify replicated susceptibility loci.

METHODS

We analyzed 1,422 LD-tagging single nucleotide polymorphisms (SNPs) in 83 genes to test single-marker and locus-wide evidence of association with bipolar disorder in the NIMH Genetics Initiative bipolar pedigrees and the Portuguese Island Collection (PIC) (N = 1,093 in 528 informative pairs). Both studies previously demonstrated significant evidence of linkage to 6q. SNPs were genotyped using an Illumina iSelect genotyping array which employs the Infinium assay. Evidence of single-marker association was assessed using the generalized disequilibrium test (GDT). Empirical estimates of gene-wide significance were obtained by permutation (via 100,000 gene-dropping simulations) of Fisher's combined test of P-values for each locus.

RESULTS

No single variant yielded significant experiment-wide evidence of association, for either the combined sample or in each subsample. Our gene-dropping simulations identified nominally significant gene-wide associations with multiple loci, of which NT5DC1 in the NIMH subsample and CCNC in the PIC were the strongest candidates. However, no one gene consistently exceeded empirical significance criteria in both independent samples or survived Bonferroni correction for the number of genes tested.

CONCLUSIONS

Using a gene-based approach to family-based association, we identified gene-wide associations with several genes, though no single locus was significantly associated with bipolar disorder in both cohorts. This suggests that chromosome 6q may harbor multiple susceptibility loci or that complex patterns of LD in this region may confound approaches based on common SNPs. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.

Cancer regression and neurological toxicity following anti-MAGE-A3 TCR gene therapy

Nine cancer patients were treated with adoptive cell therapy using autologous anti-MAGE-A3 T-cell receptors (TCR)-engineered T cells. Five patients experienced clinical regression of their cancers including 2 on-going responders. Beginning 1-2 days postinfusion, 3 patients (#'s 5, 7, and 8) experienced mental status changes, and 2 patients (5 and 8) lapsed into comas and subsequently died. Magnetic resonance imagining analysis of patients 5 and 8 demonstrated periventricular leukomalacia, and examination of their brains at autopsy revealed necrotizing leukoencephalopathy with extensive white matter defects associated with infiltration of CD3(+)/CD8(+) T cells. Patient 7, developed Parkinson-like symptoms, which resolved over 4 weeks and fully recovered. Immunohistochemical staining of patient and normal brain samples demonstrated rare positively staining neurons with an antibody that recognizes multiple MAGE-A family members. The TCR used in this study recognized epitopes in MAGE-A3/A9/A12. Molecular assays of human brain samples using real-time quantitative-polymerase chain reaction, Nanostring quantitation, and deep-sequencing indicated that MAGE-A12 was expressed in human brain (and possibly MAGE-A1, MAGE-A8, and MAGE-A9). This previously unrecognized expression of MAGE-A12 in human brain was possibly the initiating event of a TCR-mediated inflammatory response that resulted in neuronal cell destruction and raises caution for clinical applications targeting MAGE-A family members with highly active immunotherapies.

Assessment of Response to Lithium Maintenance Treatment in Bipolar Disorder: A Consortium on Lithium Genetics (ConLiGen) Report

OBJECTIVE

The assessment of response to lithium maintenance treatment in bipolar disorder (BD) is complicated by variable length of treatment, unpredictable clinical course, and often inconsistent compliance. Prospective and retrospective methods of assessment of lithium response have been proposed in the literature. In this study we report the key phenotypic measures of the "Retrospective Criteria of Long-Term Treatment Response in Research Subjects with Bipolar Disorder" scale currently used in the Consortium on Lithium Genetics (ConLiGen) study.

MATERIALS AND METHODS

Twenty-nine ConLiGen sites took part in a two-stage case-vignette rating procedure to examine inter-rater agreement [Kappa (κ)] and reliability [intra-class correlation coefficient (ICC)] of lithium response. Annotated first-round vignettes and rating guidelines were circulated to expert research clinicians for training purposes between the two stages. Further, we analyzed the distributional properties of the treatment response scores available for 1,308 patients using mixture modeling.

RESULTS

Substantial and moderate agreement was shown across sites in the first and second sets of vignettes (κ = 0.66 and κ = 0.54, respectively), without significant improvement from training. However, definition of response using the A score as a quantitative trait and selecting cases with B criteria of 4 or less showed an improvement between the two stages (ICC1 = 0.71 and ICC2 = 0.75, respectively). Mixture modeling of score distribution indicated three subpopulations (full responders, partial responders, non responders).

CONCLUSIONS

We identified two definitions of lithium response, one dichotomous and the other continuous, with moderate to substantial inter-rater agreement and reliability. Accurate phenotypic measurement of lithium response is crucial for the ongoing ConLiGen pharmacogenomic study.

Genome-wide association study meta-analysis of European and Asian-ancestry samples identifies three novel loci associated with bipolar disorder

Meta-analyses of bipolar disorder (BD) genome-wide association studies (GWAS) have identified several genome-wide significant signals in European-ancestry samples, but so far account for little of the inherited risk. We performed a meta-analysis of ∼750,000 high-quality genetic markers on a combined sample of ∼14,000 subjects of European and Asian-ancestry (phase I). The most significant findings were further tested in an extended sample of ∼17,700 cases and controls (phase II). The results suggest novel association findings near the genes TRANK1 (LBA1), LMAN2L and PTGFR. In phase I, the most significant single nucleotide polymorphism (SNP), rs9834970 near TRANK1, was significant at the P=2.4 × 10(-11) level, with no heterogeneity. Supportive evidence for prior association findings near ANK3 and a locus on chromosome 3p21.1 was also observed. The phase II results were similar, although the heterogeneity test became significant for several SNPs. On the basis of these results and other established risk loci, we used the method developed by Park et al. to estimate the number, and the effect size distribution, of BD risk loci that could still be found by GWAS methods. We estimate that >63,000 case-control samples would be needed to identify the ∼105 BD risk loci discoverable by GWAS, and that these will together explain <6% of the inherited risk. These results support previous GWAS findings and identify three new candidate genes for BD. Further studies are needed to replicate these findings and may potentially lead to identification of functional variants. Sample size will remain a limiting factor in the discovery of common alleles associated with BD.

Replication and meta-analysis of TMEM132D gene variants in panic disorder

A recent genome-wide association study in patients with panic disorder (PD) identified a risk haplotype consisting of two single-nucleotide polymorphisms (SNPs) (rs7309727 and rs11060369) located in intron 3 of TMEM132D to be associated with PD in three independent samples. Now we report a subsequent confirmation study using five additional PD case-control samples (n = 1670 cases and n = 2266 controls) assembled as part of the Panic Disorder International Consortium (PanIC) study for a total of 2678 cases and 3262 controls in the analysis. In the new independent samples of European ancestry (EA), the association of rs7309727 and the risk haplotype rs7309727-rs11060369 was, indeed, replicated, with the strongest signal coming from patients with primary PD, that is, patients without major psychiatric comorbidities (n = 1038 cases and n = 2411 controls). This finding was paralleled by the results of the meta-analysis across all samples, in which the risk haplotype and rs7309727 reached P-levels of P = 1.4e-8 and P = 1.1e-8, respectively, when restricting the samples to individuals of EA with primary PD. In the Japanese sample no associations with PD could be found. The present results support the initial finding that TMEM132D gene contributes to genetic susceptibility for PD in individuals of EA. Our results also indicate that patient ascertainment and genetic background could be important sources of heterogeneity modifying this association signal in different populations.

Identity-by-descent filtering as a tool for the identification of disease alleles in exome sequence data from distant relatives

Large-scale, deep resequencing may be the next logical step in the genetic investigation of common complex diseases. Because each individual is likely to carry many thousands of variants, the identification of causal alleles requires an efficient strategy to reduce the number of candidate variants. Under many genetic models, causal alleles can be expected to reside within identity-by-descent (IBD) regions shared by affected relatives. In distant relatives, IBD regions constitute a small portion of the genome and can thus greatly reduce the search space for causal alleles. However, the effectiveness of this strategy is unknown. We test the simulated mini-exome data set in extended pedigrees provided by Genetic Analysis Workshop 17. At the fourth- and fifth-degree level of relatedness, case-case pairs shared between 1% and 9% of the genome identical by descent. As expected, no genes were shared identical by descent by all case subjects, but 43 genes were shared by many case subjects across at least 50 replicates. We filtered variants in these genes based on population frequency, function, informativeness, and evidence of association using the family-based association test. This analysis highlighted five genes previously implicated in triglyceride, lipid, and cholesterol metabolism. Comparison with the list of true risk alleles revealed that strict IBD filtering followed by association testing of the rarest alleles was the most sensitive strategy. IBD filtering may be a useful strategy for narrowing down the list of candidate variants in exome data, but the optimal degree of relatedness of affected pairs will depend on the genetic architecture of the disease under study.

A network-based approach to prioritize results from genome-wide association studies

Genome-wide association studies (GWAS) are a valuable approach to understanding the genetic basis of complex traits. One of the challenges of GWAS is the translation of genetic association results into biological hypotheses suitable for further investigation in the laboratory. To address this challenge, we introduce Network Interface Miner for Multigenic Interactions (NIMMI), a network-based method that combines GWAS data with human protein-protein interaction data (PPI). NIMMI builds biological networks weighted by connectivity, which is estimated by use of a modification of the Google PageRank algorithm. These weights are then combined with genetic association p-values derived from GWAS, producing what we call ‘trait prioritized sub-networks.’ As a proof of principle, NIMMI was tested on three GWAS datasets previously analyzed for height, a classical polygenic trait. Despite differences in sample size and ancestry, NIMMI captured 95% of the known height associated genes within the top 20% of ranked sub-networks, far better than what could be achieved by a single-locus approach. The top 2% of NIMMI height-prioritized sub-networks were significantly enriched for genes involved in transcription, signal transduction, transport, and gene expression, as well as nucleic acid, phosphate, protein, and zinc metabolism. All of these sub-networks were ranked near the top across all three height GWAS datasets we tested. We also tested NIMMI on a categorical phenotype, Crohn’s disease. NIMMI prioritized sub-networks involved in B- and T-cell receptor, chemokine, interleukin, and other pathways consistent with the known autoimmune nature of Crohn’s disease. NIMMI is a simple, user-friendly, open-source software tool that efficiently combines genetic association data with biological networks, translating GWAS findings into biological hypotheses.

A systems approach to the biology of mood disorders through network analysis of candidate genes

Meta analysis of association data of mood disorders has shown evidence for the role of particular genes in genetic risk. Integration of association data from meta analysis with differential expression data in brains of mood disorder patients could heighten the level of support for specific genes. To identify molecular mechanisms that may be disrupted in disease, a systems approach that involves analysis of biological networks created by these selected genes was employed.Interaction networks of hierarchical groupings of selected genes were generated using the Michigan Molecular Interactions (MiMI) software. Large networks were deconvoluted into subclusters of core complexes by using a community clustering program, GLay. Network nodes were functionally annotated in DAVID Bioinformatics Resource to identify enriched pathways and functional clusters. MAPK and beta adrenergic receptor signaling pathways were significantly enriched in the ANK3 and CACNA1C network. The PBRM1 network bolstered the enrichment of chromatin remodeling and transcription regulation functional clusters. Lowering the stringency for inclusion of other genes in network seeds increased network complexity and expanded the recruitment of enriched pathways to include signaling by neurotransmitter and hormone receptors, neurotrophin, ErbB and the cell cycle. We present a strategy to interrogate mechanisms in the cellular system that might be perturbed in disease. Network analysis of meta analysis- generated candidate genes that exhibited differential expression in mood disorder brains identified signaling pathways and functional clusters that may be involved in genetic risk for mood disorders.

Interaction networks of lithium and valproate molecular targets reveal a striking enrichment of apoptosis functional clusters and neurotrophin signaling

The overall neurobiological mechanisms by which lithium and valproate stabilize mood in bipolar disorder patients have yet to be fully defined. The therapeutic efficacy and dissimilar chemical structures of these medications suggest that they perturb both shared and disparate cellular processes. To investigate key pathways and functional clusters involved in the global action of lithium and valproate, we generated interaction networks formed by well-supported drug targets. Striking functional similarities emerged. Intersecting nodes in lithium and valproate networks highlighted a strong enrichment of apoptosis clusters and neurotrophin signaling. Other enriched pathways included MAPK, ErbB, insulin, VEGF, Wnt and long-term potentiation indicating a widespread effect of both drugs on diverse signaling systems. MAPK1/3 and AKT1/2 were the most preponderant nodes across pathways suggesting a central role in mediating pathway interactions. The convergence of biological responses unveils a functional signature for lithium and valproate that could be key modulators of their therapeutic efficacy.

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.

The International Consortium on Lithium Genetics (ConLiGen): an initiative by the NIMH and IGSLI to study the genetic basis of response to lithium treatment

For more than half a decade, lithium has been successfully used to treat bipolar disorder. Worldwide, it is considered the first-line mood stabilizer. Apart from its proven antimanic and prophylactic effects, considerable evidence also suggests an antisuicidal effect in affective disorders. Lithium is also effectively used to augment antidepressant drugs in the treatment of refractory major depressive episodes and prevent relapses in recurrent unipolar depression. In contrast to many psychiatric drugs, lithium has outlasted various pharmacotherapeutic 'fashions', and remains an indispensable element in contemporary psychopharmacology. Nevertheless, data from pharmacogenetic studies of lithium are comparatively sparse, and these studies are generally characterized by small sample sizes and varying definitions of response. Here, we present an international effort to elucidate the genetic underpinnings of lithium response in bipolar disorder. Following an initiative by the International Group for the Study of Lithium-Treated Patients (www.IGSLI.org) and the Unit on the Genetic Basis of Mood and Anxiety Disorders at the National Institute of Mental Health,lithium researchers from around the world have formed the Consortium on Lithium Genetics (www.ConLiGen.org) to establish the largest sample to date for genome-wide studies of lithium response in bipolar disorder, currently comprising more than 1,200 patients characterized for response to lithium treatment. A stringent phenotype definition of response is one of the hallmarks of this collaboration. ConLiGen invites all lithium researchers to join its efforts.

Meta-analysis of genome-wide association data identifies a risk locus for major mood disorders on 3p21.1

The major mood disorders, which include bipolar disorder (BD) and major depressive disorder (MDD), are substantially heritable, but few risk loci have been identified. We performed a meta-analysis of 5 major mood disorder case-control samples, including over 13,600 unique individuals genotyped with approximately 500,000 to 1 million single nucleotide polymorphism (SNP) markers on high-density arrays. Allele-wise association results were meta-analyzed with a method that weights results by sample size. We found genome-wide significant evidence that SNPs in a region of chromosome 3p21.1were associated with major mood disorders. The SNP rs2251219 returned the smallest meta-analysis p-value, 3.63 × 10⁻⁸, with a pooled odds ratio of 0.87. Supportive results were observed in 2 out of 3 independent samples tested in a replication study. These results implicate one or more genes in this region in the etiology of major mood disorders and suggest that BD and MDD share genetic risk factors.

Two variants in Ankyrin 3 (ANK3) are independent genetic risk factors for bipolar disorder

Two recent reports have highlighted ANK3 as a susceptibility gene for bipolar disorder (BD). We first reported association between BD and the ANK3 marker rs9804190 in a genome-wide association study (GWAS) of two independent samples (Baum et al., 2008). Subsequently, a meta-analysis of GWAS data based on samples from the US and the UK reported association with a different ANK3 marker, rs10994336 (Ferreira et al., 2008). The markers lie about 340 kb apart in the gene. Here, we test both markers in additional samples and characterize the contribution of each marker to BD risk. Our previously reported findings at rs9804190, which had been based on DNA pooling, were confirmed by individual genotyping in the National Institute of Mental Health (NIMH) waves 1-4 (P=0.05; odds ratio (OR)=1.24) and German (P=0.0006; OR=1.34) samples. This association was replicated in an independent US sample known as NIMH wave 5 (466 cases, 212 controls; P=0.017; OR=1.38). A random-effects meta-analysis of all three samples was significant (P=3 x 10(-6); OR=1.32), with no heterogeneity. Individual genotyping of rs10994336 revealed a significant association in the German sample (P=0.0001; OR=1.70), and similar ORs in the NIMH 1-4 and NIMH 5 samples that were not significant at the P<0.05 level. Meta-analysis of all three samples supported an association with rs10994336 (P=1.7 x 10(-5); OR=1.54), again with no heterogeneity. There was little linkage disequilibrium between the two markers. Further analysis suggested that each marker contributed independently to BD, with no significant marker x marker interaction. Our findings strongly support ANK3 as a BD susceptibility gene and suggest true allelic heterogeneity.

A genome-wide association study implicates diacylglycerol kinase eta (DGKH) and several other genes in the etiology of bipolar disorder

The genetic basis of bipolar disorder has long been thought to be complex, with the potential involvement of multiple genes, but methods to analyze populations with respect to this complexity have only recently become available. We have carried out a genome-wide association study of bipolar disorder by genotyping over 550,000 single-nucleotide polymorphisms (SNPs) in two independent case-control samples of European origin. The initial association screen was performed using pooled DNA, and selected SNPs were confirmed by individual genotyping. While DNA pooling reduces power to detect genetic associations, there is a substantial cost saving and gain in efficiency. A total of 88 SNPs, representing 80 different genes, met the prior criteria for replication in both samples. Effect sizes were modest: no single SNP of large effect was detected. Of 37 SNPs selected for individual genotyping, the strongest association signal was detected at a marker within the first intron of diacylglycerol kinase eta (DGKH; P=1.5 x 10(-8), experiment-wide P<0.01, OR=1.59). This gene encodes DGKH, a key protein in the lithium-sensitive phosphatidyl inositol pathway. This first genome-wide association study of bipolar disorder shows that several genes, each of modest effect, reproducibly influence disease risk. Bipolar disorder may be a polygenic disease.

Gene-based SNP mapping of a psychotic bipolar affective disorder linkage region on 22q12.3: association with HMG2L1 and TOM1

Genetic linkage studies in both bipolar affective disorder (BPAD) and schizophrenia have implicated overlapping regions of chromosome 22q. We previously reported that BPAD pedigrees containing multiple members with psychotic symptoms showed suggestive linkage to chromosome 22q12.3. Now we have tested 189 single nucleotide polymorphisms (SNPs) spanning a 3 Mb region around the linkage peak for association with BPAD in 305 families, unrelated cases, and controls. SNPs were selected in or near genes, resulting in coverage at a density of 1 SNP per 6.7 kb across the 22 annotated genes in the region. The strongest signal emerged from family-based association analysis of an 11-SNP, 54 kb haplotype straddling the gene HMG2L1 and part of TOM1. A 3-marker haplotype of SNPs within TOM1 was associated with BPAD (allele-wise P = 0.0011) and with psychotic BPAD (allele-wise P = 0.00049). As hypothesized, the mean odds ratio for the risk alleles across the region was 1.39 in the psychotic but only 0.96 in the non-psychotic subset. Genotype-wise analyses yielded similar results, but the psychotic/non-psychotic distinction was more pronounced with mean odds ratios of 1.91 versus 0.8. Permutation of genotype-wise results for rs2413338 in HMG2L1 showed an empirical P = 0.037 for the difference between subsets. HMG2L1 is a negative regulator of Wnt signaling, a pathway of interest in psychotic BPAD as it is activated by both mood stabilizer and anti-psychotic medications. Further work is needed to confirm these results and uncover the functional variation underlying the association signal.

Sequence variation in DOCK9 and heterogeneity in bipolar disorder

BACKGROUND

Linkage of bipolar disorder to a broad region on chromosome 13q has been supported in several studies including a meta-analysis on genome scans. Subsequent reports have shown that variations in the DAOA (G72) locus on 13q33 display association with bipolar disorder but these may not account for all of the linkage evidence in the region.

OBJECTIVE

To identify additional susceptibility loci on 13q32-q33 by linkage disequilibrium mapping and explore the impact of phenotypic heterogeneity on association.

METHODS

In the initial phase, 98 single nucleotide polymorphism (SNPs) located on 13q32-q33 were genotyped on 285 probands with bipolar disorder and their parents were drawn from families in the NIMH Genetics Initiative consortium for bipolar disorder (NIMH1-4) and two other series. Fine scale mapping using one family series (NIMH1-2) as the test sample was targeted on a gene that displayed the highest evidence of association. A secondary analysis of familial component phenotypes of bipolar disorder was conducted.

RESULTS

Three of seven SNPs in DOCK9, a gene that encodes an activator of the Rho-GTPase Cdc42, showed significant excess allelic transmission (P=0.0477-0.00067). Fine scale mapping on DOCK9 yielded evidence of association at nine SNPs in the gene (P=0.02-0.006). Follow-up tests detected excess transmission of the same allele of rs1340 in two out of three other sets of families. The association signals were largely attributable to maternally transmitted alleles (rs1927568: P=0.000083; odds ratio=3.778). A secondary analysis of familial component phenotypes of bipolar disorder detected significant association across multiple DOCK9 markers for racing thoughts, psychosis, delusion during mania and course of illness indicators.

CONCLUSION

These results suggest that DOCK9 contributes to both risk and increased illness severity in bipolar disorder. We found evidence for the effect of phenotypic heterogeneity on association. To our knowledge this is the first report to implicate DOCK9 or the Rho-GTPase pathway in the etiology of bipolar disorder.