Skip to main content

The cytokines interleukin-6 and interferon-α induce distinct microglia phenotypes

Abstract

Background

Elevated production of the cytokines interleukin (IL)-6 or interferon (IFN)-α in the central nervous system (CNS) is implicated in the pathogenesis of neurological diseases such as neuromyelitis optica spectrum disorders or cerebral interferonopathies, respectively. Transgenic mice with CNS-targeted chronic production of IL-6 (GFAP-IL6) or IFN-α (GFAP-IFN) recapitulate important clinical and pathological features of these human diseases. The activation of microglia is a prominent manifestation found both in the human diseases and in the transgenic mice, yet little is known about how this contributes to disease pathology.

Methods

Here, we used a combination of ex vivo and in situ techniques to characterize the molecular, cellular and transcriptomic phenotypes of microglia in GFAP-IL6 versus GFAP-IFN mice. In addition, a transcriptomic meta-analysis was performed to compare the microglia response from GFAP-IL6 and GFAP-IFN mice to the response of microglia in a range of neurodegenerative and neuroinflammatory disorders.

Results

We demonstrated that microglia show stimulus-specific responses to IL-6 versus IFN-α in the brain resulting in unique and extensive molecular and cellular adaptations. In GFAP-IL6 mice, microglia proliferated, had shortened, less branched processes and elicited transcriptomic and molecular changes associated with phagocytosis and lipid processing. In comparison, microglia in the brain of GFAP-IFN mice exhibited increased proliferation and apoptosis, had larger, hyper-ramified processes and showed transcriptomic and surface marker changes associated with antigen presentation and antiviral response. Further, a transcriptomic meta-analysis revealed that IL-6 and IFN-α both contribute to the formation of a core microglia response in animal models of neurodegenerative and neuroinflammatory disorders, such as Alzheimer’s disease, tauopathy, multiple sclerosis and lipopolysaccharide-induced endotoxemia.

Conclusions

Our findings demonstrate that microglia responses to IL-6 and IFN-α are highly stimulus-specific, wide-ranging and give rise to divergent phenotypes that modulate microglia responses in neuroinflammatory and neurodegenerative diseases.

Background

Interleukin (IL)-6 and interferon (IFN)-α are cytokines that have essential roles in regulating inflammatory processes in the central nervous system (CNS) [1]. The chronic production of these cytokines is directly involved in the pathogenesis of several neuroinflammatory and neurodegenerative diseases, such as neuromyelitis optica spectrum disorder (NMOSD) and cerebral interferonopathies, respectively [1, 2]. For example, in the case of NMOSD, antibody-mediated neutralization of the IL-6 receptor ameliorates disease [3] while inhibition of the IFN-α signaling molecule JAK1 has shown promising outcomes in patients with cerebral interferonopathies [4]. These diseases are debilitating, there is no cure and additional research is required to clarify the underlying pathology to increase the efficacy and precision of new therapeutics [3, 5, 6].

Important clinical and pathological phenotypes of these diseases are recapitulated in transgenic mouse models with CNS-restricted, astrocyte-targeted production of IL-6 or IFN-α, termed glial fibrillary acidic protein (GFAP)-IL6 or GFAP-IFN mice, respectively [1, 7, 8]. GFAP-IL6 mice exhibit signs of disease such as impaired hippocampal long-term potentiation [9, 10], neuronal hyperexcitability [10, 11], progressive cognitive decline [12] and development of ataxia [2, 8]. Pathological changes in the CNS of GFAP-IL6 mice include neurodegeneration and demyelination [12,13,14], reactive astrogliosis and microgliosis [8, 12, 14], upregulation of acute-phase response proteins [8, 15, 16], proliferative angiopathy [8] and breakdown of blood–brain barrier integrity [13, 17]. This is similar to patients with NMOSD, who have increased intrathecal levels of IL-6 [18,19,20,21], reactive astrogliosis and microgliosis and destructive demyelination [22,23,24,25] and can clinically present with ataxia and seizures [26, 27]. On the other hand, GFAP-IFN mice exhibit a cerebral interferonopathy characterized by stunted growth, neuronal hyperexcitability, cognitive dysfunction, ataxia, convulsive seizures and increased mortality [28]. In the CNS, these animals show reactive astrogliosis and microgliosis, neurodegeneration, microangiopathy with aneurysms and cerebral calcification in the cerebellum and thalamus [7, 28, 29], features also commonly seen in patients with cerebral interferonopathies such as Aicardi–Goutières syndrome [30,31,32,33].

Microglia are direct responders to both IL-6 and IFN-α and the increased activation of these cells is a common feature of both IL-6- and IFN-α-mediated diseases in humans and the cytokine-transgenic mice [2]. Yet, little is known about how this feature contributes to these diverse neuropathologies. Microglia are the primary immune cells of the CNS and are highly plastic, with remarkable abilities to fine-tune their molecular and functional phenotype to different states depending on the nature, duration and context of the stimulus. In line with this ability, these cells can exist in a vast spectrum of molecular and functional phenotypes [34,35,36,37,38]. Thus, we hypothesized that microglia are both a prominent target and effector cell of IL-6 and IFN-α in the CNS. We investigated how IL-6 and IFN-α alter the molecular and cellular phenotype of microglia in the GFAP-IL6 versus GFAP-IFN mice. We also correlated the unique phenotypes of microglia in the mouse models to their transcriptomic responses. Finally, we performed a transcriptomic meta-analysis and determined that microglia from other neuropathological states exhibit IL-6 or IFN-α-like responses.

Methods

Mice

The transgenic MacGreen, GFAP-IL6 and GFAP-IFN mice were described previously [7, 8, 28, 39] and were bred and maintained under specific-pathogen-free conditions at the animal facility of the University of Sydney. Both GFAP-IL6 and GFAP-IFN mice were originally developed by I.L. Campbell at the Scripps Research Institute, La Jolla, CA, USA and breeding stock were obtained from there. For the RNA-seq and flow cytometry experiments, GFAP-IL6 and GFAP-IFN mice were crossed with MacGreen mice, which encode the enhanced green fluorescent protein (eGFP) gene under the control of the Csf1r promoter, labeling the myeloid cell compartment, including microglia and macrophages, with eGFP [39]. MacGreen and GFAP-IL6 mice were on C57BL/6 background and GFAP-IFN mice were on a mixed C57BL/6 × BALB/c background. Wildtype littermates from both GFAP-IL6 and GFAP-IFN lines were used as WT controls and no differences were found between lines. The C57BL/6 WT controls for RNA-seq were validated by performing RTPCR for key microglia genes in WT controls from both lines. Animals received food and water ad libitum. The temperature and humidity parameters in animal holding areas were set to fall between 20–24 °C and 40–70%, respectively, with light between 0545 and 1745 h. Mice were housed at a maximum density of 6 mice per cage.

For the microglia cell turnover experiment, 1-, 3- and 6-month-old mice were injected intraperitoneally with 100 mg/kg 5-bromo-2’-deoxyuridine (BrdU) in 0.9% (w/v) NaCl each day at 1500 h for 5 days to label proliferating cells. For histological analysis, mice were deeply anaesthetized with isoflurane, perfused intracardially with ice-cold phosphate-buffered saline (PBS) followed by ice-cold, neutral buffered 4% paraformaldehyde (PFA). Brains were collected, fixed in 4% neutral buffered PFA overnight at 4 °C and then paraffin-embedded. For passive tissue clearing, PFA-fixed brains were dissected into the cerebellum and cortex and were transferred to PBS. For Western blot, 1-month-old mice were euthanized with isoflurane and the cerebellum was dissected and flash frozen. For ex vivo experiments, mice were deeply anaesthetized with isoflurane, perfused intracardially with ice-cold PBS and the brain was isolated.

Histology

For histochemistry/immunohistochemistry (HC/IHC), paraffin sections (12 μm) were deparaffinized and rehydrated in graded ethanol. Antigens were unmasked with 25 mM Tris pH 8, 5 mM EDTA pH 8 and 0.05% (w/v) SDS in a vegetable steamer for 40 min. Sections were incubated in 0.3% peroxidase for 10 min and blocked in 1% goat serum with 0.1% Triton X-100 and 0.05% Tween-20 in PBS for 30 min. The primary antibodies rabbit anti-pY701-STAT1 (sc135648, Santa Cruz Biotechnology, 1:100) or rabbit anti-pY705-STAT3 (9145L, Cell Signaling Technologies, 1:100) were incubated overnight at 4 °C. Following primary antibody incubation and washing, sections were incubated with biotinylated secondary goat anti-rabbit antibodies (BA-1000, Vector Laboratories, 1:200) followed by VECTASTAIN Elite ABC Horseradish Peroxidase (HRP) Kit (PK-7200, Vector Laboratories) for 1 h at room temperature (RT). Sections were developed with 3,3’-diaminobenzidine with nickel enhancement (SK-4100, Vector Laboratories). Following washing, avidin/biotin blocking (SP-2001, Vector Laboratories) was performed. Biotinylated-tomato lectin (L0651, Sigma-Aldrich, 1:50) was incubated overnight at 4 °C. Sections were then washed and incubated with HRP-conjugated streptavidin (SA-5004, Vector Laboratories, 1:200) for 1 h at RT. Sections were developed with 3,3′-diaminobenzidine without nickel enhancement, counterstained with Mayer’s hematoxylin and mounted. Sections were viewed with a DM4000B microscope (Leica Microsystems) and imaged using a SPOT Flex 15.2 64 Mp Shifting Pixel camera and SPOT Advanced 4.5 software (Diagnostic Instruments).

For immunofluorescence, deparaffinization, rehydration, antigen retrieval and blocking were performed as described above. The primary antibodies rabbit anti-Iba1 (019–19741, Wako Pure Chemical Industries, 1:500) and rat anti-BrdU (MCA2060GA, BioRad Laboratories, 1:100) were incubated overnight at 4 °C. Following primary antibody incubation and washing, sections were incubated with secondary goat anti-rabbit IgG-AF488 (A-11034, Thermo Fisher Scientific, 1:500) and goat anti-rat IgG-AF594 (A-11007, Thermo Fisher Scientific, 1:500) for 1 h at RT. Slides were washed and then cover-slipped with Fluoroshield™ DAPI mounting media (Sigma-Aldrich). For the TUNEL assay, the primary antibody rabbit anti-Iba1 (1:500) was incubated for 2 h at RT. Following washing, sections were incubated in TUNEL reaction mixture (TMR Red, Roche) and secondary goat anti-rabbit IgG-AF488 (1:500) for 1 h at 37 °C. Slides were washed and then cover-slipped with Fluoroshield™ DAPI mounting media (Sigma-Aldrich). Fluorescent imaging was performed at the Advanced Microscopy Facility of the Bosch Institute at the University of Sydney using a Zeiss AxioScan.Z1 slide-scanning microscope at 20 × magnification and a Zeiss LSM800 confocal laser scanning microscope using 20 × Plan Apochromat NA = 0.8 air objective, 40 × Plan Apochromat NA = 1.3 oil-immersion objective, or 63 × Plan Apochromat NA = 1.4 oil-immersion objective with 405, 488 and 561 nm lasers and appropriate filters (Carl Zeiss).

Immunoblotting

Total protein for immunoblotting was isolated from the cerebella of mice in 100 mM Tris pH 7.5, 150 mM NaCl, 1 mM EDTA, 1% deoxycholic acid, 1% Triton X-100, 0.1% (w/v) SDS, 2 mM PMSF, 50 mM NaF, and 1 × Protease Inhibitor Cocktail III (539134, Merck Millipore) and 1 × Phosphatase Inhibitor Cocktail II (524625, Merck Millipore). Proteins were separated on a 10% Tris–glycine polyacrylamide gel by electrophoresis and transferred onto polyvinylidene fluoride membranes. The primary antibodies rabbit anti-pY701-STAT1 (7649, Cell Signaling Technology, 1:2000), rabbit anti-pS727-STAT1 (9177, Cell Signaling Technology, 1:2000), rabbit anti-pY705-STAT3 (9131, Cell Signaling Technology, 1:2000), rabbit anti-STAT1 (9172, Cell Signaling Technology, 1:2000), rabbit anti-STAT3 (4904, Cell Signaling Technology, 1:2000) and mouse anti-GAPDH (MAB374, Merck Millipore, 1:30,000) in 5% (w/v) skim milk powder in Tris-buffered saline and 0.1% Tween-20 were incubated with membranes overnight at 4 °C. Membranes were incubated with HRP-conjugated goat anti-rabbit IgG (7074P2, Cell Signaling Technology, 1:10,000) and peroxidase-conjugated goat anti-mouse IgG (A0168, Sigma-Aldrich, 1:10,000) for 1 h at RT. Proteins were detected by chemiluminescence (WBKLS0500, Merck Millipore) and visualized using an iBright™ 1500 (Thermo Fisher Scientific). Relative protein bands were quantified by densitometry using ImageJ software (NIH, USA) by normalizing band density to that of the corresponding GAPDH loading control.

Tissue clearing and microglia morphological analysis

PFA-fixed cerebella and cortices were incubated in 4% (w/v) acrylamide, 0.05% (w/v) bisacrylamide and 0.25% (w/v) VA-044 (Wako Pure Chemical Industries) in PBS for 48 h at 4 °C and were then incubated at 37 °C for 3 h to polymerize the hydrogel-tissue matrix. Polymerized tissues were next briefly washed, transferred into 8% (w/v) SDS in PBS and incubated at 37 °C with shaking at 225 rpm until tissues became transparent (4–5 days). Cleared tissues were then washed in PBS for 6 h at RT, with fresh changes of PBS each hour. Rabbit anti-Iba1 (019–19741, Wako Pure Chemical Industries, 1:200) in 2% goat serum with 0.1% Triton X-100 and 0.01% sodium azide in PBS was incubated with tissues at RT with shaking for 7 days. Following 6 h of washing, tissues were then incubated with secondary goat anti-rabbit IgG-AF594 (A-11037, Thermo Fisher Scientific, 1:200) at RT with shaking for 7 days. Following washing, tissues were incubated overnight in 88% (w/v) Histodenz (Sigma-Aldrich), 0.01% sodium azide and 1 μg/mL Hoechst 33342 (B2261, Sigma-Aldrich) in PBS at RT with shaking. Imaging was performed on a Zeiss LSM800 confocal laser scanning microscope using 20 × Plan Apochromat NA = 0.8 air objective with 405 and 568 nm lasers and appropriate filters (Carl Zeiss). Three-dimensional reconstructions of microglia were generated using the Filament Tracer feature in Imaris 9.2.1 (Oxford Instruments). Following surface and dendrite rendering, morphometric analysis was performed using the in-house statistics generator. The total process length, total process volume, number of branching points, number of terminal points, and total Sholl intersections were analyzed as indicators of microglia process size and complexity.

Fluorescence-activated cell sorting (FACS) of microglia

PBS-perfused brains were washed in PBS, the cerebellum was dissected and cells were mechanically dissociated by Dounce homogenization and passed through a 70 μm sieve. Cells were pelleted by centrifugation at 460g for 10 min at 4 °C before resuspension in 70% (v/v) Percoll. An equal volume of 37% (v/v) Percoll was slowly layered on top of the cell suspension to form a discontinuous Percoll gradient. Microglia and leukocytes were collected from the 37/70% (v/v) interface after centrifugation at 1825g for 25 min at 25 °C with no brake. Cells were then stained on ice with rat anti-4D4 (clone 4D4, 1:300 [40], provided by Oleg Butovsky) in 5% FBS (Thermo Fisher Scientific) and 5 mM EDTA in PBS for 25 min. Following washing, cells were then incubated on ice with goat anti-rat IgG-APC (Poly4054, Biolegend, 1:300). Cells were washed and then incubated on ice with 7AAD (51-68981E, BD Bioscience, 1:20) for 10 min. Using a BD Influx™ Cell Sorter (BD Bioscience), eGFP+ 4D4+ 7AAD cerebellar microglial cells were sorted into lysis buffer. Total RNA was isolated from FACS-sorted microglia using the RNeasy® plus micro kit (QIAGEN Inc.) according to the manufacturer’s instructions.

Library preparation for RNA-seq

The cDNA was synthesized from cerebellar microglial cell RNA at the Western Sydney University Next Generation Sequencing facility using SMARTer® Stranded Total RNA-Seq Kit—Pico Input Mammalian (Takara Bio Inc.) according to manufacturer’s instructions. The RNA-seq library fragments were amplified by PCR, using 14 enrichment cycles, before purification and RNA-seq.

