Bergman, Jeremiah

Abstract Identification of gene expression changes in post-mortem brain tissue of Alzheimer’s disease donors compared to controls have implicated numerous biological pathways for Alzheimer’s disease pathophysiology. Nonetheless, there is still limited understanding of how gene expression dysregulation underpins specific proteinopathies core to Alzheimer’s disease. Here we investigate brain transcriptomic changes in a well characterized cohort of Alzheimer’s disease donors to identify genes and networks that associate with Alzheimer’s disease endophenotypes including neuropathology measures (Braak stage, Thal phase and cerebral amyloid angiopathy score) and Alzheimer’s disease-related brain protein levels (Apolipoprotein E, Amyloid-β 40, Amyloid-β 42, tau, and phospho-Tau). Bulk transcriptome measures were collected from the temporal cortex tissue of 477 Alzheimer’s disease donors. Following quality control, transcriptome-wide association studies were performed for each endophenotype. We used weighted gene co-expression network analysis to build co-expression networks and integrated transcriptome with epigenetic and genetic data from the same donors. We detected a total of 5,740 Bonferroni significant temporal cortex gene associations with Alzheimer’s disease endophenotypes, most of which were with brain tau levels. We discovered tau-associated co-expression modules enriched in known and novel Alzheimer’s disease pathways. We found that a beneficial (or neutral) brain biochemical state of higher total tau and lower phospho-Tau are associated with increased levels of synaptic, DNA damage/repair, nucleic acid metabolism and myelin processes. In contrast, in a detrimental state of lower total and higher phospho-Tau, there is upregulation of vascular and immune, and downregulation of mitochondrial and myelin pathways. There are brain gene expression perturbations that are associated with Alzheimer’s disease endophenotypes. While some of these associations are common across multiple endophenotypes, many are distinct for different Alzheimer’s disease-related proteins. Based on these findings, we propose a hypothetical model of dynamic brain gene expression changes that track with progressive Alzheimer’s disease proteostasis. These expression changes hold potential to serve as dynamic, precision biomarkers of brain Alzheimer’s disease progression. This study demonstrates the potential of integrative multiomics and deep Alzheimer’s disease endophenotype analyses in well-characterized brain tissues to uncover with precision the complex biology of Alzheimer’s disease.

Abstract Background Progressive supranuclear palsy (PSP) is a devastating neurodegenerative disorder with rapid progression to death. Multiple efforts to develop a disease‐modifying therapy for PSP have yet to succeed. We previously reported the discovery of three potential therapeutic targets for PSP. These three genes, DDR2, KANK2, and STOM, are consistently upregulated in the brains of PSP patients and mouse models of tauopathy. Reducing their expression in a drosophila model significantly rescued the tau‐induced eye degeneration and tau seeding activity. To translate these targets into the clinic, we are developing antisense oligonucleotides (ASO) as potential therapies. Method We screened candidate ASOs in neuroglioma cell line H4. ASO hits were validated for their potency and safety. We also used iPSC‐derived neurons from PSP patients to test the potency and potential toxicity of our ASOs. To investigate target engagement, therapeutic efficacy and preliminary pharmacodynamic properties, we treated iPSC‐derived midbrain organoids with the ASO leads. Lastly, we will assess the off‐target effects of the ASOs by investigating the levels of mRNA transcripts that share similar sequences as our targets using RNAseq. Result Our in vitro screening identified ASO leads that significantly reduce the target gene expressions at mRNA and protein levels in a dose‐dependent manner. The ASO hits also exhibit low cellular toxicity in the H4 cell line. Upon treating iPSC neurons with ASO hits, we were able to validate the dose‐dependent reduction of target mRNAs. Cell viability analysis indicated low toxicity, suggesting that these ASOs are potent and safe. In midbrain organoid model, DDR2‐targeting ASO showed time‐dependent target engagement and were prioritized for further development. Conclusion We developed candidate therapeutics for PSP using ASOs that silence novel therapeutic targets discovered in our systems based cross‐species screening. We showed that our ASOs have low toxicity and confirmed target engagement at both mRNA and protein levels in vitro. The results suggest a favorable preclinical profile of our ASO leads that warrants additional investigations. Importantly, PSP shares congruent pathophysiological changes with other tauopathies, such as Alzheimer’s disease. These common perturbed pathways indicate that our ASOs can be potentially repurposed for multiple neurodegenerative tauopathies.

