Spinal cord extracts of amyotrophic lateral sclerosis spread TDP-43 pathology in cerebral organoids.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder caused by progressive loss of motor neurons and there is currently no effective therapy. Cytoplasmic mislocalization and aggregation of TAR DNA-binding protein 43 kDa (TDP-43) within the CNS is a pathological hallmark in sporadic ALS and prion-like propagation of pathogenic TDP-43 is thought to be implicated in disease progression. However, cell-to-cell transmission of pathogenic TDP-43 in the human CNS has not been confirmed experimentally. Here we used induced pluripotent stem cells (iPSCs)-derived cerebral organoids as recipient CNS tissue model that are anatomically relevant human brain. We injected postmortem spinal cord protein extracts individually from three non-ALS or five sporadic ALS patients containing pathogenic TDP-43 into the cerebral organoids to validate the templated propagation and spreading of TDP-43 pathology in human CNS tissue. We first demonstrated that the administration of spinal cord extracts from an ALS patient induced the formation of TDP-43 pathology that progressively spread in a time-dependent manner in cerebral organoids, suggesting that pathogenic TDP-43 from ALS functioned as seeds and propagated cell-to-cell to form de novo TDP-43 pathology. We also reported that the administration of ALS patient-derived protein extracts caused astrocyte proliferation to form astrogliosis in cerebral organoids, reproducing the pathological feature seen in ALS. Moreover, we showed pathogenic TDP-43 induced cellular apoptosis and that TDP-43 pathology correlated with genomic damage due to DNA double-strand breaks. Thus, our results provide evidence that patient-derived pathogenic TDP-43 can mimic the prion-like propagation of TDP-43 pathology in human CNS tissue. Our findings indicate that our assays with human cerebral organoids that replicate ALS pathophysiology have a promising strategy for creating readouts that could be used in future drug discovery efforts against ALS.

Rare-variant and polygenic analyses of amyotrophic lateral sclerosis in the French-Canadian genome.

PURPOSE: The genetic etiology of amyotrophic lateral sclerosis (ALS) includes few rare, large-effect variants and potentially many common, small-effect variants per case. The genetic risk liability for ALS might require a threshold comprised of a certain amount of variants. Here, we tested the degree to which risk for ALS was affected by rare variants in ALS genes, polygenic risk score, or both. METHODS: 335 ALS cases and 356 controls from Québec, Canada were concurrently tested by microarray genotyping and targeted sequencing of ALS genes known at the time of study inception. ALS genome-wide association studies summary statistics were used to estimate an ALS polygenic risk score (PRS). Cases and controls were subdivided into rare-variant heterozygotes and non-heterozygotes. RESULTS: Risk for ALS was significantly associated with PRS and rare variants independently in a logistic regression model. Although ALS PRS predicted a small amount of ALS risk overall, the effect was most pronounced between ALS cases and controls that were not heterozygous for a rare variant in the ALS genes surveyed. CONCLUSION: Both PRS and rare variants in ALS genes impact risk for ALS. PRS for ALS is most informative when rare variants are not observed in ALS genes.

Evolution of a Human-Specific Tandem Repeat Associated with ALS.

Tandem repeats are proposed to contribute to human-specific traits, and more than 40 tandem repeat expansions are known to cause neurological disease. Here, we characterize a human-specific 69 bp variable number tandem repeat (VNTR) in the last intron of WDR7, which exhibits striking variability in both copy number and nucleotide composition, as revealed by long-read sequencing. In addition, greater repeat copy number is significantly enriched in three independent cohorts of individuals with sporadic amyotrophic lateral sclerosis (ALS). Each unit of the repeat forms a stem-loop structure with the potential to produce microRNAs, and the repeat RNA can aggregate when expressed in cells. We leveraged its remarkable sequence variability to align the repeat in 288 samples and uncover its mechanism of expansion. We found that the repeat expands in the 3'-5' direction, in groups of repeat units divisible by two. The expansion patterns we observed were consistent with duplication events, and a replication error called template switching. We also observed that the VNTR is expanded in both Denisovan and Neanderthal genomes but is fixed at one copy or fewer in non-human primates. Evaluating the repeat in 1000 Genomes Project samples reveals that some repeat segments are solely present or absent in certain geographic populations. The large size of the repeat unit in this VNTR, along with our multiplexed sequencing strategy, provides an unprecedented opportunity to study mechanisms of repeat expansion, and a framework for evaluating the roles of VNTRs in human evolution and disease.

