Multi-ancestry GWAS reveals excitotoxicity associated with outcome after ischaemic stroke.

During the first hours after stroke onset, neurological deficits can be highly unstable: some patients rapidly improve, while others deteriorate. This early neurological instability has a major impact on long-term outcome. Here, we aimed to determine the genetic architecture of early neurological instability measured by the difference between the National Institutes of Health Stroke Scale (NIHSS) within 6 h of stroke onset and NIHSS at 24 h. A total of 5876 individuals from seven countries (Spain, Finland, Poland, USA, Costa Rica, Mexico and Korea) were studied using a multi-ancestry meta-analyses. We found that 8.7% of NIHSS at 24 h of variance was explained by common genetic variations, and also that early neurological instability has a different genetic architecture from that of stroke risk. Eight loci (1p21.1, 1q42.2, 2p25.1, 2q31.2, 2q33.3, 5q33.2, 7p21.2 and 13q31.1) were genome-wide significant and explained 1.8% of the variability suggesting that additional variants influence early change in neurological deficits. We used functional genomics and bioinformatic annotation to identify the genes driving the association from each locus. Expression quantitative trait loci mapping and summary data-based Mendelian randomization indicate that ADAM23 (log Bayes factor = 5.41) was driving the association for 2q33.3. Gene-based analyses suggested that GRIA1 (log Bayes factor = 5.19), which is predominantly expressed in the brain, is the gene driving the association for the 5q33.2 locus. These analyses also nominated GNPAT (log Bayes factor = 7.64) ABCB5 (log Bayes factor = 5.97) for the 1p21.1 and 7p21.1 loci. Human brain single-nuclei RNA-sequencing indicates that the gene expression of ADAM23 and GRIA1 is enriched in neurons. ADAM23, a presynaptic protein and GRIA1, a protein subunit of the AMPA receptor, are part of a synaptic protein complex that modulates neuronal excitability. These data provide the first genetic evidence in humans that excitotoxicity may contribute to early neurological instability after acute ischaemic stroke.

ApoE4 markedly exacerbates tau-mediated neurodegeneration in a mouse model of tauopathy.

APOE4 is the strongest genetic risk factor for late-onset Alzheimer disease. ApoE4 increases brain amyloid-ß pathology relative to other ApoE isoforms. However, whether APOE independently influences tau pathology, the other major proteinopathy of Alzheimer disease and other tauopathies, or tau-mediated neurodegeneration, is not clear. By generating P301S tau transgenic mice on either a human ApoE knock-in (KI) or ApoE knockout (KO) background, here we show that P301S/E4 mice have significantly higher tau levels in the brain and a greater extent of somatodendritic tau redistribution by three months of age compared with P301S/E2, P301S/E3, and P301S/EKO mice. By nine months of age, P301S mice with different ApoE genotypes display distinct phosphorylated tau protein (p-tau) staining patterns. P301S/E4 mice develop markedly more brain atrophy and neuroinflammation than P301S/E2 and P301S/E3 mice, whereas P301S/EKO mice are largely protected from these changes. In vitro, E4-expressing microglia exhibit higher innate immune reactivity after lipopolysaccharide treatment. Co-culturing P301S tau-expressing neurons with E4-expressing mixed glia results in a significantly higher level of tumour-necrosis factor-α (TNF-α) secretion and markedly reduced neuronal viability compared with neuron/E2 and neuron/E3 co-cultures. Neurons co-cultured with EKO glia showed the greatest viability with the lowest level of secreted TNF-α. Treatment of P301S neurons with recombinant ApoE (E2, E3, E4) also leads to some neuronal damage and death compared with the absence of ApoE, with ApoE4 exacerbating the effect. In individuals with a sporadic primary tauopathy, the presence of an ε4 allele is associated with more severe regional neurodegeneration. In individuals who are positive for amyloid-ß pathology with symptomatic Alzheimer disease who usually have tau pathology, ε4-carriers demonstrate greater rates of disease progression. Our results demonstrate that ApoE affects tau pathogenesis, neuroinflammation, and tau-mediated neurodegeneration independently of amyloid-ß pathology. ApoE4 exerts a 'toxic' gain of function whereas the absence of ApoE is protective.