Abstract Background Alzheimer's disease (AD) affects all brain cells and has complex genomic and immunological alterations. Previous research discovered missense AD risk or protective variants in microglial genes ABI3 and PLCG2, respectively. Expression levels of these genes are altered in AD and can influence microglial function. This study aims to uncover protective, and risk microglial molecular signatures associated with these variants to determine their role in microglial subtypes and states in AD by single cell expression and functional studies. Method We generated microglia‐enriched snRNAseq data from donors harboring either AD protective PLCG2 or AD risk ABI3 missense mutations, or neither mutation. After standard snRNAseq QC, we performed differential expressed gene (DEG) analysis of all microglial cells between variant carriers and non‐carriers using MAST. We defined protective signature as genes that are both down in ABI3 and up in PLCG2 mutation‐carriers. In contrast, risk signature genes are up in ABI3 and down in PLCG2 mutation‐carriers. We investigated the conservation of protective and risk signatures across multiple datasets, including scRNAseq data from iPSC‐derived microglia cells carrying PLCG2 protective variant, snRNAseq data from AD‐resilient donors, and data from AD and other diagnostic groups across multiple brain regions sourced from external datasets. Result We obtained snRNAseq profiles of 35,000 microglia from AD variant carriers. Our DEG analysis among all microglia cells provided 227 microglial protective and 293 risk signature genes defined by these variants. Using integrated analysis of multiple internal and external datasets, we further narrowed down these signatures. We determined that these high‐confidence protective signature genes are downregulated in early AD, upregulated in late AD brains and positively‐correlated with protective variant load in in vitro models. In contrast, risk signature genes are upregulated in early AD, downregulated in late AD and resilient donors. Risk signature expression is decreased with protective PLCG2 variant load, and altered with Aβ treatment in in vitro models. Conclusion Our study uncovers microglia specific protective and risk signatures associated with AD using sn/scRNAseq datasets from multiple sources and models. These findings nominate novel immune targets and pathways with implications for microglial function in health and disease, and ultimately therapeutic potential.

Abstract Background Alzheimer's disease (AD) affects all brain cells and has complex genomic and immunological alterations. Previous research discovered missense AD risk or protective variants in microglial genes ABI3 and PLCG2 , respectively. Expression levels of these genes are altered in AD and can influence microglial function. This study aims to uncover protective, and risk microglial molecular signatures associated with these variants to determine their role in microglial subtypes and states in AD by single cell expression and functional studies. Method We generated microglia‐enriched snRNAseq data from donors harboring either AD protective PLCG2 or AD risk ABI3 missense mutations, or neither mutation. After standard snRNAseq QC, we performed differential expressed gene (DEG) analysis of all microglial cells between variant carriers and non‐carriers using MAST. We defined protective signature as genes that are both down in ABI3 and up in PLCG2 mutation‐carriers. In contrast, risk signature genes are up in ABI3 and down in PLCG2 mutation‐carriers. We investigated the conservation of protective and risk signatures across multiple datasets, including scRNAseq data from iPSC‐derived microglia cells carrying PLCG2 protective variant, snRNAseq data from AD‐resilient donors, and data from AD and other diagnostic groups across multiple brain regions sourced from external datasets. Result We obtained snRNAseq profiles of 35,000 microglia from AD variant carriers. Our DEG analysis among all microglia cells provided 227 microglial protective and 293 risk signature genes defined by these variants. Using integrated analysis of multiple internal and external datasets, we further narrowed down these signatures. We determined that these high‐confidence protective signature genes are downregulated in early AD, upregulated in late AD brains and positively‐correlated with protective variant load in in vitro models. In contrast, risk signature genes are upregulated in early AD, downregulated in late AD and resilient donors. Risk signature expression is decreased with protective PLCG2 variant load, and altered with Aβ treatment in in vitro models. Conclusion Our study uncovers microglia specific protective and risk signatures associated with AD using sn/scRNAseq datasets from multiple sources and models. These findings nominate novel immune targets and pathways with implications for microglial function in health and disease, and ultimately therapeutic potential.

