A population-based study of KCNH7 p.Arg394His and bipolar spectrum disorder

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.

Adam, amigo, brain, and K channel

Voltage-dependent K+ (Kv) channels are diverse, comprising the classical Shab - Kv2, Shaker - Kv1, Shal - Kv4, and Shaw - Kv3 families. The Shaker family alone consists of Kv1.1, Kv1.2, Kv1.3, Kv1.4, Kv1.5, Kv1.6, and Kv1.7. Moreover, the Shab family comprises two functional (Kv2.1 and Kv2.2) and several "silent" alpha subunits (Kv2.3, Kv5, Kv6, Kv8, and Kv9), which do not generate K current. However, e.g., Kv8.1, via heteromerization, inhibits outward currents of the same family or even that of Shaw. This property of Kv8.1 is similar to those of designated beta subunits or non-selective auxiliary elements, including ADAM or AMIGO proteins. Kv channels and, in turn, ADAM may modulate the synaptic long-term potentiation (LTP). Prevailingly, Kv1.1 and Kv1.5 are attributed to respective brain and heart pathologies, some of which may occur simultaneously. The aforementioned channel proteins are apparently involved in several brain pathologies, including schizophrenia and seizures.

Molecular signature of excessive female aggression: study of stressed mice with genetic inactivation of neuronal serotonin synthesis

Aggression is a complex social behavior, critically involving brain serotonin (5-HT) function. The neurobiology of female aggression remains elusive, while the incidence of its manifestations has been increasing. Yet, animal models of female aggression are scarce. We previously proposed a paradigm of female aggression in the context of gene x environment interaction where mice with partial genetic inactivation of tryptophan hydroxylase-2 (Tph2+/- mice), a key enzyme of neuronal 5-HT synthesis, are subjected to predation stress resulting in pathological aggression. Using deep sequencing and the EBSeq method, we studied the transcriptomic signature of excessive aggression in the prefrontal cortex of female Tph2+/- mice subjected to rat exposure stress and food deprivation. Challenged mutants, but not other groups, displayed marked aggressive behaviors. We found 26 genes with altered expression in the opposite direction between stressed groups of both Tph2 genotypes. We identified several molecular markers, including Dgkh, Arfgef3, Kcnh7, Grin2a, Tenm1 and Epha6, implicated in neurodevelopmental deficits and psychiatric conditions featuring impaired cognition and emotional dysregulation. Moreover, while 17 regulons, including several relevant to neural plasticity and function, were significantly altered in stressed mutants, no alteration in regulons was detected in stressed wildtype mice. An interplay of the uncovered pathways likely mediates partial Tph2 inactivation in interaction with severe stress experience, thus resulting in excessive female aggression.

Locus Coeruleus Neurons' Firing Pattern Is Regulated by ERG Voltage-Gated K+ Channels

Locus coeruleus (LC) neurons, with their extensive innervations throughout the brain, control a broad range of physiological processes. Several ion channels have been characterized in LC neurons that control intrinsic membrane properties and excitability. However, ERG (ether-à-go-go-related gene) K⁺ channels that are particularly important in setting neuronal firing rhythms and automaticity have not as yet been discovered in the LC. Moreover, the neurophysiological and pathophysiological roles of ERG channels in the brain remain unclear despite their expression in several structures. By performing immunohistochemical investigations, we found that ERG-1A, ERG-1B, ERG-2 and ERG-3 are highly expressed in the LC neurons of mice. To examine the functional role of ERG channels, current-clamp recordings were performed on mouse LC neurons in brain slices under visual control. ERG channel blockade by WAY-123,398, a class III anti-arrhythmic agent, increased the spontaneous firing activity and discharge irregularity of LC neurons. Here, we have shown the presence of distinct ERG channel subunits in the LC which play an imperative role in modulating neuronal discharge patterns. Thus, we propose that ERG channels are important players behind the changes in, and/or maintenance of, LC firing patterns that are implicated in the generation of different behaviors and in several disorders.

Integration of high-resolution promoter profiling assays reveals novel, cell type-specific transcription start sites across 115 human cell and tissue types

