Scalable Nanopore sequencing of human genomes provides a comprehensive view of haplotype-resolved variation and methylation.

Long-read sequencing technologies substantially overcome the limitations of short-reads but have not been considered as a feasible replacement for population-scale projects, being a combination of too expensive, not scalable enough or too error-prone. Here we develop an efficient and scalable wet lab and computational protocol, Napu, for Oxford Nanopore Technologies long-read sequencing that seeks to address those limitations. We applied our protocol to cell lines and brain tissue samples as part of a pilot project for the National Institutes of Health Center for Alzheimer's and Related Dementias. Using a single PromethION flow cell, we can detect single nucleotide polymorphisms with F1-score comparable to Illumina short-read sequencing. Small indel calling remains difficult within homopolymers and tandem repeats, but achieves good concordance to Illumina indel calls elsewhere. Further, we can discover structural variants with F1-score on par with state-of-the-art de novo assembly methods. Our protocol phases small and structural variants at megabase scales and produces highly accurate, haplotype-specific methylation calls.

Genome-Wide Analysis of Structural Variants in Parkinson Disease.

OBJECTIVE: Identification of genetic risk factors for Parkinson disease (PD) has to date been primarily limited to the study of single nucleotide variants, which only represent a small fraction of the genetic variation in the human genome. Consequently, causal variants for most PD risk are not known. Here we focused on structural variants (SVs), which represent a major source of genetic variation in the human genome. We aimed to discover SVs associated with PD risk by performing the first large-scale characterization of SVs in PD. METHODS: We leveraged a recently developed computational pipeline to detect and genotype SVs from 7,772 Illumina short-read whole genome sequencing samples. Using this set of SV variants, we performed a genome-wide association study using 2,585 cases and 2,779 controls and identified SVs associated with PD risk. Furthermore, to validate the presence of these variants, we generated a subset of matched whole-genome long-read sequencing data. RESULTS: We genotyped and tested 3,154 common SVs, representing over 412 million nucleotides of previously uncatalogued genetic variation. Using long-read sequencing data, we validated the presence of three novel deletion SVs that are associated with risk of PD from our initial association analysis, including a 2 kb intronic deletion within the gene LRRN4. INTERPRETATION: We identified three SVs associated with genetic risk of PD. This study represents the most comprehensive assessment of the contribution of SVs to the genetic risk of PD to date. ANN NEUROL 2023;93:1012-1022.

THAP1 modulates oligodendrocyte maturation by regulating ECM degradation in lysosomes.

Mechanisms controlling myelination during central nervous system (CNS) maturation play a pivotal role in the development and refinement of CNS circuits. The transcription factor THAP1 is essential for timing the inception of myelination during CNS maturation through a cell-autonomous role in the oligodendrocyte lineage. Here, we demonstrate that THAP1 modulates the extracellular matrix (ECM) composition by regulating glycosaminoglycan (GAG) catabolism within oligodendrocyte progenitor cells (OPCs). Thap1-/- OPCs accumulate and secrete excess GAGs, inhibiting their maturation through an autoinhibitory mechanism. THAP1 controls GAG metabolism by binding to and regulating the GusB gene encoding ß-glucuronidase, a GAG-catabolic lysosomal enzyme. Applying GAG-degrading enzymes or overexpressing ß-glucuronidase rescues Thap1-/- OL maturation deficits in vitro and in vivo. Our studies establish lysosomal GAG catabolism within OPCs as a critical mechanism regulating oligodendrocyte development.

An ANXA11 P93S variant dysregulates TDP­43 and causes corticobasal syndrome.