AbstractBackgroundA complex, multicellular disease with genetic and immunological elements, Alzheimer’s disease (AD) affects millions worldwide. There has been previous research linking AD to the missense variants ABI3‐rs616338‐T and PLCG2‐rs72824905‐G, and the altered expression of these genes has been shown to disrupt microglial function. In our understanding of AD risk and resilience, limited research has been conducted on how these variants affect microglial subtypes and states in AD.MethodsWe previously identified DOCK8 protein as a target for fluorescent activated nuclei sorting (FANS) to enrich microglia nuclei from frozen human brains. Using this enrichment strategy, we generated microglia enriched snRNAseq data from temporal cortex tissue of 30 donors harboring ABI3 or PLCG2 missense mutations, or neither mutation. We introduced PLCG2 variant in either homozygote, and heterozygote forms into an AD patient derived iPSC through CRISPR/Cas9 approach, differentiated these cells into microglia cells (IMGLs) following established protocols, and generated scRNAseq data after treatment with Amyloidβ. We applied standard alignment and quality control pipelines to the snRNAseq and scRNAseq datasets in parallel. Established markers of microglia states and subtypes were utilized to annotate the clusters; analyzed for their hub gene expression using comparisons between cluster differential profiling to determine signature pathways. Differential gene expression was assessed using pseudo bulk and MAST approaches.ResultsWe obtained single nuclei transcriptome profiles of 54,000 frozen human brain nuclei, of which 35,000 are microglia, through snRNAseq, and single cell profiles of 63,000 iMGLs through scRNAseq, using our standard alignment and quality control procedures. Our differential gene expression analysis identified novel genes and pathways which were identified as part of resilience or AD risk networks in microglia subtypes and states. These alterations were validated using orthogonal experimental methods and further explored for their conservation using in vivo external datasets. Pseudotime inference will be completed to probe microglia marker gene expression dynamics, treatment conditions, along with continuous cell‐state changes.ConclusionOur study uncovers microglia subtype specific resilience and risk factors provided by AD related genetic variants using snRNAseq from frozen human brain and scRNAseq from IPSC derived microglia cells. These findings nominate novel targets and pathways with therapeutic potential.

AbstractBackgroundProgressive supranuclear palsy (PSP) is a devastating primary tauopathy with rapid progression to death. Although several therapies currently in the development pipeline show promising safety profiles and robust target engagement, few demonstrated significant efficacy in patients, underscoring the need to interrogate additional targets with novel therapeutic modalities to expand the potential therapeutic arsenal. To diversify the therapeutic avenues for PSP and related tauopathies (e.g. Alzheimer’s disease), we systematically integrated multi‐omics data from human brains of PSP and control donors with cross‐species validation to nominate high‐confidence therapeutic targets. We plan to translate our findings into safe and effective treatments for PSP using antisense oligonucleotides (ASO).MethodWe analyzed brain gene expression profiles in PSP and control individuals at bulk tissue (N = 408) and single‐cell levels (snRNAseq, N = 36). Cell‐type‐specific expression perturbations were systematically prioritized using a cross‐species validation paradigm, including rTG4510 tau mice and tau‐overexpressing Drosophila models. We conducted in vitro screening of ASO candidates against our prioritized target genes. ASO efficacy and toxicity were measured. RNAseq experiments will be performed to interrogate off‐target effects. The safety and efficacy of the lead ASOs will be assessed in PSP‐patient‐derived iPSC models.ResultWe previously reported discovery of novel genes with significant differential expression in PSP brains, characterized their brain cell‐specificity and validated them using snRNAseq data. We prioritized 21 genes using gene expression data from a mouse tauopathy model. Validation of these 21 high‐priority genes using a Drosophila tau model nominated astrocytic STOM, KANK2, and DDR2 as potential therapeutic targets for PSP. Knocking down their expression in Drosophila significantly rescues tau‐mediated neurodegenerative pathology. In vitro screening identified ASO leads that reduced the target expressions at mRNA and protein levels with low cellular toxicity.ConclusionWe developed a systems biology pipeline and nominated STOM, KANK2, and DDR2 as candidate gene targets for PSP. We identified ASOs that modulate these targets’ expression without cellular toxicity, suggesting they may be suitable as potential therapeutic candidates. Importantly, the shared pathophysiology and molecular aberrations between PSP and other tauopathies, such as Alzheimer’s disease, suggests that such therapeutic candidates may be repurposed for multiple neurodegenerative diseases, further accelerating, and streamlining the therapeutic pipeline.

