Mere responders no more, microglia nowadays are respected—if not feared—as early actors when a person is developing a neurodegenerative disease. Several presentations at last month’s Alzheimer’s Association International Conference underscored their pre-eminent role within the pathogenetic cascade of Alzheimer’s disease. Using a cell-type-specific version of a polygenic risk score, researchers reported that myriad genetic variations that are active in microglia can exacerbate amyloid deposition. In contrast, AD-related polygenic variation that is active in neurons has more impact on cognition. To chart how genetic variation alters the expression of microglial genes and influences AD and related diseases, scientists introduced the microglial genomic atlas (MiGA), the largest repository of human microglial transcriptomes to date. The atlas helped nail down associations between AD risk variants and microglial genes, which previously had been inferred from studies of other myeloid cell types.

  • Cell-type-specific polygenic risk scores implicate microglial genes in amyloid deposition.
  • The microglia genomic atlas—MiGA—clarifies how risk variants influence expression of microglial genes.
  • Neuron-specific polygenic scores associate with cognition.

Massive genome-wide association studies have uncovered dozens of genetic polymorphisms that influence the risk of Alzheimer’s disease (Feb 2021 news). Most of these individual variants impart only a smidgen of risk, but collectively, they exert a strong influence. Polygenic risk scores are calculated by tallying the effects of multiple risk variants. They have been tied to all manner of AD-related phenotypes, including plaques, tangles, brain atrophy, and cognitive decline (Escott-Price et al., 2015Jul 2016 news; Feb 2019 news). 

Lost in this polygenic fray is the question of exactly which cells these cumulative risk variants influence. After all, most, if not all, AD risk variants regulate expression in specific cell types. This was the main theme of Hyun-Sik Yang, Brigham and Women’s Hospital in Boston. Wondering how genetic variants active in different cell types contribute to different aspects of preclinical AD, Yang accessed baseline data from the A4 study, a secondary prevention trial that is testing solanezumab in people with brain amyloid who were cognitively normal at baseline. All 2,961 A4 participants had consented to genome-wide genotyping.

Yang and colleagues started by calculating an overall polygenic risk score for each participant based on the variants he or she carried. Using risk variants identified in the largest AD GWAS to date, the researchers employed a Bayesian statistical method that accounts for individual variations in genomic architecture to calculate each person’s PRS (Ge et al., 2019). The analysis incorporated contributions from thousands of genetic variants, most of which fell below a stringent threshold for genome-wide significance.

Crucially, for each variant the scientists referenced single-cell RNA sequencing data from postmortem brain samples in order to designate which cell types expressed the putative causal gene (Mathys et al., 2019). This enabled them to derive distinct polygenic risk scores for each of eight different cell types: excitatory neurons, inhibitory neurons, oligodendrocytes, oligodendrocyte progenitors, astrocytes, microglia, endothelial cells, and pericytes. Some genetic variants were active in more than one cell type. Ultimately, each A4 participant was assigned nine different PRSs: an overall one that included all variants, and one for each of the eight cell types.

Which PRS had the strongest ties to baseline amyloid levels among all the A4 participants? The microglial score was the clear winner in this regard, Yang said. Its association with baseline Aβ load was about twice as strong as either neuronal PRS, and topped the overall PRS by a third. This finding cast microglia as intimately involved in amyloid.

Microglia Drive Plaque. Among people in the A4 study, their microglia-specific PRS scores had the strongest association (y axis) with baseline amyloid burden. [Courtesy of Hyun-Sik Yang, Brigham and Women’s Hospital.]

The story was different for cognition. Using participants’ baseline scores on the preclinical AD cognitive composite (PACC) as a gauge, Yang found that more than any other PRS, a person’s excitatory neuron PRS was most strongly tied to his or her cognitive performance at study start. The higher a person’s excitatory neuron PRS, i.e., the more risk variants known to affect genes in these neurons, the lower he or she tended to perform on the test. The strength of this relationship dropped only slightly when adjusted for Aβ burden, suggesting that cumulative AD risk variants active in excitatory neurons influence cognition independently of amyloid. Yang thinks these variants could influence AD risk by weakening neuronal function. Akin to the concept of cognitive reserve, a high neuronal PRS could render neurons more susceptible to a given amount of amyloid.

Julie Williams of Cardiff University in Wales commented that Yang’s findings convincingly show that microglia are part of the causative pathway of AD, especially since A4 participants were at the earliest stages of the disease.

