Genomic view of bipolar disorder revealed by whole genome sequencing in a genetic isolate

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.

A fast linkage method for population GWAS cohorts with related individuals

Linkage analysis, a class of methods for detecting co-segregation of genomic segments and traits in families, was used to map disease-causing genes for decades before genotyping arrays and dense SNP genotyping enabled genome-wide association studies in population samples. Population samples often contain related individuals, but the segregation of alleles within families is rarely used because traditional linkage methods are computationally inefficient for larger datasets. Here, we describe Population Linkage, a novel application of Haseman-Elston regression as a method of moments estimator of variance components and their standard errors. We achieve additional computational efficiency by using modern methods for detection of IBD segments and variance component estimation, efficient preprocessing of input data, and minimizing redundant numerical calculations. We also refined variance component models to account for the biases in population-scale methods for IBD segment detection. We ran Population Linkage on four blood lipid traits in over 70,000 individuals from the HUNT and SardiNIA studies, successfully detecting 25 known genetic signals. One notable linkage signal that appeared in both was for low-density lipoprotein (LDL) cholesterol levels in the region near the gene APOE (LOD = 29.3, variance explained = 4.1%). This is the region where the missense variants rs7412 and rs429358, which together make up the ε2, ε3, and ε4 alleles each account for 2.4% and 0.8% of variation in circulating LDL cholesterol. Our results show the potential for linkage analysis and other large-scale applications of method of moments variance components estimation.

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.

Genomic and neuroimaging approaches to bipolar disorder

BACKGROUND

To date, besides genome-wide association studies, a variety of other genetic analyses (e.g. polygenic risk scores, whole-exome sequencing and whole-genome sequencing) have been conducted, and a large amount of data has been gathered for investigating the involvement of common, rare and very rare types of DNA sequence variants in bipolar disorder. Also, non-invasive neuroimaging methods can be used to quantify changes in brain structure and function in patients with bipolar disorder.

AIMS

To provide a comprehensive assessment of genetic findings associated with bipolar disorder, based on the evaluation of different genomic approaches and neuroimaging studies.

METHOD

We conducted a PubMed search of all relevant literatures from the beginning to the present, by querying related search strings.

RESULTS

ANK3, CACNA1C, SYNE1, ODZ4 and TRANK1 are five genes that have been replicated as key gene candidates in bipolar disorder pathophysiology, through the investigated studies. The percentage of phenotypic variance explained by the identified variants is small (approximately 4.7%). Bipolar disorder polygenic risk scores are associated with other psychiatric phenotypes. The ENIGMA-BD studies show a replicable pattern of lower cortical thickness, altered white matter integrity and smaller subcortical volumes in bipolar disorder.

CONCLUSIONS

The low amount of explained phenotypic variance highlights the need for further large-scale investigations, especially among non-European populations, to achieve a more complete understanding of the genetic architecture of bipolar disorder and the missing heritability. Combining neuroimaging data with genetic data in large-scale studies might help researchers acquire a better knowledge of the engaged brain regions in bipolar disorder.

Convolutional Neural Network Visualization for Identification of Risk Genes in Bipolar Disorder

BACKGROUND

Bipolar disorder (BD) is a type of chronic emotional disorder with a complex genetic structure. However, its genetic molecular mechanism is still unclear, which makes it insufficient to be diagnosed and treated.

METHODS AND RESULTS

In this paper, we proposed a model for predicting BD based on single nucleotide polymorphisms (SNPs) screening by genome-wide association study (GWAS), which was constructed by a convolutional neural network (CNN) that predicted the probability of the disease. According to the difference of GWAS threshold, two sets of data were named: group P001 and group P005. And different convolutional neural networks are set for the two sets of data. The training accuracy of the model trained with group P001 data is 96%, and the test accuracy is 91%. The training accuracy of the model trained with group P005 data is 94.5%, and the test accuracy is 92%. At the same time, we used gradient weighted class activation mapping (Grad-CAM) to interpret the prediction model, indirectly to identify high-risk SNPs of BD. In the end, we compared these high-risk SNPs with human gene annotation information.

CONCLUSION

The model prediction results of the group P001 yielded 137 risk genes, of which 22 were reported to be associated with the occurrence of BD. The model prediction results of the group P005 yielded 407 risk genes, of which 51 were reported to be associated with the occurrence of BD.

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.

A linkage and exome study of multiplex families with bipolar disorder implicates rare coding variants of ANK3 and additional rare alleles at 10q11-q21

BACKGROUND

Bipolar disorder is a highly heritable psychiatric condition for which specific genetic factors remain largely unknown. In the present study, we used combined whole-exome sequencing and linkage analysis to identify risk loci and dissect the contribution of common and rare variants in families with a high density of illness.

METHODS

Overall, 117 participants from 15 Australian extended families with bipolar disorder (72 with affective disorder, including 50 with bipolar disorder type I or II, 13 with schizoaffective disorder-manic type and 9 with recurrent unipolar disorder) underwent whole-exome sequencing. We performed genome-wide linkage analysis using MERLIN and conditional linkage analysis using LAMP. We assessed the contribution of potentially functional rare variants using a genebased segregation test.

RESULTS

We identified a significant linkage peak on chromosome 10q11-q21 (maximal single nucleotide polymorphism = rs10761725; exponential logarithm of the odds [LODexp] = 3.03; empirical p = 0.046). The linkage interval spanned 36 protein-coding genes, including a gene associated with bipolar disorder, ankyrin 3 (ANK3). Conditional linkage analysis showed that common ANK3 risk variants previously identified in genome-wide association studies - or variants in linkage disequilibrium with those variants - did not explain the linkage signal (rs10994397 LOD = 0.63; rs9804190 LOD = 0.04). A family-based segregation test with 34 rare variants from 14 genes under the linkage interval suggested rare variant contributions of 3 brain-expressed genes: NRBF2 (p = 0.005), PCDH15 (p = 0.002) and ANK3 (p = 0.014).

LIMITATIONS

We did not examine non-coding variants, but they may explain the remaining linkage signal.

CONCLUSION

Combining family-based linkage analysis with next-generation sequencing data is effective for identifying putative disease genes and specific risk variants in complex disorders. We identified rare missense variants in ANK3, PCDH15 and NRBF2 that could confer disease risk, providing valuable targets for functional characterization.

Ancestral haplotype reconstruction in endogamous populations using identity-by-descent

In this work we develop a novel algorithm for reconstructing the genomes of ancestral individuals, given genotype or sequence data from contemporary individuals and an extended pedigree of family relationships. A pedigree with complete genomes for every individual enables the study of allele frequency dynamics and haplotype diversity across generations, including deviations from neutrality such as transmission distortion. When studying heritable diseases, ancestral haplotypes can be used to augment genome-wide association studies and track disease inheritance patterns. The building blocks of our reconstruction algorithm are segments of Identity-By-Descent (IBD) shared between two or more genotyped individuals. The method alternates between identifying a source for each IBD segment and assembling IBD segments placed within each ancestral individual. Unlike previous approaches, our method is able to accommodate complex pedigree structures with hundreds of individuals genotyped at millions of SNPs.