Genome-wide association study of brain amyloid deposition as measured by Pittsburgh Compound-B (PiB)-PET imaging.

Deposition of amyloid plaques in the brain is one of the two main pathological hallmarks of Alzheimer's disease (AD). Amyloid positron emission tomography (PET) is a neuroimaging tool that selectively detects in vivo amyloid deposition in the brain and is a reliable endophenotype for AD that complements cerebrospinal fluid biomarkers with regional information. We measured in vivo amyloid deposition in the brains of ~1000 subjects from three collaborative AD centers and ADNI using 11C-labeled Pittsburgh Compound-B (PiB)-PET imaging followed by meta-analysis of genome-wide association studies, first to our knowledge for PiB-PET, to identify novel genetic loci for this endophenotype. The APOE region showed the most significant association where several SNPs surpassed the genome-wide significant threshold, with APOE*4 being most significant (P-meta = 9.09E-30; ß = 0.18). Interestingly, after conditioning on APOE*4, 14 SNPs remained significant at P < 0.05 in the APOE region that were not in linkage disequilibrium with APOE*4. Outside the APOE region, the meta-analysis revealed 15 non-APOE loci with P < 1E-05 on nine chromosomes, with two most significant SNPs on chromosomes 8 (P-meta = 4.87E-07) and 3 (P-meta = 9.69E-07). Functional analyses of these SNPs indicate their potential relevance with AD pathogenesis. Top 15 non-APOE SNPs along with APOE*4 explained 25-35% of the amyloid variance in different datasets, of which 14-17% was explained by APOE*4 alone. In conclusion, we have identified novel signals in APOE and non-APOE regions that affect amyloid deposition in the brain. Our data also highlights the presence of yet to be discovered variants that may be responsible for the unexplained genetic variance of amyloid deposition.

The TMEM106B FTLD-protective variant, rs1990621, is also associated with increased neuronal proportion.

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

Sex-specific genetic predictors of Alzheimer's disease biomarkers.

Cerebrospinal fluid (CSF) levels of amyloid-ß 42 (Aß42) and tau have been evaluated as endophenotypes in Alzheimer's disease (AD) genetic studies. Although there are sex differences in AD risk, sex differences have not been evaluated in genetic studies of AD endophenotypes. We performed sex-stratified and sex interaction genetic analyses of CSF biomarkers to identify sex-specific associations. Data came from a previous genome-wide association study (GWAS) of CSF Aß42 and tau (1527 males, 1509 females). We evaluated sex interactions at previous loci, performed sex-stratified GWAS to identify sex-specific associations, and evaluated sex interactions at sex-specific GWAS loci. We then evaluated sex-specific associations between prefrontal cortex (PFC) gene expression at relevant loci and autopsy measures of plaques and tangles using data from the Religious Orders Study and Rush Memory and Aging Project. In Aß42, we observed sex interactions at one previous and one novel locus: rs316341 within SERPINB1 (p = 0.04) and rs13115400 near LINC00290 (p = 0.002). These loci showed stronger associations among females (ß = - 0.03, p = 4.25 × 10-8; ß = 0.03, p = 3.97 × 10-8) than males (ß = - 0.02, p = 0.009; ß = 0.01, p = 0.20). Higher levels of expression of SERPINB1, SERPINB6, and SERPINB9 in PFC was associated with higher levels of amyloidosis among females (corrected p values < 0.02) but not males (p > 0.38). In total tau, we observed a sex interaction at a previous locus, rs1393060 proximal to GMNC (p = 0.004), driven by a stronger association among females (ß = 0.05, p = 4.57 × 10-10) compared to males (ß = 0.02, p = 0.03). There was also a sex-specific association between rs1393060 and tangle density at autopsy (pfemale = 0.047; pmale = 0.96), and higher levels of expression of two genes within this locus were associated with lower tangle density among females (OSTN p = 0.006; CLDN16 p = 0.002) but not males (p ≥ 0.32). Results suggest a female-specific role for SERPINB1 in amyloidosis and for OSTN and CLDN16 in tau pathology. Sex-specific genetic analyses may improve understanding of AD's genetic architecture.

