# Meta-analysis of genome-wide association studies identifies ancestry-specific associations underlying circulating total tau levels

Chloe Sarnowski, Mohsen Ghanbari, Joshua C. Bis, Mark Logue, Myriam Fornage, Aniket Mishra, Shahzad Ahmad, Alexa S. Beiser, Eric Boerwinkle, Vincent Bouteloup, et al.

► To cite this version:

Chloe Sarnowski, Mohsen Ghanbari, Joshua C. Bis, Mark Logue, Myriam Fornage, et al. Meta-analysis of genome-wide association studies identifies ancestry-specific associations underlying circulating total tau levels. Communications Biology, 2022, 5 (1), 10.1038/s42003-022-03287-y. hal-03664665

HAL Id: hal-03664665

https://hal.science/hal-03664665

Submitted on 11 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

# Meta-analysis of genome-wide association studies identifies ancestry-specific associations underlying circulating total tau levels

Chloé Sarnowski ID 1,30✉, Mohsen Ghanbari ID 2,3,30, Joshua C. Bis ID 4,30, Mark Logue 5,6,30, Myriam Fornage ID 7,30, Aniket Mishra ID 8,30, Shahzad Ahmad ID 2,9, Alexa S. Beiser 10,11,12, Eric Boerwinkle 7, Vincent Bouteloup 13, Vincent Chouraki 14, L Adrienne Cupples ID 10,11, Vincent Damotte 14, Charles S. DeCarli 15, Anita L. DeStefano 10, Luc Djoussé 16, Alison E. Fohner 17, Carol E. Franz ID 18, Tiffany F. Kautz 19, Jean-Charles Lambert ID 14, Michael J. Lyons 20, Thomas H. Mosley 21, Kenneth J. Mukamal ID 22, Matthew P. Pase ID 23,24, Eliana C. Portilla Fernandez 2, Robert A. Rissman 25, Claudia L. Satizabal 11,12,19, Ramachandran S. Vasan ID 11,26, Amber Yaqub ID 2, Stephanie Debette ID 8,27,31, Carole Dufouil ID 27,31, Lenore J. Launer ID 28,31, William S. Kremen ID 18,31, William T. Longstreth 29,31, M Arfan Ikram ID 2,31 & Sudha Seshadri 11,12,19,31

Circulating total-tau levels can be used as an endophenotype to identify genetic risk factors for tauopathies and related neurological disorders. Here, we confirmed and better characterized the association of the 17q21 MAPT locus with circulating total-tau in 14,721 European participants and identified three novel loci in 953 African American participants (4q31, 5p13, and 6q25) at $P < 5 \times 10^{-8}$. We additionally detected 14 novel loci at $P < 5 \times 10^{-7}$, specific to either Europeans or African Americans. Using whole-exome sequence data in 2,279 European participants, we identified ten genes associated with circulating total-tau when aggregating rare variants. Our genetic study sheds light on genes reported to be associated with neurological diseases including stroke, Alzheimer's, and Parkinson's (F5, MAP1B, and BCAS3), with Alzheimer's pathological hallmarks (ADAMTS12, IL15, and FHIT), or with an important function in the brain (PARD3, ELFN2, UBASH3B, SLIT3, and NSD3), and suggests that the genetic architecture of circulating total-tau may differ according to ancestry.
The protein tau is an important biomarker of neuronal injury and neurodegeneration. Alzheimer's disease (AD) and other dementias or related neurological disorders are associated with abnormal intraneuronal tau aggregates (collectively known as tauopathies)¹. Newer techniques to diagnose AD now examine CSF biomarkers to improve diagnostic certainty and aid in earlier diagnosis²⁻⁴. However, their collection is invasive and user variability can be large in the downstream quantification assays. In addition, CSF tau levels are normal or low in tauopathies like Progressive Supranuclear Palsy (PSP) and in frontotemporal dementia patients with tau mutations⁵,⁶.

Using blood biomarkers with high specificity and sensitivity for AD is ideal to lower cost, risk, and burden³. Circulating total-tau (t-tau) levels can be quantified in serum or in plasma⁷ early in AD due to blood brain barrier breakdown⁸,⁹. Particularly, they show promise as a predictive biomarker for dementia and related endophenotypes¹⁰, with higher levels in patients with dementia or mild cognitive impairment compared to controls¹¹,¹², and higher levels associated with poorer cognitive performance, and smaller hippocampal volumes¹³,¹⁴. However, elevated levels may lack diagnostic specificity for AD, and simply indicate that brain injury is common to several neurological diseases. A recent paper showed for example that higher circulating t-tau predicted a higher risk of incident stroke¹⁵. Importantly, circulating biomarkers do not need to mirror their level in CSF to be useful. Altogether, the recent literature suggests that circulating t-tau levels may be a predictive biomarker to improve risk stratification for dementia and assess AD's progression, help with enrollment of high-risk individuals into dementia prevention trials, be useful in addition to other blood biomarkers of neurodegeneration to determine cognitive improvements in clinical trials, and represent a useful biomarker for AD when added to CSF tau measures¹⁰,¹⁵⁻¹⁷.

CSF t-tau and phosphorylated tau (p-tau) levels have been used as endophenotypes in genome-wide association studies (GWAS) to detect genetic variants associated with AD risk. Similarly, circulating t-tau levels may be used as an endophenotype to identify genetic risk factors for tauopathies and related neurological disorders. Only two GWAS, conducted in the Alzheimer's Disease Neuroimaging Initiative (ADNI) study, were published for plasma t-tau or p-tau levels¹⁸,¹⁹. The modest sample size and the inclusion of only European participants has limited the statistical power to identify potential novel associations (only MAPT and APOE loci associations were statistically significant for t-tau and p-tau respectively), the ability to explore less frequent genetic variation, as well as the generalization of the findings to other ancestries. Therefore, the aim of our study was to perform large-scale meta-analyses of circulating t-tau levels, using 15,674 participants from eight studies representing two ancestries (Europeans and African Americans), to explore genetic variation underlying circulating t-tau levels and assess their overlap with known genetic determinants of neurological diseases. We detected four ancestry-specific loci at the genome-wide significance (17q21 in Europeans, and 4q31, 5p13, and 6q25 in African Americans). We identified pleiotropic associations at 17q21 and 1q24 which, combined with the detection of an enrichment of genes associated with neurological diseases or related traits, suggested that a potential overlap exists between genetic determinants of circulating t-tau levels and several neurological disorders and traits including AD and stroke.

## Results

**Populations and participants.** We included in our meta-analyses 15,674 participants from eight studies representing two major ancestries: Europeans (N = 14,721) and African Americans

(N = 953) (Table 1, 2). A description of each study is included in the Supplementary Notes 1–8. Proportion of males varied from 35% (The Cardiovascular Health Study, African Americans) to 100% (The Vietnam Era Twin Study of Aging). Mean age ranged from 49 years (3.5) for The Coronary Artery Risk Development in Young Adults, African Americans to 78 years (4.3) for The Cardiovascular Health Study, African Americans (Table 2).

**Main meta-analysis results.** Quantile-Quantile and Manhattan plots for each ancestry-specific meta-analysis of circulating t-tau levels are presented in Figs. 1, 2 and Supplementary Fig. 1. We confirmed and better characterized the strong association of genetic variants at the microtubule associated protein tau (MAPT) 17q21 locus in European participants (lead genetic variant rs242557-A, effect allele frequency = 0.38, beta = 0.20, P = 8.9 × 10⁻¹⁴³, proportion of variance explained (PVE) ~4.2%, Table 3). Stepwise model selection procedure at the 17q21 locus identified three distinct signals (rs242557-A, beta = 0.20, P = 2.3 × 10⁻¹⁴³; rs7502280-T, beta = 0.17, P = 5.7 × 10⁻³⁸, PVE ~ 1.1%, and rs2942003-T, beta = 0.16, P = 3.9 × 10⁻⁷⁸, PVE ~ 2.3%), Table 3, Supplementary Table 1 and Fig. 3. We also detected, at the genome-wide threshold (P < 5 × 10⁻⁸), three potential novel loci in African American participants at 4q31 (lead genetic variant rs111836296-T, effect allele frequency = 0.06, beta = -0.54, P = 1.7 × 10⁻⁸), 5p13 (rs74710969-T, effect allele frequency = 0.11, beta = -0.53, P = 3.4 × 10⁻⁸, PVE ~ 3%) and 6q25 (rs674432-C, effect allele frequency = 0.97, beta = 0.68, P = 1.8 × 10⁻⁸, Table 3). The variants rs111836296 and rs74710969 are extremely rare or monomorphic in European populations. We detected additional associations at P < 5 × 10⁻⁷ at eleven loci (2q23, 3p14, 3q11, 5q13, 7p21, 7p15, 7q36, 8p11, 10p11, 10q23, and 14q32) in African American participants and three loci (1q24, 3p14, and 17q23) in European participants (Supplementary Tables 2, 3). All the identified hits were ancestry specific (absence of association or monomorphic variant in one ancestry) as indicated by the high heterogeneity (Table 3 and Supplementary Table 3). Forest plots are presented in Supplementary Figs. 2–11 and indicated consistent direction of effects across studies. Regional association plots for all loci are presented in Supplementary Figs. 12–27.