Moyamoya Disease Susceptibility Gene RNF213 Regulates Endothelial Barrier Function.

BACKGROUND: Variants in the ring finger protein 213 (RNF213) gene are known to be associated with increased predisposition to cerebrovascular diseases development. Genomic studies have identified RNF213 as a major risk factor of Moyamoya disease in East Asian descendants. However, little is known about the RNF213 (ring finger protein 213) biological functions or its associated pathogenic mechanisms underlying Moyamoya disease. METHODS: To investigate RNF213 loss-of-function effect in endothelial cell, stable RNF213-deficient human cerebral endothelial cells were generated using the CRISPR-Cas9 genome editing technology. RESULTS: In vitro assays, using RNF213 knockout brain endothelial cells, showed clear morphological changes and increased blood-brain barrier permeability. Downregulation and delocalization of essential interendothelial junction proteins involved in the blood-brain barrier maintenance, such as PECAM-1 (platelet endothelial cell adhesion molecule-1), was also observed. Brain endothelial RNF213-deficient cells also showed an abnormal potential to transmigration of leukocytes and secreted high amounts of proinflammatory cytokines. CONCLUSIONS: Taken together, these results indicate that RNF213 could be a key regulator of cerebral endothelium integrity, whose disruption could be an early pathological mechanism leading to Moyamoya disease. This study also further reinforces the importance of blood-brain barrier integrity in the development of Moyamoya disease and other RNF213-associated diseases.

Investigation of the causal relationship between ALS and autoimmune disorders: a Mendelian randomization study.

BACKGROUND: Epidemiological studies have reported an association between amyotrophic lateral sclerosis (ALS) and different autoimmune disorders. This study aims to explore the causal relationship between autoimmune disorders and ALS using Mendelian randomization (MR). METHODS: To test the genetically predicted effects of liability towards immune-related outcomes on ALS risk, we used summary statistics from the largest European genome-wide association studies (GWAS) for these disorders in a two-sample MR setting. To do this, we extracted single nucleotide polymorphisms (SNPs) from the GWAS, which strongly associated with the 12 traits, and queried their effects in a large European ALS GWAS (27,265 cases and 110,881 controls). To avoid bias in our MR instruments related to the complex linkage disequilibrium structure of the human leukocyte antigen (HLA) region, we excluded SNPs within this region from the analyses. We computed inverse-variance weighted (IVW) MR estimates and undertook sensitivity analyses using MR methods robust to horizontal pleiotropy. We also performed a reverse MR analysis testing the causal effects of ALS on the above autoimmune traits. RESULTS: After applying Bonferroni correction for multiple testing, our MR analyses showed that the liability to autoimmune disorders does not affect ALS risk. Our reverse MR analysis also did not support the effects of liability to ALS on other autoimmune disorders. The results of the main IVW MR analyses were generally supported by our sensitivity MR analyses. The variance in the exposures explained by the sets of SNPs used as MR instruments ranged from 8.1 × 10-4 to 0.31. Our MR study was well-powered to detect effects as small as an odds ratio (OR) of 1.045 for ALS in the main MR and as small as an OR of 1.32 in the reverse MR. CONCLUSION: Our MR study does not support a relationship between liability to autoimmune disorders and ALS risk in the European population. The associations observed in epidemiological studies could be partly attributed to shared biology or environmental confounders.

Expanded CAG Repeats in ATXN1, ATXN2, ATXN3, and HTT in the 1000 Genomes Project.

BACKGROUND: Spinocerebellar ataxia types 1, 2, 3 and Huntington disease are neurodegenerative disorders caused by expanded CAG repeats. METHODS: We performed an in-silico analysis of CAG repeats in ATXN1, ATXN2, ATXN3, and HTT using 30× whole-=genome sequencing data of 2504 samples from the 1000 Genomes Project. RESULTS: Seven HTT-positive, 3 ATXN2-positive, 1 ATXN3-positive, and 6 possibly ATXN1-positive samples were identified. No correlation was found between the repeat sizes of the different genes. The distribution of CAG alleles varied by ethnicity. CONCLUSION: Our results suggest that there may be asymptomatic small expanded repeats in almost 0.5% of these populations. © 2020 International Parkinson and Movement Disorder Society.

Evidence for Non-Mendelian Inheritance in Spastic Paraplegia 7.