INTRODUCTION: Variants of uncertain significance (VUS) surged with affordable genetic testing, posing challenges for determining pathogenicity. We examine the pathogenicity of a novel VUS P93S in Annexin A11 (ANXA11) - an amyotrophic lateral sclerosis/frontotemporal dementia-associated gene - in a corticobasal syndrome kindred. Established ANXA11 mutations cause ANXA11 aggregation, altered lysosomal-RNA granule co-trafficking, and transactive response DNA binding protein of 43 kDa (TDP-43) mis-localization. METHODS: We described the clinical presentation and explored the phenotypic diversity of ANXA11 variants. P93S's effect on ANXA11 function and TDP-43 biology was characterized in induced pluripotent stem cell-derived neurons alongside multiomic neuronal and microglial profiling. RESULTS: ANXA11 mutations were linked to corticobasal syndrome cases. P93S led to decreased lysosome colocalization, neuritic RNA, and nuclear TDP-43 with cryptic exon expression. Multiomic microglial signatures implicated immune dysregulation and interferon signaling pathways. DISCUSSION: This study establishes ANXA11 P93S pathogenicity, broadens the phenotypic spectrum of ANXA11 mutations, underscores neuronal and microglial dysfunction in ANXA11 pathophysiology, and demonstrates the potential of cellular models to determine variant pathogenicity. HIGHLIGHTS: ANXA11 P93S is a pathogenic variant. Corticobasal syndrome is part of the ANXA11 phenotypic spectrum. Hybridization chain reaction fluorescence in situ hybridization (HCR FISH) is a new tool for the detection of cryptic exons due to TDP-43-related loss of splicing regulation. Microglial ANXA11 and related immune pathways are important drivers of disease. Cellular models are powerful tools for adjudicating variants of uncertain significance.

Coding and Noncoding Variation in LRRK2 and Parkinson's Disease Risk.

BACKGROUND: The leucine-rich repeat kinase 2 (LRRK2) gene harbors both rare highly damaging missense variants (eg, p.G2019S) and common noncoding variants (eg, rs76904798) with lower effect sizes that are associated with Parkinson's disease (PD) risk. OBJECTIVES: This study aimed to investigate in a large meta-analysis whether the LRRK2 Genome-Wide Association Study (GWAS) signal represented by rs76904798 is independently associated with PD risk from LRRK2 coding variation and whether complex linkage disequilibrium structures with p.G2019S and the 5' noncoding haplotype account for the association of LRRK2 coding variants. METHODS: We performed a meta-analysis using imputed genotypes from 17,838 patients, 13,404 proxy patients, and 173,639 healthy controls of European ancestry. We excluded carriers of p.G2019S and/or rs76904798 to clarify the role of LRRK2 coding variation in mediating disease risk and excluded carriers of relatively rare LRRK2 coding variants to assess the independence of rs76904798. We also investigated the co-inheritance of LRRK2 coding variants with p.G2019S, rs76904798, and p.N2081D. RESULTS: LRRK2 rs76904798 remained significantly associated with PD after excluding the carriers of relatively rare LRRK2 coding variants. LRRK2 p.R1514Q and p.N2081D were frequently co-inherited with rs76904798, and the allele distribution of p.S1647T significantly changed among patients after removing rs76904798 carriers. CONCLUSIONS: These data suggest that the LRRK2 coding variants previously related to PD (p.N551K, p.R1398H, p.M1646T, and p.N2081D) do not drive the 5' noncoding GWAS signal. These data, however, do not preclude the independent association of the haplotype p.N551K-p.R1398H and p.M1646T with altered disease risk. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson Movement Disorder Society. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Parkinson-Associated SNCA Enhancer Variants Revealed by Open Chromatin in Mouse Dopamine Neurons.

The progressive loss of midbrain (MB) dopaminergic (DA) neurons defines the motor features of Parkinson disease (PD), and modulation of risk by common variants in PD has been well established through genome-wide association studies (GWASs). We acquired open chromatin signatures of purified embryonic mouse MB DA neurons because we anticipated that a fraction of PD-associated genetic variation might mediate the variants' effects within this neuronal population. Correlation with >2,300 putative enhancers assayed in mice revealed enrichment for MB cis-regulatory elements (CREs), and these data were reinforced by transgenic analyses of six additional sequences in zebrafish and mice. One CRE, within intron 4 of the familial PD gene SNCA, directed reporter expression in catecholaminergic neurons from transgenic mice and zebrafish. Sequencing of this CRE in 986 individuals with PD and 992 controls revealed two common variants associated with elevated PD risk. To assess potential mechanisms of action, we screened >16,000 proteins for DNA binding capacity and identified a subset whose binding is impacted by these enhancer variants. Additional genotyping across the SNCA locus identified a single PD-associated haplotype, containing the minor alleles of both of the aforementioned PD-risk variants. Our work posits a model for how common variation at SNCA might modulate PD risk and highlights the value of cell-context-dependent guided searches for functional non-coding variation.

