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.

Biallelic variants in POPDC2 cause a novel autosomal recessive syndrome presenting with cardiac conduction defects and variable hypertrophic cardiomyopathy.

POPDC2 encodes for the Popeye domain-containing protein 2 which has an important role in cardiac pacemaking and conduction, due in part to its cAMP-dependent binding and regulation of TREK-1 potassium channels. Loss of Popdc2 in mice results in sinus pauses and bradycardia and morpholino knockdown of popdc2 in zebrafish results in atrioventricular (AV) block. We identified bi-allelic variants in POPDC2 in 4 families that presented with a phenotypic spectrum consisting of sinus node dysfunction, AV conduction defects and hypertrophic cardiomyopathy. Using homology modelling we show that the identified POPDC2 variants are predicted to diminish the ability of POPDC2 to bind cAMP. In in vitro electrophysiological studies we demonstrated that, while co-expression of wild-type POPDC2 with TREK-1 increased TREK-1 current density, POPDC2 variants found in the patients failed to increase TREK-1 current density. While patient muscle biopsy did not show clear myopathic disease, it showed significant reduction of the expression of both POPDC1 and POPDC2, suggesting that stability and/or membrane trafficking of the POPDC1-POPDC2 complex is impaired by pathogenic variants in any of the two proteins. Single-cell RNA sequencing from human hearts demonstrated that co-expression of POPDC1 and 2 was most prevalent in AV node, AV node pacemaker and AV bundle cells. Sinoatrial node cells expressed POPDC2 abundantly, but expression of POPDC1 was sparse. Together, these results concur with predisposition to AV node disease in humans with loss-of-function variants in POPDC1 and POPDC2 and presence of sinus node disease in POPDC2, but not in POPDC1 related disease in human. Using population-level genetic data of more than 1 million individuals we showed that none of the familial variants were associated with clinical outcomes in heterozygous state, suggesting that heterozygous family members are unlikely to develop clinical manifestations and therefore might not necessitate clinical follow-up. Our findings provide evidence for POPDC2 as the cause of a novel Mendelian autosomal recessive cardiac syndrome, consistent with previous work showing that mice and zebrafish deficient in functional POPDC2 display sinus and AV node dysfunction.

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.

Sequence variant affects GCSAML splicing, mast cell specific proteins, and risk of urticaria.

Urticaria is a skin disorder characterized by outbreaks of raised pruritic wheals. In order to identify sequence variants associated with urticaria, we performed a meta-analysis of genome-wide association studies for urticaria with a total of 40,694 cases and 1,230,001 controls from Iceland, the UK, Finland, and Japan. We also performed transcriptome- and proteome-wide analyses in Iceland and the UK. We found nine sequence variants at nine loci associating with urticaria. The variants are at genes participating in type 2 immune responses and/or mast cell biology (CBLB, FCER1A, GCSAML, STAT6, TPSD1, ZFPM1), the innate immunity (C4), and NF-κB signaling. The most significant association was observed for the splice-donor variant rs56043070[A] (hg38: chr1:247556467) in GCSAML (MAF = 6.6%, OR = 1.24 (95%CI: 1.20-1.28), P-value = 3.6 × 10-44). We assessed the effects of the variants on transcripts, and levels of proteins relevant to urticaria pathophysiology. Our results emphasize the role of type 2 immune response and mast cell activation in the pathogenesis of urticaria. Our findings may point to an IgE-independent urticaria pathway that could help address unmet clinical need.

Rare variants with large effects provide functional insights into the pathology of migraine subtypes, with and without aura.