Polygenic risk score of sporadic late-onset Alzheimer's disease reveals a shared architecture with the familial and early-onset forms.

OBJECTIVE: To determine whether the extent of overlap of the genetic architecture among the sporadic late-onset Alzheimer's Disease (sLOAD), familial late-onset AD (fLOAD), sporadic early-onset AD (sEOAD), and autosomal dominant early-onset AD (eADAD). METHODS: Polygenic risk scores (PRSs) were constructed using previously identified 21 genome-wide significant loci for LOAD risk. RESULTS: We found that there is an overlap in the genetic architecture among sEOAD, fLOAD, and sLOAD. The highest association of the PRS and risk (odds ratio [OR] = 2.27; P = 1.29 × 10-7) was observed in sEOAD, followed by fLOAD (OR = 1.75; P = 1.12 × 10-7) and sLOAD (OR = 1.40; P = 1.21 × 10-3). The PRS was associated with cerebrospinal fluid ptau181-Aß42 on eADAD (P = 4.36 × 10-2). CONCLUSION: Our analysis confirms that the genetic factors identified for LOAD modulate risk in sLOAD and fLOAD and also sEOAD cohorts. Specifically, our results suggest that the burden of these risk variants is associated with familial clustering and earlier onset of AD. Although these variants are not associated with risk in the eADAD, they may be modulating age at onset.

Genome-wide association study identifies four novel loci associated with Alzheimer's endophenotypes and disease modifiers.

More than 20 genetic loci have been associated with risk for Alzheimer's disease (AD), but reported genome-wide significant loci do not account for all the estimated heritability and provide little information about underlying biological mechanisms. Genetic studies using intermediate quantitative traits such as biomarkers, or endophenotypes, benefit from increased statistical power to identify variants that may not pass the stringent multiple test correction in case-control studies. Endophenotypes also contain additional information helpful for identifying variants and genes associated with other aspects of disease, such as rate of progression or onset, and provide context to interpret the results from genome-wide association studies (GWAS). We conducted GWAS of amyloid beta (Aß42), tau, and phosphorylated tau (ptau181) levels in cerebrospinal fluid (CSF) from 3146 participants across nine studies to identify novel variants associated with AD. Five genome-wide significant loci (two novel) were associated with ptau181, including loci that have also been associated with AD risk or brain-related phenotypes. Two novel loci associated with Aß42 near GLIS1 on 1p32.3 (ß = -0.059, P = 2.08 × 10-8) and within SERPINB1 on 6p25 (ß = -0.025, P = 1.72 × 10-8) were also associated with AD risk (GLIS1: OR = 1.105, P = 3.43 × 10-2), disease progression (GLIS1: ß = 0.277, P = 1.92 × 10-2), and age at onset (SERPINB1: ß = 0.043, P = 4.62 × 10-3). Bioinformatics indicate that the intronic SERPINB1 variant (rs316341) affects expression of SERPINB1 in various tissues, including the hippocampus, suggesting that SERPINB1 influences AD through an Aß-associated mechanism. Analyses of known AD risk loci suggest CLU and FERMT2 may influence CSF Aß42 (P = 0.001 and P = 0.009, respectively) and the INPP5D locus may affect ptau181 levels (P = 0.009); larger studies are necessary to verify these results. Together the findings from this study can be used to inform future AD studies.

Parkinson disease polygenic risk score is associated with Parkinson disease status and age at onset but not with alpha-synuclein cerebrospinal fluid levels.

BACKGROUND: The genetic architecture of Parkinson's Disease (PD) is complex and not completely understood. Multiple genetic studies to date have identified multiple causal genes and risk loci. Nevertheless, most of the expected genetic heritability remains unexplained. Polygenic risk scores (PRS) may provide greater statistical power and inform about the genetic architecture of multiple phenotypes. The aim of this study was to test the association between PRS and PD risk, age at onset and cerebrospinal fluid (CSF) biomarkers (α-synuclein, Aß1-42, t-tau and p-tau). METHODS: The weighted PRS was created using the genome-wide loci from Nalls et al., 2014 PD GWAs meta-analysis. The PRS was tested for association with PD status, age at onset and CSF biomarker levels in 829 cases and 432 controls of European ancestry. RESULTS: The PRS was associated with PD status (p = 5.83×10-08) and age at onset (p = 5.70×10-07). The CSF t-tau levels showed a nominal association with the PRS (p = 0.02). However, CSF α-synuclein, amyloid beta and phosphorylated tau were not found to be associated with the PRS. CONCLUSION: Our study suggests that there is an overlap in the genetic architecture of PD risk and onset, although the different loci present different weights for those phenotypes. In our dataset we found a marginal association of the PRS with CSF t-tau but not with α-synuclein CSF levels, suggesting that the genetic architecture for the CSF biomarker levels is different from that of PD risk.

