Quicksall, Zachary

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.

AbstractBackgroundNeurodegenerative disorders affect different brain regions at different time points and to varying degrees, although the underlying reasons for this variability in regional susceptibility is unknown. Brain regions have different proportions of cell types. This cell proportion variability may be further influenced by presence and degree of brain disorders, including neurodegenerative diseases. We postulate that genetic and other factors that influence brain cell type proportions may also provide insights into brain regional susceptibility to neurodegenerative diseases.MethodsUtilizing bulk transcriptomics data from 7 brain regions obtained from 3 cohorts and corresponding whole genome sequencing (WGS) data, genetic variants were evaluated for association with cell type proportions across brain regions. Utilizing cell type specific marker genes for deconvolution provided by McKenzie et.al, we estimated surrogate variables for proportion (SVP) of various cell types using BRETIGEA. We then evaluated relationship between SVP and AD as well as Age, Sex, APOE ε4 and neuropathology measures. Genetic variants from WGS were then tested for association with SVP to identify distinct and conserved patterns across 6 cell types and 7 brain regions.ResultsAnalysis of SVP by diagnosis, sex, age, APOE ε4 and neuropathology measures identified significant differences in cell proportions in various cell types and brain regions. GWAS of SVP identified 134 loci that reached genome‐wide significance (5e‐8). After conservative adjustment for multiple testing (6 cell types and 7 brain regions; 42 tests), variants proximal to SYNPR and VENTXP7 remained significant (p‐value<1.19E‐09). Variants proximal to SYNPR were associated with increased proportion of oligodendrocytes in cerebellum whereas variants near VENTXP7 were associated with reduced proportion of astrocytes in DLPFC. Variants in SYNPR have previously been associated with neurofibrillary tangles and diffuse plaque measurements while variants in VENTXP7 have been associated with measurement of neurofibrillary tangles, cognitive decline, unipolar depression, and bipolar disorder.ConclusionsIt is likely that there are genetic variants associated with cell proportion in brain, may vary by cell type and brain region, while others may act more broadly. Identification of genetic variants associated with cell proportions in different brain regions may provide insights into affected cells and regional vulnerabilities.

AbstractBackgroundTwo major risk factors of Alzheimer’s disease (AD) are advanced age and the APOE‐ε4 risk allele. However, despite having these risk factors, some individuals escape from developing AD or have slower cognitive decline. In this study we hypothesized that “resilient” individuals, defined as those who remain cognitively normal despite being >80 years‐old and carrying at least one copy of the APOE‐ε4 allele, have concerted changes in their blood transcriptomes that are associated with resilience to AD.MethodWe generated blood RNA‐sequence data from 105 APOE‐ε4 carriers from the Mayo Clinic Study of Aging (MCSA) who remained cognitively normal after age 80 (resilient), or who developed mild cognitive impairment (MCI) or AD (non‐resilient). Weighted gene co‐expression network analysis (WGCNA) was performed to evaluate the association of the blood co‐expression network with primary phenotypes of resilience and memory scores. To validate our findings, we constructed a consensus co‐expression network using both MCSA and a replication cohort of blood microarray data from 91 Alzheimer’s Disease Neuroimaging Initiative (ADNI) study participants who had a similar demographic profile. We annotated the consensus modules for their enrichment of cell types, biological processes, genes within previously published genome‐wide significant (GWS) AD‐associated loci, and genes that reside in consensus modules built from seven brain regions across three different autopsy cohorts. To identify genetic variants regulating blood consensus module co‐expression, we performed module quantitative trait locus (modQTL) analysis.ResultWe found co‐expression modules significantly (q<0.05) correlated with the memory phenotype, logical memory delayed recall, across both MCSA and ADNI cohorts. Some of these modules were enriched with blood cell types and biological processes, including telomere maintenance, mitochondrial functions, and immune response. We found 1 module significantly (q<0.05) enriched with genes near GWS AD‐associated loci and 4 modules with genes in brain consensus modules. Lastly, we found 8 modQTLs significantly (Bonferroni‐p<0.05) associated with the co‐expression of 5 different modules.ConclusionWe identified blood co‐expression modules, biological pathways and specific genetic variants associated with cognitive resilience in the presence of AD risk factors. This study provides novel insights into the potential role of peripheral blood in the development of and resilience to 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.

