Evolution of immune genes is associated with the Black Death.

Infectious diseases are among the strongest selective pressures driving human evolution1,2. This includes the single greatest mortality event in recorded history, the first outbreak of the second pandemic of plague, commonly called the Black Death, which was caused by the bacterium Yersinia pestis3. This pandemic devastated Afro-Eurasia, killing up to 30-50% of the population4. To identify loci that may have been under selection during the Black Death, we characterized genetic variation around immune-related genes from 206 ancient DNA extracts, stemming from two different European populations before, during and after the Black Death. Immune loci are strongly enriched for highly differentiated sites relative to a set of non-immune loci, suggesting positive selection. We identify 245 variants that are highly differentiated within the London dataset, four of which were replicated in an independent cohort from Denmark, and represent the strongest candidates for positive selection. The selected allele for one of these variants, rs2549794, is associated with the production of a full-length (versus truncated) ERAP2 transcript, variation in cytokine response to Y. pestis and increased ability to control intracellular Y. pestis in macrophages. Finally, we show that protective variants overlap with alleles that are today associated with increased susceptibility to autoimmune diseases, providing empirical evidence for the role played by past pandemics in shaping present-day susceptibility to disease.

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.

TCEAL1 loss-of-function results in an X-linked dominant neurodevelopmental syndrome and drives the neurological disease trait in Xq22.2 deletions.

An Xq22.2 region upstream of PLP1 has been proposed to underly a neurological disease trait when deleted in 46,XX females. Deletion mapping revealed that heterozygous deletions encompassing the smallest region of overlap (SRO) spanning six Xq22.2 genes (BEX3, RAB40A, TCEAL4, TCEAL3, TCEAL1, and MORF4L2) associate with an early-onset neurological disease trait (EONDT) consisting of hypotonia, intellectual disability, neurobehavioral abnormalities, and dysmorphic facial features. None of the genes within the SRO have been associated with monogenic disease in OMIM. Through local and international collaborations facilitated by GeneMatcher and Matchmaker Exchange, we have identified and herein report seven de novo variants involving TCEAL1 in seven unrelated families: three hemizygous truncating alleles; one hemizygous missense allele; one heterozygous TCEAL1 full gene deletion; one heterozygous contiguous deletion of TCEAL1, TCEAL3, and TCEAL4; and one heterozygous frameshift variant allele. Variants were identified through exome or genome sequencing with trio analysis or through chromosomal microarray. Comparison with previously reported Xq22 deletions encompassing TCEAL1 identified a more-defined syndrome consisting of hypotonia, abnormal gait, developmental delay/intellectual disability especially affecting expressive language, autistic-like behavior, and mildly dysmorphic facial features. Additional features include strabismus, refractive errors, variable nystagmus, gastroesophageal reflux, constipation, dysmotility, recurrent infections, seizures, and structural brain anomalies. An additional maternally inherited hemizygous missense allele of uncertain significance was identified in a male with hypertonia and spasticity without syndromic features. These data provide evidence that TCEAL1 loss of function causes a neurological rare disease trait involving significant neurological impairment with features overlapping the EONDT phenotype in females with the Xq22 deletion.

Genomic Analysis Identifies Risk Factors in Restless Legs Syndrome.

OBJECTIVE: 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. We sought to identify additional novel genetic risk factors associated with RLS susceptibility. METHODS: We performed a whole-genome sequencing and genome-wide association meta-analysis of RLS cases (n = 9,851) and controls (n = 38,957) in 3 population-based biobanks (All of Us, Canadian Longitudinal Study on Aging, and CARTaGENE). RESULTS: Genome-wide association analysis identified 9 independent risk loci, of which 8 had been previously reported, and 1 was a novel risk locus (LMX1B, rs35196838, OR 1.14, 95% CI 1.09-1.19, p value = 2.2 × 10-9). 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). INTERPRETATION: Our study expands the understanding of the genetic architecture of RLS, and highlights the contributions of common variants to this prevalent neurological disorder. ANN NEUROL 2024.

