Essential genes of the macrophage response to Staphylococcus aureus exposure

Background Although significant advances have been made in understanding the mechanisms of macrophage response to Staphylococcus aureus infection, the molecular details are still elusive. Identification of the essential genes and biological processes of macrophages that are specifically changed at different durations of S. aureus exposure is of great clinical significance. Methods We aimed to identify the significantly changed genes and biological processes of S. aureus-exposed macrophages. We systematically analyzed the macrophage gene expression profile GSE 13670 database with 8 h, 24 h or 48 h S. aureus infection. The results were further confirmed by western blot and quantitative polymerase chain reaction (qPCR) analyses. Results After 8 h of S. aureus infection, the expression of 624 genes was significantly changed. Six hundred thirteen differentially expressed genes (DEGs) were identified after 24 h of S. aureus infection. Two hundred fifty-three genes were significantly changed after 48 h of S. aureus infection. STAT1 was consistently up-regulated in these three treatments. TP53, JAK2, CEBPA, STAT3, MYC, CTNNB1 and PRKCA were only identified in the 8 h or 24 h S. aureus infection groups. CTNNB1 and PRKCA were for the first time identified as potential essential genes in S. aureus infection of macrophages. In the Gene Ontology (GO) term analysis, the defense response was shown to be the most significantly changed biological process among all processes; KEGG pathway analysis identified the JAK-STAT signaling pathway involved in early infection. Conclusions Our systematic analysis identified unique gene expression profiles and specifically changed biological processes of the macrophage response to different S. aureus exposure times. Electronic supplementary material The online version of this article (10.1186/s11658-018-0090-4) contains supplementary material, which is available to authorized users.

vaccines currently are licensed for S. aureus [3]. Long-term exposure to S. aureus results in the adaptive immune response and influences the immunological memory establishment and maintenance [4,5].
Macrophages play a critical role in innate immune responses to bacterial infection. Macrophages are primary professional phagocytes that are designed to devour and kill microbes by pinocytosis, receptor-mediated endocytosis or phagocytosis. Previous studies have identified profound gene expression profile changes in differently activated macrophages [6]. However, the underlying molecular mechanism of the macrophage response to S. aureus has still not been completely elucidated.
Recently, Koziel et al. indicated that S. aureus might induce cytoprotective mechanisms by regulating the gene expression profiles inside macrophages [7]. However, their study was mainly focused on the mechanism of cell death and apoptosis upon S. aureus infection. To further dissect the global gene expression profile change and corresponding signaling pathways in S. aureus infection-induced immune response, we systematically analyzed the macrophage gene expression profile to identify genes with significant changes in response to S. aureus exposure. We identified several essential genes and pathways during the S. aureus infection in our analysis. The results were further confirmed by western blot and qPCR analyses, providing the first molecular targets for macrophage response to S. aureus infection.