We did not observe an association of the two SNPs defining APOE in the main circulating t-tau European meta-analysis, consistent with previous finding.¹⁸ In the additional APOE4-stratified analyses, we observed similar magnitude and consistent direction of effects for the main associations identified in Europeans (Table 4).

**Overlap of circulating t-tau genetic determinants with neurological diseases and traits.** In addition to the strong associations at the MAPT locus that is known to be pleiotropic, we identified association at 1q24, a locus previously reported for stroke, Supplementary Table 4.²⁰ The lead genetic variant in our analysis (rs6686805) was in linkage disequilibrium with rs1800594, a GWAS hit for ischemic stroke. Analyses conducted with FUMA, based on the main results from the European circulating t-tau meta-analysis, identified significantly differentially expressed genes in brain cerebellar hemisphere and brain cerebellum (Supplementary Fig. 28) and enrichment of genes in gene sets reported by GWAS of neurological diseases or traits including Parkinson Disease (PD), craniofacial microsomia, intracranial volume, cognitive function, subcortical brain region volumes, and AD in APOE E4- carriers, as well as risk factors such as body mass index (Supplementary Fig. 29). Finally, a genetic risk score (GRS) based on the distinct genome-wide associations (two MAPT genetic variants, rs242557 and rs376284405) from our European meta-analysis of circulating t-tau levels (excluding the
<table><caption>Table 1 Description of the European-ancestry participants included in the meta-analysis of circulating total-tau levels.</caption><thead><tr><th></th><th>FHS (N = 6018)</th><th>RSI (N = 2169)</th><th>RSII (N = 960)</th><th>MEM1 (N = 336)</th><th>MEM2 (N = 7738)</th><th>CARDIA (N = 315)</th><th>CHS (N = 1396)</th><th>ARIC (N = 549)</th><th>VETSA (N = 754)</th><th>ADNI (N = 486)</th></tr></thead><tbody><tr><td>Men, N (%)</td><td>2790 (46)</td><td>1696 (59)</td><td>951 (44)</td><td>153 (45)</td><td>666 (37)</td><td>149 (47)</td><td>523 (37)</td><td>257 (47)</td><td>754 (100)</td><td>292 (60)</td></tr><tr><td>Age, mean (SD)</td><td>56.81 (13.98)</td><td>75.30 (6.10)</td><td>67.80 (7.10)</td><td>67.60 (8.58)</td><td>71.47 (8.56)</td><td>51.02 (3.25)</td><td>63.52 (4.38)</td><td>64.0 (60.0-67.0)</td><td>67.50 (2.51)</td><td>75.44 (6.71)</td></tr><tr><td>Age, median [25-75%]</td><td>56.1 [46.4-66.3]</td><td>74.4 [70.5-79.4]</td><td>65.5 [62.8-70.8]</td><td>68.0 (62.9-73.6)</td><td>72.2 [66.1-77.6]</td><td>52.0 [49.0-53.5]</td><td>77.91 (4.3)</td><td>77.0 (75.0-80.3)</td><td>68.1 [65.4-69.7]</td><td>75.8 [71.5-80.1]</td></tr><tr><td>Circulating t-tau, mean (SD)</td><td>4.13 (2.12)</td><td>2.60 (0.90)</td><td>2.60 (3.10)</td><td>2.10 (0.88)</td><td>2.23 (1.84)</td><td>0.54 (1.26)</td><td>0.46 (2.18)</td><td>0.41 (0.54)</td><td>2.02 (1.23)</td><td>2.88 (1.64)</td></tr><tr><td>Circulating t-tau, median [25-75%]</td><td>3.91 (3.24-4.73)</td><td>2.50 (1.90-3.0)</td><td>2.40 (1.90-2.90)</td><td>2.01 [1.55-2.60]</td><td>1.86 [1.41-2.60]</td><td>0.34 (0.20-0.57)</td><td>0.27 (0.17-0.43)</td><td>0.18 (0.09-0.34)</td><td>1.76 [1.29-2.32]</td><td>2.71 [1.93-3.45]</td></tr><tr><td>Circulating t-tau (log), mean (SD)</td><td>1.97 (0.46)</td><td>1.30 (0.50)</td><td>1.20 (0.50)</td><td>0.91 (0.81)</td><td>0.95 (0.73)</td><td>-1.54 (1.14)</td><td>-1.31 (0.83)</td><td>-1.68 (0.96)</td><td>0.82 (0.72)</td><td>0.92 (0.57)</td></tr><tr><td>Circulating t-tau (log), median [25-75%]</td><td>1.97 (1.70-2.24)</td><td>1.30 (0.90-1.60)</td><td>1.20 (0.90-1.60)</td><td>1.01 [0.63-1.38]</td><td>0.90 [0.50-1.38]</td><td>-1.56 [-2.32, -0.81]</td><td>-1.30 [-1.75, -0.85]</td><td>-1.71 [-2.34, -1.10]</td><td>0.81 [0.37-1.21]</td><td>1.00 [0.66-1.24]</td></tr></tbody></table>

FHS Framingham Heart Study, RSII and RSII The Rotterdam Study, MEMI and MEM2 The MEMENTO Study, CARDIA The Coronary Artery Risk Development in Young Adults Study, CHS The Cardiovascular Health Study, VETSA The Vietnam Era Twin Study of Aging Study, ARIC The Atherosclerosis Risk in Communities Study, ADNI The Alzheimer's Disease Neuroimaging Initiative Study.

Framingham Heart Study, FHS) was strongly associated with circulating t-tau levels (beta = 0.3, *P* = 4 × 10⁻⁹⁷, PVE = 7%) and was associated with intracranial volume (beta = 15.1, *P* = 3 × 10⁻⁴) in FHS. We did not detect significant associations with the other traits tested (Supplementary Table 5). Altogether, these findings suggest an overlap of the genetic associations of circulating t-tau levels with known genetic determinants of neurological disorders and associated traits.

**Two sample Mendelian Randomization (MR) analyses.** Using two sample MR and large GWAS summary statistics, we did not identify significant causal associations between circulating t-tau levels and AD, PD, stroke, or White Matter Hyperintensities (WMH) (Supplementary Table 6). We also tested the opposite hypothesis and did not find significant causal associations (Supplementary Table 7). Our results did not indicate significant heterogeneity or presence of directional horizontal pleiotropy, except for a few analyses (WMH or PD as exposure; stroke as outcome). We also performed power calculations for the MR where circulating t-tau levels was the exposure variable (Supplementary Table 8) that indicated that our analysis would be underpowered if the instruments had small effects on the neurological outcomes (especially for AD and PD analyses with smaller numbers of cases).

**Rare variant analyses using whole-exome sequence data.** Using SKAT (variance component test) or CMC (burden test), our rare variant aggregation tests based on whole-exome sequence data identified 10 genes (*ELFN2*, *UBASH3B*, *RUSF1*, *ZFP28*, *LCT*, *REM1*, *DELE1*, *SLIT3*, *NSD3*, and *MYO1G*) significantly (*P* ≤ 1.25 × 10⁻⁶) associated with circulating t-tau levels when aggregating rare variants (MAF ≤ 5% or MAF ≤ 1%) with high or moderate impact, including some missense and loss of function variants (Supplementary Tables 9, 10 and Supplementary Figs. 30, 31). All except one gene (*MYO1G*) were detected with SKAT at the gene level significance threshold (*P* = 1.25 × 10⁻⁶), while at least nominally significant associations were observed for most genes with CMC (Supplementary Tables 9, 10). These results indicate that rare variants in those genes were likely to have different magnitudes and directions of effects, including no effect. This is a likely scenario as the number of rare variants aggregated for each gene was somewhat large, ranging from 13 to 64 variants. Similar results were observed for the two sets of annotations tested (missense or loss of function versus high or moderate impact). This observation, combined with the fact that one set of annotations is a subset of the other, and the number of genetic variants contributing to both analyses did not differ drastically, suggested that the same variants were selected to be aggregated in both analyses. Similar results were also observed for the two MAF thresholds tested, except for *NSD3* and *SLIT3*. For these two genes, the addition of a small number of more frequent variants (1% < MAF ≤ 5%) attenuated the association.

Discussion