Genetic modifiers of risk and age at onset in GBA associated Parkinson's disease and Lewy body dementia.

Parkinson's disease is a genetically complex disorder. Multiple genes have been shown to contribute to the risk of Parkinson's disease, and currently 90 independent risk variants have been identified by genome-wide association studies. Thus far, a number of genes (including SNCA, LRRK2, and GBA) have been shown to contain variability across a spectrum of frequency and effect, from rare, highly penetrant variants to common risk alleles with small effect sizes. Variants in GBA, encoding the enzyme glucocerebrosidase, are associated with Lewy body diseases such as Parkinson's disease and Lewy body dementia. These variants, which reduce or abolish enzymatic activity, confer a spectrum of disease risk, from 1.4- to >10-fold. An outstanding question in the field is what other genetic factors that influence GBA-associated risk for disease, and whether these overlap with known Parkinson's disease risk variants. Using multiple, large case-control datasets, totalling 217 165 individuals (22 757 Parkinson's disease cases, 13 431 Parkinson's disease proxy cases, 622 Lewy body dementia cases and 180 355 controls), we identified 1691 Parkinson's disease cases, 81 Lewy body dementia cases, 711 proxy cases and 7624 controls with a GBA variant (p.E326K, p.T369M or p.N370S). We performed a genome-wide association study and analysed the most recent Parkinson's disease-associated genetic risk score to detect genetic influences on GBA risk and age at onset. We attempted to replicate our findings in two independent datasets, including the personal genetics company 23andMe, Inc. and whole-genome sequencing data. Our analysis showed that the overall Parkinson's disease genetic risk score modifies risk for disease and decreases age at onset in carriers of GBA variants. Notably, this effect was consistent across all tested GBA risk variants. Dissecting this signal demonstrated that variants in close proximity to SNCA and CTSB (encoding cathepsin B) are the most significant contributors. Risk variants in the CTSB locus were identified to decrease mRNA expression of CTSB. Additional analyses suggest a possible genetic interaction between GBA and CTSB and GBA p.N370S induced pluripotent cell-derived neurons were shown to have decreased cathepsin B expression compared to controls. These data provide a genetic basis for modification of GBA-associated Parkinson's disease risk and age at onset, although the total contribution of common genetics variants is not large. We further demonstrate that common variability at genes implicated in lysosomal function exerts the largest effect on GBA associated risk for disease. Further, these results have implications for selection of GBA carriers for therapeutic interventions.

Parkinson's disease-associated mutations in the GTPase domain of LRRK2 impair its nucleotide-dependent conformational dynamics.

Mutation in leucine-rich repeat kinase 2 (LRRK2) is a common cause of familial Parkinson's disease (PD). Recently, we showed that a disease-associated mutation R1441H rendered the GTPase domain of LRRK2 catalytically less active and thereby trapping it in a more persistently "on" conformation. However, the mechanism involved and characteristics of this on conformation remained unknown. Here, we report that the Ras of complex protein (ROC) domain of LRRK2 exists in a dynamic dimer-monomer equilibrium that is oppositely driven by GDP and GTP binding. We also observed that the PD-associated mutations at residue 1441 impair this dynamic and shift the conformation of ROC to a GTP-bound-like monomeric conformation. Moreover, we show that residue Arg-1441 is critical for regulating the conformational dynamics of ROC. In summary, our results reveal that the PD-associated substitutions at Arg-1441 of LRRK2 alter monomer-dimer dynamics and thereby trap its GTPase domain in an activated state.

The Parkinson's Disease Genome-Wide Association Study Locus Browser.