Data collection, data processing and filtering
To retrieve the human monocyte-derived macrophages (hMDMs) gene expression profile, the Gene Expression Omnibus (GEO) microarray expression dataset (http://www.ncbi. nlm.nih.gov/geo/) was searched and GSE 13670 (http://www.ncbi.nlm.nih.gov/geo/query/ acc.cgi?acc=GSE13670) was obtained for our expression analysis. Briefly, control and S. aureus-exposed macrophage cells were analyzed at different time points (8,24 and 48 h after infection) which each includes 5 samples [8]. Unprocessed data (.cel files) were collected. This dataset was Affymetrix Human Genome U133 Plus 2.0. Corresponding probe annotation files were downloaded for further analysis.
The normalization process in GCRMA [9] was used to quantify the microarray signal for our research. The package genefilter [10] in R language was employed to filter out uninformative data, such as control probe sets and other internal controls as well as to remove genes expressed uniformly close to background detection levels. The filter did not remove probe sets without Entrez gene identifiers or with identical Entrez gene identifiers.

Differentially expressed gene analysis
Three statistical comparisons were carried out as follows. Comparison 1 was conducted between the 8 h infection and control groups (hMDM_SA_8hr vs hMDM_control_8hr), Comparison 2 was made between the 24 h infection and control groups (hMDM_SA_ 24hr vs hMDM_SA_24hr), and Comparison 3 was made between 48 h infection and control groups (hMDM_SA_48hr vs hMDM_control_48hr). The package limma [11] in R language was used to identify differential expressed genes from these 3 comparisons. For those probes that had an identical Entrez gene identifier, we only kept the probe showing the greatest variance. Genes with |log2(FC)| > 1.5 and the adjusted p value < 0.01 were accepted as statistically differentially expressed. The p value was adjusted by applying Benjamini and Hochberg's (BH) false discovery rate (FDR) correction to the original p value, and the fold change threshold was set to focus on significantly differentially expressed genes.
The obtained DEGs from three comparisons were compared against each other and the common DEGs in all three comparisons or only in one comparison were identified.

Hierarchical clustering, GO and KEGG pathway analysis
To better understand the global gene expression patterns, we performed hierarchical clustering [12] to classify analyzed samples based on gene expression profiles. The DEGs were further analyzed using Gene Ontology terms (biological processes) and KEGG pathways. Heatmaps for the DEGs classified in targeted biological processes or KEGG pathways were generated by gplots [13] in the R package.
ClusterProfiler [14] in R packages was used to detect Gene Ontology categories and KEGG pathways with significant overrepresentation in DEGs compared with the whole genome. The significantly enriched biological processes were identified with an adjusted p value less than 0.01. For KEGG pathway analysis, the adjusted p value was set to less than 0.05.

Construction of biological network
Protein-protein interaction (PPI) databases were retrieved from HPRD [15], BIOGRID [16], and PIP [17] databases. Pair interactions in any of the three databases were chosen to be included in our curated PPI database. As a result, a total of 561,405 pair interactions were included in our analysis. Cytoscape [18] was used to construct the interaction network. After functional enrichment analysis, the DEGs specified in significantly altered biological processes (Gene Ontology terms) and KEGG pathways were mapped to corresponding networks for interaction analysis. In the network analysis, if the connection of a gene to other genes is more than the average connection, it suggests that this gene has multiple interactions with other genes. Therefore, this gene will be considered as a hub gene compared to other genes.
Cell culture, western blot and quantitative real-time PCR analysis RAW 264.7 cells were purchased from the Cell Bank of the Chinese Academy of Sciences and cultured following ATCC conditions. Antibodies for beta-catenin (8480) and beta-actin (3700) were purchased from Cell Signaling.
The TRIzol kit (Invitrogen USA) was used to extract total RNA from cells following the manufacturer's protocol. cDNA was synthesized using the cDNA synthesis kit (Invitrogen, USA). Quantitative real-time PCR was performed using the TaqMan system (Thermo Fisher Scientific) and 18S rRNA was used as an internal control. All TaqMan probes were purchased from Thermo Fisher Scientific. All setups were performed in triplicate. The results were derived by normalization to 18S rRNA. The p-value was set at 0.05.

Results
Differential gene expression analysis of macrophages infected by S. aureus at various time points In order to identify essentially and significantly changed genes of the macrophage response to S. aureus infection, we analyzed the gene expression profile of the macrophage database which has 3 different S. aureus exposure times (GSE 13670). At |log 2 (fold change)| > 1. 5 Table 1).
Overlapping of all the DEGs from three comparisons was performed by Limma [11] to identify the commonly changed genes from different S. aureus exposure times (Fig. 1). As shown by the Venn diagram, 103 common genes were significantly changed after 8 h, 24 h or 48 h of S. aureus infection. Among these, 94 DEGs were up-regulated and 9 DEGs were down-regulated (Additional file 1: Table S1). There were 335 unique DEGs changed only in Comparison 1, among which 187 DEGs were up-regulated and 148 were downregulated ( Fig. 1, Additional file 2: Table S2). Two hundred twenty-seven DEGs changed only in Comparison 2, among which 76 DEGs were up-regulated and 151 were downregulated ( Fig. 1, Additional file 3: Table S3). Fifty-one DEGs were changed only in Comparison 3, among which 36 DEGs were up-regulated and 15 down-regulated ( Fig. 1, Additional file 4: Table S4).

