Transcriptional landscape of human microglia implicates age, sex, and APOE-related immunometabolic pathway perturbations

Microglia have fundamental roles in health and disease; however, effects of age, sex, and genetic factors on human microglia have not been fully explored. We applied bulk and single-cell approaches to comprehensively characterize human microglia transcriptomes and their associations with age, sex, and APOE. We identified a novel microglial signature, characterized its expression in bulk tissue and single-cell microglia transcriptomes. We discovered microglial co-expression network modules associated with age, sex, and APOE-ε4 that are enriched for lipid and carbohydrate metabolism genes. Integrated analyses of modules with single-cell transcriptomes revealed significant overlap between age-associated module genes and both pro-inflammatory and disease-associated microglial clusters. These modules and clusters harbor known neurodegenerative disease genes including APOE, PLCG2, and BIN1. Meta-analyses with published bulk and single-cell microglial datasets further supported our findings. Thus, these data represent a well-characterized human microglial transcriptome resource and highlight age, sex, and APOE-related microglial immunometabolism perturbations with potential relevance in neurodegeneration.

Alzheimer's disease and progressive supranuclear palsy share similar transcriptomic changes in distinct brain regions

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.

Genome-wide analysis identifies a novel LINC-PINT splice variant associated with vascular amyloid pathology in Alzheimer's disease

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.

Plasma Biomarkers of Alzheimer's Disease in African Americans

Background/Objective:

The aim of this study was to determine if plasma concentrations of 5 surrogate markers of Alzheimer’s disease (AD) pathology and neuroinflammation are associated with disease status in African Americans.

Methods:

We evaluated 321 African Americans (159 AD, 162 controls) from the Florida Consortium for African-American Alzheimer’s Disease Studies (FCA³DS). Five plasma proteins reflecting AD neuropathology or inflammation (Aβ₄₂, tau, IL6, IL10, TNFα) were tested for associations with AD, age, sex, APOE and MAPT genotypes, and for pairwise correlations.

Results:

Plasma tau levels were higher in AD when adjusted for biological and technical covariates. APOEɛ4 was associated with lower plasma Aβ₄₂ and tau levels. Older age was associated with higher plasma Aβ₄₂, tau, and TNFα. Females had lower IL10 levels. Inflammatory proteins had strong pairwise correlations amongst themselves and with Aβ₄₂.

Conclusion:

We identified effects of demographic and genetic variants on five potential plasma biomarkers in African Americans. Plasma inflammatory biomarkers and Aβ₄₂ may reflect correlated pathologies and elevated plasma tau may be a biomarker of AD in this population.