BACKGROUND: Parkinson's disease (PD) is a neurodegenerative disease with an often complex component identifiable by genome-wide association studies. The most recent large-scale PD genome-wide association studies have identified more than 90 independent risk variants for PD risk and progression across more than 80 genomic regions. One major challenge in current genomics is the identification of the causal gene(s) and variant(s) at each genome-wide association study locus. The objective of the current study was to create a tool that would display data for relevant PD risk loci and provide guidance with the prioritization of causal genes and potential mechanisms at each locus. METHODS: We included all significant genome-wide signals from multiple recent PD genome-wide association studies including themost recent PD risk genome-wide association study, age-at-onset genome-wide association study, progression genome-wide association study, and Asian population PD risk genome-wide association study. We gathered data for all genes 1 Mb up and downstream of each variant to allow users to assess which gene(s) are most associated with the variant of interest based on a set of self-ranked criteria. Multiple databases were queried for each gene to collect additional causal data. RESULTS: We created a PD genome-wide association study browser tool (https://pdgenetics.shinyapps.io/GWASBrowser/) to assist the PD research community with the prioritization of genes for follow-up functional studies to identify potential therapeutic targets. CONCLUSIONS: Our PD genome-wide association study browser tool provides users with a useful method of identifying potential causal genes at all known PD risk loci from large-scale PD genome-wide association studies. We plan to update this tool with new relevant data as sample sizes increase and new PD risk loci are discovered. © 2020 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

The role of monogenic genes in idiopathic Parkinson's disease.

In the past two decades, mutations in multiple genes have been linked to autosomal dominant or recessive forms of monogenic Parkinson's disease (PD). Collectively, these monogenic (often familial) cases account for less than 5% of all PD, the majority being apparently sporadic cases. More recently, large-scale genome-wide association studies have identified over 40 loci that increase risk of PD. Importantly, there is overlap between monogenic and sporadic PD genes, particularly for the loci that contain the genes SNCA and LRRK2, which are mutated in monogenic dominant PD. There have also been reports of idiopathic PD cases with heterozygous variants in autosomal recessive genes suggesting that these mutations may increase risk of PD. These observations suggest that monogenic and idiopathic PD may have shared pathogenic mechanisms. Here, we focus mainly on the role of monogenic PD genes that represent pleomorphic risk loci for idiopathic PD. We also discuss the functional mechanisms that may play a role in increasing risk of disease in both monogenic and idiopathic forms.

Genetics in Parkinson disease: Mendelian versus non-Mendelian inheritance.

Parkinson's disease is a common, progressive neurodegenerative disorder, affecting 3% of those older than 75 years of age. Clinically, Parkinson's disease (PD) is associated with resting tremor, postural instability, rigidity, bradykinesia, and a good response to levodopa therapy. Over the last 15 years, numerous studies have confirmed that genetic factors contribute to the complex pathogenesis of PD. Highly penetrant mutations producing rare, monogenic forms of the disease have been discovered in singular genes such as SNCA, Parkin, DJ-1, PINK 1, LRRK2, and VPS35. Unique variants with incomplete penetrance in LRRK2 and GBA have been shown to be strong risk factors for PD in certain populations. Additionally, over 20 common variants with small effect sizes are now recognized to modulate the risk for PD. Investigating Mendelian forms of PD has provided precious insight into the pathophysiology that underlies the more common idiopathic form of disease; however, no treatment methodologies have developed. Furthermore, for identified common risk alleles, the functional basis underlying risk principally remains unknown. The challenge over the next decade will be to strengthen the findings delivered through genetic discovery by assessing the direct, biological consequences of risk variants in tandem with additional high-content, integrated datasets. This review discusses monogenic risk factors and mechanisms of Mendelian inheritance of Parkinson disease. Highly penetrant mutations in SNCA, Parkin, DJ-1, PINK 1, LRRK2 and VPS35 produce rare, monogenic forms of the disease, while unique variants within LRRK2 and GBA show incomplete penetrance and are strong risk factors for PD. Additionally, over 20 common variants with small effect sizes modulate disease risk. The challenge over the next decade is to strengthen genetic findings by assessing direct, biological consequences of risk variants in tandem with high-content, integrated datasets. This article is part of a special issue on Parkinson disease.

Systematic elucidation and in vivo validation of sequences enriched in hindbrain transcriptional control.