We apply our method to an Old Order Amish pedigree from Lancaster, Pennsylvania, whose founders came to North America from Europe during the early 18th century. The pedigree includes 1338 individuals from the past 12 generations, 394 with genotype data. The motivation for reconstruction is to understand the genetic basis of diseases segregating in the family through tracking haplotype transmission over time. Using our algorithm thread, we are able to reconstruct an average of 224 ancestral individuals per chromosome. For these ancestral individuals, on average we reconstruct 79% of their haplotypes. We also identify a region on chromosome 16 that is difficult to reconstruct—we find that this region harbors a short Amish-specific copy number variation and the gene HYDIN. thread was developed for endogamous populations, but can be applied to any extensive pedigree with the recent generations genotyped. We anticipate that this type of practical ancestral reconstruction will become more common and necessary to understand rare and complex heritable diseases in extended families.

When analyzing complex heritable traits, genomic data from many generations of an extended family increases the amount of information available for statistical inference. However, typically only genomic data from the recent generations of a pedigree are available, as ancestral individuals are deceased. In this work we present an algorithm, called thread, for reconstructing the genomes of ancestral individuals, given a complex pedigree and genomic data from the recent generations. Previous approaches have not been able to accommodate large datasets (both in terms of sites and individuals), made simplifying assumptions about pedigree structure, or did not tie reconstructed sequences back to specific individuals. We apply thread to a complex Old Order Amish pedigree of 1338 individuals, 394 with genotype data.

Inherited L1 Retrotransposon Insertions Associated With Risk for Schizophrenia and Bipolar Disorder

Studies of the genetic heritability of schizophrenia and bipolar disorder examining single nucleotide polymorphisms (SNPs) and copy number variations have failed to explain a large portion of the genetic liability, resulting in substantial missing heritability. Long interspersed element 1 (L1) retrotransposons are a type of inherited polymorphic variant that may be associated with risk for schizophrenia and bipolar disorder. We performed REBELseq, a genome wide assay for L1 sequences, on DNA from male and female persons with schizophrenia and controls (n = 63 each) to identify inherited L1 insertions and validated priority insertions. L1 insertions of interest were genotyped in DNA from a replication cohort of persons with schizophrenia, bipolar disorder, and controls (n = 2268 each) to examine differences in carrier frequencies. We identified an inherited L1 insertion in ARHGAP24 and a quadallelic SNP (rs74169643) inside an L1 insertion in SNTG2 that are associated with risk for developing schizophrenia and bipolar disorder (all odds ratios ~1.2). Pathway analysis identified 15 gene ontologies that were differentially affected by L1 burden, including multiple ontologies related to glutamatergic signaling and immune function, which have been previously associated with schizophrenia. These findings provide further evidence supporting the role of inherited repetitive genetic elements in the heritability of psychiatric disorders.

Transcriptomic signatures across human tissues identify functional rare genetic variation

Rare genetic variants are abundant across the human genome, and identifying their function and phenotypic impact is a major challenge. Measuring aberrant gene expression has aided in identifying functional, large-effect rare variants (RVs). Here, we expanded detection of genetically driven transcriptome abnormalities by analyzing gene expression, allele-specific expression, and alternative splicing from multitissue RNA-sequencing data, and demonstrate that each signal informs unique classes of RVs. We developed Watershed, a probabilistic model that integrates multiple genomic and transcriptomic signals to predict variant function, validated these predictions in additional cohorts and through experimental assays, and used them to assess RVs in the UK Biobank, the Million Veterans Program, and the Jackson Heart Study. Our results link thousands of RVs to diverse molecular effects and provide evidence to associate RVs affecting the transcriptome with human traits.

Validity of the Mood Disorder Questionnaire (MDQ) as a screening tool for bipolar spectrum disorders in anabaptist populations

The Mood Disorder Questionnaire (MDQ) is an established screening tool for bipolar spectrum disorders (BSD), but has not been validated in diverse populations and the best scoring method remains uncertain. This study assessed diagnostic validity of the MDQ among Anabaptists, an underserved population frequently involved in genetic research. 161 participants completed the MDQ and were diagnosed by a best-estimate final diagnosis (BEFD). Diagnostic agreements between alternate MDQ scoring methods and the BEFD were quantified using Cohen's Kappa (κ), sensitivity (α), and specificity (β). Scoring criteria evaluated included >7 simultaneous symptoms and at least moderate impairment, >7 simultaneous symptoms, with at least mild impairment, >7 symptoms only, with no further requirement, and three novel scoring methods that require >5 symptoms or fewer. Diagnostic agreement varied. The original method proved most specific but had the lowest κ and sensitivity. κ increased with more liberal scoring criteria, reaching a maximum under the lower-threshold symptom methods, with little loss of specificity in the 5-symptom method. Decreasing the symptom threshold below 5 conferred little or no benefit. These results support the diagnostic validity of the MDQ among this Anabaptist sample and suggest that a 5-symptom scoring method may increase diagnostic sensitivity in populations at high risk for bipolar disorder.

The genetics of bipolar disorder

Bipolar disorder (BD) is one of the most heritable mental illnesses, but the elucidation of its genetic basis has proven to be a very challenging endeavor. Genome-Wide Association Studies (GWAS) have transformed our understanding of BD, providing the first reproducible evidence of specific genetic markers and a highly polygenic architecture that overlaps with that of schizophrenia, major depression, and other disorders. Individual GWAS markers appear to confer little risk, but common variants together account for about 25% of the heritability of BD. A few higher-risk associations have also been identified, such as a rare copy number variant on chromosome 16p11.2. Large scale next-generation sequencing studies are actively searching for other alleles that confer substantial risk. As our understanding of the genetics of BD improves, there is growing optimism that some clear biological pathways will emerge, providing a basis for future studies aimed at molecular diagnosis and novel therapeutics.

Contribution of common and rare variants to bipolar disorder susceptibility in extended pedigrees from population isolates

Current evidence from case/control studies indicates that genetic risk for psychiatric disorders derives primarily from numerous common variants, each with a small phenotypic impact. The literature describing apparent segregation of bipolar disorder (BP) in numerous multigenerational pedigrees suggests that, in such families, large-effect inherited variants might play a greater role. To identify roles of rare and common variants on BP, we conducted genetic analyses in 26 Colombia and Costa Rica pedigrees ascertained for bipolar disorder 1 (BP1), the most severe and heritable form of BP. In these pedigrees, we performed microarray SNP genotyping of 838 individuals and high-coverage whole-genome sequencing of 449 individuals. We compared polygenic risk scores (PRS), estimated using the latest BP1 genome-wide association study (GWAS) summary statistics, between BP1 individuals and related controls. We also evaluated whether BP1 individuals had a higher burden of rare deleterious single-nucleotide variants (SNVs) and rare copy number variants (CNVs) in a set of genes related to BP1. We found that compared with unaffected relatives, BP1 individuals had higher PRS estimated from BP1 GWAS statistics (P = 0.001 ~ 0.007) and displayed modest increase in burdens of rare deleterious SNVs (P = 0.047) and rare CNVs (P = 0.002 ~ 0.033) in genes related to BP1. We did not observe rare variants segregating in the pedigrees. These results suggest that small-to-moderate effect rare and common variants are more likely to contribute to BP1 risk in these extended pedigrees than a few large-effect rare variants.

