Neuropathological correlates and genetic architecture of microglial activation in elderly human brain

Microglia, the resident immune cells of the brain, have important roles in brain health. However, little is known about the regulation and consequences of microglial activation in the aging human brain. Here we report that the proportion of morphologically activated microglia (PAM) in postmortem cortical tissue is strongly associated with β-amyloid, tau-related neuropathology, and the rate of cognitive decline. Effect sizes for PAM measures are substantial, comparable to that of APOE ε4, the strongest genetic risk factor for Alzheimer’s disease, and mediation models support an upstream role for microglial activation in Alzheimer’s disease via accumulation of tau. Further, we identify a common variant (rs2997325) influencing PAM that also affects in vivo microglial activation measured by [11C]-PBR28 PET in an independent cohort. Thus, our analyses begin to uncover pathways regulating resident neuroinflammation and identify overlaps of PAM’s genetic architecture with those of Alzheimer’s disease and several other traits.

T he function of immune cells in the central nervous system (CNS) has recently become a major focus in human genetics, as these cells have been implicated in susceptibility to neurodegenerative, autoimmune, and psychiatric diseases. Microglia, the brain's resident immune cells, are thought to have important roles in both tempering and exacerbating agingrelated neuropathological processes, but their precise role remains unclear as they are difficult to access in human subjects. Recently, a molecularly defined subtype of disease-associated microglia has been proposed to exist in a mouse model of Alzheimer's disease (AD) 1 . However, transcriptomic identities of isolated microglia are notoriously plastic 2 and highly susceptible to a myriad of experimental confounds 3 . Regional and temporal heterogeneity of microglia subpopulations have also been shown in human and mouse models based on both molecular and morphological characteristics.
Recent postmortem studies have shown that microglial densities in specific regions are associated with a syndromic diagnosis of both early and late-onset AD 4 , and a recent systematic review of 113 studies quantifying microglial activation in postmortem AD brain highlighted the importance of activation vs. abundance of these cells in disease 5 . However, low sample sizes, indirect measures of microglia, and lack of full antemortem and postmortem pathological assessments all limit the insights that can be drawn from the individual component studies and this systematic review. Here, we leverage two large cohort studies of cognitive aging that include antemortem longitudinal cognitive assessments and structured postmortem histopathological evaluations to characterize a postmortem measure of microglial activation, directly observed by immunohistochemical staining and light microscopy. This morphological assessment of microglial activation stage represents a clear and robust measurement of neuroinflammation that cannot be captured by a surrogate marker. We first examine how this measure relates to different agingrelated pathologies. We follow this with causal mediation analyses aimed at placing microglial activation temporally within the cascade of pathological events leading to AD. Finally, we perform genome-wide analyses to identify the genomic architecture of microglial activation and implement a high-resolution polygenic risk scoring method based on Mendelian randomization assumptions to demonstrate putatively causal effects of microglial activation on multiple human diseases and traits. Figure 1 illustrates the set of analyses performed in our study.