The goal of our study was to characterize genetic variation underlying circulating t-tau levels and to explore their overlap with known genetic determinants of neurological diseases. By performing large-scale meta-analyses in more than 15,000 participants from two major ancestries, we identified new genetic variants and genes associated with circulating t-tau levels, all associations being observed in only one ancestry (African Americans or Europeans). We identified pleiotropic signals at two regions (17q21 and 1q24) that were previously reported for plasma t-tau, AD, PD, WMH, and PSP (*MAPT*) or stroke (*F5*), respectively, and enrichment of genes associated with
<table><caption>Table 2 Description of the African American participants included in the meta-analysis of circulating total-tau levels.</caption><thead><tr><th></th><th>CARDIA (N = 111)</th><th>CHS (N = 273)</th><th>ARIC (N = 569)</th></tr></thead><tbody><tr><td>Men, N (%)</td><td>51 (46)</td><td>96 (35)</td><td>218 (38)</td></tr><tr><td>Age, mean (SD)</td><td>48.78 (3.46)</td><td>76.32 (4.93)</td><td>61.68 (4.47)</td></tr><tr><td>Age, median [25-75%]</td><td>48.0 [46.0-51.0]</td><td>75.0 [72.0-80.0]</td><td>61.0 [58.0-65.0]</td></tr><tr><td>Circulating t-tau, mean (SD)</td><td>0.57 (1.41)</td><td>0.46 (0.94)</td><td>0.50 (0.77)</td></tr><tr><td>Circulating t-tau, median [25-75%]</td><td>0.37 [0.22-0.60]</td><td>0.28 [0.18-0.46]</td><td>0.31 [0.20-0.50]</td></tr><tr><td>Circulating t-tau (log), mean (SD)</td><td>-1.54 (1.03)</td><td>-1.24 (0.85)</td><td>-1.57 (1.10)</td></tr><tr><td>Circulating t-tau (log), median [25-75%]</td><td>-1.43 [-2.24, -0.78]</td><td>-1.27 [-1.73, -0.78]</td><td>-1.69 [-2.36, -0.99]</td></tr></tbody></table>

CARDIA The Coronary Artery Risk Development in Young Adults Study, CHS The Cardiovascular Health Study, ARIC The Atherosclerosis Risk in Communities Study.

African-American meta-analysis of circulating total-tau levels

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Manhattan plot of association *P* values for the African American specific meta-analysis of GWAS of circulating total-tau levels. The -log10(*P*)-value for each single nucleotide variant on the y axis is plotted against the build 37 genomic position on the x axis (chromosomal coordinate). The dashed horizontal red line indicates the genome-wide significance threshold of *P* = 5 × 10⁻⁸ and the dashed horizontal black line indicates the threshold of *P* = 5 × 10⁻⁷.</figcaption></figure>

neurological diseases or related traits. Thus, our analyses highlighted that an overlap may exist between genetic determinants of circulating t-tau levels and several neurological disorders and traits including AD.

Our findings confirmed the importance of genes and pathways already well known to be involved in AD or other tauopathies and neurological diseases. Indeed, we first confirmed the strong association in Europeans of the 17q21 *MAPT* locus (lead genetic variant rs242557), which has been reported to be associated with circulating t-tau levels¹⁸. The *MAPT* locus has also been associated with AD, PD, and PSP, Supplementary Table 4²⁰⁻²², indicating an important role of *MAPT* in many neurodegenerative diseases. This locus has also been associated with head size²³⁻²⁵ and notably child head circumference²³, which may indicate possible effects of this inversion on brain development very early in life. We found an enrichment for gene sets reported associated with craniofacial microsomia, intracranial volume, and subcortical brain region volumes, which tend to also support this hypothesis (Supplementary Fig. 29). We also identified a

significant positive association of a GRS, constructed based on two distinct *MAPT* genetic variants (rs242557 and rs376284405), with intracranial volume in the FHS, while these variants were distinct from the ones reported by GWAS of intracranial volume at this locus (rs199525, rs8072451, and rs17689882). Despite PD is not a tauopathy, PSP and corticobasal degeneration, two PD subtypes known as Parkinson-plus syndromes, are both associated with the formation of tau deposits²⁶. Here we were also able to identify two additional and distinct signals at 17q21 (rs7502280 and rs2942003). The variant rs7502280 is located at 29 kb of the corticotropin releasing hormone receptor 1 (*CRHR1*) and at 7.3 kb of the mitogen-activated protein kinase 8 interacting protein 1 (*LOC644172*), and is a GWAS hit for relative carbohydrate intake²⁷ and sleep duration²⁸. The variant rs2942003 lies at 12 kb and tags the leucine rich repeat containing 37 member A2 (*LRRC37A2*) gene (Fig. 3). One 17q21 variant (rs439945) reported by the Parkinson Disease GWAS Consortium was found to be significantly associated with nearby gene expression probes targeting *LRRC37A* and *LRRC37A2* by a study investigating the
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Manhattan plot of association *P* values for the European ancestry specific meta-analysis of GWAS of circulating total-tau levels. The -log10(*P*)-value for each single nucleotide variant on the y axis is plotted against the build 37 genomic position on the x axis (chromosomal coordinate). The dashed horizontal red line indicates the genome-wide significance threshold of *P* = 5 × 10⁻⁸ and the dashed horizontal black line indicates the threshold of *P* = 5 × 10⁻⁷. The y axis was truncated for ease of interpretation.</figcaption></figure>

modification of gene expression in prefrontal cortex brain samples of pathologically confirmed PD cases and controls²⁹. Thus, genetic variations in both *MAPT* and *LRRC37A2* appear to be important determinants of tauopathies and neurodegenerative disorders. More research needs to be performed to understand more precisely the mechanisms underlying their contributions to other tauopathies. Particularly, the 17q21 region, a common inversion polymorphism, is complex and may affect the expression of other genes in the region that may also be involved in neurodegenerative disease pathology, possibly in a tissue-specific manner.

In addition to the *MAPT* 17q21 locus identified in Europeans, we detected three potential novel loci in African American participants (4q31, 5p13, and 6q25) at the genome-wide significance level. Two of the lead genetic variants (rs111836296 and rs74710969) were extremely rare in Europeans and lie in or tag candidate genes (*IL15* and *ADAMTS12*) linked to AD and other neurological disorders. The genetic variant rs111836296 at 4q31 lies at 6 kb and tags the interleukin 15 (*IL15*) gene (Supplementary Fig. 15). Serum IL15, a pro-inflammatory cytokine, has been studied as a possible marker of AD³⁰. The genetic variant rs74710969 at 5p13 lies in an intron of the ADAM metallopeptidase with thrombospondin type 1 motif 12 (*ADAMTS12*) gene (Supplementary Fig. 16). Previous studies have associated ADAMTSs family of secreted metalloproteases with the repair of the central nervous system, through its ability to degrade neurocan, a novel component of brain extra-cellular matrix. Alterations in this degradation processes could be associated with the pathogenesis of neurological disorders³¹. Several studies also suggest a role for *ADAMTS12* in stroke³²,³³.

We also detected 14 loci at *P* < 5 × 10⁻⁷, eleven loci in African American participants and three loci in Europeans.

Three signals identified in African Americans (3p14 and 5q13) or Europeans (17q23) lie in genes related to the tubulin-

microtubule system (*FHIT*, *MAP1B*, and *BCAS3*). The signal at 3p14 lies in the fragile histidine triad diadenosine triphosphatase (*FHIT*) gene (Supplementary Fig. 13) and the encoded protein interacts with tubulin. The signal at 5q13 lies in the microtubule associated protein 1B (*MAP1B*) gene (Supplementary Fig. 17). Proteins of this family may be involved in microtubule assembly, which is an essential step in neurogenesis. Gene knockout studies of the mouse *MAP1B* gene suggested an important role in development and function of the nervous system. Several studies are also in favor of a role of *MAP1B* in AD³⁴,³⁵. MAP1B is also a component of cortical Lewy bodies and binds alpha-synuclein filaments, which suggests that it may be involved in the pathogenesis of Lewy bodies³⁶. The signal at 17q23 lies in the BCAS3 microtubule associated cell migration factor (*BCAS3*) gene (Supplementary Fig. 27) that is highly expressed in the brain (GTEx). In mice, Rudhira, a murine WD40 domain protein that is 98% identical to *BCAS3*, has been shown to bind to microtubules and vimentin intermediate filaments to promote cell migration for angiogenic remodeling³⁷.