RABENOSYN separation-of-function mutations uncouple endosomal recycling from lysosomal degradation, causing a distinct Mendelian disorder.

Rabenosyn (RBSN) is a conserved endosomal protein necessary for regulating internalized cargo. Here, we present clinical, genetic, cellular and biochemical evidence that two distinct RBSN missense variants are responsible for a novel Mendelian disorder consisting of progressive muscle weakness, facial dysmorphisms, ophthalmoplegia and intellectual disability. Using exome sequencing, we identified recessively acting germline alleles p.Arg180Gly and p.Gly183Arg, which are both situated in the FYVE domain of RBSN. We find that these variants abrogate binding to its cognate substrate phosphatidylinositol 3-phosphate (PI3P) and thus prevent its translocation to early endosomes. Although the endosomal recycling pathway was unaltered, mutant p.Gly183Arg patient fibroblasts show accumulation of cargo tagged for lysosomal degradation. Our results suggest that these variants are separation-of-function alleles, which cause a delay in endosomal maturation without affecting cargo recycling. We conclude that distinct germline mutations in RBSN cause non-overlapping phenotypes with specific and discrete endolysosomal cellular defects.

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.

The parkin V380L variant is a genetic modifier of Machado-Joseph disease with impact on mitophagy.

Machado-Joseph disease (MJD) is an autosomal dominant neurodegenerative spinocerebellar ataxia caused by a polyglutamine-coding CAG repeat expansion in the ATXN3 gene. While the CAG length correlates negatively with the age at onset, it accounts for approximately 50% of its variability only. Despite larger efforts in identifying contributing genetic factors, candidate genes with a robust and plausible impact on the molecular pathogenesis of MJD are scarce. Therefore, we analysed missense single nucleotide polymorphism variants in the PRKN gene encoding the Parkinson's disease-associated E3 ubiquitin ligase parkin, which is a well-described interaction partner of the MJD protein ataxin-3, a deubiquitinase. By performing a correlation analysis in the to-date largest MJD cohort of more than 900 individuals, we identified the V380L variant as a relevant factor, decreasing the age at onset by 3 years in homozygous carriers. Functional analysis in an MJD cell model demonstrated that parkin V380L did not modulate soluble or aggregate levels of ataxin-3 but reduced the interaction of the two proteins. Moreover, the presence of parkin V380L interfered with the execution of mitophagy-the autophagic removal of surplus or damaged mitochondria-thereby compromising cell viability. In summary, we identified the V380L variant in parkin as a genetic modifier of MJD, with negative repercussions on its molecular pathogenesis and disease age at onset.

Defining the Genetic Landscape of Congenital Mirror Movements in 80 Affected Individuals.

BACKGROUND: Congenital mirror movements (CMM) is a rare neurodevelopmental disorder characterized by involuntary movements from one side of the body that mirror voluntary movements on the opposite side. To date, five genes have been associated with CMM, namely DCC, RAD51, NTN1, ARHGEF7, and DNAL4. OBJECTIVE: The aim of this study is to characterize the genetic landscape of CMM in a large group of 80 affected individuals. METHODS: We screened 80 individuals with CMM from 43 families for pathogenic variants in CMM genes. In large CMM families, we tested for presence of pathogenic variants in multiple affected and unaffected individuals. In addition, we evaluated the impact of three missense DCC variants on binding between DCC and Netrin-1 in vitro. RESULTS: Causal pathogenic/likely pathogenic variants were found in 35% of probands overall, and 70% with familial CMM. The most common causal gene was DCC, responsible for 28% of CMM probands and 80% of solved cases. RAD51, NTN1, and ARHGEF7 were rare causes of CMM, responsible for 2% each. Penetrance of CMM in DCC pathogenic variant carriers was 68% and higher in males than females (74% vs. 54%). The three tested missense variants (p.Ile164Thr; p.Asn176Ser; and p.Arg1343His) bind Netrin-1 similarly to wild type DCC. CONCLUSIONS: A genetic etiology can be identified in one third of CMM individuals, with DCC being the most common gene involved. Two thirds of CMM individuals were unsolved, highlighting that CMM is genetically heterogeneous and other CMM genes are yet to be discovered. © 2024 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Characterizing proteomic and transcriptomic features of missense variants in amyotrophic lateral sclerosis genes.

