Human brain single nucleus cell type enrichments in neurodegenerative diseases
==============================================================================

* Chelsea X. Alvarado
* Cory A. Weller
* Nicholas Johnson
* Hampton L. Leonard
* Andrew B. Singleton
* Xylena Reed
* Cornelis Blauewendraat
* Mike Nalls

## Abstract

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.

Keywords
*   Neurodegeneration
*   single cell expression
*   MAGMA
*   colocalization
*   Alzheimer’s Disease
*   Parkinson’s Disease
*   Lewy Body Dementia
*   Amyotrophic Lateral Sclerosis
*   Frontotemporal Dementia Lobar Degeneration
*   Progressive supranuclear palsy

## Introduction

Neurodegenerative diseases (NDD) encompass diseases characterized by progressive degeneration of cell types including neurons and glia in the central nervous system and/or peripheral nervous system. NDDs vary in their anatomic vulnerabilities and main affected neuropathologies, resulting in specific cell type dysregulation which may or may not be shared between NDDs.1,2 The complexity of biological processes and pathways involved in NDD pathogenesis has stymied progress in the understanding of disease and treatment development. While the expression of microglia in one cell type of the brain may be dysregulated, within the same disease state, other cell types may be regulated normally. Technology such as bulk RNA sequencing (RNA-seq) would be unlikely to identify these differences in expression since RNA-seq averages the expression of transcripts from all cells present in a tissue sample.

High throughput single-cell, or single-nucleus, mRNA sequencing (scRNA-seq, or snRNA-seq) technology provides a new window into understanding the functionally complex interactions between cell populations and across tissues. While bulk RNA-seq provides an assessment of expression averaged across a population of sampled cells, snRNA-seq allows for nuanced insight at the cell-type level. Instead of identifying gene set enrichment across an entire tissue, one can identify enrichment of gene expression partitioned by cell type, or changes in the composition of cell types themselves. uncertainty still limits our progress in better understanding the biological underpinnings of NDDs. In this report we leverage population-scale genome-wide association study (GWAS) data in conjunction with snRNA-seq in the brain to gain more specific insights into potential cell-type-specific mechanisms of risk within and across neurodegenerative diseases.

## Methods

### Single nucleus RNA-seq expression data

Data from adult human single nucleus brain tissues, consisting of the forebrain, midbrain, and hindbrain, was obtained from Siletti et al.3 The data was obtained by sequencing dissections from three healthy post-mortem donors. The subset of data used in this study consisted of snRNA-seq data. It included each of the 461 super clusters defined by Siletti et al. as being represented by the islands identified through a two-dimensional representation calculated by t-distributed stochastic neighbor embedding (t-SNE).3 Additionally, we calculated enrichment at the scRNA-seq data at the level of annotated cell type classes, using auto-annotation data provided by Siletti et al. in the supplementary materials of their manuscript.

### Genome-wide association study summary statistics

