The genetics of severe depression.

Genome-wide association studies (GWASs) of major depressive disorder (MDD) have recently achieved extremely large sample sizes and yielded substantial numbers of genome-wide significant loci. Because of the approach to ascertainment and assessment in many of these studies, some of these loci appear to be associated with dysphoria rather than with MDD, potentially decreasing the clinical relevance of the findings. An alternative approach to MDD GWAS is to focus on the most severe forms of MDD, with the hope that this will enrich for loci of larger effect, rendering their identification plausible, and providing potentially more clinically actionable findings. Here we review the genetics of severe depression by using clinical markers of severity including: age of onset, recurrence, degree of impairment, and treatment with ECT. There is evidence for increased family-based and Single Nucleotide Polymorphism (SNP)-based estimates of heritability in recurrent and early-onset illness as well as severe functional impariment. GWAS have been performed looking at severe forms of MDD and a few genome-wide loci have been identified. Several whole exome sequencing studies have also been performed, identifying associated rare variants. Although these findings have not yet been rigorously replicated, the elevated heritability seen in severe MDD phenotypes suggests the value of pursuing additional genome-wide interrogation of samples from this population. The challenge now is generating a cohort of adequate size with consistent phenotyping that will allow for careful and robust classifications and distinctions to be made. We are currently pursuing such a strategy in our 50-site worldwide Gen-ECT-ics consortium.

Genome-wide significant risk loci for mood disorders in the Old Order Amish founder population.

Genome-wide association studies (GWAS) of mood disorders in large case-control cohorts have identified numerous risk loci, yet pathophysiological mechanisms remain elusive, primarily due to the very small effects of common variants. We sought to discover risk variants with larger effects by conducting a genome-wide association study of mood disorders in a founder population, the Old Order Amish (OOA, n = 1,672). Our analysis revealed four genome-wide significant risk loci, all of which were associated with >2-fold relative risk. Quantitative behavioral and neurocognitive assessments (n = 314) revealed effects of risk variants on sub-clinical depressive symptoms and information processing speed. Network analysis suggested that OOA-specific risk loci harbor novel risk-associated genes that interact with known neuropsychiatry-associated genes via gene interaction networks. Annotation of the variants at these risk loci revealed population-enriched, non-synonymous variants in two genes encoding neurodevelopmental transcription factors, CUX1 and CNOT1. Our findings provide insight into the genetic architecture of mood disorders and a substrate for mechanistic and clinical studies.

Patterns of pharmacotherapy for bipolar disorder: A GBC survey.

OBJECTIVES: To understand treatment practices for bipolar disorders (BD), this study leveraged the Global Bipolar Cohort collaborative network to investigate pharmacotherapeutic treatment patterns in multiple cohorts of well-characterized individuals with BD in North America, Europe, and Australia. METHODS: Data on pharmacotherapy, demographics, diagnostic subtypes, and comorbidities were provided from each participating cohort. Individual site and regional pooled proportional meta-analyses with generalized linear mixed methods were conducted to identify prescription patterns. RESULTS: This study included 10,351 individuals from North America (n = 3985), Europe (n = 3822), and Australia (n = 2544). Overall, participants were predominantly female (60%) with BD-I (60%; vs. BD-II = 33%). Cross-sectionally, mood-stabilizing anticonvulsants (44%), second-generation antipsychotics (42%), and antidepressants (38%) were the most prescribed medications. Lithium was prescribed in 29% of patients, primarily in the Australian (31%) and European (36%) cohorts. First-generation antipsychotics were prescribed in 24% of the European versus 1% in the North American cohort. Antidepressant prescription rates were higher in BD-II (47%) compared to BD-I (35%). Major limitations were significant differences among cohorts based on inclusion/exclusion criteria, data source, and time/year of enrollment into cohort. CONCLUSIONS: Mood-stabilizing anticonvulsants, second-generation antipsychotics, and antidepressants were the most prescribed medications suggesting prescription patterns that are not necessarily guideline concordant. Significant differences exist in the prescription practices across different geographic regions, especially the underutilization of lithium in the North American cohorts and the higher utilization of first-generation antipsychotics in the European cohorts. There is a need to conduct future longitudinal studies to further explore these differences and their impact on outcomes, and to inform and implement evidence-based guidelines to help improve treatment practices in BD.

