GWAS meta-analysis reveals key risk loci in essential tremor pathogenesis

Essential tremor (ET) is a prevalent neurological disorder with a largely unknown underlying biology. In this genome-wide association study meta-analysis, comprising 16,480 ET cases and 1,936,173 controls from seven datasets, we identify 12 sequence variants at 11 loci. Evaluating mRNA expression, splicing, plasma protein levels, and coding effects, we highlight seven putative causal genes at these loci, including CA3 and CPLX1. CA3 encodes Carbonic Anhydrase III and carbonic anhydrase inhibitors have been shown to decrease tremors. CPLX1, encoding Complexin-1, regulates neurotransmitter release. Through gene-set enrichment analysis, we identify a significant association with specific cell types, including dopaminergic and GABAergic neurons, as well as biological processes like Rho GTPase signaling. Genetic correlation analyses reveals a positive association between ET and Parkinson's disease, depression, and anxiety-related phenotypes. This research uncovers risk loci, enhancing our knowledge of the complex genetics of this common but poorly understood disorder, and highlights CA3 and CPLX1 as potential therapeutic targets.

Sequence variants affecting voice pitch in humans

The genetic basis of the human vocal system is largely unknown, as are the sequence variants that give rise to individual differences in voice and speech. Here, we couple data on diversity in the sequence of the genome with voice and vowel acoustics in speech recordings from 12,901 Icelanders. We show how voice pitch and vowel acoustics vary across the life span and correlate with anthropometric, physiological, and cognitive traits. We found that voice pitch and vowel acoustics have a heritable component and discovered correlated common variants in ABCC9 that associate with voice pitch. The ABCC9 variants also associate with adrenal gene expression and cardiovascular traits. By showing that voice and vowel acoustics are influenced by genetics, we have taken important steps toward understanding the genetics and evolution of the human vocal system.

Thirty novel sequence variants impacting human intracranial volume

Intracranial volume, measured through magnetic resonance imaging and/or estimated from head circumference, is heritable and correlates with cognitive traits and several neurological disorders. We performed a genome-wide association study meta-analysis of intracranial volume (n = 79 174) and found 64 associating sequence variants explaining 5.0% of its variance. We used coding variation, transcript and protein levels, to uncover 12 genes likely mediating the effect of these variants, including GLI3 and CDK6 that affect cranial synostosis and microcephaly, respectively. Intracranial volume correlates genetically with volumes of cortical and sub-cortical regions, cognition, learning, neonatal and neurological traits. Parkinson's disease cases have greater and attention deficit hyperactivity disorder cases smaller intracranial volume than controls. Our Mendelian randomization studies indicate that intracranial volume associated variants either increase the risk of Parkinson's disease and decrease the risk of attention deficit hyperactivity disorder and neuroticism or correlate closely with a confounder.

Dose response of the 16p11.2 distal copy number variant on intracranial volume and basal ganglia

Carriers of large recurrent copy number variants (CNVs) have a higher risk of developing neurodevelopmental disorders. The 16p11.2 distal CNV predisposes carriers to e.g., autism spectrum disorder and schizophrenia. We compared subcortical brain volumes of 12 16p11.2 distal deletion and 12 duplication carriers to 6882 non-carriers from the large-scale brain Magnetic Resonance Imaging collaboration, ENIGMA-CNV. After stringent CNV calling procedures, and standardized FreeSurfer image analysis, we found negative dose-response associations with copy number on intracranial volume and on regional caudate, pallidum and putamen volumes (β = -0.71 to -1.37; P < 0.0005). In an independent sample, consistent results were obtained, with significant effects in the pallidum (β = -0.95, P = 0.0042). The two data sets combined showed significant negative dose-response for the accumbens, caudate, pallidum, putamen and ICV (P = 0.0032, 8.9 × 10-6, 1.7 × 10-9, 3.5 × 10-12 and 1.0 × 10-4, respectively). Full scale IQ was lower in both deletion and duplication carriers compared to non-carriers. This is the first brain MRI study of the impact of the 16p11.2 distal CNV, and we demonstrate a specific effect on subcortical brain structures, suggesting a neuropathological pattern underlying the neurodevelopmental syndromes.