Accurate transcription start site (TSS) annotations are essential for understanding transcriptional regulation and its role in human disease. Gene collections such as GENCODE contain annotations for tens of thousands of TSSs, but not all of these annotations are experimentally validated nor do they contain information on cell type–specific usage. Therefore, we sought to generate a collection of experimentally validated TSSs by integrating RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression (RAMPAGE) data from 115 cell and tissue types, which resulted in a collection of approximately 50 thousand representative RAMPAGE peaks. These peaks are primarily proximal to GENCODE-annotated TSSs and are concordant with other transcription assays. Because RAMPAGE uses paired-end reads, we were then able to connect peaks to transcripts by analyzing the genomic positions of the 3′ ends of read mates. Using this paired-end information, we classified the vast majority (37 thousand) of our RAMPAGE peaks as verified TSSs, updating TSS annotations for 20% of GENCODE genes. We also found that these updated TSS annotations are supported by epigenomic and other transcriptomic data sets. To show the utility of this RAMPAGE rPeak collection, we intersected it with the NHGRI/EBI genome-wide association study (GWAS) catalog and identified new candidate GWAS genes. Overall, our work shows the importance of integrating experimental data to further refine TSS annotations and provides a valuable resource for the biological community.

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.

Kv11 (ether-à-go-go-related gene) voltage-dependent K+ channels promote resonance and oscillation of subthreshold membrane potentials

Key points

Some ion channels are known to behave as inductors and make up the parallel resonant circuit in the plasma membrane of neurons, which enables neurons to respond to current inputs with a specific frequency (so‐called ‘resonant properties’).

Here, we report that heterologous expression of mouse Kv11 voltage‐dependent K⁺ channels generate resonance and oscillation at depolarized membrane potentials in HEK293 cells; expressions of individual Kv11 subtypes generate resonance and oscillation with different frequency properties.

Kv11.3‐expressing HEK293 cells exhibited transient conductance changes that opposed the current changes induced by voltage steps; this probably enables Kv11 channels to behave like an inductor.

The resonance and oscillation of inferior olivary neurons were impaired at the resting membrane potential in Kv11.3 knockout mice.

This study helps to elucidate basic ion channel properties that are crucial for the frequency responses of neurons.

Abstract

The plasma membranes of some neurons preferentially respond to current inputs with a specific frequency, and output as large voltage changes. This property is called resonance, and is thought to be mediated by ion channels that show inductor‐like behaviour. However, details of the candidate ion channels remain unclear. In this study, we mainly focused on the functional roles of Kv11 potassium (K⁺) channels, encoded by ether‐á‐go‐go‐related genes, in resonance in mouse inferior olivary (IO) neurons. We transfected HEK293 cells with long or short splice variants of Kv11.1 (Merg1a and Merg1b) or Kv11.3, and examined membrane properties using whole‐cell recording. Transfection with Kv11 channels reproduced resonance at membrane potentials depolarized from the resting state. Frequency ranges of Kv11.3‐, Kv11.1(Merg1b)‐ and Kv11.1(Merg1a)‐expressing cells were 2–6 Hz, 2–4 Hz, and 0.6–0.8 Hz, respectively. Responses of Kv11.3 currents to step voltage changes were essentially similar to those of inductor currents in the resistor–inductor–capacitor circuit. Furthermore, Kv11 transfections generated membrane potential oscillations. We also confirmed the contribution of HCN1 channels as a major mediator of resonance at more hyperpolarized potentials by transfection into HEK293 cells. The Kv11 current kinetics and properties of Kv11‐dependent resonance suggested that Kv11.3 mediated resonance in IO neurons. This finding was confirmed by the impairment of resonance and oscillation at –30 to –60 mV in Kcnh7 (Kv11.3) knockout mice. These results suggest that Kv11 channels have important roles in inducing frequency‐dependent responses in a subtype‐dependent manner from resting to depolarized membrane potentials.

Some ion channels are known to behave as inductors and make up the parallel resonant circuit in the plasma membrane of neurons, which enables neurons to respond to current inputs with a specific frequency (so‐called ‘resonant properties’).