Rare variants implicate NMDA receptor signaling and cerebellar gene networks in risk for bipolar disorder.

Bipolar disorder is an often-severe mental health condition characterized by alternation between extreme mood states of mania and depression. Despite strong heritability and the recent identification of 64 common variant risk loci of small effect, pathophysiological mechanisms remain unknown. Here, we analyzed genome sequences from 41 multiply-affected pedigrees and identified variants in 741 genes with nominally significant linkage or association with bipolar disorder. These 741 genes overlapped known risk genes for neurodevelopmental disorders and clustered within gene networks enriched for synaptic and nuclear functions. The top variant in this analysis - prioritized by statistical association, predicted deleteriousness, and network centrality - was a missense variant in the gene encoding D-amino acid oxidase (DAOG131V). Heterologous expression of DAOG131V in human cells resulted in decreased DAO protein abundance and enzymatic activity. In a knock-in mouse model of DAOG131, DaoG130V/+, we similarly found decreased DAO protein abundance in hindbrain regions, as well as enhanced stress susceptibility and blunted behavioral responses to pharmacological inhibition of N-methyl-D-aspartate receptors (NMDARs). RNA sequencing of cerebellar tissue revealed that DaoG130V resulted in decreased expression of two gene networks that are enriched for synaptic functions and for genes expressed, respectively, in Purkinje neurons or granule neurons. These gene networks were also down-regulated in the cerebellum of patients with bipolar disorder compared to healthy controls and were enriched for additional rare variants associated with bipolar disorder risk. These findings implicate dysregulation of NMDAR signaling and of gene expression in cerebellar neurons in bipolar disorder pathophysiology and provide insight into its genetic architecture.

De novo variation in bipolar disorder.

Bipolar disorder (BD) is a common, highly heritable disorder that affects 1-2% of the world's population. To date, most genetic studies of BD have focused on common gene variation, and while robustly associated loci have been identified, a substantial proportion of the heritability remains missing and could be partially attributable to rare variation. In this study, we apply a de novo paradigm in BD to identify newly arisen variants that have yet to undergo natural selection and may represent highly pathogenic variants. We performed whole genome sequencing of 97 trios of Ashkenazi Jewish descent, selecting "simplex" families with no family history of BD and an early age of onset. We found a total of 6882 de novo variants (an average of 70.9 ± 12.9 S.D. variants per trio), including 107 variants within protein-coding genes. We combined our exonic variations with the results of 79 previously published BD trios, identifying 20 loss-of-function (LoF) and 77 missense damaging de novo variants in BD. These variants showed significant enrichment for constrained genes and for genes located to the postsynaptic density (PSD) (all Bonferroni corrected p < 0.05). Pathway analyses showed enrichment in several pathways, including "Phosphoinositides (PI) and their downstream targets" (Bonferroni p = 4.2 × 10-6), a pathway prominently featured in lithium's hypothesized mechanism of action. In addition, while we found overall evidence for transmission of common variant polygenic risk of BD in our full sample (pTDT p = 2.21 × 10-4), specific trios with LoF variants showed no evidence of polygenic transmission. In sum, our findings support the de novo paradigm as a contributor to the genetic architecture of BD and provide evidence that constrained genes, as well as genes within the PSD and PI pathway harbor rare variation associated with BD.

Circadian rhythms in bipolar disorder patient-derived neurons predict lithium response: preliminary studies.

