Wang, Xue

ABSTRACT Recent development of single-cell technology across multiple omics platforms has provided new ways to obtain holistic views of cells to study disease pathobiology. Alzheimer’s disease (AD) is the most common form of dementia worldwide, yet the detailed understanding of its cellular and molecular mechanisms remains limited. In this study, we analyzed paired single-cell transcriptomic (scRNA-seq) and chromatin accessibility (scATAC-seq) data from the Seattle Alzheimer’s Disease Brain Cell Atlas (SEA-AD) Consortium to investigate the molecular mechanisms of AD at a cell-subpopulation-specific resolution focusing on glial cells. We benchmarked various multi-omics integration methods using diverse metrics and built an analytic workflow that enabled effective batch correction and cross-modality alignment, creating a unified cell state space. Through integrative analysis of 26 human brain samples, we uncovered AD-associated gene expression and pathway changes in glial subpopulations and highlighted important transcriptomic and epigenomic signatures via functional inference and interpretable machine learning paradigms, discovering the profound involvement of the Solute Carrier proteins (SLC) family genes in multiple glial cell types. We also identified glial cell-specific regulatory programs mediated by key transcription factors such as JUN and FOSL2 in astrocytes, the Zinc Finger (ZNF) family genes in microglia, and the SOX family of transcription factors in oligodendrocytes. Our study provides a comprehensive workflow and a high-resolution view of how glial regulatory programs are disrupted in AD. Our findings offer novel insights into disease-related changes in gene regulation and suggest potential targets for further research and therapy.

Abstract Blood-brain barrier (BBB) dysfunction is a key feature of Alzheimer’s disease (AD), particularly in individuals carrying the APOE-ε4 allele. This dysfunction worsens neuroinflammation and hinders the removal of toxic proteins, such as amyloid-beta (Aβ42), from the brain. In post-mortem brain tissues and in animal models, we previously reported that fibronectin accumulates at the BBB predominantly in APOE-ε4 carriers. Furthermore, we found a loss-of-function variant in the fibronectin 1 ( FN1 ) gene significantly reduces aggregated fibronectin levels and decreases AD risk among APOE-ε4 carriers. Yet, the molecular mechanisms downstream of fibronectin at the BBB remain unclear. The extracellular matrix (ECM) plays a crucial role in maintaining BBB homeostasis and orchestrating the interactions between BBB cell types, including endothelia and astrocytes. Understanding the mechanisms affecting the ECM and BBB cell types will be critical for developing effective therapies against AD, especially among APOE-ε4 carriers. Here, we demonstrate that APOE-ε4 , Aβ42, and inflammation drive the induction of FN1 expression in several models including zebrafish, mice, iPSC-derived human 3D astrocyte and 3D cerebrovascular cell cultures, and in human brains. Fibronectin accumulation disrupts astroglial-endothelial interactions and the signalling cascade between vascular endothelial growth factor (VEGF), heparin-binding epidermal growth factor (HBEGF) and Insulin-like growth factor 1 (IGF1). This accumulation of fibronectin in APOE-ε4- associated AD potentiates BBB dysfunction, which strongly implicates reducing fibronectin deposition as a potential therapeutic target for AD. Graphical abstract Accessibility text This image illustrates the effects of different APOE isoforms (ApoE-ε3 and ApoE-ε4) on blood-brain barrier (BBB) integrity, focusing on the molecular interactions between astrocytes and endothelial cells. This figure emphasizes the detrimental effects of ApoE-ε4 on BBB integrity via fibronectin accumulation and altered signaling pathways. The top section provides a schematic overview of the blood-brain barrier, highlighting astrocytes, endothelial cells, and their interface. The left panel represents the ApoE-ε3 condition: Normal fibronectin (FN1) levels support healthy interactions between astrocytes and endothelial cells. Growth factors, including VEGFA, HBEGF, and IGF1, maintain BBB integrity through their respective receptors (VEGFR and EGFR). Green arrows indicate activation of these signaling pathways. The right panel depicts the ApoE-ε4 condition: Elevated fibronectin (FN1) disrupts astrocyte-endothelium interactions. FN1 binds integrins and activates focal adhesion kinase (FAK), inhibiting VEGFA, which is required for endothelial HBEGF that in turn activates IGF1 signaling. Red symbols indicate inhibition of HBEGF, VEGFA, and IGF1 pathways, leading to BBB dysfunction. Highlights APOE-ε4 drives fibronectin deposition in Alzheimer’s, disrupting astrocyte-endothelia interactions. APOE-ε4 and fibronectin co-localize, forming aggregates at blood-brain barrier (BBB). Fibronectin alters the signaling between VEGF, IGF1, and HBEGF impairing BBB function. Reducing fibronectin restores BBB integrity and offsets APOE-ε4 pathology.