AbstractBackgroundAlzheimer’s disease (AD) is neuropathologically characterized by amyloid‐β (Aβ) plaques and tau neurofibrillary tangles often quantified by Thal phase and Braak stage, respectively. Aβ also frequently deposits in the cerebrovasculature with severity categorized by a cerebral amyloid angiopathy (CAA) score. These and related measures often show high variability within AD suggesting distinct underlying mechanisms. We hypothesize that, within the AD brain, neuropathology and levels of core AD‐related proteins are influenced by variations in DNA methylation (DNAm). To test this, we performed epigenome‐wide association studies (EWAS) using DNAm measures from the temporal cortex (TCX) and cerebellum (CER) with AD‐related neuropathologic measures (Braak, Thal, CAA) and brain biochemical levels of five proteins (apoE, Aβ40, Aβ42, tau, p‐tau).MethodsDNAm from neuropathologically‐confirmed AD cases was measured by reduced representation bisulfite sequencing (RRBS) from 471 TCX samples, 200 of which also had CER RRBS. TCX levels of five AD‐related proteins from three tissue fractions (buffer‐, detergent‐, and in‐soluble) were measured previously by ELISA (Liu 2020). CpG methylation (CpGm) was binned by a 15‐state chromatin model (Kundaje 2015) and averaged into CpGm clusters (rCpGm). rCpGms were tested for association with each AD endophenotype and expression levels of nearby genes through multi‐variable linear regression. Replication of significant rCpGms was performed in two independent datasets with AD‐related endophenotypes (Shireby 2022, De Jager 2014).ResultsOur innovative binning method demonstrated biologically relevant CpGm patterns. We found epigenome‐wide significant associations primarily with tau‐related endophenotypes including 93 that had significant and concordant associations in the replication datasets. These rCpGms also significantly associated with expression of nearby genes showing enrichment in oligodendrocyte marker genes including those related to myelination. In vitro validation of tau effects on oligodendrocyte gene expression through DNAm is ongoing.ConclusionsAlthough all endophenotypes tested are core to AD pathophysiology, our results suggest each has a distinct epigenetic architecture underlying their variability in the AD brain. In particular, we found evidence of DNAm variability associating with brain oligodendrocyte gene expression and TCX tau levels. By discovering AD brain endophenotype‐specific DNAm changes, we can identify core components of complex mechanisms revealing important biological insights into AD pathophysiology.

AbstractBackgroundInter‐cellular communication within the gliovascular unit (GVU) is critical for cerebral blood flow regulation, and maintenance of the blood‐brain‐barrier (BBB) properties. The breakdown of BBB in Alzheimer’s disease (AD) is well‐established, but precise underlying molecular changes remain unclear. Additionally, whether GVU molecular alterations observed in AD brains are also detected in blood from living patients is unknown. Further, these GVU molecular perturbations require further investigation in different model systems to identify both human brain‐specific and cross‐species conserved alterations. In this study, we investigated prioritized GVU molecules altered in AD brains for their conservation in blood and cross‐species model systems.MethodsWe performed single nucleus RNA sequencing (snRNAseq) of temporal cortex tissue in AD and control brains. We analyzed this data to detect cell‐specific GVU molecular perturbations and their interactions. We investigated molecular interactions between vascular and astrocyte clusters, the major cell types of the GVU of the BBB. To determine whether GVU transcriptional alterations detected in the brain are preserved in the blood, existing blood expression, genetic, and neuroimaging data from two longitudinal antemortem cohorts were analyzed. Using model systems, including mouse, drosophila, and zebrafish, we evaluated the cross‐species conservation of the top GVU alterations detected in AD brains.ResultsBrain snRNAseq revealed transcriptional profiles of 6,541 astrocytes and 2,210 vascular cells. The latter formed three distinct vascular clusters characterized as pericytes, endothelia and perivascular fibroblasts. We identified differentially expressed genes and their enriched pathways within these clusters and observed the highest levels of transcriptional changes within pericytes. Vascular targets that interact with astrocytic ligands have biological functions in cell signaling, angiogenesis, amyloid ß metabolism, and cytoskeletal architecture. We discovered that genetic variants influencing blood expression levels of some of the prioritized GVU genes were associated with neuroimaging burden of cerebrovascular disease in living human cohorts. Our ongoing studies in model systems revealed conservation of some of the top prioritized molecular perturbations across species.ConclusionOur findings prioritized by multiscale, cross‐tissue human data revealed GVU perturbations within interacting pericyte and astrocyte molecules, which are conserved across multiple cross‐species models. These results nominate new molecular targets and mechanistic insights for BBB disruptions in AD.