Bipolar disorder (BD) is a neuropsychiatric illness defined by recurrent episodes of mania/hypomania, depression and circadian rhythm abnormalities. Lithium is an effective drug for BD, but 30-40% of patients fail to respond adequately to treatment. Previous work has demonstrated that lithium affects the expression of "clock genes" and that lithium responders (Li-R) can be distinguished from non-responders (Li-NR) by differences in circadian rhythms. However, circadian rhythms have not been evaluated in BD patient neurons from Li-R and Li-NR. We used induced pluripotent stem cells (iPSCs) to culture neuronal precursor cells (NPC) and glutamatergic neurons from BD patients characterized for lithium responsiveness and matched controls. We identified strong circadian rhythms in Per2-luc expression in NPCs and neurons from controls and Li-R, but NPC rhythms in Li-R had a shorter circadian period. Li-NR rhythms were low amplitude and profoundly weakened. In NPCs and neurons, expression of PER2 was higher in both BD groups compared to controls. In neurons, PER2 protein levels were higher in BD than controls, especially in Li-NR samples. In single cells, NPC and neuron rhythms in both BD groups were desynchronized compared to controls. Lithium lengthened period in Li-R and control neurons but failed to alter rhythms in Li-NR. In contrast, temperature entrainment increased amplitude across all groups, and partly restored rhythms in Li-NR neurons. We conclude that neuronal circadian rhythm abnormalities are present in BD and most pronounced in Li-NR. Rhythm deficits in BD may be partly reversible through stimulation of entrainment pathways.

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.

National Network of Depression Centers' Recommendations on Harmonizing Clinical Documentation of Electroconvulsive Therapy.

ABSTRACT: Electroconvulsive therapy (ECT) is a highly therapeutic and cost-effective treatment for severe and/or treatment-resistant major depression. However, because of the varied clinical practices, there is a great deal of heterogeneity in how ECT is delivered and documented. This represents both an opportunity to study how differences in implementation influence clinical outcomes and a challenge for carrying out coordinated quality improvement and research efforts across multiple ECT centers. The National Network of Depression Centers, a consortium of 26+ US academic medical centers of excellence providing care for patients with mood disorders, formed a task group with the goals of promoting best clinical practices for the delivery of ECT and to facilitate large-scale, multisite quality improvement and research to advance more effective and safe use of this treatment modality. The National Network of Depression Centers Task Group on ECT set out to define best practices for harmonizing the clinical documentation of ECT across treatment centers to promote clinical interoperability and facilitate a nationwide collaboration that would enable multisite quality improvement and longitudinal research in real-world settings. This article reports on the work of this effort. It focuses on the use of ECT for major depressive disorder, which accounts for the majority of ECT referrals in most countries. However, most of the recommendations on clinical documentation proposed herein will be applicable to the use of ECT for any of its indications.

Host-parasite interaction associated with major mental illness.

Clinical studies frequently report that patients with major mental illness such as schizophrenia and bipolar disorder have co-morbid physical conditions, suggesting that systemic alterations affecting both brain and peripheral tissues might underlie the disorders. Numerous studies have reported elevated levels of anti-Toxoplasma gondii (T. gondii) antibodies in patients with major mental illnesses, but the underlying mechanism was unclear. Using multidisciplinary epidemiological, cell biological, and gene expression profiling approaches, we report here multiple lines of evidence suggesting that a major mental illness-related susceptibility factor, Disrupted in schizophrenia (DISC1), is involved in host immune responses against T. gondii infection. Specifically, our cell biology and gene expression studies have revealed that DISC1 Leu607Phe variation, which changes DISC1 interaction with activating transcription factor 4 (ATF4), modifies gene expression patterns upon T. gondii infection. Our epidemiological data have also shown that DISC1 607 Phe/Phe genotype was associated with higher T. gondii antibody levels in sera. Although further studies are required, our study provides mechanistic insight into one of the few well-replicated serological observations in major mental illness.

Methylomic and transcriptomic predictors of one-month exposure to cortisol in healthy individuals.