Illuminating the primary sequence encryption of enhancers is central to understanding the regulatory architecture of genomes. We have developed a machine learning approach to decipher motif patterns of hindbrain enhancers and identify 40,000 sequences in the human genome that we predict display regulatory control that includes the hindbrain. Consistent with their roles in hindbrain patterning, MEIS1, NKX6-1, as well as HOX and POU family binding motifs contributed strongly to this enhancer model. Predicted hindbrain enhancers are overrepresented at genes expressed in hindbrain and associated with nervous system development, and primarily reside in the areas of open chromatin. In addition, 77 (0.2%) of these predictions are identified as hindbrain enhancers on the VISTA Enhancer Browser, and 26,000 (60%) overlap enhancer marks (H3K4me1 or H3K27ac). To validate these putative hindbrain enhancers, we selected 55 elements distributed throughout our predictions and six low scoring controls for evaluation in a zebrafish transgenic assay. When assayed in mosaic transgenic embryos, 51/55 elements directed expression in the central nervous system. Furthermore, 30/34 (88%) predicted enhancers analyzed in stable zebrafish transgenic lines directed expression in the larval zebrafish hindbrain. Subsequent analysis of sequence fragments selected based upon motif clustering further confirmed the critical role of the motifs contributing to the classifier. Our results demonstrate the existence of a primary sequence code characteristic to hindbrain enhancers. This code can be accurately extracted using machine-learning approaches and applied successfully for de novo identification of hindbrain enhancers. This study represents a critical step toward the dissection of regulatory control in specific neuronal subtypes.

Integration of ChIP-seq and machine learning reveals enhancers and a predictive regulatory sequence vocabulary in melanocytes.

We take a comprehensive approach to the study of regulatory control of gene expression in melanocytes that proceeds from large-scale enhancer discovery facilitated by ChIP-seq; to rigorous validation in silico, in vitro, and in vivo; and finally to the use of machine learning to elucidate a regulatory vocabulary with genome-wide predictive power. We identify 2489 putative melanocyte enhancer loci in the mouse genome by ChIP-seq for EP300 and H3K4me1. We demonstrate that these putative enhancers are evolutionarily constrained, enriched for sequence motifs predicted to bind key melanocyte transcription factors, located near genes relevant to melanocyte biology, and capable of driving reporter gene expression in melanocytes in culture (86%; 43/50) and in transgenic zebrafish (70%; 7/10). Next, using the sequences of these putative enhancers as a training set for a supervised machine learning algorithm, we develop a vocabulary of 6-mers predictive of melanocyte enhancer function. Lastly, we demonstrate that this vocabulary has genome-wide predictive power in both the mouse and human genomes. This study provides deep insight into the regulation of gene expression in melanocytes and demonstrates a powerful approach to the investigation of regulatory sequences that can be applied to other cell types.

A SINE-VNTR-Alu at the LRIG2 locus is associated with proximal and distal gene expression in CRISPR and population models.

SINE-VNTR-Alu (SVA) retrotransposons represent mobile regulatory elements that have the potential to influence the surrounding genome when they insert into a locus. Evolutionarily recent mobilisation has resulted in loci in the human genome where a given retrotransposon might be observed to be present or absent, termed a retrotransposon insertion polymorphism (RIP). We previously observed that an SVA RIP ~ 2 kb upstream of LRIG2 on chromosome 1, the 'LRIG2 SVA', was associated with differences in local gene expression and methylation, and that the two were correlated. Here, we have used CRISPR-mediated deletion of the LRIG2 SVA in a cell line model to validate that presence of the retrotransposon is directly affecting local expression and provide evidence that is suggestive of a modest role for the SVA in modulating nearby methylation. Additionally, in leveraging an available Hi-C dataset we observed that the LRIG2 SVA was also involved in long-range chromatin interactions with a cluster of genes ~ 300 kb away, and that expression of these genes was to varying degrees associated with dosage of the SVA in both CRISPR cell line and population models. Altogether, these data support a regulatory role for SVAs in the modulation of gene expression, with the latter potentially involving chromatin looping, consistent with the model that RIPs may contribute to interpersonal differences in transcriptional networks.

Integration of genomic and functional approaches reveals enhancers at LMX1A and LMX1B.