Migraine is a complex neurovascular disease with a range of severity and symptoms, yet mostly studied as one phenotype in genome-wide association studies (GWAS). Here we combine large GWAS datasets from six European populations to study the main migraine subtypes, migraine with aura (MA) and migraine without aura (MO). We identified four new MA-associated variants (in PRRT2, PALMD, ABO and LRRK2) and classified 13 MO-associated variants. Rare variants with large effects highlight three genes. A rare frameshift variant in brain-expressed PRRT2 confers large risk of MA and epilepsy, but not MO. A burden test of rare loss-of-function variants in SCN11A, encoding a neuron-expressed sodium channel with a key role in pain sensation, shows strong protection against migraine. Finally, a rare variant with cis-regulatory effects on KCNK5 confers large protection against migraine and brain aneurysms. Our findings offer new insights with therapeutic potential into the complex biology of migraine and its subtypes.

A genome-wide meta-analysis identifies 50 genetic loci associated with carpal tunnel syndrome.

Carpal tunnel syndrome (CTS) is the most common entrapment neuropathy and has a largely unknown underlying biology. In a genome-wide association study of CTS (48,843 cases and 1,190,837 controls), we found 53 sequence variants at 50 loci associated with the syndrome. The most significant association is with a missense variant (p.Glu366Lys) in SERPINA1 that protects against CTS (P = 2.9 × 10-24, OR = 0.76). Through various functional analyses, we conclude that at least 22 genes mediate CTS risk and highlight the role of 19 CTS variants in the biology of the extracellular matrix. We show that the genetic component to the risk is higher in bilateral/recurrent/persistent cases than nonrecurrent/nonpersistent cases. Anthropometric traits including height and BMI are genetically correlated with CTS, in addition to early hormonal-replacement therapy, osteoarthritis, and restlessness. Our findings suggest that the components of the extracellular matrix play a key role in the pathogenesis of CTS.

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.

Rare SLC13A1 variants associate with intervertebral disc disorder highlighting role of sulfate in disc pathology.

Back pain is a common and debilitating disorder with largely unknown underlying biology. Here we report a genome-wide association study of back pain using diagnoses assigned in clinical practice; dorsalgia (119,100 cases, 909,847 controls) and intervertebral disc disorder (IDD) (58,854 cases, 922,958 controls). We identify 41 variants at 33 loci. The most significant association (ORIDD = 0.92, P = 1.6 × 10-39; ORdorsalgia = 0.92, P = 7.2 × 10-15) is with a 3'UTR variant (rs1871452-T) in CHST3, encoding a sulfotransferase enzyme expressed in intervertebral discs. The largest effects on IDD are conferred by rare (MAF = 0.07 - 0.32%) loss-of-function (LoF) variants in SLC13A1, encoding a sodium-sulfate co-transporter (LoF burden OR = 1.44, P = 3.1 × 10-11); variants that also associate with reduced serum sulfate. Genes implicated by this study are involved in cartilage and bone biology, as well as neurological and inflammatory processes.

A genome-wide meta-analysis uncovers six sequence variants conferring risk of vertigo.

Vertigo is the leading symptom of vestibular disorders and a major risk factor for falls. In a genome-wide association study of vertigo (Ncases = 48,072, Ncontrols = 894,541), we uncovered an association with six common sequence variants in individuals of European ancestry, including missense variants in ZNF91, OTOG, OTOGL, and TECTA, and a cis-eQTL for ARMC9. The association of variants in ZNF91, OTOGL, and OTOP1 was driven by an association with benign paroxysmal positional vertigo. Using previous reports of sequence variants associating with age-related hearing impairment and motion sickness, we found eight additional variants that associate with vertigo. Although disorders of the auditory and the vestibular system may co-occur, none of the six genome-wide significant vertigo variants were associated with hearing loss and only one was associated with age-related hearing impairment. Our results uncovered sequence variants associating with vertigo in a genome-wide association study and implicated genes with known roles in inner ear development, maintenance, and disease.

A meta-analysis uncovers the first sequence variant conferring risk of Bell's palsy.