Allostatic load (AL) refers to the cumulative "wear and tear" on an organism throughout its lifetime. One of the primary contributing factors to AL is prolonged exposure to stress or its primary catabolic agent cortisol. Chronic exposure to stress or cortisol is associated with numerous diseases, including cardiovascular disease, metabolic disorders, and psychiatric disorders. Therefore, a molecular marker capable of integrating a past history of cortisol exposure would be of great utility for assessing disease risk. To this end, we recruited 87 healthy males and females of European ancestry between 18 and 60 years old, extracted genomic DNA and RNA from leukocytes, and implemented a gene-centric DNA enrichment method coupled with bisulfite sequencing and RNA-Seq of total RNA for the determination of genome-wide methylation and gene transcription, respectively. Sequencing data were analyzed against awakening and bedtime cortisol data to identify differentially methylated regions (DMRs) and CpGs (DMCs) and differentially expressed genes (DEGs). Six candidate DMCs (punadjusted < 0.005) and nine DEGs (punadjusted < 0.0005) were used to construct a prediction model that could capture past 30+ days of both bedtime and awakening cortisol levels. Utilizing a cross-validation approach, we obtained a regression coefficient of R2 = 0.308 for predicting continuous awakening cortisol and an area under the curve (AUC) = 0.753 for dichotomous (high vs. low tertile) awakening cortisol, and R2 = 0.224 and AUC = 0.723 for continuous and dichotomous bedtime cortisol levels, respectively. To our knowledge, the current study represents the first attempt to identify genome-wide predictors of cortisol exposure that utilizes both methylation and transcription targets. The utility of our approach needs to be replicated in an independent cohort of samples for which similar cortisol metrics are available.

Differential responses to lithium in hyperexcitable neurons from patients with bipolar disorder.

Bipolar disorder is a complex neuropsychiatric disorder that is characterized by intermittent episodes of mania and depression; without treatment, 15% of patients commit suicide. Hence, it has been ranked by the World Health Organization as a top disorder of morbidity and lost productivity. Previous neuropathological studies have revealed a series of alterations in the brains of patients with bipolar disorder or animal models, such as reduced glial cell number in the prefrontal cortex of patients, upregulated activities of the protein kinase A and C pathways and changes in neurotransmission. However, the roles and causation of these changes in bipolar disorder have been too complex to exactly determine the pathology of the disease. Furthermore, although some patients show remarkable improvement with lithium treatment for yet unknown reasons, others are refractory to lithium treatment. Therefore, developing an accurate and powerful biological model for bipolar disorder has been a challenge. The introduction of induced pluripotent stem-cell (iPSC) technology has provided a new approach. Here we have developed an iPSC model for human bipolar disorder and investigated the cellular phenotypes of hippocampal dentate gyrus-like neurons derived from iPSCs of patients with bipolar disorder. Guided by RNA sequencing expression profiling, we have detected mitochondrial abnormalities in young neurons from patients with bipolar disorder by using mitochondrial assays; in addition, using both patch-clamp recording and somatic Ca(2+) imaging, we have observed hyperactive action-potential firing. This hyperexcitability phenotype of young neurons in bipolar disorder was selectively reversed by lithium treatment only in neurons derived from patients who also responded to lithium treatment. Therefore, hyperexcitability is one early endophenotype of bipolar disorder, and our model of iPSCs in this disease might be useful in developing new therapies and drugs aimed at its clinical treatment.

Developing Treatments for Alzheimer's and Related Disorders with Precision Medicine: A Vision.

Precision medicine, also known as personalized medicine, is concerned with finding the right treatment for the right patient at the right time. It is a way of thinking focused on parsing heterogeneity ultimately down to the level of the individual. Its main mission is to identify characteristics of heterogeneous clinical conditions so as to target tailored therapies to individuals. Precision Medicine however is not an agnostic collection of all manner of clinical, genetic and other biologic data in select cohorts. This is an important point. Simply collecting as much information as possible on individuals without applying this way of thinking should not be considered Precision Medicine.

Mutations in the pancreatic secretory enzymes CPA1 and CPB1 are associated with pancreatic cancer.