BACKGROUND: Although the typical inheritance of spastic paraplegia 7 is recessive, several reports have suggested that SPG7 variants may also cause autosomal dominant hereditary spastic paraplegia (HSP). OBJECTIVES: We aimed to conduct an exome-wide genetic analysis on a large Canadian cohort of HSP patients and controls to examine the association of SPG7 and HSP. METHODS: We analyzed 585 HSP patients from 372 families and 1175 controls, including 580 unrelated individuals. Whole-exome sequencing was performed on 400 HSP patients (291 index cases) and all 1175 controls. RESULTS: The frequency of heterozygous pathogenic/likely pathogenic SPG7 variants (4.8%) among unrelated HSP patients was higher than among unrelated controls (1.7%; OR 2.88, 95% CI 1.24-6.66, P = 0.009). The heterozygous SPG7 p.(Ala510Val) variant was found in 3.7% of index patients versus 0.85% in unrelated controls (OR 4.42, 95% CI 1.49-13.07, P = 0.005). Similar results were obtained after including only genetically-undiagnosed patients. We identified four heterozygous SPG7 variant carriers with an additional pathogenic variant in known HSP genes, compared to zero in controls (OR 19.58, 95% CI 1.05-365.13, P = 0.0031), indicating potential digenic inheritance. We further identified four families with heterozygous variants in SPG7 and SPG7-interacting genes (CACNA1A, AFG3L2, and MORC2). Of these, there is especially compelling evidence for epistasis between SPG7 and AFG3L2. The p.(Ile705Thr) variant in AFG3L2 is located at the interface between hexamer subunits, in a hotspot of mutations associated with spinocerebellar ataxia type 28 that affect its proteolytic function. CONCLUSIONS: Our results provide evidence for complex inheritance in SPG7-associated HSP, which may include recessive and possibly dominant and digenic/epistasis forms of inheritance. © 2021 International Parkinson and Movement Disorder Society.

Oligogenicity, C9orf72 expansion, and variant severity in ALS.

"Oligogenic inheritance" is used to describe cases where more than one rare pathogenic variant is observed in the same individual. While multiple variants can alter disease presentation, the necessity of multiple variants to instigate pathogenesis has not been addressed in amyotrophic lateral sclerosis (ALS). We sequenced ALS-associated genes in C9orf72-expansion-positive and negative ALS patients, alongside unaffected controls, to test the importance of oligogenicity and variant deleteriousness in ALS. We found that all groups had similar numbers of rare variants, but that variant severity was significantly higher in C9orf72-negative ALS cases, suggesting sufficiency of C9orf72 expansion to cause ALS alone.

Mutations in ATP13A2 (PARK9) are associated with an amyotrophic lateral sclerosis-like phenotype, implicating this locus in further phenotypic expansion.

BACKGROUND: Amyotrophic lateral sclerosis [1] is a genetically heterogeneous neurodegenerative disorder, characterized by late-onset degeneration of motor neurons leading to progressive limb and bulbar weakness, as well as of the respiratory muscles, which is the primary cause of disease fatality. To date, over 25 genes have been implicated as causative in ALS with C9orf72, SOD1, FUS, and TARDBP accounting for the majority of genetically positive cases. RESULTS: We identified two patients of Italian and French ancestry with a clinical diagnosis of juvenile-onset ALS who were mutation-negative in any of the known ALS causative genes. Starting with the index case, a consanguineous family of Italian origin, we performed whole-exome sequencing and identified candidate pathogenic mutations in 35 genes, 27 of which were homozygous. We next parsed all candidates against a cohort of 3641 ALS cases; only ATP13A2 was found to harbor recessive changes, in a patient with juvenile-onset ALS, similar to the index case. In vivo complementation of ATP13A2 using a zebrafish surrogate model that focused on the assessment of motor neuron morphology and cerebellar integrity confirmed the role of this gene in central and peripheral nervous system maintenance and corroborated the damaging direction of effect of the change detected in the index case of this study. CONCLUSIONS: We here expand the phenotypic spectrum associated with genetic variants in ATP13A2 that previously comprised Kufor-Rakeb syndrome, spastic paraplegia 78, and neuronal ceroid lipofuscinosis type 12 (CLN12), to also include juvenile-onset ALS, as supported by both genetic and functional data. Our findings highlight the importance of establishing a complete genetic profile towards obtaining an accurate clinical diagnosis.

Full sequencing and haplotype analysis of MAPT in Parkinson's disease and rapid eye movement sleep behavior disorder.

