{'fraction_non_alphanumeric': 0.03034142745242593, 'fraction_numerical': 0.018684052895913007, 'mean_word_length': 5.430666666666666, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 6, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
A time-resolved molecular map of the macrophage response to VSV infection

Studying the relationship between virus infection and cellular response is paradigmatic for our understanding of how perturbation changes biological systems. Immune response, in this context is a complex yet evolutionarily adapted and robust cellular change, and is experimentally amenable to molecular analysis. To visualize the full cellular response to virus infection, we performed temporal transcriptomics, proteomics, and phosphoproteomics analysis of vesicular stomatitis virus (VSV)-infected mouse macrophages. This enabled the understanding of how infection-induced changes in host gene and protein expression are coordinated with post-translational modifications by cells in time to best measure and control the infection process. The vast and complex molecular changes measured could be decomposed in a limited number of clusters within each category (transcripts, proteins, and protein phosphorylation) each with own kinetic parameter and characteristic pathways/processes, suggesting multiple regulatory options in the overall sensing and homeostatic program. Altogether, the data underscored a prevalent executive function to phosphorylation. Resolution of the molecular events affecting the RIG-I pathway, central to viral recognition, reveals that phosphorylation of the key innate immunity adaptor mitochondrial antiviral-signaling protein (MAVS) on S328/S330 is necessary for activation of type-I interferon and nuclear factor κ B (NFκB) pathways. To further understand the hierarchical relationships, we analyzed kinase–substrate relationships and found RAF1 and, to a lesser extent, ARAF to be inhibiting VSV replication and necessary for NFκB activation, and AKT2, but not AKT1, to be supporting VSV replication. Integrated analysis using the omics data revealed co-regulation of transmembrane transporters including SLC7A11, which was subsequently validated as a host factor in the VSV replication. The data sets are predicted to greatly empower future studies on the functional organization of the response of macrophages to viral challenges.


INTRODUCTION
To maintain host homeostasis and provide protection against invading pathogens, cells are constantly sensing the microenvironment through signaling machinery that, when activated, drives downstream inflammatory responses. Innate immune responses are initiated by germline encoded pattern recognition receptors such as Toll-like receptors and retinoic acid inducible gene (RIG)-like receptors expressed on the cell surface and throughout the cell cytoplasm, as well as integral membrane compartments. 1 These receptors sense pathogen-associated molecular patterns to induce strong pathogen-specific innate immune responses. Antiviral responses activated downstream from virus-specific pathogen-associated molecular patterns include the activation of type-I interferons, transcriptional activation of subsets of inflammatory genes, and a cytokine/chemokine milieu that in turn contributes to the activation of the adaptive immune response. 2 In order to activate a potent inflammatory immune response, cells employ different layers of regulation that range from chromatin remodeling and transcriptional or translational modulation to protein folding, post-translational modifications (PTMs), and localization in the appropriate subcellular compartment. 3 It is plausible that these different regulatory levels enable self-governing and failsafe mechanisms, ensuring clearance of the pathogen while minimizing damage to the host.
Several systematic studies using transcriptomics and functional RNA interference (RNAi) screens have begun to identify (i) the genome-scale regulatory networks that underlie a pathogen response; 4 (ii) a specific FOXO3-IRF7 feed-forward loop in antiviral response; 5 (iii) development of a virus-specific regulatory map; 6 (iv) identification of effector molecules of type-I interferon response 7 and several key regulators of innate immunity such as cyclic GMP-AMP synthase (cGAS) 8 and interferon induced protein with tetratricopeptide repeats 1 (IFIT1). 9 Further, mass spectrometry (MS)-based quantitative proteomics has proven to be a powerful complementary method to study innate immune responses. [10][11][12] These systems-level approaches have begun to (e) Significantly enriched biological processes (P value ⩽ 0.005) at transcriptome, proteome, and phosphoproteome. (f) The normalized mean expression changes for each time point after VSV infection with error bars depicting the extent of expression changes for the subset of genes associated with three representative biological processes that are significantly regulated at transcriptome/proteome/phosphoproteome. Heatmap depicts the expression changes of selected genes associated with the respective biological process. reveal how viruses have evolved ways to tame the host immune system by modulating the key effector proteins and signaling pathways as well disrupting homeostatic sensing mechanisms. [13][14][15] Together, these studies shows us that systems-level methodologies can assist in understanding the full spectrum of activation and regulation behind this complex machinery, and aid in defining the global and unbiased dynamics of the host response. 16 In this study, we aimed to map the complex cellular immune response of murine macrophages to vesicular stomatitis virus (VSV) infection in an unbiased approach, profiling the transcriptome, proteome, and phosphoproteome changes in a quantitative and time-resolved manner (Figure 1a). This allowed the global evaluation of the dynamics of host genes that have a role in infection, commonly referred to as the 'Infectome', and to assess the extent and the logic of the cellular response. We chose VSV, a negative-sense single-stranded RNA virus, for its wide model usage, safety, convenience, and parsimonious molecular functions, encoding five proteins. An integrated bioinformatics analysis of these data sets allowed insights into infection-induced gene expression dynamics, identification of many known and novel differentially regulated genes and proteins during the infection process, as well as numerous altered known and novel phosphosites enabling to explore potential kinase-substrate associations and perturbed signaling pathways. Overall, our data revealed that (i) molecular responses can be clustered by their kinetics and nature in modules with different properties and distinctive features-interestingly, several genes associated with these clusters had no previous innate immune annotation; (ii) phosphorylation appears to occupy a hierarchically dominant challenge-specific function; and (iii) phosphorylation of the innate immunity adapter MAVS has a key role in controlling innate immunity response pathways.

OMICS profiling reveals global insights into gene expression dynamics upon VSV infection
The cellular response of a murine macrophage cell line (RAW264.7) infected with a mutant strain of VSV (VSV-M2) was profiled over the period of time required for infection. Three time points were selected including 20 min, 3, and 6 h. VSV-M2 is not competent for suppression of cytoplasmic recognition 17 and consequently triggers a robust interferon (IFN) response (Supplementary Figure S1A,B). The upregulation of type-I IFNs, IFITs, interleukins (ILs), and chemokines at the transcript level measured using RNA sequencing (RNA-Seq) confirmed that the cells responded to the viral infection (Figure 1b). In total, we were able to quantify 19,889 transcripts, 5,828 proteins, and 1,799 phosphoproteins (available as a resource at www.infectome-map.org). The phosphoproteins included 6,526 phosphoserine, 2,188 phosphothreonine, and 817 phosphotyrosine sites. Using an absolute robust Z score threshold of 2.5, we identified 2,002 differentially regulated genes (Supplementary Figure S1C). These included 1,432 transcripts, 652 proteins, and 401 phosphoproteins (Supplementary Table S1). The differentially regulated phosphoresidues were comprised 506 phosphoserine; 102 phosphothreonine; and 46 phosphotyrosine residues. We further used a rank-based normalization method to calculate a normalized fold change that allowed the comparison of changes in the transcriptome, proteome, and phosphoproteome (Supplementary Figure S1D).
To understand the expression dynamics of transcripts, proteins, and phosphoproteins, we estimated the skewness in the distribution of expression changes upon VSV infection. We found that there is a statistically significant level of skewness in the individual distributions revealing distinct dynamics in innate immune response upon VSV infection (Figure 1c). Transcriptome distribution revealed a large downregulation at 20 min and 3 h, whereas at 6 h a profound upregulation was observed (Figure 1c). Proteome distribution showed an upregulation at 3 and 6 h, and a subtle downregulation at 20 min. Strikingly, phosphoproteome showed a strong upregulation at all three time points (Figure 1c). The overlap between the genes that are differentially expressed at transcript/protein/phosphoprotein levels, taking into account all time points, was assessed. Interestingly, we observed that only a small number of genes changed at two or three molecular levels, thus revealing a discrete molecular-level footprint (Figure 1d). We also observed that differentially regulated proteins and phosphoproteins are significantly enriched for housekeeping genes (Supplementary Figure S1E). 18 This is in line with the published observation that lipopolysaccharide (LPS)-induced changes in dendritic cells cause remodeling of housekeeping genes, primarily by protein translation or degradation. 19 We performed enrichment analysis for biological processes using genes that were regulated at the transcriptome, proteome, and phosphoproteome levels, and identified 43 significantly enriched processes (P value ⩽ 0.005; Figure 1e). We identified discrete cellular processes that were regulated at one or more molecular levels. For instance, we found genes associated with defense response regulated at the transcriptome level; lipid biosynthetic process at the proteome level; and mRNA processing at the phosphoproteome level ( Figure 1f). We also performed k-means clustering of the transcript, protein, and phosphosite measurements, and were able to segregate the regulated genes into early, late, and mid-late responders to viral infection (Supplementary Figure S1F). The mid-to-late responders at the transcript level-cluster 1 (T-C1) included expected early antiviral genes such as type-I IFN and various ILs. In addition, T-C1 also contains unique genes such as Dusp8, a member of dual specificity phosphatases that regulate MAPK signaling, 20 and MIR702, an anti-apoptotic mirtron that binds to 3'-untranslated repeat of ATF6, which is an endoplasmic reticulum stress-associated transcription factor. 21,22 The late-responder cluster at the protein level, Pr-C2, further included known IFN-induced genes such as IFITs ISG15, MDA5 and RSAD2. 2 In addition, GBP5, which is reported to promote NLRP3 inflammasome assembly upon bacterial infection, 23 and CMPK2, a mitochondrial gene involved in the nucleotide salvage pathway, 24 were also part of Pr-C2 cluster, hinting at possible involvement in antiviral responses. Phosphorylation sites in the early response cluster Ph-C2 included several innate immune-regulated genes with known sites, e.g., T7/S9 on GSK3B and several novel sites, e.g., S208 on OAS1A, S506 on IRAK2, and S3857 on UBR4. Further, we examined the overlap between phosphoproteome clusters to identify proteins with multiple differentially regulated phosphosites with varying kinetics. For instance, we found overlap between the phosphoproteome clusters Ph-C1 and Ph-C8, indicating a strong phosphorylation regulation (Supplementary Figure S1G). These included RB1, a master cell cycle regulator whose phosphorylation state determines its E2F-binding state; 25 SPAG9, a scaffold protein that is required for JNK signaling; 26 and ZC3HC1, a component of the SCF-type E3 ligase system involved in ubiquitination of cell cycle proteins. 27 To address the specificity of VSV-induced transcriptome changes with other antiviral response, we compared our 6 h transcriptome data with the transcriptome of RAW264.7 cells infected with Murine norovirus (MNV1) for 9 h. 28 There was a good correlation (Pearson correlation = 0.53) between the differentially expressed VSV-induced transcriptome at 6 h with MNV1-induced transcriptome at 9 h. The transcript levels of type-I IFNs, IFN-stimulated genes (ISGs), and proinflammatory cytokines were upregulated in both VSV and MNV1 transcriptome (Supplementary Figure S1H). To elucidate the specificity of VSV-induced phosphoproteome change, we compared our data set with LPS-induced phosphoproteome changes in RAW264.7 cells. 11 As these two data sets were generated in different technical settings, we used a robust Z score cutoff of ± 2.5 for the comparison in order to create as little bias as possible. On the basis of this threshold, there were 425 and 469 differentially regulated phosphoproteins that are LPS-and VSV-induced, respectively (Supplementary Figure S1I). There was an overlap of 137 phosphoproteins common to LPS and VSV infection. Interestingly, only 82 phosphosites from these 137 proteins were common to both LPS and VSV (Supplementary Figure S1I). Overall, whereas there is a subset of phosphoproteome that is common to LPS and VSV, and transcriptome that is common to VSV and MNV1, a large portion of the changes are distinct depending on the stimulation, indicating the presence of yet unknown differential regulatory factors and events.
A top-down strategy to identify co-regulated gene modules reveals co-regulation of transmembrane transporters and SLC7A11 as a host factor in the VSV life cycle Gene regulation during cellular perturbation events such as infection is critical for mounting an antimicrobial state while minimizing damage to the host. Systems-level approaches have aided the identification of some of these regulatory circuits  associated with innate immune response. 4,5,29 Using normalized RNA and protein fold changes of the differentially regulated genes (Supplementary Figures S1D and S2A), we constructed a gene-gene correlation matrix that allowed us to identify 12 co-regulated gene modules (Figure 2a and Supplementary Figure S2B). Gene ontology analysis showed several significantly enriched biological processes that were associated with each specific module, and some that were associated with multiple gene modules (Figure 2b). We performed enrichment analysis for transcription factor gene set based on target-binding sites using MSigDB 30 and summarized the results as a network, where nodes comprises co-regulated modules and transcription factors; moreover, edges represent a regulatory relationship (Figure 2c). We found both common and unique transcription factors associated with the co-regulated modules. Several members of the E26 transcription factor family including ETS1, ZEB1, SPI1, TCF3, ETS2, and NFE2L1 were associated with one or more modules.
In addition, such a co-regulation network could also be used for tracking the temporal kinetics of specific family or group of genes triggered by specific perturbations such as ISGs. We investigated the kinetics of ISGs 7 in the co-regulated gene modules and found four distinct groups based on changes in the transcript and protein levels after VSV infection (Supplementary Figure S2C), which could be the potential subsets of co-regulated antiviral effector genes. This network could serve as a point of departure to model host gene regulatory networks after viral infection. For instance, transcription factors ETS2 and FOXO4 that have a subset of common target genes in the co-regulated modules M1 and M7 are likely to share functional roles.
Gene ontology analysis for significantly enriched biological processes showed a group of transmembrane transporters as part of the co-regulated module M1 (Figure 2d). These included a number of genes from the understudied solute carrier family. 31 One of the highly upregulated genes in this group was an aminoacid transporter SLC7A11 that transports cysteine/glutamate. 32 SLC7A11 has also been described as a fusion-entry receptor of Kaposis's sarcoma-associated herpes virus. 33 More recently, SLC7A11 was reported to have a decisive role in reactive oxygen species-induced ferroptosis and p53-mediated tumor suppression. 34 We generated a CRISPR-Cas9-based knockout of SLC7A11 in HAP1 cells 35  Regulated novel kinases and substrates revealed by integrating publicly available kinase-substrate associations To understand the regulated kinome landscape upon VSV infection, we integrated our OMICS data with publicly available kinase-substrate interactions 36 and assembled a 638-interactioncontaining kinase-substrate network (Figure 3a) that included 68 kinases, which were regulated at the transcript/protein/phosphoprotein level and 181 substrates, which were differentially phosphorylated (Supplementary Table S2). The regulated kinases were members of nine different kinase family subgroups (Supplementary Figure S3A and B). Comparing this list with a transcription factor database, 37 we uncovered 13 transcription factors among the differentially phosphorylated substrates including Protein C-ets-2 (ETS2) at S294 and S297. To estimate the number of genes that are linked to innate immunity, we intersected this network with InnateDB. 38 We found 29 kinases and 41 substrates annotated as innate immunity genes, thus highlighting the identified differentially regulated novel kinases and substrates having a role in phosphorylation signaling upon viral infection. Clustering the phosphorylated peptides, including the kinases and substrates in this network, into seven different groups (Figures 3a and b and Supplementary Figure S3C) allowed us to identify the potential key kinase-substrate associations that showed coherent phosphorylation kinetics (Figure 3c). This kinase-substrate network also included several differentially phosphorylated kinases (Supplementary Figure S3D). On the basis of the temporal phosphorylation pattern, JAK2 and WNK1 appear to be late-responding kinases after infection, whereas AKT2, SYK, PAK1, PAK2, ABL1, ABL2, and ARAF respond early. Within this kinase-substrate network, we observed a potential link between the early responding kinases PAK2/ABL1 and the transcription factor JUN (AP-1), which showed coherent differential phosphorylation ( Figure 3d). We further investigated whether the phosphorylation kinetics of JUN correlated with the transcriptome changes of JUN target genes. To this end, we identified three clusters in the transcriptome (T-C1, T-C5, and T-C7; Supplementary Figure S1F) that showed similar temporal kinetics as compared with JUN phosphorylation. We performed motif analysis on these three clusters for enrichment of transcription factor-binding motifs using MSigDB and found statistically significant enrichment for JUN (P valueo 5.43E − 15; Figure 3e). These temporal data along with the motif analysis provide a potential functional link between the kinases PAK2/ABL1 and activation of transcription factor JUN in antiviral response; moreover, the data set provides the feasibility to identify other potential regulatory mechanisms by integrating phosphorylation events of transcription factors with temporal transcriptome data.
We sought to validate the early phosphorylation role revealed above for AKT2, a key factor in PI3K/mTOR signaling, 39 by infecting CRISPR-Cas9-based AKT2 knockout HAP1 cells with VSV-GFP. The AKT2 knockout cells showed a significant reduction in GFP-positive cells after 6 h of infection with MOI10, thereby hinting that AKT2 kinase may have a role in efficient VSV infection or replication (Figure 3f). HAP1 cells treated with the Pan-AKT inhibitors Ipatasertib and MK-2206, but not the AKT1-specific inhibitor, Perifosine, had a reduced infection after 6 h at an MOI of 10 when compared with dimethylsulfoxide-treated cells (Figure 3g and Supplementary Figure S3E). The inhibitor phenotypes support what the omics data integration indicates: AKT2, but not AKT1, was revealed to be differentially regulated by phosphorylation during virus infection.
Among the various differentially phosphorylated kinases upon VSV infection were two members of the RAF family of kinases, ARAF (A-Raf) and RAF1 (c-Raf), which were phosphorylated at S157 and S29, respectively. To investigate the roles of these kinases in regulating virus replication and immune response, we again made CRISPR-Cas9-based ARAF or RAF1 knockout in HAP1. RAF1 and, to a lesser extent, ARAF, has been linked to innate immune activity downstream from virus recognition, 40,41 and, in fact, after 4 h of infection with VSV-GFP at MOI2, the RAF1 KO cells, but not the ARAF KO, had significantly more infection compared with the WT (Figure 3h). A reduced viral MOI and infection time was used here compared with other experiments in order to limit virus turnover duration as the knockout cell lines were more susceptible VSV-mediated lysis. The increased infection phenotype was mirrored by a significant decrease in IL-6 transcription in the RAF1 KO line at the same time point (Supplementary Figure S3F). Further, there was a decrease of P65 phosphorylation in both the RAF1 and the ARAF knockout lines compared with WT at 30 min, whereas only the RAF1 KO had a pronounced and consistent loss of the phosphorylation at a 60 min post infection (Figure 3i), whereas p65 transcription levels are unchanged (Supplementary Figure S3G). Transcriptional activation of IL-6 after phosphorylation of P65 is a known and well-studied pathway post-NFκB activation. [42][43][44] Together, this indicates that RAF1 and, to a less extent, ARAF are potentially activated during early VSV infection, and may mediate a controlling immune response activated by VSV-as elucidated by mining the kinome landscape (Figure 3a).
Integrated analysis to understand the dynamics of perturbed signaling pathways To identify the signaling pathways that are perturbed upon VSV infection, we performed gene-set enrichment for signaling pathways at transcriptome, proteome, and phosphoproteome levels. This enrichment, for instance, revealed that the RIG-I signaling pathway was perturbed at all three molecular levels: p53 signaling and cytosolic DNA-sensing pathway at transcriptome and proteome levels, and the MAPK signaling pathway at transcriptome and phosphoproteome levels. In addition, we found that pathways such as JAK-STAT signaling and NOD-like receptor signaling were enriched at the transcriptome level, whereas focal adhesion and tight junction signaling was enriched at the phosphoproteome level (Supplementary Figure S4A). This substantiated the possibility that the cellular repertoire of possible responses to viral perturbation is varied and that, although some part of the program required transcriptional activation followed by protein changes, others are purely dependent on phosphorylation within this time frame. To obtain a more quantitative measure of the perturbation, we carried out Signaling Pathway Impact Analysis 45 using normalized fold changes of transcriptome, proteome, and phosphoproteome data (Supplementary Figure S1D). Signaling Pathway Impact Analysis analysis on 148 KEGG signaling pathways, taking into account the directionality and type of the interactions, showed enrichment for several pathways as seen in the classical enrichment method. In addition, we were able to derive the perturbation state of these pathways at each molecular level over time (Figure 4a and Supplementary  Table S3). For instance, tracking RIG-I signaling showed a time-dependent activation, whereas signaling during Influenza A virus infection has a sustained activation soon after infection.
The RIG-I signaling pathway is important for sensing several viruses including VSV, Influenza virus, Sendai virus, and Hepatitis C virus, 46 and initiating strong IFN response. The data resolution from the OMICS profiling was able to capture more than 95% of the genes at one or more molecular level. Thus, we created a temporal map describing the transcriptome, proteome, and phosphoproteome of RIG-I signaling upon VSV infection (Figure 4b), including novel phosphorylation sites on the key molecules involved in this signaling pathway. DHX58 (also called LGP2), a member of the RIG-like receptor family of proteins that has a virus-specific role in mediating antiviral response, 47-50 was lacking upon infection phosphorylation at S477 located in the repressor domain. We also uncovered two phosphorylation sites S645 and S648, also not previously reported, that are upregulated on another member of the RIG-like receptor family, IFIH1/MDA5. MAVS, a key adaptor protein in RIG-I signaling, and downstream IFN activation, was found to be hyperphosphorylated upon infection at S328 in our study. In addition, TBK1, a key regulator of IRF3 signaling, 51,52 was hyperphosphorylated at S716 and IRF3 at S388. Altogether, the signaling pathway impact analysis indicates that individual pathways have strong time-and molecule-level-dependent regulation; moreover, the resolution of the data sets allows temporal tracking of the key signaling pathways.
MAVS S328-S330 is required for IL-6 and IFN response Our phosphoproteomics screen revealed that a number of novel phosphosites of the key innate immune signaling proteins and clustering of these phosphosites allowed us to interrogate the subsets with different kinetics. Among the various phosphosites that followed the same temporal kinetics as IRF3 was S328 on MAVS, which was hyperphosphorylated at 3 and 6 h after infection without considerable change at the RNA or protein level (Figure 5a,b). Multiple sequence alignment of the MAVS proteins from mouse, human, and rat revealed that the amino acid S328 is specific to mouse and is replaced by Isoleucine in human and rat, whereas adjacent S330 was conserved in all these species (Supplementary Figure S5A). The region that we identified is different from a recent study 53 that described the function of a novel phosphorylation site on human MAVS. To test the role of this phosphorylation cluster S328-S330 in antiviral signaling, we reconstituted MAVS − / − immortalized primary mouse macrophages with MAVS WT, MAVS S328A, or MAVS S328A/S330A (Supplementary Figure S5B). First, we infected the reconstituted cells with VSV engineered to express GFP under its major promoter. 54 Flow cytometry quantification of the number of GFP-positive cells revealed the double-mutant S328A/S330A to have a significant increase in viral gene production while S328A had only a marginal increase (Figure 5c).
To further define the roles of the phosphosites in NFκB and IFNβ activation downstream of MAVS activation, we analyzed the mRNA expression levels of IFNβ, IL-6, and IFIT1 upon Poly(I:C) stimulation. mRNA expression of these genes was significantly reduced in both S328A and S328A/S330A mutants compared with cells reconstituted with MAVS WT (Figure 5d). Poly(I:C)-induced IFNβ protein secretion was also significantly reduced ( Figure 5e). Upon VSV infection, the double mutant showed a significant decrease in IL-6 and IFNβ transcript levels, whereas the single mutant showed a significant decrease in IL-6 but not IFNβ transcript levels (Supplementary Figure S5C,D). In addition, both mutants showed a decrease in IRF3 phosphorylation, with the double mutant having a profound effect (Supplementary Figure S5E). We further investigated whether the MAVS phosphorylation had any impact on other activation of other antiviral genes. To this end, we performed RNA-Seq-based global transcriptome analysis upon stimulation with Poly(I:C) for 3 and 9 h. We found that ISGs were part of the differentially regulated transcripts upon Poly(I:C) stimulation. Detailed analysis of this subset showed that the mRNA levels of several ISGs such as IFIT1, IFIT2, IFIT3, MX1, ISG15, CXCL10, and ISG20 were significantly reduced in both the S328A and S328A/S330A mutants (Figure 5f). This shows that the MAVS phosphorylation is critical for triggering the transcriptional rewiring during innate immune response. The phospho-mutant MAVS-expressing cells had no significant difference in induction of IL-6 mRNA expression upon LPS stimulation (Supplementary Figure S5F). Confocal microscopy revealed that the MAVS constructs were properly localized to mitochondria in our reconstituted cells (Supplementary Figure S5G). The identification of this specific MAVS phosphorylation and subsequent validation to link its functional effect emphasizes the importance of using orthogonal technologies to profile the perturbation state of a cell to understand cellular response, in addition to validating the quality of the data made available through this study.

DISCUSSION
How much of the response of mammalian cells to viral challenge is dedicated to recognition of the intrusion, measurement of the severity of the challenge, and mounting an adequate response? Are there processes that need to occur before others? Are there 'housekeeping' processes that would be triggered also by non-microbial perturbations or are most mechanisms evolved for safe-guarding homeostasis to withstand microbes? We joined the effort of several laboratories focused on gaining fundamental insights to answer questions such as these by attempting a global, unbiased survey of the events in macrophages triggered by VSV infection. [55][56][57] The analysis and validation that ensued aimed at deconvoluting the response to discern distinct functional relationships and a possible unifying logic. This study describes discrete patterns of events and reports on observations that allow to formulate distinct, testable, hypotheses. Using an integrated transcriptome, proteome, and phosphoproteome profiling, we have created a quantitative infectome map describing the gene expression dynamics of VSV-infected macrophages. This has led to identification of several known and novel differentially regulated genes that have a potential role in antiviral innate immune response. These data indicate that the temporal dynamics of the host response is distinct, depending on the molecular level, likely reflecting protein translation/degradation rates but also suggestive of specialized regulatory options, which may each be best achieved at the level of RNA, protein, or PTM. 19,[58][59][60][61] Enrichment for so-called housekeeping genes at proteome and phosphoproteome levels may be indicative of the careful orchestration of host's cellular response to perturbation scenarios such as viral infection, and a recent study 19 validated this notion that housekeeping genes are regulated at the proteome level in a LPS-induced dendritic cell response. This study extends this observation to a viral response and also shows a similar response for PTMs such as phosphorylation. The observation that a large number of genes are exclusively regulated by phosphorylation of their products underscores the substantial role of PTM in innate immune response 62 likely evolved because of the requirement of a fast response to pathogens. Analysis of other PTMs amenable to high-throughput identifications, such as acetylation and glycosylation, 63 should shed light on their role in innate immune response. The phenotypic effect of MAVS phosphorylation, which we unraveled in this study, emphasizes on the advantages on using orthogonal measurements to study perturbations.
A data set as large as the one created here requires efficient bioinformatics methods for integration and to allow to detect general principles or formulate precise hypotheses. There is a great variety of possibilities to analyze the data, and in future it should be possible to unearth may other properties that are not touched upon here. The analysis methods that we used in this study, such as clustering of genes based on temporal kinetics, quantitative pathway impact analysis, kinase-substrate interaction networks, and co-regulation networks, provide different perspective depending on the biological question and plane of intellectual focus that is being sought. Whereas the analysis presented identifies biological objects that have expression behaviors that are similar to or undistinguishable from elements known to have important roles in innate immune signaling; the data offer the basis for analyses at the systems-level perspective of regulation of antiviral signaling and response.
A general question that arises during genome-wide expression profiling of host responses at the cellular level is how specific are these changes with respect to the type of perturbation. It is known that there is a transcriptional signature for various innate immune stimuli, 4 but so far no recognizable patterns have been proposed for changes at the PTM level such as phosphorylation. We tried to address this by comparing VSV-induced phosphoproteome changes with LPS-induced changes. What emerges from this comparison is that, whereas there is small part that is common to both, a larger portion of the changes is very stimulus-specific, allowing the speculation that the various biochemical effects offered by phosphorylation (modulation of protein interfaces of the molecular machinery, intrasteric and allosteric regulation of enzymes) may provide the level of fine-tuning necessary to focus on specific challenges and responses. 58,64 In comparison, at the protein level, one can detect a large overlap of commonly regulated proteins but a plethora of phosphosites that are ligandspecific. This shows that phosphoproteome changes will have common and ligand-specific signatures. It will certainly be interesting, in the future, to systematically assess the PTM dynamics of various innate immune stimuli.
Regulation of gene expression occurs at multiple levels such as chromatin remodeling, 65 transcription, 66 post-transcriptional modifications, 67 translation, 68 and PTMs. 69 This study validates this notion that cells have multiple regulatory options under perturbation conditions and is important to profile multiple molecular layers to appreciate the cellular response in full. Without doubt, adding measurements of changes in the metabolism of infected cells should further allow to deconvolute the response. 70,71 There are a number of genes with a potential role in antiviral innate immune response identified in this study, and they could serve as a reference point of future studies, as VSV is a commonly used model virus. Two independent studies have showed that phosphorylation status of ASC at Y146 (ref. 72) and NLRC4 at S533 (ref. 73) is essential for its inflammasome activation. 74 Along these lines, the validation of the phosphosites on MAVS in this study is a clear demonstration of the significance of phosphorylation in regulating innate immune responses. In general, we envisage that the data set provided by this study together with its analytical strategy will prove to be useful to the community of scientists dedicated to the elucidation of the cellular processes involved in recognizing and responding to viral infections. The analysis strategies used in our study could be applied to characterize any immunological perturbation scenarios to get a systems-level overview. Overall, these unique data sets from this study could serve as a useful resource for generating interesting hypotheses.

Phosphopeptide enrichment and offline fractionation
Uninfected and VSV-infected RAW264.7 cells at~70% confluency were washed in cold phosphate-buffered saline (Sigma), cells were scraped, centrifuged, and the cell pellets were snap-frozen in liquid nitrogen. Frozen cell pellets were lysed in 8 M urea lysis buffer containing phosphatase and kinase inhibitors (Roche Diagnostics GmbH, Mannheim, Germany). Total protein concentration was normalized across the samples and a total of 7.7 mg protein per condition was used for tryptic digestion and phosphotyrosine enrichment. Peptide digest (100 μg) was used for whole-cell proteome and global phosphopeptide enrichment. iTRAQ 4-plex labeling was performed according to the instructions of the manufacturer. The whole-cell proteome was separated into 50 fractions and the global phosphopeptide-enriched samples were separated into 10 fractions.
Online LC-MS and data processing MS was performed on a hybrid LTQ-Orbitrap Velos mass spectrometer (ThermoFisher Scientific, Waltham, MA, USA) using the Xcalibur version 2.1.0.1140 coupled to an Agilent 1200 HPLC nanoflow system via a nanoelectrospray ion source using liquid junction (Proxeon, Odense, Denmark). The analyses of the phosphopeptides were performed in a data-dependent acquisition mode using a top 10 hybrid method with a multistage-activation collision-induced dissociation MS 2 event, followed by a higher-energy dissociation MS 2 event for each selected precursor ion. The analysis of the normal proteome was performed with a top 10 higherenergy dissociation method.
Peak list data were extracted using the Proteowizard software (http:// proteowizard.sourceforge.net/-release 3.0.3201) and searched against the SwissProt mouse database (release 2012.05) appended with common contaminants, using Mascot (version 2.3.02; MatrixScience, London, UK) and Phenyx (version 2.5.14; GeneBio, Geneva, Switzerland). Fixed modifications included cysteine carbamidomethylation and iTRAQ labels on lysine and N terminus. Variable modifications included methionine oxidation and phosphorylation of serine, threonine, and tyrosine. One tryptic missed cleavage was allowed. A two-pass search strategy was used to apply a final false discovery rate of 1% at the protein level and 0.1% at the peptide level as detailed in Bennett et al. 75 Protein-level quantitation for total cellular proteome and peptide-level quantification for phosphoproteome data were performed using isobar R package, 76 and robust Z score were calculated. Protein (for total cellular proteome) or peptide (for phosphoproteome) with an absolute Z score ⩾ 2.5 and P value ⩽ 0.05 were considered differentially expressed.

RNA-Seq
Total RNA was extracted using TRIzol reagent (Life Technologies) following the manufacturer's protocol from 1 × 10 7 cells. Quality of RNA was measured using BioRad's Experion RNA Analysis chip, and samples with a RNA integrity number greater than 7 were used for whole-transcriptome analysis. RNA (2 μg) was used for complementary DNA library synthesis using the TruSeq RNA Sample Prep kit v2 (Illumina, San Diego, CA, USA) according to the manufacturer's recommendations. RNA-Seq libraries were sequenced on an Illumina HiSeq 2000 machine with a 50 bp single-end protocol. The sequence reads were aligned to mouse reference genome build mm9 using TopHat2 (ref. 77) and Bowtie2. 78 RPKM (reads per kilobase million) and FPKM (fragments per kilobase million) values were computed using RSeQC package 79 and Cufflinks, 80 respectively. Fold changes were calculated using DEGseq R package 81 applying the MARS (MA-plot-based method with Random Sampling model) algorithm. Robust Z scores were also computed based on the fold change. Transcripts with an absolute Z score ⩾ 2.5 and false discovery rate ⩽ 0.001 were considered differentially expressed.

Bioinformatics analyses
Gene ontology biological process and KEGG pathway enrichment analysis were carried out using DAVID web service. 82 Signaling pathway impact analysis was analyzed using Signaling Pathway Impact Analysis package implemented in R statistical environment (http://www.r-project. org/). The list of kinases was downloaded from UniProt database (http:// www.uniprot.org/docs/pkinfam). Multiple sequence alignment was carried out using ClustalW2 (http://www.ebi.ac.uk/Tools/msa/clustalw2/). Networks were visualized using Cytoscape (http://www.cytoscape.org/). Pearson's correlation method was used for constructing the gene-gene correlation matrix. Optimal numbers of clusters were calculated using the statistical method clusGap, available under the cluster package in R. Transcription factor target site enrichment was performed using the GSEA software from MSigDB (version 5.0; http://www.broadinstitute.org/gsea/msigdb) and top 10 gene sets with false discovery rate q-value ⩽ 0.05 were considered significant. All statistical analyses, data processing, and visualization were performed in the R statistical environment.

Rank-based normalization
Rank-based normalization of ratio distribution was used to map proteomics and phosphoproteomics ratios to the range of observed mRNA ratios and make them comparable. mRNA ratios mr 1 , …, mr N were sorted from largest negative to largest positive. The same sorting operation was carried out with protein abundance ratios pr 1 , …, pr M . Then, the proteins were assigned ratios from the mRNA ratio distribution according to the formula pr i → mr trunc(i/M × N) , trunc(x) being the function returning the largest integer smaller than or equal to x. The same transformation was imposed to phosphoproteomics ratios ppr 1 , …, ppr K : ppr i → mr trunc(i/K × N) . CRISPR-Cas9-mediated generation of HAP1 knockout cells SLC7A11: pX459 plasmid 83 (gift from Feng Zhang lab-Addgene, Cambridge, MA, USA, #48139) was ligated with guide sequence inserts 5ʹ-CACCGCTCCAAAGGAGGTTACCTGC-3ʹ and 5ʹ-AAACGCAGGTAACCTCCTTT GGAGC-3ʹ. HAP1 cells were transiently transfected with for 24 h, followed by selection with Puromycin (5 μg/ml). Single-cell clone was generated by limiting dilution. ARAF and RAF1: lentiCRISPRv2 plasmid 84 (gift from Feng Zhang lab-Addgene #52961) was ligated with 5ʹ-CACCGCCGTGCGTTGC TTGTTGGGC-3ʹ and 5ʹ-AAACGCCCAACAAGCAACGCACGGC-3ʹ in case of ARAF; 5ʹ-CACCGAGCAGTTTGGCTATCAGCGC-3ʹ and 5ʹ-AAACGCGCTGATAGC CAAACTGCTC-3ʹ in case of RAF1. Lentiviruses were produced using secondgeneration packaging plasmids pMD2.G and psPAX2. HEK293T cells were co-transfected with packaging and lentiCRISPRv2 plasmids. Cells were washed 16 h after transfection. Virus-containing supernatants were collected 24 h after washing, filtered, and used for infection of HAP1 cells. Twenty-four hours later, a second round of infection was performed. Forty-eight hours after the last infection, HAP1 cells were selected with Puromycin (5 μg/ml). HAP1 cells deficient in ARAF, RAF1, and AKT2 were also obtained from Haplogen GmbH (Vienna, Austria) produced with CRISPR technology.

Confocal microscopy
Immortalized bone marrow-derived mouse MAVS −/− macrophages reconstituted with hemagglutinin (HA)-tagged MAVS WT or mutants were seeded on glass coverslips overnight. Before fixing, cells were stained with Mitotracker Orange (invitrogen, Carlsbad, CA, USA) according to the manufacturer's protocols. Cells were fixed in 2% formaldehyde with 0.1% triton-X114. Proteins of interest were visualized with the following: anti-HA (Santa Cruz Biotechnology, Y-11)/anti-rabbit AlexaFlour 633 (invitrogen). Slides were visualized using an LSM780 (Carl Zeiss, Oberkochen, Germany) with a GaAsP multidetector unit and two PMTs. Pinhole was set to 1 AU on each channel. Images were taken sequentially. Images were taken at × 100 and analyzed with ImageJ (NIH, open source).