Lianne Reus of Amsterdam University Medical Center also took a polygenic approach to investigate how different cell types influence AD pathogenesis. At AAIC, Reus linked cell-type-specific PRS to fluid biomarkers. She drew genotyping data from more than 617 participants in the Alzheimer’s Disease Neuroimaging Initiative (ADNI) and 918 from the European Medical Information Framework for Alzheimer’s Disease (EMIF-AD) for a total of 525 controls, 709 people with MCI, and 272 with AD. (Shaw et al., 2009Bos et al., 2018). 

To calculate cell-type-specific PRSs for each person, Reus took a more candidate-based approach. Instead of weaving thousands of variants into her analysis, Reus selected only 26. These had risen to genome-wide significance in GWAS, and had detectable levels of RNA in at least one brain-cell type in the BRAIN-RNA-Seq database (de Rojas et al., 2021Zhang et al., 2014). Reus designated a gene as cell-type specific if at least 50 percent of its expression came from a single cell type. Of the 26 variants, 10 were nonspecific, 11 were microglia-specific, four were astrocyte-specific, one was specific to neurons, and none were specific for endothelial cells or oligodendrocytes (see genes in image below). With this minimalist panel, Reus calculated four PRSs for each person: one for the 10 nonspecific loci, and one each for microglia, astrocytes, and even neurons, though this was based on only one variant. She then asked how each person’s four PRSs associated with their CSF concentration of Aβ, total tau, p-tau181, neurogranin, neurofilament light (NfL), and YKL40.

Getting Specific. Of 26 AD risk variants (right), 10 were expressed across cell types, 11 were more expressed in microglia, four in astrocytes, and one was found expressed only in neurons. [Courtesy of Lianne Reus, Amsterdam University.]

ApoE was designated as an astrocyte-specific gene. When included in the PRS for that cell type, the astrocyte polygenic scores strongly associated with CSF Aβ, total tau, p-tau181, and neurogranin. Without ApoE, a person’s astrocyte PRS, based on only three possible variants, associated only with CSF Aβ, demonstrating the power of ApoE’s connection to AD pathogenesis more broadly. The microglial PRS also associated with CSF Aβ, and nominally with total tau. Neither the neuron-specific score nor the nonspecific PRS was significantly tied to any CSF biomarker. Notably, only one variant in SLC24A4 was part of the neuronal “PRS” in this dataset.

Although Yang and Reus’s studies employed different methods to calculate PRS and to designate cell-type specificity of AD risk variants, they drew one common conclusion, namely that markers of amyloid deposition associate with variants expressed in astrocytes and/or microglia.

Glial Ties to Aβ. The astrocyte PRS linked to several CSF biomarkers (right) with high statistical significance. Without ApoE, the astrocyte PRS only linked to CSF Aβ. Microglial PRS associated with CSF Aβ and nominally with total tau. [Courtesy of Lianne Reus, Amsterdam University.]

“Cell-type-specific polygenic risk scores help to clarify the contribution of specific cell types in AD pathophysiology in different stages of the disease,” wrote Betty Tijms of Amsterdam University Medical Center and Pieter-Jelle Visser of Maastricht University in the Netherlands, who co-authored the study with Reus. “The precise cell-type-dependent interactions will probably vary between individuals, and so this approach may help to better understand interindividual variation in disease manifestation.”

Microglia Atlas Connects the Dots
While the central role of microglial variants in AD has become clear, exactly how most of them influence microglial gene expression and cellular function is unclear. That is partly because microglia are so changeable. They adapt quickly to changes in their environment. Genetic variation could therefore affect microglia differently depending on which brain region they are in, their age, or their response to a pathological condition. What’s more, because genetic variants rarely occur within the coding region of a gene, and gene expression varies dramatically by cell type, detailed transcriptomic data are needed to nail down microglia-specific relationships between variants and gene expression. This requires large numbers of human microglial samples. These are in short supply, hence researchers have used transcriptomic data from other myeloid cells, or they have made do with live microglia extracted from small numbers of people during brain surgery to chart connections between genetic variants and gene expression in the cells (Aug 2019 newsJun 2017 news; Nov 2019 news). 

Mega MiGA? Transcriptomes of 255 live microglial samples, taken from 100 human donors at autopsy, form expression profiles from four brain regions. MiGA could serve a number of purposes, including prioritization of disease genes. [Courtesy of Lopes et al., bioRxiv, 2021.]

