A novel phosphoproteomic landscape evoked in response to type I interferon in the brain and in glial cells
Journal of Neuroinflammation volume 18, Article number: 237 (2021)
Type I interferons (IFN-I) are key responders to central nervous system infection and injury and are also increased in common neurodegenerative diseases. Their effects are primarily mediated via transcriptional regulation of several hundred interferon-regulated genes. In addition, IFN-I activate several kinases including members of the MAPK and PI3K families. Yet, how changes to the global protein phosphoproteome contribute to the cellular response to IFN-I is unknown.
The cerebral phosphoproteome of mice with brain-targeted chronic production of the IFN-I, IFN-α, was obtained. Changes in phosphorylation were analyzed by ontology and pathway analysis and kinase enrichment predictions. These were verified by phenotypic analysis, immunohistochemistry and immunoblots. In addition, primary murine microglia and astrocytes, the brain's primary IFN-I-responding cells, were acutely treated with IFN-α and the global phosphoproteome was similarly analyzed.
We identified widespread protein phosphorylation as a novel mechanism by which IFN-I mediate their effects. In our mouse model for IFN-I-induced neurodegeneration, protein phosphorylation, rather than the proteome, aligned with the clinical hallmarks and pathological outcome, including impaired development, motor dysfunction and seizures. In vitro experiments revealed extensive and rapid IFN-I-induced protein phosphorylation in microglia and astrocytes. Response to acute IFN-I stimulation was independent of gene expression and mediated by a small number of kinase families. The changes in the phosphoproteome affected a diverse range of cellular processes and functional analysis suggested that this response induced an immediate reactive state and prepared cells for subsequent transcriptional responses.
Our studies reveal a hitherto unappreciated role for changes in the protein phosphorylation landscape in cellular responses to IFN-I and thus provide insights for novel diagnostic and therapeutic strategies for neurological diseases caused by IFN-I.
Inflammation of the central nervous system (CNS) is a fundamental response to a wide range of stimuli ranging from virus infection to tumors and neurodegenerative diseases. A key component of the inflammatory response is the type I interferons (IFN-I), which have both protective and detrimental effects . Importantly, increasing evidence points to a fundamental role for IFN-I in ageing and neurodegeneration [2,3,4]. In a recent study, Roy and colleagues (2020) demonstrated that IFN-I drive neuroinflammation in Alzheimer’s disease (AD) and mediate synapse loss . At its extreme, IFN-I-mediated CNS inflammation and neurodegeneration are known as “cerebral type I interferonopathies” [5, 6]. They include genetic diseases (e.g., Aicardi–Goutières Syndrome, AGS), chronic and congenital infections (e.g., HIV, toxoplasma, and cytomegalovirus) and autoinflammatory disorders (e.g., neurological manifestations of systemic lupus erythematosus, SLE) [6,7,8]. Importantly, the mechanisms of how IFN-I drive neurological diseases are unclear, making causal treatment difficult.
The IFN-I are a large family of cytokines that includes multiple IFN-α subtypes and a single IFN-β . The cellular effects of IFN-I are driven through the coordinated regulation of the expression of a large number of genes known as interferon regulated genes (IRGs). All IFN-I signal through the Janus kinase (JAK)/signal transducer and activator of transcription (STAT) pathway . Canonical activation of the pathway leads to the phosphorylation of the transcription factors STAT1 and STAT2 that form a trimolecular complex with interferon regulatory factor 9 (IRF9). This complex, termed interferon stimulated gene factor 3, is transcriptionally active and is primarily responsible for mediating the expression of IRGs [10,11,12]. However, the canonical JAK/STAT pathway may only represent the tip of the iceberg when it comes to the cellular mechanisms by which IFN-I actions are mediated. The existence of alternative IFN-I-activated signaling pathways via mitogen-activated protein kinase (MAPK) and phosphatidylinositol 3-kinase [9, 13] suggests that we understand very little of the complexity of the IFN-I system and that IFN-I may use additional, as yet undefined mechanisms to mediate cellular responses.
Advancements in mass spectrometry-based techniques provide an excellent opportunity to uncover novel networks of protein interaction and regulation on a system-wide level. Using a global phosphoproteomic approach, we determined the global phosphoproteome in the CNS of transgenic mice with CNS-restricted production of IFN-α, termed GIFN39 mice [14, 15]. These mice recapitulate the key clinical and pathological changes seen in patients with cerebral type I interferonopathies . We discovered that the changes in the phosphorylation landscape of the CNS in transgenic mice at early-stage disease aligned with disease features seen in mice with late-stage disease. The changes were largely driven by a small number of kinase families including MAPKs, cyclin-dependent kinases (CDKs), casein kinases (CKs) and calcium calmodulin kinases (CaMKs). Furthermore, the phosphoproteome of acutely IFN-α-treated microglia and astrocytes—the primary immune responding cells of the CNS—revealed an immediate reactive state that prepared these cells for the subsequent transcriptomic response. Functional comparison of the in vivo and in vitro phosphoproteomes suggested that microglia and astrocytes are prominent contributors to disease in cerebral type I interferonopathies. In summary, here we identify a novel mechanism for the regulation of the cellular response to IFN-I. Importantly, this mechanism is largely independent of altered transcription and may serve as a novel indicator for diagnostic and therapeutic purposes targeting particular kinases in a cell-type-specific manner.
Transgenic GIFN39 mice (originally obtained from the Scripps Research Institute, La Jolla, CA, USA, where they were developed by I. L. Campbell) are described previously [14, 15] and were bred and maintained in-house. Mice were maintained under specific-pathogen free conditions at the animal facility of the University of Sydney. For in vitro experiments, cells were isolated from cortices of 1–4-day old wildtype (WT) C57Bl/6 mice. For ex vivo experiments, 8- and 16-week-old mice were euthanized with isoflurane, measured, brains weighed and the cerebellum dissected and flash frozen. For histological analyses, brains were fixed in 10% neutral buffered formalin overnight and then paraffin-embedded.
Primary microglia and astrocyte purification
Mixed glial cultures and subsequent purification into microglia and astrocytes were prepared as described previously . Briefly, single-cell suspensions of papain digested cortices were plated on poly-D-lysine coated 75 cm2 flasks and incubated at 37 °C in a humidified incubator with 5% CO2 until confluent in media (Dulbecco's modified Eagle medium with high glucose, 10% fetal bovine serum, 100 U/ml penicillin and 10 μg/ml streptomycin). Microglia were collected by mechanical shaking. Astrocytes were purified by CD11b-negative selection following manufacturer’s instructions (130-049-601, Miltenyi Biotec, Macquarie Park, NSW, Australia). Flow cytometry with anti-CD11b (101211, Biolegend, San Diego, CA, USA) was used to determine astrocyte purity after sorting with an allowance of 0.5% CD11b+ contaminating cells. Purified microglia from three 75 cm2 flasks were pooled into one 75 cm2 flask and purified astrocytes were seeded at six million cells per 75 cm2 flask.
IFN-α treatment of primary microglia and astrocytes
At 80–90% confluence, microglia and astrocytes were washed with media and incubated in media with or without 1,000 U/ml IFN-α11 (130-093-130, Miltenyi Biotec) at 37 °C in a humidified incubator with 5% CO2. Treatment was stopped with two washes of ice-cold PBS and protein extracted in lysis buffer.
For both tissue and cells, phosphopeptides and non-phosphorylated peptides were obtained from the same animals or cell samples as outlined in Additional file 1: Fig. S1. For all samples, protein was extracted in urea lysis buffer (6 M urea (U6504, Sigma-Aldrich, Castle Hill, NSW, Australia), 2 M thiourea (88810, Sigma-Aldrich), 20 U/ml Benzonase (E1014, Sigma-Aldrich) with 1 × Protease Inhibitor Cocktail III (539,134, Merck, Bayswater, VIC, Australia) and 1 × Phosphatase Inhibitor Cocktail II (524,625, Merck) by hand-held homogeniser (Omni, Kennesaw GA). For microglia, protein from ten 75 cm2 flasks, and for astrocytes, protein from two 75 cm2 flasks, was pooled by methanol–chloroform precipitation and resuspended in urea lysis buffer. Cysteine residues were reduced in the presence of 1,4-dithiothreitol and subsequently alkylated with iodoacetamide. Following tryptic digestion, 250 µg of peptide was labelled with one of four isobaric tags (iTRAQ, Sciex, Mt Waverly, VIC, Australia; TMT, Thermo Fisher Scientific, Waltham, MA) as listed in Table S1. Phosphopeptides were enriched according to [17, 18], producing three peptide populations (singly and multiply phosphorylated and non-phosphorylated). Singly phosphorylated and non-phosphorylated peptides were further fractionated using offline hydrophilic interaction chromatography generating 10 fractions each. Identification and quantification were performed on an LTQ-Orbitrap (QE plus, ThermoFisher Scientific, North Ryde, NSW, Australia) mass spectrometer in data dependent mode. All experiments were performed in duplicate. Data were analyzed using Proteome Discoverer (Version 2.2, ThermoFisher Scientific) and searched using an in house MASCOT server (Version 2.4) against the UniProt Mus musculus database (database version July 2018) with the following parameters: maximum 2 missed cleavages, 20 ppm mass error (MS) and 0.2 Da mass error (MS/MS); iTRAQ was searched as a fixed modification; variable modifications included phosphorylation (Ser, Thr, Tyr), carbamidomethyl (Cys), oxidation (Met), acetylation (protein N-term) and cyclisation (Glu and Asp). False discovery rate of 0.01 was applied for non-phosphorylated peptides, which was relaxed to 0.05 for singly and multiply phosphorylated peptides, to allow for the labile nature of the modification. All experiments were performed in duplicate with moderate to strong correlation (Additional file 1: Fig. S2). Normalization of iTRAQ/TMT reporters was calculated using the sum of all intensities approach across all phosphorylated and non-phosphorylated peptide spectral matches (PSM) within a given experiment, prior to calculation of ratios using WT or 0 min IFN-α stimulated as denominator. Log2 ratios and population wide mean and standard deviation were calculated for each PSM using a total reporter ion intensity-dependent ranking for z-score calculation, adapted from . The weighted average was used to compare median z-score and fold change across duplicate labelling experiments for each phosphosite and phosphopeptide (phosphoproteome), and protein (non-phosphorylated only). For analysis of proteome, significance was defined as a median z-score ≥ 1.96 or ≤ − 1.96 (non-phosphorylated, 95% confidence interval) and ≥ 1 or ≤ − 1 (68.3% confidence interval) for phosphopeptides and phosphosites. A minimum of five PSMs per protein was used a further confidence threshold for regulated non-phosphorylated proteins. Regulation was defined with a fold change ≥ 1.5 or ≤ − 1.5.
Relative abundance of proteins were estimated by their NSAF as described in . Ontology enrichment of biological processes (GOTERM_BP_DIRECT) and pathways (KEGG_PATHWAY) was performed with Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 [21, 22]. Canonical pathways, upstream analysis and disease functions were analyzed by Ingenuity Pathway Analysis (IPA) (QIAGEN Inc., https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis) . Kinase prediction analyses utilized motif-X [24, 25] with changed parameters of foreground and background of the IPI mouse proteome and a motif length of 15; in vivo group-based prediction system (v1.0)  and PhosphoSitePlus, March 2017 . The activity score of “activity, induced” or “activity, inhibited” for each kinase was calculated by calculating the sum of the multiplication of the log2 fold change of each significant PhosphoSitePlus-annotated phosphosites on the kinase with its PSM and + 1 if activity was induced or − 1 if activity was inhibited and then divided by the sum of the PSMs of those phosphosites. This score reflects the cumulative effect of the phosphosites on the kinase to an induced or inhibited state compared with the basal status weighted on fold change and normalized to abundance. A threshold score of 0.53 reflects a log2 of a 1.5-fold change in activity and if positive, an increase, and if negative, a reduction in activity. Proteins were considered as kinases as defined in . Heatmaps, principle component analysis (PCA) and some plots were generated in R  using packages gplots , ggplot2  and VennDiagram  and the princomp function. Summary of gene ontologies were extracted using REVIGO . Sequences of mouse and human STATs were aligned with CLUSTAL O (1.2.4)  and annotated using information from PhosphoSitePlus with focus on transcription induction, inhibition, alteration and intracellular localization. IRGs were identified using the Interferome database .
Paraffin Sections (5 μm) were deparaffinized and rehydrated in graded ethanol. Antigens were unmasked with 25 mM Tris pH 9 (for GFAP), 25 mM Tris pH 8 and 0.05% SDS (for Iba1) or 10 mM citrate pH 6 in 0.05% Tween-20 (for CD3, neurofilament and cleaved caspase-3) in a vegetable steamer for 45 min. Sections were incubated in 0.03% H2O2 for 10 min and blocked in 1% goat serum with 0.1% Triton-X and 0.05% Tween-20 in PBS for 20 min. Primary antibodies against GFAP (Z0334, Agilent, Santa Clara, CA, USA; 1:1,000 dilution), Iba1 (019–19,741, Wako Pure Chemical Industries, Osaka, Japan; 1:1,000 dilution), CD3 (ab16669, Abcam, Melbourne, VIC, Australia; 1:200 dilution), cleaved-caspase-3 (9661, Cell Signaling Technologies; 1:300 dilution) and neurofilament (N0142, Sigma-Aldrich; 1:200 dilution) were incubated overnight at 4 °C. Following washing, sections were incubated with biotinylated secondary anti-rabbit and anti-mouse antibodies (BA-1000 and BA-2000, Vector Laboratories, Burlingame, CA, USA) followed by VECTASTAIN Elite ABC HRP Kit (PK-7200, Vector Laboratories). Sections were developed with 3,3′-diaminobenzidine with or without Ni enhancement (SK-4100, Vector Laboratories), prepared for dual stains or counterstained with Mayer’s Hematoxylin and mounted. For alizarin red S (ARS) post-immunostaining, sections were incubated with 2% ARS pH 4.2, destained in acetone, washed in 1:1 acetone:xylene and in xylene and then a cover slip was applied. Sections were view with the DM4000B microscope (Leica Macquarie Park, NSW, Australia) and imaged with Spot Imaging Software (Spot Imaging, Sterling Heights, MI, USA). For Iba1-ARS visualization of microglia and calcification, images were acquired at various focus depths and blended in Adobe Photoshop. Brightness and contrast were globally adjusted using Adobe Photoshop.
Total protein for immunoblotting was isolated from the cerebellum of mice or treated microglia and astrocytes in 50 mM Tris pH 7.5, 150 mM NaCl, 1 mM EDTA, 1% sodium deoxycholate, 1% Triton-X 100, 0.1% SDS and 1 × Protease Inhibitor Cocktail III and 1 × Phosphatase Inhibitor Cocktail II. Proteins were separated on a 10% Tris–glycine acrylamide gel, transferred onto polyvinylidene fluoride membranes, blocked with 5% skim milk or 5% bovine serum albumin in 137 mM NaCl, 20 mM Tris–HCl (pH 7.4) and 0.1%(v/v) Tween-20 and probed using primary and secondary antibodies. Primary antibodies used were pY701-STAT1 (7649; 1:2,000 dilution), pS727-STAT1 (9177; 1:2,000 dilution), STAT1 (9172; 1:2,000 dilution), STAT2 (4597; 1:1,000 dilution), pY705-STAT3 (9131; 1:2,000 dilution), STAT3 (4904; 1:2,000 dilution), pY694-STAT5 (9351; 1:250 dilution), pT202/204-MAPK1/3 (9101; 1:2,000 dilution), pT180/pY182-MAPK14 (9211; 1:1,000 dilution), MAPK14 (9212; 1:1,000 dilution), pT183/pY185-MAPK8/9 (9255; 1:1,000 dilution), MAPK8/9 (9252; 1:1,000 dilution), pS473-AKT (4058; 1:1,000 dilution), AKT (4691; 1:2,000 dilution) (all from Cell Signaling Technologies, Danvers, MA, USA), pY689-STAT2 (07–224; 1:1,000 dilution, Merck Millipore), GAPDH (MAB374, Merck Millipore; 1:30,000 dilution), MAPK1/3 (M5670, Sigma-Aldrich; 1:40,000 dilution) and STAT5 (sc835, Santa Cruz Biotechnology, Dallas, Texas, USA; 1:200 dilution). Secondary antibodies used were horseradish peroxidase conjugated anti-rabbit (sc2004, Santa Cruz; 1:30,000 dilution) and peroxidase conjugated anti-mouse (A0168, Sigma-Aldrich; 1:10,000 dilution). Chemiluminescence (WBKLS0500, Merck Millipore) was detected in the ChemiDoc XRS + (BioRad Laboratories, Gladesville, NSW, Australia) and relative proteins levels were quantified by densitometry using Fiji . To reduce and account for variations across membranes, membranes were incubated together with the same solutions and imaged at the same time, while identical samples were included on all membranes and normalized to each other. Source data available in (Additional file 2).
Experimental design and statistical analysis
Both genders of mice were analyzed when available in the colony and there was no previously observed effect of gender in these mice with outcome of disease [14, 15]. Only female mice were used for phosphoproteomic analyses to reduce chance of variability with mixed genders. Samples sizes were based on previous studies and are listed in figure legends. Statistical analyses of different data sets are outlined in figure legends and were performed in GraphPad Prism, version 8 (La Jolla, CA, USA). Briefly, comparisons of survival curves were made by log-rank test and comparisons between genotypes and age were made by two-way ANOVA with Tukey’s post-test. Normality of residuals and homogeneity of variance were based on QQ plots and homoscedasticity plots, respectively. A P ≤ 0.05 was considered significant.
Chronic IFN-α signaling in vivo upregulates proteins that drive an inflammatory response
To determine global changes in protein levels and changes in protein phosphorylation in the CNS of GIFN39 vs WT mice, we isolated proteins from the cerebella and performed mass spectrometry on nonenriched (proteome) and enriched phosphopeptides (phosphoproteome) as described in material and methods. The cerebellum was chosen as it has the highest level of transgene expression in the CNS . Eight-week-old GIFN39 mice were used as they show only mild inflammatory changes and minimal signs of neurodegeneration, focusing our analysis on primarily IFN-α-driven responses rather than the compounded effects from leukocyte infiltration and neurodegeneration. A total of 5,616 proteins were detected with one protein less abundant and 70 proteins more abundant (|z-score|≥ 1.96 and |fold change|≥ 1.5) in the proteome of the cerebellum of GIFN39 mice compared with WT mice (Fig. 1A and Additional file 1: Tables S1 and S2). The protein with reduced abundance was α-crystallin B chain (CRYAB) and is a non-classical chaperone preventing protein aggregation in stress responses . Yet, its role in neurodegenerative and neuroinflammatory diseases is undetermined. Of the 71 proteins, 82% are regulated by IFN-I, extrapolated using the Interferome database. The other proteins are antibody subunits and endopeptidase inhibitors. Using IPA, enriched canonical pathways consisting of the proteins with elevated abundance in the cerebellum of GIFN39 mice were associated with interferon signaling, antigen presentation, acute phase response and interestingly, but less significant, SLE signaling (Fig. 1B).
Extensive protein phosphorylation in GIFN39 mice gives deeper insights into disease processes
Analysis of the phosphoproteome of cerebella from WT and GIFN39 mice revealed 17,609 phosphosites on 4,748 proteins and a strong correlation between replicates (Additional file 1: Fig. S2 and Table S1). Comparing the phosphoproteome from GIFN39 with WT mice, 303 phosphosites on 203 proteins showed significantly increased phosphorylation and 397 phosphosites on 295 proteins decreased phosphorylation (Fig. 2A). Of these regulated phosphoproteins, 22 proteins were significantly increased in abundance, while one had reduced abundance. To dissect the contribution of increased and/or decreased phosphorylation in biological processes, functional analysis was performed on each regulated set of proteins and compared using the DAVID. A large number of processes were associated with decreased phosphorylation (Fig. 2B) reflecting the greater number of dephosphorylated phosphosites. Proteins with increased phosphorylation were associated with responses to virus and immune processes, whereas proteins with reduced phosphorylation were largely associated with neuron-related functions and endocytosis (Fig. 2C). Top processes with increased and decreased phosphorylation included regulation of the cytoskeleton and cell–cell interactions, such as the endothelial barrier. Thus, while the ontologies of the regulated proteome were associated primarily with immune responses, the regulated phosphoproteome reflected both immune responses and additional non-immune functions. Furthermore, the functional analysis of the phosphoproteome provided a greater insight into the regulated processes than the proteome.
The CNS-phosphoproteome of GIFN39 mice reflects their phenotype and disease outcome
Next, we further investigated the regulated global phosphoproteome in GIFN39 mice using IPA. By grouping phosphoproteins into protein–protein networks, IPA determined the enrichment of diseases and disorders associated with the networks and calculated an activation score. The significantly activated diseases and disorders included mortality, motor dysfunction, emotional behavior, seizures, cell death and anxiety-like behavior (Fig. 3A). By contrast, processes associated with neuronal function, such as long-term potentiation and neurite branching, were significantly reduced.
To validate these enriched diseases and disorders of the 8-week-old GIFN39 mouse, we investigated their clinical and pathological outcomes at eight and 16 weeks of age. GIFN39 mice had a significantly reduced survival compared with WT littermates (Fig. 3B) and were smaller in weight and length at both ages while remaining the same weight across 8 weeks (Fig. 3C, D). In addition, the average brain weight of GIFN39 mice was less at 16 weeks of age compared with WT mice (Fig. 3E), suggesting brain atrophy. Furthermore, at 16 weeks of age, GIFN39 mice had moderate to severe ataxia, displayed wild running and jumping and occasional convulsive seizures, as previously described [14, 15].
Pathologically, GIFN39 mice developed perivascular calcification, where some microglia were observed wrapping the deposits (Fig. 3F). Microglia were also observed to be highly ramified with increased Iba1 staining compared with WT mice, at both ages. There was increased GFAP immunoreactivity in the cerebellum of GIFN39 mice (Fig. 3G), which was also reflected in the proteomic data (Additional file 1: Table S2) and by immunoblot (Additional file 1: Fig. S3). A progressive increase in CD3+ T cell infiltrates was observed with age in the GIFN39 mice, while T cells were not observed in the cerebellum of WT mice (Fig. 3H). There was reduced neurofilament staining of Purkinje cells and dendrites in the molecular layer in 16-week-old GIFN39 mice compared with WT mice indicating neurodegeneration (Fig. 3I). In addition, there was staining for cleaved caspase-3 only observed in 16-week-old GIFN39 mice (Fig. 3J) revealing increased apoptosis. Large vascular lumens were observed in the brains of GIFN39 mice at 16 weeks of age indicating the presence of aneurysms (Fig. 3J). Together, the clinical and pathological observations in aged GIFN39 mice validate those identified by the phosphoproteome at early-stage disease.
We next probed for the phosphorylation status of various STATs as a proxy to changes in the global phosphoproteome with age. Canonical IFN-I signaling relies on Y701-STAT1 phosphorylation, which was significantly elevated compared with the cerebellum of WT mice at both ages (Additional file 1: Fig. S3). In addition, pS727-STAT1, which is considered to modulate STAT1 activity [38, 39], and total STAT1 were significantly increased in abundance and ratio in GIFN39 mice compared with WT mice. There was no detectable increase in STAT2 activation (pY689-STAT2) but total STAT2 was increased in GIFN39 mice compared with WT mice. In addition to STAT1 and STAT2, STAT3 and STAT5 have also been implicated to contribute to IFN-I signaling [9, 13]. Accordingly, STAT3 was activated (pY705-STAT3) in the cerebellum of GIFN39 mice and increased with age. However, activation of STAT5 was not detected. Thus, these findings show chronic IFN-I signaling in the brains of GIFN39 mice and further activation of alternative signaling pathways with age, possibly associated with neurodegeneration.
Extensive protein phosphorylation in GIFN39 mice is associated with only a few kinase families
The results suggest changes in the phosphoproteome is involved the progression of the disease. We next investigated which kinases drive the phosphorylation of proteins using three in silico approaches. There are no established strategies for identifying phosphatases, which mediate protein dephosphorylation, largely due to dissimilarities of substrate-motifs of targets . The first approach to identify kinases was substrate-centric, whereby enriched peptide motifs with increased phosphorylation were extracted with motif-X. Enriched motifs were matched to consensus kinase-substrate motifs to identify active kinase families. Four serine-centric motifs were enriched and associated with the MAPK/CDK (sP and sP…K), CK (s.E) and CaMK (R..s) families (Fig. 4A). The MAPK and CDK families share a common proline-directed phosphorylation site , resulting in both families matching to similar motifs. The second approach predicted kinase families based on their association to substrate motifs and the protein–protein interactome. The number of substrates linked to each kinase family was normalized to the number of phosphorylated peptides associated with a kinase family to give the percentage of phosphopeptides targeted by a kinase family. Similar to the first approach, the MAPK family was predicted to phosphorylate the largest percentage of phosphosites in the cerebellum of GIFN39 mice, followed by CDK and AKT (Fig. 4B). In addition, 27% of phosphosites were not matched to a kinase. The third approach was kinase-centric, calculating a kinase activation score. This score estimates the overall activity of a kinase from the cumulative effect of each phosphosite in relation to its abundance and effect on the kinase—whether it increased or inhibited activity, extracted from PhosphoSitePlus. Interestingly, there was a slight increase in MAPK14 activity and a slight decrease in activity from CaMK2A, CaMK2B and mTOR (Fig. 4C). However, a large number of phosphosites did not have annotations resulting in only 28 of all 250 kinases identified in our data being assigned an activation score. Although there are limitations to each approach including variable kinase-substrate specificities, multiple kinases targeting the same protein and most phosphosites having undefined roles , the three approaches identified the same few core kinase families, with the MAPK family having the most prominent presence in the regulation of the cerebellar phosphoproteome in GIFN39 mice. In support of this, immunoblots for MAPK1/3 and MAPK14 showed phosphorylation at their canonical activation sites in the cerebellum of GIFN39 mice, while phosphorylation of MAPK8/9 was not detected (Additional file 1: Fig. S3). The phosphosites for MAPK1/3 and MAPK14 were also detected in the phosphoproteome indicating basal activation. Furthermore, there was slight increased MAPK3 activation in GIFN39 mice compared with WT mice.
The proteome of primary microglia and astrocytes aligns with their expected biological roles
To determine, the contribution of specific cell types to the phosphoproteome changes in the CNS of GIFN39 mice, we next analyzed the proteome and phosphoproteome of microglia and astrocytes, which are dominant IFN-I responding CNS-resident cells [16, 43]. The phosphoproteome is sensitive to signaling with rapid global changes occurring within minutes [44,45,46,47,48]. This excluded the isolation and ex vivo analysis of cell types from the brain of mice. Hence, we investigated the response to IFN-α in vitro in primary murine microglia and astrocytes. We treated microglia or astrocytes with IFN-α for 0, 5, 15 and 30 min to determine their direct cellular response to IFN-α. The short treatment ensured a largely transcription-independent response in the absence of secondary signaling from IRGs or secreted cytokines and non-physiological responses attributed to a cultured system.
Proteomic post-hoc analysis showed high purity of our microglia and astrocyte cultures, based on the relative protein abundance of cell-type-specific markers (Additional file 1: Fig. S4A, B). Comparing the proteomes (Table S1), 3,017 proteins were common between microglia and astrocytes, while 743 and 1,240 proteins were specific to microglia and astrocytes, respectively. The biological processes and pathways associated with the common proteins mapped to fundamental cellular processes including gene transcription, translation and metabolism (Table S3). Microglia-specific proteins mapped to their known roles as primary resident immune cells: immune-related processes and inflammation. For astrocytes, the top processes and pathways included cell–cell contact, cell size and brain development. Furthermore, top ten transcriptional regulators common to both cell types were associated with homeostasis, metabolism and proliferation (Table S4). The top ten microglia-associated transcriptional regulators were involved with myeloid lineage, immune function and homeostasis, whereas astrocyte-associated transcriptional regulators were linked to cell migration and proliferation, and homeostasis. Thus, common and cell-type-specific transcriptional regulators align well with the top biological processes and pathways of the total proteome further reflecting the cell-type-specific roles of these cells in the CNS.
IFN-α induces extensive and cell-type-specific changes of the phosphoproteome in microglia and astrocytes
Investigating the IFN-α response in microglia and astrocytes by PCA of the fold changes of all phosphosites revealed distinct changes in the phosphoproteome across treatment duration and between cell types (Fig. 5A). In IFN-α-treated microglia, the overall change to the phosphoproteome was greater than in astrocytes (Additional file 1: Fig. S5). In microglia, after 5 min of IFN-α treatment, 1,480 sites on 817 proteins had significantly increased phosphorylation levels and 1,091 sites on 726 proteins had reduced phosphorylation levels, when compared with untreated cells (|z-score|≥ 1 and |fold change|≥ 1.5). Numbers decreased after 15 min before increasing again after 30 min treatment. By contrast, in astrocytes, 109 sites on 99 proteins had increased phosphorylation and 348 sites on 295 proteins had decreased phosphorylation levels after 5 min treatment. The number increased further after 15 and 30 min IFN-α treatment. In addition, although tyrosine phosphorylation is key in mediating signaling of IFN-I, only 1.2–2.6% of regulated residues were tyrosine, with greater regulation in microglia compared with astrocytes (Fig. 5B). The other two residues, serine and threonine, were most abundantly regulated. This reflects the abundance of each of these residues in the proteome determined in previous large-scale phosphoproteomic studies [42, 46]. This demonstrated that microglia and astrocytes responded rapidly and extensively to IFN-α and that the response in microglia was more wide-ranging than in astrocytes.
Next, we analyzed the functional correlates associated with the IFN-α-regulated phosphoproteome using IPA. The top 15 significant pathways in microglia were similar across all treatment timepoints and were associated with immunity, actin dynamics, proliferation, survival and apoptosis and various intracellular signal transduction and metabolism pathways (Fig. 5C). In contrast to microglia, the astrocyte response showed greater variability over time and included cell membrane interactions and remodeling, translation, proliferation, survival and apoptosis, and actin dynamics (Fig. 5D).
Microglia and astrocytes show distinct IFN-α-induced phosphorylation dynamics of the JAK/STAT pathway
To further investigate the differential signaling kinetics in microglia and astrocytes, we focused on the JAK/STAT pathway. Due to the nature of mass spectrometry, not all phosphosites can be identified owing to detection limits of peptide length, peptide physico-chemical properties, peptide abundance and use of data-dependent acquisition [49,50,51]. Nevertheless, a comprehensive and unbiased survey of the global phosphoproteome is achievable. Expectedly, pY701-STAT1, an activation marker of IFN-I signaling, was markedly increased in both cell types following IFN-α treatment (Fig. 6). pS727-STAT1 levels similarly increased with treatment time but with delayed kinetics and a lower magnitude aligning with its proposed role of modulating the activity of STAT1 [38, 39]. Together, IFN-α treatment induced the expected increase of the canonical markers of IFN-I activation with time [52,53,54] revealing the two-phase regulation of the microglia global phosphoproteome was not an artifact.
Interestingly, additional phosphosites (T699, T704 and S735) were regulated on STAT1 by IFN-α (Figs. 6A and Additional file 1: Fig. S6). These sites lie in the transactivation domain suggesting regulation of transcriptional activity of STAT1 and that pT704 may be involved in dimerization of STATs . Extrapolation of the function of these sites from annotated STATs indicate a role for modulating activity and protein-complex formation (Additional file 1: Fig. S7A). Following Y701 phosphorylation, STAT1 dimerizes with tyrosine-phosphorylated STAT2 and IRF9 to regulate gene expression. While no phosphosites were detected by mass spectrometry on STAT2, likely owing to above mentioned reasons, phosphorylation of the canonical Y689 was identified by immunoblot (Fig. 6B–D). By contrast, multiple sites were phosphorylated on IRF9, detected in the GIFN39 cerebellum and IFN-α-treated cells (Fig. 6A, S6 and S7B), which may modulate transcription though protein-complex formation as the residues lie in the C-terminal IRF-associated domain (Additional file 1: Fig. S7B). In addition to STAT1 and STAT2, phosphoproteomic analysis identified activation of STAT3, STAT5a, STAT5b and STAT6 in microglia and STAT5a and STAT5b in astrocytes in response to IFN-α (Fig. 6A). Of these, tyrosine phosphorylation of pY705-STAT3 and pY689-STAT5 was confirmed by immunoblot (Fig. 6B–D). Together, these findings reveal dichotomous signaling kinetics between microglia and astrocytes in the acute response to IFN-α. Furthermore, the data suggest novel ways in which canonical IFN-I signaling is modulated between these cell types.
Diversity of IFN-α-regulated phosphosites is modulated by a limited number of kinases
We next investigated the contributing kinases to the increased phosphorylation measured by IFN-α treatment. To identify kinases associated with the IFN-α-induced changes in phosphorylation, we used the same three in silico approaches as above. In microglia, the enriched motifs matched to substrates of MAPK/CDK, CK, CaMK and AKT (Fig. 7A). One third of enriched motifs had no clear match to a kinase. While fewer enriched motifs were identified in the treated astrocytes, the same five kinase families were matched (Fig. 7B), indicating similar core IFN-I signaling pathways. In support, the second approach similarly identified MAPK and CDK families as the key kinases contributing to IFN-α-induced phosphorylation in microglia and astrocytes (Fig. 7C). Other top predicted kinase families were CaMK, glycogen synthase kinase (GSK) and AKT. The third approach identified MAPK14 and its downstream target MAPKAPK2 to have 1.5-fold increased activity in IFN-α-treated microglia (Fig. 7D). By contrast, IFN-α treatment in astrocytes reduced activity of MAPK14. However, activation of members of the MAPK family was detected by immunoblots in both cell types following IFN-α treatment resulting in a slight but non-significant increase in MAPK14 phosphorylation (Additional file 1: Fig. S8). Moreover, there was a significantly higher abundance of these kinases in microglia compared with astrocytes (Additional file 1: Fig. S8). In addition, activity of PRKCD and CDK1 was reduced after 30 min of IFN-α treatment in microglia and astrocytes, respectively. Together, the MAPK family was a major contributor to IFN-α-induced protein phosphorylation in both microglia and astrocytes. This supports the notion seen in vivo of a core IFN-α response being largely driven by these kinases, while differences in phosphorylation kinetics and cellular response between the two cell types reflects the changes in activity of each kinase.
Astrocytes rather than microglia play a diverse role in the regulated global phosphoproteome in GIFN39 mice
We next extrapolated the contribution of direct IFN-α signaling in microglia and astrocytes in the disease to GIFN39 mice. Of note, the presence of neurons and other cells in vivo likely reduces detection of microglia- and astrocyte-specific proteomes in the cerebellum. In addition, other cytokines and mediators of inflammation in the inflamed brains of GIFN39 mice would also modulate the detected phosphoproteome. Furthermore, to investigate the contribution of IFN-α signaling by cell type in vivo would be difficult as the phosphoproteome would rapidly change during cell isolation. Hence, we restricted our comparison to functional ontologies. Although, acutely-treated cells will likely differ in their response compared with microglia and astrocytes in the brains of GIFN39 mice, comparison of the diseases and disorders in treated microglia and astrocytes suggested that IFN-α signaling in microglia and astrocytes contributed to the top diseases and disorders enriched in the regulated cerebellar phosphoproteome of GIFN39 mice (Fig. 8A). Delineating the functional contribution, the enriched biological processes of each of the regulated phosphoproteomes were compared. Overall, 46% of processes enriched in the phosphoproteome of GIFN39 mice were also enriched in IFN-α-treated cells (Fig. 8B). To summarize ontology terms, semantics were used to combine related terms : in common, was nervous system development including establishment of the endothelial intestinal barrier, cytoskeleton organization, protein autophosphorylation, endocytosis, positive regulation of protein binding, intracellular signal transduction, cell–cell adhesion and negative regulation of neuron apoptotic processes. Processes contributed by microglia include response to virus, positive regulation of microtubule polymerization, receptor mediated endocytosis and protein dephosphorylation. Processes contributed by astrocytes include eukaryotic translation initiation factor 4F complex assembly, exocytosis, substantia nigra development and endoplasmic reticulum calcium ion homeostasis. Expectedly, neuron-associated processes were exclusive to the cerebella phosphoproteome of GIFN39 mice. Together, these results indicate the cell-type-specific contribution from the direct IFN-α signaling in astrocyte, rather than microglia, have a more diverse contribution to the changes in the cerebellar phosphoproteome in GIFN39 mice.
Type I interferons have both beneficial and detrimental actions in the CNS ranging from protection against infections to mediating inflammation and neurodegeneration. Although the wide spectrum of actions of IFN-I are known to involve IRGs, here we identify a novel mechanism by which IFN-I regulate the cellular response that is independent of transcription. We demonstrate that in addition to regulating gene expression, IFN-I induce widespread changes in protein phosphorylation in vivo. These changes preceded severe disease in a mouse model for cerebral type I interferonopathies and—importantly—reflected the cardinal features of the clinical and pathological disease. Furthermore, microglia and astrocytes are key cellular components of this phosphoproteome response. Finally, in silico analysis revealed that a limited number of kinase families mediated a large proportion of the protein phosphorylation, identifying them as potential targets for future therapeutic intervention of IFN-I-mediated neuroinflammatory diseases.
Chronic production of IFN-α in young GIFN39 mice resulted in altered phosphorylation of nearly 500 proteins despite the absence of pronounced disease or overt neurodegeneration. Strikingly, the phosphoproteome accurately reflected the development of pathological features and disease outcome in GIFN39 mice. Functional analysis revealed enriched disease categories in the phosphoproteome of GIFN39 mice that are also present in humans with cerebral type I interferonopathies [56,57,58], suggesting that widespread protein phosphorylation also contributes to these diseases. Accordingly, the total number of proteins showing differential phosphorylation in GIFN39 vs WT mice was comparable to findings in mouse models of AD and several human degenerative and inflammatory diseases including AD and SLE [59,60,61]. While increased phosphorylation in the GIFN39 cerebellum was broadly associated with a response to virus infection, reduced phosphorylation was linked to neuronal functions. A similar observation has been made in AD  and various other neurodegenerative diseases [62, 63] suggesting that dysregulated protein phosphorylation plays a significant role in a wide range of neurodegenerative and neuroinflammatory diseases.
In contrast to the widespread changes in protein phosphorylation, only 71 proteins showed altered abundance in cerebella from 8-week-old GIFN39 mice compared with WT mice. This low number is probably due to a combination of factors, including the mice not yet showing prominent disease, and the relative high proportion of neurons in the cerebellum  which have a very limited response to IFN-I . Of the 71 proteins identified in our study several are associated with vascular aneurysms , such as fibronectin, which are key pathological features observed in GIFN39 mice, patients with cerebral type I interferonopathies and during IFN-I treatment [57, 66]. In addition, antigen presentation and autoimmunity were prominently featured, corroborating studies that have identified autoantibodies in type I interferonopathies [67, 68]. Of note, similar to blood mononuclear cells from patients with SLE [61, 69], a systemic type I interferonopathy, the diversity of biological processes was limited in the total proteome compared with the phosphoproteome, demonstrating the importance of posttranslational modifications for regulating cellular responses. Nevertheless, changes in the proteome—like those of the phosphoproteome—were highly specific and correlated well with key pathological features of GIFN39 mice and patients with type I interferonopathies, alike.
Microglia and astrocytes are highly responsive to IFN-I, regulating a significantly greater number of IRGs when compared with neurons [16, 43]. Here we found that acute treatment of microglia and astrocytes in vitro resulted in extensive changes in protein phosphorylation. In addition to a common IFN-I response, both cell types showed cell-type-specific responses that were reflective of their biological roles in vivo. In microglia, immunity-associated intracellular signaling was enriched, paralleling recent transcriptomic findings  and further highlighting their role as the principal immune cell of the CNS. Interestingly, SLE signaling was only enriched in microglia, suggesting a role for microglia in the CNS manifestations of this disease. Accordingly, microglia are a major source of IFN-I in CNS-lupus . By contrast, top enriched pathways in IFN-α-treated astrocytes underscored the role of these cells in blood–brain barrier regulation and integrity, and contribution to leukocyte extravasation. Accordingly, loss of IFN-I signaling in astrocytes increases blood brain barrier permeability during infection with a neurotropic virus . Together, this suggests that microglia and astrocytes are prominent contributors to the global response seen in the CNS in vivo.
The changes in protein phosphorylation in microglia and astrocytes occurred within 5 to 30 min and were thus largely independent of de novo gene expression. Furthermore, several signaling pathways enriched in the acute phosphoproteome were also identified in the transcriptome of these cells in response to 12 h treatment with IFN-α , suggesting that the acute global phosphorylation-dependent response to IFN-α guides the transcriptomic response. Thus, as phosphorylation preceded transcription and translation upon IFN-α stimulation, changes in protein phosphorylation likely prime microglia and astrocytes and instigate rapid host protection responses in a transcription-independent manner.
To address the molecular basis underlying the increased protein phosphorylation we performed in silico analysis of kinase motifs and activities. Unexpectedly, these studies revealed that only a limited set of kinases belonging to the MAPK/CDK, CK, CaMK and AKT families likely mediated the majority of the observed phosphorylation events in vivo and in vitro. These kinase families have previously been associated with the cellular response to IFN-α [53, 72]. It is tempting to speculate that differential activities and abundances of members of the MAPK family may contribute to the lower expression of several IRGs in astrocytes compared with microglia in response to IFN-α . Of note, while mass spectroscopy and immunoblots revealed no obvious differences in phosphorylation of these kinases apart from MAPK3 and MAPK14, it is conceivable that additional post translational modifications, such as acetylation of MAPK14 , also contribute to increased kinase activity. A role for the MAPK pathway in regulating cellular effects to IFN-α is further supported by the presence of S727 phosphorylation of STAT1 in GIFN39 mice and IFN-α-treated cells, which is mediated by MAPK14 [74,75,76]. It will be interesting to determine if different IFN-I subtypes activate the same kinases. While all IFN-I share the same components of the canonical signalling pathway, they bind to the receptor with different kinetics and affinities [77, 78]. It is conceivable that this could result in differences in downstream kinase activation and IFN-I-subtype specific responses . Our studies also identified a number of novel phosphosites in the canonical IFN-I signaling molecules STAT1 and IRF9. In addition to the critical Y701-STAT1, three phosphosites on STAT1 previously not reported to be IFN-α-regulated, as well as the known site S727, were seen in IFN-α-treated cells. Interestingly, T699, T704 and S735 phosphorylation was cell-type specific, suggesting that these sites may contribute to the distinct responses in microglia and astrocytes. Since these residues lie within the transactivation domain of STAT1 it is conceivable that they regulate transcriptional activity either directly or through recruitment of co-activators. Accordingly, crystal structures indicate that T704-STAT1 could be involved in STAT1 dimerization . Similarly, IRF9 was phosphorylated at multiple sites. Unlike other IRF family members, IRF9 is constitutively active  and does not require phosphorylation of the C-terminal region. However, based on sequence similarities the identified phosphosites on IRF9 could indicate the presence of a so far unappreciated signaling network that modulates interferon stimulated gene factor 3 activity.
Importantly, our results further suggest a clear dichotomy in IFN-I signaling, where the changes in the cellular phosphoproteome are mediated mostly by members of the MAPK/CDK, CK, CaMK and AKT kinase families, whereas transcriptomic changes are largely dependent on canonical JAK/STAT signaling. These findings open up new opportunities for the development of kinase-specific drugs as a mean to modulate IFN-I-mediated neuroinflammation and neurodegeneration. Accordingly, in recent years, several kinase inhibitors have been approved for the treatment of a range of degenerative diseases, immune disorders and cancers [81, 82]. For example, JAK inhibitors (e.g., baricitinib, tofacitinib) have been used to treat patients with type I interferonopathies , while LRRK2 and CDK5 inhibitors are being explored as treatments for various neurodegenerative disorders . While targeting the MAPK pathway has been less successful due to toxicity , identifying specific members, such as MAPK14, above may allow for the development of more specific inhibitors with less side effects.
This study provides evidence for a novel mechanism for IFN-I action that involves so far unappreciated and extensive changes in the protein phosphorylation landscape affecting a diverse range of cellular processes. The changes in the phosphoproteome not only preceded disease in a mouse model for cerebral type I interferonopathy but aligned with the clinical and pathological outcome. In microglia and astrocytes, IFN-α-induced protein phosphorylation was rapid and occurred in a cell-type-specific manner highlighting the biological roles of both cell types in the CNS and their potential contribution to IFN-I-induced disease. In addition to preparing the cells for the subsequent transcriptomic response, the changes in the phosphoproteome induced a reactive state in the cells. The changes are mediated, in parallel to the canonical IFN-I JAK/STAT signaling pathway, by a limited number of kinase families. These studies reveal a hitherto unappreciated role for global protein phosphorylation in response to IFN-I and provide novel insight for the development of diagnostic and therapeutic strategies for diseases of the CNS and periphery caused by IFN-I.
Availability of data and materials
The mass spectrometry proteomics data have been deposited to Zenodo (http://doi.org/10.5281/zenodo.4422055 and http://doi.org/10.5281/zenodo.4433003) and to the ProteomeXchange Consortium via the Proteomics Identification (PRIDE)  partner repository with the data set identifier PXD014443 and http://doi.org/10.6019/PXD014443. Total quantified proteins, phosphopeptides and phosphosites of WT vs GIFN39 cerebella and IFN-α-treated microglia and astrocytes (Table S1) has been deposited to Zenodo (https://doi.org/10.5281/zenodo.4287788).
Analysis of variance
Alizarin Red S
Calcium calmodulin kinase
Central nervous system
α-Crystallin B chain
Database for Annotation, Visualization and Integrated Discovery
Glyceraldehyde 3-phosphate dehydrogenase
Glial fibrillary acidic protein
Glycogen synthase kinase
Ionized calcium-binding adapter molecule 1
Type I interferons
Ingenuity Pathway Analysis
Interferon regulator factor
Interferon regulated genes
Isobaric tags for relative and absolute quantitation
Kyoto Encyclopedia of Genes and Genomes
Mitogen-activated protein kinase
MAPK-activated protein kinase 2
Normalized spectral abundance factor
Principle component analysis
Protein kinase C delta type
Peptide spectral match
Systemic lupus erythematosus
Signal transducer and activator of transcription
Paul S, Ricour C, Sommereyns C, Sorgeloos F, Michiels T. Type I interferon response in the central nervous system. Biochimie. 2007;89(6–7):770–8. https://doi.org/10.1016/j.biochi.2007.02.009.
Baruch K, Deczkowska A, David E, Castellano JM, Miller O, Kertser A, et al. Aging. Aging-induced type I interferon response at the choroid plexus negatively affects brain function. Science. 2014; 346(6205): 89–93. Doi: https://doi.org/10.1126/science.1252945
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 Investig. 2020;130(4):1912–30. https://doi.org/10.1172/jci133737.
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. https://doi.org/10.1016/j.neurobiolaging.2013.10.089.
Crow YJ. Type I interferonopathies: mendelian type I interferon up-regulation. Curr Opin Immunol. 2015;32:7–12. https://doi.org/10.1016/j.coi.2014.10.005.
McGlasson S, Jury A, Jackson A, Hunt D. Type I interferon dysregulation and neurological disease. Nat Rev Neurol. 2015;11(9):515–23. https://doi.org/10.1038/nrneurol.2015.143.
Crow YJ, Manel N. Aicardi-Goutieres syndrome and the type I interferonopathies. Nat Rev Immunol. 2015;15(7):429–40. https://doi.org/10.1038/nri3850.
Hofer MJ, Campbell IL. Type I interferon in neurological disease-the devil from within. Cytokine Growth Factor Rev. 2013;24(3):257–67. https://doi.org/10.1016/j.cytogfr.2013.03.006.
Gonzalez-Navajas JM, Lee J, David M, Raz E. Immunomodulatory functions of type I interferons. Nat Rev Immunol. 2012;12(2):125–35. https://doi.org/10.1038/nri3133.
Majoros A, Platanitis E, Kernbauer-Holzl E, Rosebrock F, Muller M, Decker T. Canonical and non-canonical aspects of JAK-STAT signaling: lessons from interferons for cytokine responses. Front Immunol. 2017;8:29. https://doi.org/10.3389/fimmu.2017.00029.
Piehler J, Thomas C, Garcia KC, Schreiber G. Structural and dynamic determinants of type I interferon receptor assembly and their functional interpretation. Immunol Rev. 2012;250(1):317–34. https://doi.org/10.1111/imr.12001.
Hertzog PJ, Williams BR. Fine tuning type I interferon responses. Cytokine Growth Factor Rev. 2013;24(3):217–25. https://doi.org/10.1016/j.cytogfr.2013.04.002.
Platanias LC. Mechanisms of type-I- and type-II-interferon-mediated signalling. Nat Rev Immunol. 2005;5(5):375–86. https://doi.org/10.1038/nri1604.
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.
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. https://doi.org/10.1016/s0006-8993(99)01328-1.
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;66(10):2058–78. https://doi.org/10.1002/glia.23460.
Engholm-Keller K, Birck P, Storling J, Pociot F, Mandrup-Poulsen T, Larsen MR. TiSH–a robust and sensitive global phosphoproteomics strategy employing a combination of TiO2, SIMAC, and HILIC. J Proteomics. 2012;75(18):5749–61. https://doi.org/10.1016/j.jprot.2012.08.007.
Mohamad M, Mitchell SJ, Wu LE, White MY, Cordwell SJ, Mach J, et al. Ultrastructure of the liver microcirculation influences hepatic and systemic insulin activity and provides a mechanism for age-related insulin resistance. Aging Cell. 2016;15(4):706–15. https://doi.org/10.1111/acel.12481.
Yang IV, Chen E, Hasseman JP, Liang W, Frank BC, Wang S, et al. Within the fold: assessing differential expression measures and reproducibility in microarray assays. Genome Biol. 2002;3(11):0062. https://doi.org/10.1186/gb-2002-3-11-research0062.
McIlwain S, Mathews M, Bereman MS, Rubel EW, MacCoss MJ, Noble WS. Estimating relative abundances of proteins from shotgun proteomics data. BMC Bioinform. 2012;13:308. https://doi.org/10.1186/1471-2105-13-308.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57. https://doi.org/10.1038/nprot.2008.211.
Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13. https://doi.org/10.1093/nar/gkn923.
Kramer A, Green J, Pollard J Jr, Tugendreich S. Causal analysis approaches in ingenuity pathway analysis. Bioinform. 2014;30(4):523–30. https://doi.org/10.1093/bioinformatics/btt703.
Chou MF, Schwartz D. Biological sequence motif discovery using motif-x. Curr Protoc Bioinformatics. 2011;13:5–24. https://doi.org/10.1002/0471250953.bi1315s35.
Schwartz D, Gygi SP. An iterative statistical approach to the identification of protein phosphorylation motifs from large-scale data sets. Nat Biotechnol. 2005;23(11):1391–8. https://doi.org/10.1038/nbt1146.
Song C, Ye M, Liu Z, Cheng H, Jiang X, Han G, et al. Systematic analysis of protein phosphorylation networks from phosphoproteomic data. Mol Cell Proteom. 2012;11(10):1070–83. https://doi.org/10.1074/mcp.M111.012625.
Hornbeck PV, Zhang B, Murray B, Kornhauser JM, Latham V, Skrzypek E. PhosphoSitePlus, 2014: mutations, PTMs and recalibrations. Nucleic Acids Res. 2015;43(Database issue):D512-20. https://doi.org/10.1093/nar/gku1267.
Caenepeel S, Charydczak G, Sudarsanam S, Hunter T, Manning G. The mouse kinome: discovery and comparative genomics of all mouse protein kinases. Proc Natl Acad Sci USA. 2004;101(32):11707–12. https://doi.org/10.1073/pnas.0306880101.
R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2017.
Warnes GR, Bolker B, Bonebakker L, Gentleman R, Liaw WHA, Lumley T, et al. gplots: Various R Programming Tools for Plotting Data. 3.0.1 ed2016.
Wickham H. ggplot2: Elegant Graphics for Data Analysis: Springer Publishing Company, Incorporated; 2009. 216 p.
Chen H. VennDiagram: Generate High-Resolution Venn and Euler Plots. 1.6.20 ed2018.
Supek F, Bosnjak M, Skunca N, Smuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS ONE. 2011;6(7): e21800. https://doi.org/10.1371/journal.pone.0021800.
Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7:539. https://doi.org/10.1038/msb.2011.75.
Rusinova I, Forster S, Yu S, Kannan A, Masse M, Cumming H, et al. Interferome v2.0: an updated database of annotated interferon-regulated genes. Nucleic Acids Res. 2013;41(Database issue): D1040–6. Doi: https://doi.org/10.1093/nar/gks1215
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.
Do Carmo S, Crynen G, Paradis T, Reed J, Iulita MF, Ducatenzeiler A, et al. Hippocampal proteomic analysis reveals distinct pathway deregulation profiles at early and late stages in a rat model of Alzheimer’s-Like Amyloid Pathology. Mol Neurobiol. 2018;55(4):3451–76. https://doi.org/10.1007/s12035-017-0580-9.
Kovarik P, Mangold M, Ramsauer K, Heidari H, Steinborn R, Zotter A, et al. Specificity of signaling by STAT1 depends on SH2 and C-terminal domains that regulate Ser727 phosphorylation, differentially affecting specific target gene expression. EMBO J. 2001;20(1–2):91–100. https://doi.org/10.1093/emboj/20.1.91.
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.
Roy J, Cyert MS. Cracking the phosphatase code: docking interactions determine substrate specificity. Sci Signal. 2009;2(100):re9. https://doi.org/10.1126/scisignal.2100re9.
Chang EJ, Begum R, Chait BT, Gaasterland T. Prediction of cyclin-dependent kinase phosphorylation substrates. PLoS ONE. 2007;2(7): e656. https://doi.org/10.1371/journal.pone.0000656.
Huttlin EL, Jedrychowski MP, Elias JE, Goswami T, Rad R, Beausoleil SA, et al. A tissue-specific atlas of mouse protein phosphorylation and expression. Cell. 2010;143(7):1174–89. https://doi.org/10.1016/j.cell.2010.12.001.
Wang J, Campbell IL. Innate STAT1-dependent genomic response of neurons to the antiviral cytokine alpha interferon. J Virol. 2005;79(13):8295–302. https://doi.org/10.1128/jvi.79.13.8295-8302.2005.
Humphrey SJ, Azimifar SB, Mann M. High-throughput phosphoproteomics reveals in vivo insulin signaling dynamics. Nat Biotechnol. 2015;33(9):990–5. https://doi.org/10.1038/nbt.3327.
Humphrey SJ, Yang G, Yang P, Fazakerley DJ, Stockli J, Yang JY, et al. Dynamic adipocyte phosphoproteome reveals that Akt directly regulates mTORC2. Cell Metab. 2013;17(6):1009–20. https://doi.org/10.1016/j.cmet.2013.04.010.
Olsen JV, Blagoev B, Gnad F, Macek B, Kumar C, Mortensen P, et al. Global, in vivo, and site-specific phosphorylation dynamics in signaling networks. Cell. 2006;127(3):635–48. https://doi.org/10.1016/j.cell.2006.09.026.
Osinalde N, Moss H, Arrizabalaga O, Omaetxebarria MJ, Blagoev B, Zubiaga AM, et al. Interleukin-2 signaling pathway analysis by quantitative phosphoproteomics. J Proteomics. 2011;75(1):177–91. https://doi.org/10.1016/j.jprot.2011.06.007.
Pinto SM, Nirujogi RS, Rojas PL, Patil AH, Manda SS, Subbannayya Y, et al. Quantitative phosphoproteomic analysis of IL-33-mediated signaling. Proteomics. 2015;15(2–3):532–44. https://doi.org/10.1002/pmic.201400303.
Ping L, Duong DM, Yin L, Gearing M, Lah JJ, Levey AI, et al. Global quantitative analysis of the human brain proteome in Alzheimer’s and Parkinson’s Disease. Scientific Data. 2018;5: 180036. https://doi.org/10.1038/sdata.2018.36.
Parker BL, Yang G, Humphrey SJ, Chaudhuri R, Ma X, Peterman S, et al. Targeted phosphoproteomics of insulin signaling using data-independent acquisition mass spectrometry. Sci Signal. 2015;8(380):rs6. https://doi.org/10.1126/scisignal.aaa3139.
Marcantonio M, Trost M, Courcelles M, Desjardins M, Thibault P. Combined enzymatic and data mining approaches for comprehensive phosphoproteome analyses: application to cell signaling events of interferon-gamma-stimulated macrophages. Mol Cell Proteom. 2008;7(4):645–60. https://doi.org/10.1074/mcp.M700383-MCP200.
Detjen KM, Welzel M, Farwig K, Brembeck FH, Kaiser A, Riecken EO, et al. Molecular mechanism of interferon alfa-mediated growth inhibition in human neuroendocrine tumor cells. Gastroenterology. 2000;118(4):735–48.
Inamura K, Matsuzaki Y, Uematsu N, Honda A, Tanaka N, Uchida K. Rapid inhibition of MAPK signaling and anti-proliferation effect via JAK/STAT signaling by interferon-alpha in hepatocellular carcinoma cell lines. Biochem Biophys Acta. 2005;1745(3):401–10. https://doi.org/10.1016/j.bbamcr.2005.06.003.
Zheng H, Hu P, Quinn DF, Wang YK. Phosphotyrosine proteomic study of interferon alpha signaling pathway using a combination of immunoprecipitation and immobilized metal affinity chromatography. Mol Cell Proteom. 2005;4(6):721–30. https://doi.org/10.1074/mcp.M400077-MCP200.
Chen X, Vinkemeier U, Zhao Y, Jeruzalmi D, Darnell JE Jr, Kuriyan J. Crystal structure of a tyrosine phosphorylated STAT-1 dimer bound to DNA. Cell. 1998;93(5):827–39.
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. https://doi.org/10.1002/acn3.213.
Ramesh V, Bernardi B, Stafa A, Garone C, Franzoni E, Abinun M, et al. Intracerebral large artery disease in Aicardi-Goutieres syndrome implicates SAMHD1 in vascular homeostasis. Dev Med Child Neurol. 2010;52(8):725–32. https://doi.org/10.1111/j.1469-8749.2010.03727.x.
Uggetti C, La Piana R, Orcesi S, Egitto MG, Crow YJ, Fazzi E. Aicardi-Goutieres syndrome: neuroradiologic findings and follow-up. AJNR Am J Neuroradiol. 2009;30(10):1971–6. https://doi.org/10.3174/ajnr.A1694.
Marttinen M, Paananen J, Neme A, Mitra V, Takalo M, Natunen T, et al. A multiomic approach to characterize the temporal sequence in Alzheimer’s disease-related pathology. Neurobiol Dis. 2019;124:454–68. https://doi.org/10.1016/j.nbd.2018.12.009.
Kempf SJ, Metaxas A, Ibanez-Vea M, Darvesh S, Finsen B, Larsen MR. An integrated proteomics approach shows synaptic plasticity changes in an APP/PS1 Alzheimer’s mouse model. Oncotarget. 2016;7(23):33627–48. https://doi.org/10.18632/oncotarget.9092.
Zhang X, Ma H, Huang J, Dai Y. Characterization of the Phosphoproteome in SLE Patients. PLoS ONE. 2012;7(12): e53129. https://doi.org/10.1371/journal.pone.0053129.
Braithwaite SP, Stock JB, Lombroso PJ, Nairn AC. Protein phosphatases and Alzheimer’s disease. Prog Mol Biol Transl Sci. 2012;106:343–79. https://doi.org/10.1016/b978-0-12-396456-4.00012-2.
Nematullah M, Hoda MN, Khan F. Protein phosphatase 2A: a double-faced phosphatase of cellular system and its role in neurodegenerative disorders. Mol Neurobiol. 2018;55(2):1750–61. https://doi.org/10.1007/s12035-017-0444-3.
Saunders A, Macosko EZ, Wysoker A, Goldman M, Krienen FM, de Rivera H, et al. Molecular diversity and specializations among the cells of the adult mouse brain. Cell. 2018;174(4):1015-30.e16. https://doi.org/10.1016/j.cell.2018.07.028.
Bakir-Gungor B, Sezerman OU. The identification of pathway markers in intracranial aneurysm using genome-wide association data from two different populations. PLoS ONE. 2013;8(3): e57022. https://doi.org/10.1371/journal.pone.0057022.
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. https://doi.org/10.1182/blood-2016-05-715987.
Sherer Y, Gorstein A, Fritzler MJ, Shoenfeld Y. Autoantibody explosion in systemic lupus erythematosus: more than 100 different antibodies found in SLE patients. Semin Arthritis Rheum. 2004;34(2):501–37.
Cuadrado E, Vanderver A, Brown KJ, Sandza A, Takanohashi A, Jansen MH, et al. Aicardi-Goutieres syndrome harbours abundant systemic and brain-reactive autoantibodies. Ann Rheum Dis. 2015;74(10):1931–9. https://doi.org/10.1136/annrheumdis-2014-205396.
Wang L, Dai Y, Qi S, Sun B, Wen J, Zhang L, et al. Comparative proteome analysis of peripheral blood mononuclear cells in systemic lupus erythematosus with iTRAQ quantitative proteomics. Rheumatol Int. 2012;32(3):585–93. https://doi.org/10.1007/s00296-010-1625-9.
Shiozawa S, Kuroki Y, Kim M, Hirohata S, Ogino T. Interferon-alpha in lupus psychosis. Arthritis Rheum. 1992;35(4):417–22.
Daniels BP, Jujjavarapu H, Durrant DM, Williams JL, Green RR, White JP, et al. Regional astrocyte IFN signaling restricts pathogenesis during neurotropic viral infection. J Clin Investig. 2017;127(3):843–56. https://doi.org/10.1172/jci88720.
Oh C, Ryoo J, Park K, Kim B, Daly MB, Cho D, et al. A central role for PI3K-AKT signaling pathway in linking SAMHD1-deficiency to the type I interferon signature. Sci Rep. 2018;8(1):84. https://doi.org/10.1038/s41598-017-18308-8.
Pillai VB, Sundaresan NR, Samant SA, Wolfgeher D, Trivedi CM, Gupta MP. Acetylation of a conserved lysine residue in the ATP binding pocket of p38 augments its kinase activity during hypertrophy of cardiomyocytes. Mol Cell Biol. 2011;31(11):2349–63. https://doi.org/10.1128/mcb.01205-10.
Decker T, Kovarik P. Serine phosphorylation of STATs. Oncogene. 2000;19(21):2628–37. https://doi.org/10.1038/sj.onc.1203481.
Goh KC, Haque SJ, Williams BR. p38 MAP kinase is required for STAT1 serine phosphorylation and transcriptional activation induced by interferons. EMBO J. 1999;18(20):5601–8. https://doi.org/10.1093/emboj/18.20.5601.
Kovarik P, Stoiber D, Eyers PA, Menghini R, Neininger A, Gaestel M, et al. Stress-induced phosphorylation of STAT1 at Ser727 requires p38 mitogen-activated protein kinase whereas IFN-gamma uses a different signaling pathway. Proc Natl Acad Sci USA. 1999;96(24):13956–61. https://doi.org/10.1073/pnas.96.24.13956.
de Weerd NA, Vivian JP, Nguyen TK, Mangan NE, Gould JA, Braniff SJ, et al. Structural basis of a unique interferon-β signaling axis mediated via the receptor IFNAR1. Nat Immunol. 2013;14(9):901–7. https://doi.org/10.1038/ni.2667.
Moraga I, Harari D, Schreiber G, Uzé G, Pellegrini S. Receptor density is key to the alpha2/beta interferon differential activities. Mol Cell Biol. 2009;29(17):4778–87. https://doi.org/10.1128/mcb.01808-08.
Fox LE, Locke MC, Lenschow DJ. Context is key: delineating the unique functions of IFNα and IFNβ in disease. Front Immunol. 2020;11: 606874. https://doi.org/10.3389/fimmu.2020.606874.
Rengachari S, Groiss S, Devos JM, Caron E, Grandvaux N, Panne D. Structural basis of STAT2 recognition by IRF9 reveals molecular insights into ISGF3 function. Proc Natl Acad Sci USA. 2018;115(4):E601–9. https://doi.org/10.1073/pnas.1718426115.
Ferguson FM, Gray NS. Kinase inhibitors: the road ahead. Nat Rev Drug Discov. 2018;17(5):353–77. https://doi.org/10.1038/nrd.2018.21.
Bhullar KS, Lagaron NO, McGowan EM, Parmar I, Jha A, Hubbard BP, et al. Kinase-targeted cancer therapies: progress, challenges and future directions. Mol Cancer. 2018;17(1):48. https://doi.org/10.1186/s12943-018-0804-2.
Sanchez GAM, Reinhardt A, Ramsey S, Wittkowski H, Hashkes PJ, Berkun Y, et al. JAK1/2 inhibition with baricitinib in the treatment of autoinflammatory interferonopathies. J Clin Investig. 2018;128(7):3041–52. https://doi.org/10.1172/jci98814.
Kontzias A, Kotlyar A, Laurence A, Changelian P, O’Shea JJ. Jakinibs: a new class of kinase inhibitors in cancer and autoimmune disease. Curr Opin Pharmacol. 2012;12(4):464–70. https://doi.org/10.1016/j.coph.2012.06.008.
Perez-Riverol Y, Csordas A, Bai J, Bernal-Llinares M, Hewapathirana S, Kundu DJ, et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. 2019;47(D1):D442. https://doi.org/10.1093/nar/gky1106.
We thank Dr. Claire L. Thompson for help with the manuscript.
B.V. is supported by the Australian Government Research Training Program scholarship. M.Y.W. is supported by Australian Research Council DECRA award and Heart Foundation (Australia) Future Leader Fellowship.
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 of practice for the care and use of animals for scientific purposes’. Ethics approval for all animal experiments was obtained from the animal ethics committee of The University of Sydney (2014/699).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Viengkhou, B., White, M.Y., Cordwell, S.J. et al. A novel phosphoproteomic landscape evoked in response to type I interferon in the brain and in glial cells. J Neuroinflammation 18, 237 (2021). https://doi.org/10.1186/s12974-021-02277-x