BACKGROUND: MAPT haplotypes are associated with PD, but their association with rapid eye movement sleep behavior disorder is unclear. OBJECTIVE: To study the role of MAPT variants in rapid eye movement sleep behavior disorder. METHODS: Two cohorts were included: (A) PD (n = 600), rapid eye movement sleep behavior disorder (n = 613) patients, and controls (n = 981); (B) dementia with Lewy bodies patients with rapid eye movement sleep behavior disorder (n = 271) and controls (n = 950). MAPT-associated variants and the entire coding sequence of MAPT were analyzed. Age-, sex-, and ethnicity-adjusted analyses were performed to examine the association between MAPT, PD, and rapid eye movement sleep behavior disorder. RESULTS: MAPT-H2 variants were associated with PD (odds ratios: 0.62-0.65; P = 0.010-0.019), but not with rapid eye movement sleep behavior disorder. In PD, the H1 haplotype odds ratio was 1.60 (95% confidence interval: 1.12-2.28; P = 0.009), and the H2 odds ratio was 0.68 (95% confidence interval: 0.48-0.96; P = 0.03). The H2/H1 haplotypes were not associated with rapid eye movement sleep behavior disorder. CONCLUSIONS: Our results confirm the protective effect of the MAPT-H2 haplotype in PD, and define its components. Furthermore, our results suggest that MAPT does not play a major role in rapid eye movement sleep behavior disorder, emphasizing different genetic background than in PD in this locus. © 2018 International Parkinson and Movement Disorder Society.

Genetic modifiers of multiple sclerosis progression, severity and onset.

The genetic contribution to clinical outcomes for multiple sclerosis (MS) has yet to be defined. We performed exome sequencing analysis in 100 MS patients presenting opposite extremes of clinical phenotype (discovery cohort), and genotyped variants of interest in 2016 MS patients (replication cohort). Linear and logistic regression analyses were used to identify significant associations with disease course, severity and onset. Our analysis of the discovery cohort nominated 38 variants in 21 genes. Replication analysis identified PSMG4 p.W99R and NLRP5 p.M459I to be associated with disease severity (p=0.002 and 0.008). CACNA1H p.R1871Q was found associated with patients presenting relapsing remitting MS at clinical onset (p=0.028) whereas NLRP5 p.M459I and EIF2AK1 p.K558R were associated with primary progressive disease (p=0.031 and 0.023). In addition, PSMG4 p.W99R and NLRP5 p.R761L were found to correlate with an earlier age at MS clinical onset, and MC1R p.R160W with delayed onset of clinical symptoms (p=0.010-0.041).

Genetic variants in IL2RA and IL7R affect multiple sclerosis disease risk and progression.

Multiple sclerosis (MS) is a common demyelinating neurodegenerative disease with a strong genetic component. Previous studies have associated genetic variants in IL2RA and IL7R in the pathophysiology of the disease. In this study, we describe the association between IL2RA (rs2104286) and IL7R (rs6897932) in the Canadian population. Genotyping 1,978 MS patients and 830 controls failed to identify any significant association between these variants and disease risk. However, stratified analysis for family history of disease and disease course identified a trend towards association for IL2RA in patients without a family history (p = 0.05; odds ratio = 0.77) and a significant association between IL7R and patients who developed progressive MS (PrMS) (p = 0.002; odds ratio = 0.73). Although not statistically significant, the effect of IL2RA (rs2104286) in patients without a family history of MS indicates that the genetic components for familial and sporadic disease are perhaps distinct. This data suggests that the onset of sporadic disease is likely determined by a large number of variants of small effect, whereas MS in patients with a family history of disease is caused by a few deleterious variants. In addition, the significant association between PrMS and rs6897932 indicates that IL7R may not be disease-causing but a determinant of disease course. Further characterization of the effect of IL2RA and IL7R genetic variants in defined MS subtypes is warranted to evaluate the effect of these genes on specific clinical outcomes and to further elucidate the mechanisms of disease onset and progression.

Copy-number variants and polygenic risk for intelligence confer risk for autism spectrum disorder irrespective of their effects on cognitive ability.