At AAIC, Katia Lopes of the Icahn School of Medicine in New York presented initial findings of transcriptomes garnered from one of the larger brain samples to date. Called the microglia genomic atlas (MiGA), the trove was built using 255 tissue samples taken from 100 human donors at autopsy. The donors included people with no signs of neurological disease, as well as those with AD, PD, and other neurological disorders. Live microglia were teased from brain samples collected from four brain regions: the medial frontal gyrus and superior temporal gyrus in the cortex, and the thalamus and subventricular zone. Researchers led by Robin Franklin at the University of Cambridge and Daniel Gaffney at the Wellcome Sanger Institute, Hinxton, both in the U.K., recently profiled transcriptomes of live microglial among 112 tissue samples acquired from 141 donors undergoing neurosurgery (Young et al., 2021). 

Sex did not significantly alter microglial profiles, but age did: The researchers detected nearly 1,700 genes whose expression in these cells changes with age. Genes upregulated in aging included those involved in lipid metabolism and immune responses, while those that flagged with age play a role in cell motility, polarity, and IL-6 cytokine signaling. Intriguingly, some of the genes that lessen with age overlapped with those reported to dampen microglia expression in people who have AD (Srinivasan et al., 2020). 

At AAIC, Lopes showed results from a quantitative trait loci (QTL) analysis. Essentially, the scientists looked for genetic polymorphisms that influenced either the expression or splicing of genes in microglia. The eQTLs modulate expression of genes that are dubbed eGenes; sQTLs influence splicing of sGenes. By combining microglial transcriptome data from all four brain regions sampled, the researchers identified 3,611 eGenes and 4,614 sGenes. Most of these were expressed across brain regions, but some were specific. For example, expression of the RNF40 gene was influenced by one eQTL, but only in microglia inhabiting the subventricular zone.

Did any of the eQTLs or sQTLs identified in MiGA influence genes associated with AD risk loci? Indeed, the researchers found 15 MiGA QTLs that influenced expression and/or splicing of genes associated with 10 AD risk loci, meaning that some loci had more than one QTL. When more than one gene associates with a given risk locus, researchers cannot always tell which gene is driving the risk. Lopes and colleagues leveraged the co-localization of QTLs with risk loci to help answer this question.

Indeed, for one such AD risk locus—ECHDC3—Lopes used data from MiGA to zero in on the gene influenced by the variation. The lead AD risk SNP resides between two genes in the locus: ECHDC3 and USP6NL. In MiGA, Lopes identified an eQTL within the locus that associated with increased expression of USP6NL, and was also nominally tied to increased AD risk.

Using fine mapping of the locus, she found three additional SNPs associated with higher AD risk. Considering the AD risk of SNPs and eQTL SNP combined, four of the five landed within a microglial enhancer. Finally, using published microglial proximity ligation assisted ChIP-Seq (PLAC-Seq) data, which essentially catches an enhancer in the act of hooking up with a gene, Lopes found that the enhancer associated with USP6NL, but not ECHDC3 (Nott et al., 2019). 

In short, three different forms of data converged on USP6NL as the causal gene. “Now that we have all this information—eQTL, co-localization with risk SNPs, and PLAC-Seq—it is exciting to overlay all of these to find the mechanism behind the altered expression,” Lopes said.

MiGA demonstrated that many neurological disease susceptibility loci are mediated through gene expression or splicing in microglia, and suggested causal variants, noted Valentina Escott-Price of Cardiff University. The atlas may therefore help further refine gene and SNP lists to make the polygenic risk scores—which rely on estimations of causal genes—more relevant to the disease and improve prediction for eventual use in people.—Jessica Shugart