Cerebrospinal fluid soluble TREM2 is higher in Alzheimer disease and associated with mutation status.

Low frequency coding variants in TREM2 are associated with increased Alzheimer disease (AD) risk, while loss of functions mutations in the gene lead to an autosomal recessive early-onset dementia, named Nasu-Hakola disease (NHD). TREM2 can be detected as a soluble protein in cerebrospinal fluid (CSF) and plasma, and its CSF levels are elevated in inflammatory CNS diseases. We measured soluble TREM2 (sTREM2) in the CSF of a large AD case-control dataset (n = 180) and 40 TREM2 risk variant carriers to determine whether CSF sTREM2 levels are associated with AD status or mutation status. We also performed genetic studies to identify genetic variants associated with CSF sTREM2 levels. CSF, but not plasma, sTREM2 was highly correlated with CSF total tau and phosphorylated-tau levels (r = 0.35, P < 1×10(-4); r = 0.40, P < 1×10(-4), respectively), but not with CSF Aß42. AD cases presented higher CSF sTREM2 levels than controls (P = 0.01). Carriers of NHD-associated TREM2 variants presented significantly lower CSF sTREM2 levels, supporting the hypothesis that these mutations lead to reduced protein production/function (R136Q, D87N, Q33X or T66M; P = 1×10(-3)). In contrast, CSF sTREM2 levels were significantly higher in R47H carriers compared to non-carriers (P = 6×10(-3)), suggesting that this variant does not impact protein expression and increases AD risk through a different pathogenic mechanism than NHD variants. In GWAS analyses for CSF sTREM2 levels the most significant signal was located on the MS4A gene locus (P = 5.45 × 10(-07)) corresponding to one of the SNPs reported to be associated with AD risk in this locus. Furthermore, SNPs involved in pathways related to virus cellular entry and vesicular trafficking were overrepresented, suggesting that CSF sTREM2 levels could be an informative phenotype for AD.

Chitinase-3-like 1 protein (CHI3L1) locus influences cerebrospinal fluid levels of YKL-40.

BACKGROUND: Alzheimer's disease (AD) pathology appears several years before clinical symptoms, so identifying ways to detect individuals in the preclinical stage is imperative. The cerebrospinal fluid (CSF) Tau/Aß42 ratio is currently the best known predictor of AD status and cognitive decline, and the ratio of CSF levels of chitinase-3-like 1 protein (CHI3L1, YKL-40) and amyloid beta (Aß42) were reported as predictive, but individual variability and group overlap inhibits their utility for individual diagnosis making it necessary to find ways to improve sensitivity of these biomarkers. METHODS: We used linear regression to identify genetic loci associated with CSF YKL-40 levels in 379 individuals (80 cognitively impaired and 299 cognitively normal) from the Charles F and Joanne Knight Alzheimer's Disease Research Center. We tested correlations between YKL-40 and CSF Tau/Aß42 ratio, Aß42, tau, and phosphorylated tau (ptau181). We used studentized residuals from a linear regression model of the log-transformed, standardized protein levels and the additive reference allele counts from the most significant locus to adjust YKL-40 values and tested the differences in correlations with CSF Tau/Aß42 ratio, Aß42, tau, and ptau181. RESULTS: We found that genetic variants on the CH13L1 locus were significantly associated with CSF YKL-40 levels, but not AD risk, age at onset, or disease progression. The most significant variant is a reported expression quantitative trait locus for CHI3L1, the gene which encodes YKL-40, and explained 12.74 % of the variance in CSF YKL-40 in our study. YKL-40 was positively correlated with ptau181 (r = 0.521) and the strength of the correlation significantly increased with the addition of genetic information (r = 0.573, p = 0.006). CONCLUSIONS: CSF YKL-40 levels are likely a biomarker for AD, but we found no evidence that they are an AD endophenotype. YKL-40 levels are highly regulated by genetic variation, and by including genetic information the strength of the correlation between YKL-40 and ptau181 levels is significantly improved. Our results suggest that studies of potential biomarkers may benefit from including genetic information.