Association of ABI3 and PLCG2 missense variants with disease risk and neuropathology in Lewy body disease and progressive supranuclear palsy

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 (n_Total = 2008, n_PSP = 1037, n_LBD-NP = 971) and Thal phase amyloid plaque scores (n_Total = 1786, n_PSP = 1018, n_LBD-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.

Deciphering cellular transcriptional alterations in Alzheimer's disease brains

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.

Divergent brain gene expression patterns associate with distinct cell-specific tau neuropathology traits in progressive supranuclear palsy

Progressive supranuclear palsy (PSP) is a neurodegenerative parkinsonian disorder characterized by tau pathology in neurons and glial cells. Transcriptional regulation has been implicated as a potential mechanism in conferring disease risk and neuropathology for some PSP genetic risk variants. However, the role of transcriptional changes as potential drivers of distinct cell-specific tau lesions has not been explored. In this study, we integrated brain gene expression measurements, quantitative neuropathology traits and genome-wide genotypes from 268 autopsy-confirmed PSP patients to identify transcriptional associations with unique cell-specific tau pathologies. We provide individual transcript and transcriptional network associations for quantitative oligodendroglial (coiled bodies = CB), neuronal (neurofibrillary tangles = NFT), astrocytic (tufted astrocytes = TA) tau pathology, and tau threads and genomic annotations of these findings. We identified divergent patterns of transcriptional associations for the distinct tau lesions, with the neuronal and astrocytic neuropathologies being the most different. We determined that NFT are positively associated with a brain co-expression network enriched for synaptic and PSP candidate risk genes, whereas TA are positively associated with a microglial gene-enriched immune network. In contrast, TA is negatively associated with synaptic and NFT with immune system transcripts. Our findings have implications for the diverse molecular mechanisms that underlie cell-specific vulnerability and disease risk in PSP.

ABI3 and PLCG2 missense variants as risk factors for neurodegenerative diseases in Caucasians and African Americans

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.

Conserved brain myelination networks are altered in Alzheimer's and other neurodegenerative diseases

INTRODUCTION

Comparative transcriptome analyses in Alzheimer's disease (AD) and other neurodegenerative proteinopathies can uncover both shared and distinct disease pathways.

METHODS

We analyzed 940 brain transcriptomes including patients with AD, progressive supranuclear palsy (PSP; a primary tauopathy), and control subjects.

RESULTS

We identified transcriptional coexpression networks implicated in myelination, which were lower in PSP temporal cortex (TCX) compared with AD. Some of these associations were retained even after adjustments for brain cell population changes. These TCX myelination network structures were preserved in cerebellum but they were not differentially expressed in cerebellum between AD and PSP. Myelination networks were downregulated in both AD and PSP, when compared with control TCX samples.

DISCUSSION

Downregulation of myelination networks may underlie both PSP and AD pathophysiology, but may be more pronounced in PSP. These data also highlight conservation of transcriptional networks across brain regions and the influence of cell type changes on these networks.

A candidate regulatory variant at the TREM gene cluster associates with decreased Alzheimer's disease risk and increased TREML1 and TREM2 brain gene expression

INTRODUCTION

We hypothesized that common Alzheimer's disease (AD)-associated variants within the triggering receptor expressed on myeloid (TREM) gene cluster influence disease through gene expression.

METHODS

Expression microarrays on temporal cortex and cerebellum from ∼400 neuropathologically diagnosed subjects and two independent RNAseq replication cohorts were used for expression quantitative trait locus analysis.

RESULTS

A variant within a DNase hypersensitive site 5' of TREM2, rs9357347-C, associates with reduced AD risk and increased TREML1 and TREM2 levels (uncorrected P = 6.3 × 10^-3 and 4.6 × 10^-2, respectively). Meta-analysis on expression quantitative trait locus results from three independent data sets (n = 1006) confirmed these associations (uncorrected P = 3.4 × 10^-2 and 3.5 × 10^-3, Bonferroni-corrected P = 6.7 × 10^-2 and 7.1 × 10^-3, respectively).

DISCUSSION

Our findings point to rs9357347 as a functional regulatory variant that contributes to a protective effect observed at the TREM locus in the International Genomics of Alzheimer's Project genome-wide association study meta-analysis and suggest concomitant increase in TREML1 and TREM2 brain levels as a potential mechanism for protection from AD.

Late-onset Alzheimer's risk variants in memory decline, incident mild cognitive impairment, and Alzheimer's disease

We tested association of nine late-onset Alzheimer's disease (LOAD) risk variants from genome-wide association studies (GWAS) with memory and progression to mild cognitive impairment (MCI) or LOAD (MCI/LOAD) in older Caucasians, cognitively normal at baseline and longitudinally evaluated at Mayo Clinic Rochester and Jacksonville (n>2000). Each variant was tested both individually and collectively using a weighted risk score. APOE-e4 associated with worse baseline memory and increased decline with highly significant overall effect on memory. CLU-rs11136000-G associated with worse baseline memory and incident MCI/LOAD. MS4A6A-rs610932-C associated with increased incident MCI/LOAD and suggestively with lower baseline memory. ABCA7-rs3764650-C and EPHA1-rs11767557-A associated with increased rates of memory decline in subjects with a final diagnosis of MCI/LOAD. PICALM-rs3851179-G had an unexpected protective effect on incident MCI/LOAD. Only APOE-inclusive risk scores associated with worse memory and incident MCI/LOAD. The collective influence of the nine top LOAD GWAS variants on memory decline and progression to MCI/LOAD appears limited. Discovery of biologically functional variants at these loci may uncover stronger effects on memory and incident disease.

Novel late-onset Alzheimer disease loci variants associate with brain gene expression

OBJECTIVE

Recent genome-wide association studies (GWAS) of late-onset Alzheimer disease (LOAD) identified 9 novel risk loci. Discovery of functional variants within genes at these loci is required to confirm their role in Alzheimer disease (AD). Single nucleotide polymorphisms that influence gene expression (eSNPs) constitute an important class of functional variants. We therefore investigated the influence of the novel LOAD risk loci on human brain gene expression.

METHODS

We measured gene expression levels in the cerebellum and temporal cortex of autopsied AD subjects and those with other brain pathologies (∼400 total subjects). To determine whether any of the novel LOAD risk variants are eSNPs, we tested their cis-association with expression of 6 nearby LOAD candidate genes detectable in human brain (ABCA7, BIN1, CLU, MS4A4A, MS4A6A, PICALM) and an additional 13 genes ±100 kb of these SNPs. To identify additional eSNPs that influence brain gene expression levels of the novel candidate LOAD genes, we identified SNPs ±100 kb of their location and tested for cis-associations.

RESULTS

CLU rs11136000 (p = 7.81 × 10(-4)) and MS4A4A rs2304933/rs2304935 (p = 1.48 × 10(-4)-1.86 × 10(-4)) significantly influence temporal cortex expression levels of these genes. The LOAD-protective CLU and risky MS4A4A locus alleles associate with higher brain levels of these genes. There are other cis-variants that significantly influence brain expression of CLU and ABCA7 (p = 4.01 × 10(-5)-9.09 × 10(-9)), some of which also associate with AD risk (p = 2.64 × 10(-2)-6.25 × 10(-5)).

CONCLUSIONS

CLU and MS4A4A eSNPs may at least partly explain the LOAD risk association at these loci. CLU and ABCA7 may harbor additional strong eSNPs. These results have implications in the search for functional variants at the novel LOAD risk loci.

Brain expression genome-wide association study (eGWAS) identifies human disease-associated variants

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<0.05, p = 7.70×10⁻⁵–1.67×10⁻⁸²). Of these, 2,089 were also significant in the temporal cortex (p = 1.85×10⁻⁵–1.70×10⁻¹⁴¹). The top cerebellar cisSNPs had 2.4-fold enrichment for human disease-associated variants (p<10⁻⁶). We identified novel cisSNP/transcript associations for human disease-associated variants, including progressive supranuclear palsy SLCO1A2/rs11568563, Parkinson's disease (PD) MMRN1/rs6532197, Paget's disease OPTN/rs1561570; and we confirmed others, including PD MAPT/rs242557, systemic lupus erythematosus and ulcerative colitis IRF5/rs4728142, and type 1 diabetes mellitus RPS26/rs1701704. In our eGWAS, there was 2.9–3.3 fold enrichment (p<10⁻⁶) of significant cisSNPs with suggestive AD–risk association (p<10⁻³) in the Alzheimer's Disease Genetics Consortium GWAS. These results demonstrate the significant contributions of genetic factors to human brain gene expression, which are reliably detected across different brain regions and pathologies. The significant enrichment of brain cisSNPs among disease-associated variants advocates gene expression changes as a mechanism for many central nervous system (CNS) and non–CNS diseases. Combined assessment of expression and disease GWAS may provide complementary information in discovery of human disease variants with functional implications. Our findings have implications for the design and interpretation of eGWAS in general and the use of brain expression quantitative trait loci in the study of human disease genetics.