AbstractBackgroundTwo main risk factors of Alzheimer’s disease (AD) are aging and APOE‐ε4. However, some individuals remain cognitively normal despite having these risk factors. They are considered “cognitively resilient”. This study aimed to identify molecular factors that confer cognitive resilience in APOE‐ε4 carriers ≥ 80 years of age and may serve as biomarkers.MethodWe applied weighted gene co‐expression network analysis (WGCNA) to generate consensus co‐expression networks from blood of participants in two antemortem cohorts, the Mayo Clinic Study of Aging (MCSA, n=105), and the Alzheimer’s Disease Neuroimaging Initiative (ADNI, n=91), using RNA‐sequencing and microarray data, respectively. We associated these networks with resilience (resilient vs non‐resilient), cognitive endophenotypes and hippocampal volume. Preservation between consensus networks from blood and those derived from postmortem brain tissues of AD and control donors from AMP‐AD (n=1174) was evaluated. We validated the human findings in four AD mouse models. Finally, machine learning models were utilized to discriminate cases (AD+mild cognitive impairment (MCI)) from controls in MCSA, ADNI and ANMerge antemortem cohorts.ResultFour consensus networks were significantly correlated with a memory phenotype (logical memory delayed recall=LMDR) and hippocampal volume in both MCSA and ADNI. Among these, blood expression module M3 was most preserved with the brain transcriptome. M3 was enriched with NDUF hub genes that are involved in the mitochondrial respiratory chain. Expression levels of M3 and many blood NDUFs had significant associations with better LMDR and hippocampal volume. In brain, NDUFs were upregulated in controls compared to AD, and their expression levels were associated with better global cognition and decreased AD neuropathology. Many NDUFs were significantly downregulated in the hippocampus or cortex of AD mice compared to wild‐types. Lastly, models that included blood NDUFs improved diagnostic accuracy of AD+MCI compared to models that only included demographic and risk variables (age, sex, APOE‐ε4) in MCSA, ADNI and ANMerge. In MCSA and ADNI, adding NDUFs’ expression to models that included established blood biomarkers (Aβ42/40, ptau181, NFL) further improved diagnostic accuracy.ConclusionOur results suggest that mitochondrial NDUFs are centrally‐linked peripheral molecular signatures that may be resilience factors against AD and serve as both therapeutic targets and novel diagnostic biomarkers.