We included data from six NDD GWAS: Amyotrophic lateral sclerosis (ALS) from van Rheenen et al.4; Alzheimer’s disease (AD) from Bellenguez, et al.5; Frontotemporal lobar degeneration (FTLD) from Pottier, et al.6; Lewy body dementia (LBD) from Chia, et al.7; Parkinson’s disease (PD) from Nalls, et al.8; and Progressive supranuclear palsy (PSP) from Höglinger, et al.9 Data for ALS, AD, LBD, and PSP were obtained from GWAS Catalog ([https://www.ebi.ac.uk/gwas](https://www.ebi.ac.uk/gwas)) and data for FTLD and PD were obtained directly from the respective authors.

### MetaBrain eQTL summary statistics

We included expression quantitative trait loci (eQTL) summary statistics from MetaBrain, a large scale eQTL meta-analysis de Klein et al.10 For colocalization analysis, we included SNPs with a reported p-value no greater than 1×10−4 for the QTL component of the study.

### Cell Type Enrichment Analyses

We conducted cell type enrichment analyses using the R package MAGMA.11,12 Cell typing for each adult human brain snRNA-seq data and disease GWAS combination. Multiple testing correction (Bonferroni method) allows for efficient identification of enriched cell types. Required input data for conducting MAGMA analysis includes formatted GWAS summary statistics and a CellTypeDataset object (CTD).

Preprocessing of disease specific GWAS summary statistics and the creation of the CTD object, which holds cell specificity data, was performed using the R packages *MungeSumstats*13 and *ECWE*14 respectively. Quality control and munging of input data was conducted using methods available through the used R packages (see supplementary methods).

### Colocalization

We conducted Bayesian colocalization analysis using the R package *coloc*15 for all pairwise combinations of putatively significant autosomal NDD GWAS (p ≤ 5×10−8) and MetaBrain eQTLs (p ≤ 1×10−4). To account for the possibility of a single SNP influencing expression at multiple loci, we iterated over significant probes (see supplementary methods). We report associations with a posterior probability of at least 90%. A summary of data used and numbers of hits is included in **Supplementary Table S1**.

### Gene Expression Summarization

We summarized expression ranks for genes of interest within the single cell adult human brain transcriptome data set adult\_human\_20221007.loom from Siletti et al.3 Using custom R scripts we converted feature counts into TPM (transcripts per million). For a given sample, feature counts were divided by maximum nonredundant intron-removed exon lengths to correct for differences in gene length. Values were then multiplied by a sample-specific constant (106 / T, where T is the sum of length-normalized counts) such that the resulting unitless vector sums to one million. We extracted exon lengths based on annotations from the GTF file used to originally annotate the single cell data (gb_pri_annot.gtf). We calculated the expression percentile rank for genes of interest using the empirical cumulative distribution function and then calculated the mean and median expression percentile rank (EPR) value for each gene for each tested cell type. In order to ease interpretation, we binned the EPR values into 3 classes: off (EPR < 10), low (10 < EPR < 90) or high (EPR > 90).

### Online methods

Additional methods details can be found in the Supplementary Materials section as well as in the Data and Code Availability section. All analyses comprising this workflow are summarized in **Figure 1**.

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/07/01/2023.06.30.23292084/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2023/07/01/2023.06.30.23292084/F1)

Figure 1: Workflow diagram.
**NDD**: Neurodegenerative Disease. **GWAS**: Genome-Wide Association Study. **eQTL**: expression Quantitative Trait Loci. **ECWE**:Expression Weighted Cell Type Enrichment.14

## Results

### Cell type enrichments identified in AD, ALS, and PD

We identified significant cell type enrichments in the adult human brain for three out of six tested diseases: AD, ALS, and PD (**Table 1; Supplementary Table 2**). In general, we found that linear regression style enrichment analysis identified more significant enriched cell types than the top 10% enrichment style (nlinear = 30, nTop 10% = 5). MAGMA. Celltyping documentation does state that using the linear regression enrichment mode results in more significant results due to overlapping cell type signatures.

View this table:
[Table 1:](http://medrxiv.org/content/early/2023/07/01/2023.06.30.23292084/T1)

Table 1: Significant Celltype enrichments across AD, ALS, and PD.
Table provides information on each significantly enriched cell type after Bonferonni correction (p < 9.2×10−4). Results are provided for supercluster and class levels as well as differing enrichment modes. **CHRP ENDO**: Choroid Plexus and Endothelial. **ENDO**: Endothelial. **MAC**:Macrophage. **MGL**: Microglia. **MONO**: Monocyte. **NK**: Natural Killer cell. **OLIGO**: Oligodendrocyte. **OPC**: Oligodendrocyte precursor cell. **PER**: Pericyte. **TCELL**: T cell. **VSMC**: Vascular smooth muscle cell.

In PD, we ran MAGMA cell typing analysis on three variations/subsets of Nalls et al. meta-GWAS8: the complete meta-GWAS; meta-GWAS excluding 23andMe data; and meta-GWAS excluding 23andMe and UK Biobank data. In the first two variations tested, no cell types at the supercluster or class level reached significant enrichment after Bonferroni correction (p < 9.2×10−4). The only GWAS variation that identified significant cell type enrichments was when the full PD meta-GWAS was used in the cell typing analysis. We identified 14 enriched cell types when using linear regression analysis at the supercluster level. At the class level, we identified 10 enriched cell types (**Supplementary Table S3**). The top enriched cell type at the supercluster and class level are *thalamic excitatory* and *monocytes* (pthalmic excitatory = 1.467 ×10−5, *β*thalmic_excitatory = 3.278×10−3, pmonocytes = 1.314×10−6, *β*monocytes = 3.178×10−3), with monocytes being the most significantly enriched cell type in PD. Genes enriched in the monocyte cell type include *FCN1, CLEC12A, S100A4, TNFRSF1B, IFI30, LYZ, CYTIP, FGR, LILRB2, KYNU (***Supplementary Tables S4-5**).

In AD, we found microglia to be the only significantly enriched cell type at the supercluster level in both linear and top 10% analyses (p < 9.2×10−4 ; plinear = 4.193×10−8, *β*linear = 3.978×10−3, *β* pTop 10% = 1.195×10−6, *β*Top 10% = 0.1278) for risk loci. At the class level, four cell types were found to be significant after Bonferroni correction and in both analyses. Significantly enriched cell types identified include macrophages, monocytes, microglia, and natural killer cells (**Supplementary Table S6**). AD was the only disease to have significant enrichments when using the *Top 10%* enrichment analyses.

In ALS, we identified one significant cell type enrichment using linear regression enrichment analysis and at the class level annotations. *Monocytes* were the only cell type to reach significance (Bonferroni-corrected threshold of p < 9.2×10−4); pMonocytes = 8.516×10−4, *β*monocytes = 1.778×10−3; **Supplementary Table S7**). We did not detect significant cell type enrichments in FTLD, LBD, and PSP. No cell types at either the supercluster or class level reached significance after MAGMA-implemented Bonferroni correction (p < 9.2×10−4; **Supplementary Table S2, Supplementary Tables S8-10**).

### Colocalization

Across all diseases tested we fine mapped a total of 205 association signals at posterior probability > 90%. This included 89 unique genes identified as harboring putative causal associations. Of these 205 associations, 20 were centered around the *HLA* region, with colocalized omic associations in the cerebellum and cortex (the latter in multiple ancestry groups) suggesting extremely complicated risk of neuroinflammation in this part of the genome. Interestingly, *GRN* is fine-mapped using colocalized expression QTL (eQTL) signals in the cerebellum for AD and the cortex for PD suggesting related but potentially different mechanisms. The *TMEM175*/*GAK* region shows multiple colocalized signals across multiple diseases (PD and LBD). Of particular interest is that *GAK* and *TMEM175* both are fine-mapped to the same SNP (rs6599388) with eQTL effects in opposite directions in LBD. At the same time, regional signals for the gene *SLC26A1* were also fine-mapped for PD, LBD and AD with the effect in PD being detected in the spinal cord while the other disease QTLs were localized to the cortex.

Variants exhibiting a colocalization posterior probability ≥ 90% are summarized in **Figure 2** and **Table 2**. Extended data is available in **Supplementary Table S11**. 30 loci were fine-mapped to a single gene per disease by leveraging QTL data (see **Table 2** for details10,16). Of these, six are known druggable genes including *ABCA1, ADAM10, CD55, FGF7, OXGR1* and *POLE* (see **Supplementary Table S12**). Mining additional data on these druggable genes from the omicSynth16 database (**Supplementary Table S13**), we note that *ABCA1* has putative functional multi-omic associations with AD and PSP in blood and brain tissues. *ADAM10* has a similar pattern of functional inferences in AD. *CD55* is shown to have multiple significant functional inferences in brain tissues mediating risk for AD and LBD. However, *FGF7* does not display any significant functional inferences in any diseases from the database query. Methylation QTLs in blood connect PSP and FTD at *OXGR1* via functional inferences using SMR, while multiple brain, blood and nerve associations connect this gene with PD risk across both expression and methylation QTLs. A similar pattern of disease and tissue associations is seen for *POLE*, although there is no significant neural tissue association for PD and PSP is also connected via blood eQTLs to the same gene. Of the fine mapped loci, 22 (Supplemental Table S16) have also been nominated elsewhere as potential therapeutic targets with likely functional impacts on neurodegenerative disease risk in the context of methylation, expression, protein or chromatin QTLs detailed in the omicSynth web application.16

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/07/01/2023.06.30.23292084/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2023/07/01/2023.06.30.23292084/F2)

Figure 2: Graphical summaries of cell type enrichment and colocalization results.
**A**. Tile plot comparison of significant cell type enrichments, by disease at each level for both MAGMA.Celltyping analysis. A Bonferonni significance line is provided on the bar chart portion of the image. **B**. Comparison of significant cell type enrichments, with celltype dendrogram, by disease at the class level for both MAGMA.Celltyping analysis. **C**. Summary of fine-mapped loci per disease. AD = Alzheimer’s Disease, ALS= Amyotrophic Lateral Sclerosis, PD = Parkinson’s Disease.

View this table:
[Table 2:](http://medrxiv.org/content/early/2023/07/01/2023.06.30.23292084/T2)

Table 2: Summary of colocalized signals.
This summary details all associations suggesting a single gene associated with multiple diseases of interest or a single gene per locus at posterior probability above 90%.

### Cell type resolution of colocalized genes

We further evaluated gene expression from snRNA-seq used in enrichment analysis. We calculated the mean and median expression percentile rank (EPR) for each gene across cells corresponding to the nominated supercluster cell types identified in our enrichment analyses (**Table 1**) and compared the aggregate mean and median values against the nominated colocalized genes (**Supplementary Figure S2, Supplementary Tables S15-16**). In order to ease interpretation we binned the mean EPR values into three categories based on the mean EPR value for each gene-cell type combination: *off, low*, and *high* (see methods).

We identified 10 colocalized genes with *high* median EPR values out of the 14 tested genes. The gene *PAM* has 6 cell type combinations (CGE interneuron, Eccentric medium spiny neuron, MGE interneuron, Hippocampal dentate gyrus, Thalamic excitatory, and Mammillary body) classified as *high* (**Figure 3**). Overall, the Mamillary body cell type had the greatest count of eight *high* EPR genes.

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/07/01/2023.06.30.23292084/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2023/07/01/2023.06.30.23292084/F3)

Figure 3: scRNA-seq expression distribution across colocalized genes.
This graphical representation visualizes the scRNA-seq expression distribution across MAGMA nominated enriched supercluster cell types. Expression percentile rank (EPR) for each cell-gene combination was binned into *high expression* (dark blue, median expression within top 10% of all genes), *low expression* (green, middle 80%) or *off* (orange, bottom 10%).

## Discussion

We identified enriched cell types in various brain regions for three (AD, ALS, and PD) out of six tested NDDs utilizing snRNA-seq data of the 461 superclusters identified by Siletti et al.3 We did not detect significant cell type enrichments, at either the supercluster or class level, in FTLD, LBD, and PSP. We speculate that due to the smaller sample sizes of each GWAS we were unable to identify any significantly enriched cell types after applying the MAGMA-implemented Bonferroni correction (p < 9.2×10−4; **Supplementary Table S2, Supplementary Tables S8-10**).

All significantly enriched cell types had an associated positive beta for risk genes, and our identified significant cell type enrichments fall in line with current knowledge on cell types implicated in various NDD pathologies which we will discuss further. Broadly, our results highlight the importance of immune-related cell types in the pathology of varying NDDs. Additionally, our results add to the growing evidence that NDDs, such as AD, ALS, and PD, are highly related to autoimmune diseases or may even be autoimmune diseases of the brain.17–20

Monocytes were the only cell type significantly enriched across the three diseases AD, ALS, and PD. It is the most significantly enriched cell type in ALS and PD (pALS=8.51×10−4, pPD = 1.31×10−6) and the third most enriched cell in AD (pAD = 9.25×10−8). Monocytes are the precursor cells to dendritic cells (DCs) and macrophages which both play important roles in immune response, neuroinflammation, and neuroimmunological response.21–23 Of the genes enriched in the monocyte cell type (previously described in our results), *KYNU* is a gene of note, being part of the tryptophan metabolic pathway which has been found to play a role in Aβ formation.19

Microglia are another significantly enriched cell type identified in our analyses, though only significantly enriched in AD. Previous literature highlights the role that microglia play in neuronal loss and neuroinflammation as well as their involvement in immune response.24–26 Microglia are known to function similarly to DCs and macrophages, derived from monocytes.25 In AD literature, microglia is implicated as an affected cell type associated with neuroinflammation.25,27 Our results highlight microglia as the second most enriched cell type in AD after Bonferroni correction (pAD = 4.19×10−8). Microglia were found to be nominally significant in LBD and PD (pLBD(Linear) = 0.0119, pPD(Linear) = 9.25×10−3, pPD(Linear) =1.96×10−03, **Supplementary Tables 2, 3, 6**). Interestingly, while not shown as enriched in our ALS analyses, microglia have been shown to play an important role in the pathogenesis of ALS in part due to their activation being neurotoxic to motor neurons.25,27–30 It is possible that due to the limited sample size in the used snRNA-seq data, the analyses were unable to detect any significant cell type enrichment in microglia for ALS.

Common risk factors across neurodegenerative diseases have always been of interest to the basic science and therapeutic industries. The association at the *SLC26A1* locus is of particular interest across multiple NDDs, as it is strongly associated with *IDUA* protein level in QTL studies, with deficiencies in this protein causing severe lysosomal storage disorders.31,32 The *HLA* region is a complex locus from a structural genetic standpoint, but also in terms of general risk of neuroinflammation so not surprisingly resulting in a number of fine-mapped associations across this locus for multiple genes.5,33 *GRN* is a positive control for this colocalization analysis effort as previous efforts show increased genetic risk in the region coincides with decreased expression of the gene in PD, AD and ALS.34

Fine-mapping efforts localizing signals to single genes within risk loci help to identify novel therapeutic targets with known biological plausibility. *ABCA1* is known to be associated with Tangier disease, which is a rare autosomal recessive disorder with low plasma levels of high-density lipoprotein (HDL) causing peripheral neuropathy.35 *ADAM10* is implicated in the formation of amyloid plaques in the brain and the processing of APP.35,36 CD55 is known to interact with viruses and cause neuroinflammation potentially leading to increased AD risk.37 *COQ8A* mutations have been shown to cause coenzyme Q10 deficiency leading to autosomal recessive ataxia, cerebellar atrophy, and progressive movement disorders.38 *C9orf72*, known to be linked to ALS, acts as a positive control here and displays a unique ALS/FTD colocalization. *ADAM10* and *C9orf72* were also nominated as potential therapeutic drug targets for neurodegenerative disease via Mendelian randomization.16

Limitations to this study generally relate to availability of data in this context. First and foremost, there is a limited amount of multi-ancestry or non-European data available for the GWAS and single-cell or QTL resources used here. This potentially introduces bias into therapeutic development and precision medicine applications. Secondly, low sample sizes for single nucleus analyses (in terms of the number of humans involved) reduces our ability to generate eQTL databases as compared to coarse methods of bulk RNA sequencing at scale.

Here we provided insights that could potentially aid in therapeutic development for NDDs. On the macro-level, we have identified cell-type level enrichments associated with disease risk in multiple neurodegenerative diseases allowing biologists and drug developers to better focus their mechanistic and therapeutic research. On the micro-level, we have used eQTL colocalization methods to narrow down the large tracts of associated loci in GWAS to potentially functional variants, metaphorically going from a neighborhood to building level resolution on a map.

## Supporting information

Supplemental Materials [[supplements/292084_file02.pdf]](pending:yes)

Supplemental Tables [[supplements/292084_file03.xlsx]](pending:yes)

## Data Availability

All GWAS summary statistics used in this report are available at [https://nih-card-ndd-smr-home-syboky.streamlit.app/About](https://nih-card-ndd-smr-home-syboky.streamlit.app/About) Single nucleus (sn) RNA-seq expression data is available at [https://github.com/linnarsson-lab/adult-human-brain](https://github.com/linnarsson-lab/adult-human-brain) which is the repository associated with [https://doi.org/10.1101/2022.10.12.511898](https://doi.org/10.1101/2022.10.12.511898) eQTL data is available at the MetaBrain browser: [https://www.metabrain.nl/](https://www.metabrain.nl/) Code used for analysis is available at [https://github.com/NIH-CARD/NDD\_Single_Cell](https://github.com/NIH-CARD/NDD_Single_Cell)

[https://github.com/NIH-CARD/NDD\_Single\_Cell](https://github.com/NIH-CARD/NDD_Single_Cell) 

[https://nih-card-ndd-smr-home-syboky.streamlit.app/About](https://nih-card-ndd-smr-home-syboky.streamlit.app/About) 

[https://github.com/linnarsson-lab/adult-human-brain](https://github.com/linnarsson-lab/adult-human-brain) 

[https://www.metabrain.nl/](https://www.metabrain.nl/) 

## Data and Code availability

Code used to generate and process our data can be accessed at our Github Repository. All data for this project are publicly available via the original publications accessed by our team. Summaries of enrichment and colocalization data are also available to browse in our community target discovery and due diligence resource omicSynth web application.16

*   Received June 30, 2023.
*   Revision received June 30, 2023.
*   Accepted July 1, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available for use under a CC0 license

## References

1.  1.Dugger BN, Dickson DW. Pathology of Neurodegenerative Diseases.Cold Spring Harb Perspect Biol. 2017;9(7):a028035. doi:10.1101/cshperspect.a028035
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImNzaHBlcnNwZWN0IjtzOjU6InJlc2lkIjtzOjExOiI5LzcvYTAyODAzNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA3LzAxLzIwMjMuMDYuMzAuMjMyOTIwODQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

2.  2.Cuevas-Diaz Duran R, Gonzalez-Orozco JC, Velasco I, Wu JQ. Single-cell and single-nuclei RNA sequencing as powerful tools to decipher cellular heterogeneity and dysregulation in neurodegenerative diseases. Front Cell Dev Biol. 2022;10:884748. doi:10.3389/fcell.2022.884748
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fcell.2022.884748&link_type=DOI) 

3.  3.Siletti K, Hodge R, Mossi Albiach A, et al. Transcriptomic diversity of cell types across the adult human brain. Published online 2022.doi:10.1101/2022.10.12.511898
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMi4xMC4xMi41MTE4OTh2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA3LzAxLzIwMjMuMDYuMzAuMjMyOTIwODQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

4.  4.Van Rheenen W, Van Der Spek RAA, Bakker MK, et al. Common and rare variant association analyses in amyotrophic lateral sclerosis identify 15 risk loci with distinct genetic architectures and neuron-specific biology. Nat Genet. 2021;53(12):1636–1648. doi:10.1038/s41588-021-00973-1
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-021-00973-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34873335&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F01%2F2023.06.30.23292084.atom) 

5.  5.Bellenguez C, Küçükali F, Jansen IE, et al. New insights into the genetic etiology of Alzheimer’s disease and related dementias. Nat Genet. 2022;54(4):412–436. doi:10.1038/s41588-022-01024-z
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-022-01024-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35379992&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F01%2F2023.06.30.23292084.atom) 

6.  6.Pottier C, Ren Y, Perkerson RB, et al. Genome-wide analyses as part of the international FTLD-TDP whole-genome sequencing consortium reveals novel disease risk factors and increases support for immune dysfunction in FTLD. Acta Neuropathol. 2019;137(6):879–899. doi:10.1007/s00401-019-01962-9
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00401-019-01962-9&link_type=DOI) 

7.  7.Chia R, Sabir MS, Bandres-Ciga S, et al. Genome sequencing analysis identifies new loci associated with Lewy body dementia and provides insights into its genetic architecture. Nat Genet. 2021;53(3):294–303.doi:10.1038/s41588-021-00785-3
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-021-00785-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33589841&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F01%2F2023.06.30.23292084.atom) 

8.  8.Nalls MA, Blauwendraat C, Vallerga CL, et al. Identification of novel risk loci,causal insights,and heritable risk for Parkinson’s disease: a meta-analysis of genome-wide association studies. Lancet Neurol. 2019;18(12):1091–1102. doi:10.1016/s1474-4422(19)30320-5
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s1474-4422(19)30320-5&link_type=DOI) 

9.  9.Höglinger GU, Melhem NM, Dickson DW, et al. Identification of common variants influencing risk of the tauopathy progressive supranuclear palsy. Nat Genet. 2011;43(7):699–705. doi:10.1038/ng.859
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.859&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21685912&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F01%2F2023.06.30.23292084.atom) 

10. 10.de Klein N, Tsai EA, Vochteloo M, et al. Brain expression quantitative trait locus and network analyses reveal downstream effects and putative drivers for brain-related diseases. Nat Genet. 2023;55(3):377–388.doi:10.1038/s41588-023-01300-6
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-023-01300-6&link_type=DOI) 

11. 11.Skene NG, Bryois J, Bakken TE, et al. Genetic identification of brain cell types underlying schizophrenia. Nat Genet. 2018;50(6):825–833. doi:10.1038/s41588-018-0129-5
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0129-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29785013&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F01%2F2023.06.30.23292084.atom) 

12. 12.de Leeuw CA, Mooij JM, Heskes T, Posthuma D. MAGMA: generalized gene-set analysis of GWAS data. PLoS Comput Biol. 2015;11(4):e1004219. doi:10.1371/journal.pcbi.1004219
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1004219&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25885710&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F01%2F2023.06.30.23292084.atom) 

13. 13.Murphy AE, Schilder BM, Skene NG. MungeSumstats: a Bioconductor package for the standardization and quality control of many GWASsummary statistics. Bioinformatics. 2021;37(23):4593–4596.doi:10.1093/bioinformatics/btab665
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btab665&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34601555&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F01%2F2023.06.30.23292084.atom) 

14. 14.Skene NG, Grant SGN. Identification of Vulnerable Cell Types in Major Brain Disorders Using Single Cell Transcriptomes and Expression Weighted Cell Type Enrichment. Front Neurosci. 2016;10:16.doi:10.3389/fnins.2016.00016
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fnins.2016.00016&link_type=DOI) 

15. 15.Giambartolomei C, Vukcevic D, Schadt EE, et al. Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet. 2014;10(5):e1004383. doi:10.1371/journal.pgen.1004383
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1004383&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24830394&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F01%2F2023.06.30.23292084.atom) 

16. 16.Alvarado CX, Makarious MB, Vitale D, et al. omicSynth: an Open Multi-omic Community Resource for Identifying Druggable Targets across Neurodegenerative Diseases. medRxiv. Published online April 12,2023:2023.04.06.23288266.doi:10.1101/2023.04.06.23288266
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMy4wNC4wNi4yMzI4ODI2NnYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDcvMDEvMjAyMy4wNi4zMC4yMzI5MjA4NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

17. 17.Lyon MS, Wosiski-Kuhn M, Gillespie R, Caress J, Milligan C. Inflammation,Immunity,and amyotrophic lateral sclerosis: I. Etiology and pathology. Muscle Nerve. 2019;59(1):10–22. doi:10.1002/mus.26289
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mus.26289&link_type=DOI) 

18. 18.Bonam SR, Muller S. Parkinson’s disease is an autoimmune disease: Areappraisal. Autoimmun Rev. 2020;19(12):102684.doi:10.1016/j.autrev.2020.102684
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.autrev.2020.102684&link_type=DOI) 

19. 19.Meier-Stephenson FS, Meier-Stephenson VC, Carter MD, et al. Alzheimer’s disease as an autoimmune disorder of innate immunity endogenously modulated by tryptophan metabolites. Alzheimers Dement. 2022;8(1). doi:10.1002/trc2.12283
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/trc2.12283&link_type=DOI) 

20. 20.Li M, Wan J, Xu Z, Tang B. The association between Parkinson’s disease and autoimmune diseases: Asystematic review and meta-analysis. Front Immunol. 2023;14:1103053. doi:10.3389/fimmu.2023.1103053
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fimmu.2023.1103053&link_type=DOI) 

21. 21.Feng Y, Li L, Sun XH. Monocytes and Alzheimer’s disease. Neurosci Bull. 2011;27(2):115–122. doi:10.1007/s12264-011-1205-3
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12264-011-1205-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21441973&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F01%2F2023.06.30.23292084.atom) 

22. 22.Wlodarczyk A, Løbner M, Cédile O, Owens T. Comparison of microglia and infiltrating CD11c+ cells as antigen presenting cells for Tcell proliferation and cytokine response.JNeuroinflammation. 2014;11(1):57.doi:10.1186/1742-2094-11-57
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1742-2094-11-57&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24666681&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F01%2F2023.06.30.23292084.atom) 

23. 23.Mammana S, Fagone P, Cavalli E, et al. The Role of Macrophages in Neuroinflammatory and Neurodegenerative Pathways of Alzheimer’s Disease,Amyotrophic Lateral Sclerosis,and Multiple Sclerosis: Pathogenetic Cellular Effectors and Potential Therapeutic Targets. Int JMol Sci. 2018;19(3):831.doi:10.3390/ijms19030831
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijms19030831&link_type=DOI) 

24. 24.Neher JJ, Neniskyte U, Zhao JW, Bal-Price A, Tolkovsky AM, Brown GC. Inhibition of microglial phagocytosis is sufficient to prevent inflammatory neuronal death.JImmunol. 2011;18(8):4973–4983.doi:10.4049/jimmunol.1003600
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTA6IjE4Ni84LzQ5NzMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNy8wMS8yMDIzLjA2LjMwLjIzMjkyMDg0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

25. 25.Gallizioli M, Miro-Mur F, Otxoa-de-Amezaga A, et al. Dendritic Cells and Microglia Have Non-redundant Functions in the Inflamed Brain with Protective Effects of Type 1 cDCs. Cell Rep. 2020;33(3):108291.doi:10.1016/j.celrep.2020.108291
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.celrep.2020.108291&link_type=DOI) 

26. 26.Kwon HS, Koh SH. Neuroinflammation in neurodegenerative disorders: the roles of microglia and astrocytes. Transl Neurodegener. 2020;9(1). doi:10.1186/s40035-020-00221-2
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s40035-020-00221-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33239064&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F01%2F2023.06.30.23292084.atom) 

27. 27.Monterey MD, Wei H, Wu X, Wu JQ. The Many Faces of Astrocytes in Alzheimer’s Disease. Front Neurol. 2021;12:619626. doi:10.3389/fneur.2021.619626
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fneur.2021.619626&link_type=DOI) 

28. 28.Zhao W, Beers DR, Appel SH. Immune-mediated Mechanisms in the Pathoprogression of Amyotrophic Lateral Sclerosis.JNeuroimmune Pharmacol. 2013;8(4):888–899. doi:10.1007/s11481-013-9489-x
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11481-013-9489-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23881705&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F01%2F2023.06.30.23292084.atom) 

29. 29.Quek H, Cun-López C, Stewart R, et al.ALSmonocyte-derived microglia-like cells reveal cytoplasmic TDP-43accumulation,DNAdamage,and cell-specific impairment of phagocytosis associated with disease progression.JNeuroinflammation. 2022;19(1). doi:10.1186/s12974-022-02421-1
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12974-022-02421-1&link_type=DOI) 

30. 30.Humphrey J, Venkatesh S, Hasan R, et al. Integrative transcriptomic analysis of the amyotrophic lateral sclerosis spinal cord implicates glial activation and suggests new risk genes. Nat Neurosci. 2023;2(1):150–162.doi:10.1038/s41593-022-01205-3
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41593-022-01205-3&link_type=DOI) 

31. 31.Yogalingam G, Guo XH, Muller VJ, et al. Identification and molecular characterization of α-L-iduronidase mutations present in mucopolysaccharidosis type Ipatients undergoing enzyme replacement therapy. Hum Mutat. 2004;24(3):199–207. doi:10.1002/humu.20081
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/humu.20081&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15300847&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F01%2F2023.06.30.23292084.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000223304800003&link_type=ISI) 

32. 32.Ruffieux H, Carayol J, Popescu R, et al. Afully joint Bayesian quantitative trait locus mapping of human protein abundance in plasma. PLoS Comput Biol. 2020;1(6):e1007882. doi:10.1371/journal.pcbi.1007882
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1007882&link_type=DOI) 

33. 33.Yu E, Ambati A, Andersen MS, et al. Fine mapping of the HLAlocus in Parkinson’s disease in Europeans. npj Parkinson’s Disease. 2021;7(1). doi:10.1038/s41531-021-00231-5
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41531-021-00231-5&link_type=DOI) 

34. 34.Nalls MA, Blauwendraat C, Sargent L, et al. Evidence for GRNconnecting multiple neurodegenerative diseases. Brain Commun. 2021;3(2):fcab095. doi:10.1093/braincomms/fcab095
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/braincomms/fcab095&link_type=DOI) 

35. 35.Brunham LR, Kang MH, Van Karnebeek C, et al. Clinical,Biochemical,and Molecular Characterization of Novel Mutations in ABCA1in Families with Tangier Disease. In: Springer Berlin Heidelberg; 2014:51–62.doi:10.1007/8904\_2014\_348
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access\_num=10.1007/8904_2014_348&link_type=DOI) 

36. 36.Manzine PR, Ettcheto M, Cano A, et al.ADAM10in Alzheimer’s disease: Pharmacological modulation by natural compounds and its role as a peripheral marker. Biomed Pharmacother. 2019;113:108661.doi:10.1016/j.biopha.2019.108661
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopha.2019.108661&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30836275&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F07%2F01%2F2023.06.30.23292084.atom) 

37. 37.Shinjyo N, Kagaya W, Pekna M. Interaction Between the Complement System and Infectious Agents - APotential Mechanistic Link to Neurodegeneration and Dementia. Front Cell Neurosci. 2021;15:710390.doi:10.3389/fncel.2021.710390
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fncel.2021.710390&link_type=DOI) 

38. 38.Zhang L, Ashizawa T, Peng D. Primary coenzyme Q10deficiency due to COQ8A gene mutations. Molecular Genetics & Genomic Medicine. 2020;8(10):e1420. doi:10.1002/mgg3.1420
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mgg3.1420&link_type=DOI)