Genetic variants that regulate gene expression levels can also influence human disease risk. Discovery of genomic loci that alter brain gene expression levels (brain expression quantitative trait loci = eQTLs) can be instrumental in the identification of genetic risk underlying both central nervous system (CNS) and non–CNS diseases. To systematically assess the role of brain eQTLs in human disease and to evaluate the influence of brain region and pathology in eQTL mapping, we performed an expression genome-wide association study (eGWAS) in 773 brain samples from the cerebellum and temporal cortex of ∼200 autopsied subjects with Alzheimer's disease (AD) and ∼200 with other brain pathologies (non–AD). We identified ∼3,000 significant associations between cisSNPs near ∼700 genes and their cerebellar transcript levels, which replicate in ADs and non–ADs. More than 2,000 of these associations were reproducible in the temporal cortex. The top cisSNPs are enriched for both CNS and non–CNS disease-associated variants. We identified novel and confirmed previous cisSNP/transcript associations for many disease loci, suggesting gene expression regulation as their mechanism of action. These findings demonstrate the reproducibility of the eQTL approach across different brain regions and pathologies, and advocate the combined use of gene expression and disease GWAS for identification and functional characterization of human disease-associated variants.

Glutathione S-transferase omega genes in Alzheimer and Parkinson disease risk, age-at-diagnosis and brain gene expression: an association study with mechanistic implications

BACKGROUND

Glutathione S-transferase omega-1 and 2 genes (GSTO1, GSTO2), residing within an Alzheimer and Parkinson disease (AD and PD) linkage region, have diverse functions including mitigation of oxidative stress and may underlie the pathophysiology of both diseases. GSTO polymorphisms were previously reported to associate with risk and age-at-onset of these diseases, although inconsistent follow-up study designs make interpretation of results difficult. We assessed two previously reported SNPs, GSTO1 rs4925 and GSTO2 rs156697, in AD (3,493 ADs vs. 4,617 controls) and PD (678 PDs vs. 712 controls) for association with disease risk (case-controls), age-at-diagnosis (cases) and brain gene expression levels (autopsied subjects).

RESULTS

We found that rs156697 minor allele associates with significantly increased risk (odds ratio = 1.14, p = 0.038) in the older ADs with age-at-diagnosis > 80 years. The minor allele of GSTO1 rs4925 associates with decreased risk in familial PD (odds ratio = 0.78, p = 0.034). There was no other association with disease risk or age-at-diagnosis. The minor alleles of both GSTO SNPs associate with lower brain levels of GSTO2 (p = 4.7 × 10-11-1.9 × 10-27), but not GSTO1. Pathway analysis of significant genes in our brain expression GWAS, identified significant enrichment for glutathione metabolism genes (p = 0.003).

CONCLUSION

These results suggest that GSTO locus variants may lower brain GSTO2 levels and consequently confer AD risk in older age. Other glutathione metabolism genes should be assessed for their effects on AD and other chronic, neurologic diseases.