Publications

Alpha-synuclein is the main protein component of Lewy bodies, the pathological hallmark of Parkinson's disease. However, genetic modifiers of cerebrospinal fluid (CSF) alpha-synuclein levels remain unknown. The use of CSF levels of amyloid beta1-42, total tau, and phosphorylated tau181 as quantitative traits in genetic studies have provided novel insights into Alzheimer's disease pathophysiology. A systematic study of the genomic architecture of CSF biomarkers in Parkinson's disease has not yet been conducted. Here, genome-wide association studies of CSF biomarker levels in a cohort of individuals with Parkinson's disease and controls (N = 1960) were performed. PD cases exhibited significantly lower CSF biomarker levels compared to controls. A SNP, proxy for APOE ε4, was associated with CSF amyloid beta1-42 levels (effect = - 0.5, p = 9.2 × 10-19). No genome-wide loci associated with CSF alpha-synuclein, total tau, or phosphorylated tau181 levels were identified in PD cohorts. Polygenic risk score constructed using the latest Parkinson's disease risk meta-analysis were associated with Parkinson's disease status (p = 0.035) and the genomic architecture of CSF amyloid beta1-42 (R2 = 2.29%; p = 2.5 × 10-11). Individuals with higher polygenic risk scores for PD risk presented with lower CSF amyloid beta1-42 levels (p = 7.3 × 10-04). Two-sample Mendelian Randomization revealed that CSF amyloid beta1-42 plays a role in Parkinson's disease (p = 1.4 × 10-05) and age at onset (p = 7.6 × 10-06), an effect mainly mediated by variants in the APOE locus. In a subset of PD samples, the APOE ε4 allele was associated with significantly lower levels of CSF amyloid beta1-42 (p = 3.8 × 10-06), higher mean cortical binding potentials (p = 5.8 × 10-08), and higher Braak amyloid beta score (p = 4.4 × 10-04). Together these results from high-throughput and hypothesis-free approaches converge on a genetic link between Parkinson's disease, CSF amyloid beta1-42, and APOE.

Apart from amyloid β deposition and tau neurofibrillary tangles, Alzheimer's disease (AD) is a neurodegenerative disorder characterized by neuronal loss and astrocytosis in the cerebral cortex. The goal of this study is to investigate genetic factors associated with the neuronal proportion in health and disease. To identify cell-autonomous genetic variants associated with neuronal proportion in cortical tissues, we inferred cellular population structure from bulk RNA-Seq derived from 1536 individuals. We identified the variant rs1990621 located in the TMEM106B gene region as significantly associated with neuronal proportion (p value = 6.40 × 10-07) and replicated this finding in an independent dataset (p value = 7.41 × 10-04) surpassing the genome-wide threshold in the meta-analysis (p value = 9.42 × 10-09). This variant is in high LD with the TMEM106B non-synonymous variant p.T185S (rs3173615; r2 = 0.98) which was previously identified as a protective variant for frontotemporal lobar degeneration (FTLD). We stratified the samples by disease status, and discovered that this variant modulates neuronal proportion not only in AD cases, but also several neurodegenerative diseases and in elderly cognitively healthy controls. Furthermore, we did not find a significant association in younger controls or schizophrenia patients, suggesting that this variant might increase neuronal survival or confer resilience to the neurodegenerative process. The single variant and gene-based analyses also identified an overall genetic association between neuronal proportion, AD and FTLD risk. These results suggest that common pathways are implicated in these neurodegenerative diseases, that implicate neuronal survival. In summary, we identified a protective variant in the TMEM106B gene that may have a neuronal protection effect against general aging, independent of disease status, which could help elucidate the relationship between aging and neuronal survival in the presence or absence of neurodegenerative disorders. Our findings suggest that TMEM106B could be a potential target for neuronal protection therapies to ameliorate cognitive and functional deficits.