Within recent years, there has been a growing number of genes associated with amyotrophic lateral sclerosis (ALS), resulting in an increasing number of novel variants, particularly missense variants, many of which are of unknown clinical significance. Here, we leverage the sequencing efforts of the ALS Knowledge Portal (3864 individuals with ALS and 7839 controls) and Project MinE ALS Sequencing Consortium (4366 individuals with ALS and 1832 controls) to perform proteomic and transcriptomic characterization of missense variants in 24 ALS-associated genes. The two sequencing datasets were interrogated for missense variants in the 24 genes, and variants were annotated with gnomAD minor allele frequencies, ClinVar pathogenicity classifications, protein sequence features including Uniprot functional site annotations, and PhosphoSitePlus post-translational modification site annotations, structural features from AlphaFold predicted monomeric 3D structures, and transcriptomic expression levels from Genotype-Tissue Expression. We then applied missense variant enrichment and gene-burden testing following binning of variation based on the selected proteomic and transcriptomic features to identify those most relevant to pathogenicity in ALS-associated genes. Using predicted human protein structures from AlphaFold, we determined that missense variants carried by individuals with ALS were significantly enriched in ß-sheets and α-helices, as well as in core, buried or moderately buried regions. At the same time, we identified that hydrophobic amino acid residues, compositionally biased protein regions and regions of interest are predominantly enriched in missense variants carried by individuals with ALS. Assessment of expression level based on transcriptomics also revealed enrichment of variants of high and medium expression across all tissues and within the brain. We further explored enriched features of interest using burden analyses and identified individual genes were indeed driving certain enrichment signals. A case study is presented for SOD1 to demonstrate proof-of-concept of how enriched features may aid in defining variant pathogenicity. Our results present proteomic and transcriptomic features that are important indicators of missense variant pathogenicity in ALS and are distinct from features associated with neurodevelopmental disorders.

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.

De Novo Pathogenic Variants in N-cadherin Cause a Syndromic Neurodevelopmental Disorder with Corpus Collosum, Axon, Cardiac, Ocular, and Genital Defects.

Cadherins constitute a family of transmembrane proteins that mediate calcium-dependent cell-cell adhesion. The extracellular domain of cadherins consists of extracellular cadherin (EC) domains, separated by calcium binding sites. The EC interacts with other cadherin molecules in cis and in trans to mechanically hold apposing cell surfaces together. CDH2 encodes N-cadherin, whose essential roles in neural development include neuronal migration and axon pathfinding. However, CDH2 has not yet been linked to a Mendelian neurodevelopmental disorder. Here, we report de novo heterozygous pathogenic variants (seven missense, two frameshift) in CDH2 in nine individuals with a syndromic neurodevelopmental disorder characterized by global developmental delay and/or intellectual disability, variable axon pathfinding defects (corpus callosum agenesis or hypoplasia, mirror movements, Duane anomaly), and ocular, cardiac, and genital anomalies. All seven missense variants (c.1057G>A [p.Asp353Asn]; c.1789G>A [p.Asp597Asn]; c.1789G>T [p.Asp597Tyr]; c.1802A>C [p.Asn601Thr]; c.1839C>G [p.Cys613Trp]; c.1880A>G [p.Asp627Gly]; c.2027A>G [p.Tyr676Cys]) result in substitution of highly conserved residues, and six of seven cluster within EC domains 4 and 5. Four of the substitutions affect the calcium-binding site in the EC4-EC5 interdomain. We show that cells expressing these variants in the EC4-EC5 domains have a defect in cell-cell adhesion; this defect includes impaired binding in trans with N-cadherin-WT expressed on apposing cells. The two frameshift variants (c.2563_2564delCT [p.Leu855Valfs∗4]; c.2564_2567dupTGTT [p.Leu856Phefs∗5]) are predicted to lead to a truncated cytoplasmic domain. Our study demonstrates that de novo heterozygous variants in CDH2 impair the adhesive activity of N-cadherin, resulting in a multisystemic developmental disorder, that could be named ACOG syndrome (agenesis of corpus callosum, axon pathfinding, cardiac, ocular, and genital defects).