LMX1A and LMX1B encode two closely related members of the LIM homeobox family of transcription factors. These genes play significant, and frequently overlapping, roles in the development of many structures in the nervous system, including the cerebellum, hindbrain, spinal cord roof plate, sensory systems and dopaminergic midbrain neurons. Little is known about the cis-acting regulatory elements (REs) that dictate their temporal and spatial expression or about the regulatory landscape surrounding them. The availability of comparative sequence data and the advent of genomic technologies such as ChIP-seq have revolutionized our capacity to identify regulatory sequences like enhancers. Despite this wealth of data, the vast majority of loci lack any significant in vivo functional exploration of their non-coding regions. We have completed a significant functional screen of conserved non-coding sequences (putative REs) scattered across these critical human loci, assaying the temporal and spatial control using zebrafish transgenesis. We first identify and describe the LMX1A paralogs lmx1a and lmx1a-like, comparing their expression during embryogenesis with that in mammals, along with lmx1ba and lmx1bb genes. Consistent with their prominent neuronal expression, 47/71 sequences selected within and flanking LMX1A and LMX1B exert spatial control of reporter expression in the central nervous system (CNS) of mosaic zebrafish embryos. Upon germline transmission, we identify CNS reporter expression in multiple independent founders for 22 constructs (LMX1A, n = 17; LMX1B, n = 5). The identified enhancers display significant overlap in their spatial control and represent only a fraction of the conserved non-coding sequences at these critical genes. Our data reveal the abundance of regulatory instruction located near these developmentally important genes.

Human brain single nucleus cell type enrichments in neurodegenerative diseases.

Background: Single-cell RNA sequencing has opened a window into clarifying the complex underpinnings of disease, particularly in quantifying the relevance of tissue- and cell-type-specific gene expression. Methods: To identify the cell types and genes important to therapeutic target development across the neurodegenerative disease spectrum, we leveraged genome-wide association studies, recent single-cell sequencing data, and bulk expression studies in a diverse series of brain region tissues. Results: We were able to identify significant immune-related cell types in the brain across three major neurodegenerative diseases: Alzheimer's disease, amyotrophic lateral sclerosis, and Parkinson's disease. Subsequently, putative roles of 30 fine-mapped loci implicating seven genes in multiple neurodegenerative diseases and their pathogenesis were identified. Conclusions: We have helped refine the genetic regions and cell types effected across multiple neurodegenerative diseases, helping focus future translational research efforts.

LRRK2: Genetic mechanisms vs genetic subtypes.

In 2004, the identification of pathogenic variants in the LRRK2 gene across several families with autosomal dominant late-onset Parkinson's disease (PD) revolutionized our understanding of the role of genetics in PD. Previous beliefs that genetics in PD was limited to rare early-onset or familial forms of the disease were quickly dispelled. Currently, we recognize LRRK2 p.G2019S as the most common genetic cause of both sporadic and familial PD, with more than 100,000 affected carriers across the globe. The frequency of LRRK2 p.G2019S is also highly variable across populations, with some regions of Asian or Latin America reporting close to 0%, contrasting to Ashkenazi Jews or North African Berbers reporting up to 13% and 40%, respectively. Patients with LRRK2 pathogenic variants are clinically and pathologically heterogeneous, highlighting the age-related variable penetrance that also characterizes LRRK2-related disease. Indeed, the majority of patients with LRRK2-related disease are characterized by a relatively mild Parkinsonism with less motor symptoms with variable presence of α-synuclein and/or tau aggregates, with pathologic pleomorphism widely described. At a functional cellular level, it is likely that pathogenic variants mediate a toxic gain-of-function of the LRRK2 protein resulting in increased kinase activity perhaps in a cell-specific manner; by contrast, some LRRK2 variants appear to be protective reducing PD risk by decreasing the kinase activity. Therefore, employing this information to define appropriate patient populations for clinical trials of targeted kinase LRRK2 inhibition strategies is very promising and demonstrates a potential future application for PD using precision medicine.

Human brain single nucleus cell type enrichments in neurodegenerative diseases.

Single cell RNA sequencing has opened a window into clarifying the complex underpinnings of disease, particularly in quantifying the relevance of tissue- and cell-type-specific gene expression. To identify the cell types and genes important to therapeutic target development across the neurodegenerative disease spectrum, we leveraged genome-wide association studies, recent single cell sequencing data, and bulk expression studies in a diverse series of brain region tissues. We were able to identify significant immune-related cell types in the brain across three major neurodegenerative diseases: Alzheimer's Disease, Amyotrophic Lateral Sclerosis, and Parkinson's Diseases. Subsequently, we identified the major role of 30 fine-mapped loci implicating seven genes in multiple neurodegenerative diseases and their pathogenesis.