Epidemiologic studies have reported inconsistent results regarding an association between Parkinson disease (PD) and cutaneous melanoma (melanoma). Identifying shared genetic architecture between these diseases can support epidemiologic findings and identify common risk genes and biological pathways. Here, we apply polygenic, linkage disequilibrium-informed methods to the largest available case-control, genome-wide association study summary statistic data for melanoma and PD. We identify positive and significant genetic correlation (correlation: 0.17, 95% CI 0.10-0.24; P = 4.09 × 10-06) between melanoma and PD. We further demonstrate melanoma and PD-inferred gene expression to overlap across tissues (correlation: 0.14, 95% CI 0.06 to 0.22; P = 7.87 × 10-04) and highlight seven genes including PIEZO1, TRAPPC2L, and SOX6 as potential mediators of the genetic correlation between melanoma and PD. These findings demonstrate specific, shared genetic architecture between PD and melanoma that manifests at the level of gene expression.

Apolipoprotein E (APOE) ε4 genotype is associated with increased risk of dementia in Parkinson's disease (PD), but the mechanism is not clear, because patients often have a mixture of α-synuclein (αSyn), amyloid-β (Aβ), and tau pathologies. APOE ε4 exacerbates brain Aβ pathology, as well as tau pathology, but it is not clear whether APOE genotype independently regulates αSyn pathology. In this study, we generated A53T αSyn transgenic mice (A53T) on Apoe knockout (A53T/EKO) or human APOE knockin backgrounds (A53T/E2, E3, and E4). At 12 months of age, A53T/E4 mice accumulated higher amounts of brainstem detergent-insoluble phosphorylated αSyn compared to A53T/EKO and A53T/E3; detergent-insoluble αSyn in A53T/E2 mice was undetectable. By immunohistochemistry, A53T/E4 mice displayed a higher burden of phosphorylated αSyn and reactive gliosis compared to A53T/E2 mice. A53T/E2 mice exhibited increased survival and improved motor performance compared to other APOE genotypes. In a complementary model of αSyn spreading, striatal injection of αSyn preformed fibrils induced greater accumulation of αSyn pathology in the substantia nigra of A53T/E4 mice compared to A53T/E2 and A53T/EKO mice. In two separate cohorts of human patients with PD, APOE ε4/ε4 individuals showed the fastest rate of cognitive decline over time. Our results demonstrate that APOE genotype directly regulates αSyn pathology independent of its established effects on Aβ and tau, corroborate the finding that APOE ε4 exacerbates pathology, and suggest that APOE ε2 may protect against αSyn aggregation and neurodegeneration in synucleinopathies.

BACKGROUND: Low frequency coding variants in TREM2 are associated with Alzheimer disease (AD) risk and cerebrospinal fluid (CSF) TREM2 protein levels are different between AD cases and controls. Similarly, TREM2 risk variant carriers also exhibit differential CSF TREM2 levels. TREM2 has three different alternative transcripts, but most of the functional studies only model the longest transcript. No studies have analyzed TREM2 expression levels or alternative splicing in brains from AD and cognitively normal individuals. We wanted to determine whether there was differential expression of TREM2 in sporadic-AD cases versus AD-TREM2 carriers vs sex- and aged-matched normal controls; and if this differential expression was due to a particular TREM2 transcript.

METHODS: We analyzed RNA-Seq data from parietal lobe brain tissue from AD cases with TREM2 variants (n = 33), AD cases (n = 195) and healthy controls (n = 118), from three independent datasets using Kallisto and the R package tximport to determine the read count for each transcript and quantified transcript abundance as transcripts per million.

RESULTS: The three TREM2 transcripts were expressed in brain cortex in the three datasets. We demonstrate for the first time that the transcript that lacks the transmembrane domain and encodes a soluble form of TREM2 (sTREM2) has an expression level around 60% of the canonical transcript, suggesting that around 25% of the sTREM2 protein levels could be explained by this transcript. We did not observe a difference in the overall TREM2 expression level between cases and controls. However, the isoform which lacks the 5' exon, but includes the transmembrane domain, was significantly lower in TREM2- p.R62H carriers than in AD cases (p = 0.007).