Whole-exome sequencing of 81 individuals from 27 multiply affected bipolar disorder families

Bipolar disorder (BD) is a highly heritable neuropsychiatric disease characterized by recurrent episodes of depression and mania. Research suggests that the cumulative impact of common alleles explains 25-38% of phenotypic variance, and that rare variants may contribute to BD susceptibility. To identify rare, high-penetrance susceptibility variants for BD, whole-exome sequencing (WES) was performed in three affected individuals from each of 27 multiply affected families from Spain and Germany. WES identified 378 rare, non-synonymous, and potentially functional variants. These spanned 368 genes, and were carried by all three affected members in at least one family. Eight of the 368 genes harbored rare variants that were implicated in at least two independent families. In an extended segregation analysis involving additional family members, five of these eight genes harbored variants showing full or nearly full cosegregation with BD. These included the brain-expressed genes RGS12 and NCKAP5, which were considered the most promising BD candidates on the basis of independent evidence. Gene enrichment analysis for all 368 genes revealed significant enrichment for four pathways, including genes reported in de novo studies of autism (p_adj < 0.006) and schizophrenia (p_adj = 0.015). These results suggest a possible genetic overlap with BD for autism and schizophrenia at the rare-sequence-variant level. The present study implicates novel candidate genes for BD development, and may contribute to an improved understanding of the biological basis of this common and often devastating disease.

A linkage and exome study implicates rare variants of KANK4 and CAP2 in bipolar disorder in a multiplex family

OBJECTIVES

Bipolar disorder (BD) is a neuropsychiatric disorder with a complex pattern of inheritance. Although many genetic studies have been conducted on BD, its genetic correlates remain uncertain. This study was aimed at  identifying the genetic underpinnings of the disorder in an Indian family, which has been under comprehensive clinical evaluation and follow-up for over 12 years.

METHODS

We analysed a four-generation family with several of its members diagnosed for BD employing a combination of genetic linkage and exome analysis.

RESULTS

We obtained suggestive LOD score for a chromosome 1 and a chromosome 6 marker (D1S410; LOD = 3.01, Ө = 0; and D6S289; LOD = 1.58, Ө = 0). Manual haplotyping of the regions encompassing these two markers helped delimit a critical genomic interval of 32.44 Mb (D1S2700-D1S435; chromosome 1p31.1-13.2) and another of 10.34 Mb (D6S470-D6S422; chromosome 6p22.3-22.2). We examined the exomic sequences corresponding to these two intervals and found rare variants, NM_181712.4: c.2461G>T (p.Asp821Tyr) in KANK4 at 1p31.1-13.2; and NM_006366:c.-93G>A, in the 5' UTR of CAP2 at 6p22.3-22.2.

CONCLUSIONS

Our studysuggests involvement of KANK4 or CAP2 or both in BD in this family. Further analysis of these two genes in BD patients and functional evaluation of the allelic variants identified are suggested.

Using linkage studies combined with whole-exome sequencing to identify novel candidate genes for familial colorectal cancer

Colorectal cancer (CRC) is a complex disorder for which the majority of the underlying germline predisposition factors remain still unidentified. Here, we combined whole‐exome sequencing (WES) and linkage analysis in families with multiple relatives affected by CRC to identify candidate genes harboring rare variants with potential high‐penetrance effects. Forty‐seven affected subjects from 18 extended CRC families underwent WES. Genome‐wide linkage analysis was performed under linear and exponential models. Suggestive linkage peaks were identified on chromosomes 1q22–q24.2 (maxSNP = rs2134095; LODlinear = 2.38, LODexp = 2.196), 7q31.2–q34 (maxSNP = rs6953296; LODlinear = 2.197, LODexp = 2.149) and 10q21.2–q23.1 (maxSNP = rs1904589; LODlinear = 1.445, LODexp = 2.195). These linkage signals were replicated in 10 independent sets of random markers from each of these regions. To assess the contribution of rare variants predicted to be pathogenic, we performed a family‐based segregation test with 89 rare variants predicted to be deleterious from 78 genes under the linkage intervals. This analysis showed significant segregation of rare variants with CRC in 18 genes (weighted p‐value > 0.0028). Protein network analysis and functional evaluation were used to suggest a plausible candidate gene for germline CRC predisposition. Etiologic rare variants implicated in cancer germline predisposition may be identified by combining traditional linkage with WES data. This approach can be used with already available NGS data from families with several sequenced members to further identify candidate genes involved germline predisposition to disease. This approach resulted in one candidate gene associated with increased risk of CRC but needs evidence from further studies.

What's new?

Inherited genetic factors are thought to account for more than one‐third of colorectal cancer (CRC) cases. Most predisposing genetic factors, however, remain unidentified. Here, genome‐wide linkage analysis using whole‐exome sequencing (WES) data was performed in families with marked CRC aggregation. The combined linkage‐sequencing approach identified possible linkage peaks on chromosomes 1q22‐q24.2, 7q31.2‐q34, and 10q21.2‐q23.1. Analyses of potentially pathogenic variants revealed significant segregation of rare variants in 18 genes, while functional analyses identified a plausible candidate gene for germline CRC predisposition. The findings underscore the utility of linkage analysis employing WES for the discovery of candidate genes for disease predisposition.

Clinical and genetic validity of quantitative bipolarity

Research has yet to provide a comprehensive understanding of the genetic basis of bipolar disorder (BP). In genetic studies, defining the phenotype by diagnosis may miss risk-allele carriers without BP. The authors aimed to test whether quantitatively detected subclinical symptoms of bipolarity identifies a heritable trait that infers risk for BP. The Quantitative Bipolarity Scale (QBS) was administered to 310 Old Order Amish or Mennonite individuals from multigenerational pedigrees; 110 individuals had psychiatric diagnoses (20 BP, 61 major depressive disorders (MDD), 3 psychotic disorders, 26 other psychiatric disorders). Familial aggregation of QBS was calculated using the variance components method to derive heritability and shared household effects. The QBS score was significantly higher in BP subjects (31.5 ± 3.6) compared to MDD (16.7 ± 2.0), other psychiatric diagnoses (7.0 ± 1.9), and no psychiatric diagnosis (6.0 ± 0.65) (all p < 0.001). QBS in the whole sample was significantly heritable (h² = 0.46 ± 0.15, p < 0.001) while the variance attributed to the shared household effect was not significant (p = 0.073). When subjects with psychiatric illness were removed, the QBS heritability was similar (h² = 0.59 ± 0.18, p < 0.001). These findings suggest that quantitative bipolarity as measured by QBS can separate BP from other psychiatric illnesses yet is significantly heritable with and without BP included in the pedigrees suggesting that the quantitative bipolarity describes a continuous heritable trait that is not driven by a discrete psychiatric diagnosis. Bipolarity trait assessment may be used to supplement the diagnosis of BP in future genetic studies and could be especially useful for capturing subclinical genetic contributions to a BP phenotype.