RNA-seq

Paired-end, 126 bp read length RNA-seq was performed on a HiSeq 2500 (Illumina) at the Western Sydney University Next Generation Sequencing facility.

RNA-seq analysis

RNA-seq data were processed using a custom pipeline designed by Boris Guennewig [41,42,43]. The quality of the input data was assessed using FastQC (version 0.11.3) [44] and reads were mapped to the mm10 reference genome using the STAR aligner (version 2.5.2a) [45]. Potential transcripts were identified using StringTie (version 1.3.3b) [46] and known GENCODE genes were quantified using RSEM (version 1.3.0) [47]. For differential expression analysis, RSEM count data were imported into the R project environment [48]. Reads with primarily zero counts were filtered out to obtain the genes for downstream analyses. Outlier samples were identified using principal component analysis (PCA) and hierarchical clustering. Reads were normalized using the trimmed mean of means (TMM) and differential expression analysis was performed in R using the edgeR package [49]. A differentially expressed gene (DEG) was called significant if its false discovery rate (FDR) was ≤ 0.05 after Benjamini–Hochberg correction. Significantly enriched biological processes were identified using the gene ontology (GO) enrichment tool WebGestalt [50] and clustering of biological processes was visualized using Enrichment Map [51]. Two-way analysis was performed using the DEG lists generated by edgeR. Genes that were at least twofold differentially expressed by GFAP-IL6 versus WT microglia and/or genes at least twofold differentially expressed by GFAP-IFN versus WT microglia were plotted. A constant of 0.125 was added to the normalized reads prior to the fold-change calculation to avoid infinite fold-change values for genes with zero counts across a group. “Core response” genes were differentially expressed in both GFAP-IL6 versus WT microglia and GFAP-IFN versus WT microglia, with no significant difference in expression between GFAP-IL6 versus GFAP-IFN microglia. “IL-6-skewed” genes were differentially expressed in GFAP-IL6 versus WT microglia, but not GFAP-IFN versus WT microglia, or were differentially expressed by both GFAP-IL6 versus WT and GFAP-IFN versus WT microglia but were differentially expressed by GFAP-IL6 versus GFAP-IFN microglia, or were genes that were significantly upregulated by GFAP-IL6 versus WT microglia and significantly downregulated by GFAP-IFN versus WT microglia. “IFN-α-skewed” genes were differentially expressed in GFAP-IFN versus WT microglia, but not GFAP-IL6 versus WT microglia, or were differentially expressed by both GFAP-IFN versus WT and GFAP-IL6 versus WT microglia but were also differentially expressed by GFAP-IFN versus GFAP-IL6 microglia, or were genes that were significantly upregulated by GFAP-IFN versus WT microglia and significantly downregulated by GFAP-IL6 versus WT microglia. Significantly enriched biological processes were identified and visualized as above.

RNA-seq meta-analysis

In addition to our RNA-seq data, RNA-seq FASTQ files were downloaded from the sequence read archive (SRA) and processed using the custom pipeline designed by Boris Guennewig. The files downloaded from SRA included WT Clec7a, APP-PS1 Clec7a and APP-PS1 Clec7a+ microglia (GSE102563, n = 6 per group) [52]; non-transgenic (Non-Tg) and hMAPT-P301S microglia (GSE93180, n = 6 per group) [53]; unmanipulated, EAE CD11c and EAE CD11c+ microglia (GSE78809, n = 3 per group) [54]; and vehicle and LPS microglia (GSE75246, n = 5 per group) [55]. Data processing and DEG analysis was performed as described above and was performed separately for each study. Normalized gene reads were then Z-score transformed within each study before the data from the different studies were combined. Z-score transformation has been used previously to account for the “batch” effect of combining data from different studies [53]. The pheatmap package was used in R [56] to perform clustering analysis and generate the heatmap. For each gene cluster, significantly enriched biological processes were identified using WebGestalt [50]. For the genes in each cluster, we calculated log2-fold changes by comparing each sample group with its respective control, as well as the median log2-fold change for each comparison. A constant of 0.125 was added to the normalized reads prior to the fold-change calculation to avoid infinite fold-change values for genes with zero counts across a group.

Microglia cDNA synthesis and real-time PCR (RTPCR) analysis

Following FACS-sorting of microglia from the whole brain and RNA isolation as above, cDNA was synthesized using RevertAid First Strand cDNA synthesis kit (Thermo Fisher Scientific) according to the manufacturer’s instructions. For each sample, we used cDNA from the same number of cells to normalize the RTPCR. Accordingly, the cDNA was diluted such that RTPCR was performed on cDNA from 250 microglial cells. The RTPCR was set up using SensiFAST™ SYBR® Lo-ROX (Meridian Bioscience) and 400 nM primer pairs and was performed with a 7500 Fast Real-Time PCR System (Thermo Fisher Scientific) using the ΔΔCt setting with the cycle program: 95 °C for 2 min, then 40 cycles of 95 °C for 5 s then 60 °C for 30 s, followed by melt curve analysis. The CT for each gene of interest was normalized to the CT of the housekeeping gene 18S rRNA. Primer sequences: Apoe (forward: TGTGGGCCGTGCTGTTGGTC; reverse: GCCTGCTCCCAGGGTTGGTTG) [57], Axl (forward: TGAGCCAACCGTGGAAAGAG; reverse: AGGCCACCTTATGCCGATCTA) [58], B2m (forward: GTGACCCTGGTCTTTCTGGT; reverse: GTATGTTCGGCTTCCCATTC) [59], C4b (forward: GACAAGGCACCTTCAGAACC; reverse: CAGCAGCTTAGTCAGGGTTACA), Cst3 (forward: CGCTCCTTGCTGTTCCTGCT; reverse: TGCCCTTGTTGTACTCGCTCAC) [60], Ctsb (forward: AGACCTGCTTACTTGCTGTG; reverse: GGAGGGATGGTGTATGGTAAG) [61], Fn1 (forward: ACCGACAGTGGTGTGGTCTA; reverse: CACCATAAGTCTGGGTCACG) [62], H2-D1 (forward: TCCGAGATTGTAAAGCGTGAAGA; reverse: GAACCCAAGCTCACAGGGAA) [63], Slc39a14 (forward: GGAACCCTCTACTCCAACGC; reverse: ATGGTTATGCCCGTGATGGT) [64] and 18S (forward: CACGGCCGGTACAGTGAAAC; reverse: AGAGGAGCGAGCGACCAA) [65].

Flow cytometry

Following ex vivo isolation, single-cell suspensions of microglia and leukocytes from the brain were incubated at 4 °C for 30 min with rat anti-CD16/32-BB700 (2.4G2, BD Bioscience, 1:200) and LIVE/DEAD™ Fixable Blue Dead Cell Stain (Thermo Fisher Scientific, 1:500). Cells were washed and then stained at 4 °C for 30 min with fluorophore-conjugated antibodies: rat anti-B220-BUV661 (RA3-6B2, BD Bioscience, 1:200), rat anti-CD115-AF594 (AFS98, Biolegend, 1:200), rat anti-CD11b-BUV737 (M1/70, BD Bioscience, 1:200), Armenian hamster anti-CD11c-APC-R700 (N418, BD Bioscience, 1:200), hamster anti-CD3e-PE/Cy5.5 (145-2C11, Thermo Fisher Scientific, 1:200), rat anti-CD4-BV570 (RM4-5, Biolegend, 1:200), rat anti-CD45-APC/Cy7 (30-F11, BD Bioscience, 1:200), rat anti-CD62L-BV650 (MEL-14, Biolegend, 1:200), mouse anti-CD64-APC (X54-5/7.1, BD Bioscience, 1:200), rat anti-CD69-BV785 (H1.2F3, BD Bioscience, 1:200), hamster anti-CD80-PE/CF594 (16-10A1, BD Bioscience, 1:200), rat anti-CD86-BV605 (GL1, Biolegend, 1:200), rat anti-CD8a-BUV805 (53–6.7, BD Bioscience, 1:200), rat anti-F4/80-BUV395 (T45-2342, BD Bioscience, 1:200), rat anti-Ly6C-PE/Cy7 (HK1.4, Biolegend, 1:200), rat anti-Ly6G-BUV563 (1A8, BD Bioscience, 1:200), rat anti-MHC-II (I-A/I-E)-BV510 (M5/114.15.2, BD Bioscience, 1:200), mouse anti-NK1.1-PE/Cy5 (PK136, Biolegend, 1:200), rat anti-SCA-1-BV711 (D7, Biolegend, 1:200) and rabbit anti-TMEM119-PE (106–6, Abcam, 1:200). Cells were then fixed in 4% PFA prior to analysis performed on an LSR-X flow cytometer (BD Bioscience) and data were analyzed with FlowJo software version 10 (BD Bioscience).

Alternatively, cells were incubated at 4 °C for 30 min with unconjugated rat anti-CD16/32 (2.4G2, BD Bioscience, 1:200) and LIVE/DEAD™ Fixable Blue Dead Cell stain (Thermo Fisher Scientific, 1:500). Following washing, cells were stained at 4 °C for 25 min with rat anti-FCRLS (1:300 [66], provided by Oleg Butovsky) or rat anti-4D4 (1:300 [40]), followed by staining at 4 °C for 20 min with secondary goat anti-rat IgG-APC (Poly4054, Biolegend, 1:300). Cells were washed and then incubated with mouse anti-MHC-I (H-2 Kb/H-2Db)-PE (28-8-6, Biolegend, 1:100) and rat anti-CD11b-PerCP/Cy5.5 (M1/70, BD Bioscience, 1:200) at 4 °C for 30 min. Cells were then fixed in 4% PFA prior to analysis performed on an LSR-II flow cytometer (BD Bioscience) and data were analyzed with FlowJo software version 10 (BD Bioscience).

Computational analysis of flow cytometry data