To evaluate whether germline variants in genes encoding pancreatic secretory enzymes contribute to pancreatic cancer susceptibility, we sequenced the coding regions of CPB1 and other genes encoding pancreatic secretory enzymes and known pancreatitis susceptibility genes (PRSS1, CPA1, CTRC, and SPINK1) in a hospital series of pancreatic cancer cases and controls. Variants in CPB1, CPA1 (encoding carboxypeptidase B1 and A1), and CTRC were evaluated in a second set of cases with familial pancreatic cancer and controls. More deleterious CPB1 variants, defined as having impaired protein secretion and induction of endoplasmic reticulum (ER) stress in transfected HEK 293T cells, were found in the hospital series of pancreatic cancer cases (5/986, 0.5%) than in controls (0/1,045, P = 0.027). Among familial pancreatic cancer cases, ER stress-inducing CPB1 variants were found in 4 of 593 (0.67%) vs. 0 of 967 additional controls (P = 0.020), with a combined prevalence in pancreatic cancer cases of 9/1,579 vs. 0/2,012 controls (P < 0.01). More ER stress-inducing CPA1 variants were also found in the combined set of hospital and familial cases with pancreatic cancer than in controls [7/1,546 vs. 1/2,012; P = 0.025; odds ratio, 9.36 (95% CI, 1.15-76.02)]. Overall, 16 (1%) of 1,579 pancreatic cancer cases had an ER stress-inducing CPA1 or CPB1 variant, compared with 1 of 2,068 controls (P < 0.00001). No other candidate genes had statistically significant differences in variant prevalence between cases and controls. Our study indicates ER stress-inducing variants in CPB1 and CPA1 are associated with pancreatic cancer susceptibility and implicate ER stress in pancreatic acinar cells in pancreatic cancer development.

Defining Major Depressive Disorder Cohorts Using the EHR: Multiple Phenotypes Based on ICD-9 Codes and Medication Orders.

BACKGROUND: Major Depressive Disorder (MDD) is one of the most common mental illnesses and a leading cause of disability worldwide. Electronic Health Records (EHR) allow researchers to conduct unprecedented large-scale observational studies investigating MDD, its disease development and its interaction with other health outcomes. While there exist methods to classify patients as clear cases or controls, given specific data requirements, there are presently no simple, generalizable, and validated methods to classify an entire patient population into varying groups of depression likelihood and severity. METHODS: We have tested a simple, pragmatic electronic phenotype algorithm that classifies patients into one of five mutually exclusive, ordinal groups, varying in depression phenotype. Using data from an integrated health system on 278,026 patients from a 10-year study period we have tested the convergent validity of these constructs using measures of external validation, including patterns of psychiatric prescriptions, symptom severity, indicators of suicidality, comorbidity, mortality, health care utilization, and polygenic risk scores for MDD. RESULTS: We found consistent patterns of increasing morbidity and/or adverse outcomes across the five groups, providing evidence for convergent validity. LIMITATIONS: The study population is from a single rural integrated health system which is predominantly white, possibly limiting its generalizability. CONCLUSION: Our study provides initial evidence that a simple algorithm, generalizable to most EHR data sets, provides categories with meaningful face and convergent validity that can be used for stratification of an entire patient population.

Targeted sequencing of the LRRTM gene family in suicide attempters with bipolar disorder.

Glutamatergic signaling is the primary excitatory neurotransmission pathway in the brain, and its relationship to neuropsychiatric disorders is of considerable interest. Our previous attempted suicide genome-wide association study, and numerous studies investigating gene expression, genetic variation, and DNA methylation have implicated aberrant glutamatergic signaling in suicide risk. The glutamatergic pathway gene LRRTM4 was an associated gene identified in our attempted suicide genome-wide association study, with association support seen primarily in females. Recent evidence has also shown that glutamatergic signaling is partly regulated by sex-related hormones. The LRRTM gene family encodes neuronal leucine-rich transmembrane proteins that localize to and promote glutamatergic synapse development. In this study, we sequenced the coding and regulatory regions of all four LRRTM gene members plus a large intronic region of LRRTM4 in 476 bipolar disorder suicide attempters and 473 bipolar disorder nonattempters. We identified two male-specific variants, one female- and five male-specific haplotypes significantly associated with attempted suicide in LRRTM4. Furthermore, variants within significant haplotypes may be brain expression quantitative trait loci for LRRTM4 and some of these variants overlap with predicted hormone response elements. Overall, these results provide supporting evidence for a sex-specific association of genetic variation in LRRTM4 with attempted suicide.