Current understanding of bipolar disorder: Toward integration of biological basis and treatment strategies

Biological studies of bipolar disorder initially focused on the mechanism of action for antidepressants and antipsychotic drugs, and the roles of monoamines (e.g., serotonin, dopamine) have been extensively studied. Thereafter, based on the mechanism of action of lithium, intracellular signal transduction systems, including inositol metabolism and intracellular calcium signaling, have drawn attention. Involvement of intracellular calcium signaling has been supported by genetics and cellular studies. Elucidation of the neural circuits affected by calcium signaling abnormalities is critical, and our previous study suggested a role of the paraventricular thalamic nucleus. The genetic vulnerability of mitochondria causes calcium dysregulation and results in the hyperexcitability of serotonergic neurons, which are suggested to be susceptible to oxidative stress. Efficacy of anticonvulsants, animal studies of candidate genes, and studies using induced pluripotent stem cell-derived neurons have suggested a relation between bipolar disorder and the hyperexcitability of neurons. Recent genetic findings suggest the roles of polyunsaturated acids. At the systems level, social rhythm therapy targets circadian rhythm abnormalities, and cognitive behavioral therapy may target emotion/cognition (E/C) imbalance. In the future, pharmacological and psychosocial treatments may be combined and optimized based on the biological basis of each patient, which will realize individualized treatment.

Health Needs Assessment of Five Pennsylvania Plain Populations

We performed a health needs assessment for five Plain communities in Pennsylvania from a random sample of households, comparing them to the general population of Pennsylvania adults. Plain respondents were more likely to drink well water, as likely to eat fruit and vegetables and much more likely to drink raw milk and be exposed to agricultural chemicals. Plain respondents were less likely to receive screening exams compared to the general population and there was variation from settlement to settlement in whether respondents had a regular doctor, whether they received preventive screenings or had their children vaccinated, with Mifflin County Amish generally lowest in these and Plain Mennonites highest. Plain respondents reported good physical and mental health compared to the general population but Groffdale Mennonite respondents had a high proportion of diagnoses of depression and were more likely to be receiving treatment for a mental health condition. Most Plain respondents would want a spouse tested for genetic disease with Mifflin County Amish least in favor of these tests. Despite their geographic and genetic isolation, the health of Plain communities in Pennsylvania is similar to that of other adults in the state.

Twin-based heritability of actimetry traits

There is a critical need for phenotypes with substantial heritability that can be used as endophenotypes in behavioral genetic studies. Activity monitoring, called actimetry, has potential as a means of assessing sleep and circadian rhythm traits that could serve as endophenotypes relevant to a range of psychopathologies. This study examined a range of actimetry traits for heritability using a classic twin design. The sample consisted of 195 subjects from 45 monozygotic (MZ) and 50 dizygotic (DZ) twin pairs aged 16-40 years. Subjects wore both a research-grade actimeter (GENEActiv) and a consumer-oriented device (FitBit) for 2 weeks. Sleep and circadian traits were extracted from GENEActiv data using PennZzz and ChronoSapiens software programs. Sleep statistics for a limited number of FitBit-collected traits were generated by its accompanying mobile app. Broad sense heritability was computed on a set of 33 MZ and 38 DZ twin pairs with complete data using both OpenMX and SOLAR software. These analyses yielded a large number of actimetry-derived traits, 20 of which showed high heritability (h² > 0.6), seven of which remain significant after Bonferroni correction. These results indicate that actimetry enables assessing a range of phenotypes with substantial heritability that may be useful as endophenotypes for genetic studies.

Whole-exome sequencing identifies variants associated with structural MRI markers in patients with bipolar disorders

BACKGROUND

Bipolar disorder (BD) is one of the most heritable psychiatric disorders. A growing number of whole-exome sequencing (WES) studies for BD has been performed, however, no research has examined the association between single nucleotide variants (SNVs) from WES and structural magnetic resonance imaging (MRI) data.

METHODS

We sequenced whole-exomes in 53 patients with BD and 82 healthy control participants at an initial discovery stage and investigated the impacts of SNVs in risk genes from WES analysis on the cortical gray-matter thickness and integrity of white matter tracts and in the following stage. Cortical thickness and white matter integrity were investigated using the FreeSurfer and TRACULA (Tracts Constrained by UnderLying Anatomy).

RESULTS

We identified 122 BD-related genes including KMT2C, AHNAK, CDH23, DCHS1, FRAS1, MACF1 and RYR3 and observed 27 recurrent copy number alteration regions including gain on 8p23.1 and loss on 15q11.1 - q11.2. Among them, single nucleotide polymorphism (SNP) rs4639425 in KMT2C gene, which regulates histone H3 lysine 4 (H3K4) methylation involved in chromatin remodeling, was associated with widespread alterations of white matter integrity including the cingulum, uncinate fasciculus, cortico-spinal tract, and superior longitudinal fasciculus.

LIMITATION

The small sample size of patients with BD in the genome data may cause our study to be underpowered when searching for putative rare mutations.

CONCLUSION

This study first combined a WES approach and neuroimaging findings in psychiatric disorders. We postulate the rs4639425 may be associated with BD-related microstructural changes of white matter tracts.

Rediscovering the value of families for psychiatric genetics research

As it is likely that both common and rare genetic variation are important for complex disease risk, studies that examine the full range of the allelic frequency distribution should be utilized to dissect the genetic influences on mental illness. The rate limiting factor for inferring an association between a variant and a phenotype is inevitably the total number of copies of the minor allele captured in the studied sample. For rare variation, with minor allele frequencies of 0.5% or less, very large samples of unrelated individuals are necessary to unambiguously associate a locus with an illness. Unfortunately, such large samples are often cost prohibitive. However, by using alternative analytic strategies and studying related individuals, particularly those from large multiplex families, it is possible to reduce the required sample size while maintaining statistical power. We contend that using whole genome sequence (WGS) in extended pedigrees provides a cost-effective strategy for psychiatric gene mapping that complements common variant approaches and WGS in unrelated individuals. This was our impetus for forming the "Pedigree-Based Whole Genome Sequencing of Affective and Psychotic Disorders" consortium. In this review, we provide a rationale for the use of WGS with pedigrees in modern psychiatric genetics research. We begin with a focused review of the current literature, followed by a short history of family-based research in psychiatry. Next, we describe several advantages of pedigrees for WGS research, including power estimates, methods for studying the environment, and endophenotypes. We conclude with a brief description of our consortium and its goals.

Exome sequencing in large, multiplex bipolar disorder families from Cuba

