- Open Access
Microglia RAGE exacerbates the progression of neurodegeneration within the SOD1G93A murine model of amyotrophic lateral sclerosis in a sex-dependent manner
Journal of Neuroinflammation volume 18, Article number: 139 (2021)
Burgeoning evidence highlights seminal roles for microglia in the pathogenesis of neurodegenerative diseases including amyotrophic lateral sclerosis (ALS). The receptor for advanced glycation end products (RAGE) binds ligands relevant to ALS that accumulate in the diseased spinal cord and RAGE has been previously implicated in the progression of ALS pathology.
We generated a novel mouse model to temporally delete Ager from microglia in the murine SOD1G93A model of ALS. Microglia Ager deficient SOD1G93A mice and controls were examined for changes in survival, motor function, gliosis, motor neuron numbers, and transcriptomic analyses of lumbar spinal cord. Furthermore, we examined bulk-RNA-sequencing transcriptomic analyses of human ALS cervical spinal cord.
Transcriptomic analysis of human cervical spinal cord reveals a range of AGER expression in ALS patients, which was negatively correlated with age at disease onset and death or tracheostomy. The degree of AGER expression related to differential expression of pathways involved in extracellular matrix, lipid metabolism, and intercellular communication. Microglia display increased RAGE immunoreactivity in the spinal cords of high AGER expressing patients and in the SOD1G93A murine model of ALS vs. respective controls. We demonstrate that microglia Ager deletion at the age of symptomatic onset, day 90, in SOD1G93A mice extends survival in male but not female mice. Critically, many of the pathways identified in human ALS patients that accompanied increased AGER expression were significantly ameliorated by microglia Ager deletion in male SOD1G93A mice.
Our results indicate that microglia RAGE disrupts communications with cell types including astrocytes and neurons, intercellular communication pathways that divert microglia from a homeostatic to an inflammatory and tissue-injurious program. In totality, microglia RAGE contributes to the progression of SOD1G93A murine pathology in male mice and may be relevant in human disease.
Amyotrophic lateral sclerosis (ALS), a progressive, fatal, neurodegenerative disease, is characterized by the inexorable death of motor neurons. Affected patients experience progressive muscle atrophy, motor function decline, and eventual paralysis. ALS rapidly progresses in most patients, resulting in an estimated survival of 2–4 years after onset of symptoms . Greater than 90% of ALS patients have no known family history of the disease and hence the non-familial form of the disease is known as sporadic ALS (sALS). A smaller subset of patients (5–10%) is affected by familial ALS (fALS) . Mutations in superoxide dismutase 1 (SOD1) account for approximately 20% of fALS patients . Regardless of the underlying disease-related mutation, ALS is characterized by the accumulation of proteinaceous aggregates, which is associated with a pro-damage inflammatory and oxidative stress state involving multiple cell types. Ultimately, these processes collectively augur motor neuron death [2, 3]. Mice carrying multiple copies of a human SOD1G93A transgene display hallmarks of ALS, including accumulation of proteinaceous aggregates, neuroinflammation, and progressive motor function decline with eventual paralysis and death in all SOD1G93A transgenic mice [4, 5]. While ALS results in motor neuron death, several landmark studies have elucidated critical roles for non-neuronal cells during disease progression [6,7,8,9,10,11]. Altogether, these considerations illustrate that ALS pathology is driven by dysfunction across a myriad of cell types whose intercommunications drive processes that irreparably damage neurons.
Spatial transcriptomic analysis of SOD1G93A mouse spinal cords suggests that microglial dysfunction is evident by post-natal day 30 (P30), with astrocyte dysfunction becoming apparent later, by P70. These findings suggest that microglial dysfunction may precede and contribute to astrocyte dysfunction, which precedes neuronal death . Critically, these key findings implicating glial cells in ALS pathobiology were also observed in spatial transcriptomics analyses of human ALS spinal cord . Studies in SOD1G93A mice affirmed mediating roles for glial cells in impacting survival in SOD1G93A mice; when the myeloid and lymphoid compartment of these mice was replaced with a wild-type (WT) bone marrow transplant within 24 h of birth, the recipient SOD1G93A mice demonstrated improved survival compared to those SOD1G93A mice receiving SOD1G93A bone marrow . Moreover, an expanding number of genes with known ALS-linked mutations mediate myeloid function, including C9ORF72, TARDP, and OPTN [13, 14]. Specifically, C9ORF72 regulates phagosome and lysosome pathways , while TARDP regulates phagocytosis . OPTN deficiency induces an inflammatory profile in microglia . It has previously been shown that reducing myeloid inflammation increases lifespan of SOD1G93A mice [7, 18]. In a distinct inducible murine model of ALS, microglia were required for the clearance of protein aggregates . Microglia from SOD1G93A mice display unique transcriptomic signatures reminiscent of both protective and inflammatory macrophages . Altogether, these findings suggest that microglia may exert opposing roles in ALS, mediating stage-dependent alterations in protective vs. damage-provoking functions. In fact, accruing evidence suggests that microglia may undergo step-wise activation toward a dysfunctional inflammatory disease-associated microglia (DAM) phenotype, at least in the context of Alzheimer’s disease (AD) [21, 22]. Triggering receptor expressed on myeloid cells 2 (TREM2) has been implicated in the step-wise activation of microglia toward this DAM state; however, it is established that additional factors synergize or act independently of TREM2 in this “activation pathway” and that TREM2 plays complex pro- and anti-inflammatory effects in immune cells . The identity of such additional factors is thus essential to discover, as they may unveil new therapeutic targets for neurodegenerative disorders.
The receptor for advanced glycation end products (RAGE) is an immunoglobulin (Ig)-type transmembrane receptor expressed by many cell types and is upregulated in patient and murine ALS tissues [24,25,26]. RAGE binds a diverse set of ligands, particularly those noted as damage associated molecular patterns (DAMPs), including several upregulated in ALS patient and SOD1G93A mouse spinal cord, such as S100 calcium binding protein B, carboxymethyllysine (CML)-advanced glycation end product (AGE), and high-mobility group box 1 (HMGB1) [25, 27,28,29]. Spatial transcriptomic analyses recently identified an enrichment of the “AGE-RAGE pathway in diabetic complications” Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway in both murine and human ALS spinal cord within glia-related modules of gene expression .
While the discrete role of microglia RAGE in the context of ALS has never been studied, RAGE activation on myeloid cells, including microglia, stimulates NF-κB activity, thereby promoting pro-inflammatory cytokine expression and generation of reactive oxygen species [30,31,32,33,34]. Ligand engagement with RAGE induces a conformational change in its cytoplasmic tail that signals, at least in part, through diaphanous related formin 1 (DIAPH1) [35, 36]. Microglia have been shown to induce motor neuron death and drive pathology in SOD1G93A mice via activation of NF-κB . We hypothesize that microglia RAGE contributes to oxidative stress and exacerbation of pro-inflammatory gene expression, processes which divert microglia from their homeostatic and protective functions, thereby redirecting them to a phenotype that contributes to neurodegeneration.
In the current study, analysis of human ALS patient RNA-sequencing data uncovered for the first time that a spectrum of expression of AGER, the gene encoding RAGE, was evident in ALS cervical spinal cord tissue, and that this was related to alterations in pathways involving lipid metabolism, cellular microenvironment, and intercellular communication. Furthermore, patients with high AGER expression displayed increased levels of RAGE in microglia in the cervical spinal cord. In a murine model, we demonstrate a significant survival benefit and improvement in motor function performance in male but not female microglia Ager-deficient SOD1G93A Agerfl/fl Cx3cr1Cre/+ mice (microglia Ager deletion) relative to SOD1G93A Ager+/+ Cx3cr1Cre/+ mice (Cre-controls; microglia express Ager) subjected to deletion of microglia Ager at 3 months of age. Remarkably, deletion of microglia Ager from SOD1G93A mice modulated the lumbar spinal cord transcriptome in a manner such that a number of pathways related to increased AGER expression in human ALS patients were significantly ameliorated. Altogether, these data suggest that RAGE may contribute to the pathological re-programming of innate microglia functions in the ALS spinal cord.
Materials and methods
All experiments were performed on C57BL/6J mice intercrossed with the following strains: B6.129P2(Cg)-CX3CR1tm2.1(cre/ERT2)Litt/WganJ (JAX Stock No: 021160), AgerFlox/Flox mice (C56BL6J background), and B6.Cg-Tg(SOD1G93A)1Gur/J (JAX Stock No:004435) [4, 37,38,39]. Note that all mice employed in this study were extensively backcrossed into or generated directly in C57BL/6J background. In all mice, copy number was analyzed and any mice with alterations were not enrolled in the study . Three mice over the course of the study were censored due to non-ALS deaths. All mice were maintained under pathogen-free conditions. Unless otherwise noted, all mice were housed in a temperature (19–23 °C) and humidity (30–70% relative humidity)-regulated environment with 12-h light/dark cycle, lights on/off at 6:30, and received standard chow food pellets (Lab Diets, Cat: 5053) and water ad libitum. Tamoxifen (TAM) was used to delete Ager from Cx3cr1-expressing cells at 90 days of age. All animals, irrespective of genotype, were administered 0.2 μg of TAM (Sigma, Cat: T5648) dissolved in corn oil (Sigma-Aldrich, Cat: C8267) once daily every day for a total of five intraperitoneal injections except for the female control mice used in validation of Ager deletion, which received five intraperitoneal injections administered every other day. Onset of disease was identified by the age at which the animal reached maximum weight. The humane endpoint was determined by 20% weight loss from maximum weight, or the inability of the animal to right itself within 15 s of being placed on its side. When a mouse reached the humane endpoint, by either criterion, it was immediately sacrificed and date of death was recorded for survival analyses. Mice were weighed 2–3 times weekly, same time of day, with calibrated balances.
Measurement of motor function was performed using the Hanging Wire test as detailed previously . Briefly, mice were placed on a cage wire lid, inverted 50 cm above padding and observed for up to 1 min. The amount of time until the mouse fell to the padding was recorded. All mice were trained to the procedure for 1 week prior to data acquisition (8 weeks of age). All tests were performed twice-weekly until the mice were no longer able to perform the test. In all cases, the operator was naïve to experimental code. Triplicate measurements were taken and mean value was recorded for each session and weekly means were calculated. Animals were given 30-s breaks between test replicates.
Murine tissue staining
Mice were anesthetized with Ketamine and Xylazine and then underwent rapid cardiac perfusion with PBS, followed by 4% paraformaldehyde (PFA). The spinal column was rapidly removed, and the lumbar spinal cord was retrieved. Lumbar spinal cords were washed in PBS and then drop-fixed in 4% PFA for an additional 1 h at 4 °C. Gastrocnemius tissue was collected and drop-fixed 16 h in 1% PFA at 4 °C. After fixation, tissue was incubated with 15% sucrose in PBS for 24 h followed by incubation with 30% sucrose in PBS for an additional 24 h. Cyroprotected tissue was frozen in optimal cutting temperature (OCT) compound (FisherScientific, Cat: 23-730-571) and kept at − 80 °C until sectioning. Furthermore, 8-μm-thick serial sections were collected on Superfrost PLUS slides (FisherScientific, Cat: 22-037-246) using a Microm cryostat (ThermoFisher, Model: HM550). Sections were washed 3× with PBS for 5 min and then permeabilized with 0.2% Triton-X 100 in PBS for 10 min and washed 3× with PBS. Blocking was conducted with Serum-Free blocking buffer (Dako, Cat: X090930-2) for 1 h at RT. All primary antibodies were diluted in Antibody Diluent (Dako, Cat: S3022) and applied overnight at 4 °C. Subsequently, slides were washed 3× with PBS, then incubated with secondary antibodies diluted in Antibody Diluent for 1 h at RT. For slides containing DAPI, slides were washed and placed in a 1 μg/mL DAPI (Invitrogen, Cat: D1306) solution for 5 min at RT and washed 3× with PBS before mounting with fluorescent mounting media (Dako, Cat: S302380-2). Primary antibodies used: 1 μg/mL Mouse anti-RAGE (Millipore, Cat: MAB5328), 0.5 μg/mL Rat anti-CD11B (M1/70) (Invitrogen, Cat:14-0112-82), 0.25 μg/mL Rat anti-GFAP (2.2B10) (Invitrogen, Cat: 13-0300), 5 μg/mL Rat anti-CLEC7A (Invivogen, Cat: mabg-mdect), 1 μg/mL Mouse anti-AMIGO2 (G-7) (SantaCruz, Cat: sc-373699), 5 μg/mL Mouse anti-NeuN (Millipore, Cat: MAB377), 5 μg/mL Rat anti-CD68 (Abcam, Cat:ab53444), 5 μg/mL Rat anti-F4/80 (Abcam, Cat: ab6640), and 3 μg/mL Chicken anti-MAP2 (Abcam, Cat:ab5392). Secondary antibodies utilized: Donkey anti-Rat Alexa Fluor 488 (Invitrogen, Cat: A-21208), Donkey anti-Mouse Alexa Fluor 546 (Invitrogen, Cat: A10036), Donkey anti-Rat Alexa Fluor 594 (Invitrogen, Cat: A-21209), Donkey anti-Mouse Alexa Fluor 488 (Invitrogen, Cat: A21203), and Donkey anti-Chicken Alexa Fluor 647 (Jackson ImmunoResearch, Cat: 703-605-155). All secondary antibodies were used at 1 μg/mL. All experiments included negative controls by omission of primary antibody.
Human tissue staining
De-identified paraffin-embedded tissue sections from sporadic ALS patient and control cervical spinal cord tissue were provided by the Target ALS Multicenter Postmortem Tissue Core (www.targetals.org). Then, 3 × 5 min washes with Clear-Rite 3 (ThermoFisher, Cat: 6901TS) were used to de-paraffinize the sections. Slides were rehydrated in a series of 5 min EtOH washes (2 × 100%, 1 × 90%, 1 × 70%) and then washed 3× with ddH20. Antigen retrieval was performed by steaming in Epitope Retrieval Solution (IHC World, Cat: IW-1100) for 1 h then washed 3× with PBS. Sections were permeabilized for 10 min with 0.4% Triton X100 in PBS then washed 3× with PBS. Autofluorescence was reduced by incubating with 1× TrueBlack in 70% ethanol (Biotium, Cat: 23007) for 30 s, then washed 3× with large volumes of PBS. Blocking was conducted with 5% Donkey Serum (SigmaAldrich, Cat: D9663) in PBS overnight at 4 °C. All primary antibodies were diluted in 2.5% Donkey Serum in PBS and applied for 24 h at 4 °C. Subsequently, slides were washed 3× with PBS, then incubated in secondary antibodies diluted in 2.5% Donkey Serum in PBS for 1 h at RT. Slides were then washed 3× with PBS before mounting with fluorescent mounting media. Primary antibodies utilized: 5 μg/mL Rabbit anti-IBA1 (Wako, Cat: 013-27691), 4 μg/mL Goat anti-RAGE (R&D, Cat: AF1179), 5 μg/mL Mouse anti-GFAP (BD Biosciences, Cat: 556330), and 5 μg/mL Mouse anti-NeuN (Millipore, Cat: MAB377). Secondary antibodies utilized: Donkey anti-rabbit Alexa Fluor 488 (Invitrogen, Cat: A21206), Donkey anti-mouse Alexa Fluor 488 (Invitrogen, Cat: A21203), Donkey anti-Goat Alexa Fluor 594 (Invitrogen, Cat: A32758), and Donkey anti-Goat Alexa Fluor 647 (Invitrogen, Cat: A32849). All secondary antibodies were used at 1 μg/mL. As above, all experiments included negative controls with omission of primary antibodies.
Imaging and quantification
Multicolor wide-field images were taken on a Leica 5500B microscope at × 20 or × 40 magnification, as indicated. All microscope settings were kept identical for each experiment. For mouse tissue imaging: 2–4 images/tissue slice and 3–4 tissue slices/sample were collected for analysis. For human tissue imaging: 2–4 images/tissue slice and 3 tissues slices/patient were collected for analysis, four regions of interest per image,150 × 150 μm in size, were selected as to exclude auto-fluorescent blood vessels for each image before further analysis. All analyses were performed with the Fiji distribution of ImageJ (NIH) .
Cell number analyses
Cell numbers that displayed overlap of the designated antibodies with DAPI were manually counted per image and averaged per mouse. All analysis was completed by a naïve experimenter blinded to all experimental codes.
Cell area analysis
To quantify the positive area of each stain in end-stage murine and human tissues, the images underwent background removal with the rolling ball radius set to 50 within ImageJ. Images were then subjected to automated thresholding with the optimal thresholding algorithm of each signal being selected by a naïve experimenter. The following thresholding algorithms were utilized: Li for: CD11B in Wild-Type vs SOD1G93A mice; Moments used for: RAGE, c-Type lectin domain containing 7a protein (or dectin-1) (CLEC7A), glial fibrillary acidic protein (GFAP), microtubule-associated protein 2 (MAP2), and adhesion molecule with Ig like domain 2 (AMIGO2); Triangle used for: CD11B, CLEC7A, on 120 day old tissue, and F4/80, CD68, and neuronal nuclei antigen (NeuN); and Otsu for: human ionized calcium-binding adapter molecule 1 (IBA1), and human GFAP [43,44,45]. Calculation of positive area was calculated per μm2 and the average of 4–8 images per sample was used for statistical analysis.
BV2 mouse microglia-like cells were obtained as a generous gift from Dr. Colin K. Combs (University of Nebraska). BV2 cells were grown in growth medium (Dulbecco’s modified Eagle’s medium (DMEM) (Thermofisher, Cat:11885-084) containing 10% (v/v) heat-inactivated fetal bovine serum (FBS) (Corning, Cat:35-010-CV)) with 100 units/mL penicillin and 100 μg/mL streptomycin (Thermofisher, Cat:10378016). The cells were incubated at 37 °C under 5% CO2 and 95% relative humidity.
BV2 cells were grown to 50% confluency in normal growth medium. Growth medium was replaced with fresh growth medium supplemented with 5 μg/mL polybrene (Sigma, Cat:TR-1003) and 1.5 × 106 lentiviral particles of either shAger (Sigma, Cat:TRCN0000071745) or non-targeting shRNA control (Sigma, Cat: SHC002H) lentivirus for 8 h. When cells reached 90% confluency, shRNA expressing cells were selected using growth medium with 4 μg/mL puromycin (Thermofisher, Cat: A1113803) and were maintained under puromycin selection for at least three passages before use.
Treatment of BV2 cells
BV2 cells were plated in a 6-well plate at 3 × 105 cells per well and cultured in growth medium containing 1% FBS overnight. Cells were treated with 300 μg/mL RAGE ligand CML-AGE or human serum albumin control for 24 h before lysing in Qiazol reagent (QIAGEN Inc, Cat: 79306) and frozen until RNA was isolated . For RAGE inhibitor studies, cells were first pre-treated with either 10 μM RAGE inhibitor (a C11 analog) or 0.1% v/v dimethyl sulfoxide (DMSO) in 1% FBS containing growth medium for 2 h before treatment .
CD11B+ cell isolation via autoMACS
Mice were deeply anesthetized with Ketamine and Xylazine and then cardiac-perfused with 50 mL of ice-cold PBS. Subsequently, mice were decapitated and the brain was removed and collected into ice-cold tubes containing PBS. Brain tissue was dissociated using the Adult Brain Tissue Dissociation Kit (Miltenyi, Cat: 130-107-677) as per manufacturer’s guidelines with the gentleMACS dissociator (Miltenyi, Cat:130-093-235) and strained with 70 μm filters (Miltenyi, Cat: 130-098-462). Homogenates underwent debris removal as per the manufacturer’s guidelines (Miltenyi, Cat:130-109-398). Samples were then incubated with CD11B-microbeads (Miltenyi, Cat: 130-097-142) as per manufacturer’s guidelines, then washed with Running Buffer (Miltenyi, Cat: 130-091-221) and centrifuged at 500×g for 5 min. Pellets were resuspended with 2 mL Running Buffer and subjected to the autoMACS Pro Separator (Miltenyi, Cat:130-092-545) “Possel” selection protocol. Positive fractions were collected, centrifuged at 500×g for 5 min and then re-suspended in 200 μL Qiazol lysis reagent (Qiagen, Cat: 79306) and flash frozen until RNA was isolated.
Tissue harvest and RNA isolation
Mice were anesthetized with Ketamine and Xylazine. The spinal column was rapidly removed, and the lumbar spinal cord was retrieved. Lumbar spinal cords were flash frozen and stored at − 80 °C until processed. Frozen tissue was kept on dry ice until homogenized in Qiazol reagent (Qiagen Inc, Cat: 79306). Crude homogenate was mixed with chloroform, at a ratio of 70 μL chloroform per 350 μL homogenate, to facilitate removal of lipids, and allowed to equilibrate at room temperature for 2 min. Samples were centrifuged at 13,000×g for 15 min at 4 °C. The aqueous phase containing the RNA was transferred to a new tube. Then the following steps were completed as detailed by the manufacturer’s instructions using an RNeasy Mini Kit (Qiagen, Cat: 74104) with on-column DNAse digestion (Qiagen, Cat:79254). RNA concentration was determined via a NanoDrop spectrophotometer (ThermoFisher, Model: ND-1000). RNA integrity (RIN) for RNA-sequencing samples was measured using RNA 6000 Pico Kit in a 2100 Bioanalyzer (Agilent).
Real-time quantitative polymerase chain reaction
One microgram of RNA was used with the iScript cDNA synthesis kit (BioRad, Cat:1708890) as per the manufacturer’s instructions to generate cDNA. Taqman gene expression assays were used to evaluate Ager (Life Technologies, Cat: Mm01134790_g1), 18s rRNA (Applied Biosystems, Cat:4310893E), Hprt (Life Technologies, Cat:Mm03024075_m1), Il1a (Life Technologies, Cat:Mm00439620_m1), and Malat1 (Life Technologies, Cat: Mm01227912_s1) levels using the Taqman Fast Universal PCR Master Mix (Life Technologies, Cat: 4367846) with the 7500 Fast Real-Time PCR System (Applied Biosystems, Cat: LS4351106). Expression values were calculated using the ΔΔCt method relative to 18s rRNA or Hprt.
Five hundred nanograms of RNA per biological replicate (1.5 μg total) for each condition was pooled before proceeding with the mouse cytokine cDNA plate array (Signosis, Cat: AP-1141) as per the manufacturer’s instructions. Resulting luminescence was quantified using a SpectraMax M5 microplate reader (Molecular Devices, San Jose, CA) and normalized to 18s rRNA.
Murine RNA sequencing
High-quality lumbar spinal cord RNA samples, as confirmed by Bioanalyzer with RIN values ranging from 7.9 to 8.7, free of DNase and RNase were prepared for sequencing using the TruSeqStranded mRNA library prep kit (Illumina, Cat: 20020594) and the NovaSeq 6000 SP Reagent Kit v1.5 (Illumina, Cat: 20028401). High-throughput RNA sequencing (RNA-seq) was completed using an Illumina NovaSeq 6000 sequencer performed by the NYU Grossman School of Medicine Genome Technology Core.
Murine RNA-seq data analysis
The resulting 49–61 M read pairs per sample were processed following standard quality control practices to remove low-quality reads and adapter sequence contamination (~ 15% of reads) [47,48,49]. Remaining high-quality read pairs were aligned to the mouse genome (mm10) using STAR 2.7.3a with a mean of 88% uniquely mapping . Read pair counts per gene were summed with the featureCounts function in subread 1.6.3 using the GENCODE M25 annotation release [50,51,52,53]. Read counts were normalized using the trimmed mean of M values (TMM) method in edgeR 3.24.3 within the R environment 3.5.1 [54,55,56]. Differential expression was analyzed using edgeR via the exactTest function [55, 57]. Significantly differentially expressed genes were determined by a cut-off of adjusted p value < 0.05 and subjected to over representation analysis of KEGG pathways using the R package clusterProfiler v3.16.1 [58,59,60,61]. Significantly differentially expressed genes with false-discovery rate (FDR) < 0.05 and corresponding log2 fold changes were used as input for all ingenuity pathway analysis (IPA) including graphical summary, canonical pathways, upstream regulators, and causal network analyses (QIAGEN Inc.) . Adjustments were made for multiple testing to control the FDR at 0.05 .
Human RNA-seq data analysis
Raw RNA-seq data of cervical spinal cord and de-identified metadata were obtained from:The Target ALS Multicentered Postmortem Tissue Core, the New York Genome Center for Genomics of Neurodegenerative Disease, Amyotrophic Lateral Sclerosis Association, and TOW Foundation (www.targetals.org). Illumina paired-end 100-bp read data were downloaded from TargetALS (mean of 43 M reads per sample) and processed following standard quality control practices to remove low-quality reads and adapter sequence contamination (~ 4% of reads) [47,48,49]. Remaining high-quality read pairs were aligned to the human genome (hg19) using STAR 2.6.1d with a mean of 95% uniquely mapping. Read pair counts per gene were summed with the featureCounts function in subread 1.6.3 using the GENCODE 30 annotation release [50,51,52,53]. Only samples with RIN values of ≥ 6 were used for all analyses. Read counts were normalized using the TMM method on log2 counts per million in edgeR 3.24.3 within the R environment 3.5.1 [54,55,56]. Differential expression was analyzed using generalized linear models via edgeR, with either (1) diagnosis of ALS spectrum motor neuron disease (n = 76) or non-neurological control (n = 11) as categorical variables; or (2) AGER as a continuous variable in ALS patients (n = 76) to control and test for differential expression related to diagnosis or increasing AGER expression while controlling for sex and testing for any potential interactions within the model . Significant differentially expressed genes with log2 fold changes with absolute values of ≥ 0.5 used as the input for evaluation of overrepresentation of KEGG gene sets and for all IPA analyses as per the murine data analysis described above. Competitive gene set testing via CAMERA, and rotational gene set testing via ROAST were both completed using edgeR 3.24.3 to evaluate KEGG and Gene Ontology (GO) gene sets [55, 57,58,59,60, 65,66,67,68]. Testing for linear correlation was completed by first normalizing AGER expression using the TMM method and generating counts per million mapped reads (CPM) with edgeR 3.24.3 as above. A linear model was generated with normalized AGER CPM as the predictor of either age at death/tracheostomy or age at onset within the base R environment 3.5.1.
Data are shown as mean ±SEM or as indicated. Normality of the data was assessed using the Shapiro-Wilk’s normality test. If normality assumption was met, data were subsequently evaluated by independent two-sample t tests, two-way ANOVA with post-hoc Tukey’s test, one way ANOVA with post-hoc Holm-Šídák multiple comparisons test, or mixed effects analysis with Geisser-Greenhouse correction with post-hoc Holm-Sidak’s multiple comparisons test, as indicated. Non-parametric Mann–Whitney tests were implemented instead to assess differences if normality assumption was violated. Survival data were visualized by Kaplan–Meier curve and the Logrank test (Mantel–Cox) was used to evaluate differences in the survival distributions between groups. All analyses were performed with GraphPad Prism 9 (GraphPad Software, San Diego, CA) and R 3.6.1 using R package “nlme” [56, 69]. p values < 0.05 were used to denote statistical significance.
Evaluation of AGER expression in human amyotrophic lateral sclerosis spinal cord
Previous studies have documented increased RAGE expression in human ALS spectrum motor neuron disease (ALS) spinal cord [25, 26]; however, detailed assessments of roles for RAGE in ALS are lacking . To address this key point, we obtained bulk RNA-seq data from ALS and control patient cervical spinal cord tissue from the Target ALS Multicentered Postmortem Tissue Core, the New York Genome Center for Genomics of Neurodegenerative Disease, Amyotrophic Lateral Sclerosis Association and TOW Foundation. We limited our analysis to RNA-seq samples with RIN ≥ 6. In total, 76 ALS and 11 non-neurological control cervical spinal cord data met these criteria and were eligible for analysis.
Expression of the gene encoding RAGE, AGER, was not significantly different (FDR = 0.78) between ALS patients and non-neurological controls (Fig. 1A). However, when we examined gene set differences in expression between the ALS patients and non-neurological controls with rotational gene set testing (ROAST) and competitive gene set testing (CAMERA) analyses, we observed significant enrichment of the “AGE-RAGE signaling pathway in diabetic complications” which was visualized by a barcode enrichment plot for the pathway (Fig. 1B, see Additional Files 1 and 2, Supplemental Table 1.1-1.3). Altogether, these analyses suggest upregulation of the AGE-RAGE signaling pathway gene set in ALS patients relative to control patients.
Further analysis revealed that there was a spectrum of AGER expression (normalized counts) in the cervical spinal cord across ALS patients (Fig. 1A). Thus, it was logical to determine if the level of AGER correlated with available ALS patient phenotypic metadata. The ALS patients had a median onset of 60.5 years of age (range 32–80), median age at death of 65 years (range 32–80), and a median disease duration (onset to death) of 36mo (range 6–156 months), and the majority had no identifiable family history of ALS (see Additional Files 1 and 2, Supplemental Table1.4). We found normalized AGER counts per million (CPM) was negatively correlated with the age at disease onset (coefficient = − 2.449, p = 0.07, adjusted r2 = 0.032) and the age at death or tracheostomy (coefficient = − 2.609, p = 0.037, adjusted r2 = 0.045) in the ALS patients (Fig. 1C, D). We obtained a randomly-selected and available subset of tissue sections from the patients displaying highest vs. lowest 10% of AGER RNA values, that is, the extremes (Fig. 1E). We demonstrated that the difference in AGER RNA was also present at the protein level (Fig. 1F, G).
We next sought to examine genes and pathways that may be modulated with increasing AGER expression in ALS patients. We tested for differential expression using AGER as the predictive variable and found that there were many significantly differentially expressed genes associated with AGER expression (see Additional Files 1 and 2, Supplemental Table 1.5). Importantly, there were no differential genes related to AGER expression that differed significantly by sex. We narrowed subsequent analyses to differential genes with ≥ 0.5 |Log2| fold change; it is important to note that Log2 fold changes represent the degree to which each gene changes per change in AGER expression. KEGG pathway enrichment results indicate upregulation of cell-cell communication and extracellular matrix remodeling pathways with the extent of AGER expression (see Additional Files 1 and 2, Supplemental Table 1.6). Ingenuity pathway analysis (IPA) indicated several enriched canonical pathways of genes including hepatic fibrosis/hepatic stellate cell activation (Table 1). Altogether, these data support roles for RAGE in human ALS.
However, as these analyses are from bulk RNA-seq data of cervical spinal cord tissue, it was not possible to discern which cell type(s) may be expressing RAGE and if this too differs across the spectrum of AGER-expressing ALS patients. To address this point, we sought to co-localize RAGE expression with several cell type markers within the subset of patients with the extremes of RAGE expression. We performed a series of immunohistochemistry (IHC) experiments to determine which cell types in the ALS cervical spinal cord were expressing RAGE and if that differed between these ALS patients. We found that there were no differences in microglia (ionized calcium-binding adapter molecule 1, IBA1) or astrocyte (glial fibrillary acidic protein, GFAP) staining area in either the ventral horn gray matter or associated white matter between high and low AGER patients (Fig. 2A–C). However, there were significantly higher levels of RAGE overlap with IBA1 in both areas, but not with GFAP in either area, in high vs. low AGER patients (Fig. 2A, D, E). Furthermore, there was no difference in the overlap of RAGE with a neuronal marker, neuronal nuclei antigen (NeuN), within the ventral horn between high and low AGER patients (Fig. 2A, F). However, additional RAGE signal was detected, which was not localized to any of these cell-type markers and the cellular or extracellular nature remains unclear. These data support the recent spatial transcriptomic analyses that indicated enrichment of the AGE-RAGE pathway in a glia module associated with disease progression in the SOD1G93A murine ALS model and in a corresponding module in human ALS patient tissue . Altogether, these data suggest that AGER expression may be associated with (1) ALS pathology and (2) glial perturbation.
SOD1 G93A mice exhibit increased RAGE-expressing microglia during pathology progression
Based on these analyses in human ALS and the accumulating evidence that microglia are dysfunctional early in the SOD1G93A murine model, which may promote/propagate astrocytic and neuronal dysfunction and that RAGE is known to regulate in vitro microglia-like cell responses to several ligands increased in SOD1G93A mice [25, 33, 34, 70,71,72,73], we sought to examine if microglia expressed RAGE in a prototypic murine model of ALS, the SOD1G93A model. There are currently no mouse models of “sporadic” ALS; however, murine models of familial ALS provide a means to model and test fundamental mechanisms of motor degeneration and reduced survival in these animals vs. their unaffected littermate controls. Hence, we began by performing IHC in SOD1G93A and WT mice (all in the C57BL/6J background) lumbar spinal cord tissue at age 120 days. The overlap area of RAGE with the myeloid marker, integrin subunit alpha m (CD11B), was higher in SOD1G93A mice vs. littermate WT control mice (Fig. 3A, B). While RAGE is expressed in multiple cell-types in spinal cord; the findings that RAGE overlap with microglia is increased in SOD1G93A mice and in a subset of human patients with higher AGER expression, alongside the spatial transcriptomic analyses, implicating the AGE-RAGE pathway alteration in glia, suggested it was logical to probe potential roles for microglia RAGE in ALS-like pathology in the SOD1G93A mouse model.
Microglia Ager deletion extends survival in male SOD1 G93A mice
As accruing evidence has suggested stage-specific roles for microglia in the context of ALS and other neurodegenerative diseases, we sought to assess whether microglia RAGE affects the pathological progression of SOD1G93A mice . We employed a tamoxifen (TAM)-inducible model in which microglia expressed Ager during development and early life and administered TAM to all mice at the age of 90 days to induce Ager deletion in SOD1G93A Agerfl/fl Cx3cr1Cre/+, and SOD1G93A Ager+/+ Cx3cr1Cre/+ mice to evaluate potential roles for RAGE in the progression of pathology (Fig. 4A). We confirmed knock-down of Ager expression in primary CD11B cell isolates from central nervous system (CNS) tissues (Fig. 4B). The area of RAGE overlap with IBA1, but not GFAP, was significantly lower in the SOD1G93A Agerfl/fl Cx3cr1Cre/+ lumbar spinal cord tissue relative to Cre-expressing controls at the end of the study, suggesting Ager knock-down in microglia was maintained (Fig. 4C, D, Supplemental Figure 1A-B). The area of RAGE overlap with the pan neuronal marker microtubule-associated protein 2 (MAP2) at the age of 120 days did not differ between the two genotypes, suggesting that neuronal Ager was not modulated by this approach in this study (Supplemental Figure 1C-D).
As Cx3cr1 Cre-expressing mice, regardless of the administration of TAM, or not, are heterozygous for the Cx3cr1 locus, we evaluated SOD1G93A Ager+/+ Cx3cr1Cre/+ (microglia Ager-expressing, Cx3cr1 Cre-expressing) and SOD1G93A Agerfl/fl Cx3cr1Cre/+ (microglia Ager deficient, Cx3cr1-Cre expressing) mice . As expected, there were no differences in the age of disease onset between any groups, as defined as the age at which the mouse reached maximum weight (Supplemental Figure 2A-B, see Additional file 1; p > 0.05). We found that male SOD1G93A Agerfl/fl Cx3cr1Cre/+ mice displayed significantly longer survival relative to SOD1G93A Ager+/+ Cx3cr1Cre/+ mice, median lifespan of 159 (range 149–175 days) days vs. 152 days (range 136–171 days), respectively (Fig. 5A; p < 0.05). We did not observe differences in lifespan between the female mice groups (Fig. 5B; p > 0.05). Altogether, these data suggested that microglia RAGE may contribute to pathology progression in male SOD1G93A mice.
Male SOD1 G93A Ager fl/fl Cx3cr1 Cre/+ mice displayed reduced pathology during disease progression
We next investigated if microglia Ager contributed to the progressive motor function decline experienced by SOD1G93A mice [4, 5, 38]. Mixed effects analysis of motor function data indicated a significant interaction between time and genotype in male but not female mice (p = 0.038) suggesting that male but not female SOD1G93A Agerfl/fl Cx3cr1Cre/+ mice had time-dependent protection in motor function (Fig. 6A, B). However, there were no significant differences in motor function at any one time point after multiple corrections in either male or female mice (Fig. 6A, B).
As SOD1G93A mice experience pronounced weight loss with disease progression, we next examined if microglia Ager expression modulated this weight loss [4, 5, 38]. Mixed effects analysis of body weight data normalized to weight at tamoxifen administration indicated significant interactions between time and genotype in male and female mice (p = 0.0067, p = 0.0134 respectively) suggesting that in SOD1G93A Agerfl/fl Cx3cr1Cre/+ mice there are time dependent alterations in body weight (Fig. 6C, D). However, there were no significant differences in normalized body-weight measures at any single time point after multiple corrections in male or female mice (Fig. 6C, D). As female mice did not display any significant differences in survival, or motor function, these mice were not investigated further (Figs. 5B and 6B, D).
While no significant differences in motor function or body weight were noted at any single time point, there were time-dependent reductions in the rate of motor function decline and weight in male SOD1G93A Agerfl/fl Cx3cr1Cre/+ mice suggesting a slower rate of disease progression. This led us to consider the possibility that microglia Ager deletion may have affected the numbers of surviving motor neurons within the lumbar ventral horn. We examined the number of surviving motor neurons by IHC at the age of 120 days in male mice and found significantly higher numbers of motor neurons, as labeled by NeuN, DAPI, and choline acetyltransferase (ChAT), in the lumbar ventral horn in the SOD1G93A Agerfl/fl Cx3cr1Cre/+ mice relative to the Cre-recombinase expressing controls (Fig. 6E, F). Altogether, these data suggest that deletion of microglia Ager may reduce/delay neuron death in male SOD1G93A mice. These findings led us to perform additional experiments to provide insight into the mechanisms by which this protection may occur.
Microglia Ager deletion does not impact skeletal muscle macrophages
Next, we considered that deletion of Ager from Cx3cr1 expressing cells at 90 days of age may have modulated skeletal muscle pathology either by direct or indirect effects on skeletal muscle macrophages and/or peripheral monocyte-derived macrophages. We evaluated end-stage and day 120 gastrocnemius muscle for macrophage content, as labeled by F4/80 and CD68, which revealed no differences between groups at either time point (Supplemental Figure 3A-D, 4A-D). Altogether, these data suggest the beneficial effects of microglia Ager reduction were likely restricted to the nervous system.
Male SOD1 G93A Ager fl/fl Cx3cr1 Cre/+ mice lumbar spinal cord exhibits transcriptomic alterations suggestive of improved homeostatic function
Hence, to uncover cell intrinsic and cell-cell communication pathway mechanisms underlying the benefits of microglia Ager deletion in male SOD1G93A mice, we performed RNA-seq on isolated lumbar spinal cord tissue from male SOD1G93A Agerfl/fl Cx3cr1Cre/+ mice and SOD1G93A Ager+/+ Cx3cr1Cre/+ mice at humane endpoint (end-stage). We identified 78 differentially expressed genes between the two genotypes (Fig. 7, Supplemental Table 1.7, see Additional Files 1 and 2). Overrepresentation analysis of the differential gene list indicated enrichment for several KEGG Pathways (Table 2, Supplemental Table 1.8, see Additional Files 1 and 2). Enriched KEGG pathway gene sets included “Cardiac muscle contraction,” “Calcium signaling pathway,” “PPAR signaling pathway,” “Regulation of actin cytoskeleton,” and “Cholesterol metabolism”. Furthermore, canonical pathway analysis using IPA indicated enrichment of several pathways within the differential gene list including “Calcium signaling,” “Actin cytoskeleton signaling,” and “Integrin signaling” among others (Table 2, Supplemental Table 1.9, see Additional Files 1 and 2) . Altogether, these data suggest genotype-dependent effects on cell-cell crosstalk mechanisms, as exemplified by alterations of lipid, metabolic, and integrin signaling gene sets.
To further address this concept, we analyzed potential causal networks that may be modulated in our dataset, which could explain the observed transcriptomic changes . As Ager was deleted solely in Cx3cr1-expressing cells, we limited our analysis to potential cytokines that could originate from microglia and cause transcriptomic alterations across multiple cell types which would be present in the bulk RNA-sequencing data. We identified several significant putative causal networks belonging to numerous cytokine families, including interleukin-1 (IL1), interleukin-3, interferon-α (IFN-α), and interferon-β (IFN-β) (Table 3, Supplemental Table 1.10, see Additional Files 1 and 2).
To evaluate potential roles of microglia RAGE in mediating cytokine expression changes, we turned to an in vitro system. We first conducted a commercially available cytokine screen to evaluate RAGE-dependent effects of RAGE ligand carboxymethyllysine (CML)-AGE, a ligand increased in SOD1G93A spinal cords and in human patient tissues, on BV2 microglia-like cells (Supplemental 5A) [25, 28, 29, 46]. An interesting candidate discovered by this approach was Il1a, as this was also predicted by the RNA-seq causal network analysis (Supplemental Figure 5A, Table 3, Supplemental Table 1.10). In fact, in validation experiments, CML-AGE significantly induced Il1a expression, which was significantly reduced by pre-treatment with a RAGE inhibitor in BV2 cells (Supplemental Figure 5B) . Furthermore, lentiviral transduction of BV2 cells with short hairpin RNA to significantly reduce Ager expression in BV2 cells prevented the CML-AGE associated increase in Il1a expression (Supplemental Figure 5C-D). Interestingly, expression of Malat1, a differentially expressed inflammation-associated lncRNA in the microglia Ager-deficient spinal cord, exhibited a RAGE-dependent increase in BV2 cells in response to CML-AGE (Supplemental Figure 5E, Supplemental Table 1.7) .
Collectively, these results point to significant modulation of intrinsic microglia inflammation. In one or more cell types, prompted by microglia Ager deletion, fundamental changes in general cellular health are observed, as exemplified by alterations in actin cytoskeleton and calcium related pathways, both of which would be predicted to fundamentally alter cell intrinsic and intercellular communication properties in the spinal cord. To verify the implications of these transcriptomic alterations, we examined these points specifically in end-stage tissues.
Microglia Ager deletion reduces the accumulation of damage-associated microglia in male SOD1 G93A mice
The RNA-seq data analysis suggested that microglia Ager deletion resulted in reduction in IL1 and IFN signaling (Table 3), both of which were previously suggested to be dysfunctional in ALS models and in patients [75,76,77,78,79]. In fact, in vitro experiments supported microglia RAGE affecting Il1a and Malat1 expression (Supplemental Figure 3). Beyond the number and density of microglia, their specific gene expression patterns aid in designation of these cells as pro-damage vs. homeostatic phenotype. Specifically, the c-type lectin domain containing 7a protein (CLEC7A) is an established marker of pro-damage/disease-associated microglia, which has been shown to be highly upregulated in SOD1G93A microglia . By IHC, at end-stage, we found that there was a significant reduction in CLEC7A+ area, and CLEC7A+ cell number in SOD1G93A Agerfl/fl Cx3cr1Cre/+ vs. SOD1G93A Ager+/+ Cx3cr1Cre/+ mice (Fig. 8A–C), suggesting that an attenuated or altered DAM phenotype was induced by deletion of microglia Ager. These key findings raise the possibility that RAGE contributes to a cell-intrinsic transition of microglia from a homeostatic to a dysfunctional phenotype.
Microglia Ager deletion reduces accumulation of reactive astrocytes at the end-stage of disease in SOD1 G93A mice
The RNA-seq results indicated reduction of a putative causal network involving IL1 (Table 3). In this context, accumulating evidence suggests that microglia-secreted molecules complement component 1q (C1q), IL1α, and TNF can induce astrocyte reactivity and promote neurotoxicity [70, 71]. Prompted by this consideration, we thus investigated if the SOD1G93A Agerfl/fl Cx3cr1Cre/+ mice displayed alterations in astrocytes. At end-stage, we found that GFAP+ area was significantly lower in SOD1G93A Agerfl/fl Cx3cr1Cre/+ mice relative to Cre-expressing controls (Fig. 9A, B). Recent work has suggested distinct reactive astrocyte phenotypes and identified markers of those states [70, 80]. AMIGO2, adhesion molecule with Ig-like domain 2, was proposed as a marker of “A1” reactive inflammatory astrocytes induced specifically by C1q, IL1α and TNF . Although it is acknowledged that this classification likely does not illuminate the breadth of astrocyte properties and contributions to ALS operative in vivo, we nevertheless examined if the GFAP alterations were concomitant alongside alterations in “A1” astrocytes to begin to define if microglia RAGE might impact astrocyte gene expression. Indeed, we found that GFAP+ AMIGO2+ overlap area was significantly reduced in SOD1G93A Agerfl/fl Cx3cr1Cre/+ mice vs. the Cre-expressing controls (Fig. 9A, C). Altogether, these data suggest that microglia Ager deletion may reduce astrocytic dysfunction likely through altered cytokine expression.
Microglia Ager deletion reduces microgliosis at an earlier stage in male SOD1 G93A mice
We next considered that it was important to also examine microglia and expression of CLEC7A at a time point within the progression phase, but short of end-stage analysis. It is important to note that microglial CLEC7A expression increases over time in this model, and not all microglia may have notable expression of CLEC7A at 120 days of age . Thus, we addressed this point by examining the expression of both CD11B and CLEC7A at day 120. The results indicated significant reductions in the amount of CD11B+ cells and CD11B+ area in SOD1G93A mice devoid of microglia Ager (Fig. 10A–C). Similar to the end-stage tissue analysis, we observed reduced CLEC7A+ area at day 120 (Fig. 10D, E). Altogether, these data suggested that RAGE expression on microglia may contribute to the promotion of microgliosis (Fig. 10A–D), and may affect CLEC7A expression at a stage within the microglia phenotypic transition but prior to frank end-stage tissue pathologies (Figs. 8A–C and 10D, E). As such, we surmised that the overall RAGE-dependent mechanisms in microglia in SOD1G93A mice were likely through alterations in microglia cell intrinsic and cell-cell communication pathways between microglia and other cell-types.
Finally, we considered the possibility that the processes affected by microglia Ager deletion in SOD1G93A mice might relate to the processes accompanied by increasing human AGER expression in ALS patients. Accordingly, we compared the results of both RNA-seq analyses and found that half of the ingenuity canonical pathways identified in the human data analysis overlapped between the two data sets: “Hepatic fibrosis/Hepatic Stellate Cell Activation,” “Atherosclerosis Signaling,” and “Agranulocyte Adhesion and Diapedesis“ (Tables 1 and 2, Supplemental Figure 6-8). Collectively, although it is acknowledged that the human and mouse spinal cord sequencing experiments were not identically designed, these analyses nevertheless suggest that RAGE may modulate extracellular matrix composition, cell-cell communication and lipid metabolism in ALS spinal cord tissues in patients and in SOD1G93A mice.
Our study sought to uncover potential RAGE-dependent roles in ALS by utilizing human patient cervical spinal cord RNA-seq data and the SOD1G93A mouse model of ALS-like pathology. Although RAGE expression was not significantly different between human ALS patients relative to non-neurological controls, this might have been accounted for, in part, by the unexpected observation that ALS patients displayed a range of AGER mRNA in the spinal cord. Accordingly, our analysis revealed upregulation of the “AGE-RAGE pathway in diabetic complications” within cervical spinal cord RNA-seq data of ALS vs. control patients. We found that the amount of AGER negatively correlated with age at onset and age at death or tracheostomy in ALS patients’ cervical spinal cord. If and to what degree the stratification of ALS patients by spinal cord AGER expression may be useful in prediction of prognosis, putative responsiveness to pharmacological interventions and/or the identification of ALS predictive biomarkers, for example, is an important subject for future investigation.
Furthermore, analysis of the consequences of varied degrees of AGER expression across patients’ transcriptomic analyses, using AGER expression as a continuous variable, indicated enrichment in pathways involved in lipid metabolism, extracellular matrix, and cell-cell communication. In fact, microglia displayed increased RAGE protein overlap within both the anterior white matter and the ventral horn in high AGER patients relative to low AGER patients. Of note, our data in human ALS did not suggest sex-dependent differences in AGER-related gene expression patterns. In parallel, the present work in the murine model indicated that there was a significant increase in overlap of RAGE with microglia in SOD1G93A mice at the age of 120 days.
Beyond microglia, other cell-types express RAGE in the spinal cord of both SOD1G93A and human ALS patients. In fact, the roles of RAGE expression in astrocytes are not clear either in homeostasis or in neurodegeneration. RAGE expression in motor neurons has been posited to play a role in the cell-death pathways induced by conditioned media from SOD1G93A-expressing astrocytes [81, 82]. A recent report indicated that global constitutive Ager deletion caused a decrease in survival of SOD1G93A mice, while administering a CNS-permeable RAGE antagonist, FPS-ZM1, starting at 60 days of age, reduced gliosis, and increased motor neuron number at 130 days; but did not alter survival in either sex . However, the authors reported a potential sex-dependent impact on the progression of weight loss and time to death following 10% weight loss . In a separate report, another global constitutive Ager deleted SOD1G93A mouse line exhibited increased survival; however, this Ager deleted mouse line has been reported to harbor a large genomic duplication, which may have confounded the results and explained the discrepancy between the two studies [83, 84]. However, both studies suggested RAGE inhibition reduced gliosis and improved motor function. Treatment of SOD1G93A mice with soluble RAGE (sRAGE), which sequesters RAGE ligands, reduced motor pathology, and extended lifespan in SOD1G93A mice . As there was no significant blood-brain permeability of sRAGE when administered peripherally , the mechanism by which sRAGE exerted protective effects remains elusive. It is possible that circulating soluble factors, modulated by sRAGE treatment, entered the CNS and exerted beneficial effects. Altogether, these studies suggest RAGE could exert complex protective and deleterious effects in a time- and location-dependent manner but cell-types responsible for any of these effects are unclear and warrant further study.
Our findings revealed that deletion of microglia Ager improved survival and motor function in male but not female SOD1G93A mice. In this context, in the mouse model employed in our study, the “knock-in” of the Cre recombinase into the Cx3cr1 locus results in an obligate heterozygous deletion of this gene in all mice, that is fully independent of tamoxifen administration. Thus, all Cre-expressing mice in our study are susceptible to the effects of hemizygous deletion of Cx3cr1. After the initiation of the present work, recent studies have independently implicated Cx3cr1 in the pathogenesis of several neurodegenerative diseases, including ALS. For example, heterozygous deletion of Cx3cr1 induces changes within microglia reminiscent of aged cells [86,87,88]. In fact, in a model of Alzheimer disease (AD), heterozygosity of Cx3cr1 in male mice resulted in reductions of pathology . In contrast, homozygous global deletion of Cx3cr1 in male but not female SOD1G93A mice reduced survival . Altogether, these considerations led us to conclude that it was critical to directly compare our findings in microglia Ager-deleted mice to the Cx3cr1-Cre-control mice to account for the obligate allelic loss of Cx3cr1. As we only observed differences in the male mice, we focused further investigation into these groups to determine what may drive the observed phenotypic differences.
Skeletal muscle macrophage content was not significantly different between the genotypes suggesting the beneficial effects of microglia Ager deletion were likely restricted to the nervous system. Accordingly, transcriptomic analysis of male SOD1G93A Agerfl/fl Cx3cr1Cre/+ and Cre-control lumbar spinal cord tissues indicated microglia Ager deletion reduced gene expression in pathways related to lipid metabolism, actin cytoskeleton, extracellular matrix, and cell-cell communication. Predicted cytokine regulators of these effects included cytokines, such IL1α, C-X-C motif chemokine ligand 1 (CXCL1), C-X-C motif chemokine ligand 9 (CXCL9), TNF superfamily member 10 (TNFSF10), and IFNs. Several of these cytokines are known to be dysregulated in the context of ALS in murine models and in patients [75,76,77]. IFN was recently linked to the acquisition of the damage/disease-associated (DAM) microglia phenotype . DAM are known to accumulate in the SOD1G93A mouse model and are thought to both provoke and ameliorate disease-related pathways [20, 21, 23]. The transition of microglia to this dysfunctional state is theorized to occur in at least several stages with later stages being TREM2-dependent . Critically, CLEC7A, a marker of DAM highly upregulated over time by SOD1G93A microglia , was significantly downregulated via microglia Ager deletion at day 120 and in end-stage tissues. Furthermore, the number of CLEC7A+ microglia reduced by microglia Ager deletion suggested that microglia Ager deletion reduces the activation/acquisition of the DAM phenotype or, alternatively, shifts microglia to a novel phenotype. In fact, BV2 microglia-like cells treated with CML-AGE, a canonical RAGE ligand present in the ALS spinal cord, demonstrated a RAGE-dependent increase in Il1a and the DEG identified by RNA-seq Malat1. Malat1 has been linked to the regulation of pro-inflammatory responses of peripheral macrophages. Malat1-deficient macrophages display reduced LPS responsiveness and altered fibrotic phenotypes . Altogether, these considerations suggest that microglia RAGE may amplify this dysfunctional state, at least in male mice.
Recent work has implicated microglia-sourced IL1α, TNF, and C1q in inducing astrocyte dysfunction and promoting their neurotoxicity in neurodegeneration [70, 90]. In this context, and the considerations that we observe less DAM and that IL1α activity was predicted to be downregulated by microglia Ager deletion, we examined the impact of microglia Ager deletion on the accumulation of GFAP+ “pan” reactive astrocytes as well as the overlap of GFAP with an “A1” reactive astrocyte marker, AMIGO2 [70, 80]. Consistent with the predicted downregulation of IL1α activity in the transcriptomic data and reductions in DAM, we observed lower GFAP+ area and lower overlap area of GFAP+ and AMIGO2+. These data indicate that microglia RAGE signaling may contribute to transformation of astrocytes to a more neurotoxic phenotype.
As we had performed bulk RNA-seq on the murine lumbar spinal cords at sacrifice, we compared the results between the analysis of the ALS patients and the murine dataset. This exercise revealed that half of the ingenuity canonical pathways overlapped between the two analyses and suggests that microglia RAGE may modulate extracellular matrix composition, cell-cell communication, and fatty acid metabolism in ALS and that this may influence disease pathology.
It is important to note that we are unable to distinguish if microglia Ager deletion interacted with potential Cx3cr1-mediated sex-dependent impacts on pathology, or if microglia RAGE has sex-dependent roles in microglia. In support of the latter, male microglia exhibited an enrichment of the RAGE receptor binding Gene Ontology term relative to female microglia isolated from cortical tissues . Furthermore, recent work indicated sex-dependent actions of HMGB1, a known RAGE ligand, on spinal cord microglia in a model of hypersensitivity . Additionally, recent work indicated unexplained sex-dependent effects of FPS-ZM1 treatment on several measures of progression in SOD1G93A mice . Yet, the transcriptomic signature we observed in human ALS patients related to AGER expression was independent of sex. However, the human patient data analyzed were obtained almost entirely from sporadic ALS patients with no known family history of illness; thus, the generalizability of our sex-dependent findings in this model of familial ALS is unclear. Finally, it is also possible that the observed sex-dependent effect in mice is a result of a complex interaction between roles of RAGE within microglia that were altered due to the developmental loss of one allele of Cx3cr1. Further study in distinct models of ALS will be required to discern if microglia RAGE has sex-dependent effects in the pathogenesis of ALS.
Hence, although we may not be able to fully disentangle the effects of the Cx3cr1 heterozygosity on survival from potential roles for RAGE in pathology the SOD1G93A mouse model, the direct comparison to the Cre-control reflects the most rigorous design for this work. Importantly, it is known that monocytes also express Cx3cr1; however, these cells replenish from Cx3cr1 deficient pre-cursors within a few days, and monocyte recruitment into the SOD1G93A spinal cord has been a point of contention within the field. Recent tracing studies have suggested that recruitment is limited until late states of pathology . While we cannot entirely rule out early incipient roles for monocytic RAGE in the context of SOD1G93A pathology, we predict that the largest benefits were from microglia Ager reduction within the CNS. Hence, future studies require testing the RAGE hypothesis in distinct models of ALS, and in a model employing Cre-recombinase mice in which expression of native microglia genes is not affected, and is more specific to microglia, such as the recently developed Tmem119-2A-CreERT2 mouse model .
Altogether, these data demonstrate that microglia RAGE expression may contribute to microglia dysfunction and ultimately to the disruption of communication with multiple cell types including astrocytes and neurons in the diseased ALS spinal cord. Our findings that motor neuron numbers were higher in lumbar spinal cords of male mice devoid of microglia Ager vs. controls at day 120 indicate that such communications emitted from microglia cues were, ultimately, aligned with neuronal survival in the pre-morbid/end-stage state. Collectively, our findings implicate microglia and RAGE in human ALS and provide mechanistic evidence that microglia RAGE contributes to pathology in male SOD1G93A mice.
In summary, our findings suggest microglia RAGE may impact human ALS and contribute to ALS-like pathology in male SOD1G93A mice. Male SOD1G93A mice bearing microglia Ager deficiency from the age of 3 months exhibited reduced gliosis, neuronal and motor function loss, and reduced dysfunctional transcriptomic signatures in lumbar spinal cord tissue. Critically, the observation of a spectrum of AGER and RAGE expression in human ALS spinal cord unveils a myriad of areas for research with respect to implications for ALS vulnerability, severity of disease, and, perhaps, responsiveness to RAGE-directed therapeutics. In this very context, the observed overlap between the human and murine dataset analyses indicate that the pathways altered in mice may be relevant in human patient pathologies. In conclusion, these data provide novel evidence that microglia-specific RAGE expression is a disease-modifying factor in ALS, potentially through modulation of the molecular switch from homeostatic microglia to dysfunctional microglia.
Availability of data and materials
Murine raw FastQ sequencing data and corresponding normalized count data are uploaded for public access on the NCBI GEO database under accession number: GSE160402. All raw Human RNA-seq data from Target ALS samples are publicly available via The Target ALS Multicentered Postmortem Tissue Core, the New York Genome Center for Genomics of Neurodegenerative Disease, Amyotrophic Lateral Sclerosis Association and TOW Foundation (www.targetals.org). Access can be requested by emailing ALSData@nyugenome.org.
Advanced glycation end product
Amyotrophic lateral sclerosis
Adhesion molecule with Ig like domain 2
Complement component 1q
Integrin subunit alpha m
c-Type lectin domain containing 7a protein (or dectin-1)
Central nervous system
C-X-C motif chemokine ligand 1
C-X-C motif chemokine ligand 9
Diaphanous related formin 1
Familial amyotrophic lateral sclerosis
Glial fibrillary acidic protein
High-mobility group box 1
Ionized calcium-binding adapter molecule 1
Ingenuity pathway analysis
- MAP 2:
Microtubule-associated protein 2
Neuronal nuclei antigen
Receptor for advanced glycation end products
Sporadic amyotrophic lateral sclerosis
Superoxide dismutase 1
Triggering receptor expressed on myeloid cells 2
NINDS. Amyotrophic Lateral Sclerosis (ALS) Fact Sheet. Bethesda: NIH Neurological Institute; 2013. No. 16-916
Chen H, Kankel MW, Su SC, Han SWS, Ofengeim D. Exploring the genetics and non-cell autonomous mechanisms underlying ALS/FTLD. Cell Death Differ. 2018;25(4):646–60.
Nguyen DKH, Thombre R, Wang J. Autophagy as a common pathway in amyotrophic lateral sclerosis. Neurosci Lett. 2019;697:34–48. https://doi.org/10.1016/j.neulet.2018.04.006. Epub 2018 Apr 4.
Gurney ME, Pu H, Chiu AY, Dal Canto MC, Polchow CY, Alexander DD, et al. Motor neuron degeneration in mice that express a human Cu,Zn superoxide dismutase mutation. Science. 1994;264(5166):1772–5. https://doi.org/10.1126/science.8209258.
Dal Canto MC, Gurney ME. Neuropathological changes in two lines of mice carrying a transgene for mutant human Cu,Zn SOD, and in mice overexpressing wild type human SOD: a model of familial amyotrophic lateral sclerosis (FALS). Brain Res. 1995;676(1):25–40. https://doi.org/10.1016/0006-8993(95)00063-V.
Haidet-Phillips AM, Hester ME, Miranda CJ, Meyer K, Braun L, Frakes A, et al. Astrocytes from familial and sporadic ALS patients are toxic to motor neurons. Nat Biotechnol. 2011;29(9):824–8. https://doi.org/10.1038/nbt.1957.
Frakes AE, Ferraiuolo L, Haidet-Phillips AM, Schmelzer L, Braun L, Miranda CJ, et al. Microglia induce motor neuron death via the classical NF-kappaB pathway in amyotrophic lateral sclerosis. Neuron. 2014;81(5):1009–23. https://doi.org/10.1016/j.neuron.2014.01.013.
Boillee S, Yamanaka K, Lobsiger CS, Copeland NG, Jenkins NA, Kassiotis G, et al. Onset and progression in inherited ALS determined by motor neurons and microglia. Science. 2006;312(5778):1389–92. https://doi.org/10.1126/science.1123511.
Beers DR, Henkel JS, Xiao Q, Zhao W, Wang J, Yen AA, et al. Wild-type microglia extend survival in PU.1 knockout mice with familial amyotrophic lateral sclerosis. Proc Natl Acad Sci U S A. 2006;103(43):16021–6. https://doi.org/10.1073/pnas.0607423103.
Miller SJ, Glatzer JC, Hsieh YC, Rothstein JD. Cortical astroglia undergo transcriptomic dysregulation in the G93A SOD1 ALS mouse model. J Neurogenet. 2018;32(4):322–35. https://doi.org/10.1080/01677063.2018.1513508.
Philips T, Rothstein JD. Glial cells in amyotrophic lateral sclerosis. Exp Neurol. 2014;262 Pt B:111–20.
Maniatis S, Aijo T, Vickovic S, Braine C, Kang K, Mollbrink A, et al. Spatiotemporal dynamics of molecular pathology in amyotrophic lateral sclerosis. Science. 2019;364(6435):89–93. https://doi.org/10.1126/science.aav9776.
Hardiman O, Al-Chalabi A, Chio A, Corr EM, Logroscino G, Robberecht W, et al. Amyotrophic lateral sclerosis. Nat Rev Dis Prim. 2017;3(1):17085. https://doi.org/10.1038/nrdp.2017.85.
Li Q, Barres BA. Microglia and macrophages in brain homeostasis and disease. Nat Rev Immunol. 2018;18(4):225–42. https://doi.org/10.1038/nri.2017.125.
O'Rourke JG, Bogdanik L, Yanez A, Lall D, Wolf AJ, Muhammad AK, et al. C9orf72 is required for proper macrophage and microglial function in mice. Science. 2016;351(6279):1324–9. https://doi.org/10.1126/science.aaf1064.
Paolicelli RC, Jawaid A, Henstridge CM, Valeri A, Merlini M, Robinson JL, et al. TDP-43 Depletion in microglia promotes amyloid clearance but also induces synapse loss. Neuron. 2017;95(2):297–308 e6. https://doi.org/10.1016/j.neuron.2017.05.037.
Ito Y, Ofengeim D, Najafov A, Das S, Saberi S, Li Y, et al. RIPK1 mediates axonal degeneration by promoting inflammation and necroptosis in ALS. Science. 2016;353(6299):603–8. https://doi.org/10.1126/science.aaf6803.
Bartlett R, Sluyter V, Watson D, Sluyter R, Yerbury JJ. P2X7 antagonism using Brilliant Blue G reduces body weight loss and prolongs survival in female SOD1(G93A) amyotrophic lateral sclerosis mice. PeerJ. 2017;5:e3064. https://doi.org/10.7717/peerj.3064.
Spiller KJ, Restrepo CR, Khan T, Dominique MA, Fang TC, Canter RG, et al. Microglia-mediated recovery from ALS-relevant motor neuron degeneration in a mouse model of TDP-43 proteinopathy. Nat Neurosci. 2018;21(3):329–40. https://doi.org/10.1038/s41593-018-0083-7.
Chiu IM, Morimoto ET, Goodarzi H, Liao JT, O'Keeffe S, Phatnani HP, et al. A neurodegeneration-specific gene-expression signature of acutely isolated microglia from an amyotrophic lateral sclerosis mouse model. Cell Rep. 2013;4(2):385–401. https://doi.org/10.1016/j.celrep.2013.06.018.
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. https://doi.org/10.1016/j.cell.2017.05.018.
Deczkowska A, Keren-Shaul H, Weiner A, Colonna M, Schwartz M, Amit I. Disease-associated microglia: a universal immune sensor of neurodegeneration. Cell. 2018;173(5):1073–81. https://doi.org/10.1016/j.cell.2018.05.003.
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–81 e9. https://doi.org/10.1016/j.immuni.2017.08.008.
MacLean M, Derk J, Ruiz HH, Juranek JK, Ramasamy R, Schmidt AM. The receptor for advanced glycation end products (RAGE) and DIAPH1: implications for vascular and neuroinflammatory dysfunction in disorders of the central nervous system. Neurochem Int. 2019;126:154–64. https://doi.org/10.1016/j.neuint.2019.03.012.
Juranek JK, Daffu GK, Wojtkiewicz J, Lacomis D, Kofler J, Schmidt AM. Receptor for advanced glycation end products and its inflammatory ligands are upregulated in amyotrophic lateral sclerosis. Front Cell Neurosci. 2015;9:485.
Casula M, Iyer AM, Spliet WG, Anink JJ, Steentjes K, Sta M, et al. Toll-like receptor signaling in amyotrophic lateral sclerosis spinal cord tissue. Neuroscience. 2011;179:233–43. https://doi.org/10.1016/j.neuroscience.2011.02.001.
Lo Coco D, Veglianese P, Allievi E, Bendotti C. Distribution and cellular localization of high mobility group box protein 1 (HMGB1) in the spinal cord of a transgenic mouse model of ALS. Neurosci Lett. 2007;412(1):73–7. https://doi.org/10.1016/j.neulet.2006.10.063.
Kikuchi S, Shinpo K, Ogata A, Tsuji S, Takeuchi M, Makita Z, et al. Detection of N epsilon-(carboxymethyl)lysine (CML) and non-CML advanced glycation end-products in the anterior horn of amyotrophic lateral sclerosis spinal cord. Amyotroph Lateral Scler Other Motor Neuron Disord. 2002;3(2):63–8. https://doi.org/10.1080/146608202760196020.
Kaufmann E, Boehm BO, Sussmuth SD, Kientsch-Engel R, Sperfeld A, Ludolph AC, et al. The advanced glycation end-product N epsilon-(carboxymethyl)lysine level is elevated in cerebrospinal fluid of patients with amyotrophic lateral sclerosis. Neurosci Lett. 2004;371(2-3):226–9. https://doi.org/10.1016/j.neulet.2004.08.071.
Yeh CH, Sturgis L, Haidacher J, Zhang XN, Sherwood SJ, Bjercke RJ, et al. Requirement for p38 and p44/p42 mitogen-activated protein kinases in RAGE-mediated nuclear factor-kappaB transcriptional activation and cytokine secretion. Diabetes. 2001;50(6):1495–504. https://doi.org/10.2337/diabetes.50.6.1495.
Shen C, Ma Y, Zeng Z, Yin Q, Hong Y, Hou X, et al. RAGE-specific inhibitor FPS-ZM1 attenuates AGEs-induced neuroinflammation and oxidative stress in rat primary microglia. Neurochem Res. 2017;42(10):2902–11. https://doi.org/10.1007/s11064-017-2321-x.
Manigrasso MB, Pan J, Rai V, Zhang J, Reverdatto S, Quadri N, et al. Small molecule inhibition of ligand-stimulated RAGE-DIAPH1 signal transduction. Sci Rep. 2016;6(1):22450. https://doi.org/10.1038/srep22450.
Bianchi R, Kastrisianaki E, Giambanco I, Donato R. S100B protein stimulates microglia migration via RAGE-dependent up-regulation of chemokine expression and release. J Biol Chem. 2011;286(9):7214–26. https://doi.org/10.1074/jbc.M110.169342.
Bianchi R, Giambanco I, Donato R. S100B/RAGE-dependent activation of microglia via NF-kappaB and AP-1 Co-regulation of COX-2 expression by S100B, IL-1beta and TNF-alpha. Neurobiol Aging. 2010;31(4):665–77. https://doi.org/10.1016/j.neurobiolaging.2008.05.017.
Hudson BI, Kalea AZ, Del Mar AM, Harja E, Boulanger E, D'Agati V, et al. Interaction of the RAGE cytoplasmic domain with diaphanous-1 is required for ligand-stimulated cellular migration through activation of Rac1 and Cdc42. J Biol Chem. 2008;283(49):34457–68. https://doi.org/10.1074/jbc.M801465200.
Rai V, Maldonado AY, Burz DS, Reverdatto S, Yan SF, Schmidt AM, et al. Signal transduction in receptor for advanced glycation end products (RAGE): solution structure of C-terminal rage (ctRAGE) and its binding to mDia1. J Biol Chem. 2012;287(7):5133–44. https://doi.org/10.1074/jbc.M111.277731.
Parkhurst CN, Yang G, Ninan I, Savas JN, Yates JR 3rd, Lafaille JJ, et al. Microglia promote learning-dependent synapse formation through brain-derived neurotrophic factor. Cell. 2013;155(7):1596–609. https://doi.org/10.1016/j.cell.2013.11.030.
Wooley CM, Sher RB, Kale A, Frankel WN, Cox GA, Seburn KL. Gait analysis detects early changes in transgenic SOD1(G93A) mice. Muscle Nerve. 2005;32(1):43–50. https://doi.org/10.1002/mus.20228.
Hurtado del Pozo C, Ruiz HH, Arivazhagan L, Aranda JF, Shim C, Daya P, et al. A Receptor of the immunoglobulin superfamily regulates adaptive thermogenesis. Cell Rep. 2019;28(3):773–91.e7.
Scott S, Kranz JE, Cole J, Lincecum JM, Thompson K, Kelly N, et al. Design, power, and interpretation of studies in the standard murine model of ALS. Amyotroph Lateral Scler. 2008;9(1):4–15. https://doi.org/10.1080/17482960701856300.
Juranek JK, Daffu GK, Geddis MS, Li H, Rosario R, Kaplan BJ, et al. Soluble RAGE treatment delays progression of amyotrophic lateral sclerosis in SOD1 mice. Front Cell Neurosci. 2016;10:117.
Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9(7):676–82. https://doi.org/10.1038/nmeth.2019.
Otsu N. A threshold selection method from Gray-level histograms. IEEE Trans Syst Man Cybern. 1979;9(1):62–6. https://doi.org/10.1109/TSMC.1979.4310076.
Tsai WH. Moment-preserving thresholding—a new approach. Comput Vision Graph. 1985;29(3):377–93. https://doi.org/10.1016/0734-189X(85)90133-1.
Zack GW, Rogers WE, Latt SA. Automatic measurement of sister chromatid exchange frequency. J Histochem Cytochem. 1977;25(7):741–53. https://doi.org/10.1177/25.7.70454.
Senatus L, López-Díez R, Egaña-Gorroño L, Liu J, Hu J, Daffu G, Li Q, Rahman K, Vengrenyuk Y, Barrett TJ, Dewan MZ, Guo L, Fuller D, Finn AV, Virmani R, Li H, Friedman RA, Fisher EA, Ramasamy R, Schmidt AM. RAGE impairs murine diabetic atherosclerosis regression and implicates IRF7 in macrophage inflammation and cholesterol metabolism. JCI Insight. 2020;5(13):e137289. https://doi.org/10.1172/jci.insight.137289.
Wingett S, Andrews S. FastQ Screen: a tool for multi-genome mapping and quality control [version 2; peer review: 4 approved]. F1000Research. 2018;7:1338.
Andrews S. FastQC: a quality control tool for high throughput sequence data; 2010.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. https://doi.org/10.1093/bioinformatics/btu170.
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. https://doi.org/10.1093/bioinformatics/bts635.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30. https://doi.org/10.1093/bioinformatics/btt656.
Liao Y, Smyth GK, Shi W. The Subread aligner: fast, accurate and scalable read mapping by seed-and-vote. Nucleic Acids Res. 2013;41(10):e108. https://doi.org/10.1093/nar/gkt214.
Frankish A, Diekhans M, Ferreira AM, Johnson R, Jungreis I, Loveland J, et al. GENCODE reference annotation for the human and mouse genomes. Nucleic Acids Res. 2019;47(D1):D766–D73. https://doi.org/10.1093/nar/gky955.
Robinson MD, Oshlack A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010;11(3):R25. https://doi.org/10.1186/gb-2010-11-3-r25.
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. https://doi.org/10.1093/bioinformatics/btp616.
Team RC. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018.
Nikolayeva O, Robinson MD. edgeR for differential RNA-seq and ChIP-seq analysis: an application to stem cell biology. Methods Mol Biol. 2014;1150:45–79. https://doi.org/10.1007/978-1-4939-0512-6_3.
Kanehisa M, Sato Y, Furumichi M, Morishima K, Tanabe M. New approach for understanding genome variations in KEGG. Nucleic Acids Res. 2019;47(D1):D590–D5. https://doi.org/10.1093/nar/gky962.
Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004;32(Database issue):D277–80. https://doi.org/10.1093/nar/gkh063.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30. https://doi.org/10.1093/nar/28.1.27.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7. https://doi.org/10.1089/omi.2011.0118.
Kramer A, Green J, Pollard J Jr, Tugendreich S. Causal analysis approaches in ingenuity pathway analysis. Bioinformatics. 2014;30(4):523–30. https://doi.org/10.1093/bioinformatics/btt703.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300. https://doi.org/10.1111/j.2517-6161.1995.tb02031.x.
Lun AT, Chen Y, Smyth GK. It's DE-licious: a recipe for differential expression analyses of RNA-seq experiments using quasi-likelihood methods in edgeR. Methods Mol Biol. 2016;1418:391–416. https://doi.org/10.1007/978-1-4939-3578-9_19.
Wu D, Smyth GK. Camera: a competitive gene set test accounting for inter-gene correlation. Nucleic Acids Res. 2012;40(17):e133. https://doi.org/10.1093/nar/gks461.
Wu D, Lim E, Vaillant F, Asselin-Labat ML, Visvader JE, Smyth GK. ROAST: rotation gene set tests for complex microarray experiments. Bioinformatics. 2010;26(17):2176–82. https://doi.org/10.1093/bioinformatics/btq401.
The Gene Ontology C. The gene ontology resource: 20 years and still GOing strong. Nucleic Acids Res. 2019;47(D1):D330–D8.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat Genet. 2000;25(1):25–9. https://doi.org/10.1038/75556.
Pinheiro J BD, DebRoy S, Sarkar D, R Core Team. nlme: linear and nonlinear mixed effects models. 2020.
Liddelow SA, Guttenplan KA, Clarke LE, Bennett FC, Bohlen CJ, Schirmer L, et al. Neurotoxic reactive astrocytes are induced by activated microglia. Nature. 2017;541(7638):481–7. https://doi.org/10.1038/nature21029.
Joshi AU, Minhas PS, Liddelow SA, Haileselassie B, Andreasson KI, Dorn GW 2nd, et al. Fragmented mitochondria released from microglia trigger A1 astrocytic response and propagate inflammatory neurodegeneration. Nat Neurosci. 2019;22(10):1635–48. https://doi.org/10.1038/s41593-019-0486-0.
Chen J, Sun Z, Jin M, Tu Y, Wang S, Yang X, et al. Inhibition of AGEs/RAGE/Rho/ROCK pathway suppresses non-specific neuroinflammation by regulating BV2 microglial M1/M2 polarization through the NF-kappaB pathway. J Neuroimmunol. 2017;305:108–14. https://doi.org/10.1016/j.jneuroim.2017.02.010.
Serrano A, Donno C, Giannetti S, Peric M, Andjus P, D'Ambrosi N, et al. The astrocytic S100B protein with its receptor RAGE is aberrantly expressed in SOD1(G93A) models, and its inhibition decreases the expression of proinflammatory genes. Mediat Inflamm. 2017;2017:1626204.
Cui H, Banerjee S, Guo S, Xie N, Ge J, Jiang D, Zörnig M, Thannickal VJ, Liu G. Long noncoding RNA Malat1 regulates differential activation of macrophages and response to lung injury. JCI Insight. 2019;4(4):e124522. https://doi.org/10.1172/jci.insight.124522.
Hensley K, Floyd RA, Gordon B, Mou S, Pye QN, Stewart C, et al. Temporal patterns of cytokine and apoptosis-related gene expression in spinal cords of the G93A-SOD1 mouse model of amyotrophic lateral sclerosis. J Neurochem. 2002;82(2):365–74. https://doi.org/10.1046/j.1471-4159.2002.00968.x.
Meissner F, Molawi K, Zychlinsky A. Mutant superoxide dismutase 1-induced IL-1beta accelerates ALS pathogenesis. Proc Natl Acad Sci U S A. 2010;107(29):13046–50. https://doi.org/10.1073/pnas.1002396107.
Hu Y, Cao C, Qin XY, Yu Y, Yuan J, Zhao Y, et al. Increased peripheral blood inflammatory cytokine levels in amyotrophic lateral sclerosis: a meta-analysis study. Sci Rep. 2017;7(1):9094. https://doi.org/10.1038/s41598-017-09097-1.
Ye J, Cheung J, Gerbino V, Ahlsen G, Zimanyi C, Hirsh D, et al. Effects of ALS-associated TANK binding kinase 1 mutations on protein-protein interactions and kinase activity. Proc Natl Acad Sci U S A. 2019;116(49):24517–26. https://doi.org/10.1073/pnas.1915732116.
Liu J, Gao L, Zang D. Elevated Levels of IFN-gamma in CSF and serum of patients with amyotrophic lateral sclerosis. PLoS One. 2015;10(9):e0136937. https://doi.org/10.1371/journal.pone.0136937.
Zamanian JL, Xu L, Foo LC, Nouri N, Zhou L, Giffard RG, et al. Genomic analysis of reactive astrogliosis. J Neurosci. 2012;32(18):6391–410. https://doi.org/10.1523/JNEUROSCI.6221-11.2012.
Kim MJ, Vargas MR, Harlan BA, Killoy KM, Ball LE, Comte-Walters S, Gooz M, Yamamoto Y, Beckman JS, Barbeito L, Pehar M. Nitration and Glycation Turn Mature NGF into a Toxic Factor for Motor Neurons: A Role for p75NTR and RAGE Signaling in ALS. Antioxid Redox Signal. 2018;28(18):1587–602. https://doi.org/10.1089/ars.2016.6966. Epub 2017 Jun 26.
Liu L, Killoy KM, Vargas MR, Yamamoto Y, Pehar M. Effects of RAGE inhibition on the progression of the disease in hSOD1(G93A) ALS mice. Pharmacol Res Perspect. 2020;8(4):e00636. https://doi.org/10.1002/prp2.636.
Lee JD, McDonald TS, Fung JNT, Woodruff TM. Absence of receptor for advanced glycation end product (RAGE) reduces inflammation and extends survival in the hSOD1(G93A) mouse model of amyotrophic lateral sclerosis. Mol Neurobiol. 2020;57(10):4143–55. https://doi.org/10.1007/s12035-020-02019-9.
Bartling B, Zunkel K, Al-Robaiy S, Dehghani F, Simm A. Gene doubling increases glyoxalase 1 expression in RAGE knockout mice. Biochim Biophys Acta Gen Subj. 1864;2020(1):129438.
Deane R, Du Yan S, Submamaryan RK, LaRue B, Jovanovic S, Hogg E, et al. RAGE mediates amyloid-beta peptide transport across the blood-brain barrier and accumulation in brain. Nat Med. 2003;9(7):907–13. https://doi.org/10.1038/nm890.
Liu C, Hong K, Chen H, Niu Y, Duan W, Liu Y, et al. Evidence for a protective role of the CX3CL1/CX3CR1 axis in a model of amyotrophic lateral sclerosis. Biol Chem. 2019;400(5):651–61. https://doi.org/10.1515/hsz-2018-0204.
Hickman SE, Allison EK, Coleman U, Kingery-Gallagher ND, El Khoury J. Heterozygous CX3CR1 deficiency in microglia restores neuronal beta-amyloid clearance pathways and slows progression of Alzheimer's like-disease in PS1-APP mice. Front Immunol. 2019;10:2780. https://doi.org/10.3389/fimmu.2019.02780.
Gyoneva S, Hosur R, Gosselin D, Zhang B, Ouyang Z, Cotleur AC, Peterson M, Allaire N, Challa R, Cullen P, Roberts C, Miao K, Reynolds TL, Glass CK, Burkly L, Ransohoff RM. Cx3cr1-deficient microglia exhibit a premature aging transcriptome. Life Sci Alliance. 2019;2(6):e201900453. https://doi.org/10.26508/lsa.201900453.
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. https://doi.org/10.1172/JCI133737.
Guttenplan KA, Weigel MK, Adler DI, Couthouis J, Liddelow SA, Gitler AD, et al. Knockout of reactive astrocyte activating factors slows disease progression in an ALS mouse model. Nat Commun. 2020;11(1):3753. https://doi.org/10.1038/s41467-020-17514-9.
Guneykaya D, Ivanov A, Hernandez DP, Haage V, Wojtas B, Meyer N, et al. Transcriptional and translational differences of microglia from male and female brains. Cell Rep. 2018;24(10):2773–83 e6. https://doi.org/10.1016/j.celrep.2018.08.001.
Agalave NM, Rudjito R, Farinotti AB, Khoonsari PE, Sandor K, Nomura Y, et al. Sex-dependent role of microglia in disulfide high mobility group box 1 protein-mediated mechanical hypersensitivity. Pain. 2021;162(2):446–58. https://doi.org/10.1097/j.pain.0000000000002033.
Chiot A, Zaidi S, Iltis C, Ribon M, Berriat F, Schiaffino L, et al. Modifying macrophages at the periphery has the capacity to change microglial reactivity and to extend ALS survival. Nat Neurosci. 2020;23(11):1339–51. https://doi.org/10.1038/s41593-020-00718-z.
Kaiser T, Feng G. Tmem119-EGFP and Tmem119-CreERT2 Transgenic Mice for Labeling and Manipulating Microglia. eNeuro. 2019;6(4):ENEURO.0448-18.2019. https://doi.org/10.1523/ENEURO.0448-18.2019.
The authors gratefully acknowledge the expert assistance of Ms. Latoya Woods of the Diabetes Research Program, Department of Medicine, NYU Grossman School of Medicine, in the preparation of this manuscript. We thank the NYU Grossman School of Medicine’s Genome Technology Center; which is partially funded by the Cancer Center Support grant P30CA016087 at NYU Langone Health’s Cancer Center. We acknowledge the Target ALS postmortem tissue core for providing us with human cervical spinal cord sections. We acknowledge the Target ALS Multicentered Postmortem Tissue Core, the New York Genome Center for Genomics of Neurodegenerative Disease, Amyotrophic Lateral Sclerosis Association and TOW Foundation for access to the human RNA-seq data and corresponding metadata.
Funding for this project included: Department of Defense award: W81XWH-17-1-0346 to A.M.S; F31NS120424 from the National Institutes of Health National Institute of Neurological Disorders and Stroke to M.M; and support from the Diabetes Research Program, NYU Grossman School of Medicine.
Ethics approval and consent to participate
All mouse experiments were performed under protocols approved by the New York University Grossman School of Medicine Institutional Animal Care and Use Committee (IACUC) in accordance with international and NIH guidelines. All human data, and tissue used in this study was provided by TargetALS in a de-identified manner.
Consent for publication
A.M.S is an inventor on awarded (US 9,353,3078, US 9364472) and pending (US 16/094270, Europe 177864360.0) patent applications assigned to NYU Grossman Medical Center.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1.
Figures and figure legend for Supplemental Figure 1–8. Also, table legends for Supplemental Tables 1.1–1.10.
Additional file 2.
Excel workbook containing Supplemental Tables 1.1–1.10.
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.
About this article
Cite this article
MacLean, M., Juranek, J., Cuddapah, S. et al. Microglia RAGE exacerbates the progression of neurodegeneration within the SOD1G93A murine model of amyotrophic lateral sclerosis in a sex-dependent manner. J Neuroinflammation 18, 139 (2021). https://doi.org/10.1186/s12974-021-02191-2