Search

Abstract Determining the genetic architecture of Alzheimer’s disease pathologies can enhance mechanistic understanding and inform precision medicine strategies. Here, we perform a genome-wide association study of cortical tau quantified by positron emission tomography in 3046 participants from 12 independent studies. The CYP1B1-RMDN2 locus is associated with tau deposition. The most significant signal is at rs2113389, explaining 4.3% of the variation in cortical tau, while APOE4 rs429358 accounts for 3.6%. rs2113389 is associated with higher tau and faster cognitive decline. Additive effects, but no interactions, are observed between rs2113389 and diagnosis, APOE4, and amyloid beta positivity. CYP1B1 expression is upregulated in AD. rs2113389 is associated with higher CYP1B1 expression and methylation levels. Mouse model studies provide additional functional evidence for a relationship between CYP1B1 and tau deposition but not amyloid beta. These results provide insight into the genetic basis of cerebral tau deposition and support novel pathways for therapeutic development in AD.

Abstract Progressive supranuclear palsy (PSP) is a neurodegenerative movement and cognitive disorder characterized by abnormal accumulation of the microtubule-associated protein tau in the brain. Biochemically, inclusions in PSP are enriched for tau proteoforms with four microtubule-binding domain repeats (4R), an isoform that arises from alternative tau pre-mRNA splicing. While preferential aggregation and reduced degradation of 4R tau protein is thought to play a role in inclusion formation and toxicity, an alternative hypothesis is that altered expression of tau mRNA isoforms plays a causal role. This stems from the observation that PSP is associated with common variation in the tau gene (MAPT) at the 17q21.31 locus which contains low copy number repeats flanking a large recurrent genomic inversion. The complex genomic structural changes at the locus give rise to two dominant haplotypes, termed H1 and H2, that have the potential to markedly influence gene expression. Here, we explored haplotype-dependent differences in gene expression using a bulk RNA-seq dataset derived from human post-mortem brain tissue from PSP (n = 84) and controls (n = 77) using a rigorous computational pipeline, including alternative pre-mRNA splicing. We found 3579 differentially expressed genes in the temporal cortex and 10,011 in the cerebellum. We also found 7214 differential splicing events in the temporal cortex and 18,802 in the cerebellum. In the cerebellum, total tau mRNA levels and the proportion of transcripts encoding 4R tau were significantly increased in PSP compared to controls. In the temporal cortex, the proportion of reads that expressed 4R tau was increased in cases compared to controls. 4R tau mRNA levels were significantly associated with the H1 haplotype in the temporal cortex. Further, we observed a marked haplotype-dependent difference in KANSL1 expression that was strongly associated with H1 in both brain regions. These findings support the hypothesis that sporadic PSP is associated with haplotype-dependent increases in 4R tau mRNA that might play a causal role in this disorder.

Abstract To uncover molecular changes underlying blood-brain-barrier dysfunction in Alzheimer’s disease, we performed single nucleus RNA sequencing in 24 Alzheimer’s disease and control brains and focused on vascular and astrocyte clusters as main cell types of blood-brain-barrier gliovascular-unit. The majority of the vascular transcriptional changes were in pericytes. Of the vascular molecular targets predicted to interact with astrocytic ligands, SMAD3, upregulated in Alzheimer’s disease pericytes, has the highest number of ligands including VEGFA, downregulated in Alzheimer’s disease astrocytes. We validated these findings with external datasets comprising 4,730 pericyte and 150,664 astrocyte nuclei. Blood SMAD3 levels are associated with Alzheimer’s disease-related neuroimaging outcomes. We determined inverse relationships between pericytic SMAD3 and astrocytic VEGFA in human iPSC and zebrafish models. Here, we detect vast transcriptome changes in Alzheimer’s disease at the gliovascular-unit, prioritize perturbed pericytic SMAD3-astrocytic VEGFA interactions, and validate these in cross-species models to provide a molecular mechanism of blood-brain-barrier disintegrity in Alzheimer’s disease.