Bell's palsy is the most common cause of unilateral facial paralysis and is defined as an idiopathic and acute inability to control movements of the facial muscles on the affected side. While the pathogenesis remains unknown, previous studies have implicated post-viral inflammation and resulting compression of the facial nerve. Reported heritability estimates of 4-14% suggest a genetic component in the etiology and an autosomal dominant inheritance has been proposed. Here, we report findings from a meta-analysis of genome-wide association studies uncovering the first unequivocal association with Bell's palsy (rs9357446-A; P = 6.79 × 10-23, OR = 1.23; Ncases = 4714, Ncontrols = 1,011,520). The variant also confers risk of intervertebral disc disorders (P = 2.99 × 10-11, OR = 1.04) suggesting a common pathogenesis in part or a true pleiotropy.

Sequence Variants in TAAR5 and Other Loci Affect Human Odor Perception and Naming.

Olfactory receptor (OR) genes in humans form a special class characterized by unusually high DNA sequence diversity, which should give rise to differences in perception and behavior. In the largest genome-wide association study to date based on olfactory testing, we investigated odor perception and naming with smell tasks performed by 9,122 Icelanders, with replication in a separate sample of 2,204 individuals. We discovered an association between a low-frequency missense variant in TAAR5 and reduced intensity rating of fish odor containing trimethylamine (p.Ser95Pro, pcombined = 5.6 × 10-15). We demonstrate that TAAR5 genotype affects aversion to fish odor, reflected by linguistic descriptions of the odor and pleasantness ratings. We also discovered common sequence variants in two canonical olfactory receptor loci that associate with increased intensity and naming of licorice odor (trans-anethole: lead variant p.Lys233Asn in OR6C70, pcombined = 8.8 × 10-16 and pcombined = 1.4 × 10-9) and enhanced naming of cinnamon (trans-cinnamaldehyde; intergenic variant rs317787-T, pcombined = 5.0 × 10-17). Together, our results show that TAAR5 genotype variation influences human odor responses and highlight that sequence diversity in canonical OR genes can lead to enhanced olfactory ability, in contrast to the view that greater tolerance for mutations in the human OR repertoire leads to diminished function.

GraphTyper2 enables population-scale genotyping of structural variation using pangenome graphs.

Analysis of sequence diversity in the human genome is fundamental for genetic studies. Structural variants (SVs) are frequently omitted in sequence analysis studies, although each has a relatively large impact on the genome. Here, we present GraphTyper2, which uses pangenome graphs to genotype SVs and small variants using short-reads. Comparison to the syndip benchmark dataset shows that our SV genotyping is sensitive and variant segregation in families demonstrates the accuracy of our approach. We demonstrate that incorporating public assembly data into our pipeline greatly improves sensitivity, particularly for large insertions. We validate 6,812 SVs on average per genome using long-read data of 41 Icelanders. We show that GraphTyper2 can simultaneously genotype tens of thousands of whole-genomes by characterizing 60 million small variants and half a million SVs in 49,962 Icelanders, including 80 thousand SVs with high-confidence.

A PRPH splice-donor variant associates with reduced sural nerve amplitude and risk of peripheral neuropathy.

Nerve conduction (NC) studies generate measures of peripheral nerve function that can reveal underlying pathology due to axonal loss, demyelination or both. We perform a genome-wide association study of sural NC amplitude and velocity in 7045 Icelanders and find a low-frequency splice-donor variant in PRPH (c.996+1G>A; MAF = 1.32%) associating with decreased NC amplitude but not velocity. PRPH encodes peripherin, an intermediate filament (IF) protein involved in cytoskeletal development and maintenance of neurons. Through RNA and protein studies, we show that the variant leads to loss-of-function (LoF), as when over-expressed in a cell line devoid of other IFs, it does not allow formation of the normal filamentous structure of peripherin, yielding instead punctate protein inclusions. Recall of carriers for neurological assessment confirms that from an early age, homozygotes have significantly lower sural NC amplitude than non-carriers and are at risk of a mild, early-onset, sensory-negative, axonal polyneuropathy.

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-8), 47% less corpus callosum (CC) volume (ß = -2.4SD, P = 5.5 × 10-10) and lower brain-wide fractional anisotropy (P = 6.7 × 10-4). 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.