Truncating Mutations in UBAP1 Cause Hereditary Spastic Paraplegia.

The diagnostic gap for rare neurodegenerative diseases is still considerable, despite continuous advances in gene identification. Many novel Mendelian genes have only been identified in a few families worldwide. Here we report the identification of an autosomal-dominant gene for hereditary spastic paraplegia (HSP) in 10 families that are of diverse geographic origin and whose affected members all carry unique truncating changes in a circumscript region of UBAP1 (ubiquitin-associated protein 1). HSP is a neurodegenerative disease characterized by progressive lower-limb spasticity and weakness, as well as frequent bladder dysfunction. At least 40% of affected persons are currently undiagnosed after exome sequencing. We identified pathological truncating variants in UBAP1 in affected persons from Iran, USA, Germany, Canada, Spain, and Bulgarian Roma. The genetic support ranges from linkage in the largest family (LOD = 8.3) to three confirmed de novo mutations. We show that mRNA in the fibroblasts of affected individuals escapes nonsense-mediated decay and thus leads to the expression of truncated proteins; in addition, concentrations of the full-length protein are reduced in comparison to those in controls. This suggests either a dominant-negative effect or haploinsufficiency. UBAP1 links endosomal trafficking to the ubiquitination machinery pathways that have been previously implicated in HSPs, and UBAP1 provides a bridge toward a more unified pathophysiology.

A new approach for rare variation collapsing on functional protein domains implicates specific genic regions in ALS.

Large-scale sequencing efforts in amyotrophic lateral sclerosis (ALS) have implicated novel genes using gene-based collapsing methods. However, pathogenic mutations may be concentrated in specific genic regions. To address this, we developed two collapsing strategies: One focuses rare variation collapsing on homology-based protein domains as the unit for collapsing, and the other is a gene-level approach that, unlike standard methods, leverages existing evidence of purifying selection against missense variation on said domains. The application of these two collapsing methods to 3093 ALS cases and 8186 controls of European ancestry, and also 3239 cases and 11,808 controls of diversified populations, pinpoints risk regions of ALS genes, including SOD1, NEK1, TARDBP, and FUS While not clearly implicating novel ALS genes, the new analyses not only pinpoint risk regions in known genes but also highlight candidate genes as well.

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.

Diagnostic Yield of Whole Exome Sequencing for Adults with Ataxia: a Brazilian Perspective.