Here, we report that heterologous expression of mouse Kv11 voltage‐dependent K⁺ channels generate resonance and oscillation at depolarized membrane potentials in HEK293 cells; expressions of individual Kv11 subtypes generate resonance and oscillation with different frequency properties.

Kv11.3‐expressing HEK293 cells exhibited transient conductance changes that opposed the current changes induced by voltage steps; this probably enables Kv11 channels to behave like an inductor.

The resonance and oscillation of inferior olivary neurons were impaired at the resting membrane potential in Kv11.3 knockout mice.

This study helps to elucidate basic ion channel properties that are crucial for the frequency responses of neurons.

Antidepressive effect of an inward rectifier K+ channel blocker peptide, tertiapin-RQ

Renal outer medullary K⁺ channel, ROMK (Kir1.1, kcnj1) is expressed in the kidney and brain, but its role in the central nervous system remains unknown. Recent studies suggested an involvement of the ROMK channel in mental diseases. Tertiapin (TPN) is a European honey bee venom peptide and is reported to selectively block the ROMK channel. Here, we have chemically synthesized a series of mutated TPN peptides, including TPN-I8R and -M13Q (TPN-RQ), reported previously, and examined their blocking activity on the ROMK channel. Among 71 peptides tested, TPN-RQ was found to block the ROMK channel most effectively. Whole-cell patch-clamp recordings showed the essential roles of two disulfide bonds and the circular structure for the blockade activity. To examine the central role, we injected TPN-RQ intracerebroventricularly and examined the effects on depression- and anxiety-like behaviors in mice. TPN-RQ showed an antidepressive effect in tail-suspension and forced swim tests. The injection of TPN-RQ also enhanced the anxiety-like behavior in the elevated plus-maze and light/dark box tests and impaired spontaneous motor activities in balance beam and wheel running tests. Administration of TPM-RQ suppressed the anti-c-Fos immunoreactivity in the lateral septum, without affecting immunoreactivity in antidepressant-related nuclei, e.g. the dorsal raphe nucleus and locus coeruleus. TPN-RQ may exert its antidepressive effects through a different mechanism from current drugs.

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.

Kv2.1 voltage-gated potassium channels in developmental perspective

Kv2.1 voltage-gated potassium channels consist of two types of α-subunits: (a) electrically-active Kcnb1 α-subunits and (b) silent or modulatory α-subunits plus β-subunits that, similar to silent α-subunits, also regulate electrically-active subunits. Voltage-gated potassium channels were traditionally viewed, mainly by electrophysiologists, as regulators of the electrical activity of the plasma membrane in excitable cells, a role that is performed by transmembrane protein domains of α-subunits that form the electric pore. Genetic studies revealed a role for this region of α-subunits of voltage-gated potassium channels in human neurodevelopmental disorders, such as epileptic encephalopathy. The N- and C-terminal domains of α-subunits interact to form the cytoplasmic subunit of heterotetrameric potassium channels that regulate electric pores. Subsequent animal studies revealed the developmental functions of Kcnb1-containing voltage-gated potassium channels and illustrated their role during brain development and reproduction. These functions of potassium channels are discussed in this review in the context of regulatory interactions between electrically-active and regulatory subunits.

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.

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.

Tale of tail current

The largest biomass of channel proteins is located in unicellular organisms and bacteria that have no organs. However, orchestrated bidirectional ionic currents across the cell membrane via the channels are important for the functioning of organs of organisms, and equally concern both fauna or flora. Several ion channels are activated in the course of action potentials. One of the hallmarks of voltage-dependent channels is a 'tail current' - deactivation as observed after prior and sufficient activation predominantly at more depolarized potentials e.g. for Kv while upon hyperpolarization for HCN α subunits. Tail current also reflects the timing of channel closure that is initiated upon termination of stimuli. Finally, deactivation of currents during repolarization could be a selective estimate for given channel as in case of HERG, if dedicated long and more depolarized 'tail pulse' is used. Since from a holding potential of e.g. -70 mV are often a family of outward K⁺ currents comprising I_A and I_K are simultaneously activated in native cells.

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.