Construction of biological network
Heatmaps or PPI networks were constructed and visualized for these significantly changed genes commonly or specifically changed in these three comparisons (Figs. 2, 3, 4 and 5). Several hub genes which are key members of significant pathways/networks were identified. Among 103 common DEGs, STAT1 is the hub gene in the constituted network ( Fig. 2a and b). STAT1 is involved in cell growth regulation and antiviral and immune defense, such as inflammation and innate and adaptive immunity, antiproliferative responses and tumor suppression, and it participates in crosstalk with other signal transduction pathways [19]. The expression of STAT1 was significantly increased among the common genes in our analysis. For DEGs changed only in Comparison 1, the biological network was more complicated and the key genes in the network were TP53, JAK2, CEBPA, STAT3, MYC and CTNNB1 ( Fig. 3a and b). For DEGs changed only in Comparison 2, the key gene was PRKCA ( Fig. 4a and b). No biological network was obtained for DEGs changed only in Comparison 3 (Fig. 5).

GO and KEGG pathway analysis of DEGs in macrophages infected by S. aureus with different durations of treatment
To better understand the protein and biological modules involved in the process of macrophages exposed to S. aureus, we used ClusterProfiler [14] to perform GO categories and KEGG pathways enrichment analysis with DEGs at different time points. Seven hundred forty-two biological processes and 25 KEGG pathways were enriched in the analysis of these common DEGs. Among them were significantly changed processes or pathways related to infection, such as defense response, immune system process, signal transduction and JAK-STAT signaling pathway (Table 2). For the S. aureus 8 h treatment group (Comparison 1), 338 biological processes and 9 KEGG pathways were over-represented. Consistent with the results of analysis of common DEGs, these significantly enriched processes and pathways are also related to infection, including defense response, immune system process, myeloid cell differentiation and symbiosis, encompassing mutualism through parasitism, JAK-STAT signaling pathway, apoptosis and cell cycle (Table 3). These pathways and processes might play major roles in the early stage of S. aureus infection. For 24 h treatment of the S. aureus group (Comparison 2), 62 biological processes and 7 KEGG pathways were over-represented. The included significant GO biological processes were "response to stress", "response to stimulus" and "cell activation" (Table 4). No significant KEGG pathway was identified for Comparison 2. No GO biological processes or KEGG pathways were significantly enriched after the analysis of the DEGs in Comparison 3. Cell culture, quantitative real-time PCR and western blot analyses To confirm the results obtained from bioinformatic analysis, quantitative real-time PCR was employed to detect the mRNA expression of key genes after different S. aureus exposure times in RAW 264.7 murine macrophages. As shown in Fig. 6a, the expression level of the HUB gene STAT1 was significantly increased (2.5, 3.8, and 4.5 fold respectively) at all exposure time settings compared to the control group. While a 2.2-fold increase was observed for beta-catenin (CTNNB1) after 8 h of S. aureus exposure, PRKCA was downregulated (~45% compared to control) at the 24 h time point. These results were consistent with previous analyses. Considering the early change of beta-catenin and the enriched inflammation signaling in murine macrophages after S. aureus exposure and the effect of beta-catenin on inflammation signaling [20], we are curious whether beta-catenin was involved in the inflammatory response with cytokine production regulation in a general condition (without S. aureus exposure). Here, several typical cytokines, including IL-12, IL-6, IL-4, IFNγ, and TNFα, were examined by qPCR analysis after knockdown of beta-catenin in RAW 264.7 murine macrophages in the general situation (without S. aureus exposure). Figure 6d indicates that the knockdown of beta-catenin by siRNA leads to a reduced expression level for IL-12 (55% compared to control), IL-6 (45% compared to control), IL-4 (65% compared to control), IFNγ (40% compared to control), and TNFα (45% compared to control). In contrast, over-expression of CTNNB1 in macrophages leads to the increased expression of IL-12, IL-6, IL-4, IFNγ, and TNFα (Fig. 6e).