MAP1B mutations cause intellectual disability and extensive white matter deficit

Discovery of coding variants in genes that confer risk of neurodevelopmental disorders is an important step towards understanding the pathophysiology of these disorders. Whole-genome sequencing of 31,463 Icelanders uncovers a frameshift variant (E712KfsTer10) in microtubule-associated protein 1B (MAP1B) that associates with ID/low IQ in a large pedigree (genome-wide corrected P = 0.022). Additional stop-gain variants in MAP1B (E1032Ter and R1664Ter) validate the association with ID and IQ. Carriers have 24% less white matter (WM) volume (β = −2.1SD, P = 5.1 × 10⁻⁸), 47% less corpus callosum (CC) volume (β = −2.4SD, P = 5.5 × 10⁻¹⁰) and lower brain-wide fractional anisotropy (P = 6.7 × 10⁻⁴). In summary, we show that loss of MAP1B function affects general cognitive ability through a profound, brain-wide WM deficit with likely disordered or compromised axons.

Intellectual disability (ID) is characterized by an intelligence quotient of below 70 and impaired adaptive skills. Here, analyzing whole genome sequences from 31,463 Icelanders, Walters et al. identify variants in MAP1B associated with ID and extensive brain-wide white matter deficits.

Truncating mutations in RBM12 are associated with psychosis

Thus far, a handful of highly penetrant mutations conferring risk of psychosis have been discovered. Here we used whole-genome sequencing and long-range phasing to investigate an Icelandic kindred containing ten individuals with psychosis (schizophrenia, schizoaffective disorder or psychotic bipolar disorder). We found that all affected individuals carry RBM12 (RNA-binding-motif protein 12) c.2377G>T (P = 2.2 × 10^-4), a nonsense mutation that results in the production of a truncated protein lacking a predicted RNA-recognition motif. We replicated the association in a Finnish family in which a second RBM12 truncating mutation (c.2532delT) segregates with psychosis (P = 0.020). c.2377G>T is not fully penetrant for psychosis; however, we found that carriers unaffected by psychosis resemble patients with schizophrenia in their non-psychotic psychiatric disorder and neuropsychological test profile (P = 0.0043) as well as in their life outcomes (including an increased chance of receiving disability benefits, P = 0.011). As RBM12 has not previously been linked to psychosis, this work provides new insight into psychiatric disease.

15q11.2 CNV affects cognitive, structural and functional correlates of dyslexia and dyscalculia

Several copy number variants have been associated with neuropsychiatric disorders and these variants have been shown to also influence cognitive abilities in carriers unaffected by psychiatric disorders. Previously, we associated the 15q11.2(BP1-BP2) deletion with specific learning disabilities and a larger corpus callosum. Here we investigate, in a much larger sample, the effect of the 15q11.2(BP1-BP2) deletion on cognitive, structural and functional correlates of dyslexia and dyscalculia. We report that the deletion confers greatest risk of the combined phenotype of dyslexia and dyscalculia. We also show that the deletion associates with a smaller left fusiform gyrus. Moreover, tailored functional magnetic resonance imaging experiments using phonological lexical decision and multiplication verification tasks demonstrate altered activation in the left fusiform and the left angular gyri in carriers. Thus, by using convergent evidence from neuropsychological testing, and structural and functional neuroimaging, we show that the 15q11.2(BP1-BP2) deletion affects cognitive, structural and functional correlates of both dyslexia and dyscalculia.

CNVs conferring risk of autism or schizophrenia affect cognition in controls