Association Study of KCNH7 Polymorphisms and Individual Responses to Risperidone Treatment in Schizophrenia

Risperidone has been used to treat the symptoms of schizophrenia and to reduce its relapse. However, the responses to treatment show great variability among patients. The potassium channel has been reported as an effective target for antipsychotics. KCNH7, a member of the voltage-gated K⁺ channel Kv11 family, is primarily expressed in the brain. Here, we assessed the genetic association of KCNH7 with risperidone responses in 393 schizophrenia patients. The patients were treated with risperidone for 6 weeks. The reduction rates of Positive and Negative Syndrome Scale (PANSS) scores were determined to quantify drug response. We also examined the associations between six single-nucleotide polymorphisms (SNPs) of KCNH7 and the risperidone responses for a total of 6 weeks. The SNP rs77699177 (C > T) in the KCNH7 gene intron was significantly associated with the treatment response reflected by the PANSS reduction rate (CC, 55.8 ± 23.0; TC, 70.9 ± 20.3, P = 0.000110), indicating that patients with the TC genotype have better efficacy for antipsychotic therapy. The rs2241240 SNP also showed a significant association with treatment responses after 6 weeks of treatment (P = 0.00256). The findings indicate that the voltage-gated K⁺ channel KCNH7 is a potential functional marker for the identification of the response to risperidone treatment in schizophrenia patients. Note: The study was registered under clinical trial number ChiCTR-RNC-09000522 (http://www.chictr.org/).

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.

Mapping of breakpoints in balanced chromosomal translocations by shallow whole-genome sequencing points to EFNA5, BAHD1 and PPP2R5E as novel candidates for genes causing human Mendelian disorders

BACKGROUND

Mapping the breakpoints in de novo balanced chromosomal translocations (BCT) in symptomatic individuals provides a unique opportunity to identify in an unbiased way the likely causative genetic defect and thus find novel human disease candidate genes. Our aim was to fine-map breakpoints of de novo BCTs in a case series of nine patients.

METHODS

Shallow whole-genome mate pair sequencing (SGMPS) together with long-range PCR and Sanger sequencing. In one case (BCT disrupting BAHD1 and RET) cDNA analysis was used to verify expression of a fusion transcript in cultured fibroblasts.

RESULTS

In all nine probands 11 disrupted genes were found, that is, EFNA5, EBF3, LARGE, PPP2R5E, TXNDC5, ZNF423, NIPBL, BAHD1, RET, TRPS1 and SLC4A10. Five subjects had translocations that disrupted genes with so far unknown (EFNA5, BAHD1, PPP2R5E, TXNDC5) or poorly delineated impact on the phenotype (SLC4A10, two previous reports of BCT disrupting the gene). The four genes with no previous disease associations (EFNA5, BAHD1, PPP2R5E, TXNDC5), when compared with all human genes by a bootstrap test, had significantly higher pLI (p<0.017) and DOMINO (p<0.02) scores indicating enrichment in genes likely to be intolerant to single copy damage. Inspection of individual pLI and DOMINO scores, and local topologically associating domain structure suggested that EFNA5, BAHD1 and PPP2R5E were particularly good candidates for novel disease loci. The pathomechanism for BAHD1 may involve deregulation of expression due to fusion with RET promoter.

CONCLUSION

SGMPS in symptomatic carriers of BCTs is a powerful approach to delineate novel human gene-disease associations.

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.

Spinal muscular atrophy within Amish and Mennonite populations: Ancestral haplotypes and natural history

We correlate chromosome 5 haplotypes and SMN2 copy number with disease expression in 42 Mennonite and 14 Amish patients with spinal muscular atrophy (SMA). A single haplotype (A1) with 1 copy of SMN2 segregated among all Amish patients. SMN1 deletions segregated on four different Mennonite haplotypes that carried 1 (M1a, M1b, M1c) or 2 (M2) copies of SMN2. DNA microsatellite and microarray data revealed structural similarities among A1, M1a, M1b, and M2. Clinical data were parsed according to both SMN1 genotype and SMN2 copy number (2 copies, n = 44; 3 copies, n = 9; or 4 copies, n = 3). No infant with 2 copies of SMN2 sat unassisted. In contrast, all 9 Mennonites with the M1a/M2 genotype (3 copies of SMN2) sat during infancy at a median age of 7 months, and 5 (56%) walked and dressed independently at median ages of 18 and 36 months, respectively. All are alive at a median age of 11 (range 2–31) years without ventilatory support. Among 13 Amish and 26 Mennonite patients with 2 copies of SMN2 who did not receive feeding or ventilatory support, A1/A1 as compared to M1a/M1a genotype was associated with earlier clinical onset (p = 0.0040) and shorter lifespan (median survival 3.9 versus 5.7 months, p = 0.0314). These phenotypic differences were not explained by variation in SMN1 deletion size or SMN2 coding sequence, which were conserved across haplotypes. Distinctive features of SMA within Plain communities provide a population-specific framework to study variations of disease expression and the impact of disease-modifying therapies administered early in life.

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.

Ether-à-go-go K+ channels: effective modulators of neuronal excitability

Mammalian ether-à-go-go (EAG) channels are voltage-gated K⁺ channels. They are encoded by the KCNH gene family and divided into three subfamilies, eag (Kv10), erg (eag-related gene; Kv11) and elk (eag-like; Kv12). All EAG channel subtypes are expressed in the brain where they effectively modulate neuronal excitability. This Topical Review describes the biophysical properties of each of the EAG channel subtypes, their function in neurons and the neurological diseases induced by EAG channel mutations. In contrast to the function of erg currents in the heart, where they contribute to repolarization of the cardiac action potential, erg currents in neurons are involved in the maintenance of the resting potential, setting of action potential threshold and frequency accommodation. They can even support high frequency firing by preventing a depolarization-induced Na⁺ channel block. EAG channels are modulated differentially, e.g. eag channels by intracellular Ca²⁺ , erg channels by extracellular K⁺ and GPCRs, and elk channels by changes in pH. So far, only currents mediated by erg channels have been recorded in neurons with the help of selective blockers. Neuronal eag and elk currents have not been isolated due to the lack of suitable channel blockers. However, findings in KO mice indicate a physiological role of eag1 currents in synaptic transmission and an involvement of elk2 currents in cognitive performance. Human eag1 and eag2 gain-of-function mutations underlie syndromes associated with epileptic seizures.

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.

Arthropod toxins acting on neuronal potassium channels

Arthropod venoms are a rich mixture of biologically active compounds exerting different physiological actions across diverse phyla and affecting multiple organ systems including the central nervous system. Venom compounds can inhibit or activate ion channels, receptors and transporters with high specificity and affinity providing essential insights into ion channel function. In this review, we focus on arthropod toxins (scorpions, spiders, bees and centipedes) acting on neuronal potassium channels. A brief description of the K⁺ channels classification and structure is included and a compendium of neuronal K⁺ channels and the arthropod toxins that modify them have been listed. This article is part of the Special Issue entitled 'Venom-derived Peptides as Pharmacological Tools.'

Effects of Temperature on Heteromeric Kv11.1a/1b and Kv11.3 Channels

Kv11.1 channels are crucial in cardiac physiology, and there is increasing evidence of physiological roles of different Kv11 channels outside the heart. The HERG (human Kv11.1a) channel has previously been shown to carry substantially more current at elevated temperatures, and we have now comparably investigated the temperature dependence of neuronal Kv11.3 channels and the more ubiquitous heteromeric Kv11.1a/1b channels. Transiently expressed rat Kv11 channels were studied at 21°C, 30°C, and 35°C. At near-physiological temperature, the maximal sustained outward current density was almost three times the mean value obtained at room temperature for Kv11.1a/1b, and increased by ∼150% for Kv11.3. For both channels, reduced inactivation contributed to the current increase at higher temperature. Elevated temperature moved Kv11.1a/1b isochronal activation curves to more negative potentials, but shifted the potential of half-maximal Kv11.3 channel activation to more depolarized values and reduced its voltage sensitivity. Thus, increased temperature stabilized the open state over the closed state of Kv11.1a/1b channels and exerted the opposite effect on Kv11.3 channel activation. Both Kv11 channels exhibited an overall high temperature sensitivity of most gating parameters, with remarkably high Q10 factors of ∼5 for the rate of Kv11.1a/1b activation. The Q10 factors for Kv11.3 gating were more uniform, but still higher for activation than for inactivation kinetics. The results demonstrate that characteristic differences between Kv11.1a/1b and Kv11.3 determined at room temperature do not necessarily apply to physiological conditions. The data provided here can aid in the design of models that will enhance our understanding of the role of Kv11 currents in excitable cells.

Potassium channel gene associations with joint processing speed and white matter impairments in schizophrenia

Patients with schizophrenia show decreased processing speed on neuropsychological testing and decreased white matter integrity as measured by diffusion tensor imaging, two traits shown to be both heritable and genetically associated indicating that there may be genes that influence both traits as well as schizophrenia disease risk. The potassium channel gene family is a reasonable candidate to harbor such a gene given the prominent role potassium channels play in the central nervous system in signal transduction, particularly in myelinated axons. We genotyped members of the large potassium channel gene family focusing on putatively functional single nucleotide polymorphisms (SNPs) in a population of 363 controls, 194 patients with schizophrenia spectrum disorder (SSD) and 28 patients with affective disorders with psychotic features who completed imaging and neuropsychological testing. We then performed three association analyses using three phenotypes - processing speed, whole-brain white matter fractional anisotropy (FA) and schizophrenia spectrum diagnosis. We extracted SNPs showing an association at a nominal P value of <0.05 with all three phenotypes in the expected direction: decreased processing speed, decreased FA and increased risk of SSD. A single SNP, rs8234, in the 3' untranslated region of voltage-gated potassium channel subfamily Q member 1 (KCNQ1) was identified. Rs8234 has been shown to affect KCNQ1 expression levels, and KCNQ1 levels have been shown to affect neuronal action potentials. This exploratory analysis provides preliminary data suggesting that KCNQ1 may contribute to the shared risk for diminished processing speed, diminished white mater integrity and increased risk of schizophrenia.

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.

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.

Targeted massively parallel sequencing of autism spectrum disorder-associated genes in a case control cohort reveals rare loss-of-function risk variants

Background

Autism spectrum disorder (ASD) is highly heritable, yet genome-wide association studies (GWAS), copy number variation screens, and candidate gene association studies have found no single factor accounting for a large percentage of genetic risk. ASD trio exome sequencing studies have revealed genes with recurrent de novo loss-of-function variants as strong risk factors, but there are relatively few recurrently affected genes while as many as 1000 genes are predicted to play a role. As such, it is critical to identify the remaining rare and low-frequency variants contributing to ASD.

Methods

We have utilized an approach of prioritization of genes by GWAS and follow-up with massively parallel sequencing in a case-control cohort. Using a previously reported ASD noise reduction GWAS analyses, we prioritized 837 RefSeq genes for custom targeting and sequencing. We sequenced the coding regions of those genes in 2071 ASD cases and 904 controls of European white ancestry. We applied comprehensive annotation to identify single variants which could confer ASD risk and also gene-based association analysis to identify sets of rare variants associated with ASD.

Results

We identified a significant over-representation of rare loss-of-function variants in genes previously associated with ASD, including a de novo premature stop variant in the well-established ASD candidate gene RBFOX1. Furthermore, ASD cases were more likely to have two damaging missense variants in candidate genes than controls. Finally, gene-based rare variant association implicates genes functioning in excitatory neurotransmission and neurite outgrowth and guidance pathways including CACNAD2, KCNH7, and NRXN1.

Conclusions

We find suggestive evidence that rare variants in synaptic genes are associated with ASD and that loss-of-function mutations in ASD candidate genes are a major risk factor, and we implicate damaging mutations in glutamate signaling receptors and neuronal adhesion and guidance molecules. Furthermore, the role of de novo mutations in ASD remains to be fully investigated as we identified the first reported protein-truncating variant in RBFOX1 in ASD. Overall, this work, combined with others in the field, suggests a convergence of genes and molecular pathways underlying ASD etiology.

Electronic supplementary material

The online version of this article (doi:10.1186/s13229-015-0034-z) contains supplementary material, which is available to authorized users.

Rare variants in neuronal excitability genes influence risk for bipolar disorder

We sequenced the genomes of 200 individuals from 41 families multiply affected with bipolar disorder (BD) to identify contributions of rare variants to genetic risk. We initially focused on 3,087 candidate genes with known synaptic functions or prior evidence from genome-wide association studies. BD pedigrees had an increased burden of rare variants in genes encoding neuronal ion channels, including subunits of GABAA receptors and voltage-gated calcium channels. Four uncommon coding and regulatory variants also showed significant association, including a missense variant in GABRA6. Targeted sequencing of 26 of these candidate genes in an additional 3,014 cases and 1,717 controls confirmed rare variant associations in ANK3, CACNA1B, CACNA1C, CACNA1D, CACNG2, CAMK2A, and NGF. Variants in promoters and 5' and 3' UTRs contributed more strongly than coding variants to risk for BD, both in pedigrees and in the case-control cohort. The genes and pathways identified in this study regulate diverse aspects of neuronal excitability. We conclude that rare variants in neuronal excitability genes contribute to risk for BD.

Whole genome/exome sequencing in mood and psychotic disorders

Recent developments in DNA sequencing technologies have allowed for genetic studies using whole genome or exome analysis, and these have been applied in the study of mood and psychotic disorders, including bipolar disorder, depression, schizophrenia, and schizoaffective disorder. In this review, the current situation, recent findings, methodological problems, and future directions of whole genome/exome analysis studies of these disorders are summarized. Whole genome/exome studies of bipolar disorder have included pedigree analysis and case-control studies, demonstrating the role of previously implicated pathways, such as calcium signaling, cyclic adenosine monophosphate response element binding protein (CREB) signaling, and potassium channels. Extensive analysis of trio families and case-control studies showed that de novo mutations play a role in the genetic architecture of schizophrenia and indicated that mutations in several molecular pathways, including chromatin regulation, activity-regulated cytoskeleton, post-synaptic density, N-methyl-D-aspartate receptor, and targets of fragile X mental retardation protein, are associated with this disorder. Depression is a heterogeneous group of diseases and studies using exome analysis have been conducted to identify rare mutations causing Mendelian diseases that accompany depression. In the near future, clarification of the genetic architecture of bipolar disorder and schizophrenia is expected. Identification of causative mutations using these new technologies will facilitate neurobiological studies of these disorders.

Toward a Deeper Understanding of the Genetics of Bipolar Disorder

Bipolar disorder is a common, complex psychiatric disorder characterized by mania and depression. The disease aggregates in families, but despite much effort, it has been difficult to delineate the basic genetic model or identify specific genetic risk factors. Not only single gene Mendelian transmission and common variant hypotheses but also multivariate threshold models and oligogenic quasi-Mendelian modes of inheritance have dominated the discussion at times. Almost complete sequence information of the human genome and falling sequencing costs now offer the opportunity to test these models in families in which the disorder is transmitted over several generations. Exome-wide sequencing studies have revealed an astonishing number of rare and potentially damaging mutations in brain-expressed genes that could have contributed to the disease manifestation. However, the statistical analysis of these data has been challenging, because genetic risk factors displayed a high degree of dissimilarity across families. This scenario is not unique to bipolar disorder, but similar results have also been found in schizophrenia, a potentially related psychiatric disorder. Recently, our group has published data which supported an oligogenic genetic model of transmission in a family with bipolar disorder. In this family, three affected siblings shared rare, damaging mutations in multiple genes, which were linked to stress response pathways. These pathways are also the target for drugs frequently used to treat bipolar disorder. This article discusses these findings in the context of previously proclaimed disease models and suggests future research directions, including biological confirmation and phenotype stratification as an approach to disease heterogeneity.