AbstractBackgroundSystemic inflammation plays a pivotal role in many chronic diseases including Alzheimer’s disease (AD). Assessing the composition of immune pathways in neurodegenerative diseases can contribute to precision medicine. Using publicly available transcriptomic data, we sought to elucidate transcriptional networks pertinent to inflammatory pathways across brain regions and peripheral blood in AD/mild cognitive impairment (MCI) and peripheral blood in Parkinson’s disease (PD).MethodFor the AD/MCI vs. control dataset, we analyzed bulk‐RNAseq collected from 6 brain regions of donors from ROSMAP, Mayo Clinic, and Mount Sinai School of Medicine (MSSM) brain banks available from the AMP‐AD consortium. Ante‐mortem, blood RNAseq expression data was retrieved from the AMP‐AD Emory Vascular cohort and Mayo Clinic Study of Aging (MCSA). We also collected blood‐derived microarray expression data from the Alzheimer’s Disease Neuroimaging Initiative (ADNI). For the PD vs Control dataset, blood‐derived bulk‐RNAseq from the PDBP and PPMI cohorts were available through the AMP‐PD consortium. Following quality control, normalization, and residual generation to account for biological and technical variables, co‐expression network modules and their enriched pathways were identified using WGCNA within each dataset. Module/trait correlation tests for aging and diagnosis (cases [AD/MCI or PD] vs control) phenotypes were evaluated. Gene Ontology enrichment analyses were executed to identify enriched pathways and brain or blood cell types within the modules. Modules were tested for preservation across cohorts.ResultWe identified conserved immune signatures across brain regions and cohorts. Modules involved in immune response were preserved across all cohorts. Blood consensus modules involved in immune response were preserved in the brain and vice versa. Some immune modules were associated with AD/MCI, PD, and/or aging. Brain immune modules are significantly associated with aging and/or AD. Significant correlations (q<0.05) with PD diagnosis were present. In the MCSA and Emory vascular cohorts, there were no significant (q<0.05) associations between modules and diagnosis, while in ADNI there were nominal (p<0.05) associations.ConclusionPreserved transcriptional immune networks were identified across blood and brain and across two neurodegenerative diseases. Expanding gene co‐expression network analyses to other diseases and integrating additional omics measures and phenotypes can further strengthen these findings to unravel the immune signatures across complex diseases.

AbstractBackgroundA hallmark of Alzheimer’s disease is amyloid‐β (Aβ) deposition in the brain parenchyma as senile plaques which Aß frequently accompanies in the cerebral vasculature as cerebral amyloid angiopathy (CAA). Previous studies demonstrated that the cerebrovascular system is essential for Aβ clearance and interactions exist between cerebrovascular disease and AD. Neuroimaging data, including measures on small vessel disease, amyloid deposits, and microhemorrhages (MCH), can reflect both cerebrovascular disease and AD pathology, thereby serving as intermediate phenotypes of both conditions. We aim to identify genes and co‐expression networks in blood associated with these neuroimaging endophenotypes and the biological pathways enriched in these perturbed transcriptional networks.MethodWe conducted differential gene expression and Weighted Gene Co‐Expression Network analyses (WGCNA) on blood RNA sequencing data of 407 participants from the longitudinally followed Mayo Clinic Study of Aging cohort. The neuroimaging phenotypes include MCH, infarction, and amyloid deposits measured by Pittsburgh compound B positron emission tomography (PiB PET).ResultWGCNA analysis showed a nominally significant association (p<0.05) between several module eigengenes (MEs) and MCH and PiB PET. One module, M42, remained significantly associated with PiB PET after Bonferroni correction for all 43 modules tested. Pathway enrichment analysis of M42 identified two significant pathways (FDR‐adjusted p<0.05) including phosphodiesterase I activity and IgE binding. These pathways contain genes such as MS4A2 encoding the high‐affinity IgE receptor involved in the immune response. Similarly, ENPP2 in the phosphodiesterase I activity pathway represents an additional candidate gene previously implicated in Alzheimer’s disease.ConclusionWe generated co‐expression networks and identified one network module with a statistically significant association between brain amyloid deposits and blood transcriptome data. Further investigations of this network and its enriched pathways are necessary to understand their roles in amyloidosis, CAA, and neurodegeneration. Future studies will include seeking replication of these results in other datasets. Such studies can help elucidate the biological mechanisms of pathogenesis of cerebrovascular damage and amyloid deposits. In addition, in a clinical setting, blood levels of transcripts discovered in this study may serve as predictive or diagnostic biomarkers with mechanistic implications.