CONCLUSION: Using bulk RNA-Seq data from three different cohorts, we were able to quantify the expression level of the three TREM2 transcripts, demonstrating: (1) all three transcripts of them are highly expressed in the human cortex, (2) that up to 25% of the sTREM2 may be due to the expression of a specific isoform and not TREM2 cleavage; and (3) that TREM2 risk variants do not affect expression levels, suggesting that the effect of the TREM2 variants on CSF levels occurs at post-transcriptional level.

BACKGROUND: Alzheimer's disease (AD) is the most common form of dementia. This neurodegenerative disorder is associated with neuronal death and gliosis heavily impacting the cerebral cortex. AD has a substantial but heterogeneous genetic component, presenting both Mendelian and complex genetic architectures. Using bulk RNA-seq from the parietal lobes and deconvolution methods, we previously reported that brains exhibiting different AD genetic architecture exhibit different cellular proportions. Here, we sought to directly investigate AD brain changes in cell proportion and gene expression using single-cell resolution.

METHODS: We generated unsorted single-nuclei RNA sequencing data from brain tissue. We leveraged the tissue donated from a carrier of a Mendelian genetic mutation, PSEN1 p.A79V, and two family members who suffer from sporadic AD, but do not carry any autosomal mutations. We evaluated alternative alignment approaches to maximize the titer of reads, genes, and cells with high quality. In addition, we employed distinct clustering strategies to determine the best approach to identify cell clusters that reveal neuronal and glial cell types and avoid artifacts such as sample and batch effects. We propose an approach to cluster cells that reduces biases and enable further analyses.

RESULTS: We identified distinct types of neurons, both excitatory and inhibitory, and glial cells, including astrocytes, oligodendrocytes, and microglia, among others. In particular, we identified a reduced proportion of excitatory neurons in the Mendelian mutation carrier, but a similar distribution of inhibitory neurons. Furthermore, we investigated whether single-nuclei RNA-seq from the human brains recapitulate the expression profile of disease-associated microglia (DAM) discovered in mouse models. We also determined that when analyzing human single-nuclei data, it is critical to control for biases introduced by donor-specific expression profiles.

CONCLUSION: We propose a collection of best practices to generate a highly detailed molecular cell atlas of highly informative frozen tissue stored in brain banks. Importantly, we have developed a new web application to make this unique single-nuclei molecular atlas publicly available.

Infantile globoid cell leukodystrophy (GLD, Krabbe disease) is a fatal demyelinating disorder caused by a deficiency in the lysosomal enzyme galactosylceramidase (GALC). GALC deficiency leads to the accumulation of the cytotoxic glycolipid, galactosylsphingosine (psychosine). Complementary evidence suggested that psychosine is synthesized via an anabolic pathway. Here, we show instead that psychosine is generated catabolically through the deacylation of galactosylceramide by acid ceramidase (ACDase). This reaction uncouples GALC deficiency from psychosine accumulation, allowing us to test the long-standing "psychosine hypothesis." We demonstrate that genetic loss of ACDase activity (Farber disease) in the GALC-deficient mouse model of human GLD (twitcher) eliminates psychosine accumulation and cures GLD. These data suggest that ACDase could be a target for substrate reduction therapy (SRT) in Krabbe patients. We show that pharmacological inhibition of ACDase activity with carmofur significantly decreases psychosine accumulation in cells from a Krabbe patient and prolongs the life span of the twitcher (Twi) mouse. Previous SRT experiments in the Twi mouse utilized l-cycloserine, which inhibits an enzyme several steps upstream of psychosine synthesis, thus altering the balance of other important lipids. Drugs that directly inhibit ACDase may have a more acceptable safety profile due to their mechanistic proximity to psychosine biogenesis. In total, these data clarify our understanding of psychosine synthesis, confirm the long-held psychosine hypothesis, and provide the impetus to discover safe and effective inhibitors of ACDase to treat Krabbe disease.