Bipolar disorder (BD) is a major psychiatric illness affecting around 1% of the global population. BD is characterized by recurrent manic and depressive episodes, and has an estimated heritability of around 70%. Research has identified the first BD susceptibility genes. However, the underlying pathways and regulatory networks remain largely unknown. Research suggests that the cumulative impact of common alleles with small effects explains only around 25-38% of the phenotypic variance for BD. A plausible hypothesis therefore is that rare, high penetrance variants may contribute to BD risk. The present study investigated the role of rare, nonsynonymous, and potentially functional variants via whole exome sequencing in 15 BD cases from two large, multiply affected families from Cuba. The high prevalence of BD in these pedigrees renders them promising in terms of the identification of genetic risk variants with large effect sizes. In addition, SNP array data were used to calculate polygenic risk scores for affected and unaffected family members. After correction for multiple testing, no significant increase in polygenic risk scores for common, BD-associated genetic variants was found in BD cases compared to healthy relatives. Exome sequencing identified a total of 17 rare and potentially damaging variants in 17 genes. The identified variants were shared by all investigated BD cases in the respective pedigree. The most promising variant was located in the gene SERPING1 (p.L349F), which has been reported previously as a genome-wide significant risk gene for schizophrenia. The present data suggest novel candidate genes for BD susceptibility, and may facilitate the discovery of disease-relevant pathways and regulatory networks.

Genetic pleiotropy between mood disorders, metabolic, and endocrine traits in a multigenerational pedigree

Bipolar disorder (BD) is a mental disorder characterized by alternating periods of depression and mania. Individuals with BD have higher levels of early mortality than the general population, and a substantial proportion of this is due to increased risk for comorbid diseases. To identify the molecular events that underlie BD and related medical comorbidities, we generated imputed whole-genome sequence data using a population-specific reference panel for an extended multigenerational Old Order Amish pedigree (n = 394), segregating BD and related disorders. First, we investigated all putative disease-causing variants at known Mendelian disease loci present in this pedigree. Second, we performed genomic profiling using polygenic risk scores (PRS) to establish each individual's risk for several complex diseases. We identified a set of Mendelian variants that co-occur in individuals with BD more frequently than their unaffected family members, including the R3527Q mutation in APOB associated with hypercholesterolemia. Using PRS, we demonstrated that BD individuals from this pedigree were enriched for the same common risk alleles for BD as the general population (β = 0.416, p = 6 × 10^-4). Furthermore, we find evidence for a common genetic etiology between BD risk and polygenic risk for clinical autoimmune thyroid disease (p = 1 × 10^-4), diabetes (p = 1 × 10^-3), and lipid traits such as triglyceride levels (p = 3 × 10^-4) in the pedigree. We identify genomic regions that contribute to the differences between BD individuals and unaffected family members by calculating local genetic risk for independent LD blocks. Our findings provide evidence for the extensive genetic pleiotropy that can drive epidemiological findings of comorbidities between diseases and other complex traits.

An examination of multiple classes of rare variants in extended families with bipolar disorder

Bipolar disorder (BD) is a complex psychiatric condition with high heritability, the genetic architecture of which likely comprises both common variants of small effect and rare variants of higher penetrance, the latter of which are largely unknown. Extended families with high density of illness provide an opportunity to map novel risk genes or consolidate evidence for existing candidates, by identifying genes carrying pathogenic rare variants. We performed whole-exome sequencing (WES) in 15 BD families (117 subjects, of whom 72 were affected), augmented with copy number variant (CNV) microarray data, to examine contributions of multiple classes of rare genetic variants within a familial context. Linkage analysis and haplotype reconstruction using WES-derived genotypes enabled exclusion of false-positive single-nucleotide variants (SNVs), CNV inheritance estimation, de novo variant identification and candidate gene prioritization. We found that rare predicted pathogenic variants shared among ≥3 affected relatives were overrepresented in postsynaptic density (PSD) genes (P = 0.002), with no enrichment in unaffected relatives. Genome-wide burden of likely gene-disruptive variants was no different in affected vs. unaffected relatives (P = 0.24), but correlated significantly with age of onset (P = 0.017), suggesting that a high disruptive variant burden may expedite symptom onset. The number of de novo variants was no different in affected vs. unaffected offspring (P = 0.89). We observed heterogeneity within and between families, with the most likely genetic model involving alleles of modest effect and reduced penetrance: a possible exception being a truncating X-linked mutation in IRS4 within a family-specific linkage peak. Genetic approaches combining WES, CNV and linkage analyses in extended families are promising strategies for gene discovery.

From the Psychiatrist's Couch to Induced Pluripotent Stem Cells: Bipolar Disease in a Dish

Bipolar disease (BD) is one of the major public health burdens worldwide and more people are affected every year. Comprehensive genetic studies have associated thousands of single nucleotide polymorphisms (SNPs) with BD risk; yet, very little is known about their functional roles. Induced pluripotent stem cells (iPSCs) are powerful tools for investigating the relationship between genotype and phenotype in disease-relevant tissues and cell types. Neural cells generated from BD-specific iPSCs are thought to capture associated genetic risk factors, known and unknown, and to allow the analysis of their effects on cellular and molecular phenotypes. Interestingly, an increasing number of studies on BD-derived iPSCs report distinct alterations in neural patterning, postmitotic calcium signaling, and neuronal excitability. Importantly, these alterations are partly normalized by lithium, a first line treatment in BD. In light of these exciting findings, we discuss current challenges to the field of iPSC-based disease modelling and future steps to be taken in order to fully exploit the potential of this approach for the investigation of BD and the development of new therapies.

Rare Risk Variants Identification by Identity-by-Descent Mapping and Whole-Exome Sequencing Implicates Neuronal Development Pathways in Schizophrenia and Bipolar Disorder

Schizophrenia (SCZ) and bipolar disorder (BPD) are highly heritable disorders with an estimated co-heritability of 68%. Hundreds of common alleles have been implicated, but recently a role for rare, high-penetrant variants has been also suggested in both disorders. This study investigated a familial cohort of SCZ and BPD patients from a closed population sample, where the high recurrence of the disorders and the homogenous genetic background indicate a possible enrichment in rare risk alleles. A total of 230 subjects (161 cases, 22 unaffected relatives, and 47 controls) were genetically investigated through an innovative strategy that integrates identity-by-descent (IBD) mapping and whole-exome sequencing (WES). IBD analysis allowed to track high-risk haplotypes (IBD_risk) shared exclusively by multiple patients from different families and possibly carrying the most penetrant alleles. A total of 444 non-synonymous sequence variants, of which 137 disruptive, were identified in IBD_risk haplotypes by WES. Interestingly, gene sets previously implicated in SCZ (i.e., post-synaptic density (PSD) proteins, voltage-gated calcium channels (VGCCs), and fragile X mental retardation protein (FMRP) targets) were found significantly enriched in genes carrying IBD_risk variants. Further, IBD_risk variants were preferentially affecting genes involved in the extracellular matrix (ECM) biology and axon guidance processes which appeared to be functionally connected in the pathway-derived meta-network analysis. Results thus confirm rare risk variants as key factors in SCZ and BPD pathogenesis and highlight a role for the development of neuronal connectivity in the etiology of both disorders.

Family studies to find rare high risk variants in migraine

INTRODUCTION

Migraine has long been known as a common complex disease caused by genetic and environmental factors. The pathophysiology and the specific genetic susceptibility are poorly understood. Common variants only explain a small part of the heritability of migraine. It is thought that rare genetic variants with bigger effect size may be involved in the disease. Since migraine has a tendency to cluster in families, a family approach might be the way to find these variants. This is also indicated by identification of migraine-associated loci in classical linkage-analyses in migraine families. A single migraine study using a candidate-gene approach was performed in 2010 identifying a rare mutation in the TRESK potassium channel segregating in a large family with migraine with aura, but this finding has later become questioned. The technologies of next-generation sequencing (NGS) now provides an affordable tool to investigate the genetic variation in the entire exome or genome. The family-based study design using NGS is described in this paper. We also review family studies using NGS that have been successful in finding rare variants in other common complex diseases in order to argue the promising application of a family approach to migraine.