Furthermore, two signals (1q24 and 10p11) lie in genes reported to be associated with AD or stroke or have an important function in the brain (*F5*, and *PARD3*). The signal at 1q24 identified in Europeans lies in the coagulation factor 5 (*F5*) gene (Supplementary Fig. 24). The lead genetic variant rs6686805 is in linkage disequilibrium with rs1800594, a GWAS hit for blood protein levels and ischemic stroke (Supplementary Tables 4, 11),³⁸⁻⁴⁰ and with rs6030, a missense variant in *F5*. A rare protective variant in *F5* (rs2027885) has been reported to be associated with AD in African Americans⁴¹ and with hippocampal atrophy⁴². The signal at 10p11 identified in African Americans lies in the par-3 family cell polarity regulator (*PARD3*) gene (Supplementary Fig. 21) that is required for establishment of neuronal polarity and normal axon formation in cultured hippocampal neurons⁴³,⁴⁴. Par3 regulates microtubule stability and
<table><caption>Table 3 Lead genetic variants passing the genome-wide significance threshold ($P < 5 \times 10^{-8}$) in the meta-analysis of GWAS of circulating total-tau levels in European-ancestry participants (N = 14,721) or African American participants (N = 953).</caption><thead><tr><th rowspan="3">rsid</th><th rowspan="3">Chr637Pos</th><th rowspan="3">Alleles</th><th colspan="10">Europeans (N = 14,721)</th><th colspan="4">Multi-ancestry</th></tr><tr><th colspan="2">SE</th><th colspan="2">Beta</th><th colspan="2">EAF</th><th colspan="2">P</th><th colspan="2">P<sub>Q</sub></th><th colspan="2">P<sub>RE2</sub></th><th colspan="2">Gene</th></tr><tr><th>P</th><th>SE</th><th>P</th><th>SE</th><th>P</th><th>SE</th><th>P</th><th>SE</th><th>P</th><th>SE</th><th>P</th><th>SE</th><th>P</th><th>SE</th></tr></thead><tbody><tr><td>rs111836296*</td><td>4142551581</td><td>T/C</td><td>--</td><td>--</td><td>--</td><td>--</td><td>0.06</td><td>-0.54</td><td>0.10</td><td>1.7E-08</td><td>0</td><td>0.37</td><td>--</td><td>--</td><td>intergenic</td></tr><tr><td>rs74710969*</td><td>523681842</td><td>T/G</td><td>--</td><td>--</td><td>--</td><td>0.11</td><td>-0.53</td><td>0.10</td><td>3.4E-08</td><td>0</td><td>0.88</td><td>--</td><td>--</td><td>ADAMTS2</td></tr><tr><td>rs674432</td><td>615801153</td><td>C/G</td><td>0.83</td><td>-0.001</td><td>0.01</td><td>0.97</td><td>0.68</td><td>0.12</td><td>1.3E-08</td><td>0</td><td>0.32</td><td>7.6E-06</td><td>96.84</td><td>1.8E-08</td></tr><tr><td>rs7502280*</td><td>1743670221</td><td>T/G</td><td>0.88</td><td>0.17</td><td>0.01</td><td>--</td><td>0.45</td><td>0</td><td>0</td><td>--</td><td>--</td><td>0.04</td><td>1.8E-142</td><td>76.46</td><td>0.04</td></tr><tr><td>rs29257</td><td>174499712</td><td>A/G</td><td>0.88</td><td>0.20</td><td>0.01</td><td>0.28</td><td>0.05</td><td>0.07</td><td>0.50</td><td>76.2</td><td>0</td><td>1.00</td><td>1.5E-77</td><td>0.00</td><td>0.57</td></tr><tr><td>rs2942003</td><td>1744976704</td><td>T/G</td><td>0.34</td><td>0.16</td><td>0.01</td><td>0.86</td><td>0.09</td><td>0.13</td><td>0.51</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.00</td><td>intergenic</td></tr></tbody></table>

EAF effect allele frequency, Alleles: effect (Alternate) allele/non-effect (Reference) allele.

P<sub>RE2</sub>: I-square heterogeneity statistic; P<sub>Q</sub>: Cochran's Q statistic; P value.

rs111836296 on chr4 and rs74710969 on chr5 are extremely rare in Europeans and thus they were not included in the multi-ancestry meta-analysis; rs7502280 is not present in the 1000 Genomes reference panel and thus was not present in the African American GWAS results.

P<sub>RE2</sub>: Han and Eskin's random effects model to detect associations under heterogeneity⁵⁹.

P<sub>RE2</sub>: Han and Eskin's random effects model to detect associations under heterogeneity⁵⁹.

organization, crucial for neuronal polarization⁴⁵. Moreover, atypical protein kinase C (aPKC) in complex with PAR-3/PAR-6 negatively regulates microtubule affinity-regulating kinases, which in turn causes dephosphorylation of microtubule-associated proteins, such as tau, leading to the assembly of microtubules and elongation of axons⁴⁶. Par3 also regulates APP processing and trafficking⁴⁷, polarized convergence between APP and BACE1 in hippocampal neurons⁴⁸, and retrograde endosome-to- trans-Golgi network trafficking of BACE1 along with aPKC⁴⁹. Brain regulatory marks (promoter and enhancer) are reported at the lead variant rs12245909 (HaploReg v4.1), which may suggest a functional role of this variant in the brain.

We looked up the main distinct lead genetic variants from the published ADNI GWAS of circulating tau levels¹⁸ in our European meta-analysis (excluding ADNI). We confirmed the strong association of the MAPT rs242557-A genetic variant (Table 5 and Supplementary Table 12). However, we did not find evidence of association for the three other loci that were detected at $P < 10^{-5}$ in the original ADNI GWAS, suggesting that these signals may have been false positives.

Among the 10 genes identified when leveraging whole exome sequence data and aggregating rare variants with high or moderate impact, four have a function relevant to the brain (ELFN2, UBASH3B, SLIT3, and NSD3). Interestingly, SLIT3 and NSD3 associations were more impacted by the choice of the MAF threshold to select rare variants to aggregate. For both genes, results were only significant with SKAT when using a MAF ≤ 1%, indicating that only rarer variations were contributing to the associations. The extracellular leucine rich repeat and fibronectin type III domain containing 2 gene (ELFN2), is overexpressed in the brain. The encoded protein is a postsynaptic adhesion molecule that selectively binds with group III metabotropic glutamate receptors⁵⁰,⁵¹. ELFN1, a protein of the same family, has been reported to be associated with neuropsychiatric disorders (attention deficit hyperactivity disorder, post-traumatic stress disorder, and epilepsy). Distinct neuronal expression patterns are reported for ELFN1 and ELFN2⁵¹. The ubiquitin associated and SH3 domain containing B gene (UBASH3B) is overexpressed in the brain. The encoded protein is a phosphatase, and the concerted action of protein kinases and phosphatases represents a critical signaling event controlling synaptic functions and higher-order brain functions, such as learning and memory⁵². The slit guidance ligand 3 gene (SLIT3) encodes an axon guidance molecule expressed by motor neurons⁵³,⁵⁴. SLIT3 may also play a role in essential tremor disease pathogenesis⁵⁵. The nuclear receptor binding SET domain protein 3 gene (NSD3) is highly expressed in the brain. The encoded protein is a SET domain-containing methyltransferase, an epigenetic regulator that is selectively expressed in primary microglia⁵⁶. Follow-up studies are needed to characterize the potential role of these four genes in tauopathies.

A summary of the neurological traits reported in the GWAS catalog for genetic variants in the main genes identified in the meta-analyses of circulating t-tau levels (IL15, FHIT, ADAMTS12, PARD3, F5, BCAS3, UBASH3B, and SLIT3) and described above is available in Supplementary Table 11.

By performing meta-analyses separately in African Americans and European-ancestry participants, we were able to identify ancestry-specific associations for circulating t-tau levels. The lead variants at the three loci identified at the genome-wide threshold in African American participants were extremely rare in European populations. In addition, most loci identified in African American participants were driven by the largest study, ARIC. Two of the findings identified at the genome-wide threshold in African Americans were low frequency variants, with no linkage disequilibrium support (Supplementary Figs. 15, 18), and with
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Regional association plot at the MAPT 17q21 locus. Genetic variants are plotted with their *P* values (-log10 values, left y axis) as a function of the build 37 genomic position (x axis). Estimated recombination rates (right y axis) reflect the local linkage disequilibrium (LD) structure around the top distinct associated genetic variants (red and blue diamonds, denoting rs242557 and rs2942003, respectively), identified using stepwise model selection procedure, and their correlated proxies. Genetic variants in LD with rs242557 are indicated with circles according to a light to dark red scale, from r² = 0 to 1. Genetic variants in LD with rs2942003 are indicated with triangles according to a light to dark blue scale, from r² = 0 to 1. Gray crosses (X) represent SNPs with missing LD. The third distinct genetic variant (rs7502280), indicated in black on the plot, was not present in the reference panel to calculate the LD with the other variants in the region. Reference panel used was the 1000 Genomes November 2014 European population.</figcaption></figure>

