Microglial gene signature reveals loss of homeostatic microglia associated with neurodegeneration of Alzheimer’s disease
Acta Neuropathologica Communications volume 9, Article number: 1 (2021)
Microglia-mediated neuroinflammation has been implicated in the pathogenesis of Alzheimer’s disease (AD). Although microglia in aging and neurodegenerative disease model mice show a loss of homeostatic phenotype and activation of disease-associated microglia (DAM), a correlation between those phenotypes and the degree of neuronal cell loss has not been clarified. In this study, we performed RNA sequencing of microglia isolated from three representative neurodegenerative mouse models, AppNL-G-F/NL-G-F with amyloid pathology, rTg4510 with tauopathy, and SOD1G93A with motor neuron disease by magnetic activated cell sorting. In parallel, gene expression patterns of the human precuneus with early Alzheimer’s change (n = 11) and control brain (n = 14) were also analyzed by RNA sequencing. We found that a substantial reduction of homeostatic microglial genes in rTg4510 and SOD1G93A microglia, whereas DAM genes were uniformly upregulated in all mouse models. The reduction of homeostatic microglial genes was correlated with the degree of neuronal cell loss. In human precuneus with early AD pathology, reduced expression of genes related to microglia- and oligodendrocyte-specific markers was observed, although the expression of DAM genes was not upregulated. Our results implicate a loss of homeostatic microglial function in the progression of AD and other neurodegenerative diseases. Moreover, analyses of human precuneus also suggest loss of microglia and oligodendrocyte functions induced by early amyloid pathology in human.
Alzheimer’s disease (AD) is the most common neurodegenerative disease that causes dementia, neuropathologically characterized by the accumulation of amyloid β (Aβ), phosphorylated Tau, and neuronal dystrophy and loss . Neuroinflammation is defined as an inflammatory response within the central nervous system (CNS), and it is mediated by activation of the innate immune system of the brain in response to inflammatory challenges, including misfolded protein aggregates that often accumulate in lesions of neurodegenerative diseases [22, 43]. Microglia are the resident innate immune cells of the CNS, and are key players to mediate neuroinflammation, playing critical roles in the recognition and clearance of Aβ in AD [9, 15, 27]. The activation phenotype of microglia was previously classified by the expression pattern of cytokines in analogy of activated macrophages: the proinflammatory “classical” activation phenotype (M1) and the anti-inflammatory “alternative” activated phenotype (M2) . However, this simplistic view of microglial phenotypes does not adequately reflect the complex physiology of microglia . The progression of neurodegenerative disease induces the loss of microglial homeostatic molecules and functions , leading to chronically progressive neuroinflammation . In addition, recent studies demonstrated that a common disease-associated microglia (DAM) or “neurodegenerative” phenotype, defined by a small set of upregulated genes, was observed in neurodegenerative diseases including AD, amyotrophic lateral sclerosis (ALS), and frontotemporal dementia, and aging [10, 20, 24]. However, it remains unclear whether the loss of homeostatic function in microglia or the DAM phenotype is correlated with the degree of neuronal cell loss, and whether DAM is beneficial or detrimental to neurodegenerative diseases.
Several transcriptomic studies using single nucleus analysis examined the molecular aspects of neuroinflammation in the prefrontal or entorhinal cortex of AD [11, 26, 49]; however, neuroinflammatory alteration in the precuneus, which is vulnerable to Aβ deposition in preclinical AD, has yet to be examined. The precuneus is medially located in the parietal lobe of the cerebral cortex, and it is a component of the default mode network, which is implicated in episodic memory retrieval, and displays high metabolic activity during the baseline resting state . Amyloid PET studies demonstrated that the precuneus is one of the brain regions where Aβ accumulation preferentially starts in preclinical AD [30, 47]. Therefore, it is important to uncover the neuroinflammatory aspects of the precuneus at the early amyloid pathology stage to better understand microglial response in early AD.
In this study, we first performed comparative gene expression analysis of isolated microglia from the three mouse models of neurodegenerative diseases: AppNL-G-F/NL-G-F mice that display an amyloid pathology , rTg4510 mice with tauopathy , and SOD1G93A mice with motor neurodegeneration . Despite robust neuroinflammation with microglial responses in all mouse models, AppNL-G-F/NL-G-F mice do not show neuronal death, whereas rTg4510 and SOD1G93A mice show a substantial loss of neurons. We found that most homeostatic microglial genes were downregulated in rTg4510 and SOD1G93A mice, and were correlated with the degree of neuronal cell loss. In contrast, DAM genes were uniformly upregulated in all disease models, and this alteration was not correlated with neurodegeneration. Moreover, precunei of individuals with early AD pathology show downregulation of some microglial genes linked to homeostatic functions and genes linked to oligodendrocyte function.
Materials and methods
Postmortem human brain tissues
The postmortem brains from 25 individuals (non-AD = 14 and AD = 11) were obtained by autopsy with informed consent and diagnosed by a neuropathologist in the Brain Bank for Aging Research, Tokyo Metropolitan Geriatric Hospital and Institute of Gerontology. All subjects with early AD pathology did not have apparent family history of dementia, therefore, were recognized as sporadic cases. The subjects were neuropathologically grouped according to the neurofibrillary tangle staging of Braak and Braak . The tissues were dissected from precuneus. A variation in the ratio of gray and white matter of dissected precuneus tissues, which may cause a substantial effect on relative levels of oligodendrocyte-derived mRNAs, was negligible, since we confirmed there was no changes in the OLIG1/2 mRNA levels in non-AD and AD brain as internal controls. The use of human postmortem brain tissue was approved by the Ethics Committee of Research Institute of Environmental Medicine, Nagoya University and Tokyo Metropolitan Institute. Brain tissues for RNA preparation were immediately frozen using liquid nitrogen and stored at − 80 °C before use.
Heterozygous App+/NL-G-F mice (C57BL/6-App < tm3(NL-G-F)Tcs >), carrying App gene with humanized Aβ sequence (G676R, F681Y, R684H), Swedish (KM670/671NL), Beyreuther/Iberian (I716F), and Arctic (E693G) mutations, were previously established by a knock-in strategy . Homozygous AppNL-G-F/NL-G-F and wild-type mice were obtained by crossbreeding, and were maintained as inbred lines. Tg(Camk2a-tTA)1Mmay Fgf14Tg(tetO-MAPT*P301L)4510Kha/J mice were previously established . A parental mutant tau responder line, conditionally expressing the 4R0N isoform of human P301L mutant tau, in the FVB/N strain (Clea Inc., Tokyo, Japan), and a tTA activator line, under the control of CaMKII promoter, in the 129 + ter/SV strain (Clea Inc.) were crossbred to generate rTg4510 mice . tTA activator line was used as control for rTg4510 mice. Transgenic mice expressing the inherited ALS-linked human SOD1G93A gene (B6.Cg-Tg (SOD1*G93A)1Gur/J) on the C57BL/6 background were obtained from Jackson Laboratory (Bar Harbor, ME, USA) . Genotyping of mice was performed as previously described [18, 23, 37].
All mice were maintained under a standard specific pathogen-free environment (12 h light–dark-cycle; 23 ± 1 °C; 50 ± 5% humidity) with free access to food and water throughout experiments. Animals were treated in compliance with the guidelines established by the Institutional Animal Care and Use Committee of Nagoya University and National Institutes for Quantum and Radiological Science and Technology.
Microglia isolation from brain and spinal cord
Magnetic-activated cell sorting (MACS) of microglia in brain or spinal cord is performed as described elsewhere . In brief, the cerebral cortex or spinal cord, dissected from mice transcardially perfused with phosphate-buffered saline (PBS), was dissociated at 37 °C for 15 min using the Neural Tissue Dissociation Kit-Postnatal Neurons (Miltenyi Biotec, Bergisch-Gladbach, Germany) by the gentle MACS Dissociator (Miltenyi Biotec). For isolation of microglia, myelin debris was removed by using Myelin Removal Beads II (Miltenyi Biotec). Purified cells were incubated with anti-CD16/CD32 antibodies (Thermo Fisher Scientific, Waltham, MA, USA) for blocking Fc receptors, and then incubated with anti-CD11b microBeads (Miltenyi Biotec) for isolating microglia. CD11b-positive microglia were isolated by magnetic cell sorting through an LS column (Miltenyi Biotec).
For the RNA-seq of the mouse samples, total RNA was extracted from MACS-isolated microglia of each models using an RNeasy Mini Kit (Qiagen, Hilden, Germany). The RNAs were sampled from cerebral cortices of 8-month-old AppNL-G-F/NL-G-F and 7-month-old rTg4510 mice and lumbar spinal cords of 5-month-old SODG93A mice together with the corresponding wild-type or control mice, respectively. For the RNA-seq of the human samples, total RNA was prepared from precuneus of frozen postmortem brain using mirVana™ miRNA Isolation Kit (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer instructions. The total RNA was qualified by using Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Libraries were prepared by using TruSeq mRNA or TruSeq Stranded mRNA (Illumina, San Diego, CA, USA), and, from these libraries, 151-nt paired-end reads were sequenced on the HiSeq X Ten with the HiSeq X Reagent Kits (Illumina) and the NovaSeq 6000 with the NovaSeq Reagent Kits (Illumina).
In silico analysis of the RNA-seq data
The reference genome assemblies and gene annotations, the mouse genome version UCSC mm10 and the human genome version UCSC hg38, were retrieved from iGenomes (https://support.illumina.com/sequencing/sequencing_software/igenome.html). The Illumina adapter sequences, as well as low quality bases (quality score < 20), were trimmed from 3′-ends of the sequencing reads with the Trim Galore v0.5.0 (https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/). The reads were qualified with FastQC v0.11.8 (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) before and after trimming. The processed reads were mapped on the reference genome assembly with HISAT2 v2.1.0 . Mapping reads prepared with the stranded kit was processed by employing strand specificity information. The alignments on rDNA regions were removed with BEDTools v2.25.0 . From mapping data, expressions of the individual annotated genes were quantified with StringTie v1.3.5 . Heatmaps of the gene expression profiles of target gene sets were created with the pheatmap package implemented in R (https://cran.r-project.org/web/packages/pheatmap/).
Differential expression analysis was performed using the Gene Spring Ver. 14.9.1 software from Agilent Technologies (Santa Clara, CA, USA) according to the manufacture’s protocol. Relative similarities of gene expression profiles of the samples were represented by a plot of the principal component analysis (PCA). The PCA analysis was performed with the prcomp function in R by using CPM values of the individual genes.
Quantification of mRNA levels by real-time PCR
Total RNA was extracted from MACS-sorted microglia using the RNeasy Micro Kit (Qiagen) according to the manufacturer’s instructions. Complementary DNA (cDNA) from MACS sorted cells was generated and amplified from 2.5 or 5 ng of total RNA by using the PrimeScript™ RT reagent Kit (Perfect Real Time) (TaKaRa Bio, Kusatsu, Japan) and 1/50 of the yield was amplified with the SYBR Premix Ex Taq II (Tli RNaseH Plus) (TaKaRa Bio) using the Thermal Cycler Dice Real Time System II or III (TaKaRa Bio). The thermocycle protocol was as follows: 1 cycle at 95 °C for 30 s, 40 cycles at 95 °C for 5 s and 60 °C for 30 s, and a dissociation stage of 95 °C for 15 s, 60 °C for 30 s, and 95 °C for 15 s. Actin was used for normalization. The primers used for real-time RT-PCR were as follows: human ACTB primer, forward, 5′-CTGGAACGGTGAAGGTGACA-3′; reverse, 5′-CGGCCACATTGTGAACTTTG-3′; human HBEGF primer, forward, 5′-GGACCCATGTCTTCGGAAAT-3′; reverse, 5′-CCCATGACACCTCTCTCCAT-3′; mouse Actb primer, forward, 5′-CGGACTCATCGTACTCCTGCTT-3′; reverse, 5′-TTGGCCTCACTGTCCACCTT-3′; mouse P2ry12 primer, forward, 5′-CATTGACCGCTACCTGAAGACC-3′; reverse, 5′-GCCTCCTGTTGGTGAGAATCATG-3′; mouse Sall1 primer, forward, 5′-TGTCAAGTTCCCAGAAATGTTCCA-3′; reverse, 5′-ATGCCGCCGTTCTGAATGA-3′; mouse Apoe primer, forward, 5′-GAACCGCTTCTGGGATTACCTG-3′; reverse, 5′- GCCTTTACTTCCGTCATAGTGTC-3′; mouse Itgax (Cd11c) primer, forward, 5′-CTGGATAGCCTTTCTTCTGCTG-3′; reverse, 5′-GCACACTGTGTCCGAACTCA-3′; mouse Hbegf primer, forward, 5′-ACCAGTGGAGAATCCCCTATAC-3′; reverse, 5′-GCCAAGACTGTAGTGTGGTCA-3′.
Immunofluorescence analysis was performed as described previously . In brief, mice were deeply anesthetized and perfused intracardially with PBS and 4% paraformaldehyde in PBS. The brains or spinal cords were dissected, post-fixed with the same fixative, and cryoprotected with 30% sucrose containing PBS. Twenty-micrometer-thick coronal brain sections or spinal cord transverse sections were fixed with 4% paraformaldehyde in PBS for 5 min and then permeabilized with 0.1% Triton X−100/PBS for 10 min. After an incubation in blocking solution (5% goat or donkey serum/PBS) for 1 h, sections were incubated with a combination of following antibodies: rabbit anti-Iba-1 (#019–19741, 1:500; Wako, Osaka, Japan), goat anti-AIF-1/Iba1 (#NB100-1028, 1:250; Novus Biologicals, CO, USA), mouse anti-human SOD1 (#M062-3, 1:500; MBL, Nagoya, Japan), mouse anti-AT8 (#MN1020, 1:500, Invitrogen, CA, USA), rabbit anti-Tmem119 (#ab209064, 1:500; abcam, Cambridge, UK), or mouse anti-ApoE (#sc-390925, 1:100, Santa Cruz Biotechnology, Inc, CA, USA) at 4 °C overnight. After washing with PBS, sections were incubated with fluorescent-conjugated anti-rabbit, anti-mouse, or anti-goat IgGs (1:1000; Thermo Fisher Scientific) at room temperature for 1 h. For thioflavin-S staining, after immunostaining with primary and secondary antibodies, sections were incubated with 0.02% thioflavin-S (#T1892, Sigma) at room temperature for 8 min, followed by incubation with 50% ethanol at room temperature for 1 min twice. After washing in PBS, sections were mounted on slides and analyzed with a confocal microscope (LSM700, Carl Zeiss, Oberkochen, Germany).
For RNA sequencing data, differential expression analysis between two groups with replicates was performed with edgeR v3.24.3  implemented in R v3.5.1, and q-value was calculated for multiple testing correction of the p-values with the q-value package v2.24.1  in R. TPM (Transcripts Per Million) and CPM (Counts Per Million) values were computed with the StringTie and edgeR, respectively. A statistical comparison of the fold change values of gene expression quantities between two mutant experiments were performed with Wilcoxon’s signed rank test, followed by a multiple testing correction with the Bonferroni–Holm method . The representative log fold change values of the experiments for this purpose were computed with the edgeR. For quantitative PCR, differences between two groups were analyzed by a two-tailed Student’s t test with the Bonferroni method.
Altered gene expression profiles of isolated microglia in mouse models of neurodegenerative diseases
To compare microglial phenotypes in different neurodegenerative diseases, we analyzed three representative model mice: AppNL-G-F/NL-G-F mice with amyloid pathology , rTg4510 mice with tauopathy , and SOD1G93A mice modeling inherited ALS  (Fig. 1a). All mouse models exhibited accumulation of disease-causing proteins at 2.5 ~ 3 months of age followed by neuroinflammation/gliosis (Fig. 1b). Despite age-dependent typical amyloid pathology and neuroinflammation, AppNL-G-F/NL-G-F mice did not show neuronal loss or neurofibrillary tangles (NFTs) . In contrast, rTg4510 and SOD1G93A mice exhibited robust neurodegeneration. Therefore, we used the AppNL-G-F/NL-G-F mouse as a mild model and rTg4510 and SOD1G93A mice as severe models of neurodegenerative diseases. We analyzed gene expression in microglia isolated by magnetic-activated cell sorting (MACS) from cerebral cortex of AppNL-G-F/NL-G-F and rTg4510 mice or lumbar spinal cord of SOD1G93A mice at the middle-to-late disease stage (Fig. 1a, b).
We aimed to characterize gene expression profiles in microglia isolated from each mouse model via RNA sequencing (RNA-seq) (Fig. 1a). Analysis of expression profiles revealed 3318, 2631, and 3106 differentially expressed genes in AppNL-G-F/NL-G-F, rTg4510, and SOD1G93A microglia, respectively (Fig. 1c, d, Additional file 1: Tables S1a–c). We then performed principal component analysis (PCA) by employing normalized read counts of the genes for individual samples. Based on PCA plot, we found that the first principal component (PC1), accounting for 88.8% of the total variance, solely separated the samples into four categories comprised of the three neurogenerative and one normal phenotype (Fig. 1e). Furthermore, the four clusters were located in order of the degree of neuronal cell loss along the PC1, suggesting that gene expression profiles gradually change in accordance with the degree of neuronal cell loss.
Decrease in expression of homeostatic microglial genes in mouse models of tauopathy and SOD1-ALS
We first focused on the homeostatic microglial genes, as diseased microglia lose homeostatic functions in an apolipoprotein E (APOE)- and TREM2-dependent manner . To clarify whether the loss of homeostatic microglial genes correlates with the degree of neuronal cell loss, we investigated expression levels of the 68 homeostatic genes , expressions of which are all enriched in mature microglia, among the three models. As mentioned earlier, AppNL-G-F/NL-G-F mice show neuroinflammation and an amyloid pathology without neuronal cell death, whereas rTg4510 and SOD1G93A mice show robust neuroinflammation with neuronal cell death. A heat map of the 68 gene expression profiles illustrates changes in expression of homeostatic microglial genes associated with the degree of neuronal cell loss (Fig. 2a). P2ry12 and Sall1 exhibited a typical trend of the changes of gene expression level among the three models; a non-significant change in AppNL-G-F/NL-G-F microglia, but substantial reductions in rTg4510 and SOD1G93A microglia (Fig. 2a–c). Moreover, the degree of reduction in expression levels of homeostatic microglial genes significantly correlated with the degree of neuronal cell loss (Fig. 2d and Additional file 1: Table S1d). Furthermore, although a loss of TMEM119 immunoreactivity, another representative homeostatic microglial gene, was restricted to amyloid plaque-associated microglia in AppNL-G-F/NL-G-F mice (Fig. 2e), TMEM119 immunoreactivity was decreased in microglia of rTg4510 and SOD1G93A mice (Fig. 2f, g). These results indicate that decreased expression of microglial homeostatic genes correlates with the degree of neuronal cell loss.
No association between changes in DAM genes and the degree of neuronal cell loss
To examine whether DAM are common in mouse models of neurodegenerative diseases, and whether this activation is also associated with the degree of neuronal cell loss, we compared expression levels of 162 DAM genes among three mouse models (FC > 1.5, q < 0.05)  (Fig. 3a). A heat map shows that almost all DAM genes including Itgax (Fig. 3c) were uniformly upregulated in microglia of the three mouse models, while only 12 of these genes (7.41%), including Apoe, Axl, and Cybb, exhibited upregulation in accord with the degree of neuronal cell loss (Fig. 3a–d and Additional file 1: Table S1e). Thus, the degree of increased expression of DAM genes was not associated with the degree of neuronal cell loss (Fig. 3d and Additional file 1: Table S1e). Although upregulation of ApoE was limited in plaque-associated microglia of AppNL-G-F/NL-G-F mice, immunoreactivity of ApoE was increased in microglia of all three models (Fig. 3e–g). These results suggest that neurodegenerative diseases have upregulation of most DAM genes in common, but levels of expression were not correlated with the degree of neuronal cell loss.
Differentially expressed genes for various CNS cell types in human precuneus with AD pathology
To uncover differentially expressed genes in human brain with an early amyloid pathology, we created a gene expression profile of precuneus derived from the individuals neuropathologically diagnosed with early AD and controls. All subjects with early AD did not have apparent family history of dementia, therefore, were recognized as sporadic cases. The precuneus is located on the medial side of the parietal cortex, which is vulnerable to early amyloid deposition in patients with AD [30, 47]. Data for Braak neuropathological stages, neuropathological diagnosis, Clinical Dementia Ratings (CDR), APOE genotypes, and RNA quality of human brain samples are described in Additional file 1: Table S1f. Postmortem subjects were identified and selected based on the Braak neuropathological staging as follows: 14 non-AD brains were scored 0–A for senile plaque (SP) and 0–II for NFT, and 11 AD brains were C for SP and III–IV for NFT (Fig. 4a and Additional file 1: Table S1f). We found 643 deregulated genes, consisting of 127 upregulated and 516 downregulated, in AD brains compared with the controls (|FC| > 1.5, q < 0.05) (Additional file 1: Table S1g).
In AD and other neurodegenerative diseases, neuroinflammation and neurodegeneration lead to changes in cell population of brain lesions. Therefore, we examined expression levels of representative CNS markers for each cell-type: neurons, astrocytes, microglia, and oligodendrocytes in the precuneus of AD brains and controls (Fig. 4b, Additional file 1: Table S1h). We generated a list of genes enriched in microglia under physiological conditions from previous studies [3, 5, 8] and compared gene expression profiles of AD precuneus. Although levels of the representative activation markers for microglia (AIF1, CD68, and LGALS3) were unchanged, expressions of BCO2, CAPN3, small G protein-associated genes (RASGRP3 and RAPGEF5), PACC1 (also known as TMEM206), and P2RX7 were significantly decreased in AD precuneus compared with controls (Fig. 4c). Moreover, expression of the oligodendrocyte markers, MBP, MAG, CLDN11, MOG, and CNP were significantly decreased in AD precuneus (Fig. 4c). No differences in the expression of representative neuronal and astrocytic genes were observed in the precuneus of AD brains (Fig. 4c, Additional file 1: Table S1h). These findings indicate that the early amyloid pathology induces moderate dysfunction of microglia and oligodendrocytes in AD precuneus.
Subtle alteration of DAM gene expression in AD precuneus
Reactive microgliosis with neuroinflammation is one of the neuropathological hallmarks of AD brain [9, 15, 40, 43]. Thus, we next investigated whether DAM genes are deregulated in human precuneus with early AD pathology (Fig. 5a). Among the 162 DAM genes, we found only 8 genes were deregulated in AD precuneus, and unexpectedly, all of them (APBB2, ARAP2, DHCR7, ENPP2, MYO1E, CDS22, KCNJ2, and SLC44A1) were downregulated (Fig. 5b). In addition, expression levels of representative DAM genes (ITGAX, CST7, APOE, CSF1, and AXL) were not altered in the precuneus of AD brain (Fig. 5b). Contrary to the upregulation of DAM genes in mouse models of AD, the expression of DAM genes was marginally altered in human sporadic AD precuneus.
Expression levels of the genes defined as AD risk factors in microglia of mouse models of AD and human sporadic AD precuneus
We evaluated expression levels of 54 genes [5, 13], which were previously defined as risk factors for AD by genome-wide association study, in human sporadic AD brains and isolated microglia from AppNL-G-F/NL-G-F and rTg4510 mice (Fig. 6a and Additional file 1: Table S1i). The previous study identified GRIN2B, INPP5D and PSEN1 as significantly altered genes in the temporal or frontal cortex of AD brain . We focused on those risk genes, first confirming their expressions in microglia by the transcriptome database , and then examined their expression levels in human AD precuneus and AD model mice (cut-off TPM > 5, q < 0.05; Fig. 6b). We found no upregulated genes, but the expression level of HBEGF (log2FC = −0.843, q = 0.00861) was significantly decreased in human AD precuneus (Fig. 6b). Further, of 54 risk genes with significantly altered expression, we found 19 and 14 deregulated genes in cortical microglia of AppNL-G-F/NL-G-F and rTg4510 mice, respectively. Comparison of AD risk genes among the three groups (human AD, AppNL-G-F/NL-G-F and rTg4510) revealed that 12 genes (Apoe, Scimp, Pld3, Psen1, Psen2, Sqstm1, Trem2, Treml2, Ap2a2, Cass4, Lmo4 and Hbegf) were commonly altered in AppNL-G-F/NL-G-F and rTg4510 cortical microglia, and only one gene (HBEGF/Hbegf) was significantly decreased in AD precuneus as well as AppNL-G-F/NL-G-F and rTg4510 cortical microglia (Fig. 6b and Additional file 1: Table S1i). In addition, HBEGF/Hbegf expression levels in human precuneus and mouse microglia of either models were validated by quantitative PCR (Fig. 6c, d), and were consistent with the RNA-seq data (Fig. 6a, b and Additional file 1: Table S1i). These results indicate that expression changes of AD risk genes were robust in microglia of the mouse models of AD, but minor in the precuneus of individuals with early AD pathology.
Profiling of commonly altered genes in human sporadic AD precuneus and microglia of the mouse models of AD
To understand the neuroinflammatory nature of early amyloid pathology in AD, we generated expression profiles of commonly altered genes in the precuneus of human sporadic AD brain and cortical microglia isolated from AppNL-G-F/NL-G-F and/or rTg4510 mice (Fig. 6e). We compared gene profiles of AD precunei (|FC|> 1.2, q < 0.05) and deregulated genes of cortical microglia of AppNL-G-F/NL-G-F or rTg4510 mice (|FC| > 1.5, q < 0.05, cut-off TPM > 5). We found that deregulated 32 genes, consisting of 8 upregulated and 24 downregulated, were common among AD precuneus, AppNL-G-F/NL-G-F and rTg4510 (Fig. 6f, g, Additional file 1: Tables S1j and k). The eight upregulated genes include a chemokine (CXCL10), and interferon-induced genes (STAT1, IFIT3, ISG15), indicating neuroinflammatory changes in AD precuneus and AD mice (Table S1j). In human precuneus of AD brains and AppNL-G-F/NL-G-F mice, we found 44 deregulated genes (3 upregulated and 41 downregulated), including TNFSF10 and SLC6A8, whereas in AD precuneus and rTg4510 mice, we found 13 deregulated genes (2 upregulated and 11 downregulated) including PTPRD. These results suggest that some deregulated genes are common between microglia of different AD mouse models and human sporadic AD precuneus, and that AppNL-G-F/NL-G-F and rTg4510 cortical microglia may represent different neuroinflammatory aspects relevant to AD pathologies.
Our gene expression analysis of isolated microglia from mice of three neurodegenerative disease models including AD revealed that downregulation of homeostatic microglial genes and robust upregulation of DAM genes. More importantly, the extent of downregulation of homeostatic microglial genes was correlated with the degree of neuronal cell loss. In addition, a gene expression profile of human precuneus of early AD pathology shows reduced expression of some microglial homeostatic genes without induction of DAM genes, suggesting that the microglial signature in human early AD has limited similarity to that found in AD mice.
Previous studies reported a loss of homeostatic microglial genes in neurodegenerative diseases such as AD, ALS, and multiple sclerosis [20, 24, 50]. In contrast to their findings, we found that rTg4510 and SOD1G93A mice exhibited more robust downregulation in homeostatic microglial genes than that of AppNL-G-F/NL-G-F mice, which exhibit neuroinflammation without neuronal loss (Fig. 2). In addition, rTg4510 mice show brain atrophy because of neurodegeneration, and SOD1G93A mice exhibit fatal motor neurodegeneration. Of particular interest, changes in homeostatic microglial genes include P2ry12, Sall1, and Tmem119 (Fig. 2a–d). P2RY12 regulates microglial activation via extracellular nucleotides [3, 14]. SALL1 inhibits a reactive microglia phenotype and promotes a physiological surveilling phenotype , and TMEM119 is expressed in microglia-derived cells and is absent in macrophages [3, 46]. Although a loss of Tmem119 was limited to Aβ plaque-associated microglia in AppNL-G-F/NL-G-F mice, more robust decreases in Tmem119 were observed in rTg4510 and SOD1G93A microglia. Although there is a possibility that a loss of homeostatic microglial genes might be attributed to overexpression of transgenes, our result indicates that a loss of unique microglial homeostatic genes correlates with severity of neurodegeneration, and that decreased homeostatic microglial markers may be one of the hallmarks of progressive neuronal loss.
Consistent with the previous studies [3, 20, 24], we found that DAM genes were uniformly upregulated in all three neurodegenerative mouse models. However, we also found that upregulation of most DAM genes was not correlated with the degree of neuronal cell loss (Fig. 3a, c, d). As a rare exception, Apoe was correlated with the degree of neuronal cell loss (Fig. 3b), and the change in Apoe expression was negatively associated with altered expression of homeostatic microglial genes (Fig. 2). The ApoE4 triggers induction of DAM with impaired homeostatic functions and is responsible for driving neurodegeneration [24, 39]. This finding indicates that ApoE signaling may be involved in the loss of homeostatic microglial function. Combined with previous studies, our findings suggest that most DAM genes do not directly accelerate progression of neurodegeneration, but a few DAM genes, such as APOE may be involved in severity of neurodegeneration. Since the relationship between APOE genotypes and expressions of DAM and homeostatic microglial genes cannot be determined due to the limitation of APOE e4 carriers in our study, further studies are required to clarify the potential role of APOE e4 in human microglial phenotypes.
Despite its importance as a vulnerable region of amyloid deposition in early AD, gene expression in the precuneus of patients with AD has not been well explored in the past. Our study provides the first comprehensive gene expression profile of precuneus at the early stage of AD pathology (Braak NFT stage III–IV). Although the slight elevation of astrocytic and microglia activation markers in early AD precuneus was not statistically significant, the upregulation of chemokine and proinflammatory genes (CXCL10, STAT1, ISG15, and ITIF3) indicates neuroinflammatory changes in early AD precuneus. We also found downregulation of several microglial markers and DAM genes in early AD precuneus (Figs. 4, 5). In particular, the unaltered state of most DAM genes in human AD precuneus was in striking contrast with the data from AD mouse models. One possible interpretation is that the expression of genes of human AD and mouse models of AD is discordant, as demonstrated by a study using 5XFAD mice . Another possibility is that microglial function may be suppressed at the early stage of amyloid pathology. A recent study pointed out the low sensitivity of single nucleus RNA-seq to detect DAM genes in human postmortem brains . This low sensitivity may be attributed to redistribution of DAM mRNAs to the cytosol or instability of DAM mRNAs. As our study is based on the RNA-seq of whole brain tissues, cytosolic mRNAs can be detected, therefore, the latter possibility namely, instability of DAM mRNAs in human postmortem brain, remains. Further comparative gene expression analyses using the brain samples with early and advanced neuropathology will enable a conclusion regarding the significance of homeostatic microglia and DAM in human AD.
Another striking finding was the robust downregulation of genes linked to oligodendrocytes (MBP, MAG, CLDN11, MOG, and CNP) in early AD precuneus. In recent years, increasing evidence suggests that white matter abnormalities are also an important component of AD , and oligodendrocyte abnormalities in advanced AD brain were also reported by a recent transcriptomic study . Compared with the microglia-mediated neuroinflammation, the role of oligodendrocytes in AD has been underestimated. Our data of the deregulation of oligodendrocytes and microglia in early AD precuneus suggest that both microglia and oligodendrocyte dysfunction may play a role in development of early AD pathology.
Although many risk genes for AD have been linked to microglia functions , we found downregulation of HBEGF/Hbegf was the only common change among the AD risk genes in human precuneus of early AD and AD model mice. A previous study reported an association between the single nucleotide polymorphism rs77493189 in HBEGF and late-onset AD . Moreover, forebrain-specific Hbegf knockout mice exhibited neurotrophic and growth factor imbalances as well as impaired memory function and synaptic plasticity . In this study, we performed comparative analysis of the gene expression profiles of sporadic AD precuneus and microglia isolated from familial AD models. Therefore, a small number of risk genes commonly deregulated in human precuneus and mouse microglia may be attributed to the different etiology between sporadic and familial AD. Despite of the limitation in our study, our results suggest that supplementation of HBEGF may be a viable therapeutic target for AD.
Comparative analysis of the gene expression profiles of AD precuneus and microglia isolated from AppNL-G-F/NL-G-F and rTg4510 mouse cortices revealed that the number of commonly altered genes with AD precuneus were greater in AppNL-G-F/NL-G-F microglia (76 genes) compared with rTg4510 microglia (45 genes). AppNL-G-F/NL-G-F mice specifically shared 44 deregulated genes with AD precuneus, and TNFSF10, one of the upregulated genes, is implicated in neuroinflammation and Aβ accumulation . In contrast, SLC6A8 is one of the 41 downregulated genes, and it mediates inflammation by regulating creatine uptake . We also found that PTPRD was decreased in both AD precuneus and rTg4510 mice; an association was reported between the rs560380 polymorphism in PTPRD and NFT burden . These results suggest that AppNL-G-F/NL-G-F and rTg4510 cortical microglia represent different aspects of AD pathologies.
Microglial gene signature in mice revealed that a loss of homeostatic microglia function is associated with the degree of neuronal cell loss. In humans, our evaluation of the precuneus of early AD pathology also suggests a loss of microglia and oligodendrocyte function induced by early amyloid pathology. Results from the present study indicate a correlation between glial phenotypes and severity of neurodegeneration, and also provide important resources to better understand the role of glial dysfunction in AD progression.
amyotrophic lateral sclerosis
adaptor related protein complex 2 subunit alpha 2
amyloid precursor protein
AXL receptor tyrosine kinase
cas scaffold protein family member 4
counts per million
C-X-C motif chemokine ligand 10
cytochrome b-245 beta chain
disease associated microglia
heparin binding EGF like growth factor
interferon induced protein with tetratricopeptide repeats 3
ISG15 ubiquitin-like modifier
integrin subunit alpha X
LIM domain only 4
positron emission tomography
purinergic receptor P2Y12
phospholipase D family member 3
spalt like transcription factor 1
SLP adaptor and CSK interacting membrane protein
solute carrier family 6 member 8
signal transducer and activator of transcription 1
transmembrane protein 119
TNF superfamily member 10
transcripts per million
triggering receptor expressed on myeloid cells 2
triggering receptor expressed on myeloid cells like 2
Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B 57:289–300. https://doi.org/10.2307/2346101
Braak HBE (1991) Neuropathological stageing of Alzheimer-related changes. Acta Neuropathol 82:239–259. https://doi.org/10.1007/bf00308809
Butovsky O, Jedrychowski MP, Moore CS, Cialic R, Lanser AJ, Gabriely G, Koeglsperger T, Dake B, Wu PM, Doykan CE et al (2014) Identification of a unique TGF-beta-dependent molecular and functional signature in microglia. Nat Neurosci 17:131–143. https://doi.org/10.1038/nn.3599
Cantarella G, Di Benedetto G, Puzzo D, Privitera L, Loreto C, Saccone S, Giunta S, Palmeri A, Bernardini R (2015) Neutralization of TNFSF10 ameliorates functional outcome in a murine model of Alzheimer’s disease. Brain 138:203–216. https://doi.org/10.1093/brain/awu318
Castillo E, Leon J, Mazzei G, Abolhassani N, Haruyama N, Saito T, Saido T, Hokama M, Iwaki T, Ohara T et al (2017) Comparative profiling of cortical gene expression in Alzheimer’s disease patients and mouse models demonstrates a link between amyloidosis and neuroinflammation. Sci Rep 7:17762. https://doi.org/10.1038/s41598-017-17999-3
Cavanna AE, Trimble MR (2006) The precuneus: a review of its functional anatomy and behavioural correlates. Brain 129:564–583. https://doi.org/10.1093/brain/awl004
Chibnik LB, White CC, Mukherjee S, Raj T, Yu L, Larson EB, Montine TJ, Keene CD, Sonnen J, Schneider JA et al (2018) Susceptibility to neurofibrillary tangles: role of the PTPRD locus and limited pleiotropy with other neuropathologies. Mol Psychiatry 23:1521–1529. https://doi.org/10.1038/mp.2017.20
Chiu IM, Morimoto ET, Goodarzi H, Liao JT, O’Keeffe S, Phatnani HP, Muratet M, Carroll MC, Levy S, Tavazoie S et al (2013) A neurodegeneration-specific gene-expression signature of acutely isolated microglia from an amyotrophic lateral sclerosis mouse model. Cell Rep 4:385–401. https://doi.org/10.1016/j.celrep.2013.06.018
De Strooper B, Karran E (2016) The cellular phase of Alzheimer’s disease. Cell 164:603–615. https://doi.org/10.1016/j.cell.2015.12.056
Friedman BA, Srinivasan K, Ayalon G, Meilandt WJ, Lin H, Huntley MA, Cao Y, Lee SH, Haddick PCG, Ngu H et al (2018) Diverse brain myeloid expression profiles reveal distinct microglial activation states and aspects of Alzheimer’s disease not evident in mouse models. Cell Rep 22:832–847. https://doi.org/10.1016/j.celrep.2017.12.066
Grubman A, Chew G, Ouyang JF, Sun G, Choo XY, McLean C, Simmons RK, Buckberry S, Vargas-Landin DB, Poppe D et al (2019) A single-cell atlas of entorhinal cortex from individuals with Alzheimer’s disease reveals cell-type-specific gene expression regulation. Nat Neurosci 22:2087–2097. https://doi.org/10.1038/s41593-019-0539-4
Gurney ME, Pu H, Chiu AY, Dal Canto MC, Polchow CY, Alexander DD, Caliendo J, Hentati A, Kwon YW, Deng HX (1994) Motor neuron degeneration in mice that express a human Cu, Zn superoxide dismutase mutation. Science 264:1772–1775. https://doi.org/10.1126/science.8209258
Hansen DV, Hanson JE, Sheng M (2018) Microglia in Alzheimer’s disease. J Cell Biol 217:459–472. https://doi.org/10.1083/jcb.201709069
Haynes SE, Hollopeter G, Yang G, Kurpius D, Dailey ME, Gan WB, Julius D (2006) The P2Y12 receptor regulates microglial activation by extracellular nucleotides. Nat Neurosci 9:1512–1519. https://doi.org/10.1038/nn1805
Heneka MT, Carson MJ, Khoury JE, Landreth GE, Brosseron F, Feinstein DL, Jacobs AH, Wyss-Coray T, Vitorica J, Ransohoff RM et al (2015) Neuroinflammation in Alzheimer’s disease. Lancet Neurol 14:388–405. https://doi.org/10.1016/s1474-4422(15)70016-5
Holtman IR, Raj DD, Miller JA, Schaafsma W, Yin Z, Brouwer N, Wes PD, Moller T, Orre M, Kamphuis W et al (2015) Induction of a common microglia gene expression signature by aging and neurodegenerative conditions: a co-expression meta-analysis. Acta Neuropathol Commun 3:31. https://doi.org/10.1186/s40478-015-0203-5
Holtman IR, Skola D, Glass CK (2017) Transcriptional control of microglia phenotypes in health and disease. J Clin Invest 127:3220–3229. https://doi.org/10.1172/JCI90604
Ishikawa A, Tokunaga M, Maeda J, Minamihisamatsu T, Shimojo M, Takuwa H, Ono M, Ni R, Hirano S, Kuwabara S et al (2018) In vivo visualization of tau accumulation, microglial activation, and brain atrophy in a mouse model of Tauopathy rTg4510. J Alzheimers Dis 61:1037–1052. https://doi.org/10.3233/JAD-170509
Ji L, Zhao X, Zhang B, Kang L, Song W, Zhao B, Xie W, Chen L, Hu X (2019) Slc6a8-mediated creatine uptake and accumulation reprogram macrophage polarization via regulating cytokine responses. Immunity 51(272–284):e277. https://doi.org/10.1016/j.immuni.2019.06.007
Keren-Shaul H, Spinrad A, Weiner A, Matcovitch-Natan O, Dvir-Szternfeld R, Ulland TK, David E, Baruch K, Lara-Astaiso D, Toth B et al (2017) A unique microglia type associated with restricting development of Alzheimer’s disease. Cell 169(1276–1290):e1217. https://doi.org/10.1016/j.cell.2017.05.018
Kim D, Paggi JM, Park C, Bennett C, Salzberg SL (2019) Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol 37:907–915. https://doi.org/10.1038/s41587-019-0201-4
Komine O, Yamanaka K (2015) Neuroinflammation in motor neuron disease. Nagoya J Med Sci 77:537–549
Komine O, Yamashita H, Fujimori-Tonou N, Koike M, Jin S, Moriwaki Y, Endo F, Watanabe S, Uematsu S, Akira S et al (2018) Innate immune adaptor TRIF deficiency accelerates disease progression of ALS mice with accumulation of aberrantly activated astrocytes. Cell Death Differ 25:2130–2146. https://doi.org/10.1038/s41418-018-0098-3
Krasemann S, Madore C, Cialic R, Baufeld C, Calcagno N, El Fatimy R, Beckers L, O’Loughlin E, Xu Y, Fanek Z et al (2017) The TREM2-APOE pathway drives the transcriptional phenotype of dysfunctional microglia in neurodegenerative diseases. Immunity 47(566–581):e569. https://doi.org/10.1016/j.immuni.2017.08.008
Liu JZ, Erlich Y, Pickrell JK (2017) Case-control association mapping by proxy using family history of disease. Nat Genet 49:325–331. https://doi.org/10.1038/ng.3766
Mathys H, Davila-Velderrain J, Peng Z, Gao F, Mohammadi S, Young JZ, Menon M, He L, Abdurrob F, Jiang X et al (2019) Single-cell transcriptomic analysis of Alzheimer’s disease. Nature 570:332–337. https://doi.org/10.1038/s41586-019-1195-2
Mosher KI, Wyss-Coray T (2014) Microglial dysfunction in brain aging and Alzheimer’s disease. Biochem Pharmacol 88:594–604. https://doi.org/10.1016/j.bcp.2014.01.008
Nasrabady SE, Rizvi B, Goldman JE, Brickman AM (2018) White matter changes in Alzheimer’s disease: a focus on myelin and oligodendrocytes. Acta Neuropathol Commun 6:22. https://doi.org/10.1186/s40478-018-0515-3
Oyagi A, Moriguchi S, Nitta A, Murata K, Oida Y, Tsuruma K, Shimazawa M, Fukunaga K, Hara H (2011) Heparin-binding EGF-like growth factor is required for synaptic plasticity and memory formation. Brain Res 1419:97–104. https://doi.org/10.1016/j.brainres.2011.09.003
Palmqvist S, Scholl M, Strandberg O, Mattsson N, Stomrud E, Zetterberg H, Blennow K, Landau S, Jagust W, Hansson O (2017) Earliest accumulation of beta-amyloid occurs within the default-mode network and concurrently affects brain connectivity. Nat Commun 8:1214. https://doi.org/10.1038/s41467-017-01150-x
Perry VHHC (2014) Microglial priming in neurodegenerative disease. Nat Rev Neurol 10:217–224. https://doi.org/10.1038/nrneurol.2014.38
Pertea M, Pertea GM, Antonescu CM, Chang TC, Mendell JT, Salzberg SL (2015) StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol 33:290–295. https://doi.org/10.1038/nbt.3122
Querfurth HW, LaFerla FM (2010) Alzheimer’s Disease. N Engl J Med 362:329–344. https://doi.org/10.1056/NEJMra0909142
Quinlan AR, Hall IM (2010) BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26:841–842. https://doi.org/10.1093/bioinformatics/btq033
Ransohoff RM (2016) A polarizing question: do M1 and M2 microglia exist? Nat Neurosci 19:987–991. https://doi.org/10.1038/nn.4338
Robinson MD, McCarthy DJ, Smyth GK (2010) edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26:139–140. https://doi.org/10.1093/bioinformatics/btp616
Saito T, Matsuba Y, Mihira N, Takano J, Nilsson P, Itohara S, Iwata N, Saido TC (2014) Single App knock-in mouse models of Alzheimer’s disease. Nat Neurosci 17:661–663. https://doi.org/10.1038/nn.3697
Santacruz K, Lewis J, Spires T, Paulson J, Kotilinek L, Ingelsson M, Guimaraes A, DeTure M, Ramsden M, McGowan E et al (2005) Tau suppression in a neurodegenerative mouse model improves memory function. Science 309:476–481. https://doi.org/10.1126/science.1113694
Shi Y, Yamada K, Liddelow SA, Smith ST, Zhao L, Luo W, Tsai RM, Spina S, Grinberg LT, Rojas JC et al (2017) ApoE4 markedly exacerbates tau-mediated neurodegeneration in a mouse model of tauopathy. Nature 549:523–527. https://doi.org/10.1038/nature24016
Skokowa J, Cario G, Uenalan M, Schambach A, Germeshausen M, Battmer K, Zeidler C, Lehmann U, Eder M, Baum C et al (2006) LEF-1 is crucial for neutrophil granulocytopoiesis and its expression is severely reduced in congenital neutropenia. Nat Med 12:1191–1197. https://doi.org/10.1038/nm1474
Sobue A, Ito N, Nagai T, Shan W, Hada K, Nakajima A, Murakami Y, Mouri A, Yamamoto Y, Nabeshima T et al (2018) Astroglial major histocompatibility complex class I following immune activation leads to behavioral and neuropathological changes. Glia 66:1034–1052. https://doi.org/10.1002/glia.23299
Storey JD, Tibshirani R (2003) Statistical significance for genomewide studies. Proc Natl Acad Sci USA 100:9440–9445. https://doi.org/10.1073/pnas.1530509100
Streit WJ, Mrak RE, Griffin WS (2004) Microglia and neuroinflammation: a pathological perspective. J Neuroinflammation 1:14. https://doi.org/10.1186/1742-2094-1-14
Tang Y, Le W (2016) Differential roles of M1 and M2 Microglia in neurodegenerative diseases. Mol Neurobiol 53:1181–1194. https://doi.org/10.1007/s12035-014-9070-5
Thrupp N, Sala Frigerio C, Wolfs L, Skene NG, Fattorelli N, Poovathingal S, Fourne Y, Matthews PM, Theys T, Mancuso R et al (2020) Single-nucleus RNA-Seq is not suitable for detection of microglial activation genes in humans. Cell Rep 32:108189. https://doi.org/10.1016/j.celrep.2020.108189
van Wageningen TA, Vlaar E, Kooij G, Jongenelen CAM, Geurts JJG, van Dam AM (2019) Regulation of microglial TMEM119 and P2RY12 immunoreactivity in multiple sclerosis white and grey matter lesions is dependent on their inflammatory environment. Acta Neuropathol Commun 7:206. https://doi.org/10.1186/s40478-019-0850-z
Wolk DA, Price JC, Saxton JA, Snitz BE, James JA, Lopez OL, Aizenstein HJ, Cohen AD, Weissfeld LA, Mathis CA et al (2009) Amyloid imaging in mild cognitive impairment subtypes. Ann Neurol 65:557–568. https://doi.org/10.1002/ana.21598
Zhang Y, Chen K, Sloan SA, Bennett ML, Scholze AR, O’Keeffe S, Phatnani HP, Guarnieri P, Caneda C, Ruderisch N et al (2014) An RNA-sequencing transcriptome and splicing database of glia, neurons, and vascular cells of the cerebral cortex. J Neurosci 34:11929–11947. https://doi.org/10.1523/JNEUROSCI.1860-14.2014
Zhou Y, Song WM, Andhey PS, Swain A, Levy T, Miller KR, Poliani PL, Cominelli M, Grover S, Gilfillan S et al (2020) Human and mouse single-nucleus transcriptomics reveal TREM2-dependent and TREM2-independent cellular responses in Alzheimer’s disease. Nat Med 26:131–142. https://doi.org/10.1038/s41591-019-0695-9
Zrzavy T, Hametner S, Wimmer I, Butovsky O, Weiner HL, Lassmann H (2017) Loss of ‘homeostatic’ microglia and patterns of their activation in active multiple sclerosis. Brain 140:1900–1913. https://doi.org/10.1093/brain/awx113
We thank Center for Animal Research and Education (CARE), Nagoya University for a technical support on animal experiments.
This work was partly supported by Grants-in-Aid for Scientific Research 18H06101, 20K16488 (to AS), and 16H06277 (CoBiA, to SM) from the Japan Society for the Promotion of Science (JSPS), the Strategic Research Program for Brain Sciences from AMED under Grant Number 20dm0107135 (KY) and JP20dm0107106 (SM), Brain Mapping by Integrated Neurotechnologies for Disease Studies from AMED under Grant Number 20dm0207072 (MH), Takeda Science Foundation (KY), and The Hori Science and Arts Foundation (KY).
Availability of data and materials
Data, material and software information supporting the conclusions of this article is included within the article and its additional files. The raw RNA-seq datasets analyzed in the current study are available from the corresponding author on reasonable request, and will be deposited to a public data repository.
Ethics approval and consent to participate
The experiments using human brains were approved by the Ethics Committees at Nagoya University (Approval No. #328). The experiments using genetically modified mice were approved by the Animal Care and Use Committee and the recombinant DNA experiment committee of Nagoya University (approval numbers RIEM19269 and RIEM19275, and #143, respectively).
The authors report no biomedical financial interests or potential conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Assembled multiple supplementary tables containing data for RNA sequencing and clinical and neuropathological data of human subjects. Table S1a. A list of differentially expressed genes in isolated-microglia from AppNL-G-F/NL-G-F mice. Table S1b. A list of differentially expressed genes in isolated-microglia from rTg4510 mice. Table S1c. A list of differentially expressed genes in isolated-microglia from SOD1G93A mice. Table S1d. Changes in homeostatic microglial markers in each mouse model of neurodegenerative disease. Table S1e. Changes in DAM markers in each mouse model of neurodegenerative disease. Table S1f. Clinical and neuropathological data of human subjects. Table S1g. A list of differentially expressed genes in the precuneus of AD brain. Table S1h. A list of differentially expressed genes enriched in neurons in the precuneus of AD brain. Table S1i. Changes in the genes defined as risk factors for AD in the AD precuneus and microglia isolated from AppNL-G-F/NL-G-F and rTg4510 mice. Table S1j. A list of commonly upregulated genes in the precuneus of AD brain and isolated-microglia from AppNL-G-F/NL-G-F mice and/or rTg4510 mice. Table S1k. A list of commonly downregulated genes in the precuneus of AD brain and isolated-microglia from AppNL-G-F/NL-G-F mice and/or rTg4510 mice.
About this article
Cite this article
Sobue, A., Komine, O., Hara, Y. et al. Microglial gene signature reveals loss of homeostatic microglia associated with neurodegeneration of Alzheimer’s disease. acta neuropathol commun 9, 1 (2021). https://doi.org/10.1186/s40478-020-01099-x