METHOD

PubMed was searched to find studies that looked for rare genetic variants in common complex diseases through a family-based design using NGS, excluding studies looking for de-novo mutations, or using a candidate-gene approach and studies on cancer. All issues from Nature Genetics and PLOS genetics 2014, 2015 and 2016 (UTAI June) were screened for relevant papers. Reference lists from included and other relevant papers were also searched. For the description of the family-based study design using NGS an in-house protocol was used.

RESULTS

Thirty-two successful studies, which covered 16 different common complex diseases, were included in this paper. We also found a single migraine study. Twenty-three studies found one or a few family specific variants (less than five), while other studies found several possible variants. Not all of them were genome wide significant. Four studies performed follow-up analyses in unrelated cases and controls and calculated odds ratios that supported an association between detected variants and risk of disease. Studies of 11 diseases identified rare variants that segregated fully or to a large degree with the disease in the pedigrees.

CONCLUSION

It is possible to find rare high risk variants for common complex diseases through a family-based approach. One study using a family approach and NGS to find rare variants in migraine has already been published but with strong limitations. More studies are under way.

The Sardinian Puzzle: Concentration of Major Psychoses and Suicide in the Same Sub-Regions Across One Century

BACKGROUND

Sardinia, the second largest Mediterranean island has long been considered a privileged observatory for the study of several medical conditions. The peculiar epidemiology of mood disorders and suicide across Sardinian sub-regions has long intrigued clinicians and researchers.

OBJECTIVE

The principal aim of the present study was to test whether the geographical distribution of suicides committed in Sardinian over the last three decades are comparable with the geographical origin of patients hospitalized up to half a century ago.

METHOD

The distribution of the municipalities of origin of the patients hospitalized in Sardinia between 1901 and 1964 for schizophrenia, bipolar disorder, and depression was reanalyzed and compared with the distribution of municipalities where suicides were committed between 1980 and 2013. Data were also analyzed by the altitude above the sea level and by the population size of the municipalities.

RESULTS

There was a significant variation of hospitalization and suicide rates across Sardinian sub-regions. The sub-regions of origin of the patients hospitalized for schizophrenia and bipolar disorder correlated with each other (P = 0.047). Both hospitalizations and suicides were more incident in municipalities with a higher altitude and a smaller population size. The incidence of hospitalizations and suicides correlated significantly with each other both at the municipality (P = 1.86 x 10-7) and at the sub-region level (P = 1.71 x 10-7).

CONCLUSION

The present study confirms the peculiar geographical distribution of major psychoses and suicide in Sardinia. The two phenomena appear to have been correlated for as long as one century.

A novel relationship for schizophrenia, bipolar and major depressive disorder Part 3: Evidence from chromosome 3 high density association screen

Familial clustering of schizophrenia (SCZ), bipolar disorder (BPD), and major depressive disorder (MDD) was systematically reported (Aukes et al, Genet Med 2012, 14, 338-341) and convergent evidence from genetics, symptomatology, and psychopharmacology imply that there are intrinsic connections between these three major psychiatric disorders, for example, any two or even three of these disorders could co-exist in some families. A total of 60, 838 single-nucleotide polymorphisms (SNPs) on chromosome 3 were genotyped by Affymetrix Genome-Wide Human SNP array 6.0 on 119 SCZ, 253 BPD (type-I), 177 MDD patients and 1,000 controls. The population of Shandong province was formed in 14 century and believed that it belongs to homogenous population. Associated SNPs were systematically revealed and outstanding susceptibility genes (CADPS, GRM7,KALRN, LSAMP, NLGN1, PRICKLE2, ROBO2) were identified. Unexpectedly, flanking genes for the associated SNPs distinctive for BPD and/or MDD were replicated in an enlarged cohort of 986 SCZ patients. The evidence from this chromosome 3 analysis supports the notion that both of bipolar and MDD might be subtypes of schizophrenia rather than independent disease entity. Also, a similar finding was detected on chromosome 5, 6, 7, and 8 (Chen et al. Am J Transl Res 2017;9 (5):2473-2491; Curr Mol Med 2016;16(9):840-854; Behav Brain Res 2015;293:241-251; Mol Neurobiol 2016. doi: 10.1007/s12035-016-0102-1). Furthermore, PRICKLE2 play an important role in the pathogenesis of three major psychoses in this population.

A population-specific reference panel empowers genetic studies of Anabaptist populations

Genotype imputation is a powerful strategy for achieving the large sample sizes required for identification of variants underlying complex phenotypes, but imputation of rare variants remains problematic. Genetically isolated populations offer one solution, however population-specific reference panels are needed to assure optimal imputation accuracy and allele frequency estimation. Here we report the Anabaptist Genome Reference Panel (AGRP), the first whole-genome catalogue of variants and phased haplotypes in people of Amish and Mennonite ancestry. Based on high-depth whole-genome sequence (WGS) from 265 individuals, the AGRP contains >12 M high-confidence single nucleotide variants and short indels, of which ~12.5% are novel. These Anabaptist-specific variants were more deleterious than variants with comparable frequencies observed in the 1000 Genomes panel. About 43,000 variants showed enriched allele frequencies in AGRP, consistent with drift. When combined with the 1000 Genomes Project reference panel, the AGRP substantially improved imputation, especially for rarer variants. The AGRP is freely available to researchers through an imputation server.

Gene-wide Association Study Reveals RNF122 Ubiquitin Ligase as a Novel Susceptibility Gene for Attention Deficit Hyperactivity Disorder

Attention Deficit Hyperactivity Disorder (ADHD) is a common childhood-onset neurodevelopmental condition characterized by pervasive impairment of attention, hyperactivity, and/or impulsivity that can persist into adulthood. The aetiology of ADHD is complex and multifactorial and, despite the wealth of evidence for its high heritability, genetic studies have provided modest evidence for the involvement of specific genes and have failed to identify consistent and replicable results. Due to the lack of robust findings, we performed gene-wide and pathway enrichment analyses using pre-existing GWAS data from 607 persistent ADHD subjects and 584 controls, produced by our group. Subsequently, expression profiles of genes surpassing a follow-up threshold of P-value < 1e-03 in the gene-wide analyses were tested in peripheral blood mononucleated cells (PBMCs) of 45 medication-naive adults with ADHD and 39 healthy unrelated controls. We found preliminary evidence for genetic association between RNF122 and ADHD and for its overexpression in adults with ADHD. RNF122 encodes for an E3 ubiquitin ligase involved in the proteasome-mediated processing, trafficking, and degradation of proteins that acts as an essential mediator of the substrate specificity of ubiquitin ligation. Thus, our findings support previous data that place the ubiquitin-proteasome system as a promising candidate for its involvement in the aetiology of ADHD.