Differentially methylated regions in bipolar disorder and suicide.

The addition of a methyl group to, typically, a cytosine-guanine dinucleotide (CpG) creates distinct DNA methylation patterns across the genome that can regulate gene expression. Aberrant DNA methylation of CpG sites has been associated with many psychiatric disorders including bipolar disorder (BD) and suicide. Using the SureSelectXT system, Methyl-Seq, we investigated the DNA methylation status of CpG sites throughout the genome in 50 BD individuals (23 subjects who died by suicide and 27 subjects who died from other causes) and 31 nonpsychiatric controls. We identified differentially methylated regions (DMRs) from three analyses: (a) BD subjects compared to nonpsychiatric controls (BD-NC), (b) BD subjects who died by suicide compared to nonpsychiatric controls (BDS-NC), and (c) BDS subjects compared to BD subjects who died from other causes (BDS-BDNS). One DMR from the BDS-NC analysis, located in ARHGEF38, was significantly hypomethylated (23.4%) in BDS subjects. This finding remained significant after multiple testing (PBootstrapped = 9.0 × 10-3 ), was validated using pyrosequencing, and was more significant in males. A secondary analysis utilized Ingenuity Pathway Analysis to identify enrichment in nominally significant DMRs. This identified an association with several pathways including axonal guidance signaling, calcium signaling, ß-adrenergic signaling, and opioid signaling. Our comprehensive study provides further support that DNA methylation alterations influence the risk for BD and suicide. However, further investigation is required to confirm these associations and identify their functional consequences.

Vascular Pathology and Trajectories of Late-Life Major Depressive Disorder in Secondary Psychiatric Care.

OBJECTIVE: To examine 5-year trajectories of psychiatrist-treated late-life major depressive disorder (MDD), and evaluate whether previous vascular pathology is associated with more severe trajectories of late-life MDD. METHODS: Data were obtained from nationally representative civil, psychiatric, hospital, and prescription registers in Denmark. The sample included 11,092 older adults (≥60 years) who received their first diagnosis of MDD in a psychiatric facility in Denmark between 2000 and 2007. Trajectories of inpatient or outpatient contact at psychiatric hospitals for MDD over the 5-year period following index MDD diagnosis were modeled using latent class growth analysis. Measures of vascular disease (stroke, heart disease, vascular dementia) and vascular risk factors (hypertension, diabetes) were defined based on medication prescriptions and hospital-based diagnoses. Other predictors included demographic characteristics and characteristics of the index MDD diagnosis. RESULTS: The final model included 4 trajectories with consistently low (66% of the sample), high decreasing (19%), consistently high (9%), and moderate fluctuating (6%) probabilities of contact at a psychiatric hospital for MDD during the 5-year period following the index MDD diagnosis. We found no significant associations between any form of vascular pathology and trajectory class membership. Relative to the consistently low class, older age, greater severity and >12 months of prior antidepressant medication use predicted membership in the other three classes. CONCLUSIONS: A notable proportion (34%) of individuals diagnosed with MDD in late-life require secondary psychiatric treatment for extended time periods. We did not find evidence that vascular pathology predicts hospital contact trajectories in secondary-treated late-life MDD.

The Pharmacogenomics of Bipolar Disorder study (PGBD): identification of genes for lithium response in a prospective sample.