Inhibition of p38α MAPK restores neuronal p38γ MAPK and ameliorates synaptic degeneration in a mouse model of DLB/PD.

Alterations in the p38 mitogen-activated protein kinases (MAPKs) play an important role in the pathogenesis of dementia with Lewy bodies (DLB) and Parkinson's disease (PD). Activation of the p38α MAPK isoform and mislocalization of the p38γ MAPK isoform are associated with neuroinflammation and synaptic degeneration in DLB and PD. Therefore, we hypothesized that p38α might be associated with neuronal p38γ distribution and synaptic dysfunction in these diseases. To test this hypothesis, we treated in vitro cellular and in vivo mouse models of DLB/PD with SKF-86002, a compound that attenuates inflammation by inhibiting p38α/ß, and then investigated the effects of this compound on p38γ and neurodegenerative pathology. We found that inhibition of p38α reduced neuroinflammation and ameliorated synaptic, neurodegenerative, and motor behavioral deficits in transgenic mice overexpressing human α-synuclein. Moreover, treatment with SKF-86002 promoted the redistribution of p38γ to synapses and reduced the accumulation of α-synuclein in mice overexpressing human α-synuclein. Supporting the potential value of targeting p38 in DLB/PD, we found that SKF-86002 promoted the redistribution of p38γ in neurons differentiated from iPS cells derived from patients with familial PD (carrying the A53T α-synuclein mutation) and healthy controls. Treatment with SKF-86002 ameliorated α-synuclein-induced neurodegeneration in these neurons only when microglia were pretreated with this compound. However, direct treatment of neurons with SKF-86002 did not affect α-synuclein-induced neurotoxicity, suggesting that SKF-86002 treatment inhibits α-synuclein-induced neurotoxicity mediated by microglia. These findings provide a mechanistic connection between p38α and p38γ as well as a rationale for targeting this pathway in DLB/PD.

An ANXA11 P93S variant dysregulates TDP-43 and causes corticobasal syndrome.

As genetic testing has become more accessible and affordable, variants of uncertain significance (VUS) are increasingly identified, and determining whether these variants play causal roles in disease is a major challenge. The known disease-associated Annexin A11 (ANXA11) mutations result in ANXA11 aggregation, alterations in lysosomal-RNA granule co-trafficking, and TDP-43 mis-localization and present as amyotrophic lateral sclerosis or frontotemporal dementia. We identified a novel VUS in ANXA11 (P93S) in a kindred with corticobasal syndrome and unique radiographic features that segregated with disease. We then queried neurodegenerative disorder clinic databases to identify the phenotypic spread of ANXA11 mutations. Multi-modal computational analysis of this variant was performed and the effect of this VUS on ANXA11 function and TDP-43 biology was characterized in iPSC-derived neurons. Single-cell sequencing and proteomic analysis of iPSC-derived neurons and microglia were used to determine the multiomic signature of this VUS. Mutations in ANXA11 were found in association with clinically diagnosed corticobasal syndrome, thereby establishing corticobasal syndrome as part of ANXA11 clinical spectrum. In iPSC-derived neurons expressing mutant ANXA11, we found decreased colocalization of lysosomes and decreased neuritic RNA as well as decreased nuclear TDP-43 and increased formation of cryptic exons compared to controls. Multiomic assessment of the P93S variant in iPSC-derived neurons and microglia indicates that the pathogenic omic signature in neurons is modest compared to microglia. Additionally, omic studies reveal that immune dysregulation and interferon signaling pathways in microglia are central to disease. Collectively, these findings identify a new pathogenic variant in ANXA11, expand the range of clinical syndromes caused by ANXA11 mutations, and implicate both neuronal and microglia dysfunction in ANXA11 pathophysiology. This work illustrates the potential for iPSC-derived cellular models to revolutionize the variant annotation process and provides a generalizable approach to determining causality of novel variants across genes.