Results
Active microglia discriminate pathological AD. The characteristics of ROS/MAP participants with microglial count data are presented in Table 1. We first performed pairwise Spearman correlations of each individual microglial density measurement followed by hierarchical clustering (Supplementary Figure 1), finding that stage I microglial densities were more similar between regions, whereas stage II and III microglial densities were more highly correlated within cortical and subcortical regions separately.
Following this observation and prior reports of the presence of morphologically defined active microglia in much smaller samples of AD brains (n range = 34-106) 5 , we performed benchmarking discriminatory analyses for each microglial density phenotype (see Fig. 2a for illustrative examples of staged microglia) in two steps. First, Welch t-tests comparing mean measures of microglial density between individuals with a postmortem pathological diagnosis of AD vs. non-AD found that total microglial density as measured in any region was not discriminative of AD status (0.33 < p < 0.68; Fig. 2b). On the other hand, stage III microglial density was different between AD and non-AD subjects, though this was only true in cortical regions (midfrontal (MF) p = 1.5 × 10 −8 , Cohen's d [95%CI] = 0.80 [0.52,1.08]; inferior temporal (IT) p = 6.4 × 10 −9 , Cohen's d [95% CI] = 0.84[0.55, 1.12]). Further, a stronger association was observed with the proportion of stage III microglia density relative to total microglia (the proportion of active microglia, or PAM) (MF p = 1.8 × 10 −10 , Cohen's d [95%CI] = 0.91[0.63,1.19]; IT p = 1.5 × 10 −11 , Cohen's d [95%CI] = 0.99[0.70,1.28]), confirming that morphologically activated microglia rather than the total number of microglia is most important for the accumulation of AD-related pathology in aging 6,7 . This is consistent with earlier work in mice demonstrating that microglial activation rather than proliferation mediates neurodegeneration 8 .
Second, logistic regression modeling of pathological AD in the same sample confirmed our t-test results, finding that models including PAM outperformed other models with a maximum AUC of 0.795 for IT and 0.792 for MF ( Fig. 2c; pathological distributions between pathological AD and non-AD subjects are described in Supplementary Data 9). Notably, the effects of both cortical PAM measures were independent of and improved model performance to a greater extent than the major APOE ε4 genetic risk factor for AD. In the case of IT, the inclusion of PAM increased model performance over the co-variate only model by 18% (AUC = 0.745 vs. 0.565), whereas including APOE ε4 status only yielded a 9.6% increase (AUC = 0.661  [2.46,7.15]). Importantly, neither APOE ε4 nor ε2 status were related to either cortical PAM measure, whether or not co-variates were included (all p > 0.1).
PAM effects on neuropathology and cognitive decline. Having found robust effects of PAM on pathologically defined AD, and given the high degree of neuropathological heterogeneity found in our cohort 9,10 , we sought to identify whether effects of PAM were specific to β-amyloid (Aβ) and tau (the defining pathologic characteristics of AD) or were also associated with other neuropathological features commonly observed in aged individuals. Using robust regression across 18 pathologies, we found that cortical PAM measures were associated to the greatest extent with total Aβ load and neuritic amyloid plaque count, and, to a lesser but still significant extent, with paired helical filament (PHF) tau, neurofibrillary tangles, and diffuse plaques (Fig. 3a). All relationships were in the positive direction, whereby an increase in PAM paralleled an increase in pathology. Notably, we found no significant associations of either subcortical PAM with any pathology, nor did we find association of any PAM measure with pathologies not related to Aβ or tau accumulation (Fig. 3a). Tests of person-specific linear trajectories of cognitive decline revealed significant associations of IT PAM with global cognitive decline, as well as with decline in all five cognitive sub-domains (Fig. 3b). MF PAM was also associated with global cognitive decline, and four of the five cognitive sub-domains. Similar to our neuropathological findings, there were no significant associations of subcortical PAM with cognitive decline.
To test whether associations were missed due to a lack of regional specificity in pathological measures, we performed posthoc analyses of all PAM measures against detailed regional pathology data (totaling 95 regional measures of amyloid, tau, Lewy body, and infarct pathology in both cortical and subcortical regions). These analyses mirrored the brain-wide results, finding exclusively amyloid-related and tau-related associations with cortical PAM measures; there were no significant associations for either subcortical PAM measure with any regional phenotypic measure (Supplementary Figure 3).
Mediation modeling of PAM. We next investigated whether activated microglia might contribute to or be the result of accumulating AD pathologies. Mediation analyses found no evidence for direct or mediator effects of PAM on cognitive decline in the presence of PHFtau; rather, these analyses pointed toward indirect effects of cortical PAM on cognitive decline via PHFtau (Supplementary Data 2). Both direct and indirect effects of PAM were found for PHFtau formation, and thus the sum of evidence across our models suggests a synergistic involvement of PAM and Aβ load in affecting global cognitive decline via their effects on PHFtau (Supplementary Figure 4). Put simply, our data suggest    Genetic architecture of microglial activation. Given the imperfect correlation between MF and IT PAM measures (Spearman ρ = 0.66), we performed two separate genome-wide association studies (GWAS). All significant and suggestive GWAS results are listed in Table 2 (details in Supplementary Datas 4 and 5). For IT PAM, a single locus on chromosome 1 reached genome-wide significance (rs183093970; linear regression β = 0.154, SE = 0.024, p = 5.47 × 10 −10 ) (Fig. 4a, b). However, while the 191 kb region encompassed by this lead SNP contains three independent signals, all three have low minor allele frequency (MAF) (0.02 > MAF > 0.015). Notably, this association was driven by only seven individuals in our sample carrying minor alleles tagging this haplotype and should therefore be considered cautiously. Beyond this genome-wide significant locus, 27 additional independent regions were associated at p < 1 × 10 −5 , and mapping based on position and combined eQTL evidence identified a total of 52 candidate genes as possible functional targets of these variants (Fig. 4c). For MF PAM, a different locus on chromosome 1 reached genome-wide significance (top SNP: rs2997325 T linear regression β = 0.039, SE = 0.0066, p = 1.88 × 10 −8 ; Fig. 4d, e, g). Beyond this lead SNP, 11 additional regions surpassed our suggestive threshold for association, which mapped to a total of 26 genes (Fig. 4f). In contrast to rs183093970, rs2997325 is a relatively common variant (MAF = 0.37), lies 8.9 kb 3′ of a long intergenic non-coding (Linc) RNA (RP11-170N11.1, LINC01361), and influences the expression of LINC01361 in multiple tissues with combined eQTL mapping evidence of p = 5.35 × 10 −12 12 . However, LINC01361 is not measured in our cortical RNAseq data. Q-value analysis revealed no significant overlap in genomic loci implicated by the two PAM GWAS, though many lead SNPs in the MF GWAS did reach nominal significance (p < 0.05) in the IT GWAS, and vice versa (Supplementary Data 4).
Gene enrichment analyses were performed separately on the two sets of mapped genes. For the 52 genes from the IT PAM GWAS, none of the 30 general tissue types analyzed in GTEx showed Bonferroni significant enrichment (Fig. 4j), though in more fine-grained analyses of 53 tissues, sigmoid colon was significantly enriched for differential expression of this gene set (Supplementary Figure 6b). Enrichment for functional gene categories and diseases found over 600 significant results, primarily relating to immunologic signatures (Supplementary Data 6). For the list of 26 mapped genes from the MF PAM GWAS, seven of the 30 general tissue types analyzed in GTEx showed Bonferroni significant enrichment (Fig. 4k). In contrast to      Both analyses co-varied for genotype platform, age at death, sex, postmortem interval, APOE ε4 status, and top three EIGENSTRAT principal components. g Regional association plot highlighting the MF PAM genome-wide significant locus surrounding rs2997325 (p = 1.88 × 10 −8 , n = 225). The color of each dot represents degree of linkage disequilibrium at that SNP based on 1000 Genomes Phase 3 reference data. The combined annotationdependent depletion (CADD) score is plotted below the regional association plot on a scale from 0 to 20, where 20 indicates a variant with highest predicted deleteriousness. The RegulomeDB score is plotted below the CADD score, and summarizes evidence for effects on regulatory elements of each plotted SNP (5 = transcription factor-binding site or DNAase peak, 6 = other motif altered, 7 = no evidence). Below the RegulomeDB plot is an eQTL plot showing -log 10 (p-values) for association of each plotted SNP with the expression of the mapped gene RP11-170N11.  Fig. 4i), consistent with our finding that the same allele increased MF PAM (Fig. 4h).
Role of activated microglia across human traits. Having generated genome-wide profiles of genetic risk for both cortical PAM measures, we used a high-resolution polygenic scoring-based method to test for overlap in the genomic underpinnings of microglial activation and AD, which would suggest a causal link between microglial activation and AD susceptibility. Secondarily, we tested 28 other brain and immune-related traits with publicly available GWAS data to more broadly assess the role of microglia in susceptibility to human disease. For MF PAM (Fig. 5a, b), AD showed the strongest polygenic association (p = 1.8 × 10 −10 , r 2 = 7.3 × 10 −4 ), and five other traits had optimal evidence for cogenetic regulation at a corrected threshold, with educational attainment showing the second strongest effect (p = 1.1 × 10 −5 , r 2 = 6.2 × 10 −5 ). For IT PAM, AD susceptibility (p = 4.9 × 10 −13 , r 2 = 9.4 × 10 −4 ) was also the most strongly associated, and 12 other traits met our significance threshold with educational attainment (p = 8.2 × 10 −7 , r 2 = 7.8 × 10 −5 ) also demonstrating the second strongest effect. These analyses provide evidence that the genetic predisposition to having activated microglia also contributes to making an individual more likely to develop AD. This result therefore goes beyond simply finding an enrichment of AD genes that happen to be expressed in microglia: they provide evidence that genomic propensity for the active microglia state is causally related to AD risk. The association with educational attainment is intriguing and consistent with our understanding that microglia play an important role in sculpting the developing brain by pruning synapses 13 . Likewise, the association of IT PAM with schizophrenia extends the narrative of the involvement of microglia in this neuropsychiatric disease 14 . To be thorough, we repeated these analyses in the reverse direction, asking whether genetic susceptibility for each of these 29 traits influences microglial activation. In these analyses, we found no association for AD, suggesting that the causal chain of events most likely flows from genetic risk → microglial activation → AD. We did find that primary sclerosing cholangitis (PSC) risk influenced both PAM measures while a few other traits demonstrated single associations (Fig. 5c, d). The PSC and inflammatory bowel disease associations may represent shared architecture relating to activated myeloid cells, since microglia and peripheral macrophages (which are implicated in these two inflammatory diseases) share many molecular functions.

