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/ ).

A proteogenomic view of Parkinson's disease causality and heterogeneity.

The pathogenesis and clinical heterogeneity of Parkinson's disease (PD) have been evaluated from molecular, pathophysiological, and clinical perspectives. High-throughput proteomic analysis of cerebrospinal fluid (CSF) opened new opportunities for scrutinizing this heterogeneity. To date, this is the most comprehensive CSF-based proteomics profiling study in PD with 569 patients (350 idiopathic patients, 65 GBA + mutation carriers and 154 LRRK2 + mutation carriers), 534 controls, and 4135 proteins analyzed. Combining CSF aptamer-based proteomics with genetics we determined protein quantitative trait loci (pQTLs). Analyses of pQTLs together with summary statistics from the largest PD genome wide association study (GWAS) identified 68 potential causal proteins by Mendelian randomization. The top causal protein, GPNMB, was previously reported to be upregulated in the substantia nigra of PD patients. We also compared the CSF proteomes of patients and controls. Proteome differences between GBA + patients and unaffected GBA + controls suggest degeneration of dopaminergic neurons, altered dopamine metabolism and increased brain inflammation. In the LRRK2 + subcohort we found dysregulated lysosomal degradation, altered alpha-synuclein processing, and neurotransmission. Proteome differences between idiopathic patients and controls suggest increased neuroinflammation, mitochondrial dysfunction/oxidative stress, altered iron metabolism and potential neuroprotection mediated by vasoactive substances. Finally, we used proteomic data to stratify idiopathic patients into "endotypes". The identified endotypes show differences in cognitive and motor disease progression based on previously reported protein-based risk scores.Our findings not only contribute to the identification of new therapeutic targets but also to shape personalized medicine in CNS neurodegeneration.

MS4A4A modifies the risk of Alzheimer disease by regulating lipid metabolism and immune response in a unique microglia state.

Genome-wide association studies (GWAS) have identified many modifiers of Alzheimer disease (AD) risk enriched in microglia. Two of these modifiers are common variants in the MS4A locus (rs1582763: protective and rs6591561: risk) and serve as major regulators of CSF sTREM2 levels. To understand their functional impact on AD, we used single nucleus transcriptomics to profile brains from carriers of these variants. We discovered a "chemokine" microglial subpopulation that is altered in MS4A variant carriers and for which MS4A4A is the major regulator. The protective variant increases MS4A4A expression and shifts the chemokine microglia subpopulation to an interferon state, while the risk variant suppresses MS4A4A expression and reduces this subpopulation of microglia. Our findings provide a mechanistic explanation for the AD variants in the MS4A locus. Further, they pave the way for future mechanistic studies of AD variants and potential therapeutic strategies for enhancing microglia resilience in AD pathogenesis.

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.

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.

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.

Examination of the Effect of Rare Variants in TREM2, ABI3, and PLCG2 in LOAD Through Multiple Phenotypes.

BACKGROUND: Rare variants in PLCG2 (p.P522R), ABI3 (p.S209F), and TREM2 (p.R47H, p.R62H) have been associated with late onset Alzheimer's disease (LOAD) risk in Caucasians. After the initial report, several studies have found positive results in cohorts of different ethnic background and with different phenotype. OBJECTIVE: In this study, we aim to evaluate the association of rare coding variants in PLCG2, ABI3, and TREM2 with LOAD risk and their effect at different time points of the disease. METHODS: We used a European American cohort to assess the association of the variants prior onset (using CSF Aß42, tau, and pTau levels, and amyloid imaging as endophenotypes) and after onset (measured as rate of memory decline). RESULTS: We confirm the association with LOAD risk of TREM2 p.R47H, p.R62H and ABI3 p.S209F variants, and the protective effect of PLCG2 p.P522R. In addition, ABI3 and TREM2 gene-sets showed significant association with LOAD risk. TREM2 p.R47H and PLCG2 p.P522R variants were also statistically associated with increase of amyloid imaging and AD progression, respectively. We did not observe any association of ABI3 p.S209F with any of the other AD endophenotypes. CONCLUSION: The results of this study highlight the importance of including biomarkers and alternative phenotypes to better understand the role of novel candidate genes with the disease.

Genome-wide association study of rate of cognitive decline in Alzheimer's disease patients identifies novel genes and pathways.

INTRODUCTION: Variability exists in the disease trajectories of Alzheimer's disease (AD) patients. We performed a genome-wide association study to examine rate of cognitive decline (ROD) in patients with AD. METHODS: We tested for interactions between genetic variants and time since diagnosis to predict the ROD of a composite cognitive score in 3946 AD cases and performed pathway analysis on the top genes. RESULTS: Suggestive associations (P < 1.0 × 10-6 ) were observed on chromosome 15 in DNA polymerase-γ (rs3176205, P = 1.11 × 10-7 ), chromosome 7 (rs60465337,P = 4.06 × 10-7 ) in contactin-associated protein-2, in RP11-384F7.1 on chromosome 3 (rs28853947, P = 5.93 × 10-7 ), family with sequence similarity 214 member-A on chromosome 15 (rs2899492, P = 5.94 × 10-7 ), and intergenic regions on chromosomes 16 (rs4949142, P = 4.02 × 10-7 ) and 4 (rs1304013, P = 7.73 × 10-7 ). Significant pathways involving neuronal development and function, apoptosis, memory, and inflammation were identified. DISCUSSION: Pathways related to AD, intelligence, and neurological function determine AD progression, while previously identified AD risk variants, including the apolipoprotein (APOE) ε4 and ε2 variants, do not have a major impact.

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.

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.

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.

Integrative system biology analyses of CRISPR-edited iPSC-derived neurons and human brains reveal deficiencies of presynaptic signaling in FTLD and PSP.

Mutations in the microtubule-associated protein tau (MAPT) gene cause autosomal dominant frontotemporal lobar degeneration with tau inclusions (FTLD-tau). MAPT p.R406W carriers present clinically with progressive memory loss and neuropathologically with neuronal and glial tauopathy. However, the pathogenic events triggered by the expression of the mutant tau protein remain poorly understood. To identify the genes and pathways that are dysregulated in FTLD-tau, we performed transcriptomic analyses in induced pluripotent stem cell (iPSC)-derived neurons carrying MAPT p.R406W and CRISPR/Cas9-corrected isogenic controls. We found that the expression of the MAPT p.R406W mutation was sufficient to create a significantly different transcriptomic profile compared with that of the isogeneic controls and to cause the differential expression of 328 genes. Sixty-one of these genes were also differentially expressed in the same direction between MAPT p.R406W carriers and pathology-free human control brains. We found that genes differentially expressed in the stem cell models and human brains were enriched for pathways involving gamma-aminobutyric acid (GABA) receptors and pre-synaptic function. The expression of GABA receptor genes, including GABRB2 and GABRG2, were consistently reduced in iPSC-derived neurons and brains from MAPT p.R406W carriers. Interestingly, we found that GABA receptor genes, including GABRB2 and GABRG2, are significantly lower in symptomatic mouse models of tauopathy, as well as in brains with progressive supranuclear palsy. Genome wide association analyses reveal that common variants within GABRB2 are associated with increased risk for frontotemporal dementia (P < 1 × 10-3). Thus, our systems biology approach, which leverages molecular data from stem cells, animal models, and human brain tissue can reveal novel disease mechanisms. Here, we demonstrate that MAPT p.R406W is sufficient to induce changes in GABA-mediated signaling and synaptic function, which may contribute to the pathogenesis of FTLD-tau and other primary tauopathies.

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.

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.

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.

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.

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.