Introduction: Rare copy number variants (CNVs) and polygenic risk for intelligence (PRS-IQ) both confer susceptibility for autism spectrum disorder (ASD) but have opposing effects on cognitive ability. The field has struggled to disentangle the effects of these two classes of genomic variants on cognitive ability from their effects on ASD susceptibility, in part because previous studies did not include controls with cognitive measures. We aim to investigate the impact of these genomic variants on ASD risk while adjusting for their known effects on cognitive ability. Methods: In a cohort of 8,426 subjects with ASD and 169,804 controls with cognitive assessments, we found that rare coding CNVs and PRS-IQ increased ASD risk, even after adjusting for their effects on cognitive ability. Results: Bottom decile PRS-IQ and CNVs both decreased cognitive ability but had opposing effects on ASD risk. Models combining both classes of variants showed that the effects of rare CNVs and PRS-IQ on ASD risk and cognitive ability were largely additive, further suggesting that susceptibility for ASD is conferred independently from its effects on cognitive ability. Despite imparting mostly additive effects on ASD risk, rare CNVs and PRS-IQ showed opposing effects on core and associated features and developmental history among subjects with ASD. Discussion: Our findings suggest that cognitive ability itself may not be the factor driving the underlying liability for ASD conferred by these two classes of genomic variants. In other words, ASD risk and cognitive ability may be two distinct manifestations of CNVs and PRS-IQ. This study also highlights the challenge of understanding how genetic risk for ASD maps onto its dimensional traits.

Molecular epidemiology of viral haemorrhagic septicaemia virus (VHSV) in British Columbia, Canada, reveals transmission from wild to farmed fish.

Viral haemorrhagic septicaemia virus (VHSV) is a fish pathogen found throughout the Northern Hemisphere and is capable of infecting and causing mortality in numerous marine and freshwater hosts. In the coastal waters of British Columbia, Canada, the virus has been detected for 20 yr with many occurrences of mass mortalities among populations of Pacific herring Clupea pallasii (Valenciennes) and sardine Sardinops sagax as well as detections among cultured Atlantic Salmo salar and Chinook Oncorhynchus tshawytscha salmon. We compared nucleotide sequence of the full glycoprotein (G) gene coding region (1524 nt) of 63 VHSV isolates sampled during its recorded presence from 1993 to 2011 from 6 species and a total of 29 sites. Phylogenetic analysis showed that all isolates fell into sub-lineage IVa within the major VHSV genetic group IV. Of the 63 virus isolates, there were 42 unique sequences, each of which was ephemeral, being repeatedly detected at most only 1 yr after its initial detection. Multiple sequence types were revealed during single viral outbreak events, and genetic heterogeneity was observed within isolates from individual fish. Moreover, phylogenetic analysis revealed a close genetic linkage between VHSV isolates obtained from pelagic finfish species and farmed salmonids, providing evidence for virus transmission from wild to farmed fish.

Genomic analysis identifies risk factors in restless legs syndrome.

Restless legs syndrome (RLS) is a neurological condition that causes uncomfortable sensations in the legs and an irresistible urge to move them, typically during periods of rest. The genetic basis and pathophysiology of RLS are incompletely understood. Here, we present a whole-genome sequencing and genome-wide association meta-analysis of RLS cases (n = 9,851) and controls (n = 38,957) in three population-based biobanks (All of Us, Canadian Longitudinal Study on Aging, and CARTaGENE). Genome-wide association analysis identified nine independent risk loci, of which eight had been previously reported, and one was a novel risk locus (LMX1B, rs35196838, OR = 1.14, 95% CI = 1.09-1.19, p-value = 2.2 × 10-9). A genome-wide, gene-based common variant analysis identified GLO1 as an additional risk gene (p-value = 8.45 × 10-7). Furthermore, a transcriptome-wide association study also identified GLO1 and a previously unreported gene, ELFN1. A genetic correlation analysis revealed significant common variant overlaps between RLS and neuroticism (rg = 0.40, se = 0.08, p-value = 5.4 × 10-7), depression (rg = 0.35, se = 0.06, p-value = 2.17 × 10-8), and intelligence (rg = -0.20, se = 0.06, p-value = 4.0 × 10-4). Our study expands the understanding of the genetic architecture of RLS and highlights the contributions of common variants to this prevalent neurological disorder.

Copy-number variants and polygenic risk for intelligence confer risk for autism spectrum disorder irrespective of their effects on cognitive ability.