<table><caption>**Table 4** Results of the APOE4-stratified analyses for the lead genetic variants in each locus passing the threshold of *P* < 5 × 10⁻⁷ in the European meta-analysis of GWAS of circulating total-tau levels.</caption><thead><tr><th rowspan="2">rsid</th><th rowspan="2">Chr</th><th rowspan="2">Build 37 Pos (bp)</th><th rowspan="2">Eff</th><th rowspan="2">NEff</th><th rowspan="2">EAF</th><th colspan="2">Main analysis<br>(N = 14,721)</th><th colspan="2">APOE4 carriers<br>(N = 3640)</th><th colspan="2">APOE4 Non-<br>carriers<br>(N = 10,574)</th><th rowspan="2">Gene</th></tr><tr><th>Beta</th><th>P</th><th>Beta</th><th>P</th><th>Beta</th><th>P</th></tr></thead><tbody><tr><td>rs6686805</td><td>1</td><td>169,512,643</td><td>A</td><td>C</td><td>0.67</td><td>-0.04</td><td>3.4E-07</td><td>-0.02</td><td>0.21</td><td>-0.04</td><td>5.1E-07</td><td>F5</td></tr><tr><td>rs139843727</td><td>3</td><td>66,316,022</td><td>A</td><td>C</td><td>0.99</td><td>0.22</td><td>2.9E-07</td><td>0.08</td><td>0.44</td><td>0.24</td><td>6.5E-07</td><td>SLC25A26</td></tr><tr><td>rs7502280</td><td>17</td><td>43,670,221</td><td>T</td><td>G</td><td>0.88</td><td>0.17</td><td>8.0E-38</td><td>0.16</td><td>2.2E-09</td><td>0.18</td><td>1.4E-30</td><td>intergenic</td></tr><tr><td>rs242557</td><td>17</td><td>44,019,712</td><td>A</td><td>G</td><td>0.38</td><td>0.20</td><td>8.9E-143</td><td>0.19</td><td>1.3E-37</td><td>0.20</td><td>1.4E-105</td><td>MAPT</td></tr><tr><td>rs2942003</td><td>17</td><td>44,576,704</td><td>T</td><td>G</td><td>0.34</td><td>0.16</td><td>1.9E-78</td><td>0.17</td><td>3.5E-22</td><td>0.16</td><td>3.4E-58</td><td>intergenic</td></tr><tr><td>rs4968553</td><td>17</td><td>59,428,962</td><td>C</td><td>G</td><td>0.16</td><td>-0.05</td><td>4.6E-07</td><td>-0.05</td><td>0.02</td><td>-0.06</td><td>1.7E-07</td><td>BCAS3</td></tr></tbody></table>

EAF effect allele frequency, *Eff* effect (alternate) allele, *Neff* non-effect (reference) allele.

The association of the two SNPs defining APOE in the main European meta-analysis were: rs429358-T (Beta = -0.02, *P* = 0.10) and rs7412-T (Beta = 0.01, *P* = 0.49).

only two of the three cohorts that contributed to the meta-analysis. Caution is thus needed regarding the interpretation of these findings as such results are typically seen in GWAS of admixed populations with small sample sizes and could be driven by a few outliers. We also found that the strong association of the MAPT locus with circulating t-tau levels was specific to European-ancestry participants. This result is consistent with the recent finding from the Florida Consortium for African American Alzheimer's Disease studies⁵⁷. The majority of loci identified in European-ancestry participants were driven by the two largest studies, FHS and RSI. The fact that we did not detect an

association in the African American participants for the novel loci detected at *P* < 5 × 10⁻⁷ in the European-ancestry participants may be due to a lack of power because of the limited sample size of this subgroup. However, our multi-ancestry meta-analysis showed that the hits identified were ancestry specific (Supplementary Table 3). The high heterogeneity observed across ancestry suggests that the genetic architecture of circulating t-tau levels may differ between European-ancestry and African American populations.

We explored the potential pleiotropy of loci previously reported for several neurological disorders with circulating t-tau
<table><caption>Table 5 Look-up of the main hits (P < 10⁻⁵) from the published ADNI GWAS of circulating tau levels (Chen et al., 2017) in our European meta-analysis of GWAS of circulating total-tau levels (seven studies excluding ADNI).</caption><thead><tr><th rowspan="2">rsid</th><th rowspan="2">Chr</th><th rowspan="2">Build 37 Pos (bp)</th><th rowspan="2">Eff</th><th rowspan="2">NEff</th><th rowspan="2">EAF</th><th rowspan="2">Main analysis (N = 14,235)</th><th colspan="2">APOE4 carriers (N = 3401)</th><th colspan="2">APOE4 Non-carriers (N = 10,327)</th><th rowspan="2">Gene</th></tr><tr><th>Beta</th><th>P</th><th>Beta</th><th>P</th><th>Beta</th><th>P</th></tr></thead><tbody><tr><td>rs2187213</td><td>6</td><td>162,634,337</td><td>A</td><td>G</td><td>0.35</td><td>0.0005</td><td>0.95</td><td>0.005</td><td>0.77</td><td>-0.0006</td><td>0.95</td><td>PARK2</td></tr><tr><td>rs7047280</td><td>9</td><td>23,297,808</td><td>T</td><td>C</td><td>0.60</td><td>-0.0009</td><td>0.90</td><td>-0.008</td><td>0.56</td><td>0.003</td><td>0.71</td><td>ELAVL2</td></tr><tr><td>rs7072793</td><td>10</td><td>6,106,266</td><td>T</td><td>C</td><td>0.59</td><td>0.008</td><td>0.22</td><td>0.02</td><td>0.17</td><td>0.002</td><td>0.82</td><td>IL2RA</td></tr><tr><td>rs242557</td><td>17</td><td>44,019,712</td><td>A</td><td>G</td><td>0.36</td><td>0.20</td><td>6.4E-136</td><td>0.19</td><td>3.0E-35</td><td>0.20</td><td>1.4E-101</td><td>MAPT</td></tr></tbody><tfoot><tr><td colspan="12">EAF effect allele frequency: Eff effect (alternate) allele, Neff non-effect (reference) allele.</td></tr></tfoot></table>

levels by performing a look-up of these loci in our European meta-analysis. We also used MR analyses to evaluate the potential causal associations between circulating t-tau levels with several neurological disorders and traits, but we did not identify significant causal associations. Limitations of these analyses are the availability of large, complete, and publicly available GWAS summary statistics, and the strength of the MR instruments due to the limited number of associations in the European circulating t-tau meta-analysis and the limited number of AD and PD cases.

Strengths of our study are the large sample size, with the inclusion of eight population-based cohorts with ancestral diversity, with genotype data and information on circulating t-tau levels measured with ultra-sensitive assays. We leveraged large imputations reference panels (1000 Genomes and the Haplotype Reference Consortium) to study common genetic variations complemented with whole exome sequence data to explore less frequent genetic variations. Some limitations include the modest sample size of the African American sample that has limited our ability to confirm the GWAS findings and perform secondary analyses in this subgroup, such as stratification on *APOE4* status, and the fact that the contributed studies only had circulating t-tau measurement available, and not specifically phosphorylated tau. The fact that the *APOE* locus was not associated with circulating t-tau levels may suggest that circulating t-tau and phosphorylated tau have different genetic architectures. Replication of the African American potential novel loci and the less common variant association identified in Europeans is needed to confirm our findings but the availability of samples with circulating t-tau levels and genetic data is limited.

In conclusion, our large multi-ancestry meta-analysis identified new genetic variants and loci associated with circulating t-tau levels. Notably, our study revealed that the genetic architecture underlying circulating t-tau levels might differs between African American and European-ancestry populations and that genetic variation underlying circulating t-tau levels may overlap with known genetic determinants of neurological disorders. To better understand how these variants may contribute to AD and other tauopathies, further investigations of these findings will be necessary, including cohorts with a broader ancestral diversity, biological experiments, functional and omic studies, and animal models.

## Methods

**Populations and participants.** We included in our multi-ancestry and ancestry-specific meta-analyses of total-tau participants from eight studies: seven cohorts from the Cohorts for heart and aging research in genomic epidemiology (CHARGE) consortium (the Framingham Heart Study (FHS), the Rotterdam Study (RSI and RSII), the MEMENTO Study, the Coronary Artery Risk Development in Young Adults (CARDIA) Study, the Cardiovascular Health Study (CHS), the Vietnam Era Twin Study of Aging (VETSA) Study, and the Atherosclerosis Risk in Communities (ARIC) Study), and the ADNI Study. All participants included in this study provided written informed consent for genetic testing and analyses. Study-specific information including study description, and detailed information about genotyping and imputations and GWAS analysis is included in the Supplementary Notes 1–8.

**Tau quantification.** Circulating (in plasma or in serum) t-tau levels were quantified using the Human Total Tau kit on the Simoa™ HD-1 analyzer (ADNI, plasma), the Simoa™ Tau 2.0 Kit and the Simoa™ HD-1 analyzer (FHS, plasma), the Simoa™ Tau 2.0 Kit (MEMENTO_1, plasma), the Simoa™ Human Neurology 4-Plex A assay with the Simoa™ HD-X analyzer (the Atherosclerosis Risk in Communities study, MEMENTO_2 and the Coronary Artery Risk Development in Young Adults, plasma and the Cardiovascular Health Study, serum), the Simoa™ Human Neurology 3-Plex A assay with the Simoa™ HD-1 analyzer (the Rotterdam Study, plasma) and high throughput bioassays platforms or single analyte assays using the Simoa™ HD-X or Fujirebio analyzer (the Vietnam Era Twin Study of Aging, plasma).