Discussion
Microglial activation is a well-known phenomenon that has been implicated in a myriad of pathological processes. However, a majority of studies have been conducted either in small numbers of human subjects with limited clinicopathologic data or in murine model systems whose relevance to human disease is unclear given (1) that they do not recapitulate the events and conditions observed in the aging human brain and (2) the emerging understanding of significant differences between aging human and murine microglia 15,16 . Our study provides several important insights. First, microglial activation is not a general feature of the AD brain: we found it to be elevated in cortical but not in subcortical structures. Further work now needs to be conducted to sample a wider array of regions in individuals with a range of pathologic burden and those without pathology so that we can better understand the extent to which these regional differences may be related to the nature of the brain region itself or the extent to which that region is affected by amyloid or tau pathology in a given individual. Currently, the data are too sparse to be able to interpret our results beyond simply noting that PAM in subcortical regions do not appear to relate to the burden of AD pathology. Second, other common neuropathologies of older age do not appear to be significantly associated with activated microglia in the regions tested and therefore do not confound observed associations with Aβ and tau pathologies. Third, our modeling suggests that microglial activation leads to cognitive decline indirectly via the accumulation of PHFtau. This histologybased result is consistent with transcriptome-based findings we have recently reported 17 .
We also describe the previously unmapped genetic landscape of microglial activation measured in postmortem human brain; we found no effect of the APOE locus despite reports that APOE may be a ligand for important microglial receptors such as TREM2 18,19 . However, a discovery GWAS returned two significant results: the uncommon rs183093970 variant, whose role remains to be validated, and rs2997325, whose role we confirmed by assessing the related phenotype of in vivo microglial PET imaging. The effect of the rs2997325 T allele in increasing microglial activation and [ 11 C]-PBR28 binding should be validated more extensively, but, given its high frequency and strong effect, it may be a clinically relevant biomarker requiring attention in the many studies evaluating the utility of other microgliatargeted ligands in a clinical setting to diagnose or stage neurological diseases. While little is known about the biology of LINC01362 that appears to be influenced by rs2997325, several of the suggestive loci implicated by our PAM GWAS were found within or near genes of functional significance. For example, the Lymphocyte antigen 75 gene (LY75), codes for CD205, a dendritic cell surface receptor that interacts with MHC class I molecules 20 and plays an important role in T cell function 21 . Further investigation of cortical PAM genetic architecture is warranted to extend our initial set of observations; replication of associations with both histologic measures of microglial activation and TSPO imaging in populations at risk of AD are needed to assess the utility of these measures as biomarkers that may inform the diagnosis and/or management of aging-related cognitive decline. In addition, in vitro knock-down and DNA editing experiments, combined with amyloid and tau-sensitive assays, will be needed to explore the role of these variants in biological events related to AD.
There are important limitations to consider in our study. First, while we were able to test for some regional-specific effects of PAM on pathology, we did not always have data on pathology in the exact tissue samples in which microglia were counted. Therefore, tightly coupled pathology-microglial associations, such as those known to exist with acute infarction 22 , may have been missed. Whether or not microglial activation is a regionspecific phenomenon in aging is an unresolved question; both global and focal distributions have been reported in the aging brain, albeit using different measurements 23,24 . Reassuringly, a recent postmortem investigation of microglial activation measured morphologically in the brains of 11 late onset AD subjects vs. 12 age-matched controls also found increases in microglial activation in cortical but not sub-cortical regions 4 . Second, we have a limited sample size for genome-wide analyses, although it was sufficient to discover the rs2997325 variant which has a strong effect on both microglial activation and in vivo microglial imaging. Third, our moderate sample size also means that we cannot exclude the possibility that activated microglia may have weak, undetected effects on non-AD pathologies. Finally, we acknowledge that our chosen methodologies for quantifying microglia have intrinsic biases and limitations. For example, the manual identification of stained microglia carries a degree of subjectivity and, while [ 11 C]-PBR28 is a well-validated biomarker of microglial activation in humans, experimental confounds and the potential for non-specific binding are important to consider in any molecular imaging experiment. Further, we recognize that infiltrating macrophages present in brain tissue samples may represent a source of error in our counts since they are largely indistinguishable from microglia once they are activated. showing -log 10 (p-values) for genetic score (GRS) associations of inferior temporal cortex (IT) PAM and midfrontal cortex (MF) PAM across all p-value thresholds ranging from 0 to 0.5, with 5.0 × 10 −5 regular increments (10,000 total scores). Each PAM genetic score was mapped onto the summary statistics of each published GWAS trait and tested for significance. The dotted line represents a corrected statistical significance threshold of p = 8.6 × 10 −4 , corrected for 29 GWAS traits and two PAM measures. The width of each violin represents the density of PAM polygenic scores associated with each trait at a given significance. For example, a bottom-skewed violin (e.g. educational attainment (Edu)) indicates that a majority of scores across the tested set of thresholds tended to achieve greater significance, whereas a top-skewed violin (e.g. bipolar disorder (BPD)) indicates that a majority of tested scores tended toward lower significance. Peaks achieving at least one score above corrected statistical significance are labeled for their respective GWAS trait. b Network plot illustrating significant results for PAM → GWAS trait analyses, where an arrow indicates a directional effect of a peak PAM GRS on a GWAS trait at corrected significance (edges are labeled with % variance explained in GWAS trait by optimal PAM GRS, and edge thickness is proportional to the -log(p-value) of the association

ARTICLE
Reassuringly, recent single-cell analyses performed by our group have shown that these cells represent a negligible minority (~1%) of myeloid lineage cells present in brain tissue from similar subjects 25 . These same single-cell sequencing experiments will be important for identifying molecular subpopulations of microglia that may correspond to those activation states captured by the PAM phenotype. Despite these limitations, our polygenic analyses revealed significant links that clarify and expand the roles of activated microglia in human diseases. This is clearest for AD where we refine the narrative of the involvement of myeloid cells in AD susceptibility by showing that the proportion of activated microglia has a causal role in AD. A similar result for cognitive attainment is intriguing as it may be capturing microglia's role in synaptic pruning during development and memory consolidation; however, this is confounded by the fact that this enrichment may also be capturing the role of microglia in age-related cognitive decline that we describe in this report. The enrichment for Crohn's disease is unlikely to represent a causal role of microglia in this disease; rather, this result captures the likely overlap of genetic architecture between activated microglia and other activated myeloid cells in the periphery that are known to be involved in Crohn's disease.
The relevance of recent reports of pathologically important subtypes of microglia in murine model systems, such as diseaseassociated microglia (DAM) 1 and dark microglia 26 , remains to be demonstrated in the aging human brain since mouse and human microglia diverge significantly in functional molecular signatures with age 27 . By studying activated microglia in the target organ in humans, we have found several key insights, including (1) the observation that they may accelerate cognitive decline through an effect on PHFtau accumulation, which enables us to design mechanistic studies, (2) the discovery and validation of a chromosome 1 locus (among several other suggestive plausible loci) that provides a biological foundation for dissecting the mechanisms of microglial activation, and (3) the effect of rs2997325 on [ 11 C]-PBR28 binding which may need to be accounted for in human imaging studies of in vivo microglial activation. Thus, microglial activation is a central component of AD susceptibility, and we have begun to elaborate its place in the causal chain of events leading to increased accumulation of tau pathology and subsequent cognitive decline, as well as regulatory mechanisms that influence this activation.

Methods
Study subjects. All antemortem cognitive and postmortem data analyzed in this study were gathered as part of the Religious Orders Study and Memory and Aging Project (ROS/MAP) [28][29][30] , two longitudinal cohort studies of the elderly, one from across the United States and the other from the greater Chicago area. All subjects were recruited free of dementia (mean age at entry = 78 ± 8.7 (SD) years), agreed to annual clinical and neurocognitive evaluation, and signed an Anatomical Gift Act allowing for brain autopsy at time of death. In vivo [ 11 C]-PBR28 PET imaging acquisitions were collected on data from the Indiana Memory and Aging Study (IMAS), an ongoing neuroimaging and biomarker study based at the Indiana University School of Medicine including elderly subjects at multiple levels of cognitive impairment 31 . Written informed consent was obtained from all ROS/ MAP and IMAS participants. Study protocols were approved by each site's Institutional Review Board. Full methods can be found in Supplementary Methods.
Genomics. Genotype array data for 2067 ROS/MAP subjects were imputed using the Michigan Imputation Server (Haplotype Reference Consortium reference v1.1). For IMAS, data were imputed using IMPUTE v2.2 (1000 Genomes phase 1 reference). APOE (rs429358 and rs7412) genotyping was carried out separately using standard protocols.
Tanscriptomics. RNA sequencing of postmortem dorsolateral prefrontal cortical tissue from 538 ROS/MAP subjects were available for analysis at the time of study. Following sequencing and standard quality control 32 , The Speakeasy algorithm 33 was used to cluster expressed genes into functionally cohesive gene modules, which have been extensively validated for robustness and pathophysiological relevance 16,34 . Average values of gene expression for each of 47 modules were used as quantitative outcomes.
Proteomics. Selected reaction monitoring-based (SRM) quantitative proteomics was used to measure the abundance of 67 proteins (Supplementary Data 3) in frozen dorsolateral prefrontal cortical tissue from 807 ROS/MAP participants according to a standard protocol 35,36 .
Postmortem neuropathology. A total of 18 disease-related and age-related neuropathologies were measured brain-wide (n = 985). Pathologies included multiple validated measures of Aβ peptides, neuritic and diffuse plaques, hyperphosphorylated tau protein, neurofibrillary tangles, micro and macro cerebral infarcts, cerebral atherosclerosis, TDP43 proteinopathy, and hippocampal sclerosis, among others (see Supplementary Methods) 37 . A subset of up to 225 brain samples were also evaluated for the presence of microglia at three stages of activation in four regions (midfrontal (MF) cortex, inferior temporal (IT) cortex, ventral medial caudate (VM), and posterior putamen (PPUT)), based on morphology: stage I (thin ramified processes), stage II (plump cytoplasm and thicker processes), and stage III (appearance of macrophages).
Antemortem cognitive decline. A total of 1932 subjects with genomic data also had longitudinal cognitive performance data available at the time of study (the full list of cognitive tests can be found in Supplementary Data 10).
In vivo TSPO PET imaging. Scans of in vivo [ 11 C]-PBR28 binding were assessed in 27 subjects (n CN = 13, n MCI = 7, n AD = 7) as part of the IMAS study 31 .
Statistical analysis. Regression analyses were performed in R (v3.3.3) 38 . PAM was calculated by the following formula: where r represents each of four regions and S1, S2, and S3 represent microglial densities measured in region r at stage I, II, and III, respectively. To address potential concerns related to PAM distributions and model fitting, we performed extensive model validation and sensitivity analyses (Supplementary Methods; Supplementary Figure 2; Supplementary Data 1). Differences in model performance are reported as change in C-index, expressed as a percent, which corresponds to the area under the receiver operating characteristics curve (ΔAUC). Iterative re-weighted least-squares robust regression was used and model validation was performed using the .632+ bootstrap method 39 . Models were corrected using the Bonferroni procedure. Causal mediation modeling for identifying direct and indirect effects of PAM on cognitive decline was performed using the 'mediation' R package.
GWAS were performed in PLINK 40 (v1.90) using imputed genotype dosages and additive linear models, co-varying for age at death, postmortem interval (PMI), sex, genotype batch, and the first three EIGENSTRAT 41 principal components. Significance thresholds of p < 2.5 × 10 −8 and p < 1.0 × 10 −5 (two-sided) were deemed genome-wide significant and suggestive, respectfully. Post-processing of GWAS results was conducted using the full complement of state-of-the-art tools available through the recently released Functional Mapping and Annotation of Genome-Wide Association Studies platform (FUMAGWAS; http://fuma.ctglab.nl/) 42 . Analyses included positional and eQTL-based gene mapping, assessment of Combined Annotation Dependent Depletion (CADD v1.3) 43 and RegulomeDB scores 44 , tissue expression specificity using GTEx v7, comprehensive gene set enrichment analyses, and mapping of gene targets to the DrugBank database 45 .
Statistical Parametric Mapping version 8 (SPM8) was used for imaging analysis. Freesurfer (v5.1) was used to define subject-specific regions of interest and average standardized uptake value ratios (SUVRs) were extracted for three bilateral ROIs. As six ROI values were strongly correlated and the sample size is moderate for genetic analysis of this type (n = 27), we carried out a multivariate analysis by performing genetic association analysis using six phenotypes simultaneously to minimize the number of test performed and increase the statistical power. We used the software for correlated phenotype analysis (SCOPA) program 46 to perform a single omnibus genetic test for all correlated phenotypes, modeling rs2997325 genotype additively.
To assess causal relationships between the genetic determinants of cortical PAM (based on our PAM GWAS results) and 29 brain and immune-related traits (and vice versa), we used a genetic risk score-based method 47  Reporting summary. Further information on experimental design is available in the Nature Research Reporting Summary linked to this article.