In a small fraction of patients with schizophrenia or autism, alleles of copy-number variants (CNVs) in their genomes are probably the strongest factors contributing to the pathogenesis of the disease. These CNVs may provide an entry point for investigations into the mechanisms of brain function and dysfunction alike. They are not fully penetrant and offer an opportunity to study their effects separate from that of manifest disease. Here we show in an Icelandic sample that a few of the CNVs clearly alter fecundity (measured as the number of children by age 45). Furthermore, we use various tests of cognitive function to demonstrate that control subjects carrying the CNVs perform at a level that is between that of schizophrenia patients and population controls. The CNVs do not all affect the same cognitive domains, hence the cognitive deficits that drive or accompany the pathogenesis vary from one CNV to another. Controls carrying the chromosome 15q11.2 deletion between breakpoints 1 and 2 (15q11.2(BP1-BP2) deletion) have a history of dyslexia and dyscalculia, even after adjusting for IQ in the analysis, and the CNV only confers modest effects on other cognitive traits. The 15q11.2(BP1-BP2) deletion affects brain structure in a pattern consistent with both that observed during first-episode psychosis in schizophrenia and that of structural correlates in dyslexia.

Variant in the sequence of the LINGO1 gene confers risk of essential tremor

We identified a marker in LINGO1 showing genome-wide significant association (P = 1.2 x 10(-9), odds ratio = 1.55) with essential tremor. LINGO1 has potent, negative regulatory influences on neuronal survival and is also important in regulating both central-nervous-system axon regeneration and oligodendrocyte maturation. Increased axon integrity observed in Lingo1 mouse [corrected] knockout models highlights the potential role of LINGO1 in the pathophysiology of ET [corrected]

Disruption of the neurexin 1 gene is associated with schizophrenia

Deletions within the neurexin 1 gene (NRXN1; 2p16.3) are associated with autism and have also been reported in two families with schizophrenia. We examined NRXN1, and the closely related NRXN2 and NRXN3 genes, for copy number variants (CNVs) in 2977 schizophrenia patients and 33 746 controls from seven European populations (Iceland, Finland, Norway, Germany, The Netherlands, Italy and UK) using microarray data. We found 66 deletions and 5 duplications in NRXN1, including a de novo deletion: 12 deletions and 2 duplications occurred in schizophrenia cases (0.47%) compared to 49 and 3 (0.15%) in controls. There was no common breakpoint and the CNVs varied from 18 to 420 kb. No CNVs were found in NRXN2 or NRXN3. We performed a Cochran-Mantel-Haenszel exact test to estimate association between all CNVs and schizophrenia (P = 0.13; OR = 1.73; 95% CI 0.81-3.50). Because the penetrance of NRXN1 CNVs may vary according to the level of functional impact on the gene, we next restricted the association analysis to CNVs that disrupt exons (0.24% of cases and 0.015% of controls). These were significantly associated with a high odds ratio (P = 0.0027; OR 8.97, 95% CI 1.8-51.9). We conclude that NRXN1 deletions affecting exons confer risk of schizophrenia.

Whole-genome analysis of histone H3 lysine 4 and lysine 27 methylation in human embryonic stem cells

We mapped Polycomb-associated H3K27 trimethylation (H3K27me3) and Trithorax-associated H3K4 trimethylation (H3K4me3) across the whole genome in human embryonic stem (ES) cells. The vast majority of H3K27me3 colocalized on genes modified with H3K4me3. These commodified genes displayed low expression levels and were enriched in developmental function. Another significant set of genes lacked both modifications and was also expressed at low levels in ES cells but was enriched for gene function in physiological responses rather than development. Commodified genes could change expression levels rapidly during differentiation, but so could a substantial number of genes in other modification categories. SOX2, POU5F1, and NANOG, pluripotency-associated genes, shifted from modification by H3K4me3 alone to colocalization of both modifications as they were repressed during differentiation. Our results demonstrate that H3K27me3 modifications change during early differentiation, both relieving existing repressive domains and imparting new ones, and that colocalization with H3K4me3 is not restricted to pluripotent cells.

A study of the relationships between oligonucleotide properties and hybridization signal intensities from NimbleGen microarray datasets