Abstract Progressive supranuclear palsy (PSP) is a neurodegenerative parkinsonian disorder characterized by cell-type-specific tau lesions in neurons and glia. Prior work uncovered transcriptome changes in human PSP brains, although their cell-specificity is unknown. Further, systematic data integration and experimental validation platforms to prioritize brain transcriptional perturbations as therapeutic targets in PSP are currently lacking. In this study, we combine bulk tissue (n = 408) and single nucleus RNAseq (n = 34) data from PSP and control brains with transcriptome data from a mouse tauopathy and experimental validations in Drosophila tau models for systematic discovery of high-confidence expression changes in PSP with therapeutic potential. We discover, replicate, and annotate thousands of differentially expressed genes in PSP, many of which reside in glia-enriched co-expression modules and cells. We prioritize DDR2, STOM, and KANK2 as promising therapeutic targets in PSP with striking cross-species validations. We share our findings and data via our interactive application tool PSP RNAseq Atlas ( https://rtools.mayo.edu/PSP_RNAseq_Atlas/ ). Our findings reveal robust glial transcriptome changes in PSP, provide a cross-species systems biology approach, and a tool for therapeutic target discoveries in PSP with potential application in other neurodegenerative diseases.

Abstract Background The apolipoprotein E (APOE) gene is the strongest genetic risk factor for Alzheimer’s disease (AD); however, how it modulates brain homeostasis is not clear. The apoE protein is a major lipid carrier in the brain transporting lipids such as cholesterol among different brain cell types. Methods We generated three-dimensional (3-D) cerebral organoids from human parental iPSC lines and its isogenic APOE-deficient (APOE −/− ) iPSC line. To elucidate the cell-type-specific effects of APOE deficiency in the cerebral organoids, we performed scRNA-seq in the parental and APOE −/− cerebral organoids at Day 90. Results We show that APOE deficiency in human iPSC-derived cerebral organoids impacts brain lipid homeostasis by modulating multiple cellular and molecular pathways. Molecular profiling through single-cell RNA sequencing revealed that APOE deficiency leads to changes in cellular composition of isogenic cerebral organoids likely by modulating the eukaryotic initiation factor 2 (EIF2) signaling pathway as these events were alleviated by the treatment of an integrated stress response inhibitor (ISRIB). APOE deletion also leads to activation of the Wnt/β-catenin signaling pathway with concomitant decrease of secreted frizzled-related protein 1 (SFRP1) expression in glia cells. Importantly, the critical role of apoE in cell-type-specific lipid homeostasis was observed upon APOE deletion in cerebral organoids with a specific upregulation of cholesterol biosynthesis in excitatory neurons and excessive lipid accumulation in astrocytes. Relevant to human AD, APOE4 cerebral organoids show altered neurogenesis and cholesterol metabolism compared to those with APOE3. Conclusions Our work demonstrates critical roles of apoE in brain homeostasis and offers critical insights into the APOE4-related pathogenic mechanisms.

Abstract Neurogenesis, crucial for brain resilience, is reduced in Alzheimer’s disease (AD) that induces astroglial reactivity at the expense of the pro-neurogenic potential, and restoring neurogenesis could counteract neurodegenerative pathology. However, the molecular mechanisms promoting pro-neurogenic astroglial fate despite AD pathology are unknown. In this study, we used APP/PS1dE9 mouse model and induced Nerve growth factor receptor (Ngfr) expression in the hippocampus. Ngfr, which promotes neurogenic fate of astroglia during the amyloid pathology-induced neuroregeneration in zebrafish brain, stimulated proliferative and neurogenic outcomes. Histological analyses of the changes in proliferation and neurogenesis, single-cell transcriptomics, spatial proteomics, and functional knockdown studies showed that the induced expression of Ngfr reduced the reactive astrocyte marker Lipocalin-2 (Lcn2), which we found was sufficient to reduce neurogenesis in astroglia. Anti-neurogenic effects of Lcn2 was mediated by Slc22a17, blockage of which recapitulated the pro-neurogenicity by Ngfr. Long-term Ngfr expression reduced amyloid plaques and Tau phosphorylation. Postmortem human AD hippocampi and 3D human astroglial cultures showed elevated LCN2 levels correlate with reactive gliosis and reduced neurogenesis. Comparing transcriptional changes in mouse, zebrafish, and human AD brains for cell intrinsic differential gene expression and weighted gene co-expression networks revealed common altered downstream effectors of NGFR signaling, such as PFKP, which can enhance proliferation and neurogenesis in vitro when blocked. Our study suggests that the reactive non-neurogenic astroglia in AD can be coaxed to a pro-neurogenic fate and AD pathology can be alleviated with Ngfr. We suggest that enhancing pro-neurogenic astroglial fate may have therapeutic ramifications in AD.

Abstract Despite decades of genetic studies on late-onset Alzheimer’s disease, the underlying molecular mechanisms remain unclear. To better comprehend its complex etiology, we use an integrative approach to build robust predictive (causal) network models using two large human multi-omics datasets. We delineate bulk-tissue gene expression into single cell-type gene expression and integrate clinical and pathologic traits, single nucleotide variation, and deconvoluted gene expression for the construction of cell type-specific predictive network models. Here, we focus on neuron-specific network models and prioritize 19 predicted key drivers modulating Alzheimer’s pathology, which we then validate by knockdown in human induced pluripotent stem cell-derived neurons. We find that neuronal knockdown of 10 of the 19 targets significantly modulates levels of amyloid-beta and/or phosphorylated tau peptides, most notably JMJD6. We also confirm our network structure by RNA sequencing in the neurons following knockdown of each of the 10 targets, which additionally predicts that they are upstream regulators of REST and VGF. Our work thus identifies robust neuronal key drivers of the Alzheimer’s-associated network state which may represent therapeutic targets with relevance to both amyloid and tau pathology in Alzheimer’s disease.

Abstract Chromatin immunoprecipitation (ChIP) is an antibody-based approach that is frequently utilized in chromatin biology and epigenetics. The challenge in experimental variability by unpredictable nature of usable input amounts from samples and undefined antibody titer in ChIP reaction still remains to be addressed. Here, we introduce a simple and quick method to quantify chromatin inputs and demonstrate its utility for normalizing antibody amounts to the optimal titer in individual ChIP reactions. For a proof of concept, we utilized ChIP-seq validated antibodies against the key enhancer mark, acetylation of histone H3 on lysine 27 (H3K27ac), in the experiments. The results indicate that the titration-based normalization of antibody amounts improves assay outcomes including the consistency among samples both within and across experiments for a broad range of input amounts.

Abstract Background Alzheimer’s disease (AD) is neuropathologically characterized by amyloid-beta (Aβ) plaques and neurofibrillary tangles. The main protein components of these hallmarks include Aβ40, Aβ42, tau, phosphor-tau, and APOE. We hypothesize that genetic variants influence the levels and solubility of these AD-related proteins in the brain; identifying these may provide key insights into disease pathogenesis. Methods Genome-wide genotypes were collected from 441 AD cases, imputed to the haplotype reference consortium (HRC) panel, and filtered for quality and frequency. Temporal cortex levels of five AD-related proteins from three fractions, buffer-soluble (TBS), detergent-soluble (Triton-X = TX), and insoluble (Formic acid = FA), were available for these same individuals. Variants were tested for association with each quantitative biochemical measure using linear regression, and GSA-SNP2 was used to identify enriched Gene Ontology (GO) terms. Implicated variants and genes were further assessed for association with other relevant variables. Results We identified genome-wide significant associations at seven novel loci and the APOE locus. Genes and variants at these loci also associate with multiple AD-related measures, regulate gene expression, have cell-type specific enrichment, and roles in brain health and other neuropsychiatric diseases. Pathway analysis identified significant enrichment of shared and distinct biological pathways. Conclusions Although all biochemical measures tested reflect proteins core to AD pathology, our results strongly suggest that each have unique genetic architecture and biological pathways that influence their specific biochemical states in the brain. Our novel approach of deep brain biochemical endophenotype GWAS has implications for pathophysiology of proteostasis in AD that can guide therapeutic discovery efforts focused on these proteins.

Abstract Background The S209F variant of Abelson Interactor Protein 3 (ABI3) increases risk for Alzheimer’s disease (AD), but little is known about its function in relation to AD pathogenesis. Methods Here, we use a mouse model that is deficient in Abi3 locus to study how the loss of function of Abi3 impacts two cardinal neuropathological hallmarks of AD—amyloid β plaques and tau pathology. Our study employs extensive neuropathological and transcriptomic characterization using transgenic mouse models and adeno-associated virus-mediated gene targeting strategies. Results Analysis of bulk RNAseq data confirmed age-progressive increase in Abi3 levels in rodent models of AD-type amyloidosis and upregulation in AD patients relative to healthy controls. Using RNAscope in situ hybridization, we localized the cellular distribution of Abi3 in mouse and human brains, finding that Abi3 is expressed in both microglial and non-microglial cells. Next, we evaluated Abi3 −/− mice and document that both Abi3 and its overlapping gene, Gngt2, are disrupted in these mice. Using multiple transcriptomic datasets, we show that expression of Abi3 and Gngt2 are tightly correlated in rodent models of AD and human brains, suggesting a tight co-expression relationship. RNAseq of the Abi3-Gngt2−/− mice revealed upregulation of Trem2, Plcg2, and Tyrobp, concomitant with induction of an AD-associated neurodegenerative signature, even in the absence of AD-typical neuropathology. In APP mice, loss of Abi3-Gngt2 resulted in a gene dose- and age-dependent reduction in Aβ deposition. Additionally, in Abi3-Gngt2 −/− mice, expression of a pro-aggregant form of human tau exacerbated tauopathy and astrocytosis. Further, using in vitro culture assays, we show that the AD-associated S209F mutation alters the extent of ABI3 phosphorylation. Conclusions These data provide an important experimental framework for understanding the role of Abi3-Gngt2 function and early inflammatory gliosis in AD. Our studies also demonstrate that inflammatory gliosis could have opposing effects on amyloid and tau pathology, highlighting the unpredictability of targeting immune pathways in AD.

Abstract Risk stratification for hospitalized adults with COVID-19 is essential to inform decisions about individual patients and allocation of resources. So far, risk models for severe COVID outcomes have included age but have not been optimized to best serve the needs of either older or younger adults. Models also need to be updated to reflect improvements in COVID-19 treatments. This retrospective study analyzed data from 6906 hospitalized adults with COVID-19 from a community health system across five states in the western United States. Risk models were developed to predict mechanical ventilation illness or death across one to 56 days of hospitalization, using clinical data available within the first hour after either admission with COVID-19 or a first positive SARS-CoV-2 test. For the seven-day interval, models for age ≥ 18 and

Abstract Genetics play an important role in late-onset Alzheimer’s Disease (AD) etiology and dozens of genetic variants have been implicated in AD risk through large-scale GWAS meta-analyses. However, the precise mechanistic effects of most of these variants have yet to be determined. Deeply phenotyped cohort data can reveal physiological changes associated with genetic risk for AD across an age spectrum that may provide clues to the biology of the disease. We utilized over 2000 high-quality quantitative measurements obtained from blood of 2831 cognitively normal adult clients of a consumer-based scientific wellness company, each with CLIA-certified whole-genome sequencing data. Measurements included: clinical laboratory blood tests, targeted chip-based proteomics, and metabolomics. We performed a phenome-wide association study utilizing this diverse blood marker data and 25 known AD genetic variants and an AD-specific polygenic risk score (PGRS), adjusting for sex, age, vendor (for clinical labs), and the first four genetic principal components; sex-SNP interactions were also assessed. We observed statistically significant SNP-analyte associations for five genetic variants after correction for multiple testing (for SNPs in or near NYAP1, ABCA7, INPP5D, and APOE), with effects detectable from early adulthood. The ABCA7 SNP and the APOE2 and APOE4 encoding alleles were associated with lipid variability, as seen in previous studies; in addition, six novel proteins were associated with the e2 allele. The most statistically significant finding was between the NYAP1 variant and PILRA and PILRB protein levels, supporting previous functional genomic studies in the identification of a putative causal variant within the PILRA gene. We did not observe associations between the PGRS and any analyte. Sex modified the effects of four genetic variants, with multiple interrelated immune-modulating effects associated with the PICALM variant. In post-hoc analysis, sex-stratified GWAS results from an independent AD case–control meta-analysis supported sex-specific disease effects of the PICALM variant, highlighting the importance of sex as a biological variable. Known AD genetic variation influenced lipid metabolism and immune response systems in a population of non-AD individuals, with associations observed from early adulthood onward. Further research is needed to determine whether and how these effects are implicated in early-stage biological pathways to AD. These analyses aim to complement ongoing work on the functional interpretation of AD-associated genetic variants.

Resources reporting RNA editing sites from brain tissue have been published. Here, the authors provide an atlas of RNA editing events found in the aged and Alzheimer’s disease human brain tissue resulting in changes at protein level.

Abstract Cerebral amyloid angiopathy (CAA) contributes to accelerated cognitive decline in Alzheimer’s disease (AD) dementia and is a common finding at autopsy. The APOEε4 allele and male sex have previously been reported to associate with increased CAA in AD. To inform biomarker and therapeutic target discovery, we aimed to identify additional genetic risk factors and biological pathways involved in this vascular component of AD etiology. We present a genome-wide association study of CAA pathology in AD cases and report sex- and APOE-stratified assessment of this phenotype. Genome-wide genotypes were collected from 853 neuropathology-confirmed AD cases scored for CAA across five brain regions, and imputed to the Haplotype Reference Consortium panel. Key variables and genome-wide genotypes were tested for association with CAA in all individuals and in sex and APOEε4 stratified subsets. Pathway enrichment was run for each of the genetic analyses. Implicated loci were further investigated for functional consequences using brain transcriptome data from 1,186 samples representing seven brain regions profiled as part of the AMP-AD consortium. We confirmed association of male sex, AD neuropathology and APOEε4 with increased CAA, and identified a novel locus, LINC-PINT, associated with lower CAA amongst APOEε4-negative individuals (rs10234094-C, beta = −3.70 [95% CI −0.49—−0.24]; p = 1.63E-08). Transcriptome profiling revealed higher LINC-PINT expression levels in AD cases, and association of rs10234094-C with altered LINC-PINT splicing. Pathway analysis indicates variation in genes involved in neuronal health and function are linked to CAA in AD patients. Further studies in additional and diverse cohorts are needed to assess broader translation of our findings.

Abstract Introduction Passive immunotherapies targeting Aβ continue to be evaluated as Alzheimer’s disease (AD) therapeutics, but there remains debate over the mechanisms by which these immunotherapies work. Besides the amount of preexisting Aβ deposition and the type of deposit (compact or diffuse), there is little data concerning what factors, independent of those intrinsic to the antibody, might influence efficacy. Here we (i) explored how constitutive priming of the underlying innate activation states by Il10 and Il6 might influence passive Aβ immunotherapy and (ii) evaluated transcriptomic data generated in the AMP-AD initiative to inform how these two cytokines and their receptors’ mRNA levels are altered in human AD and an APP mouse model. Methods rAAV2/1 encoding EGFP, Il6 or Il10 were delivered by somatic brain transgenesis to neonatal (P0) TgCRND8 APP mice. Then, at 2 months of age, the mice were treated bi-weekly with a high-affinity anti-Aβ1–16 mAb5 monoclonal antibody or control mouse IgG until 6 months of age. rAAV mediated transgene expression, amyloid accumulation, Aβ levels and gliosis were assessed. Extensive transcriptomic data was used to evaluate the mRNA expression levels of IL10 and IL6 and their receptors in the postmortem human AD temporal cortex and in the brains of TgCRND8 mice, the later at multiple ages. Results Priming TgCRND8 mice with Il10 increases Aβ loads and blocks efficacy of subsequent mAb5 passive immunotherapy, whereas priming with Il6 priming reduces Aβ loads by itself and subsequent Aβ immunotherapy shows only a slightly additive effect. Transcriptomic data shows that (i) there are significant increases in the mRNA levels of Il6 and Il10 receptors in the TgCRND8 mouse model and temporal cortex of humans with AD and (ii) there is a great deal of variance in individual mouse brain and the human temporal cortex of these interleukins and their receptors. Conclusions The underlying immune activation state can markedly affect the efficacy of passive Aβ immunotherapy. These results have important implications for ongoing human AD immunotherapy trials, as they indicate that underlying immune activation states within the brain, which may be highly variable, may influence the ability for passive immunotherapy to alter Aβ deposition.

Alzheimer’s disease (AD) is typically associated with hippocampal and cortical pathology, although hippocampal sparing and limbic predominant forms exist. The authors use transcriptomic analysis and neuropathology to identify genes associated with selective hippocampal vulnerability in AD.

Abstract SORL1 is strongly associated with both sporadic and familial forms of Alzheimer’s disease (AD), but a lack of information about alternatively spliced transcripts currently limits our understanding of the role of SORL1 in AD. Here, we describe a SORL1 transcript (SORL1-38b) characterized by inclusion of a novel exon (E38b) that encodes a truncated protein. We identified E38b-containing transcripts in several brain regions, with the highest expression in the cerebellum and showed that SORL1-38b is largely located in neuronal dendrites, which is in contrast to the somatic distribution of transcripts encoding the full-length SORLA protein (SORL1-fl). SORL1-38b transcript levels were significantly reduced in AD cerebellum in three independent cohorts of postmortem brains, whereas no changes were observed for SORL1-fl. A trend of lower 38b transcript level in cerebellum was found for individuals carrying the risk variant at rs2282649 (known as SNP24), although not reaching statistical significance. These findings suggest synaptic functions for SORL1-38b in the brain, uncovering novel aspects of SORL1 that can be further explored in AD research.

Abstract Background There are challenges in generating mRNA-Seq data from whole-blood derived RNA as globin gene and rRNA are frequent contaminants. Given the abundance of erythrocytes in whole blood, globin genes comprise some 80% or more of the total RNA. Therefore, depletion of globin gene RNA and rRNA are critical steps required to have adequate coverage of reads mapping to the reference transcripts and thus reduce the total cost of sequencing. In this study, we directly compared the performance of probe hybridization (GLOBINClear Kit and Globin-Zero Gold rRNA Removal Kit) and RNAse-H enzymatic depletion (NEBNext® Globin & rRNA Depletion Kit and Ribo-Zero Plus rRNA Depletion Kit) methods from 1 μg of whole blood-derived RNA on mRNA-Seq profiling. All RNA samples were treated with DNaseI for additional cleanup before the depletion step and were processed for poly-A selection for library generation. Results Probe hybridization revealed a better overall performance than the RNAse-H enzymatic depletion method, detecting a higher number of genes and transcripts without 3′ region bias. After depletion, samples treated with probe hybridization showed globin genes at 0.5% (±0.6%) of the total mapped reads; the RNAse-H enzymatic depletion had 3.2% (±3.8%). Probe hybridization showed more junction reads and transcripts compared with RNAse-H enzymatic depletion and also had a higher correlation (R > 0.9) than RNAse-H enzymatic depletion (R > 0.85). Conclusion In this study, our results showed that 1 μg of high-quality RNA from whole blood could be routinely used for transcriptional profiling analysis studies with globin gene and rRNA depletion pre-processing. We also demonstrated that the probe hybridization depletion method is better suited to mRNA sequencing analysis with minimal effect on RNA quality during depletion procedures.

The limited understanding of the temporal molecular changes in late-onset Alzheimer’s disease hinder the development of therapeutic treatment. The authors use manifold learning to develop a molecular model for disease progression from RNASeq data from human postmortem brain samples.

APOE4 is a strong genetic risk factor for late-onset Alzheimer’s disease. Here, the authors show that APOE4 is associated with AD features in hiPSCs-derived cerebral organoids. Isogenic conversion of APOE4 to APOE3 attenuates the AD-associated phenotype.

Abstract Missense variants ABI3_rs616338-T and PLCG2_rs72824905-G were previously associated with elevated or reduced risk of Alzheimer’s disease (AD), respectively. Despite reports of associations with other neurodegenerative diseases, there are few studies of these variants in purely neuropathologically diagnosed cohorts. Further, the effect of these mutations on neurodegenerative disease pathologies is unknown. In this study, we tested the effects of ABI3_rs616338-T and PLCG2_rs72824905-G on disease risk in autopsy cohorts comprised of 973 patients diagnosed neuropathologically with Lewy body disease (LBD-NP) and 1040 with progressive supranuclear palsy (PSP), compared to 3351 controls. LBD-NP patients were further categorized as high, intermediate and low likelihood of clinical dementia with Lewy bodies (DLB-CL) based on DLB Consortium criteria. We also tested for association with both Braak neurofibrillary tau tangle (nTotal = 2008, nPSP = 1037, nLBD-NP = 971) and Thal phase amyloid plaque scores (nTotal = 1786, nPSP = 1018, nLBD-NP = 768). Additionally, 841 PSP patients had quantitative tau neuropathology measures that were assessed for genetic associations. There was no statistically significant association with disease risk for either LBD-NP or PSP in our study. LBD intermediate category disease risk was significantly associated with ABI3_rs616338-T (OR = 2.65, 95% CI 1.46–4.83, p = 0.001). PLCG2_rs72824905-G was associated with lower Braak stage (ß = − 0.822, 95% CI − 1.439 to − 0.204, p = 0.009). This effect was more pronounced in the PSP (ß = − 0.995, 95% CI − 1.773 to − 0.218, p = 0.012) than LBD-NP patients (ß = − 0.292, 95% CI − 1.283 to 0.698, p = 0.563). PLCG2_rs72824905-G also showed association with reduced quantitative tau pathology for each lesion type and overall tau burden in PSP (ß = − 0.638, 95% CI − 1.139 to − 0.136, p = 0.013). These findings support a role for PLCG2_rs72824905-G in suppressing tau neuropathology. ABI3_rs616338-T may influence disease risk specifically in the LBD-NP intermediate category comprised of patients with diffuse neocortical or limbic LB, concurrently with moderate or high AD neuropathology, respectively. Our study provides a potential mechanism of action for the missense PLCG2 variant and suggests a differential disease risk effect for ABI3 in a distinct LBD-NP neuropathologic category.

Abstract Large-scale brain bulk-RNAseq studies identified molecular pathways implicated in Alzheimer’s disease (AD), however these findings can be confounded by cellular composition changes in bulk-tissue. To identify cell intrinsic gene expression alterations of individual cell types, we designed a bioinformatics pipeline and analyzed three AD and control bulk-RNAseq datasets of temporal and dorsolateral prefrontal cortex from 685 brain samples. We detected cell-proportion changes in AD brains that are robustly replicable across the three independently assessed cohorts. We applied three different algorithms including our in-house algorithm to identify cell intrinsic differentially expressed genes in individual cell types (CI-DEGs). We assessed the performance of all algorithms by comparison to single nucleus RNAseq data. We identified consensus CI-DEGs that are common to multiple brain regions. Despite significant overlap between consensus CI-DEGs and bulk-DEGs, many CI-DEGs were absent from bulk-DEGs. Consensus CI-DEGs and their enriched GO terms include genes and pathways previously implicated in AD or neurodegeneration, as well as novel ones. We demonstrated that the detection of CI-DEGs through computational deconvolution methods is promising and highlight remaining challenges. These findings provide novel insights into cell-intrinsic transcriptional changes of individual cell types in AD and may refine discovery and modeling of molecular targets that drive this complex disease.

Abstract The discovery of mutations associated with familial forms of Alzheimer’s disease (AD), has brought imperative insights into basic mechanisms of disease pathogenesis and progression and has allowed researchers to create animal models that assist in the elucidation of the molecular pathways and development of therapeutic interventions. Position 717 in the amyloid precursor protein (APP) is a hotspot for mutations associated with autosomal dominant AD (ADAD) and the valine to isoleucine amino acid substitution (V717I) at this position was among the first ADAD mutations identified, spearheading the formulation of the amyloid cascade hypothesis of AD pathogenesis. While this mutation is well described in multiple kindreds and has served as the basis for the generation of widely used animal models of disease, neuropathologic data on patients carrying this mutation are scarce. Here we present the detailed clinical and neuropathologic characterization of an APP V717I carrier, which reveals important novel insights into the phenotypic variability of ADAD cases. While age at onset, clinical presentation and widespread parenchymal beta-amyloid (Aβ) deposition are in line with previous reports, our case also shows widespread and severe cerebral amyloid angiopathy (CAA). This patient also presented with TDP-43 pathology in the hippocampus and amygdala, consistent with limbic predominant age-related TDP-43 proteinopathy (LATE). The APOE ε2/ε3 genotype may have been a major driver of the prominent vascular pathology seen in our case. These findings highlight the importance of neuropathologic examinations of genetically determined AD cases and demonstrate striking phenotypic variability in ADAD cases.

Vast numbers of differentially expressed genes and perturbed networks have been identified in Alzheimer’s disease (AD), however, neither disease nor brain region specificity of these transcriptome alterations has been explored. Using RNA-Seq data from 231 temporal cortex and 224 cerebellum samples from patients with AD and progressive supranuclear palsy (PSP), a tauopathy, we identified a striking correlation in the directionality and magnitude of gene expression changes between these 2 neurodegenerative proteinopathies. Further, the transcriptomic changes in AD and PSP brains ware highly conserved between the temporal and cerebellar cortices, indicating that highly similar transcriptional changes occur in pathologically affected and grossly less affected, albeit functionally connected, areas of the brain. Shared up- or downregulated genes in AD and PSP are enriched in biological pathways. Many of these genes also have concordant protein changes and evidence of epigenetic control. These conserved transcriptomic alterations of 2 distinct proteinopathies in brain regions with and without significant gross neuropathology have broad implications. AD and other neurodegenerative diseases are likely characterized by common disease or compensatory pathways with widespread perturbations in the whole brain. These findings can be leveraged to develop multifaceted therapies and biomarkers that address these common, complex, and ubiquitous molecular alterations in neurodegenerative diseases.

Genetic studies have shifted to sequencing-based rare variants discovery after decades of success in identifying common disease variants by Genome-Wide Association Studies using Single Nucleotide Polymorphism chips. Sequencing-based studies require large sample sizes for statistical power and therefore often inadvertently introduce batch effects because samples are typically collected, processed, and sequenced at multiple centers. Conventionally, batch effects are first detected and visualized using Principal Components Analysis and then controlled by including batch covariates in the disease association models. For sequencing-based genetic studies, because all variants included in the association analyses have passed sequencing-related quality control measures, this conventional approach treats every variant as equal and ignores the substantial differences still remaining in variant qualities and characteristics such as genotype quality scores, alternative allele fractions (fraction of reads supporting alternative allele at a variant position) and sequencing depths. In the Alzheimer's Disease Sequencing Project (ADSP) exome dataset of 9,904 cases and controls, we discovered hidden variant-level differences between sample batches of three sequencing centers and two exome capture kits. Although sequencing centers were included as a covariate in our association models, we observed differences at the variant level in genotype quality and alternative allele fraction between samples processed by different exome capture kits that significantly impacted both the confidence of variant detection and the identification of disease-associated variants. Furthermore, we found that a subset of top disease-risk variants came exclusively from samples processed by one exome capture kit that was more effective at capturing the alternative alleles compared to the other kit. Our findings highlight the importance of additional variant-level quality control for large sequencing-based genetic studies. More importantly, we demonstrate that automatically filtering out variants with batch differences may lead to false negatives if the batch discordances come largely from quality differences and if the batch-specific variants have better quality.

Summary: We present a consensus atlas of the human brain transcriptome in Alzheimer’s disease (AD), based on meta-analysis of differential gene expression in 2,114 postmortem samples. We discover 30 brain coexpression modules from seven regions as the major source of AD transcriptional perturbations. We next examine overlap with 251 brain differentially expressed gene sets from mouse models of AD and other neurodegenerative disorders. Human-mouse overlaps highlight responses to amyloid versus tau pathology and reveal age- and sex-dependent expression signatures for disease progression. Human coexpression modules enriched for neuronal and/or microglial genes broadly overlap with mouse models of AD, Huntington’s disease, amyotrophic lateral sclerosis, and aging. Other human coexpression modules, including those implicated in proteostasis, are not activated in AD models but rather following other, unexpected genetic manipulations. Our results comprise a cross-species resource, highlighting transcriptional networks altered by human brain pathophysiology and identifying correspondences with mouse models for AD preclinical studies.

Abstract Background Rare coding variants ABI3_rs616338-T and PLCG2_rs72824905-G were identified as risk or protective factors, respectively, for Alzheimer’s disease (AD). Methods We tested the association of these variants with five neurodegenerative diseases in Caucasian case-control cohorts: 2742 AD, 231 progressive supranuclear palsy (PSP), 838 Parkinson’s disease (PD), 306 dementia with Lewy bodies (DLB) and 150 multiple system atrophy (MSA) vs. 3351 controls; and in an African-American AD case-control cohort (181 AD, 331 controls). 1479 AD and 1491 controls were non-overlapping with a prior report. Results Using Fisher’s exact test, there was significant association of both ABI3_rs616338-T (OR = 1.41, p = 0.044) and PLCG2_rs72824905-G (OR = 0.56, p = 0.008) with AD. These OR estimates were maintained in the non-overlapping replication AD-control analysis, albeit at reduced significance (ABI3_rs616338-T OR = 1.44, p = 0.12; PLCG2_rs72824905-G OR = 0.66, p = 0.19). None of the other cohorts showed significant associations that were concordant with those for AD, although the DLB cohort had suggestive findings (Fisher’s test: ABI3_rs616338-T OR = 1.79, p = 0.097; PLCG2_rs72824905-G OR = 0.32, p = 0.124). PLCG2_rs72824905-G showed suggestive association with pathologically-confirmed MSA (OR = 2.39, p = 0.050) and PSP (OR = 1.97, p = 0.061), although in the opposite direction of that for AD. We assessed RNA sequencing data from 238 temporal cortex (TCX) and 224 cerebellum (CER) samples from AD, PSP and control patients and identified co-expression networks, enriched in microglial genes and immune response GO terms, and which harbor PLCG2 and/or ABI3. These networks had higher expression in AD, but not in PSP TCX, compared to controls. This expression association did not survive adjustment for brain cell type population changes. Conclusions We validated the associations previously reported with ABI3_rs616338-T and PLCG2_rs72824905-G in a Caucasian AD case-control cohort, and observed a similar direction of effect in DLB. Conversely, PLCG2_rs72824905-G showed suggestive associations with PSP and MSA in the opposite direction. We identified microglial gene-enriched co-expression networks with significantly higher levels in AD TCX, but not in PSP, a primary tauopathy. This co-expression network association appears to be driven by microglial cell population changes in a brain region affected by AD pathology. Although these findings require replication in larger cohorts, they suggest distinct effects of the microglial genes, ABI3 and PLCG2 in neurodegenerative diseases that harbor significant vs. low/no amyloid ß pathology.

Abstract Background Single nucleotide polymorphisms (SNPs) inherited as one of two common haplotypes at the transmembrane protein 106B (TMEM106B) locus are associated with the risk of multiple neurodegenerative diseases, including frontotemporal lobar degeneration with pathological inclusions of TDP-43. Among the associated variants, rs3173615 (encoding p.T185S) is the only coding variant; however, non-coding variants may also contribute to disease risk. It has been reported that the risk haplotype is associated with higher levels of TMEM106B and increased levels of TMEM106B cause cytotoxicity; however, the precise mechanism through which TMEM106B haplotypes contribute to neurodegeneration is unclear. Methods We utilized RNA sequencing data derived from temporal cortex (TCX) and cerebellum (CER) from 312 North American Caucasian subjects neuropathologically diagnosed with Alzheimer’s disease, progressive supranuclear palsy, pathological aging or normal controls to analyze transcriptome signatures associated with the risk (TT) and protective (SS) TMEM106B haplotypes. In cohorts matched for disease phenotype, we used Analysis of Variance (ANOVA) to identify differentially expressed genes and Weighted Gene Co-expression Network Analysis (WGCNA) to identify gene networks associated with the risk and protective TMEM106B haplotypes. Results A total of 110 TCX and 116 CER samples were included in the analyses. When comparing TT to SS carriers, we detected 593 differentially expressed genes in TCX and 7 in CER. Gene co-expression network analyses further showed that in both TCX and CER the SS haplotype was positively correlated with gene networks involved in synaptic transmission, whereas the TT haplotype was positively correlated with gene networks enriched for immune response. Gene expression patterns of 5 cell-type-specific markers revealed significantly reduced expression of the neuronal marker and relative increases in all other cell markers in TT as compared to SS carriers in TCX with a similar but non-significant trend in CER. Conclusions By comparing the common TMEM106B risk and protective haplotypes we identified significant and partly conserved transcriptional differences across TCX and CER and striking changes in cell-type composition, especially in TCX. These findings illustrate the profound effect of TMEM106B haplotypes on brain health and highlight the importance to better understand TMEM106B’s function and dysfunction in the context of neurodegenerative diseases.

Abstract Background After decades of identifying risk factors using array-based genome-wide association studies (GWAS), genetic research of complex diseases has shifted to sequencing-based rare variants discovery. This requires large sample sizes for statistical power and has brought up questions about whether the current variant calling practices are adequate for large cohorts. It is well-known that there are discrepancies between variants called by different pipelines, and that using a single pipeline always misses true variants exclusively identifiable by other pipelines. Nonetheless, it is common practice today to call variants by one pipeline due to computational cost and assume that false negative calls are a small percent of total. Results We analyzed 10,000 exomes from the Alzheimer’s Disease Sequencing Project (ADSP) using multiple analytic pipelines consisting of different read aligners and variant calling strategies. We compared variants identified by using two aligners in 50,100, 200, 500, 1000, and 1952 samples; and compared variants identified by adding single-sample genotyping to the default multi-sample joint genotyping in 50,100, 500, 2000, 5000 and 10,000 samples. We found that using a single pipeline missed increasing numbers of high-quality variants correlated with sample sizes. By combining two read aligners and two variant calling strategies, we rescued 30% of pass-QC variants at sample size of 2000, and 56% at 10,000 samples. The rescued variants had higher proportions of low frequency (minor allele frequency [MAF] 1–5%) and rare (MAF

Abstract Background Simultaneous consideration of two neuropathological traits related to Alzheimer’s disease (AD) has not been attempted in a genome-wide association study. Methods We conducted genome-wide pleiotropy analyses using association summary statistics from the Beecham et al. study (PLoS Genet 10:e1004606, 2014) for AD-related neuropathological traits, including neuritic plaque (NP), neurofibrillary tangle (NFT), and cerebral amyloid angiopathy (CAA). Significant findings were further examined by expression quantitative trait locus and differentially expressed gene analyses in AD vs. control brains using gene expression data. Results Genome-wide significant pleiotropic associations were observed for the joint model of NP and NFT (NP + NFT) with the single-nucleotide polymorphism (SNP) rs34487851 upstream of C2orf40 (alias ECRG4, P = 2.4 × 10−8) and for the joint model of NFT and CAA (NFT + CAA) with the HDAC9 SNP rs79524815 (P = 1.1 × 10−8). Gene-based testing revealed study-wide significant associations (P ≤ 2.0 × 10−6) for the NFT + CAA outcome with adjacent genes TRAPPC12, TRAPPC12-AS1, and ADI1. Risk alleles of proxy SNPs for rs79524815 were associated with significantly lower expression of HDAC9 in the brain (P = 3.0 × 10−3), and HDAC9 was significantly downregulated in subjects with AD compared with control subjects in the prefrontal (P = 7.9 × 10−3) and visual (P = 5.6 × 10−4) cortices. Conclusions Our findings suggest that pleiotropy analysis is a useful approach to identifying novel genetic associations with complex diseases and their endophenotypes. Functional studies are needed to determine whether ECRG4 or HDAC9 is plausible as a therapeutic target.

Abstract Background While age and the APOE ε4 allele are major risk factors for Alzheimer’s disease (AD), a small percentage of individuals with these risk factors exhibit AD resilience by living well beyond 75 years of age without any clinical symptoms of cognitive decline. Methods We used over 200 “AD resilient” individuals and an innovative, pedigree-based approach to identify genetic variants that segregate with AD resilience. First, we performed linkage analyses in pedigrees with resilient individuals and a statistical excess of AD deaths. Second, we used whole genome sequences to identify candidate SNPs in significant linkage regions. Third, we replicated SNPs from the linkage peaks that reduced risk for AD in an independent dataset and in a gene-based test. Finally, we experimentally characterized replicated SNPs. Results Rs142787485 in RAB10 confers significant protection against AD (p value = 0.0184, odds ratio = 0.5853). Moreover, we replicated this association in an independent series of unrelated individuals (p value = 0.028, odds ratio = 0.69) and used a gene-based test to confirm a role for RAB10 variants in modifying AD risk (p value = 0.002). Experimentally, we demonstrated that knockdown of RAB10 resulted in a significant decrease in Aβ42 (p value = 0.0003) and in the Aβ42/Aβ40 ratio (p value = 0.0001) in neuroblastoma cells. We also found that RAB10 expression is significantly elevated in human AD brains (p value = 0.04). Conclusions Our results suggest that RAB10 could be a promising therapeutic target for AD prevention. In addition, our gene discovery approach can be expanded and adapted to other phenotypes, thus serving as a model for future efforts to identify rare variants for AD and other complex human diseases.

A Correction to this paper has been published: https://doi.org/10.1038/s41467-020-20261-6

An amendment to this paper has been published and can be accessed via the original article.

Correction The original version of this article [1] unfortunately contained a typographical error. The ‘Alzheimer’s Disease Neuroimaging Initiative’ was erroneously included as ‘Alzheimer’s Disease Neuroimaging Initative’ in the author list of the article.

Leucine rich repeat transmembrane protein 3 (LRRTM3) is member of a synaptic protein family. LRRTM3 is a nested gene within α-T catenin (CTNNA3) and resides at the linkage peak for late-onset Alzheimer's disease (LOAD) risk and plasma amyloid β (Aβ) levels. In-vitro knock-down of LRRTM3 was previously shown to decrease secreted Aβ, although the mechanism of this is unclear. In SH-SY5Y cells overexpressing APP and transiently transfected with LRRTM3 alone or with BACE1, we showed that LRRTM3 co-localizes with both APP and BACE1 in early endosomes, where BACE1 processing of APP occurs. Additionally, LRRTM3 co-localizes with APP in primary neuronal cultures from Tg2576 mice transduced with LRRTM3-expressing adeno-associated virus. Moreover, LRRTM3 co-immunoprecipitates with both endogenous APP and overexpressed BACE1, in HEK293T cells transfected with LRRTM3. SH-SY5Y cells with knock-down of LRRTM3 had lower BACE1 and higher CTNNA3 mRNA levels, but no change in APP. Brain mRNA levels of LRRTM3 showed significant correlations with BACE1, CTNNA3 and APP in ∼400 humans, but not in LRRTM3 knock-out mice. Finally, we assessed 69 single nucleotide polymorphisms (SNPs) within and flanking LRRTM3 in 1,567 LOADs and 2,082 controls and identified 8 SNPs within a linkage disequilibrium block encompassing 5'UTR-Intron 1 of LRRTM3 that formed multilocus genotypes (MLG) with suggestive global association with LOAD risk (p = 0.06), and significant individual MLGs. These 8 SNPs were genotyped in an independent series (1,258 LOADs and 718 controls) and had significant global and individual MLG associations in the combined dataset (p = 0.02-0.05). Collectively, these results suggest that protein interactions between LRRTM3, APP and BACE1, as well as complex associations between mRNA levels of LRRTM3, CTNNA3, APP and BACE1 in humans might influence APP metabolism and ultimately risk of AD.

Genetic variants that modify brain gene expression may also influence risk for human diseases. We measured expression levels of 24,526 transcripts in brain samples from the cerebellum and temporal cortex of autopsied subjects with Alzheimer's disease (AD, cerebellar n=197, temporal cortex n=202) and with other brain pathologies (non-AD, cerebellar n=177, temporal cortex n=197). We conducted an expression genome-wide association study (eGWAS) using 213,528 cisSNPs within ± 100 kb of the tested transcripts. We identified 2,980 cerebellar cisSNP/transcript level associations (2,596 unique cisSNPs) significant in both ADs and non-ADs (q

Abstract Background Alzheimer’s disease (AD) is characterized by the presence of neurofibrillary tangles made of hyperphosphorylated tau and senile plaques composed of beta-amyloid. These pathognomonic deposits have been implicated in the pathogenesis, although the molecular mechanisms and consequences remain undetermined. UFM1 is an important, but understudied ubiquitin-like protein that is covalently attached to substrates. UFMylation has recently been identified as major modifier of tau aggregation upon seeding in experimental models. However, potential alterations of the UFM1 pathway in human AD brain have not been investigated yet. Methods Here we used frontal and temporal cortex samples from individuals with or without AD to measure the protein levels of the UFMylation pathway in human brain. We used multivariable regression analyses followed by Bonferroni correction for multiple testing to analyze associations of the UFMylation pathway with neuropathological characteristics, primary biochemical measurements of tau and additional biochemical markers from the same cases. We further studied associations of the UFMylation cascade with cellular stress pathways using Spearman correlations with bulk RNAseq expression data and functionally validated these interactions using gene-edited neurons that were generated by CRISPR-Cas9. Results Compared to controls, human AD brain had increased protein levels of UFM1. Our data further indicates that this increase mainly reflects conjugated UFM1 indicating hyperUFMylation in AD. UFMylation was strongly correlated with pathological tau in both AD-affected brain regions. In addition, we found that the levels of conjugated UFM1 were negatively correlated with soluble levels of the deUFMylation enzyme UFSP2. Functional analysis of UFM1 and/or UFSP2 knockout neurons revealed that the DNA damage response as well as the unfolded protein response are perturbed by changes in neuronal UFM1 signaling. Conclusions There are marked changes in the UFMylation pathway in human AD brain. These changes are significantly associated with pathological tau, supporting the idea that the UFMylation cascade might indeed act as a modifier of tau pathology in human brain. Our study further nominates UFSP2 as an attractive target to reduce the hyperUFMylation observed in AD brain but also underscores the critical need to identify risks and benefits of manipulating the UFMylation pathway as potential therapeutic avenue for AD.