BACKGROUND: Bipolar disorder is a serious and common psychiatric disorder characterized by manic and depressive mood switches and a relapsing and remitting course. The cornerstone of clinical management is stabilization and prophylaxis using mood-stabilizing medications to reduce both manic and depressive symptoms. Lithium remains the gold standard of treatment with the strongest data for both efficacy and suicide prevention. However, many patients do not respond to this medication, and clinically there is a great need for tools to aid the clinician in selecting the correct treatment. Large genome wide association studies (GWAS) investigating retrospectively the effect of lithium response are in the pipeline; however, few large prospective studies on genetic predictors to of lithium response have yet been conducted. The purpose of this project is to identify genes that are associated with lithium response in a large prospective cohort of bipolar patients and to better understand the mechanism of action of lithium and the variation in the genome that influences clinical response. METHODS/DESIGN: This study is an 11-site prospective non-randomized open trial of lithium designed to ascertain a cohort of 700 subjects with bipolar I disorder who experience protocol-defined relapse prevention as a result of treatment with lithium monotherapy. All patients will be diagnosed using the Diagnostic Interview for Genetic Studies (DIGS) and will then enter a 2-year follow-up period on lithium monotherapy if and when they exhibit a score of 1 (normal, not ill), 2 (minimally ill) or 3 (mildly ill) on the Clinical Global Impressions of Severity Scale for Bipolar Disorder (CGI-S-BP Overall Bipolar Illness) for 4 of the 5 preceding weeks. Lithium will be titrated as clinically appropriate, not to exceed serum levels of 1.2 mEq/L. The sample will be evaluated longitudinally using a wide range of clinical scales, cognitive assessments and laboratory tests. On relapse, patients will be discontinued or crossed-over to treatment with valproic acid (VPA) or treatment as usual (TAU). Relapse is defined as a DSM-IV manic, major depressive or mixed episode or if the treating physician decides a change in medication is clinically necessary. The sample will be genotyped for GWAS. The outcome for lithium response will be analyzed as a time to event, where the event is defined as clinical relapse, using a Cox Proportional Hazards model. Positive single nucleotide polymorphisms (SNPs) from past genetic retrospective studies of lithium response, the Consortium on Lithium Genetics (ConLiGen), will be tested in this prospective study sample; a meta-analysis of these samples will then be performed. Finally, neurons will be derived from pluripotent stem cells from lithium responders and non-responders and tested in vivo for response to lithium by gene expression studies. SNPs in genes identified in these cellular studies will also be tested for association to response. DISCUSSION: Lithium is an extraordinarily important therapeutic drug in the clinical management of patients suffering from bipolar disorder. However, a significant proportion of patients, 30-40 %, fail to respond, and there is currently no method to identify the good lithium responders before initiation of treatment. Converging evidence suggests that genetic factors play a strong role in the variation of response to lithium, but only a few genes have been tested and the samples have largely been retrospective or quite small. The current study will collect an entirely unique sample of 700 patients with bipolar disorder to be stabilized on lithium monotherapy and followed for up to 2 years. This study will produce useful information to improve the understanding of the mechanism of action of lithium and will add to the development of a method to predict individual response to lithium, thereby accelerating recovery and reducing suffering and cost. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT01272531 Registered: January 6, 2011.

A hybrid likelihood model for sequence-based disease association studies.

In the past few years, case-control studies of common diseases have shifted their focus from single genes to whole exomes. New sequencing technologies now routinely detect hundreds of thousands of sequence variants in a single study, many of which are rare or even novel. The limitation of classical single-marker association analysis for rare variants has been a challenge in such studies. A new generation of statistical methods for case-control association studies has been developed to meet this challenge. A common approach to association analysis of rare variants is the burden-style collapsing methods to combine rare variant data within individuals across or within genes. Here, we propose a new hybrid likelihood model that combines a burden test with a test of the position distribution of variants. In extensive simulations and on empirical data from the Dallas Heart Study, the new model demonstrates consistently good power, in particular when applied to a gene set (e.g., multiple candidate genes with shared biological function or pathway), when rare variants cluster in key functional regions of a gene, and when protective variants are present. When applied to data from an ongoing sequencing study of bipolar disorder (191 cases, 107 controls), the model identifies seven gene sets with nominal p-values < 0.05, of which one MAPK signaling pathway (KEGG) reaches trend-level significance after correcting for multiple testing.