AbstractBackgroundGiven the overlap of symptoms between Alzheimer’s disease (AD) and other dementias, it is essential to develop widely accessible biomarkers that will ensure an accurate diagnosis of AD to enable its proper treatment. African Americans (AA) are twice more likely to develop AD than non‐Hispanic whites (NHW) and yet they remain underrepresented in AD research. The few studies conducted in AA suggest that there are different factors that influence AD in this population; therefore, it is possible that different biomarkers and/or therapies may be needed to effectively treat AD in AA. In this study we are testing the hypothesis that genetic variants that show association with AD‐risk in AA and with plasma levels of transcript and proteins may serve as AD biomarkers which are more accessible and less costly than cerebrospinal fluid or neuroimaging biomarkers.MethodsIn 189 AD cases and 183 cognitively unimpaired (CU) AA participants of the Florida Consortium of African American Alzheimer’s Disease Studies (FCA3DS), we performed targeted DNA sequencing of 10 AD‐associated loci. Genetic variants were tested for association with plasma transcripts and plasma total tau protein levels previously measured in this cohort using a nanoString® custom panel (Reddy J, et al. 2022) and Simoa assays (Deniz K, et al. 2021), respectively. In receiver operating characteristic (ROC) analyses, the most significant expression/protein quantitative trait loci (eQTL/pQTL) were added sequentially to a base model that included only age and sex to identify eQTL/pQTL that improved the accuracy to correctly classify AD cases vs. CU.ResultsOf the 5112 variants identified in the targeted sequencing, 71 were nominally associated (p<0.05) with plasma transcript/protein levels and with AD‐risk. ROC analysis of age, sex, eQTL/pQTL in APOE, ABCA7, AKAP9, EPHA1, SORL1, and TLR4, and plasma levels of APP, ABCA7, AKAP9, CD14, CLU, and APOE‐ɛ4 dosage achieved 77.6% area under the curve (AUC) to discriminate AD vs. CU, a 19.9% improvement over a model that only included age and sex.ConclusionsPlasma transcript levels and eQTL/pQTL are promising diagnostic biomarkers for AD that may improve accessibility and reduce costs for a more accurate diagnosis of AD.

AbstractBackgroundRecent studies increased the number of genes that are associated with the co‐existence of vascular pathologies and Alzheimer’s disease (AD). However, the biological functions and the mechanism of action by which they contribute to the disease pathology are still to be further elucidated. Therefore, animal models that allow streamlined large‐throughput functional screens in a biologically relevant manner are essential. Towards this goal, we generated adult zebrafish models of amyloid toxicity to functionally analyze candidate genes and how they relate to disease perturbations.MethodThrough cerebroventricular injection of human amyloid‐beta42 monomers into the cerebrospinal fluid of the adult zebrafish brain, we inflict acute amyloid toxicity. By performing immunohistological analyses, single cell transcriptomics, gene editing for generating gene variants or full knockouts, transient functional knockdown of candidate genes, pharmacological intervention of selected signalling pathways and comparing the findings to human AD brain transcriptomics datasets, we determined the similarities and differences of the effects of amyloid pathology on the vasculature in zebrafish in relation to the human conditions.ResultsWe found that amyloid pathology in zebrafish brain shows histological and molecular similarities to human brains. Integration of single cell transcriptomics showed remarkable similarities in neurons to amyloid in humans and zebrafish. We identified several molecular mechanisms that underlie this pathology‐induced cellular response and translated to mammalian models including rodent brains. By performing genome‐wide association in AD patients, we identified that the FMNL2 links cerebrovascular disease and AD, by regulating the astroglial and blood vessel interactions and controlling the efficient clearance of toxic proteins from the brain. Our new single cell transcriptomics analyses of the gliovascular niche cells of the zebrafish brain revealed potential cellular interactions and the pathways perturbed in vascular disease and AD.ConclusionWe propose zebrafish as a useful model for functional investigation of vascular pathology and AD‐related genes identified in clinical studies by providing in vivo biological knowledge on disease mechanisms, on which novel drug development strategies can be based.

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.