Single-nucleus RNA-sequencing of autosomal dominant Alzheimer disease and risk variant carriers.

Genetic studies of Alzheimer disease (AD) have prioritized variants in genes related to the amyloid cascade, lipid metabolism, and neuroimmune modulation. However, the cell-specific effect of variants in these genes is not fully understood. Here, we perform single-nucleus RNA-sequencing (snRNA-seq) on nearly 300,000 nuclei from the parietal cortex of AD autosomal dominant (APP and PSEN1) and risk-modifying variant (APOE, TREM2 and MS4A) carriers. Within individual cell types, we capture genes commonly dysregulated across variant groups. However, specific transcriptional states are more prevalent within variant carriers. TREM2 oligodendrocytes show a dysregulated autophagy-lysosomal pathway, MS4A microglia have dysregulated complement cascade genes, and APOEε4 inhibitory neurons display signs of ferroptosis. All cell types have enriched states in autosomal dominant carriers. We leverage differential expression and single-nucleus ATAC-seq to map GWAS signals to effector cell types including the NCK2 signal to neurons in addition to the initially proposed microglia. Overall, our results provide insights into the transcriptional diversity resulting from AD genetic architecture and cellular heterogeneity. The data can be explored on the online browser ( http://web.hararilab.org/SNARE/ ).

Chi3l1/YKL-40 is controlled by the astrocyte circadian clock and regulates neuroinflammation and Alzheimer's disease pathogenesis.

Regulation of glial activation and neuroinflammation are critical factors in the pathogenesis of Alzheimer's disease (AD). YKL-40, a primarily astrocytic protein encoded by the gene Chi3l1, is a widely studied cerebrospinal fluid biomarker that increases with aging and early in AD. However, the function of Chi3l1/YKL-40 in AD is unknown. In a cohort of patients with AD, we observed that a variant in the human CHI3L1 gene, which results in decreased CSF YKL-40 expression, was associated with slower AD progression. At baseline, Chi3l1 deletion in mice had no effect on astrocyte activation while modestly promoting microglial activation. In a mouse APP/PS1 model of AD, Chi3l1 deletion decreased amyloid plaque burden and increased periplaque expression of the microglial lysosomal marker CD68, suggesting that Chi3l1 may suppress glial phagocytic activation and promote amyloid accumulation. Accordingly, Chi3l1 knockdown increased phagocytosis of zymosan particles and of ß-amyloid peptide in both astrocytes and microglia in vitro. We further observed that expression of Chi3l1 is regulated by the circadian clock, as deletion of the core clock proteins BMAL1 or CLOCK/NPAS2 strongly suppresses basal Chi3l1 expression, whereas deletion of the negative clock regulators PER1/PER2 increased Chi3l1 expression. Basal Chi3l1 mRNA was nonrhythmic because of a long mRNA half-life in astrocytes. However, inflammatory induction of Chi3l1 was gated by the clock. Our findings reveal Chi3l1/YKL-40 as a modulator of glial phagocytic activation and AD pathogenesis in both mice and humans and suggest that the astrocyte circadian clock regulates inflammatory Chi3l1 induction.

TREM2 brain transcript-specific studies in AD and TREM2 mutation carriers.

BACKGROUND: Low frequency coding variants in TREM2 are associated with Alzheimer disease (AD) risk and cerebrospinal fluid (CSF) TREM2 protein levels are different between AD cases and controls. Similarly, TREM2 risk variant carriers also exhibit differential CSF TREM2 levels. TREM2 has three different alternative transcripts, but most of the functional studies only model the longest transcript. No studies have analyzed TREM2 expression levels or alternative splicing in brains from AD and cognitively normal individuals. We wanted to determine whether there was differential expression of TREM2 in sporadic-AD cases versus AD-TREM2 carriers vs sex- and aged-matched normal controls; and if this differential expression was due to a particular TREM2 transcript. METHODS: We analyzed RNA-Seq data from parietal lobe brain tissue from AD cases with TREM2 variants (n = 33), AD cases (n = 195) and healthy controls (n = 118), from three independent datasets using Kallisto and the R package tximport to determine the read count for each transcript and quantified transcript abundance as transcripts per million. RESULTS: The three TREM2 transcripts were expressed in brain cortex in the three datasets. We demonstrate for the first time that the transcript that lacks the transmembrane domain and encodes a soluble form of TREM2 (sTREM2) has an expression level around 60% of the canonical transcript, suggesting that around 25% of the sTREM2 protein levels could be explained by this transcript. We did not observe a difference in the overall TREM2 expression level between cases and controls. However, the isoform which lacks the 5' exon, but includes the transmembrane domain, was significantly lower in TREM2- p.R62H carriers than in AD cases (p = 0.007). CONCLUSION: Using bulk RNA-Seq data from three different cohorts, we were able to quantify the expression level of the three TREM2 transcripts, demonstrating: (1) all three transcripts of them are highly expressed in the human cortex, (2) that up to 25% of the sTREM2 may be due to the expression of a specific isoform and not TREM2 cleavage; and (3) that TREM2 risk variants do not affect expression levels, suggesting that the effect of the TREM2 variants on CSF levels occurs at post-transcriptional level.

An atlas of cortical circular RNA expression in Alzheimer disease brains demonstrates clinical and pathological associations.

Parietal cortex RNA-sequencing (RNA-seq) data were generated from individuals with and without Alzheimer disease (AD; ncontrol = 13; nAD = 83) from the Knight Alzheimer Disease Research Center (Knight ADRC). Using this and an independent (Mount Sinai Brain Bank (MSBB)) AD RNA-seq dataset, cortical circular RNA (circRNA) expression was quantified in the context of AD. Significant associations were identified between circRNA expression and AD diagnosis, clinical dementia severity and neuropathological severity. It was demonstrated that most circRNA-AD associations are independent of changes in cognate linear messenger RNA expression or estimated brain cell-type proportions. Evidence was provided for circRNA expression changes occurring early in presymptomatic AD and in autosomal dominant AD. It was also observed that AD-associated circRNAs co-expressed with known AD genes. Finally, potential microRNA-binding sites were identified in AD-associated circRNAs for miRNAs predicted to target AD genes. Together, these results highlight the importance of analyzing non-linear RNAs and support future studies exploring the potential roles of circRNAs in AD pathogenesis.

A single-nuclei RNA sequencing study of Mendelian and sporadic AD in the human brain.

BACKGROUND: Alzheimer's disease (AD) is the most common form of dementia. This neurodegenerative disorder is associated with neuronal death and gliosis heavily impacting the cerebral cortex. AD has a substantial but heterogeneous genetic component, presenting both Mendelian and complex genetic architectures. Using bulk RNA-seq from the parietal lobes and deconvolution methods, we previously reported that brains exhibiting different AD genetic architecture exhibit different cellular proportions. Here, we sought to directly investigate AD brain changes in cell proportion and gene expression using single-cell resolution. METHODS: We generated unsorted single-nuclei RNA sequencing data from brain tissue. We leveraged the tissue donated from a carrier of a Mendelian genetic mutation, PSEN1 p.A79V, and two family members who suffer from sporadic AD, but do not carry any autosomal mutations. We evaluated alternative alignment approaches to maximize the titer of reads, genes, and cells with high quality. In addition, we employed distinct clustering strategies to determine the best approach to identify cell clusters that reveal neuronal and glial cell types and avoid artifacts such as sample and batch effects. We propose an approach to cluster cells that reduces biases and enable further analyses. RESULTS: We identified distinct types of neurons, both excitatory and inhibitory, and glial cells, including astrocytes, oligodendrocytes, and microglia, among others. In particular, we identified a reduced proportion of excitatory neurons in the Mendelian mutation carrier, but a similar distribution of inhibitory neurons. Furthermore, we investigated whether single-nuclei RNA-seq from the human brains recapitulate the expression profile of disease-associated microglia (DAM) discovered in mouse models. We also determined that when analyzing human single-nuclei data, it is critical to control for biases introduced by donor-specific expression profiles. CONCLUSION: We propose a collection of best practices to generate a highly detailed molecular cell atlas of highly informative frozen tissue stored in brain banks. Importantly, we have developed a new web application to make this unique single-nuclei molecular atlas publicly available.

The MS4A gene cluster is a key modulator of soluble TREM2 and Alzheimer's disease risk.

Soluble triggering receptor expressed on myeloid cells 2 (sTREM2) in cerebrospinal fluid (CSF) has been associated with Alzheimer's disease (AD). TREM2 plays a critical role in microglial activation, survival, and phagocytosis; however, the pathophysiological role of sTREM2 in AD is not well understood. Understanding the role of sTREM2 in AD may reveal new pathological mechanisms and lead to the identification of therapeutic targets. We performed a genome-wide association study (GWAS) to identify genetic modifiers of CSF sTREM2 obtained from the Alzheimer's Disease Neuroimaging Initiative. Common variants in the membrane-spanning 4-domains subfamily A (MS4A) gene region were associated with CSF sTREM2 concentrations (rs1582763; P = 1.15 × 10-15); this was replicated in independent datasets. The variants associated with increased CSF sTREM2 concentrations were associated with reduced AD risk and delayed age at onset of disease. The single-nucleotide polymorphism rs1582763 modified expression of the MS4A4A and MS4A6A genes in multiple tissues, suggesting that one or both of these genes are important for modulating sTREM2 production. Using human macrophages as a proxy for microglia, we found that MS4A4A and TREM2 colocalized on lipid rafts at the plasma membrane, that sTREM2 increased with MS4A4A overexpression, and that silencing of MS4A4A reduced sTREM2 production. These genetic, molecular, and cellular findings suggest that MS4A4A modulates sTREM2. These findings also provide a mechanistic explanation for the original GWAS signal in the MS4A locus for AD risk and indicate that TREM2 may be involved in AD pathogenesis not only in TREM2 risk-variant carriers but also in those with sporadic disease.

Genome-wide association study for variants that modulate relationships between cerebrospinal fluid amyloid-beta 42, tau, and p-tau levels.

BACKGROUND: A relationship quantitative trait locus exists when the correlation between multiple traits varies by genotype for that locus. Relationship quantitative trait loci (rQTL) are often involved in gene-by-gene (G×G) interactions or gene-by-environmental interactions, making them a powerful tool for detecting G×G. METHODS: We performed genome-wide association studies to identify rQTL between tau and Aß42 and ptau and Aß42 with over 3000 individuals using age, gender, series, APOE ε2, APOE ε4, and two principal components for population structure as covariates. Each significant rQTL was separately screened for interactions with other loci for each trait in the rQTL model. Parametric bootstrapping was used to assess significance. RESULTS: We found four significant tau/Aß42 rQTL from three unique locations and six ptau/Aß42 rQTL from five unique locations. G×G screens with these rQTL produced four significant G×G interactions (one Aß42, two ptau, and one tau) with four rQTL where each second locus was from a unique location. On follow-up, rs1036819 and rs74025622 were associated with Alzheimer's disease (AD) case/control status; rs15205 and rs79099429 were associated with rate of decline. CONCLUSIONS: The two most significant rQTL (rs8027714 and rs1036819) for ptau/Aß42 are on different chromosomes and both are strong hits for pelvic organ prolapse. While diseases of the nervous system can cause pelvic organ prolapse, it is unlikely related to the ptau/Aß42 relationship but may suggest that these two loci share a pathway. In addition to a ptau/Aß42 rQTL and association with AD case/control status, rs1036819 is a strong rQTL for case/control status/Aß42 and for tau/Aß42. It resides in the ZFAT gene, which is related to autoimmune thyroid disease. For tau, rs9817620 interacts with the tau/Aß42 rQTL rs74025622. It is in the CHL1 gene, which is a neural cell adhesion molecule and may be involved in signal transduction pathways. CHL1 is related to BACE1, which is a ß-secretase enzyme that initiates production of the ß-amyloid peptide involved in AD and is a primary drug target. Overall, there are numerous loci that affect the relationship between these important AD endophenotypes and some are due to interactions with other loci. Some affect the risk of AD and/or rate of progression.

Evaluation of Gene-Based Family-Based Methods to Detect Novel Genes Associated With Familial Late Onset Alzheimer Disease.

Gene-based tests to study the combined effect of rare variants on a particular phenotype have been widely developed for case-control studies, but their evolution and adaptation for family-based studies, especially studies of complex incomplete families, has been slower. In this study, we have performed a practical examination of all the latest gene-based methods available for family-based study designs using both simulated and real datasets. We examined the performance of several collapsing, variance-component, and transmission disequilibrium tests across eight different software packages and 22 models utilizing a cohort of 285 families (N = 1,235) with late-onset Alzheimer disease (LOAD). After a thorough examination of each of these tests, we propose a methodological approach to identify, with high confidence, genes associated with the tested phenotype and we provide recommendations to select the best software and model for family-based gene-based analyses. Additionally, in our dataset, we identified PTK2B, a GWAS candidate gene for sporadic AD, along with six novel genes (CHRD, CLCN2, HDLBP, CPAMD8, NLRP9, and MAS1L) as candidate genes for familial LOAD.

Assessment of the Genetic Architecture of Alzheimer's Disease Risk in Rate of Memory Decline.

Many genetic studies for Alzheimer's disease (AD) have been focused on the identification of common genetic variants associated with AD risk and not on other aspects of the disease, such as age at onset or rate of dementia progression. There are multiple approaches to untangling the genetic architecture of these phenotypes. We hypothesized that the genetic architecture of rate of progression is different than the risk for developing AD dementia. To test this hypothesis, we used longitudinal clinical data from ADNI and the Knight-ADRC at Washington University, and we calculated PRS (polygenic risk score) based on the IGAP study to compare the genetic architecture of AD risk and dementia progression. Dementia progression was measured by the change of Clinical Dementia Rating Sum of Boxes (CDR)-SB per year. Out of the 21 loci for AD risk, no association with the rate of dementia progression was found. The PRS rate was significantly associated with the rate of dementia progression (ß= 0.146, p = 0.03). In the case of rare variants, TREM2 (ß= 0.309, p = 0.02) was also associated with the rate of dementia progression. TREM2 variant carriers showed a 23% faster rate of dementia compared with non-variant carriers. In conclusion, our results indicate that the recently identified common and rare variants for AD susceptibility have a limited impact on the rate of dementia progression in AD patients.

Genetic variants associated with Alzheimer's disease confer different cerebral cortex cell-type population structure.

BACKGROUND: Alzheimer's disease (AD) is characterized by neuronal loss and astrocytosis in the cerebral cortex. However, the specific effects that pathological mutations and coding variants associated with AD have on the cellular composition of the brain are often ignored. METHODS: We developed and optimized a cell-type-specific expression reference panel and employed digital deconvolution methods to determine brain cellular distribution in three independent transcriptomic studies. RESULTS: We found that neuronal and astrocyte relative proportions differ between healthy and diseased brains and also among AD cases that carry specific genetic risk variants. Brain carriers of pathogenic mutations in APP, PSEN1, or PSEN2 presented lower neuron and higher astrocyte relative proportions compared to sporadic AD. Similarly, the APOE ε4 allele also showed decreased neuronal and increased astrocyte relative proportions compared to AD non-carriers. In contrast, carriers of variants in TREM2 risk showed a lower degree of neuronal loss compared to matched AD cases in multiple independent studies. CONCLUSIONS: These findings suggest that genetic risk factors associated with AD etiology have a specific imprinting in the cellular composition of AD brains. Our digital deconvolution reference panel provides an enhanced understanding of the fundamental molecular mechanisms underlying neurodegeneration, enabling the analysis of large bulk RNA-sequencing studies for cell composition and suggests that correcting for the cellular structure when performing transcriptomic analysis will lead to novel insights of AD.