Previous studies using whole exome sequencing (WES) have shown that a significant proportion of adult patients with undiagnosed ataxia in European and North American cohorts have a known genetic cause. Little is known about the diagnostic yield of WES in non-Caucasian ataxic populations. Herein, we used WES to investigate a Brazilian cohort of 76 adult patients with idiopathic ataxia previously screened for trinucleotide expansions in known ataxia genes. We collected clinical and radiological data from each patient. WES was performed following standard procedures. Only variants labeled as pathogenic or likely pathogenic according to American college of medical genetics and genomics (ACMG) criteria were retrieved. We determined the diagnostic yield of WES for the whole cohort and also for subgroups defined according to presence or not of pyramidal signs, peripheral neuropathy, and cerebellar atrophy. There were 41 women and 35 men. Mean age at testing was 48 years. Pyramidal signs, peripheral neuropathy, tremor, and cerebellar atrophy were found in 38.1%, 13.1%, 10.5%, and 68.3% of all subjects, respectively. Diagnostic yield of WES was 35.5%. Thirty-six distinct mutations were found in 20 different genes, determining the diagnosis of 18 autosomal recessive and 9 autosomal dominant ataxias. SACS and SPG7 were the most frequently found underlying genes. WES performed better in the subgroup with vs the subgroup without spasticity (p = 0.005). WES was diagnostic in 35.5% of cases of the Brazilian cohort of ataxia cases. These results have implications for diagnosis, genetic counseling and eventually treatment.

Targeted sequencing of Parkinson's disease loci genes highlights SYT11, FGF20 and other associations.

Genome-wide association studies (GWAS) have identified numerous loci associated with Parkinson's disease. The specific genes and variants that drive the associations within the vast majority of these loci are unknown. We aimed to perform a comprehensive analysis of selected genes to determine the potential role of rare and common genetic variants within these loci. We fully sequenced 32 genes from 25 loci previously associated with Parkinson's disease in 2657 patients and 3647 controls from three cohorts. Capture was done using molecular inversion probes targeting the exons, exon-intron boundaries and untranslated regions (UTRs) of the genes of interest, followed by sequencing. Quality control was performed to include only high-quality variants. We examined the role of rare variants (minor allele frequency < 0.01) using optimized sequence Kernel association tests. The association of common variants was estimated using regression models adjusted for age, sex and ethnicity as required in each cohort, followed by a meta-analysis. After Bonferroni correction, we identified a burden of rare variants in SYT11, FGF20 and GCH1 associated with Parkinson's disease. Nominal associations were identified in 21 additional genes. Previous reports suggested that the SYT11 GWAS association is driven by variants in the nearby GBA gene. However, the association of SYT11 was mainly driven by a rare 3' UTR variant (rs945006601) and was independent of GBA variants (P = 5.23 × 10-5 after exclusion of all GBA variant carriers). The association of FGF20 was driven by a rare 5' UTR variant (rs1034608171) located in the promoter region. The previously reported association of GCH1 with Parkinson's disease is driven by rare non-synonymous variants, some of which are known to cause dopamine-responsive dystonia. We also identified two LRRK2 variants, p.Arg793Met and p.Gln1353Lys, in 10 and eight controls, respectively, but not in patients. We identified common variants associated with Parkinson's disease in MAPT, TMEM175, BST1, SNCA and GPNMB, which are all in strong linkage disequilibrium with known GWAS hits in their respective loci. A common coding PM20D1 variant, p.Ile149Val, was nominally associated with reduced risk of Parkinson's disease (odds ratio 0.73, 95% confidence interval 0.60-0.89, P = 1.161 × 10-3). This variant is not in linkage disequilibrium with the top GWAS hits within this locus and may represent a novel association. These results further demonstrate the importance of fine mapping of GWAS loci, and suggest that SYT11, FGF20, and potentially PM20D1, BST1 and GPNMB should be considered for future studies as possible Parkinson's disease-related genes.

Genetic architecture and adaptations of Nunavik Inuit.