**Genotyping and imputation.** Studies used the densest imputation reference panel available to them at the time of analyses, either 1000 Genomes or the Haplotype
Reference Consortium. A description of the reference panel used by each study is provided in Supplementary Tables 13–15.

**Genome-wide association studies (GWAS) and quality control.** Each study evaluated the association of single-nucleotide genetic variants with log2-transformed circulating t-tau levels under an additive model. Analyses were adjusted for age, sex, and additional study specific covariates to control for population structure (including principal components). Studies with both Europeans and African Americans analyzed each ancestry separately. The minimum sample size for each ancestry/phenotype combination for inclusion in this study was fixed to 100. Additional stratified analyses according to *APOE4* carriers vs non *APOE4* carriers) were performed in European participants only, given the modest sample size of the African American sample. Study-specific GWAS results were filtered based on an imputation quality score greater or equal to 0.30 and a minor allele count greater or equal to 20. The minor allele frequency threshold to include a genetic variant in the meta-analyses is indicated for each study in Supplementary Tables 13–15.

**Multi-ancestry and ancestry-specific meta-analyses.** GWAS results across all studies were meta-analysed by ancestry with METAL using an inverse variance-weighted average method⁵⁸. For each meta-analysis, we selected results for which two third of the studies contributed. Results from the European and the African American meta-analyses, and for which at least two studies contributed to each meta-analysis, were then meta-analysed using Metasoft⁵⁹,⁶⁰.

**Conditional and joint association analysis based on summary statistics.** Conditional and joint association analysis in loci associated with circulating t-tau levels at the genome-wide threshold ($P < 5 \times 10^{-8}$) were performed using the Genome-wide Complex Trait Analysis⁶¹,⁶² based on the European meta-analysis results. A stepwise model selection procedure was used to select distinct associated genetic variants ($P < 5 \times 10^{-8}$). The FHS Haplotype Reference Consortium impu- ted data was used as the reference panel with unrelated participants only.