Well-defined relationships between oligonucleotide properties and hybridization signal intensities (HSI) can aid chip design, data normalization and true biological knowledge discovery. We clarify these relationships using the data from two microarray experiments containing over three million probes from 48 high-density chips. We find that melting temperature (T_m) has the most significant effect on HSI while length for the long oligonucleotides studied has very little effect. Analysis of positional effect using a linear model provides evidence that the protruding ends of probes contribute more than tethered ends to HSI, which is further validated by specifically designed match fragment sliding and extension experiments. The impact of sequence similarity (SeqS) on HSI is not significant in comparison with other oligonucleotide properties. Using regression and regression tree analysis, we prioritize these oligonucleotide properties based on their effects on HSI. The implications of our discoveries for the design of unbiased oligonucleotides are discussed. We propose that isothermal probes designed by varying the length is a viable strategy to reduce sequence bias, though imposing selection constraints on other oligonucleotide properties is also essential.

Induced pluripotent stem cell lines derived from human somatic cells

Somatic cell nuclear transfer allows trans-acting factors present in the mammalian oocyte to reprogram somatic cell nuclei to an undifferentiated state. We show that four factors (OCT4, SOX2, NANOG, and LIN28) are sufficient to reprogram human somatic cells to pluripotent stem cells that exhibit the essential characteristics of embryonic stem (ES) cells. These induced pluripotent human stem cells have normal karyotypes, express telomerase activity, express cell surface markers and genes that characterize human ES cells, and maintain the developmental potential to differentiate into advanced derivatives of all three primary germ layers. Such induced pluripotent human cell lines should be useful in the production of new disease models and in drug development, as well as for applications in transplantation medicine, once technical limitations (for example, mutation through viral integration) are eliminated.

Dynamics of Yeast Myosin I

Cortical actin patches are dynamic structures required for endocytosis in yeast. Recent studies have shown that components of cortical patches localize to the plasma membrane in a precisely orchestrated manner, and their movements at and away from the plasma membrane may define the endocytic membrane invagination and vesicle scission events, respectively. Here, through live-cell imaging, we analyze the dynamics of the highly conserved class I unconventional myosin, Myo5, which also localizes to cortical patches and is known to be involved in endocytosis and actin nucleation. Myo5 exhibits a pattern of dynamic localization different from all cortical patch components analyzed to date. Myo5 associates with cortical patches only transiently and remains stationary during its brief cortical lifespan. The peak of Myo5 association with cortical patches immediately precedes the fast movement of Arp2/3 complex-associated structures away from the plasma membrane, thus correlating precisely with the proposed vesicle scission event. To further test the role of Myo5, we generated a temperature-sensitive mutant myo5 allele. In the myo5 mutant cells, Myo5 exhibits a significantly extended cortical lifespan as a result of a general impairment of Myo5 function, and Arp2 patches exhibit an extended slow-movement phase prior to the fast movement toward the cell interior. The myo5 mutant cells are defective in fluid-phase endocytosis and exhibit an increased number of invaginations on the membrane. Based on these results, we hypothesize that the myosin I motor protein facilitates the membrane fusion/vesicle scission event of endocytosis.

A two-tiered mechanism by which Cdc42 controls the localization and activation of an Arp2/3-activating motor complex in yeast

The establishment of cell polarity in budding yeast involves assembly of actin filaments at specified cortical domains. Elucidation of the underlying mechanism requires an understanding of the machinery that controls actin polymerization and how this machinery is in turn controlled by signaling proteins that respond to polarity cues. We showed previously that the yeast orthologue of the Wiskott-Aldrich Syndrome protein, Bee1/Las17p, and the type I myosins are key regulators of cortical actin polymerization. Here, we demonstrate further that these proteins together with Vrp1p form a multivalent Arp2/3-activating complex. During cell polarization, a bifurcated signaling pathway downstream of the Rho-type GTPase Cdc42p recruits and activates this complex, leading to local assembly of actin filaments. One branch, which requires formin homologues, mediates the recruitment of the Bee1p complex to the cortical site where the activated Cdc42p resides. The other is mediated by the p21-activated kinases, which activate the motor activity of myosin-I through phosphorylation. Together, these findings provide insights into the essential processes leading to polarization of the actin cytoskeleton.