Rare copy number variants (CNVs) and polygenic risk for intelligence (PRS-IQ) both confer risk for autism spectrum disorder (ASD) but have opposing effects on cognitive ability. The field has struggled to disentangle the effects of these two classes of genomic variants on cognitive ability from their effects on ASD risk, in part because previous studies did not include controls with cognitive measures. We aim to investigate the impact of these genomic variants on ASD risk while adjusting for their known effects on cognitive ability. In a cohort of 8,426 subjects with ASD and 169,804 controls with cognitive assessments, we found that rare coding CNVs and PRS-IQ increased ASD risk, even after adjusting for their effects on cognitive ability. Bottom decile PRS-IQ and CNVs both decreased cognitive ability but had opposing effects on ASD risk. Models combining both classes of variants showed that the effects of rare CNVs and PRS-IQ on ASD risk and cognitive ability were largely additive, further suggesting that risk for ASD is conferred independently from its effects on cognitive ability. Despite imparting mostly additive effects on ASD risk, rare CNVs and PRS-IQ showed opposing effects on core and associated features and developmental history among subjects with ASD. Our findings suggest that cognitive ability itself may not be the factor driving the underlying risk for ASD conferred by these two classes of genomic variants. In other words, ASD risk and cognitive ability may be two distinct manifestations of CNVs and PRS-IQ. This study also highlights the challenge of understanding how genetic risk for ASD maps onto its dimensional traits.

Transcriptome-wide association study reveals increased neuronal FLT3 expression is associated with Tourette's syndrome.

Tourette's Syndrome (TS) is a neurodevelopmental disorder that is characterized by motor and phonic tics. A recent TS genome-wide association study (GWAS) identified a genome-wide significant locus. However, determining the biological mechanism of GWAS signals remains difficult. To characterize effects of expression quantitative trait loci (eQTLs) in TS and understand biological underpinnings of the disease. Here, we conduct a TS transcriptome-wide association study (TWAS) consisting of 4819 cases and 9488 controls. We demonstrate that increased expression of FLT3 in the dorsolateral prefrontal cortex (DLPFC) is associated with TS. We further show that there is global dysregulation of FLT3 across several brain regions and probabilistic causal fine-mapping of the TWAS signal prioritizes FLT3 with a posterior inclusion probability of 0.849. After, we proxy the expression with 100 lymphoblastoid cell lines, and demonstrate that TS cells has a 1.72 increased fold change compared to controls. A phenome-wide association study also points toward FLT3 having links with immune-related pathways such as monocyte count. We further identify several splicing events in MPHOSPH9, CSGALNACT2 and FIP1L1 associated with TS, which are also implicated in immune function. This analysis of expression and splicing begins to explore the biology of TS GWAS signals.

Questioning the Association of the STMN2 Dinucleotide Repeat With Amyotrophic Lateral Sclerosis.

Objectives: Recently, the number of dinucleotide CA repeats in an intron of the STMN2 gene was reported to be associated with an increased risk for amyotrophic lateral sclerosis (ALS). Therefore, we sought to replicate this observation in an independent group of ALS patients and a much larger control group. Methods: Here, we used whole-genome sequencing and tested the STMN2 CA repeat in a case-control cohort of the European genetic background and in genomes from various populations in the gnomAD cohort to attempt to replicate this proposed association. Results: We find that repeats well above the previously reported pathogenic threshold of 19 are commonly observed in unaffected individuals across different populations. Furthermore, we did not observe an association between longer STMN2 CA repeats and ALS phenotype. Discussion: In summary, our results do not support a role of STMN2 CA repeats toward ALS risk. As TDP-43 aggregation is central to ALS pathogenesis, lowered expression of STMN2 could be used as a biomarker for ALS. Therefore, a variant associated both with the risk for ALS and the level of STMN2 expression would be clinically useful. However, for a variant to be actionable, it must be strongly replicated in independent cohorts and exceed the rigorous statistical thresholds applied.

Transcriptomic effects of propranolol and primidone converge on molecular pathways relevant to essential tremor.

Essential tremor (ET) is one of the most common movement disorders, affecting nearly 5% of individuals over 65 years old. Despite this, few genetic risk loci for ET have been identified. Recent advances in pharmacogenomics have previously been useful to identify disease related molecular targets. Notably, gene expression has proven to be quite successful for the inference of drug response in cell models. We sought to leverage this approach in the context of ET where many patients are responsive to two drugs: propranolol and primidone. In this study, cerebellar DAOY and neural progenitor cells were treated for 5 days with clinical concentrations of propranolol and primidone, after which RNA-sequencing was used to identify convergent differentially expressed genes across treatments. Propranolol was found to affect the expression of genes previously associated with ET and other movement disorders such as TRAPPC11. Pathway enrichment analysis of these convergent drug-targeted genes identified multiple terms related to calcium signaling, endosomal sorting, axon guidance, and neuronal morphology. Furthermore, genes targeted by ET drugs were enriched within cell types having high expression of ET-related genes in both cortical and cerebellar tissues. Altogether, our results highlight potential cellular and molecular mechanisms associated with tremor reduction and identify relevant genetic biomarkers for drug-responsiveness in ET.