Molecular circadian rhythm shift due to bright light exposure before bedtime is related to subthreshold bipolarity

This study examined the link between circadian rhythm changes due to bright light exposure and subthreshold bipolarity. Molecular circadian rhythms, polysomnography, and actigraphy data were studied in 25 young, healthy male subjects, divided into high and low mood disorder questionnaire (MDQ) score groups. During the first 2 days of the study, the subjects were exposed to daily-living light (150 lux) for 4 hours before bedtime. Saliva and buccal cells were collected 5 times a day for 2 consecutive days. During the subsequent 5 days, the subjects were exposed to bright light (1,000 lux), and saliva and buccal cell samples were collected in the same way. Molecular circadian rhythms were analyzed using sine regression. Circadian rhythms of cortisol (F = 16.956, p < 0.001) and relative PER1/ARNTL gene expression (F = 122.1, p < 0.001) showed a delayed acrophase in both groups after bright light exposure. The high MDQ score group showed a significant delay in acrophase compared to the low MDQ score group only in salivary cortisol (F = 8.528, p = 0.008). The high MDQ score group showed hypersensitivity in cortisol rhythm shift after bright light exposure, suggesting characteristic molecular circadian rhythm changes in the high MDQ score group may be related to biological processes downstream from core circadian clock gene expression.

Commonalities in Development of Pure Breeds and Population Isolates Revealed in the Genome of the Sardinian Fonni's Dog

The island inhabitants of Sardinia have long been a focus for studies of complex human traits due to their unique ancestral background and population isolation reflecting geographic and cultural restriction. Population isolates share decreased genomic diversity, increased linkage disequilibrium, and increased inbreeding coefficients. In many regions, dogs and humans have been exposed to the same natural and artificial forces of environment, growth, and migration. Distinct dog breeds have arisen through human-driven selection of characteristics to meet an ideal standard of appearance and function. The Fonni's Dog, an endemic dog population on Sardinia, has not been subjected to an intensive system of artificial selection, but rather has developed alongside the human population of Sardinia, influenced by geographic isolation and unregulated selection based on its environmental adaptation and aptitude for owner-desired behaviors. Through analysis of 28 dog breeds, represented with whole-genome sequences from 13 dogs and ∼170,000 genome-wide single nucleotide variants from 155 dogs, we have produced a genomic illustration of the Fonni's Dog. Genomic patterns confirm within-breed similarity, while population and demographic analyses provide spatial identity of Fonni's Dog to other Mediterranean breeds. Investigation of admixture and fixation indices reveals insights into the involvement of Fonni's Dogs in breed development throughout the Mediterranean. We describe how characteristics of population isolates are reflected in dog breeds that have undergone artificial selection, and are mirrored in the Fonni's Dog through traditional isolating factors that affect human populations. Lastly, we show that the genetic history of Fonni's Dog parallels demographic events in local human populations.

Health Needs Assessment of Plain Populations in Lancaster County, Pennsylvania

We performed a health needs assessment for three Plain communities in Lancaster County, Pennsylvania from a random sample of households. Compared with the general population of adults, Plain respondents were more likely to be married, to have children, and they had large families; they were more likely to drink well water, to eat fruit and vegetables, to drink raw milk, and to live on a farm. Plain respondents had better physical and mental health and were less likely to have been diagnosed with various medical conditions compared with the general population of adults in Lancaster County but Old Order Mennonite respondents were more likely to have been diagnosed compared with Old Order Amish respondents. Plain respondents usually have a regular doctor and often receive preventive care but Old Order Mennonite respondents were more likely to have a regular doctor, to receive preventive care, to have had their children vaccinated, and to receive routine dental care compared with Old Order Amish respondents. Despite their relative geographic and genetic isolation, and despite the small, relative differences noted, the health of Plain communities in Lancaster County is similar to that of other adults in the County.

Exome Sequencing of Familial Bipolar Disorder

IMPORTANCE

Complex disorders, such as bipolar disorder (BD), likely result from the influence of both common and rare susceptibility alleles. While common variation has been widely studied, rare variant discovery has only recently become feasible with next-generation sequencing.

OBJECTIVE

To utilize a combined family-based and case-control approach to exome sequencing in BD using multiplex families as an initial discovery strategy, followed by association testing in a large case-control meta-analysis.

DESIGN, SETTING, AND PARTICIPANTS

We performed exome sequencing of 36 affected members with BD from 8 multiplex families and tested rare, segregating variants in 3 independent case-control samples consisting of 3541 BD cases and 4774 controls.

MAIN OUTCOMES AND MEASURES

We used penalized logistic regression and 1-sided gene-burden analyses to test for association of rare, segregating damaging variants with BD. Permutation-based analyses were performed to test for overall enrichment with previously identified gene sets.

RESULTS

We found 84 rare (frequency <1%), segregating variants that were bioinformatically predicted to be damaging. These variants were found in 82 genes that were enriched for gene sets previously identified in de novo studies of autism (19 observed vs. 10.9 expected, P = .0066) and schizophrenia (11 observed vs. 5.1 expected, P = .0062) and for targets of the fragile X mental retardation protein (FMRP) pathway (10 observed vs. 4.4 expected, P = .0076). The case-control meta-analyses yielded 19 genes that were nominally associated with BD based either on individual variants or a gene-burden approach. Although no gene was individually significant after correction for multiple testing, this group of genes continued to show evidence for significant enrichment of de novo autism genes (6 observed vs 2.6 expected, P = .028).

CONCLUSIONS AND RELEVANCE

Our results are consistent with the presence of prominent locus and allelic heterogeneity in BD and suggest that very large samples will be required to definitively identify individual rare variants or genes conferring risk for this disorder. However, we also identify significant associations with gene sets composed of previously discovered de novo variants in autism and schizophrenia, as well as targets of the FRMP pathway, providing preliminary support for the overlap of potential autism and schizophrenia risk genes with rare, segregating variants in families with BD.

High-throughput sequencing of the synaptome in major depressive disorder

Major depressive disorder (MDD) is among the leading causes of worldwide disability. Despite its significant heritability, large-scale genome-wide association studies (GWASs) of MDD have yet to identify robustly associated common variants. Although increased sample sizes are being amassed for the next wave of GWAS, few studies have as yet focused on rare genetic variants in the study of MDD. We sequenced the exons of 1742 synaptic genes previously identified by proteomic experiments. PLINK/SEQ was used to perform single variant, gene burden and gene set analyses. The GeneMANIA interaction database was used to identify protein-protein interaction-based networks. Cases were selected from a familial collection of early-onset, recurrent depression and were compared with screened controls. After extensive quality control, we analyzed 259 cases with familial, early-onset MDD and 334 controls. The distribution of association test statistics for the single variant and gene burden analyses were consistent with the null hypothesis. However, analysis of prioritized gene sets showed a significant association with damaging singleton variants in a Cav2-adaptor gene set (odds ratio=2.6; P=0.0008) that survived correction for all gene sets and annotation categories tested (empirical P=0.049). In addition, we also found statistically significant evidence for an enrichment of rare variants in a protein-based network of 14 genes involved in actin polymerization and dendritic spine formation (nominal P=0.0031). In conclusion, we have identified a statistically significant gene set and gene network of rare variants that are over-represented in MDD, providing initial evidence that calcium signaling and dendrite regulation may be involved in the etiology of depression.