Computational analysis of flow cytometry data was performed using the Spectre R package [67] (package publicly available: https://github.com/ImmuneDynamics/Spectre). Live CD45+ and/or CD11b+ cells were manually gated and exported as CSV-channel files. Keywords denoting the sample and group names were then added to the samples, before the samples were merged into a single data table. The Flow Self-Organizing Maps (FlowSOM) algorithm [68] was then run on the merged dataset to cluster the dataset, where every cell is assigned to a specific cluster. The data were downsampled so that the relative number of cells in each sample was represented proportionally. Subsequently, the downsampled data were analyzed by the dimensionality reduction algorithm Uniform Manifold Approximation and Projection (UMAP) [69] for cellular visualization. Following FlowSOM clustering and dimensionality reduction with UMAP, summary tables containing expression level, cell frequency and cell number data of both the large FlowSOM and smaller UMAP datasets were exported and the total cell number of each cell type cluster was calculated.

Statistics

Results are presented as individual values per mouse and mean ± SEM. The specific statistical tests used to determine significance are indicated in the figure legends. Statistical calculations were performed using Prism version 9 (GraphPad Software). For all data comparisons, a p-value < 0.05 was considered statistically significant.

Results

Activation of specific cytokine signaling pathways in the brain of GFAP-IL6 and GFAP-IFN mice

Both IL-6 and IFN-α signal via the JAK/STAT signal transduction pathway [70, 71], with canonical IL-6 signaling principally mediated through tyrosine 705 phosphorylated STAT3 (pY705-STAT3), whereas both tyrosine 701 and serine 727 phosphorylation of STAT1 (pY701-STAT1 and pS727-STAT1) are required for maximal activation of the IFN-α response [72, 73]. To assess stimulus-specific responses, we first performed immunoblots for phosphorylated STAT3 and STAT1 on cerebellum from GFAP-IL6 and GFAP-IFN mice (Additional file 1: Figs. S1, S2) since transgene expression of both cytokines is highest in this region compared with other areas of the brain [7, 8, 15, 74]. In agreement with previous reports [17, 75], GFAP-IL6 mice had high levels of pY705-STAT3 and low levels of pY701- and pS727-STAT1 in the cerebellum, while the cerebellum of GFAP-IFN mice had low levels of pY705-STAT3 and high levels of pY701- and pS727-STAT1 (Additional file 1: Fig. S1a). We next performed dual-label immunohistochemistry/histochemistry for pY705-STAT3 or pY701-STAT1 with tomato lectin, a microglia marker. In the cerebellum of WT mice, pY705-STAT3 and pY701-STAT1 were not detectable (Additional file 1: Fig. S1b, c). By contrast, since all resident cell types in the CNS are capable of responding to IL-6 (via trans-signaling) and IFN-α [2], the cerebellum of GFAP-IL6 and GFAP-IFN mice contained lectin-positive microglia and other, lectin-negative CNS-resident cells, with nuclear pY705-STAT3 or pY701-STAT1, respectively (Additional file 1: Fig. S1d–g). Importantly, GFAP-IL6 microglia had strong nuclear staining for pY705-STAT3, while nuclear pY701-STAT1 was not detected (Additional file 1: Fig. S1d, e). On the other hand, GFAP-IFN microglia had strong nuclear staining for pY701-STAT1, while nuclear pY705-STAT3 was not detected (Additional file 1: Fig. S1f, g). These findings highlight the remarkable signaling specificity of microglia in response to the neuroinflammation induced by chronic IL-6 versus IFN-α production, despite the pleiotropic effects of these cytokines and the secondary inflammatory factors which are induced by the chronic neuroinflammation in these animals. Together, our findings are consistent with microglia mounting stimulus-specific responses to transgene-driven production of IL-6 versus IFN-α. Further, while both GFAP-IL6 and GFAP-IFN mice exhibited robust microgliosis, with increased lectin binding compared with WT (Additional file 1: Fig. S1b-g), these cells also appeared to exhibit dramatic differences in number and morphology. Therefore, we next aimed to precisely dissect the molecular and cellular changes made by microglia in response to chronic IL-6 versus IFN-α production in the brain.

Microglia of GFAP-IL6 and GFAP-IFN mice have distinct turnover patterns

Given the apparent stimulus-specific cytokine modulation of microglia cell numbers in the brains of GFAP-IL6 versus GFAP-IFN mice, we questioned whether there were changes to microglia turnover in response to chronic IL-6 versus IFN-α production. Indeed, the number of microglia in the cerebellum, cortex and hippocampus of the GFAP-IL6 mice was increased compared with both WT and GFAP-IFN mice at all ages studied (Fig. 1a–i). Although transgene expression of IL-6 is lower in the cortex and hippocampus than the cerebellum [7, 8, 15, 74], GFAP-IL6 mice still had significantly increased numbers of microglia in these regions compared with WT mice (Fig. 1h–i). Notably, GFAP-IL6 mice had significantly greater numbers of microglia compared with WT mice at all ages. However, the number of microglia decreased with age, with significantly fewer cells in the cortex and hippocampus of 3-month-old mice and significantly fewer cells in the cerebellum and cortex of 6-month-old mice compared with 1-month-old GFAP-IL6 mice. In contrast, microglia numbers in the cerebellum of GFAP-IFN mice were slightly, but not significantly, increased compared with WT mice and these numbers remained largely unchanged at all ages studied. The density of microglia in discrete brain regions also differed by genotype. Significantly greater numbers of microglia were seen in the cerebellum of GFAP-IL6 mice than in the cortex and hippocampus (Fig. 1g–i). By contrast, the density of microglia was roughly equal in the cerebellum, cortex and hippocampus of WT and GFAP-IFN mice, with the exception of the hippocampus of 3-month-old GFAP-IFN mice, which had significantly less microglia compared with the cerebellum and cortex.

Fig. 1
figure 1

Microglia in the brain of GFAP-IL6 versus GFAP-IFN mice have unique turnover patterns. a–c Representative immunofluorescence images (Iba1+ microglia, green; BrdU+, red; DAPI, blue) from the cerebellum of 1-month-old WT (a), GFAP-IL6 (b) and GFAP-IFN (c) mice. d–f Representative images (Iba1+ microglia, green; TUNEL+, red; DAPI, blue) from the hippocampus of 6-month-old WT (d), GFAP-IL6 (e) and GFAP-IFN (f) mice. Scale bars, 20 μm. g–o Quantification of the total number of Iba1+ microglia per mm2 (g-i), the number of Iba1+BrdU+ microglia per mm2 (j–l) and the number of Iba1+TUNEL+ apoptotic microglia per section (m–o) in the cerebellum (g, j, m), cortex (h, k, n) and hippocampus (i, l, o) at 1, 3 and 6 months of age. n = 3–5 mice/group. Graphs show individual values per mouse and mean ± SEM. *, p < 0.05 compared with WT of the same age; ^, p < 0.05 compared with GFAP-IL6 of the same age; #, p < 0.05 compared with 1-month-old of the same genotype; x, p < 0.05 compared with 3-month-old of the same genotype using two-way ANOVA with Tukey’s post-test

To determine the basis for the differences in microglia numbers between mouse lines and age, we assessed proliferation and apoptosis rates. In all three brain regions of WT mice, the numbers of proliferative microglia as evaluated by BrdU incorporation and immunostaining (Fig. 1j–l) and apoptotic microglia as evaluated by TUNEL staining were very low (Fig. 1m–o). There were significantly increased numbers of BrdU+ microglia in the cerebellum of GFAP-IL6 mice as compared with WT mice at all ages (Fig. 1j). Interestingly, at later ages, the cerebellum of GFAP-IL6 mice contained significantly fewer BrdU+ microglia as compared with 1-month-old GFAP-IL6 mice. There were also low numbers of BrdU+ microglia in the cortex and hippocampus (Fig. 1k, l), consistent with lower numbers of total microglia in these regions in GFAP-IL6 mice (Fig. 1g–i). In addition, there was no significant increase in the number of TUNEL+ microglia in the brain of GFAP-IL6 mice (Fig. 1m–o). On the other hand, while the cerebellum of 1-month-old GFAP-IFN mice had comparable numbers of BrdU+ microglia to 1-month-old GFAP-IL6 mice, there was a small, but not statistically significant, increase in the number of BrdU+ microglia in the cerebellum of GFAP-IFN mice at later ages (Fig. 1j). Furthermore, in the cortex and hippocampus, GFAP-IFN mice had significantly increased numbers of BrdU+ microglia at all ages examined compared with WT mice (Fig. 1k, l). In addition, GFAP-IFN mice had a progressive increase in the number of TUNEL+ microglia, particularly in the hippocampus (Fig. 1m–o). Taken together, these findings indicate that the increase in microglia number and density in GFAP-IL6 mice is accounted for by increased proliferation but not apoptosis. On the other hand, despite increased proliferation, the number and density of microglia in GFAP-IFN mice are relatively unchanged, likely due in part to increased apoptosis.

Microglia in the brains of GFAP-IL6 and GFAP-IFN mice exhibit unique morphological changes in response to chronic IL-6 versus IFN-α production

Microglia morphology in GFAP-IL6 and GFAP-IFN mice appeared distinct (Fig. 1a–f, Additional file 1: Fig. S1b–g). To quantify morphological differences of microglia in our mouse models, we generated three-dimensional reconstructions of Iba1-immunostained microglia from the cerebellum and cortex of GFAP-IL6 and GFAP-IFN mice (Fig. 2a, b). Compared with WT microglia, the total process length of cerebellar and cortical GFAP-IL6 microglia, as well as the number of branching points, terminal points and the total number of Sholl intersections, were reduced (Fig. 2c, d). This was most pronounced in the cerebellum (Fig. 2c). On the other hand, in GFAP-IFN mice, cerebellar and cortical microglial cell total process length and the number of branching points, terminal points and the total number of Sholl intersections were significantly increased compared with both WT and GFAP-IL6 microglial cells at all ages (Fig. 2c, d). Taken together, the cytokine environments induced by chronic IL-6 versus IFN-α production in the brain of GFAP-IL6 versus GFAP-IFN mice cause microglia to adopt distinct morphological states.

Fig. 2
figure 2

Microglia have distinct morphologies in GFAP-IL6 versus GFAP-IFN mice. a-b Representative three-dimensional reconstructions of Iba1+ microglia from the cerebellum (a) and cortex (b) of 1- and 6-month-old WT, GFAP-IL6 and GFAP-IFN mice. c-d Imaris automated quantification of the total process length, branching points, terminal points and total Sholl intersections of cerebellar (c) and cortical (d) microglia. n = 2–3 mice/group, n = 15–31 cells/genotype. Graphs show individual values per cell and mean ± SEM. *, p < 0.05 compared with WT of the same age; ^, p < 0.05 compared with GFAP-IL6 of the same age; #, p < 0.05 compared with 1-month-old of the same genotype; x, p < 0.05 compared with 3-month-old of the same genotype using two-way ANOVA with Tukey’s post-test

Divergent transcriptional landscapes arise in microglia from GFAP-IL6 versus GFAP-IFN mice

To examine the effects of chronic IL-6 versus IFN-α signaling in the brain on the microglial cell phenotype on the transcriptional level, GFAP-IL6 and GFAP-IFN mice were crossed with MacGreen mice to label the myeloid compartment including microglia with eGFP [39]. We further used the microglia-specific 4D4 antibody, which does not bind to recruited monocytes/macrophages [40], to isolate and purify microglia from the cerebellum of 1-month-old mice by FACS of live dual-labeled eGFP+ 4D4+ cells and performed RNA-seq (Fig. 3a). We focused on this early age as the direct effects of IL-6 versus IFN-α are more discernible, since neuropathological changes including inflammation and neurodegeneration that arise because of chronic cytokine signaling are modest at this age. The purity of the microglial cell preparation was confirmed by robust expression of microglia signature genes, such as Cx3cr1, Csf1r, Hexb, Olfml3 and Tmem119, while expression of genes specific for astrocytes, oligodendrocytes, neurons, endothelial cells, pericytes, T cells, B cells, granulocytes and monocytes was not detectable (Fig. 3b). We next performed principal component analysis (PCA) on all genes that passed the expression level criteria to examine the overall differences in the transcriptional landscape of cerebellar microglia (Fig. 3c). PCA showed that cerebellar microglia from WT, GFAP-IL6 and GFAP-IFN mice all had a high degree of separation from one another, indicative of highly divergent transcriptomes in these cells. Compared with WT microglia, microglia from GFAP-IL6 mice had 445 upregulated and 439 downregulated genes, with upregulation of genes classically associated with an IL-6-response, including Spp1, Fn1, Socs3, Saa3 and Apoe (Fig. 3d, Additional file 2: Table S1). The IL-6-regulated differentially expressed genes (DEGs) were enriched in gene ontology (GO) biological processes such as immune response, myeloid cell differentiation, motility, cell turnover, phagocytosis, metabolism and antigen presentation (Fig. 3e). On the other hand, in response to the chronic production of IFN-α in the brain of GFAP-IFN mice, cerebellar microglia had 869 upregulated and 680 downregulated genes compared with WT cells, with upregulation of genes associated with responsiveness to IFN-α, including Ifit3, Stat2, B2m, H2-Q7, Oas2 and Slfn5 (Fig. 3f, Additional file 2: Table S1). The IFN-α-regulated DEGs were enriched in biological processes including immune response, viral response, response to IFN, response to cytokine, antigen presentation, regulation of T cell cytotoxicity, metabolism and cell cycle (Fig. 3g).

Fig. 3
figure 3

Cerebellar microglia regulate distinct subsets of genes in response to the cytokine environments induced by chronic IL-6 versus IFN-α signaling. a Microglia were isolated from the cerebellum of 1-month-old MacGreen-WT, -GFAP-IL6 and -GFAP-IFN mice and purified by FACS of live eGFP+ 4D4+ cells. RNA was isolated and reverse transcribed into cDNA, which was then amplified by PCR and sequenced. b Fragments per kilobase of transcript per million mapped reads (FPKM) of cell type-specific genes for microglia, other CNS-resident cells and peripheral leukocytes. c PCA of RNA-seq datasets of cerebellar microglia from WT, GFAP-IL6 and GFAP-IFN mice. d MA plot (representing log-ratio (M) on the y-axis and mean average (A) on the x-axis) showing transcripts differentially expressed by GFAP-IL6 cerebellar microglia compared with WT microglia. e Enrichment map of enriched GO biological processes by WebGestalt generated from the IL-6-regulated DEGs. f MA plot showing transcripts differentially expressed by GFAP-IFN cerebellar microglia compared with WT microglia. g Enrichment map of enriched GO biological processes by WebGestalt generated from the IFN-α-regulated DEGs. For d, f, the number of significantly (FDR < 0.05) upregulated and downregulated genes are indicated. For e, g, nodes in enrichment maps are significantly enriched in GO lists (FDR < 0.05) and were used to name clusters. n = 3 mice/group

We next asked which DEGs are commonly regulated to a similar degree in response to IL-6 or IFN-α and which genes are IL-6- versus IFN-α-skewed. For this, we generated a two-way fold-change plot using genes at least twofold differentially expressed by GFAP-IL6 versus WT microglia and/or genes at least twofold differentially expressed by GFAP-IFN versus WT microglia (Fig. 4a) and performed GO analysis on the upregulated DEGs (Fig. 4b-d). There were no GOs that were significantly enriched by downregulated genes alone. Cerebellar microglia upregulated 144 genes and downregulated 143 genes to a similar degree in response to chronic production of IL-6 or IFN-α, with regulation of genes associated with microglia activation (upregulation of C4b, Bhlhe40, Ccl2, Ccl12 and downregulation of P2ry12) (Fig. 4a, Additional file 3: Table S2). Interestingly, microglia from GFAP-IFN mice had a more pronounced and more extensive response than cells from GFAP-IL6 mice, since the overlapping core response genes comprised only 20.5% of the genes regulated by microglia in response to the IFN-α-induced cytokine environment, while core response genes comprised 37.0% of the genes regulated by microglia in response to the IL-6-induced cytokine environment (Additional file 1: Fig. S3). The core response genes upregulated by microglia from both transgenic mice were significantly enriched in processes including immune response, Fc-receptor signaling, antigen processing and presentation, IL-10 production, nitric oxide synthase biosynthesis, apoptosis, catabolism and translation (Fig. 4b). Since these biological processes are predominantly related to the modulation of the microglia immune response, these common transcriptional changes likely represent the overlapping target genes activated by both IL-6 and IFN-α, as well as the milieus induced by these cytokines within the brain tissue.

Fig. 4
figure 4

Microglia in GFAP-IL6 versus GFAP-IFN mice acquire unique transcriptional programs in addition to a common set of core response genes. a Two-way fold-change plot of differentially expressed genes in GFAP-IL6 versus WT microglia and GFAP-IFN versus WT microglia to identify core response genes (pink), as well as IL-6-skewed (orange) and IFN-α-skewed (blue) genes. b Enrichment map of top 100 significantly enriched GO biological processes by WebGestalt generated from the DEGs that are commonly upregulated by IL-6 and IFN-α. c Enrichment map of top 100 significantly enriched GO biological processes by WebGestalt generated from the DEGs that are upregulated and IL-6-skewed. d Enrichment map of top 100 enriched GO biological processes by WebGestalt generated from the DEGs that are upregulated and IFN-α-skewed. For b–d, nodes in enrichment maps are significantly enriched in GO lists (FDR < 0.05) and were used to name clusters

Two-way analysis also identified 252 upregulated and 237 downregulated genes that were exclusively regulated by IL-6 or were skewed such that they were significantly regulated by IL-6 compared with IFN-α and included Spp1, Apoe, Saa3, Socs3, Fn1, Slc39a14 and Ccnb2 (increased expression) and Ddit3 and Itgae (decreased expression) (Fig. 4a, Additional file 3: Table S2). IL-6-skewed genes that were upregulated by GFAP-IL6 microglia were enriched for biological processes including leukocyte activation, response to cytokine, response to endoplasmic reticulum (ER) stress, motility, neuron death, differentiation, cell signaling and metabolism (Fig. 4c). Consistent with the presence of proliferating microglia in the brain of GFAP-IL6 mice (Fig. 1), genes regulated by GFAP-IL6 microglia were also enriched in processes associated with proliferation (Ccnb2, Cdc42, Fn1, Spp1 and Vegfa) (Fig. 4a, c). In addition, genes that were associated with the microglial cell response to IL-6 were enriched in functional processes including phagocytosis and lipid processing (Apoe, Apoc1, Abca1, Ch25h and Spp1) and import of extracellular material such as iron (Slc39a14 and Trf).

On the other hand, we identified 675 upregulated and 435 downregulated genes that were exclusively regulated by IFN-α or were skewed such that they were significantly regulated by IFN-α compared with IL-6 and included Axl, Cdkn1a, Ifit3, Irf7, H2-D1, Stat2, Tap1 and Usp18 (increased expression) and Cd14 (decreased expression) (Fig. 4a, Additional file 3: Table S2). IFN-α-skewed genes that were upregulated by GFAP-IFN microglia were enriched for biological processes including response to cytokine, response to IFN, cell cycle, metabolism, response to virus, immune response, regulation of T cell cytotoxicity and antigen processing and presentation (Fig. 4d). Microglia from GFAP-IFN mice, but not GFAP-IL6 mice, had high levels of mRNA for MHC-I genes H2-D1, H2-K1, H2-Q4, H2-Q6, H2-Q7, H2-T10, H2-T22, H2-T23, as well as other genes associated with antigen processing, such as B2m, Tap1 and Tapbp (Fig. 4a, d). Enhanced MHC-I expression may be induced by expression of the NOD-like receptor Nlrc5, the master regulator of MHC-I gene expression [76], which was also upregulated exclusively by GFAP-IFN microglia. Interestingly, consistent with the presence of proliferating and apoptotic microglia in the brain of GFAP-IFN mice (Fig. 1), biological processes associated with both cell cycle transition and apoptosis were significantly enriched by IFN-α-skewed genes. IFN-α-skewed genes included proliferation promoters (Ccnd2, Ccnf, Cdc23, Ube2c and Kifc1), as well as proliferation inhibitors (Cdkn1a, Cdkn2c, Cdkn2d, E2f7 and E2f8) (Fig. 4a, d). Furthermore, GFAP-IFN microglia also expressed IFN-α-skewed genes that induce apoptosis (Casp4, Tnfsf10, Pidd1, Shisa5, Ifit2, Ifit3, Oas1h and Rnasel), as well as genes that promote survival (Axl, Adar and Apip).

Expression of selected core response, IL-6- and IFN-α-skewed genes, further analyzed by RTPCR of microglia from 1-, 3- and 6-month-old mice, for the most part, validated the RNA-seq data (Additional file 1: Fig. S4). Furthermore, the expression of some core response and IL-6- and IFN-α-skewed genes was further enhanced in older GFAP-IL6 and GFAP-IFN mice. In summary, microglia adopt a common transcriptional program and express overlapping target genes in response to the cytokine environments induced by both IL-6 and IFN-α, but also express stimulus-specific IL-6- or IFN-α-skewed genes.

Microglia in the brain of GFAP-IL6 versus GFAP-IFN mice have distinct surface marker profiles related to the biological processes identified by transcriptomic analysis

Having identified that microglia regulate unique genes associated with distinct functions in response to chronic IL-6 versus IFN-α production, such as phagocytosis and antigen presentation, we next examined whether these cells showed altered surface marker proteins associated with these cell functions. Since many of these functional markers are not specific to microglia, high-dimensional flow cytometry was performed on brain leukocytes from MacGreen-WT, MacGreen-GFAP-IL6 and MacGreen-GFAP-IFN mice and UMAP dimensionality reduction was first used to discern the global leukocyte landscape including microglia (Fig. 5, Additional file 1: Figs. S5, S6). UMAP analysis identified a large cluster of microglial cells as well as several leukocyte subpopulations present in the brain of the mice (Fig. 5a). Compared with WT, GFAP-IL6 mice had significantly increased numbers of conventional dendritic cells (cDCs) in the brain (Additional file 1: Fig. S6). On the other hand, compared with WT and GFAP-IL6 mice, there were increased numbers of Ly6Clow monocytes, Ly6Chigh monocytes, NK cells, CD4+ T cells and CD8+ T cells present in the brain parenchyma of GFAP-IFN mice at all ages studied. In addition, the number of dendritic cells increased in the brain of GFAP-IFN mice with age and compared with WT and GFAP-IL6 mice, the number of cDCs in the brain of GFAP-IFN mice was significantly increased at 3 months of age and the number of plasmacytoid dendritic cells (pDCs) was significantly increased at 3 and 6 months of age.

Fig. 5
figure 5

Distinct global leukocyte landscapes in the brain of GFAP-IL6 versus GFAP-IFN mice. a UMAP plot of entire dataset (brains of MacGreen-WT, -GFAP-IL6 and -GFAP-IFN mice at 1, 3 and 6 months of age) labeled with FlowSOM cluster identities. b UMAP plots of the 6-month-old dataset split into WT, GFAP-IL6 and GFAP-IFN mice. c UMAP plots of the same dataset colored by the expression of TMEM119, CD11b, CD16/32 or CD64. d Cluster overlay of CD16/32 versus CD64 levels on the surface of microglia from the brains of WT, GFAP-IL6 and GFAP-IFN mice at 1, 3 and 6 months of age. n = 3–5 mice/group

Dimensionality reduction separated microglia by genotype, as WT, GFAP-IL6 and GFAP-IFN microglia occupied largely distinct areas of the microglia cluster, while genotype-dependent divergent clustering was not observed for other cell types (Fig. 5b). This suggests that microglia elicit robust, homogenous responses to chronic IL-6 and IFN-α production which drive these cells to adopt divergent states. Separation of microglial cells from WT, GFAP-IL6 and GFAP-IFN mice was based on differing surface levels of TMEM119, CD11b, CD16/32 and CD64 proteins (Fig. 5c). TMEM119 levels were lower on GFAP-IL6 microglia compared with WT and GFAP-IFN microglia. By contrast, levels of CD11b and CD16/32 were elevated on the surface of GFAP-IL6 microglia compared with WT and GFAP-IFN microglia. Additionally, CD64 levels were increased on microglia from both GFAP-IL6 and GFAP-IFN mice, with levels highest on the surface of GFAP-IL6 microglia. Two-way analysis of CD16/32 and CD64 levels demonstrated that microglia had differential responses in GFAP-IL6 versus GFAP-IFN mice (Fig. 5d). Compared with 1-month-old WT microglia, GFAP-IL6 microglia were CD16/32hi CD64hi, while GFAP-IFN microglia had lower levels of CD16/32 than GFAP-IL6 microglia and were CD16/32int CD64hi. Microglia in the brains of 3- and 6-month-old GFAP-IL6 mice remained CD16/32hi CD64hi, however, cells in the brain of older GFAP-IFN mice had reduced levels of CD64 compared with 1-month-old GFAP-IFN mice.

To further examine the microglial cell surface profile, we next used the data from Fig. 5 and gated for microglia (live eGFP+ CD45low CD11b+ Ly6C TMEM119+ cells; Additional file 1: Fig. S7a) and quantified the median fluorescence intensity (MFI) of microglia signature and myeloid lineage markers (Fig. 6a–e), as well as the percentage of microglia that were positive for these markers (Fig. 6i-k, m–o). In a separate experiment, we gated microglia (live eGFP+ CD11b+ 4D4+ cells or live eGFP+ CD11b+ FCRLS+ cells; Additional file 1: Fig. S7b-c) and quantified the percentage of MHC-I+ cells and the MFI of 4D4, FCRLS and MHC-I (Fig. 6f-h, l). GFAP-IL6 microglia had reduced TMEM119 levels compared with both WT and GFAP-IFN mice, while GFAP-IFN microglia had elevated levels of TMEM119 compared with WT and GFAP-IL6 mice (Fig. 6a). Levels of microglia-specific 4D4 were slightly reduced on 1-month-old GFAP-IL6 and GFAP-IFN microglia compared with WT, however, at older ages, microglia from all three genotypes had comparable 4D4 levels (Fig. 6f). On the other hand, levels of the scavenger receptor FCRLS were significantly elevated on GFAP-IL6 and GFAP-IFN microglia compared with WT at all ages studied, however, these levels were highest on GFAP-IL6 microglia, being approximately 1.4-fold higher than that in GFAP-IFN mice (Fig. 6g).

Fig. 6
figure 6

Microglia in the brain of GFAP-IL6 versus GFAP-IFN mice have unique surface marker expression profiles. a–e From 1-, 3- and 6-month-old MacGreen-WT, -GFAP-IL6 and -GFAP-IFN mice, microglia were gated (live eGFP+ CD45low CD11b+ Ly6C TMEM119+ cells) and the median fluorescence intensity (MFI) was quantified for a TMEM119, b CD16/32, c CD64, d CD11b and e SCA-1 levels. f–h In a separate experiment, microglia were gated (live eGFP+ CD11b+ 4D4+ cells or live eGFP+ CD11b+ FCRLS+ cells) and the MFI was quantified for f 4D4, g FCRLS and h MHC-I levels. i-o Percentages of microglia positive for i CD64, j CD11c, k SCA-1, l MHC-I, m MHC-II, n CD80 and o CD86. n = 3–5 mice/group. Graphs show individual values per mouse and mean ± SEM. *, p < 0.05 compared with WT; ^, p < 0.05 compared with GFAP-IL6; #, p < 0.05 compared with 1-month-old of same genotype; x, p < 0.05 compared with 3-month-old of same genotype using two-way ANOVA with Tukey’s post-test. † N.B. The results for 6-month-old microglia are not shown due to high autofluorescence interference and are therefore not reliable

Next, we examined the levels of several myeloid functional markers on the surface of microglia from GFAP-IL6 and GFAP-IFN mice. As demonstrated above, compared with WT, both GFAP-IL6 and GFAP-IFN microglia had significantly elevated levels of CD16/32 at all ages studied (Fig. 6b). However, surface CD16/32 levels on GFAP-IL6 microglia were 2.3-fold higher than those from GFAP-IFN mice and continued to increase with age. Due to low levels, only half of the microglia from the brain of WT mice were detected as CD64+ (Fig. 6i). By contrast, surface CD64 levels were significantly elevated to comparable levels on GFAP-IL6 and GFAP-IFN microglia at 1 month of age (Fig. 6c) and almost all microglia in the brain of these mice were CD64+ positive (Fig. 6i). While CD64 levels on GFAP-IL6 microglia remained elevated at all ages studied, levels on GFAP-IFN microglia gradually decreased with age, such that microglia from 6-month-old GFAP-IFN mice had CD64 levels 1.5-fold lower than GFAP-IL6 microglia (Fig. 6c). At all ages examined, microglia in the brain of GFAP-IL6 mice had significantly higher levels of CD11b compared with WT and GFAP-IFN mice (Fig. 6d). While CD11c+ microglia were barely detected in the brain of WT mice (Fig. 6j), the brain of both GFAP-IL6 and GFAP-IFN mice had a similarly sized CD11c+ microglia subpopulation that accounted for approximately 7% of microglia at 1 month of age, 10% of microglia at 3 months of age and 15% of microglia at 6 months of age.

GFAP-IFN microglia had increased surface levels of markers associated with IFN-α signaling, antigen presentation and lymphocyte co-stimulation, such as stem cell antigen-1 (SCA-1), MHC-I, MHC-II, CD80 and CD86 (Fig. 6e, h, k, l, m–o). Consistent with a robust microglial cell IFN-I response, virtually all microglia in the brain of GFAP-IFN mice were SCA-1+ (Fig. 6k). Although almost all microglia from all three genotypes were MHC-I+ (Fig. 6l), MHC-I levels were elevated by GFAP-IL6 and GFAP-IFN microglia compared with WT microglia (Fig. 6h). At all ages examined, MHC-I levels were highest on the surface of GFAP-IFN microglia and were 2.3-fold higher than levels of MHC-I on the surface of GFAP-IL6 microglia. The brain of WT mice had small numbers of MHC-II+ microglia (Fig. 6m). While 0.3% and 0.9% of 1- and 3-month-old GFAP-IL6 microglia, respectively, were MHC-II+, this increase was not significant. By contrast, the MHC-II+ microglial cell population was significantly increased in GFAP-IFN mice compared with WT, with 1.2% and 3% of microglia positive for MHC-II at 1 and 3 months of age, respectively. Additionally, the number of MHC-II+ microglia in the brain of 3-month-old GFAP-IFN mice was 2.5-fold higher than in 1-month-old GFAP-IFN mice. Only 0.5% and 0.2% of microglia were CD80+ in the brains of 1-month-old WT and GFAP-IL6 mice, respectively (Fig. 6n). By contrast, the proportion of microglia that were CD80+ in the brain of GFAP-IFN mice was 8.9-fold and 22-fold higher, respectively, at 1 month of age. Compared with 1-month-old GFAP-IFN mice, 3-month-old GFAP-IFN mice had a CD80+ microglia population that was approximately 4.2-fold smaller. Similarly, while a very small number (0.2–1%) of WT and GFAP-IL6 microglia were CD86+, at 1 and 3 months of age, there was a significant increase (30- and 5-fold, respectively) in the number of CD86+ microglia in the brain of GFAP-IFN mice (Fig. 6o). In summary, our findings demonstrate that the global leukocyte landscape is divergent in the brain of WT, GFAP-IL6 and GFAP-IFN mice and microglia adopt distinct surface marker profiles in these animals.

IL-6- and IFN-α-like microglia responses are present in distinct neuropathological states

We next asked whether microglia from other neuropathological states exhibited responses that aligned with those observed in the GFAP-IL6 and GFAP-IFN mice. We selected microglia RNA-seq studies reported for a number of different neurological diseases in mice, including Alzheimer’s disease (AD) [52], tauopathy [53], experimental autoimmune encephalomyelitis (EAE) [54] and lipopolysaccharide (LPS) endotoxemia [55]. In some disease models, microglia were isolated into subsets based on expression of disease-associated markers on their cell surface, such as Clec7a in AD [52] and CD11c in EAE [54]. Unfortunately, at the time of analysis, there were no appropriate RNA-seq datasets from purified microglia from mice during viral infection that we could include for comparison in our meta-analysis. Hierarchical clustering separated the samples in our gene expression matrix into a large cluster of microglia samples from “control” groups and a large cluster of microglia samples from “disease” groups (Fig. 7a, Additional file 1: Fig. S8). Additionally, there were 1,759 genes differentially expressed in at least 4 different conditions analyzed, which were separated into 22 co-regulated clusters by hierarchical clustering (Additional file 1: Fig. S8, Additional file 4: Table S3).

Fig. 7
figure 7

Meta-analysis of microglia gene expression datasets identified core and cytokine-specific co-regulated gene clusters. a Heatmap of z-scores (within-study-normalized) for 1,759 genes differentially expressed in at least 4 comparisons. Hierarchical clustering identified 22 clusters of co-regulated genes. b Summary of gene set changes from the meta-analysis, which identified upregulated danger response genes (clusters 2 and 9), downregulated danger response genes (clusters 1, 8 and 22), chronic response genes (cluster 7), IL-6-response genes (cluster 19) and IFN-response genes (clusters 5, 12 and 14). Differential expression was calculated by comparing each condition with its respective control. The median log2-fold change for each gene cluster is shown

Using hierarchical clustering and visual inspection, we identified gene clusters 1, 2, 8, 9 and 22 that were dysregulated in almost every disease state and termed them "danger response genes". They included upregulated genes such as Axl, Bhlhe40, C1qa and Spp1 and downregulated genes such as Gpr34, Klk9, P2ry12 and Rnf216 (Fig. 7b, Additional file 1: Fig. S9a–d, Additional file 4: Table S3). While APP-PS1 Clec7a microglia, which are not associated with Aβ-plaques [52], did not differentially regulate these genes, clusters 2 and 9 were broadly upregulated and clusters 1, 8 and 22 were broadly downregulated by microglia in the other neurological disease states analyzed, including APP-PS1 Clec7a+, hMAPT-301S, EAE CD11c, EAE CD11c+ and LPS microglia. Interestingly, 50% of the genes identified by Friedman et al. in a microglia-specific cluster were identified as downregulated universal danger genes in our meta-analysis (Additional file 1: Fig. S10a, b) [53]. We also identified a cluster (cluster 7), which included 125 genes that were upregulated by microglia in response to chronic stimuli, such as neurodegeneration and autoimmunity, but were unchanged or decreased in microglia responding to acute LPS-induced endotoxemia; these genes were termed "chronic response genes" and included genes such as Apoe, Clec7a, Gas6, Itgax and Siglecf (Fig. 7b, Additional file 1: Fig. S9e-f, Additional 1: Table S3). While microglia not associated with pathology in APP-PS1 mice (APP-PS1 Clec7a microglia) did not upregulate the chronic response genes, these genes were upregulated by APP-PS1 Clec7a+, hMAPT-P301S, EAE CD11c and EAE CD11c+ microglia. While many of the danger response and chronic response genes were also dysregulated by GFAP-IL6 and GFAP-IFN microglia, fold-change analysis demonstrated that the distribution of gene expression was more variable than in microglia from other neuropathological conditions (Additional file 1: Fig. S9b, d, f). This is likely the result of the modest neuropathological features of disease in these mice at 1 month of age, compared with the more pronounced perturbation of microglia in mice at later stages of disease in the other neuropathological states. However, many of the danger response and chronic response genes upregulated by microglia in various different neuropathological states were classified by us as either core response genes commonly regulated by IL-6 and IFN-α, or IL-6- or IFN-α-skewed, with genes that are regulated by microglia in response to IL-6 and/or IFN-α comprising 36.1%, 35.8% and 30.4% of the danger response gene clusters 2 and 9, clusters 1, 8 and 22 and the chronic response gene cluster 7, respectively (Additional file 1: Fig. S10c–e).

In addition to core transcriptional programs, meta-analysis also identified cytokine-specific gene signatures that were only regulated in certain disease states (Fig. 7, Additional file 1: Fig. S11, Additional file 4: Table S3). Cluster 19, termed the "IL-6-response gene" cluster, consisted of 22 genes, including Ccnb2, C5ar1, Hs3st3b1 and Id2, which were highly upregulated in microglia from GFAP-IL6 mice and mice with LPS-induced endotoxemia (Fig. 7b, Additional file 1: Fig. S11a, e), with the median log2 expression of these genes upregulated 2.1- and 1.9-fold, respectively, compared with their respective controls. Similarly, these genes were modestly upregulated in both CD11c and CD11c+ microglia during EAE, with the median log2 expression of these genes upregulated 0.9- and 0.8-fold respectively. By contrast, in neurodegenerative conditions such as AD and tauopathy, microglial cell expression of the IL-6-response cluster genes were unchanged or only slightly increased.

Hierarchical clustering also identified three clusters, clusters 5, 12 and 14, which included a total of 179 genes; these genes were termed "IFN-response genes" and included Adar, H2-Q7, Ifit1, Irf7, Tnfsf10 and Gbp2 (Fig. 7b, Additional file 1: Fig. S11b-d, f–h, Additional file 4: Table S3). In addition to GFAP-IFN microglia, most of the IFN-response cluster genes were highly upregulated in both CD11c and CD11c+ microglia during EAE and in microglia during LPS-induced endotoxemia. On the other hand, regulation of the IFN-response genes by microglia from GFAP-IL6 mice and in neurodegenerative states, such as Clec7a+ cells from APP-PS1 mice and microglia from hMAPT-301S mice, was more variable, with only modest upregulation of IFN-response genes in some clusters and unchanged or downregulated expression in others. Taken together, in distinct neuropathological states, microglia uniquely regulate genes associated with IL-6 or IFN-α signaling, which are components of the universal danger response genes and, in chronic states of perturbation, the chronic microglia response.

Discussion

Here we examined the phenotypic changes of microglia in response to chronic IL-6 versus IFN-α signaling in the brain. In response to the respective cytokine milieu, microglia undergo distinct molecular and cellular adaptations as they fine-tune their phenotypes in the brain of GFAP-IL6 versus GFAP-IFN mice. Their adaptations correlated with the expression of unique IL-6- or IFN-α-skewed genes in addition to a core transcriptional program. Further, transcriptomic meta-analysis demonstrated that both IL-6 and IFN-α contribute to the formation of a core microglial cell transcriptional response in a wide range of neuropathological conditions. Finally, our findings suggest that IFN-α signaling has a more extensive effect on the microglial cell transcriptional landscape in disease compared with IL-6.

Chronic IL-6 and IFN-α production induces distinct manifestations of neuroinflammation in GFAP-IL6 and GFAP-IFN mice

The distinct activation of the signature transcription factors STAT3 and STAT1 in the brains of GFAP-IL6 and GFAP-IFN mice further highlighted the specific nature of the response by microglia to IL-6 and IFN-α, respectively, and is consistent with previous reports [17, 75]. Interestingly, serine phosphorylation of STAT1, which is required for the maximal activation of this transcription factor [72, 73], was barely detectable in GFAP-IL6 brains further indicating the minor role that STAT1 activation has in the CNS of these mice [75]. By contrast, low levels of STAT3 activation by IFN-α were also previously observed in the CNS of GFAP-IFN mice [75] and may play a role in the modulation of the IFN-I response [77]. The direct comparison of the activation of IL-6 and IFN-α signaling pathways in the CNS of GFAP-IL6 and GFAP-IFN mice demonstrates the remarkable specificity by which these pathways are activated. Although these cytokines are pleiotropic, the activation of STAT3 versus STAT1 is unique in the CNS of GFAP-IL6 and GFAP-IFN mice. Furthermore, the chronic production of IL-6 and IFN-α causes progressive, destructive diseases that are associated with the production of secondary mediators of inflammation in these mice [7, 8, 12, 15, 16, 28, 29], however, the specific activation of STAT3 versus STAT1 suggests that IL-6 and IFN-α induce distinct manifestations of neuroinflammation.

The dynamics of microglial cell proliferation and apoptosis are divergently altered by IL-6- versus IFN-α-induced neuroinflammation

Significantly increased numbers of microglia have been reported in the brain of adult GFAP-IL6 mice [14, 78,79,80,81]. Our findings here in juvenile and adult mice revealed that this increase was the consequence of significantly augmented proliferation of these cells in the brain while apoptosis was unchanged. This is consistent with reports that IL-6 can stimulate the proliferation of microglia in vitro [82] and in vivo [83, 84]. However, the total number of microglia in GFAP-IL6 mice, as well as the number of BrdU+ microglia, decreased from as early as 1 month of age correlating with IL-6 mRNA levels peaking in 1- and 3-month-old GFAP-IL6 mice before declining by 8 and 12 months of age [14], possibly due to astrocyte degeneration and loss [13]. The notion that chronic IL-6 signaling directly promotes microglia proliferation is also supported by our observation of greater numbers of proliferating microglia present in brain regions of GFAP-IL6 mice shown to have high levels of IL-6, such as the cerebellum, compared with lower numbers of proliferating cells in regions with lower levels of IL-6, such as the cortex and hippocampus [8, 15, 74]. Since apoptosis did not appear to substantially contribute to the turnover of microglia in older GFAP-IL6 mice, the mechanism by which microglial cell numbers decreased in the cerebellum of aged animals remains unknown but may be the result of microglial cell egress to other regions of the brain or due to activation of alternative cell death pathways.

In contrast to IL-6, IFN-α exhibits anti-proliferative effects in bone marrow-derived macrophages [85, 86] and elevated maternal levels of the type I IFN, IFN-β, contribute to arrested microglial cell proliferation in the CNS of newborn mice [87]. It was therefore surprising to observe increased numbers of BrdU+ microglia in all three brain regions of GFAP-IFN mice independent of age. In vesicular stomatitis virus-mediated encephalitis, astroglial and neuronal, but not microglial IFNAR signaling is required for the proliferation of microglia, suggesting that astrocytes and neurons respond to IFN-I and in turn regulate the ability of microglia to proliferate during viral encephalitis [88]. Thus, it is conceivable that in GFAP-IFN mice, astrocytes and neurons respond to high levels of IFN-α by producing factors which promote the proliferation of microglia. Importantly, TUNEL staining indicated that the increase in microglial cell number due to proliferation was offset by the loss of these cells due, at least in part, to increased apoptosis. Induction of apoptosis is a well-known action of IFN-I [89,90,91] and microglia from GFAP-IFN mice had significantly increased expression of a number of genes linked to the apoptosis-promoting function of IFN-α, including 2′-5′-oligoadenylate synthase (OAS), RNase L [92, 93], TRAIL (encoded by Tnfsf10) [94], interferon-induced protein with tetratricopeptide repeats (IFIT)2, IFIT3 [95] and the non-canonical inflammasome activator caspase-4 (encoded by Casp4) [96]. Thus, the cytokine environment induced by IFN-α in GFAP-IFN mice promotes both microglia proliferation and death, overall resulting in relatively stable cell numbers.

While the impact of microgliosis has been the subject of considerable attention, the functional outcomes of changes in microglial cell turnover in neurological diseases have largely been overlooked [97]. In mice intraperitoneally administered with LPS, there is transient, robust microglial cell proliferation in brain regions with substantial neuroinflammation [97]. Blockade of microglia proliferation in these animals prolonged the signs of LPS-induced sickness, demonstrating that the transient increase in microglia number is protective and these cells attenuate the sickness response [98]. Similarly, proliferating microglial cells are neuroprotective in cerebral ischemia and their selective ablation exacerbates neuroinflammation and ischemic damage [99]. On the other hand, the expansion of the microglial cell population during prion disease is detrimental and contributes to disease progression [100]. Together, these findings demonstrate that the regulation of the dynamics of microglia turnover is a critical component of the microglial cell response to neuroinflammation and can be helpful or harmful depending on the context.

Neuroinflammation in GFAP-IL6 versus GFAP-IFN mice induces unique morphological changes in microglia

In addition to cell numbers, in the brain of GFAP-IL6 versus GFAP-IFN mice, microglia exhibited stimulus-specific morphological changes. Microglia from GFAP-IL6 mice had shorter processes with less branching and complexity, similar to earlier reports in these animals of microglia with larger somata and processes that occupy a smaller area of the brain [101]. Importantly, the microglia in GFAP-IL6 mice had a similar morphology to microglia in the brains of patients with NMOSD [24] and animals with pathologies associated with elevated IL-6, such as acute LPS-induced endotoxemia [102, 103], stroke [104] and ischemic stroke and reperfusion [105]. On the other hand, microglia from GFAP-IFN mice were hypertrophied, had longer processes and were hyper-ramified, comparable to microglia in the brain of mice reported during aging [106, 107], chronic stress [108, 109], IFN-α-induced depression [110] and following chronic, CNS-targeted production of IFN-β [107]. Furthermore, in neurodegenerative diseases such as AD, concurrent activation of IFN-I signaling [111,112,113] and the production of cellular injury- and death-associated factors and other AD-related proinflammatory factors, including IL-1β, IL-6 and TNF [112, 114], may contribute to the heterogenous appearance of classically “reactive” [112, 115] or hypertrophied microglia [111, 115, 116]. These findings further demonstrate that activated microglia do not follow a linear program towards an amoeboid morphology, but rather, microglia exist in a morphological spectrum that depends on the nature, duration and context of the stimuli [34]. Since the morphologies of microglia in GFAP-IL6 versus GFAP-IFN mice were similar to those observed in human diseases associated with elevated IL-6 versus IFN-α levels, it is conceivable that the microglial cell morphological changes observed in these diseases in humans reflect direct microglia responses to these cytokines.

The distinct molecular signatures induced by chronic production of IL-6 versus IFN-α are associated with unique microglia activities

Paralleling the stimulus-specific differences in microglia turnover and morphology, we also observed stimulus-specific expression of genes and production of proteins in microglia from GFAP-IL6 versus GFAP-IFN mice. Correlating with progressive demyelination, proliferative angiopathy and iron accumulation observed in the brain of GFAP-IL6 mice [8, 13, 117], microglia expressed genes and proteins related to phagocytosis and processing of damaged myelin and other cellular detritus, as well as angiogenesis and the regulation and metabolism of iron. One of the most highly upregulated genes in GFAP-IL6 microglia, Apoe, is a core component of a neurodegenerative microglia phenotype (MGnD), also known as disease-associated microglia (DAM) transcriptional signature [52, 118]. Microglial cell secretion of APOE protein has been suggested to be a key player in the clearance of dead cells, extracellular debris and apoptotic neurons [119, 120], thereby enabling microglia to suppress damage to healthy neurons in neurodegenerative states [119]. However, APOE has been demonstrated to be an intrinsic regulator of a microglia phenotype associated with neurodegeneration and may contribute to disease [52, 121]. GFAP-IL6 microglia also uniquely upregulated Vegfa, which promotes angiogenesis and may contribute to the proliferative angiopathy that occurs in the brain of these animals [8]. VEGF-A also drives chemotaxis and proliferation of microglia and other cells [122], while LPS-activated microglia co-cultured with endothelial cells produce VEGF-A and promote angiogenesis [123]. Microglial cell expression of Vegfa in response to IL-6 may therefore induce both microglial cell proliferation and contribute to pathological angiogenesis observed in the brain of GFAP-IL6 mice [8]. Aberrant deposition of iron is a prominent feature of neurodegenerative diseases [124, 125] and progressively accumulates in the cerebellum of GFAP-IL6 mice, likely as a result of chronic leakage of the blood–brain barrier [117]. Increased microglial cell expression of IL-6-skewed genes associated with iron ion uptake and transport, such as Slc39a14 and transferrin (Trf) [126], may therefore impart these cells with the ability to sense and take up iron in the brain of GFAP-IL6 mice. However, since iron accumulates in the brain of older GFAP-IL6 mice, this function is likely impaired, is not sufficient to restrain iron accumulation as these animals age, or iron remains in the brain following uptake by microglia.

In contrast to GFAP-IL6 mice, GFAP-IFN mice have unique brain pathology. The microglia in these animals predominantly regulated genes and surface marker proteins that were associated with antiviral immunity, with enrichment for processes related to antigen processing and presentation and regulation of T cell cytotoxicity. Microglia from GFAP-IFN mice had significantly increased levels of MHC-I and MHC-II, as well as the T cell co-stimulation molecules CD80 and CD86. The upregulation of antigen presenting-related markers, together with the hyper-ramified morphology of GFAP-IFN microglia, are indicative of the acquisition of an antigen-presenting cell phenotype, likely allowing these cells to communicate with, interact with and present antigen to peripheral leukocytes that progressively infiltrate the brain of these mice. These findings overlap with the transcriptomic profile of primary microglia stimulated with IFN-α in vitro [127] and are also consistent with heightened antiviral activity in the brain of GFAP-IFN mice, as indicated by improved survival and reduced immune pathology in lymphocytic choriomeningitis virus and herpes simplex virus-1-infected animals [7, 128].

Microglia responses in neuroinflammatory and neurodegenerative diseases comprise IL-6- and IFN-α-like responses

Meta-analysis of the transcriptional landscapes of microglia in diverse neuropathological diseases including AD, tauopathy, EAE, LPS endotoxemia, as well as GFAP-IL6 and GFAP-IFN mice, identified a core microglial cell transcriptional response, termed danger response genes, which were similarly regulated in most of the disease states studied. These genes included Axl, Bhlhe40, Cd274, Ctsb, Ctss, Spp1 and others, which were previously identified as widely regulated by microglia in many different conditions [53, 57, 129]. The upregulation of these genes allow microglia to react to perturbation, enabling them to phagocytose damaged cells and other material (Axl, C1qa, C1qb) [130,131,132,133,134], interact and communicate with recruited lymphocytes (B2m, Cd274, H2-D1) [135, 136], process and metabolize myelin debris and other detritus (Ch25h, Lpl) [137, 138], produce secondary inflammatory mediators (Ccl3, Ctsb, Ctsc, Ctss, Trim25) [139,140,141,142,143,144] and other functions. Furthermore, Bhlhe40 was identified as a danger response gene and is suggested to be a putative transcription factor which regulates microglial cell activity in response to disturbances in their local environment [53]. In addition, danger genes universally downregulated by perturbed microglia included genes which distinguish microglia from myeloid cells in peripheral tissues. This is consistent with increasing evidence demonstrating that during neurological disease, microglia decrease expression of homeostatic signature genes and increase expression of neuroinflammatory genes in order to react to perturbation [34]. Chronic response genes that were upregulated in states of chronic neuroinflammation, but not in response to acute stimuli, included Apoe, Cd83, Chst2, Clec7a, Itgax and Itgb2 and were similar to clusters of “primed” or “neurodegenerative” genes in other meta-analysis studies [53, 57]. Further, the chronic response of microglia to nucleic acid-positive amyloid, but not acute exposure to IFN-β, induced the expression of Apoe and Clec7a [112]. Importantly, a considerable number of the danger and chronic response genes were classified by us as IL-6- and/or IFN-α-regulated and comprised 36.1%, 35.8% and 30.4% of the danger response gene clusters 2 and 9, clusters 1, 8 and 22 and the chronic response gene cluster 7, respectively. These findings indicate that both IL-6 and IFN-α signaling, likely in combination with other inflammatory factors, have an important role in the formation of the core microglial cell transcriptional response that arises as a result of a wide range of pathogenic stimuli and/or loss of homeostasis.

We also identified clusters of cytokine-specific genes that were regulated in specific disease states. In both CD11c and CD11c+ microglia during EAE and in microglia during acute LPS-induced endotoxemia, IL-6- and IFN-response genes were highly enriched. However, these genes were upregulated to a lesser degree by microglia in neurodegenerative mouse models, such as APP-PS1 and hMAPT-P301S mice. Microglia have been shown to upregulate IFN-response genes in a number of distinct neuropathological conditions, including AD and LPS endotoxemia [53, 112, 145]. However, our meta-analysis identified multiple IFN-response gene clusters. Cluster 5 genes included the classical IFN-regulated genes Adar, Ifit1, Irf7 and Stat1 and these were upregulated to varying degrees by microglia in every disease state studied. On the other hand, microglia from GFAP-IFN mice and mice with EAE expressed high levels of IFN-response cluster 12 and 14 genes, while microglia in LPS endotoxemia only expressed high levels of cluster 14 genes and expression of cluster 12 and 14 genes was not changed or was instead downregulated in neurodegenerative conditions. Although, to our knowledge, comparative analysis of the levels of IL-6 or IFN-α in the CNS in multiple sclerosis (MS)/EAE versus AD has not been performed in humans or mice, our findings are indicative of a stronger, more extensive response made by microglia to IL-6 and IFN-I, likely due to higher levels of these cytokines in the CNS, during autoimmunity and acute insult compared with neurodegenerative diseases. Furthermore, the identification of one small cluster of IL-6-regulated genes and three large clusters of IFN-regulated genes, in addition to the larger overlap of IFN-skewed genes with danger and chronic response genes compared with IL-6-skewed genes, suggested that microglia elicit a more pronounced, extensive response to IFN-α compared with IL-6. In line with this, GFAP-IFN microglia had a larger total number of regulated genes and a smaller overlap of commonly regulated genes compared with GFAP-IL6 cells. These findings suggest that IFN-α may induce a stronger, dominant transcriptional response in disease states where both IL-6 and IFN-α are present, which may conceal or dilute the transcriptional effects of IL-6. This phenomenon should be considered when interpreting transcriptional profiles from neuroinflammatory diseases in which both IL-6 and IFN-α are present at elevated levels. For example, transcriptional analyses ascribe IFN-α with an important role in NMOSD [146, 147] but may under-appreciate the critical role of IL-6 signaling in these diseases. A prominence of both of these cytokines in NMOSD is consistent with clinical evidence, where IFN-β therapy exacerbates disease severity and increases relapse rates [148,149,150,151,152], in contrast to patients with MS where IFN-β ameliorates disease [153]. On the other hand, treatment with the IL-6R antibody tocilizumab significantly ameliorates disease in patients with NMOSD [148, 154,155,156,157].

Limitations and future directions

Redundancy in the downstream effects of different cytokines must be considered when interpreting the results of this study. For example, chronic production of IL-6, IFN-α or TNF in the CNS of mice all promote a neuroinflammatory state which predisposes these animals to seizures and neurodegeneration [2, 158, 159]. Furthermore, transgenic production of IL-6 or IFN-α in the CNS is also associated with production of secondary mediators which are likely to contribute to the effects observed in microglia. Increased expression of genes which encode for IL-1β, TNF and other inflammatory cytokines is detected in the CNS of GFAP-IL6 [14] and GFAP-IFN mice (Phillip West, unpublished data). Although these caveats must be considered when interpreting the effects of these cytokines, the cytokine environments induced by IL-6 versus IFN-α are distinct. Accordingly, while GFAP-IL6, GFAP-IFN and other transgenic lines including GFAP-TNF mice exhibit clinical and molecular phenotypes with certain overlapping features, the differences are more prominent and can be attributed to the unique actions of each cytokine [160, 161]. Furthermore, in this study we observed stark differences in the signaling pathways activated by IL-6 versus IFN-α, as well as divergent alterations in the microglia molecular and cellular phenotype. However, in order to attribute the microglia responses to the direct actions of IL-6 versus IFN-α, loss-of-function experiments are required and are the focus of an ongoing study. Finally, although we did not determine the actual cytokine levels in the CNS of GFAP-IL6 and GFAP-IFN mice, the levels of cerebellar mRNA for transgenic IL-6 versus IFN-α were comparable (Phillip West, unpublished data).

The pharmacological ablation of microglia is an emerging strategy used to discern the functional significance of these cells in health and disease [2]. Accordingly, microglia ablation has enabled researchers to elucidate important actions of microglia, including maintenance of circadian rhythm [162], regulation of satiety and metabolism [163], and restraint of neuronal hyperexcitability [164]. In addition, the ablation of microglia is protective in certain neuropathological contexts but is harmful in others [165, 166], suggesting diverse functions of these cells in disease. These findings highlight the value of depleting microglia to dissect their functional significance in GFAP-IL6 and GFAP-IFN mice and this has been the approach of a further study (manuscript in preparation).

Furthermore, the production of GFAP by astrocytes is highly variable and is universally regulated in response to perturbation of the CNS [167]. Astrogliosis is readily apparent in the brain of GFAP-IL6 and GFAP-IFN mice and these cells produce increased levels of GFAP [8, 28]. Since the reactivity of astrocytes is potentially associated with the differential regulation of cytokine production in GFAP-IL6 and GFAP-IFN mice, the activation of astrocytes should also be considered when interpreting the functional state of microglia in GFAP-IL6 and GFAP-IFN mice. Detailed analysis of the molecular and cellular phenotype of astrocytes in these animals is warranted in the future and may potentially provide insight into their interaction with microglia, as well as their regulation of the microglia phenotype, during disease progression.

Conclusions

Taken together, we demonstrate here that microglia are a major target and effector cell of IL-6 and IFN-α in vivo and respond to chronic production of these cytokines in the brain by making wide-ranging cellular, molecular and transcriptional adaptations. Importantly, microglia showed stimulus-specific responses, giving rise to unique and distinct microglia phenotypes. These findings further highlight the exquisite responsiveness of microglia to altered cytokine signaling and pathology in the CNS. In addition, the identified IL-6- and IFN-response genes here are also substantial components of the transcriptional response of microglia reported in a wide range of neuroinflammatory and neurodegenerative disorders, indicating that both cytokines likely contribute to the phenotype and function of microglia in these diseases. The stimulus-specific responses of microglia reflect known functional and pathological roles of IL-6 versus IFN-α in the brain and likely impart these cells with their unique functions in IL-6 versus IFN-α-mediated neurological disease.

Availability of data and materials

The authors declare that all data supporting the findings of this study are available within the article and its Additional Information files. The raw RNA-seq data have been deposited in the European Nucleotide Archive under the entry code PRJEB46781.

Abbreviations

AD:

Alzheimer’s disease

ANOVA:

Analysis of variance

BPs:

Biological processes

BrdU:

5-Bromo-2’-deoxyuridine

CNS:

Central nervous system

DEGs:

Differentially expressed genes

EAE:

Experimental autoimmune encephalomyelitis

eGFP:

Enhanced green fluorescent protein

FACS:

Fluorescence-activated cell sorting

FDR:

False discovery rate

FlowSOM:

Flow Self-Organizing Maps

FPKM:

Fragments per kilobase of transcript per million mapped reads

GFAP:

Glial fibrillary acidic protein

GO:

Gene ontology

IFN:

Interferon

IL:

Interleukin

LPS:

Lipopolysaccharide

MS:

Multiple sclerosis

NMOSD:

Neuromyelitis optica spectrum disorder

PCA:

Principal component analysis

pS:

Phospho-serine

pY:

Phospho-tyrosine

STAT:

Signal transducer and activator of transcription

TUNEL:

Terminal deoxynucleotidyl transferase dUTP nick end labeling

UMAP:

Uniform Manifold Approximation and Projection

WT:

Wild type

References

  1. Hofer MJ, Campbell IL. Immunoinflammatory diseases of the central nervous system - the tale of two cytokines. Br J Pharmacol. 2016;173(4):716–28.

    Article  CAS  PubMed  Google Scholar 

  2. West PK, Viengkhou B, Campbell IL, Hofer MJ. Microglia responses to interleukin-6 and type I interferons in neuroinflammatory disease. Glia. 2019;67(10):1821–41.

    PubMed  Google Scholar 

  3. Wallach AI, Tremblay M, Kister I. Advances in the treatment of neuromyelitis optica spectrum disorder. Neurol Clin. 2021;39(1):35–49.

    Article  PubMed  Google Scholar 

  4. Crow YJ, Neven B, Frémond ML. JAK inhibition in the type I interferonopathies. J Allergy Clin Immunol. 2021;148(4):991–3.

    Article  CAS  PubMed  Google Scholar 

  5. Ma X, Kermode AG, Hu X, Qiu W. NMOSD acute attack: Understanding, treatment and innovative treatment prospect. J Neuroimmunol. 2020;348:577387.

    Article  CAS  PubMed  Google Scholar 

  6. Crow YJ, Shetty J, Livingston JH. Treatments in Aicardi-Goutières syndrome. Dev Med Child Neurol. 2020;62(1):42–7.

    Article  PubMed  Google Scholar 

  7. Akwa Y, Hassett DE, Eloranta ML, Sandberg K, Masliah E, Powell H, et al. Transgenic expression of IFN-alpha in the central nervous system of mice protects against lethal neurotropic viral infection but induces inflammation and neurodegeneration. J Immunol. 1998;161(9):5016–26.

    CAS  PubMed  Google Scholar 

  8. Campbell IL, Abraham CR, Masliah E, Kemper P, Inglis JD, Oldstone MB, et al. Neurologic disease induced in transgenic mice by cerebral overexpression of interleukin 6. Proc Natl Acad Sci U S A. 1993;90(21):10061–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Bellinger FP, Madamba SG, Campbell IL, Siggins GR. Reduced long-term potentiation in the dentate gyrus of transgenic mice with cerebral overexpression of interleukin-6. Neurosci Lett. 1995;198(2):95–8.

    Article  CAS  PubMed  Google Scholar 

  10. Steffensen SC, Campbell IL, Henriksen SJ. Site-specific hippocampal pathophysiology due to cerebral overexpression of interleukin-6 in transgenic mice. Brain Res. 1994;652(1):149–53.

    Article  CAS  PubMed  Google Scholar 

  11. Samland H, Huitron-Resendiz S, Masliah E, Criado J, Henriksen SJ, Campbell IL. Profound increase in sensitivity to glutamatergic- but not cholinergic agonist-induced seizures in transgenic mice with astrocyte production of IL-6. J Neurosci Res. 2003;73(2):176–87.

    Article  CAS  PubMed  Google Scholar 

  12. Heyser CJ, Masliah E, Samimi A, Campbell IL, Gold LH. Progressive decline in avoidance learning paralleled by inflammatory neurodegeneration in transgenic mice expressing interleukin 6 in the brain. Proc Natl Acad Sci U S A. 1997;94(4):1500–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Brett FM, Mizisin AP, Powell HC, Campbell IL. Evolution of neuropathologic abnormalities associated with blood-brain barrier breakdown in transgenic mice expressing interleukin-6 in astrocytes. J Neuropathol Exp Neurol. 1995;54(6):766–75.

    Article  CAS  PubMed  Google Scholar 

  14. Chiang CS, Stalder A, Samimi A, Campbell IL. Reactive gliosis as a consequence of interleukin-6 expression in the brain: studies in transgenic mice. Dev Neurosci. 1994;16(3–4):212–21.

    Article  CAS  PubMed  Google Scholar 

  15. Barnum SR, Jones JL, Muller-Ladner U, Samimi A, Campbell IL. Chronic complement C3 gene expression in the CNS of transgenic mice with astrocyte-targeted interleukin-6 expression. Glia. 1996;18(2):107–17.

    Article  CAS  PubMed  Google Scholar 

  16. Carrasco J, Hernandez J, Gonzalez B, Campbell IL, Hidalgo J. Localization of metallothionein-I and -III expression in the CNS of transgenic mice with astrocyte-targeted expression of interleukin 6. Exp Neurol. 1998;153(2):184–94.

    Article  CAS  PubMed  Google Scholar 

  17. Campbell IL, Erta M, Lim SL, Frausto R, May U, Rose-John S, et al. Trans-signaling is a dominant mechanism for the pathogenic actions of interleukin-6 in the brain. J Neurosci. 2014;34(7):2503–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Uzawa A, Mori M, Arai K, Sato Y, Hayakawa S, Masuda S, et al. Cytokine and chemokine profiles in neuromyelitis optica: significance of interleukin-6. Mult Scler. 2010;16(12):1443–52.

    Article  CAS  PubMed  Google Scholar 

  19. Uzawa A, Mori M, Masuda H, Ohtani R, Uchida T, Sawai S, et al. Interleukin-6 analysis of 572 consecutive CSF samples from neurological disorders: A special focus on neuromyelitis optica. Clin Chim Acta. 2017;469:144–9.

    Article  CAS  PubMed  Google Scholar 

  20. Uzawa A, Mori M, Sawai S, Masuda S, Muto M, Uchida T, et al. Cerebrospinal fluid interleukin-6 and glial fibrillary acidic protein levels are increased during initial neuromyelitis optica attacks. Clin Chim Acta. 2013;421:181–3.

    Article  CAS  PubMed  Google Scholar 

  21. Matsushita T, Tateishi T, Isobe N, Yonekawa T, Yamasaki R, Matsuse D, et al. Characteristic cerebrospinal fluid cytokine/chemokine profiles in neuromyelitis optica, relapsing remitting or primary progressive multiple sclerosis. PLoS ONE. 2013;8(4):e61835.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Popescu BF, Lennon VA, Parisi JE, Howe CL, Weigand SD, Cabrera-Gomez JA, et al. Neuromyelitis optica unique area postrema lesions: nausea, vomiting, and pathogenic implications. Neurology. 2011;76(14):1229–37.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Popescu BF, Lucchinetti CF. Pathology of demyelinating diseases. Annu Rev Pathol. 2012;7:185–217.

    Article  CAS  PubMed  Google Scholar 

  24. Saji E, Arakawa M, Yanagawa K, Toyoshima Y, Yokoseki A, Okamoto K, et al. Cognitive impairment and cortical degeneration in neuromyelitis optica. Ann Neurol. 2013;73(1):65–76.

    Article  PubMed  Google Scholar 

  25. Guo Y, Weigand SD, Popescu BF, Lennon VA, Parisi JE, Pittock SJ, et al. Pathogenic implications of cerebrospinal fluid barrier pathology in neuromyelitis optica. Acta Neuropathol. 2017;133(4):597–612.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. McKeon A, Lennon VA, Lotze T, Tenenbaum S, Ness JM, Rensel M, et al. CNS aquaporin-4 autoimmunity in children. Neurology. 2008;71(2):93–100.

    Article  CAS  PubMed  Google Scholar 

  27. Chan KH, Tse CT, Chung CP, Lee RL, Kwan JS, Ho PW, et al. Brain involvement in neuromyelitis optica spectrum disorders. Arch Neurol. 2011;68(11):1432–9.

    Article  PubMed  Google Scholar 

  28. Campbell IL, Krucker T, Steffensen S, Akwa Y, Powell HC, Lane T, et al. Structural and functional neuropathology in transgenic mice with CNS expression of IFN-alpha. Brain Res. 1999;835(1):46–61.

    Article  CAS  PubMed  Google Scholar 

  29. Kavanagh D, McGlasson S, Jury A, Williams J, Scolding N, Bellamy C, et al. Type I interferon causes thrombotic microangiopathy by a dose-dependent toxic effect on the microvasculature. Blood. 2016;128(24):2824–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Crow YJ, Chase DS, Lowenstein Schmidt J, Szynkiewicz M, Forte GM, Gornall HL, et al. Characterization of human disease phenotypes associated with mutations in TREX1, RNASEH2A, RNASEH2B, RNASEH2C, SAMHD1, ADAR, and IFIH1. Am J Med Genet A. 2015;167(2):296–312.

    Article  CAS  Google Scholar 

  31. Aicardi J, Goutieres F. A progressive familial encephalopathy in infancy with calcifications of the basal ganglia and chronic cerebrospinal fluid lymphocytosis. Ann Neurol. 1984;15(1):49–54.

    Article  CAS  PubMed  Google Scholar 

  32. Klok MD, Bakels HS, Postma NL, van Spaendonk RM, van der Knaap MS, Bugiani M. Interferon-alpha and the calcifying microangiopathy in Aicardi-Goutieres syndrome. Ann Clin Transl Neurol. 2015;2(7):774–9.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Goutieres F, Aicardi J, Barth PG, Lebon P. Aicardi-Goutieres syndrome: an update and results of interferon-alpha studies. Ann Neurol. 1998;44(6):900–7.

    Article  CAS  PubMed  Google Scholar 

  34. Dubbelaar ML, Kracht L, Eggen BJL, Boddeke E. The kaleidoscope of microglial phenotypes. Front Immunol. 2018;9:1753.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Cho CE, Damle SS, Wancewicz EV, Mukhopadhyay S, Hart CE, Mazur C, et al. A modular analysis of microglia gene expression, insights into the aged phenotype. BMC Genomics. 2019;20(1):164.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Stratoulias V, Venero JL, Tremblay ME, Joseph B. Microglial subtypes: diversity within the microglial community. Embo j. 2019;38(17):e101997.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Tay TL, Bechade C, D’Andrea I, St-Pierre MK, Henry MS, Roumier A, et al. Microglia gone rogue: impacts on psychiatric disorders across the lifespan. Front Mol Neurosci. 2018;10:421.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Lier J, Streit WJ, Bechmann I. Beyond activation: characterizing microglial functional phenotypes. Cells. 2021;10:9.

    Article  CAS  Google Scholar 

  39. Sasmono RT, Oceandy D, Pollard JW, Tong W, Pavli P, Wainwright BJ, et al. A macrophage colony-stimulating factor receptor-green fluorescent protein transgene is expressed throughout the mononuclear phagocyte system of the mouse. Blood. 2003;101(3):1155–63.

    Article  CAS  PubMed  Google Scholar 

  40. Butovsky O, Siddiqui S, Gabriely G, Lanser AJ, Dake B, Murugaiyan G, et al. Modulating inflammatory monocytes with a unique microRNA gene signature ameliorates murine ALS. J Clin Invest. 2012;122(9):3063–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Guennewig B, Bitar M, Obiorah I, Hanks J, O’Brien EA, Kaczorowski DC, et al. THC exposure of human iPSC neurons impacts genes associated with neuropsychiatric disorders. Transl Psychiatry. 2018;8(1):89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Guennewig B, Lim J, Marshall L, McCorkindale AN, Paasila PJ, Patrick E, et al. Defining early changes in Alzheimer’s disease from RNA sequencing of brain regions differentially affected by pathology. Sci Rep. 2021;11(1):4865.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Tsatsanis A, McCorkindale AN, Wong BX, Patrick E, Ryan TM, Evans RW, et al. The acute phase protein lactoferrin is a key feature of Alzheimer’s disease and predictor of Aβ burden through induction of APP amyloidogenic processing. Mol Psychiatry. 2021;8:90.

    Google Scholar 

  44. Andrews S. FastQC: a quality control tool for high throughput sequence data. 2015. http://www.bioinformatics.babraham.ac.uk/projects/fastqc.

  45. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.

    Article  CAS  PubMed  Google Scholar 

  46. Pertea M, Pertea GM, Antonescu CM, Chang TC, Mendell JT, Salzberg SL. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015;33(3):290–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. R Core Team. R: A language and environment for statistical computing 2017. https://www.R-project.org/.

  49. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.

    Article  CAS  PubMed  Google Scholar 

  50. Liao Y, Wang J, Jaehnig EJ, Shi Z, Zhang B. WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res. 2019;47(W1):W199-w205.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Merico D, Isserlin R, Stueker O, Emili A, Bader GD. Enrichment map: a network-based method for gene-set enrichment visualization and interpretation. PLoS ONE. 2010;5(11):e13984.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Krasemann S, Madore C, Cialic R, Baufeld C, Calcagno N, El Fatimy R, et al. The TREM2-APOE pathway drives the transcriptional phenotype of dysfunctional microglia in neurodegenerative diseases. Immunity. 2017;47(3):566-81e9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Friedman BA, Srinivasan K, Ayalon G, Meilandt WJ, Lin H, Huntley MA, et al. Diverse brain myeloid expression profiles reveal distinct microglial activation states and aspects of Alzheimer’s disease not evident in mouse models. Cell Rep. 2018;22(3):832–47.

    Article  CAS  PubMed  Google Scholar 

  54. Wlodarczyk A, Holtman IR, Krueger M, Yogev N, Bruttger J, Khorooshi R, et al. A novel microglial subset plays a key role in myelinogenesis in developing brain. Embo j. 2017;36(22):3292–308.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Srinivasan K, Friedman BA, Larson JL, Lauffer BE, Goldstein LD, Appling LL, et al. Untangling the brain’s neuroinflammatory and neurodegenerative transcriptional responses. Nat Commun. 2016;7:11295.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Kolde R. pheatmap: pretty heatmaps 2019. http://www.cran.r-project.org/package=pheatmap/.

  57. Holtman IR, Raj DD, Miller JA, Schaafsma W, Yin Z, Brouwer N, et al. Induction of a common microglia gene expression signature by aging and neurodegenerative conditions: a co-expression meta-analysis. Acta Neuropathol Commun. 2015;3:31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Bauer T, Zagorska A, Jurkin J, Yasmin N, Koffel R, Richter S, et al. Identification of Axl as a downstream effector of TGF-beta1 during Langerhans cell differentiation and epidermal homeostasis. J Exp Med. 2012;209(11):2033–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Besnard V, Dagher R, Madjer T, Joannes A, Jaillet M, Kolb M, et al. Identification of periplakin as a major regulator of lung injury and repair in mice. JCI Insight. 2018;3:5.

    Article  Google Scholar 

  60. Xu Y, Schnorrer P, Proietto A, Kowalski G, Febbraio MA, Acha-Orbea H, et al. IL-10 controls cystatin C synthesis and blood concentration in response to inflammation through regulation of IFN regulatory factor 8 expression. J Immunol. 2011;186(6):3666–73.

    Article  CAS  PubMed  Google Scholar 

  61. Moon HY, Becke A, Berron D, Becker B, Sah N, Benoni G, et al. Running-induced systemic cathepsin B secretion is associated with memory function. Cell Metab. 2016;24(2):332–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Fierro-Fernandez M, Miguel V, Marquez-Exposito L, Nuevo-Tapioles C, Herrero JI, Blanco-Ruiz E, et al. MiR-9-5p protects from kidney fibrosis by metabolic reprogramming. Faseb j. 2020;34(1):410–31.

    Article  CAS  PubMed  Google Scholar 

  63. Grabert K, Michoel T, Karavolos MH, Clohisey S, Baillie JK, Stevens MP, et al. Microglial brain region-dependent diversity and selective regional sensitivities to aging. Nat Neurosci. 2016;19(3):504–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Mercadante CJ, Prajapati M, Conboy HL, Dash ME, Herrera C, Pettiglio MA, et al. Manganese transporter Slc30a10 controls physiological manganese excretion and toxicity. J Clin Invest. 2019;129(12):5442–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Artuz CM, Knights AJ, Funnell APW, Gonda TJ, Ravid K, Pearson RCM, et al. Partial reprogramming of heterologous cells by defined factors to generate megakaryocyte lineage-restricted biomolecules. Biotechnol Rep (Amst). 2018;20:e00285.

    Article  Google Scholar 

  66. Butovsky O, Jedrychowski MP, Moore CS, Cialic R, Lanser AJ, Gabriely G, et al. Identification of a unique TGF-beta-dependent molecular and functional signature in microglia. Nat Neurosci. 2014;17(1):131–43.

    Article  CAS  PubMed  Google Scholar 

  67. Ashhurst TM, Marsh-Wakefield F, Putri GH, Spiteri AG, Shinko D, Read MN, et al. Integration, exploration, and analysis of high-dimensional single-cell cytometry data using Spectre. Cytometry A. 2021;9:78.

    Google Scholar 

  68. Van Gassen S, Callebaut B, Van Helden MJ, Lambrecht BN, Demeester P, Dhaene T, et al. FlowSOM: using self-organizing maps for visualization and interpretation of cytometry data. Cytometry A. 2015;87(7):636–45.

    Article  PubMed  Google Scholar 

  69. McInnes L, Healy J, Melville J. UMAP: Uniform Manifold Approximation and Projection for dimension reduction. arXiv. 2018:1802.03426.

  70. Heinrich PC, Behrmann I, Muller-Newen G, Schaper F, Graeve L. Interleukin-6-type cytokine signalling through the gp130/Jak/STAT pathway. Biochem J. 1998;334(Pt 2):297–314.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Platanias LC. Mechanisms of type-I- and type-II-interferon-mediated signalling. Nat Rev Immunol. 2005;5(5):375–86.

    Article  CAS  PubMed  Google Scholar 

  72. Wen Z, Zhong Z, Darnell JE Jr. Maximal activation of transcription by Stat1 and Stat3 requires both tyrosine and serine phosphorylation. Cell. 1995;82(2):241–50.

    Article  CAS  PubMed  Google Scholar 

  73. Pilz A, Ramsauer K, Heidari H, Leitges M, Kovarik P, Decker T. Phosphorylation of the Stat1 transactivating domain is required for the response to type I interferons. EMBO Rep. 2003;4(4):368–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Quintana A, Muller M, Frausto RF, Ramos R, Getts DR, Sanz E, et al. Site-specific production of IL-6 in the central nervous system retargets and enhances the inflammatory response in experimental autoimmune encephalomyelitis. J Immunol. 2009;183(3):2079–88.

    Article  CAS  PubMed  Google Scholar 

  75. Sanz E, Hofer MJ, Unzeta M, Campbell IL. Minimal role for STAT1 in interleukin-6 signaling and actions in the murine brain. Glia. 2008;56(2):190–9.

    Article  PubMed  Google Scholar 

  76. Vijayan S, Sidiq T, Yousuf S, van den Elsen PJ, Kobayashi KS. Class I transactivator, NLRC5: a central player in the MHC class I pathway and cancer immune surveillance. Immunogenetics. 2019;71(3):273–82.

    Article  CAS  PubMed  Google Scholar 

  77. Tsai MH, Pai LM, Lee CK. Fine-tuning of type I interferon response by STAT3. Front Immunol. 2019;10:1448.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Gyengesi E, Liang H, Millington C, Sonego S, Sirijovski D, Gunawardena D, et al. Investigation into the effects of tenilsetam on markers of neuroinflammation in GFAP-IL6 mice. Pharm Res. 2018;35(1):22.

    Article  CAS  PubMed  Google Scholar 

  79. Gyengesi E, Rangel A, Ullah F, Liang H, Niedermayer G, Asgarov R, et al. Chronic microglial activation in the GFAP-IL6 mouse contributes to age-dependent cerebellar volume loss and impairment in motor function. Front Neurosci. 2019;13:303.

    Article  PubMed  PubMed Central  Google Scholar 

  80. Recasens M, Almolda B, Pérez-Clausell J, Campbell IL, González B, Castellano B. Chronic exposure to IL-6 induces a desensitized phenotype of the microglia. J Neuroinflammation. 2021;18(1):31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Almolda B, Villacampa N, Manders P, Hidalgo J, Campbell IL, Gonzalez B, et al. Effects of astrocyte-targeted production of interleukin-6 in the mouse on the host response to nerve injury. Glia. 2014;62(7):1142–61.

    Article  PubMed  Google Scholar 

  82. Streit WJ, Hurley SD, McGraw TS, Semple-Rowland SL. Comparative evaluation of cytokine profiles and reactive gliosis supports a critical role for interleukin-6 in neuron-glia signaling during regeneration. J Neurosci Res. 2000;61(1):10–20.

    Article  CAS  PubMed  Google Scholar 

  83. Fattori E, Lazzaro D, Musiani P, Modesti A, Alonzi T, Ciliberto G. IL-6 expression in neurons of transgenic mice causes reactive astrocytosis and increase in ramified microglial cells but no neuronal damage. Eur J Neurosci. 1995;7(12):2441–9.

    Article  CAS  PubMed  Google Scholar 

  84. Savarin C, Hinton DR, Valentin-Torres A, Chen Z, Trapp BD, Bergmann CC, et al. Astrocyte response to IFN-γ limits IL-6-mediated microglia activation and progressive autoimmune encephalomyelitis. J Neuroinflammation. 2015;12:79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Hertzog PJ, Hwang SY, Kola I. Role of interferons in the regulation of cell proliferation, differentiation, and development. Mol Reprod Dev. 1994;39(2):226–32.

    Article  CAS  PubMed  Google Scholar 

  86. Hwang SY, Hertzog PJ, Holland KA, Sumarsono SH, Tymms MJ, Hamilton JA, et al. A null mutation in the gene encoding a type I interferon receptor component eliminates antiproliferative and antiviral responses to interferons alpha and beta and alters macrophage responses. Proc Natl Acad Sci U S A. 1995;92(24):11284–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Ben-Yehuda H, Matcovitch-Natan O, Kertser A, Spinrad A, Prinz M, Amit I, et al. Maternal type-I interferon signaling adversely affects the microglia and the behavior of the offspring accompanied by increased sensitivity to stress. Mol Psychiatry. 2020;25(5):1050–67.

    Article  CAS  PubMed  Google Scholar 

  88. Chhatbar C, Detje CN, Grabski E, Borst K, Spanier J, Ghita L, et al. Type I interferon receptor signaling of neurons and astrocytes regulates microglia activation during viral encephalitis. Cell Rep. 2018;25(1):118-29.e4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Swiecki M, Wang Y, Vermi W, Gilfillan S, Schreiber RD, Colonna M. Type I interferon negatively controls plasmacytoid dendritic cell numbers in vivo. J Exp Med. 2011;208(12):2367–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Cheng L, Yu H, Li G, Li F, Ma J, Li J, et al. Type I interferons suppress viral replication but contribute to T cell depletion and dysfunction during chronic HIV-1 infection. JCI Insight. 2017;2:12.

    Article  Google Scholar 

  91. van Boxel-Dezaire AH, Zula JA, Xu Y, Ransohoff RM, Jacobberger JW, Stark GR. Major differences in the responses of primary human leukocyte subsets to IFN-beta. J Immunol. 2010;185(10):5888–99.

    Article  PubMed  Google Scholar 

  92. Li G, Xiang Y, Sabapathy K, Silverman RH. An apoptotic signaling pathway in the interferon antiviral response mediated by RNase L and c-Jun NH2-terminal kinase. J Biol Chem. 2004;279(2):1123–31.

    Article  CAS  PubMed  Google Scholar 

  93. Rusch L, Zhou A, Silverman RH. Caspase-dependent apoptosis by 2’,5’-oligoadenylate activation of RNase L is enhanced by IFN-beta. J Interferon Cytokine Res. 2000;20(12):1091–100.

    Article  CAS  PubMed  Google Scholar 

  94. Peteranderl C, Herold S. The impact of the interferon/TNF-related apoptosis-inducing ligand signaling axis on disease progression in respiratory viral infection and beyond. Front Immunol. 2017;8:313.

    PubMed  PubMed Central  Google Scholar 

  95. Mears HV, Sweeney TR. Better together: the role of IFIT protein-protein interactions in the antiviral response. J Gen Virol. 2018;99(11):1463–77.

    Article  CAS  PubMed  Google Scholar 

  96. Casson CN, Yu J, Reyes VM, Taschuk FO, Yadav A, Copenhaver AM, et al. Human caspase-4 mediates noncanonical inflammasome activation against gram-negative bacterial pathogens. Proc Natl Acad Sci U S A. 2015;112(21):6688–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  97. Furube E, Kawai S, Inagaki H, Takagi S, Miyata S. Brain region-dependent heterogeneity and dose-dependent difference in transient microglia population increase during lipopolysaccharide-induced inflammation. Sci Rep. 2018;8(1):2203.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  98. Torii K, Takagi S, Yoshimura R, Miyata S. Microglial proliferation attenuates sickness responses in adult mice during endotoxin-induced inflammation. J Neuroimmunol. 2022;365:577832.

    Article  CAS  PubMed  Google Scholar 

  99. Lalancette-Hébert M, Gowing G, Simard A, Weng YC, Kriz J. Selective ablation of proliferating microglial cells exacerbates ischemic injury in the brain. J Neurosci. 2007;27(10):2596–605.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  100. Gómez-Nicola D, Fransen NL, Suzzi S, Perry VH. Regulation of microglial proliferation during chronic neurodegeneration. J Neurosci. 2013;33(6):2481–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  101. Ullah F, Asgarov R, Venigalla M, Liang H, Niedermayer G, Münch G, et al. Effects of a solid lipid curcumin particle formulation on chronic activation of microglia and astroglia in the GFAP-IL6 mouse model. Sci Rep. 2020;10(1):2365.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  102. Hanamsagar R, Alter MD, Block CS, Sullivan H, Bolton JL, Bilbo SD. Generation of a microglial developmental index in mice and in humans reveals a sex difference in maturation and immune reactivity. Glia. 2017;65(9):1504–20.

    Article  PubMed  PubMed Central  Google Scholar 

  103. Madore C, Joffre C, Delpech JC, De Smedt-Peyrusse V, Aubert A, Coste L, et al. Early morphofunctional plasticity of microglia in response to acute lipopolysaccharide. Brain Behav Immun. 2013;34:151–8.

    Article  CAS  PubMed  Google Scholar 

  104. Kluge MG, Kracht L, Abdolhoseini M, Ong LK, Johnson SJ, Nilsson M, et al. Impaired microglia process dynamics post-stroke are specific to sites of secondary neurodegeneration. Glia. 2017;65(12):1885–99.

    Article  PubMed  Google Scholar 

  105. Morrison HW, Filosa JA. A quantitative spatiotemporal analysis of microglia morphology during ischemic stroke and reperfusion. J Neuroinflammation. 2013;10:4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  106. Raj DD, Jaarsma D, Holtman IR, Olah M, Ferreira FM, Schaafsma W, et al. Priming of microglia in a DNA-repair deficient model of accelerated aging. Neurobiol Aging. 2014;35(9):2147–60.

    Article  CAS  PubMed  Google Scholar 

  107. Deczkowska A, Matcovitch-Natan O, Tsitsou-Kampeli A, Ben-Hamo S, Dvir-Szternfeld R, Spinrad A, et al. Mef2C restrains microglial inflammatory response and is lost in brain ageing in an IFN-I-dependent manner. Nat Commun. 2017;8(1):717.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  108. Hinwood M, Tynan RJ, Charnley JL, Beynon SB, Day TA, Walker FR. Chronic stress induced remodeling of the prefrontal cortex: structural re-organization of microglia and the inhibitory effect of minocycline. Cereb Cortex. 2013;23(8):1784–97.

    Article  PubMed  Google Scholar 

  109. Hellwig S, Brioschi S, Dieni S, Frings L, Masuch A, Blank T, et al. Altered microglia morphology and higher resilience to stress-induced depression-like behavior in CX3CR1-deficient mice. Brain Behav Immun. 2016;55:126–37.

    Article  PubMed  Google Scholar 

  110. Zheng LS, Kaneko N, Sawamoto K. Minocycline treatment ameliorates interferon-alpha- induced neurogenic defects and depression-like behaviors in mice. Front Cell Neurosci. 2015;9:5.

    PubMed  PubMed Central  Google Scholar 

  111. Akiyama H, Ikeda K, Katoh M, McGeer EG, McGeer PL. Expression of MRP14, 27E10, interferon-alpha and leukocyte common antigen by reactive microglia in postmortem human brain tissue. J Neuroimmunol. 1994;50(2):195–201.

    Article  CAS  PubMed  Google Scholar 

  112. Roy ER, Wang B, Wan YW, Chiu G, Cole A, Yin Z, et al. Type I interferon response drives neuroinflammation and synapse loss in Alzheimer disease. J Clin Invest. 2020;130(4):1912–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  113. Yamada T, Horisberger MA, Kawaguchi N, Moroo I, Toyoda T. Immunohistochemistry using antibodies to alpha-interferon and its induced protein, MxA, in Alzheimer’s and Parkinson’s disease brain tissues. Neurosci Lett. 1994;181(1–2):61–4.

    Article  CAS  PubMed  Google Scholar 

  114. Taylor JM, Minter MR, Newman AG, Zhang M, Adlard PA, Crack PJ. Type-1 interferon signaling mediates neuro-inflammatory events in models of Alzheimer’s disease. Neurobiol Aging. 2014;35(5):1012–23.

    Article  CAS  PubMed  Google Scholar 

  115. Walker DG, Lue LF. Immune phenotypes of microglia in human neurodegenerative disease: challenges to detecting microglial polarization in human brains. Alzheimers Res Ther. 2015;7(1):56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  116. Mattiace LA, Davies P, Dickson DW. Detection of HLA-DR on microglia in the human brain is a function of both clinical and technical factors. Am J Pathol. 1990;136(5):1101–14.

    CAS  PubMed  PubMed Central  Google Scholar 

  117. Castelnau PA, Garrett RS, Palinski W, Witztum JL, Campbell IL, Powell HC. Abnormal iron deposition associated with lipid peroxidation in transgenic mice expressing interleukin-6 in the brain. J Neuropathol Exp Neurol. 1998;57(3):268–82.

    Article  CAS  PubMed  Google Scholar 

  118. Keren-Shaul H, Spinrad A, Weiner A, Matcovitch-Natan O, Dvir-Szternfeld R, Ulland TK, et al. A unique microglia type associated with restricting development of Alzheimer’s disease. Cell. 2017;169(7):1276-90.e17.

    Article  CAS  PubMed  Google Scholar 

  119. Hansen DV, Hanson JE, Sheng M. Microglia in Alzheimer’s disease. J Cell Biol. 2018;217(2):459–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  120. Atagi Y, Liu CC, Painter MM, Chen XF, Verbeeck C, Zheng H, et al. Apolipoprotein E is a ligand for triggering receptor expressed on myeloid cells 2 (TREM2). J Biol Chem. 2015;290(43):26043–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  121. Shi Y, Yamada K, Liddelow SA, Smith ST, Zhao L, Luo W, et al. ApoE4 markedly exacerbates tau-mediated neurodegeneration in a mouse model of tauopathy. Nature. 2017;549(7673):523–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  122. Forstreuter F, Lucius R, Mentlein R. Vascular endothelial growth factor induces chemotaxis and proliferation of microglial cells. J Neuroimmunol. 2002;132(1–2):93–8.

    Article  CAS  PubMed  Google Scholar 

  123. Ding X, Gu R, Zhang M, Ren H, Shu Q, Xu G, et al. Microglia enhanced the angiogenesis, migration and proliferation of co-cultured RMECs. BMC Ophthalmol. 2018;18(1):249.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  124. Routhe LJ, Andersen IK, Hauerslev LV, Issa II, Moos T, Thomsen MS. Astrocytic expression of ZIP14 (SLC39A14) is part of the inflammatory reaction in chronic neurodegeneration with iron overload. Glia. 2020;68(9):1810–23.

    Article  PubMed  Google Scholar 

  125. Bishop GM, Dang TN, Dringen R, Robinson SR. Accumulation of non-transferrin-bound iron by neurons, astrocytes, and microglia. Neurotox Res. 2011;19(3):443–51.

    Article  CAS  PubMed  Google Scholar 

  126. Liuzzi JP, Aydemir F, Nam H, Knutson MD, Cousins RJ. Zip14 (Slc39a14) mediates non-transferrin-bound iron uptake into cells. Proc Natl Acad Sci U S A. 2006;103(37):13612–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  127. Li W, Viengkhou B, Denyer G, West PK, Campbell IL, Hofer MJ. Microglia have a more extensive and divergent response to interferon-alpha compared with astrocytes. Glia. 2018.

  128. Carr DJ, Veress LA, Noisakran S, Campbell IL. Astrocyte-targeted expression of IFN-alpha1 protects mice from acute ocular herpes simplex virus type 1 infection. J Immunol. 1998;161(9):4859–65.

    CAS  PubMed  Google Scholar 

  129. Hammond TR, Dufort C, Dissing-Olesen L, Giera S, Young A, Wysoker A, et al. Single-cell RNA sequencing of microglia throughout the mouse lifespan and in the injured brain reveals complex cell-state changes. Immunity. 2019;50(1):253-71.e6.

    Article  CAS  PubMed  Google Scholar 

  130. Zagórska A, Través PG, Lew ED, Dransfield I, Lemke G. Diversification of TAM receptor tyrosine kinase function. Nat Immunol. 2014;15(10):920–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  131. Fourgeaud L, Través PG, Tufail Y, Leal-Bailey H, Lew ED, Burrola PG, et al. TAM receptors regulate multiple features of microglial physiology. Nature. 2016;532(7598):240–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  132. Huang Y, Happonen KE, Burrola PG, O’Connor C, Hah N, Huang L, et al. Microglia use TAM receptors to detect and engulf amyloid β plaques. Nat Immunol. 2021;22(5):586–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  133. Fraser DA, Pisalyaput K, Tenner AJ. C1q enhances microglial clearance of apoptotic neurons and neuronal blebs, and modulates subsequent inflammatory cytokine production. J Neurochem. 2010;112(3):733–43.

    Article  CAS  PubMed  Google Scholar 

  134. Hong S, Beja-Glasser VF, Nfonoyim BM, Frouin A, Li S, Ramakrishnan S, et al. Complement and microglia mediate early synapse loss in Alzheimer mouse models. Science. 2016;352(6286):712–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  135. Chauhan P, Lokensgard JR. Glial cell expression of PD-L1. Int J Mol Sci. 2019;20:7.

    Google Scholar 

  136. Goddery EN, Fain CE, Lipovsky CG, Ayasoufi K, Yokanovich LT, Malo CS, et al. Microglia and perivascular macrophages act as antigen presenting cells to promote CD8 T cell infiltration of the brain. Front Immunol. 2021;12:726421.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  137. Nugent AA, Lin K, van Lengerich B, Lianoglou S, Przybyla L, Davis SS, et al. TREM2 regulates microglial cholesterol metabolism upon chronic phagocytic challenge. Neuron. 2020;105(5):837-54.e9.

    Article  CAS  PubMed  Google Scholar 

  138. Poliani PL, Wang Y, Fontana E, Robinette ML, Yamanishi Y, Gilfillan S, et al. TREM2 sustains microglial expansion during aging and response to demyelination. J Clin Invest. 2015;125(5):2161–70.

    Article  PubMed  PubMed Central  Google Scholar 

  139. Kang SS, Ebbert MTW, Baker KE, Cook C, Wang X, Sens JP, et al. Microglial translational profiling reveals a convergent APOE pathway from aging, amyloid, and tau. J Exp Med. 2018;215(9):2235–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  140. Pelisch N, Rosas Almanza J, Stehlik KE, Aperi BV, Kroner A. CCL3 contributes to secondary damage after spinal cord injury. J Neuroinflammation. 2020;17(1):362.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  141. Lowry JR, Klegeris A. Emerging roles of microglial cathepsins in neurodegenerative disease. Brain Res Bull. 2018;139:144–56.

    Article  CAS  PubMed  Google Scholar 

  142. Nakanishi H. Microglial cathepsin B as a key driver of inflammatory brain diseases and brain aging. Neural Regen Res. 2020;15(1):25–9.

    Article  PubMed  Google Scholar 

  143. Liu Y, Liu K, Huang Y, Sun M, Tian Q, Zhang S, et al. TRIM25 promotes TNF-α-induced NF-κB activation through potentiating the K63-linked ubiquitination of TRAF2. J Immunol. 2020;204(6):1499–507.

    Article  CAS  PubMed  Google Scholar 

  144. Torre S, Polyak MJ, Langlais D, Fodil N, Kennedy JM, Radovanovic I, et al. USP15 regulates type I interferon response and is required for pathogenesis of neuroinflammation. Nat Immunol. 2017;18(1):54–63.

    Article  CAS  PubMed  Google Scholar 

  145. Sala Frigerio C, Wolfs L, Fattorelli N, Thrupp N, Voytyuk I, Schmidt I, et al. The major risk factors for Alzheimer’s disease: age, sex, and genes modulate the microglia response to Aβ plaques. Cell Rep. 2019;27(4):1293-306.e6.

    Article  CAS  PubMed  Google Scholar 

  146. Oji S, Nicolussi EM, Kaufmann N, Zeka B, Schanda K, Fujihara K, et al. Experimental neuromyelitis optica induces a type I interferon signature in the spinal cord. PLoS ONE. 2016;11(3):e0151244.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  147. Wlodarczyk A, Khorooshi R, Marczynska J, Holtman IR, Burton M, Jensen KN, et al. Type I interferon-activated microglia are critical for neuromyelitis optica pathology. Glia. 2020;9:89.

    Google Scholar 

  148. Harmel J, Ringelstein M, Ingwersen J, Mathys C, Goebels N, Hartung HP, et al. Interferon-beta-related tumefactive brain lesion in a Caucasian patient with neuromyelitis optica and clinical stabilization with tocilizumab. BMC Neurol. 2014;14:247.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  149. Kim SH, Kim W, Li XF, Jung IJ, Kim HJ. Does interferon beta treatment exacerbate neuromyelitis optica spectrum disorder? Mult Scler. 2012;18(10):1480–3.

    Article  CAS  PubMed  Google Scholar 

  150. Palace J, Leite MI, Nairne A, Vincent A. Interferon beta treatment in neuromyelitis optica: increase in relapses and aquaporin 4 antibody titers. Arch Neurol. 2010;67(8):1016–7.

    Article  PubMed  Google Scholar 

  151. Shimizu Y, Yokoyama K, Misu T, Takahashi T, Fujihara K, Kikuchi S, et al. Development of extensive brain lesions following interferon beta therapy in relapsing neuromyelitis optica and longitudinally extensive myelitis. J Neurol. 2008;255(2):305–7.

    Article  PubMed  Google Scholar 

  152. Uzawa A, Mori M, Hayakawa S, Masuda S, Kuwabara S. Different responses to interferon beta-1b treatment in patients with neuromyelitis optica and multiple sclerosis. Eur J Neurol. 2010;17(5):672–6.

    Article  CAS  PubMed  Google Scholar 

  153. Jakimovski D, Kolb C, Ramanathan M, Zivadinov R, Weinstock-Guttman B. Interferon β for multiple sclerosis. Cold Spring Harb Perspect Med. 2018;8:11.

    Article  CAS  Google Scholar 

  154. Araki M, Matsuoka T, Miyamoto K, Kusunoki S, Okamoto T, Murata M, et al. Efficacy of the anti-IL-6 receptor antibody tocilizumab in neuromyelitis optica: a pilot study. Neurology. 2014;82(15):1302–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  155. Ringelstein M, Ayzenberg I, Harmel J, Lauenstein AS, Lensch E, Stogbauer F, et al. Long-term therapy with interleukin 6 receptor blockade in highly active neuromyelitis optica spectrum disorder. JAMA Neurol. 2015;72(7):756–63.

    Article  PubMed  Google Scholar 

  156. Araki M, Aranami T, Matsuoka T, Nakamura M, Miyake S, Yamamura T. Clinical improvement in a patient with neuromyelitis optica following therapy with the anti-IL-6 receptor monoclonal antibody tocilizumab. Mod Rheumatol. 2013;23(4):827–31.

    Article  CAS  PubMed  Google Scholar 

  157. Ayzenberg I, Kleiter I, Schroder A, Hellwig K, Chan A, Yamamura T, et al. Interleukin 6 receptor blockade in patients with neuromyelitis optica nonresponsive to anti-CD20 therapy. JAMA Neurol. 2013;70(3):394–7.

    Article  PubMed  Google Scholar 

  158. Vezzani A, Balosso S, Ravizza T. The role of cytokines in the pathophysiology of epilepsy. Brain Behav Immun. 2008;22(6):797–803.

    Article  CAS  PubMed  Google Scholar 

  159. Kramer K, Schaudien D, Eisel UL, Herzog S, Richt JA, Baumgärtner W, et al. TNF-overexpression in Borna disease virus-infected mouse brains triggers inflammatory reaction and epileptic seizures. PLoS ONE. 2012;7(7):e41476.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  160. Campbell IL, Stalder AK, Chiang CS, Bellinger R, Heyser CJ, Steffensen S, et al. Transgenic models to assess the pathogenic actions of cytokines in the central nervous system. Mol Psychiatry. 1997;2(2):125–9.

    Article  CAS  PubMed  Google Scholar 

  161. Campbell IL. Transgenic mice and cytokine actions in the brain: bridging the gap between structural and functional neuropathology. Brain Res Brain Res Rev. 1998;26(2–3):327–36.

    Article  CAS  PubMed  Google Scholar 

  162. Sominsky L, Dangel T, Malik S, De Luca SN, Singewald N, Spencer SJ. Microglial ablation in rats disrupts the circadian system. Faseb j. 2021;35(2):e21195.

    Article  CAS  PubMed  Google Scholar 

  163. De Luca SN, Sominsky L, Soch A, Wang H, Ziko I, Rank MM, et al. Conditional microglial depletion in rats leads to reversible anorexia and weight loss by disrupting gustatory circuitry. Brain Behav Immun. 2019;77:77–91.

    Article  PubMed  Google Scholar 

  164. Badimon A, Strasburger HJ, Ayata P, Chen X, Nair A, Ikegami A, et al. Negative feedback control of neuronal activity by microglia. Nature. 2020;9:78.

    Google Scholar 

  165. Han J, Zhu K, Zhang XM, Harris RA. Enforced microglial depletion and repopulation as a promising strategy for the treatment of neurological disorders. Glia. 2019;67(2):217–31.

    Article  PubMed  Google Scholar 

  166. Green KN, Crapser JD, Hohsfield LA. To kill a microglia: a case for CSF1R inhibitors. Trends Immunol. 2020;41(9):771–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  167. Liddelow SA, Barres BA. Reactive astrocytes: production, function, and therapeutic potential. Immunity. 2017;46(6):957–67.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Sue Ling Lim for technical assistance.

Funding

P.K.W. was the recipient of an Australian Government Research Training Program scholarship. O.B. is supported by BrightFocus Foundation 2020A016806; NIH/NINDS R01 NS088137, NIH/NIA R01 AG051812 and R01 AG054672; NIH/NEI R01 EY027921; National Health and Medical Research Council Australia (RG180378).

Author information

Authors and Affiliations

Authors

Contributions

PKW, ILC and MJH conceived the study. PKW, ILC and MJH designed the experiments. PKW, TMA, BV, EH and SRJ performed the experiments. BV developed an optimized passive tissue clearing technique. PKW analyzed experimental results. PKW, ANM and BG processed and analyzed the sequencing data. PKW and TMA designed, performed and analyzed flow cytometry. OB provided key resources. PKW, ILC and MJH wrote the manuscript. All authors contributed to the review and editing of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Markus J. Hofer.

Ethics declarations

Ethics approval and consent to participate

All animal experiments were performed in compliance with the NSW Animal Research Act and its associated regulations and the 2013 NHMRC Australian code for the care and use of animals for scientific purposes. Ethical approval for the use of all mice was obtained from the animal ethics committee of the University of Sydney (2014/699, 2015/758 and 2018/1428).

Consent for publication

Not applicable.

Competing interests

B.G. is a director of Pacific Analytics PTY LTD & SMRTR PTY LTD, Australia, which had no role in the design, execution, analysis or preparation of the manuscript. O.B. has collaborations with Sanofi, GSK, Regulus Therapeutics; research funding from Sanofi, GSK, miRagen Therapeutics, honoraria for lectures, consultancy: Camp4. The other authors declare no other competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Supplementary figures 1–11.

Additional file 2.

Table S1. Differential gene expression analysis and top 100 enriched biological processes (BPs) of microglia from WT, GFAP-IL6 and GFAP-IFN mice from Fig. 3.

Additional file 3. Table S2

Reads (FPKM) for common, IL-6- and IFN-α-skewed genes and top 100 enriched BPs by upregulated genes from Fig. 4.

Additional file 4. Table S3.

Reads (FPKM) and log2 fold-changes from meta-analysis and the top 20 enriched BPs by gene clusters of interest in Fig. 7 and Fig. S8.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

West, P.K., McCorkindale, A.N., Guennewig, B. et al. The cytokines interleukin-6 and interferon-α induce distinct microglia phenotypes. J Neuroinflammation 19, 96 (2022). https://doi.org/10.1186/s12974-022-02441-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12974-022-02441-x

Keywords