AbstractBackgroundAlzheimer’s disease (AD) is a complex, multicellular disease with immunological and genetic components. Previous studies demonstrated associations between AD and the missense variants ABI3‐rs616338‐T and PLCG2‐rs72824905‐G. The perturbed expression of these genes has been demonstrated to disrupt microglial function. To date, the transcriptomic profiles of variant carriers in the context of microglial subtype and state have not been explored in human brain. This represents a knowledge gap in our understanding of AD risk and resilience. We have developed and optimized a single nucleus RNAseq (snRNAseq) approach to study the impact of ABI3 and PLCG2 variants on microglial enriched nuclei from frozen human brain.MethodDOCK8 was identified as a nuclear microglial marker based on previous brain snRNAseq data. DOCK8 based microglial enrichment using fluorescent activated nuclei sorting (FANS) was validated using RT‐qPCR and snRNAseq in human cerebellum and temporal cortex (TCX), respectively. We are applying this enrichment and snRNAseq approach to TCX tissue from a cohort of 30 donors harboring ABI3 or PLCG2 missense mutations, or neither. Standard snRNAseq alignment and quality control pipelines will be applied, and nuclei annotated based on putative microglial markers. Differential gene expression will be assessed using pseudo bulk and MAST approaches, for association with ABI3 and PLCG2 genotype, diagnosis groups, risk factors, and pathology scores.ResultDOCK8 has several commercially available antibodies, making it an ideal microglia‐selection marker. Using RT‐qPCR in microglia enriched nuclei, we demonstrated an enrichment of microglia marker genes and depletion of other brain cell type markers. This is further supported by snRNA‐seq data from DOCK8‐sorted TCX tissue that yielded a 5x increase in microglia proportion (56% vs. 11% in unselected nuclei), enabling analyses of microglia sub‐populations. Leveraging this approach, we will investigate affected genes and pathways and uncover the resilience and risk mechanisms within microglia sub‐populations in ABI3 and PLCG2 variant‐carriers.ConclusionOur study presents a novel method for enriching microglial nuclei from frozen human brains and enables transcriptomic profiling of microglial subpopulations in ABI3 and PLCG2 missense variant‐carriers. We expect these insights will reveal novel genes and pathways that may play a key role in microglial function in AD.

AbstractBackgroundGene expression changes have been identified in post‐mortem brain tissue of AD cases when compared with controls, implicating various genes and pathways. Such findings provide important insights into molecular dysregulation associated with this disease and nominate therapeutic targets. Here we investigate brain gene and co‐expression network changes in a cohort of AD cases to identify associations with degree of AD neuropathology including Braak stage (Tau), Thal phase (Amyloid beta) and brain tissue levels of five AD‐related proteins (apoE, Aβ40, Aβ42, tau, and p‐tau).MethodGene expression measures were collected from temporal cortex (TCX) tissue of 477 AD cases, available through the Mayo Clinic Brain Bank, using RNA sequencing (Illumina TruSeq mRNA) performed on the HiSeq4000. Following quality control 456 cases remained for analysis. Braak and Thal were provided by an experienced neuropathologist and TCX levels of five AD‐related proteins from three fractions (buffer‐soluble, detergent‐soluble, and insoluble) were measured by ELISA (Liu et al. 2020). Weighted gene coexpression network analysis (WGCNA) R package was used to build co‐expression networks. Genes and networks were tested for association with neuropathology and AD‐related protein levels, adjusting for relevant covariates.ResultA total of 57,353 genes were profiled, of which 19,044 were reliably expressed above background levels and formed 20 co‐expression modules. 36 genes were associated with Braak stage, 10 with Thal phase and up to 9,000 with the AD‐related protein fractions. Brain levels of pTau had the greatest number of associations, with an overrepresentation of genes involved in synaptic signalling and immune response gene ontology biological processes. 16 of the coexpression modules associated with at least one trait after Bonferroni correction for number of tests.ConclusionWe have identified a number of genes and networks associated with different aspects of AD neuropathology and associated proteins, further lending support to the involvement of synaptic and immune biology in AD. However it is not clear whether observed changes cause disease or are a consequence of the neurodegenerative process. Comparison of these findings across different proteins and fractions is expected to reveal common and distinct vulnerabilities associated with amyloid and tau, and different stages of protein aggregation.