Genetics of Bipolar Disorder: Recent Update and Future Directions

Although genetic studies of Bipolar Disorder have been pursued for decades, it has only been in the last several years that clearly replicated findings have emerged. These findings, typically of modest effects, point to a polygenic genetic architecture consisting of multiple common and rare susceptibility variants. While larger genome-wide association studies are ongoing, the advent of whole exome and genome sequencing should lead to the identification of rare, and potentially more penetrant, variants. Progress along both fronts will provide novel insights into the biology of Bipolar Disorder and help usher in a new era of personalized medicine and improved treatments.

Transcriptomic Analysis of Induced Pluripotent Stem Cells Derived from Patients with Bipolar Disorder from an Old Order Amish Pedigree

Fibroblasts from patients with Type I bipolar disorder (BPD) and their unaffected siblings were obtained from an Old Order Amish pedigree with a high incidence of BPD and reprogrammed to induced pluripotent stem cells (iPSCs). Established iPSCs were subsequently differentiated into neuroprogenitors (NPs) and then to neurons. Transcriptomic microarray analysis was conducted on RNA samples from iPSCs, NPs and neurons matured in culture for either 2 weeks (termed early neurons, E) or 4 weeks (termed late neurons, L). Global RNA profiling indicated that BPD and control iPSCs differentiated into NPs and neurons at a similar rate, enabling studies of differentially expressed genes in neurons from controls and BPD cases. Significant disease-associated differences in gene expression were observed only in L neurons. Specifically, 328 genes were differentially expressed between BPD and control L neurons including GAD1, glutamate decarboxylase 1 (2.5 fold) and SCN4B, the voltage gated type IV sodium channel beta subunit (-14.6 fold). Quantitative RT-PCR confirmed the up-regulation of GAD1 in BPD compared to control L neurons. Gene Ontology, GeneGo and Ingenuity Pathway Analysis of differentially regulated genes in L neurons suggest that alterations in RNA biosynthesis and metabolism, protein trafficking as well as receptor signaling pathways may play an important role in the pathophysiology of BPD.

Disruption of sonic hedgehog signaling in Ellis-van Creveld dwarfism confers protection against bipolar affective disorder

Ellis-van Creveld syndrome, an autosomal recessively inherited chondrodysplastic dwarfism, is frequent among Old Order Amish of Pennsylvania. Decades of longitudinal research on bipolar affective disorder (BPAD) revealed cosegregation of high numbers of EvC and Bipolar I (BPI) cases in several large Amish families descending from the same pioneer. Despite the high prevalence of both disorders in these families, no EvC individual has ever been reported with BPI. The proximity of the EVC gene to our previously reported chromosome 4p16 BPAD locus with protective alleles, coupled with detailed clinical observations that EvC and BPI do not occur in the same individuals, led us to hypothesize that the genetic defect causing EvC in the Amish confers protection from BPI. This hypothesis is supported by a significant negative association of these two disorders when contrasted with absence of disease (P=0.029, Fisher's exact test, two-sided, verified by permutation to estimate the null distribution of the test statistic). As homozygous Amish EVC mutations causing EvC dwarfism do so by disrupting sonic hedgehog (Shh) signaling, our data implicate Shh signaling in the underlying pathophysiology of BPAD. Understanding how disrupted Shh signaling protects against BPI could uncover variants in the Shh pathway that cause or increase risk for this and related mood disorders.

On the Teneurin track: a new synaptic organization molecule emerges

To achieve proper synaptic development and function, coordinated signals must pass between the pre- and postsynaptic membranes. Such transsynaptic signals can be comprised of receptors and secreted ligands, membrane associated receptors, and also pairs of synaptic cell adhesion molecules. A critical open question bridging neuroscience, developmental biology, and cell biology involves identifying those signals and elucidating how they function. Recent work in Drosophila and vertebrate systems has implicated a family of proteins, the Teneurins, as a new transsynaptic signal in both the peripheral and central nervous systems. The Teneurins have established roles in neuronal wiring, but studies now show their involvement in regulating synaptic connections between neurons and bridging the synaptic membrane and the cytoskeleton. This review will examine the Teneurins as synaptic cell adhesion molecules, explore how they regulate synaptic organization, and consider how some consequences of human Teneurin mutations may have synaptopathic origins.

Copy number variants encompassing Mendelian disease genes in a large multigenerational family segregating bipolar disorder

BACKGROUND

Bipolar affective disorder (BP) is a common, highly heritable psychiatric disorder characterized by periods of depression and mania. Using dense SNP genotype data, we characterized CNVs in 388 members of an Old Order Amish Pedigree with bipolar disorder. We identified CNV regions arising from common ancestral mutations by utilizing the pedigree information. By combining this analysis with whole genome sequence data in the same individuals, we also explored the role of compound heterozygosity.

RESULTS

Here we describe 541 inherited CNV regions, of which 268 are rare in a control population of European origin but present in a large number of Amish individuals. In addition, we highlight a set of CNVs found at higher frequencies in BP individuals, and within genes known to play a role in human development and disease. As in prior reports, we find no evidence for an increased burden of CNVs in BP individuals, but we report a trend towards a higher burden of CNVs in known Mendelian disease loci in bipolar individuals (BPI and BPII, p = 0.06).

CONCLUSIONS

We conclude that CNVs may be contributing factors in the phenotypic presentation of mood disorders and co-morbid medical conditions in this family. These results reinforce the hypothesis of a complex genetic architecture underlying BP disorder, and suggest that the role of CNVs should continue to be investigated in BP data sets.

Variants in Ion Channel Genes Link Phenotypic Features of Bipolar Illness to Specific Neurobiological Process Domains

Recent advances in genome-wide association studies are pointing towards a major role for voltage-gated ion channels in neuropsychiatric disorders and, in particular, bipolar disorder (BD). The phenotype of BD is complex, with symptoms during mood episodes and deficits persisting between episodes. We have tried to elucidate the common neurobiological mechanisms associated with ion channel signaling in order to provide a new perspective on the clinical symptoms and possible endophenotypes seen in BD patients. We propose a model in which the multiple variants in genes coding for ion channel proteins would perturb motivational circuits, synaptic plasticity, myelination, hypothalamic-pituitary-adrenal axis function, circadian neuronal rhythms, and energy regulation. These changes in neurobiological mechanisms would manifest in endophenotypes of aberrant reward processing, white matter hyperintensities, deficits in executive function, altered frontolimbic connectivity, increased amygdala activity, increased melatonin suppression, decreased REM latency, and aberrant myo-inositol/ATP shuttling. The endophenotypes result in behaviors of poor impulse control, motivational changes, cognitive deficits, abnormal stress response, sleep disturbances, and energy changes involving different neurobiological process domains. The hypothesis is that these disturbances start with altered neural circuitry during development, following which multiple environmental triggers may disrupt the neuronal excitability balance through an activity-dependent molecular process, resulting in clinical mood episodes.