Discussion
S. aureus is one of the most frequently isolated pathogens, with significant morbidity and mortality. The modulation of macrophages' functions by S. aureus has a significant impact on the immune response to bacterium infection [7]. Therefore, understanding the molecular mechanism of S. aureus infection leading to the cell immune response and apoptosis is of great clinical significance.
In the current study, we analyzed the gene expression profiles of macrophages exposed to S. aureus for different durations. Interestingly, we identified more DEGs during the early stage of S. aureus infection (624 DEGs and 613 DEGs in Compassion 1 and Comparison 2, respectively) compared to the later stage (253 DEGs in Comparison 3), suggesting that the macrophages underwent more significant gene expression profile changes within 24 h after the S. aureus exposure.
We identified eight potential essential genes in all three comparisons or specific in only one comparison. Among them, the hub gene STAT1 was the only one consistently upregulated from all 3 comparisons, indicating that its expression is constitutively elevated during different S. aureus infection times. As far as we know, this is the first molecular evidence that STAT1 is elevated in S. aureus infected human monocyte-derived macrophages. The increased STAT1 expression after S. aureus exposure was further confirmed by qPCR analysis in murine macrophages (RAW 264.7) at different time points. Similarly, a previous study indicated that STAT1 was involved in Cryptococcus neoformans infection in BALB/c mice [21]. The STAT1-mediated signal transduction pathway played an important role in elicitation of the classical macrophage phenotype during C. neoformans challenge [22]. Whether the expression of STAT1 is directly related to the absolute number of macrophages is not clear. S. aureus infection stimulated the expression of STAT1, which suggested that the STAT1 signaling pathway was activated to rapidly and tightly regulate the macrophage response to the infectious challenges.  The mRNA expression of six key genes was significantly changed after 8 h of S. aureus infection. Among them the expression of JAK2, STAT3 and CTNNB1 was upregulated, while TP53, CEBPA and MYC were consistently down-regulated. As a key component of the JAK2/STAT3/Snail pathway, JAK2 is associated with macrophages on tumor infiltration involved in epithelial-mesenchymal transition [23]. STAT3 mediates the expression of a variety of genes in response to cell stimuli, and plays a key role in many cellular processes, such as cell growth and apoptosis [24]. STAT3 and IL-10 play a Fig. 6 In macrophages with different S. aureus exposure times (a) STAT1 was significantly increased (8,24 and 48 h), b beta-catenin (CTNNB1) was increased in Comparison 1 (8 h), c PRKCA was down-regulated in Comparison 2 (24 h). d 48 h after CTNNB1 knockdown, the expression levels for IL-12, IL-6, IL-4, IFNγ, and TNFα were reduced in macrophages by qPCR analysis. e IL-12, IL-6, IL-4, IFNγ, and TNFα expression levels were increased after CTNNB1 over-expression in macrophages key role in driving immune dysregulation and severe immunodeficiency [25]. A recent study suggested that the absence of TP53 in endometrial cells initiates chronic inflammation, and a TP53 mutant in endometrial cancer cells induces normal macrophages to express genes that are involved in the inflammatory reaction through signal pathways [26]. In our study, after 24 h of S. aureus infection, the expression of TP53 in macrophages was down-regulated. MYC was down-regulated in our analyses; several studies have already suggested that MYC played a role in apoptosis, the cell cycle and cellular transformation [27,28]. The expression of CEBPA was also down-regulated during 8 h of S. aureus treatment. Its encoded protein can modulate the expression of genes involved in cell cycle regulation [29]. PRKCA is a serine-and threonine-specific protein kinase, which is activated by phosphatidylserine in a calcium-dependent manner. PKC family members phosphorylate a wide variety of protein targets and are known to be involved in diverse cellular signaling pathways. PRKCA was down-regulated after 24 h of S. aureus infection. To our knowledge, this is the first time CTNNB1 and PRKCA have been identified as essential genes in the relevant studies. Suppression of PRKCA might play a role in the infection response, which is similarly reported in Brucella and some other intracellular pathogens, such as Salmonella, Leishmania, and Legionella infected macrophages [30].
Our KEGG pathways analysis identified the JAK-STAT signaling pathway as the most significant process in the early stage (8 h) of S. aureus infection of macrophages. The JAK-STAT pathway is one of several important intracellular signaling pathways responsible for the activation of macrophages, the inflammatory response and inhibition of apoptosis [8,31]. Similar to our results, Zhu et al. also found that the JAK-STAT signaling pathway was activated during the phagosome maturation of macrophage exposure to S. aureus [32].
The Wnt/beta-catenin pathway might also be involved in the macrophage inflammatory response to S. aureus. As shown in Fig. 6e, the qPCR analysis in murine macrophages suggested that the elevated expression of CTNNB1 led to upregulation of IL-12, IL-6, IL-4, IFNγ, and TNFα expression levels, further highlighting the importance of CTNNB1 in signaling/activation of the inflammatory response.