**Overlap of circulating t-tau genetic determinants with neurological diseases and traits.** We extracted from the GWAS Catalog²⁰ (https://www.ebi.ac.uk/gwas/downloads/summary-statistics) all reported associations for AD, t-tau (in plasma or CSF), stroke, PSP, Parkinson's Disease (PD), and White Matter Hyperintensities (WMH), and looked them up in our European meta-analysis results. We also used the FUMA GWAS platform (Functional Mapping and Annotation of Genome- Wide Association Studies, https://fuma.ctglab.nl/) using as input the summary statistics from the European circulating t-tau meta-analysis to leverage functional, biological information to prioritize genes and check expression patterns and shared molecular functions between these genes⁶³,⁶⁴. Finally, we performed a GRS asso- ciation analysis in the FHS. We first re-ran the European circulating t-tau levels meta-analysis without FHS. We then used the Genome-wide Complex Trait Analysis to identify distinct genome-wide associations on chromosome 17. We computed the GRS using the distinct variants identified by the Genome-wide Complex Trait Analysis for all FHS participants using the Haplotype Reference Consortium imputed genotypes and weights from the new meta-analysis. We tested the association of the GRS with incident AD (140 cases, 2775 controls) and stroke (149 cases, 3461 controls), and with four brain MRI phenotypes (hippo- campal volume, white matter hyperintensities, total brain volume, and intracranial volume, N ranging from 3489 to 4310). Details regarding trait measurements and definitions in the FHS have been published elsewhere⁶⁵⁻⁶⁷. We used logistic or linear mixed-effects model, adjusted for age at baseline or at MRI and sex, while accounting for relatedness. For the brain MRI analyses, we excluded participants with dementia, stroke, large brain infarcts, tumor or any other finding that could have affected the scan and additionally adjusted the hippocampal volume, white matter hyperintensities, and total brain volume analyses for intracranial volume. If a significant association was detected ($P = 0.05/6 = 0.008$), an additional adjust- ment for *APOE4* was performed.

**Two sample mendelian randomization (MR) analyses.** We used the TwoSam- pleMR R package in R⁶⁸,⁶⁹ to assess the causal association of circulating t-tau levels with AD, PD, stroke, and WMH using publicly available large European GWAS summary statistics. We selected large GWAS with independent samples from the ones included our meta-analysis. Briefly, we used an arbitrary threshold of $P < 5 \times 10^{-6}$ in the European meta-analysis of circulating t-tau levels to select the genetic variants to be used in the MR and performing clumping to select distinct variants based on 1000 Genomes European linkage disequilibrium reference panel. We then extracted these SNPs from the outcome GWAS based on summary sta- tistics from the GWAS catalog (https://www.ebi.ac.uk/gwas/downloads/summary- statistics) or the IEU GWAS database (IGD, https://gwas.mrcieu.ac.uk/) and per- formed harmonization of the alleles. As MAPT locus is known to be pleiotropic, we conducted the mendelian randomization using different methods and carefully check for presence of heterogeneity and horizontal pleiotropy. For this analysis, we conducted power calculations for a continuous exposure and a PVE on the exposure of 8%, using the power analysis calculator https://sb452.shinyapps.io/

/power/. Finally, we tested the opposite hypothesis that AD, PD, stroke, and WMH are causally associated with circulating t-tau levels.

**Rare variant analyses using whole exome sequence data.** To explore the association of rare variations with circulating t-tau levels, we selected the two largest studies (FHS and RSI) to perform rare-variant aggregation tests based on whole exome sequence data from the Cohorts for heart and aging research in genomic epidemiology (CHARGE) consortium⁷⁰. A total of 2279 participants were included in the analyses. Information on sequencing and quality control is provided in the Supplementary Notes 9 and 10. To make sure that the same allele was coded as the effect allele for FHS and RSI, we used the effect (alternate) allele from a consensus SNP info file from the Cohorts for heart and aging research in genomic epidemiology (CHARGE) consortium⁷¹. Annotations of the exome variants was performed with dbNSFP⁷². We selected variants with a MAF ≤ 1% or 5%, and (1) missense and loss of function variants only or (2) variants with high or moderate impact from Ensembl Variant Effect Predictor⁷³ including missense and loss of function variants. The analyses were performed using the R package seqMeta (http://cran.r-project.org/web/packages/seqMeta/). Each cohort used the seqMeta prepScores function to generate single variant score statistics and genotype covariance matrices for all variants. Results were then meta- analyzed using the skatMeta and burdenMeta functions. We used a Bonferroni correction for the number of genes included in the analyses and the number of tests ($P = 0.05/20,000/2 = 1.25 \times 10^{-6}$) and filtered the results based on a cumulative minor allele count of 30, that accounts for the number of genetic variants per gene.

**Reporting summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## Data availability

All data supporting the findings of this study are available either within the main article or the supplementary information. Summary statistics from the ancestry-specific meta- analyses of circulating levels of total-tau have been deposited and are publicly accessible on the GWAS catalog FTP (study accession numbers GCST90095138, and GCST90095139). Genome-wide summary statistics for complex disorders used in the secondary analyses were downloaded from public repositories (GWAS catalog: https:// www.ebi.ac.uk/gwas/downloads/summary-statistics; and IEU GWAS database: https:// gwas.mrcieu.ac.uk/).

Received: 3 August 2021; Accepted: 17 March 2022;
Published online: 08 April 2022

## References

1. Schraen-Maschke, S. et al. Tau as a biomarker of neurodegenerative diseases. *Biomark. Med.* **2**, 363–384 (2008).

2. Blennow, K. & Hampel, H. CSF markers for incipient Alzheimer's disease. *Lancet Neurol.* **2**, 605–613 (2003).

3. Olsson, B. et al. CSF and blood biomarkers for the diagnosis of Alzheimer's disease: a systematic review and meta-analysis. *Lancet Neurol.* **15**, 673–684 (2016).

4. Zwan, M. D. et al. Use of amyloid-PET to determine cutpoints for CSF markers: a multicenter study. *Neurology* **86**, 50–58 (2016).

5. Wagshal, D. et al. Divergent CSF tau alterations in two common tauopathies: Alzheimer's disease and progressive supranuclear palsy. *J. Neurol. Neurosurg. Psychiatry* **86**, 244–250 (2015).

6. Rosso, S. M. et al. Total tau and phosphorylated tau 181 levels in the cerebrospinal fluid of patients with frontotemporal dementia due to P301L and G272V tau mutations. *Arch. Neurol.* **60**, 1209–1213 (2003).

7. Zetterberg, H. & Burnham, S. C. Blood-based molecular biomarkers for Alzheimer's disease. *Mol. Brain* **12**, 26-019-0448-1 (2019).

8. Nation, D. A. et al. Blood-brain barrier breakdown is an early biomarker of human cognitive dysfunction. *Nat. Med.* **25**, 270–276 (2019).

9. Deane, R., Bell, R. D., Sagare, A. & Zlokkovic, B. V. Clearance of amyloid-beta peptide across the blood-brain barrier: implication for therapies in Alzheimer's disease. *CNS Neurol. Disord. Drug Targets* **8**, 16–30 (2009).

10. Pase, M. P. et al. Assessment of Plasma Total Tau Level as a predictive biomarker for dementia and related endophenotypes. *JAMA Neurol.* **76**, 598–606 (2019).

11. Dage, J. L. et al. Levels of tau protein in plasma are associated with neurodegeneration and cognitive function in a population-based elderly cohort. *Alzheimers Dement* **12**, 1226–1234 (2016).

12. Zetterberg, H. et al. Plasma tau levels in Alzheimer's disease. *Alzheimers Res. Ther.* **5**, 9 (2013).

13. Mattsson, N. et al. Plasma tau in Alzheimer disease. *Neurology* **87**, 1827–1835 (2016).
14. Mielke, M. M. et al. Association of Plasma Total Tau Level with cognitive decline and risk of mild cognitive impairment or dementia in the mayo clinic study on aging. *JAMA Neurol.* **74**, 1073–1080 (2017).

15. Pase, M. P. et al. Plasma total-tau as a biomarker of stroke risk in the community. *Ann. Neurol.* **86**, 463–467 (2019).

16. Fossati, S. et al. Plasma tau complements CSF tau and P-tau in the diagnosis of Alzheimer's disease. *Alzheimers Dement. (Amst.)* **11**, 483-492 (2019).

17. Nam, E., Lee, Y. B., Moon, C. & Chang, K. A. Serum Tau proteins as potential biomarkers for the assessment of Alzheimer's disease progression. *Int. J. Mol. Sci.* **21**, https://doi.org/10.3390/ijms21145007 (2020).

18. Chen, J. et al. Genome-wide association study identifies MAPT locus influencing human plasma tau levels. *Neurology* **88**, 669–676 (2017).

19. Lord, J. et al. A genome-wide association study of plasma phosphorylated tau181. *Neurobiol. Aging* **106**, 304.e1–304.e3 (2021).

20. MacArthur, J. et al. The new NHGRI-EBI Catalog of published genome-wide association studies (GWAS Catalog). *Nucleic Acids Res.* **45**, D896–D901 (2017).

21. Hoglinger, G. U. et al. Identification of common variants influencing risk of the tauopathy progressive supranuclear palsy. *Nat. Genet.* **43**, 699–705 (2011).

22. Jun, G. et al. A novel Alzheimer disease locus located near the gene encoding tau protein. *Mol. Psychiatry* **21**, 108–117 (2016).

23. Adams, H. H. et al. Novel genetic loci underlying human intracranial volume identified through genome-wide association. *Nat. Neurosci.* **19**, 1569–1582 (2016).

24. Hibar, D. P. et al. Common genetic variants influence human subcortical brain structures. *Nature* **520**, 224–229 (2015).

25. Shin, J. et al. Global and regional development of the human cerebral cortex: molecular architecture and occupational aptitudes. *Cereb. Cortex* **30**, 4121–4139 (2020).

26. Zhang, X. et al. Tau pathology in Parkinson's disease. *Front. Neurol.* **9**, 809 (2018).

27. Meddens, S. F. W. et al. Genomic analysis of diet composition finds novel loci and associations with health and lifestyle. *Mol. Psychiatry* **26**, 2056–2069 (2021).

28. Doherty, A. et al. GWAS identifies 14 loci for device-measured physical activity and sleep duration. *Nat. Commun.* **9**, 5257 (2018).

29. Latourelle, J. C., Dumitriu, A., Hadzi, T. C., Beach, T. G. & Myers, R. H. Evaluation of Parkinson disease risk variants as expression-QTLs. *PLoS ONE* **7**, e46199 (2012).

30. Bishnoi, R. J., Palmer, R. F. & Royall, D. R. Serum interleukin (IL)-15 as a biomarker of Alzheimer's disease. *PLoS ONE* **10**, e0117282 (2015).

31. Fontanil, T. et al. Neurocan is a new substrate for the ADAMTS12 metalloprotease: potential implications in neuropathies. *Cell. Physiol. Biochem.* **52**, 1003–1016 (2019).

32. Arning, A. et al. A genome-wide association study identifies a gene network of ADAMTS genes in the predisposition to pediatric stroke. *Blood* **120**, 5231–5236 (2012).

33. Witten, A. et al. ADAMTS12, a new candidate gene for pediatric stroke. *PLoS ONE* **15**, e0237928 (2020).

34. Gevorkian, G. et al. Amyloid-beta peptide binds to microtubule-associated protein 1B (MAP1B). *Neurochem. Int.* **52**, 1030–1036 (2008).

35. Hasegawa, M., Arai, T. & Ihara, Y. Immunochemical evidence that fragments of phosphorylated MAP5 (MAP1B) are bound to neurofibrillary tangles in Alzheimer's disease. *Neuron* **4**, 909–918 (1990).

36. Jensen, P. H. et al. Microtubule-associated protein 1B is a component of cortical Lewy bodies and binds alpha-synuclein filaments. *J. Biol. Chem.* **275**, 21500–21507 (2000).

37. Joshi, D. & Inamdar, M. S. Rudhira/BCAS3 couples microtubules and intermediate filaments to promote cell migration for angiogenic remodeling. *Mol. Biol. Cell* **30**, 1437–1450 (2019).

38. Cheng, Y. C. et al. Genome-Wide Association analysis of young-onset stroke identifies a locus on chromosome 10q25 Near HABP2. *Stroke* **47**, 307–316 (2016).

39. Sun, B. B. et al. Genomic atlas of the human plasma proteome. *Nature* **558**, 73–79 (2018).

40. Hinds, D. A. et al. Genome-wide association analysis of self-reported events in 6135 individuals and 252 827 controls identifies 8 loci associated with thrombosis. *Hum. Mol. Genet.* **25**, 1867–1874 (2016).

41. Logue, M. W. et al. Targeted Sequencing of Alzheimer disease genes in African Americans implicates novel risk variants. *Front. Neurosci.* **12**, 592 (2018).

42. Melville, S. A. et al. Multiple loci influencing hippocampal degeneration identified by genome scan. *Ann. Neurol.* **72**, 65–75 (2012).

43. Khazaei, M. R. & Puschel, A. W. Phosphorylation of the par polarity complex protein Par3 at serine 962 is mediated by aurora a and regulates its function in neuronal polarity. *J. Biol. Chem.* **284**, 33571–33579 (2009).

44. Chen, X. et al. Rare Deleterious PARD3 Variants in the aPKC-Binding Region are implicated in the pathogenesis of human cranial neural tube defects via disrupting apical tight junction formation. *Hum. Mutat.* **38**, 378–389 (2017).

45. Chen, S. et al. Regulation of microtubule stability and organization by mammalian Par3 in specifying neuronal polarity. *Dev. Cell.* **24**, 26–40 (2013).

46. Chen, Y. M. et al. Microtubule affinity-regulating kinase 2 functions downstream of the PAR-3/PAR-6/atypical PKC complex in regulating hippocampal neuronal polarity. *Proc. Natl Acad. Sci. U. S. A.* **103**, 8534–8539 (2006).

47. Sun, M., Asghar, S. Z. & Zhang, H. The polarity protein Par3 regulates APP trafficking and processing through the endocytic adaptor protein Numb. *Neurobiol. Dis.* **93**, 1–11 (2016).

48. Sun, M., Huang, C., Wang, H. & Zhang, H. Par3 regulates polarized convergence between APP and BACE1 in hippocampal neurons. *Neurobiol. Aging* **77**, 87–93 (2019).

49. Sun, M. & Zhang, H. Par3 and aPKC regulate BACE1 endosome-to-TGN trafficking through PACS1. *Neurobiol. Aging* **60**, 129–140 (2017).

50. Dunn, H. A., Zucca, S., Dao, M., Orlandi, C. & Martemyanov, K. A. ELFN2 is a postsynaptic cell adhesion molecule with essential roles in controlling group III mGluRs in the brain and neuropsychiatric behavior. *Mol. Psychiatry* **24**, 1902–1919 (2019).

51. Matsunaga, H. & Aruga, J. Trans-Synaptic regulation of metabotropic glutamate receptors by Elfn proteins in health and disease. *Front. Neural Circuits* **15**, 634875 (2021).

52. Smidak, R. et al. Quantitative proteomics reveals protein kinases and phosphatases in the individual phases of contextual fear conditioning in the C57BL/6J mouse. *Behav. Brain Res.* **303**, 208–217 (2016).

53. Punnamoottil, B., Rinkwitz, S., Giacomotto, J., Svahn, A. J. & Becker, T. S. Motor neuron-expressed microRNAs 218 and their enhancers are nested within introns of Slit2/3 genes. *Genesis* **53**, 321–328 (2015).

54. Brose, K. et al. Slit proteins bind Robo receptors and have an evolutionarily conserved role in repulsive axon guidance. *Cell* **96**, 795–806 (1999).

55. Odgerel, Z. et al. Whole genome sequencing and rare variant analysis in essential tremor families. *PLoS ONE* **14**, e0220512 (2019).

56. Das, A. et al. Transcriptome sequencing reveals that LPS-triggered transcriptional responses in established microglia BV2 cell lines are poorly representative of primary microglia. *J. Neuroinflammation* **13**, 182 (2016).

57. Deniz, K. et al. Plasma Biomarkers of Alzheimer's Disease in African Americans. *J. Alzheimers Dis.* **79**, 323–334 (2021).

58. Willer, C. J., Li, Y. & Abecasis, G. R. METAL: fast and efficient meta-analysis of genomewide association scans. *Bioinformatics* **26**, 2190–2191 (2010).

59. Han, B. & Eskin, E. Random-effects model aimed at discovering associations in meta-analysis of genome-wide association studies. *Am. J. Hum. Genet.* **88**, 586–598 (2011).

60. Kang, E. Y. et al. ForestPMPlot: a flexible tool for visualizing heterogeneity between studies in meta-analysis. *G3 (Bethesda)* **6**, 1793–1798 (2016).

61. Yang, J., Lee, S. H., Goddard, M. E. & Visscher, P. M. GCTA: a tool for genome-wide complex trait analysis. *Am. J. Hum. Genet.* **88**, 76–82 (2011).

62. Yang, J. et al. Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits. *Nat. Genet.* **44**, 369–375 (2012).

63. Watanabe, K., Taskesen, E., van Bochoven, A. & Posthuma, D. Functional mapping and annotation of genetic associations with FUMA. *Nat. Commun.* **8**, 1826 (2017).

64. Watanabe, K., Umicevic Mirkov, M., de Leeuw, C. A., van den Heuvel, M. P. & Posthuma, D. Genetic mapping of cell type specificity for complex traits. *Nat. Commun.* **10**, 3222 (2019).

65. Malik, R. et al. Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes. *Nat. Genet.* **50**, 524–537 (2018).

66. Satizabal, C., Beiser, A. S. & Seshadri, S. Incidence of dementia over three decades in the Framingham Heart Study. *N. Engl. J. Med.* **375**, 93–94 (2016).

67. Sarnowski, C. et al. Whole genome sequence analyses of brain imaging measures in the Framingham Study. *Neurology* **90**, e188–e196 (2018).

68. Hemani, G., Tilling, K. & Davey Smith, G. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. *PLoS Genet.* **13**, e1007081 (2017).

69. Hemani, G. et al. The MR-Base platform supports systematic causal inference across the human phenome. *Elife* **7**, https://doi.org/10.7554/eLife.34408 (2018).

70. Psaty, B. M. et al. Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium: Design of prospective meta-analyses of genome-wide association studies from 5 cohorts. *Circ. Cardiovasc. Genet.* **2**, 73–80 (2009).

71. Grove, M. L. et al. Best practices and joint calling of the HumanExome BeadChip: the CHARGE Consortium. *PLoS ONE* **8**, e68095 (2013).

72. Liu, X., Jian, X. & Boerwinkle, E. dbNSFP: a lightweight database of human nonsynonymous SNPs and their functional predictions. *Hum. Mutat.* **32**, 894–899 (2011).

73. McLaren, W. et al. The ensembl variant effect predictor. *Genome Biol.* **17**, 122 (2016).
## Acknowledgements

We thank all the participating cohorts from the CHARGE consortium and ADNI participants for contributing to this study. Detailed study-specific acknowledgements are provided in the Supplementary Notes 11 and 12.

## Author contributions

C.S., E.C.P.F., J.C.B., M.L., M.F., and A.M. conducted the cohort specific GWAS analysis. C.S., S.A., and E.C.P.F. conducted the exome sequence analyses. A.Y., V.C., and V.D. performed secondary analyses. V.B., and R.A.R. were involved in cohort-specific tau quantification. C.S., M.G., J.C.B., M.F., A.S.B., E.B., L.A.C., J.C.L., C.S.D., A.L.D., L.D., T.H.M., K.J.M., M.P.P., C.L.S., R.S.V., S.D., C.E.F., M.J.L., C.D., L.J.L., W.S.K., W.T.L., M.A.I., and S.S. designed and conceived the study. C.S., M.G., M.F., S.A., C.S.D., A.E.F., T.F.K., K.J.M., M.P.P., C.L.S., S.D., W.S.K., W.T.L., M.A.I., and S.S. wrote and edited the manuscript. All authors reviewed and approved the manuscript.

## Competing interests

C.S.D. reports personal fees from Novartis Pharmaceutical outside the submitted work. All the other authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s42003-022-03287-y.

**Correspondence** and requests for materials should be addressed to Chloé. Sarnowski.

**Peer review information** *Communications Biology* thanks the anonymous reviewers for their contribution to the peer review of this work. Primary Handling Editor: George Inglis. Peer reviewer reports are available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022

¹Department of Epidemiology, Human Genetics, and Environmental Sciences, University of Texas Health Science Center at Houston, Houston, TX, USA. ²Department of Epidemiology, Erasmus University Medical Center, Rotterdam, The Netherlands. ³Department of Genetics, School of Medicine, Mashhad University of Medical Sciences, Mashhad, Iran. ⁴Cardiovascular Health Research Unit, Department of Medicine, University of Washington, Seattle, WA, USA. ⁵National Center for PTSD, Behavioral Sciences Division, VA Boston Healthcare System, Boston, MA, USA.

⁶Department of Psychiatry and Biomedical Genetics, Boston University School of Medicine, Boston, MA, USA. ⁷University of Texas Health Sciences Center at Houston, Houston, TX, USA. ⁸University of Bordeaux, Inserm, Bordeaux Population Health Research Center, team VINTAGE, UMR 1219, F-33000 Bordeaux, France. ⁹Division of Systems Biomedicine and Pharmacology, Leiden Academic Centre for Drug Research, Leiden University, Leiden, The Netherlands. ¹⁰Department of Biostatistics, Boston University School of Public Health, Boston, MA, USA. ¹¹Boston University and the NHLBI's Framingham Heart Study, Boston, MA, USA. ¹²Department of Neurology, Boston University School of Medicine, Boston, MA, USA. ¹³Centre Inserm U1219 Bordeaux Population Health, CIC1401-EC, Institut de Santé Publique, d'Epidémiologie et de Développement, Université de Bordeaux, CHU de Bordeaux, Pôle Santé Publique, Bordeaux, France. ¹⁴Univ. Lille, Inserm, CHU Lille, Institut Pasteur de Lille, U1167 - RID-AGE- LabEx DISTALZ - Risk factors and molecular determinants of aging diseases, F-59000 Lille, France.

¹⁵Department of Neurology and Center for Neuroscience, University of California at Davis, Davis, CA, USA. ¹⁶Department of Medicine, Division of Aging, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA. ¹⁷Institute of Public Health Genetics and Department of Epidemiology and Cardiovascular Health Research Unit, University of Washington, Seattle, WA, USA. ¹⁸Department of Psychiatry and Center for Behavior Genetics of Aging, University of California, San Diego, La Jolla, CA, USA. ¹⁹Glenn Biggs Institute for Alzheimer's & Neurodegenerative Diseases, University of Texas Health Sciences Center, San Antonio, TX, USA. ²⁰Department of Psychology and Brain Sciences, Boston University, Boston, MA, USA. ²¹University of Mississippi Medical Center, Jackson, MS, USA. ²²Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, MA, USA. ²³Turner Institute for Brain and Mental Health, Monash University, Melbourne, VIC, Australia. ²⁴Harvard T.H. Chan School of Public Health, Harvard University, Boston, MA, USA. ²⁵Department of Neurosciences, University of California, San Diego, La Jolla, CA, USA. ²⁶Preventive Medicine & Epidemiology, Boston University School of Medicine, Boston, MA, USA. ²⁷Bordeaux University Hospital, Department of Neurology, Institute for Neurodegenerative Diseases, Bordeaux, France. ²⁸National Institute on Aging, Baltimore, MD, USA.

²⁹Departments of Neurology and Epidemiology, University of Washington, Seattle, WA, USA. ³⁰These authors contributed equally: Chloé Sarnowski, Mohsen Ghanbari, Joshua C. Bis, Mark Logue, Myriam Fornage, Aniket Mishra. ³¹These authors jointly supervised this work: Stephanie Debette, Carole Dufouil, Lenore J. Launer, William S. Kremen, William T. Longstreth, M Arfan Ikram, Sudha Seshadri. Data used in preparation of this article were obtained from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database (adni.loni.usc.edu). As such, the investigators within the ADNI contributed to the design and implementation of ADNI and/or provided data but did not participate in analysis or writing of this report. A complete listing of ADNI investigators can be found at: http://adni.loni.usc.edu/wp-content/uploads/how_to_apply/ADNI_Acknowledgement_List.pdf. ³²email: Chloe.Sarnowski@uth.tmc.edu