Comments

  1. The Polygenic Risk Score (PRS), capturing AD risk across the whole genome, does not work uniformly well in all datasets. For example, we have shown in our earlier work that SNPs with AD association at p<0.5, and recently at 0.1, improve AD risk prediction over and above the genome-wide significant/suggestive SNPs (Escott-Price et al., 2017Leonenko et al., 2021). Others have shown best p-value thresholds for SNP selection to be p<5x10-8, or p<1x10-5 (de Rojas et al., 2021; Zang et al., 2020). These differences could be explained by the specifics of the datasets, including age of the participants, clinical or pathological assessments of cases and/or controls, etc.

    There is a need in the field to understand these differences. Dissection of the AD risk by genes involved in different mechanisms of disease development and progression could help. One way forward is to look at the biological pathways suggested by GWAS studies (e.g., Kunkle et al., 2019). However, the definition of disease pathways is based mostly on the different functional categories defined by, for example, gene ontology. This can be noisy because there is little expert scrutiny, inclusion thresholds are low, and almost all AD genes are implicated in more than one pathway (Koopmans et al., 2019). Our analyses also show low predictability by PRS limited to these gene sets (Bellou et al., 2020). 

    At AAIC there were a number of presentations, including by Hyun-Sik Yang and Lianne Reus, suggesting that refining the PRS by genes expressed in relevant tissues, for example microglia or neuron-specific genes, may lead to improved prediction of a particular aspect of the disease. The work of Katia Lopes introducing the microglia genomic atlas, showing that many neurological disease susceptibility loci are mediated through gene expression or splicing in microglia, and suggesting causal variants, may help to further refine gene and SNP lists to make the PRS even more relevant to the disease and therefore improve prediction.

    References:

    . Polygenic risk score analysis of pathologically confirmed Alzheimer disease. Ann Neurol. 2017 Aug;82(2):311-314. Epub 2017 Aug 9 PubMed.

    . Identifying individuals with high risk of Alzheimer's disease using polygenic risk scores. Nat Commun. 2021 Jul 23;12(1):4506. PubMed.

    . Common variants in Alzheimer's disease and risk stratification by polygenic risk scores. Nat Commun. 2021 Jun 7;12(1):3417. PubMed. Correction.

    . Risk prediction of late-onset Alzheimer's disease implies an oligogenic architecture. Nat Commun. 2020 Sep 23;11(1):4799. PubMed.

    . Genetic meta-analysis of diagnosed Alzheimer's disease identifies new risk loci and implicates Aβ, tau, immunity and lipid processing. Nat Genet. 2019 Mar;51(3):414-430. Epub 2019 Feb 28 PubMed.

    . Age-dependent effect of APOE and polygenic component on Alzheimer's disease. Neurobiol Aging. 2020 Sep;93:69-77. Epub 2020 Apr 30 PubMed.

    . SynGO: An Evidence-Based, Expert-Curated Knowledge Base for the Synapse. Neuron. 2019 Jul 17;103(2):217-234.e4. Epub 2019 Jun 3 PubMed.

Make a Comment

To make a comment you must login or register.

References

News Citations

  1. Massive GWAS Meta-Analysis Digs Up Trove of Alzheimer’s Genes
  2. Are Early Harbingers of Alzheimer’s Scattered Across the Genome?
  3. Multi-Gene Score Predicts Cognitive Decline Independently of Brain Imaging
  4. AD Genetic Risk Tied to Changes in Microglial Gene Expression
  5. What Makes a Microglia? Tales from the Transcriptome
  6. Cell-Specific Enhancer Atlas Centers AD Risk in Microglia. Again.

Paper Citations

  1. . Common polygenic variation enhances risk prediction for Alzheimer's disease. Brain. 2015 Dec;138(Pt 12):3673-84. Epub 2015 Oct 21 PubMed.
  2. . Polygenic prediction via Bayesian regression and continuous shrinkage priors. Nat Commun. 2019 Apr 16;10(1):1776. PubMed.
  3. . Single-cell transcriptomic analysis of Alzheimer's disease. Nature. 2019 Jun;570(7761):332-337. Epub 2019 May 1 PubMed.
  4. . Cerebrospinal fluid biomarker signature in Alzheimer's disease neuroimaging initiative subjects. Ann Neurol. 2009 Apr;65(4):403-13. PubMed.
  5. . The EMIF-AD Multimodal Biomarker Discovery study: design, methods and cohort characteristics. Alzheimers Res Ther. 2018 Jul 6;10(1):64. PubMed.
  6. . Common variants in Alzheimer's disease and risk stratification by polygenic risk scores. Nat Commun. 2021 Jun 7;12(1):3417. PubMed. Correction.
  7. . An RNA-sequencing transcriptome and splicing database of glia, neurons, and vascular cells of the cerebral cortex. J Neurosci. 2014 Sep 3;34(36):11929-47. PubMed.
  8. . A map of transcriptional heterogeneity and regulatory variation in human microglia. Nat Genet. 2021 Jun;53(6):861-868. Epub 2021 Jun 3 PubMed.
  9. . Alzheimer's Patient Microglia Exhibit Enhanced Aging and Unique Transcriptional Activation. Cell Rep. 2020 Jun 30;31(13):107843. PubMed.
  10. . Brain cell type-specific enhancer-promoter interactome maps and disease-risk association. Science. 2019 Nov 29;366(6469):1134-1139. Epub 2019 Nov 14 PubMed.

Further Reading

Primary Papers

  1. . Atlas of genetic effects in human microglia transcriptome across brain regions, aging and disease pathologies. BioRxiv, October 28, 2020 bioRxiv.