AbstractBackgroundLoss of function mutations in PINK1 and PRKN are the most common causes of early‐onset Parkinson disease (PD). The encoded ubiquitin (Ub) kinase PINK1 and the E3 Ub ligase PRKN together regulate the autophagic degradation of dysfunctional mitochondria, termed mitophagy. During this cytoprotective process, PINK1 and PRKN identify and jointly decorate damaged mitochondria with phosphorylated Ub (pS65‐Ub) that serves as the mitophagy tag. Accumulating evidence suggest a pathogenic role for impaired mitophagy in PD and related neurodegenerative conditions. This is in line with our previous finding of increased levels of the mitophagy marker pS65‐Ub in the hippocampus of postmortem Lewy body disease (LBD) brain.MethodWe performed a genome‐wide association study (GWAS) of hippocampal pS65‐Ub level to identify novel regulators of the PINK1‐PRKN‐directed mitophagy in LBD. In our two‐stage GWAS, we obtained 1,012 autopsy confirmed LBD samples (754 discovery stage, 258, replication stage) from Caucasian participants of European ancestry. Using an established, mostly automated workflow, hippocampal sections were immunostained for pS65‐Ub, scanned, and quantified with unbiased algorithms.ResultOur findings identified that APOE rs429358 (β: 0.50, 95% CI: 0.41 to 0.69, p = 8.67×10−25) was very strongly associated with pS65‐Ub levels in LBD brain. The APOE4 genotype is known to exacerbate Aβ and tau deposition. We did observe strong correlations of pS65‐Ub levels with both tau tangle and senile plaque density in the entire LBD cohort. Nevertheless, the increased pS65‐Ub levels in APOE4‐carrying LBD brains were also independent of Aβ and tau deposition in a subgroup (N = 175) with medium levels of these pathologies. To further explore APOE4 effects on mitophagy outside disease context, this observation was also confirmed in APOE‐targeted replacement mice and in induced pluripotent stem cell‐derived astrocytes.ConclusionOur findings nominate a novel mitophagy regulator in LBD brain and highlight a strong association of APOE4 with mitophagy alteration. Our results suggested that APOE4 may mediate mitophagy defects through two mechanisms, by aggravating Alzheimer’s‐type pathologies and also independent of them. With APOE4 being the strongest known risk factor for Alzheimer’s disease and dementia with Lewy bodies, our findings suggest a common mechanistic link underscoring the importance of mitochondrial quality control.

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, the identification of which may reveal important biological insights into AD pathophysiology. 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 clusters from the temporal cortex (TCX) and cerebellum (CER) with AD‐related neuropathologic measures (Braak, Thal, and CAA) and brain biochemical levels of five proteins (apoE, Aβ40, Aβ42, tau, and p‐tau).MethodsDNAm from neuropathologically‐confirmed AD cases available through the Mayo Clinic brain bank was measured by reduced representation bisulfite sequencing (RRBS) from 471 brain samples, 200 of which had both TCX and CER measures. Neuropathology was scored by an experienced neuropathologist. TCX levels of five AD‐related proteins from three tissue fractions (buffer‐, detergent‐, and in‐soluble) were measured by ELISA (Liu 2020). CpG methylation (CpGm) was binned by a 15‐state chromatin model (Kundaje 2015) and averaged. Averaged CpGm clusters were tested for their association with each endophenotype through multi‐variable linear regression. The influence of CpGm clusters on expression levels of nearby genes was also investigated through multi‐variable regression.ResultsOur innovative binning method demonstrated expected and biologically relevant CpGm patterns. We identified multiple epigenome‐wide significant CpGm clusters unique to each AD‐related endophenotype and brain region. Comparisons between the TCX and CER further revealed brain region‐specific enrichment of distinct CpGm clusters. Significant CpGm clusters were further characterized by their influence on gene expression.ConclusionsAlthough all neuropathologic and biochemical endophenotypes tested are core to AD pathophysiology, our results suggest each has a distinct epigenetic architecture underlying their variability in the AD brain. By discovering brain region‐specific DNAm changes that contribute to each neuropathologic or biochemical outcome, we can identify core components of complex mechanisms that individually and jointly contribute to 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.

AbstractBackgroundAlzheimer’s disease (AD) is a progressive neurodegenerative disorder neuropathologically characterized by amyloid‐β (Aβ) plaques and tau neurofibrillary tangles often measured by Thal phase and Braak stage, respectively. Frequently, Aβ also deposits in the brain vasculature with severity categorized by a cerebral amyloid angiopathy (CAA) score. Assessment of these and related endophenotypes can reveal biological insights into the phenotypic variability in AD. We hypothesize that, amongst AD patients, severity of neuropathology and brain levels of AD‐related proteins are influenced by variation in DNA methylation (DNAm), the identification of which may uncover biological pathways for AD pathogenesis. To test this, we performed epigenome‐wide association studies (EWAS) using DNAm from the temporal cortex (TCX) and cerebellum (CER) with AD‐related neuropathology measures (Braak, Thal, and CAA) and brain levels of five proteins (apoE, Aβ40, Aβ42, tau, and p‐tau).MethodsDNAm from autopsy‐confirmed AD cases available through the Mayo Clinic Brain Bank was measured by reduced representation bisulfite sequencing (RRBS) from 471 brain samples, 200 of which had both TCX and CER measures. Bisulfite converted reads were sequenced on the Illumina High Seq4000 and quality filtered. Neuropathology scoring was performed by an experienced neuropathologist. TCX levels of five AD‐related proteins from three fractions (buffer‐soluble, detergent‐soluble, and insoluble) were measured by ELISA (Liu et al. 2020). Individual CpG methylation ratios were tested for association with each neuropathological and biochemical measure using the R package CpGassoc.ResultsWe identified multiple epigenome‐wide significant CpGs unique to each AD‐related measure. Comparisons between the TCX and CER revealed region specific CpG associations. The effects of sex, age at death, APOE genotype, and presence of genetic variants for each significant CpG association were also identified.ConclusionsAlthough all neuropathological and biochemical measures tested reflect core AD pathologies, our results demonstrate unique associations of epigenetic factors with these individual phenotypes, suggesting that their variability may be governed by distinct epigenetic processes. Discovery of DNAm changes that underlie specific neuropathological or biochemical outcomes can enhance our understanding of AD pathophysiology and may lead to the identification of novel gene or pathway targets for the different facets of this complex condition.