The Canadian Inuit have a distinct population background that may entail particular implications for the health of its individuals. However, the number of genetic studies examining this Inuit population is limited, and much remains to be discovered in regard to its genetic characteristics. In this study, we generated whole-exome sequences and genomewide genotypes for 170 Nunavik Inuit, a small and isolated founder population of Canadian Arctic indigenous people. Our study revealed the genetic background of Nunavik Inuit to be distinct from any known present-day population. The majority of Nunavik Inuit show little evidence of gene flow from European or present-day Native American peoples, and Inuit living around Hudson Bay are genetically distinct from those around Ungava Bay. We also inferred that Nunavik Inuit have a small effective population size of 3,000 and likely split from Greenlandic Inuit ∼10.5 kya. Nunavik Inuit went through a bottleneck at approximately the same time and might have admixed with a population related to the Paleo-Eskimos. Our study highlights population-specific genomic signatures in coding regions that show adaptations unique to Nunavik Inuit, particularly in pathways involving fatty acid metabolism and cellular adhesion (CPNE7, ICAM5, STAT2, and RAF1). Subsequent analyses in selection footprints and the risk of intracranial aneurysms (IAs) in Nunavik Inuit revealed an exonic variant under weak negative selection to be significantly associated with IA (rs77470587; P = 4.6 × 10-8).

GCH1 mutations in hereditary spastic paraplegia.

GCH1 mutations have been associated with dopa-responsive dystonia (DRD), Parkinson's disease (PD) and tetrahydrobiopterin (BH4 )-deficient hyperphenylalaninemia B. Recently, GCH1 mutations have been reported in five patients with hereditary spastic paraplegia (HSP). Here, we analyzed a total of 400 HSP patients (291 families) from different centers across Canada by whole exome sequencing (WES). Three patients with heterozygous GCH1 variants were identified: monozygotic twins with a p.(Ser77_Leu82del) variant, and a patient with a p.(Val205Glu) variant. The former variant is predicted to be likely pathogenic and the latter is pathogenic. The three patients presented with childhood-onset lower limb spasticity, hyperreflexia and abnormal plantar responses. One of the patients had diurnal fluctuations, and none had parkinsonism or dystonia. Phenotypic differences between the monozygotic twins were observed, who responded well to levodopa treatment. Pathway enrichment analysis suggested that GCH1 shares processes and pathways with other HSP-associated genes, and structural analysis of the variants indicated a disruptive effect. In conclusion, GCH1 mutations may cause HSP; therefore, we suggest a levodopa trial in HSP patients and including GCH1 in the screening panels of HSP genes. Clinical differences between monozygotic twins suggest that environmental factors, epigenetics, and stochasticity could play a role in the clinical presentation.

Genetic, Structural, and Functional Evidence Link TMEM175 to Synucleinopathies.

OBJECTIVE: The TMEM175/GAK/DGKQ locus is the 3rd strongest risk locus in genome-wide association studies of Parkinson disease (PD). We aimed to identify the specific disease-associated variants in this locus, and their potential implications. METHODS: Full sequencing of TMEM175/GAK/DGKQ followed by genotyping of specific associated variants was performed in PD (n = 1,575) and rapid eye movement sleep behavior disorder (RBD) patients (n = 533) and in controls (n = 1,583). Adjusted regression models and a meta-analysis were performed. Association between variants and glucocerebrosidase (GCase) activity was analyzed in 715 individuals with available data. Homology modeling, molecular dynamics simulations, and lysosomal localization experiments were performed on TMEM175 variants to determine their potential effects on structure and function. RESULTS: Two coding variants, TMEM175 p.M393T (odds ratio [OR] = 1.37, p = 0.0003) and p.Q65P (OR = 0.72, p = 0.005), were associated with PD, and p.M393T was also associated with RBD (OR = 1.59, p = 0.001). TMEM175 p.M393T was associated with reduced GCase activity. Homology modeling and normal mode analysis demonstrated that TMEM175 p.M393T creates a polar side-chain in the hydrophobic core of the transmembrane, which could destabilize the domain and thus impair either its assembly, maturation, or trafficking. Molecular dynamics simulations demonstrated that the p.Q65P variant may increase stability and ion conductance of the transmembrane protein, and lysosomal localization was not affected by these variants. INTERPRETATION: Coding variants in TMEM175 are likely to be responsible for the association in the TMEM175/GAK/DGKQ locus, which could be mediated by affecting GCase activity. ANN NEUROL 2020;87:139-153.