Next Article in Journal
In Vitro Selection of a Single-Stranded DNA Molecular Recognition Element against Atrazine
Next Article in Special Issue
Adoptive Immunotherapy Strategies with Cytokine-Induced Killer (CIK) Cells in the Treatment of Hematological Malignancies
Previous Article in Journal
Artesunate Induces Apoptosis of Bladder Cancer Cells by miR-16 Regulation of COX-2 Expression
Previous Article in Special Issue
Transgene IL-6 Enhances DC-Stimulated CTL Responses by Counteracting CD4+25+Foxp3+ Regulatory T Cell Suppression via IL-6-Induced Foxp3 Downregulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Data Mining of Atherosclerotic Plaque Transcriptomes Predicts STAT1-Dependent Inflammatory Signal Integration in Vascular Disease

by
Krzysztof Sikorski
1,*,†,
Joanna Wesoly
2 and
Hans A. R. Bluyssen
1
1
Department of Human Molecular Genetics, Adam Mickiewicz University in Poznan, Poznan 61-614, Poland
2
Laboratory of High-Throughput Technologies, Institute of Molecular Biology and Biotechnology, Adam Mickiewicz University, Umultowska 89, Poznan 61-614, Poland
*
Author to whom correspondence should be addressed.
Current address: Institute of Clinical Medicine, Department of Immunology, University of Oslo, Sognsvannsveien 20, Oslo 0027, Norway.
Int. J. Mol. Sci. 2014, 15(8), 14313-14331; https://doi.org/10.3390/ijms150814313
Submission received: 12 June 2014 / Revised: 29 July 2014 / Accepted: 1 August 2014 / Published: 18 August 2014
(This article belongs to the Special Issue Mechanism of Action and Applications of Cytokines in Immunotherapy)

Abstract

:
Atherosclerotic plaque development involves multiple extra- and intra-cellular signals engaging cells from the immune system and from the vasculature. Pro-inflammatory pathways activated by interferon gamma (IFNγ) and toll-like receptor 4 (TLR4) ligands are profoundly involved in plaque formation and have been shown to involve cross-talk in all atheroma-interacting cell types leading to increased activation of signal transducer and activator of transcription-1 (STAT1) and elevated expression of pro-inflammatory mediators. Here we demonstrate that in Gene Expression Omnibus repository (GEO) deposited microarray datasets, obtained from human coronary and carotid atherosclerotic plaques, a significant increase in expression of pro-inflammatory and immunomodulatory genes can be detected. Moreover, increased expression of multiple chemokines, adhesion molecules and matrix-remodeling molecules was commonly detected in both plaque types and correlated with the presence of putative STAT1 binding sites in their promoters, suggesting strong involvement of STAT1 in plaque development. We also provide evidence to suggest that STAT1-nuclear factor kappa-light-chain-enhancer of activated B cells (NFκB) or STAT1-interferon-regulated factor (IRF) regulatory modules are over-represented in the promoters of these inflammatory genes, which points to a possible contribution of IFNγ and TLR4 cross-talk in the process of atherogenesis. Finally, a subset of these genes encodes for secreted proteins that could serve as a basis of a non-invasive diagnostic assay. The results of our in silico analysis in vitro provide potential evidence that STAT1-dependent IFNγ-TLR4 cross-talk plays a crucial role in coronary and carotid artery plaque development and identifies a STAT1-dependent gene signature that could represent a novel diagnostic tool to monitor and diagnose plaque progression in human atherosclerosis.

1. Introduction

Developing atherosclerotic plaque integrates multiple extra- and intracellular signals recognized by immune as well as vascular cells [1,2,3]. Activation of multiple signaling pathways in all cell types present in the atheroma leads to cross-talk at the level of cell-cell interaction but also at that of transcription factor activation. For instance, endothelial and smooth muscle cells respond to damage and pathogen associated molecular patterns, which activate Toll-like receptors (e.g., TLR4) [4]. The same ligands activate TH1 lymphocytes to release pro-inflammatory interferon (IFN)γ, which also acts on vascular cells [5]. Consequently, cell activation leads to release of inflammatory cytokines and chemokines and presentation of adhesion molecules, allowing leukocyte infiltration of the vessel wall and plaque progression [6]. It has been shown in many studies that IFNγ plays a key role in atherogenesis [7,8]. IFNγ activates an inflammatory transcriptional program through the canonical janus-kinase-signal transducer and activator of transcription (JAK)-STAT signaling pathway. This involves early activation of STAT1 (Signal Transducer and Activator of Transcription-1) and, at a later stage, of IRFs (Interferon Regulatory Factors), to generate a prolonged interferon response [9]. The immediate IFNγ response involves binding of STAT1 dimers to the Gamma Activated Sequence (GAS) element in the promoters of genes such as interferon-regulated factor 1 (IRF1), guanylate binding proteins (GBP) and intercellular adhesion molecule 1 (ICAM1). In contrast, the IRF-mediated response to IFNγ utilizes the Interferon Stimulated Response Element (ISRE) to induce expression of chemokine genes such as chemokine (C-C motif) ligand 5 CCL5 [10] or C-X-C motif chemokine 10 (CXCL10) [11].
Similar to IFNγ, TLR4 expression was detected in human and mouse atherosclerotic plaques. Moreover, circulating monocytes from acute coronary syndrome and coronary arteriosclerotic patients exhibit elevated TLR4 expression [12]. Mice deficient in TLR4 show reduced atherosclerosis proving that Toll-like receptor-dependent signaling participates in disease development [13]. Signaling of TLR4 involves two major pathways: 1. MyD88 (Myeloid Differentiation factor 88)-dependent NFκB (Nuclear Factor kappa B) activation and TRAM (TLR4 adaptor molecule)-dependent IRF3 activation. NFκB induces expression of many inflammatory genes, such as inducible nitric oxide synthase (iNOS), tumor necrosis factor alpha (TNFα) and interleukin 6 (IL6). IRF3, on the other hand, induces production of IFNβ, which in an auto- and paracrine manner activates STAT1 [14]. Thus, IFNγ and TLR4 signaling pathways utilize common transcription factors, including STAT1 and IRFs.
Over the years, cross-talk has been shown to exist between IFNγ and lipopolysaccharide (LPS) in the different cell types building the atherosclerotic plaque, resulting in increased expression of inflammatory mediators [15]. This cross-talk encompasses a complex, multi-layered mechanism relying on increased activation of STAT1 as well as on interactions of STAT1 with other transcription factors (i.e., IRF1 and NFκB), resulting in increased expression of genes such as CXCL10 and ICAM1 [11,16,17]. The interactions can be either direct at the protein level or indirect at the level of promoter binding sites. The latter requires specific regulatory modules containing STAT1 binding elements and NFκB or IRF binding sites in close proximity. Similar cross-talk phenomena have been proven to exist for other cytokine combinations (e.g., TNFα and IFNγ, and IL1β with IFNγ) and to particularly contribute to inflammatory gene expression. For instance, IFNγ and TNFα synergistically induced CXCL9 [18], ICAM1 [19] and iNOS [20], which depended on GAS and NFκB elements in their promoters and involved interaction of bound STAT1 and NFκB factors with CREB-binding protein and enhanced recruitment of RNA polymerase II [18]. Alternatively, the human IDO1 gene and murine Tap1 and Lmp2 genes were shown to possess combined ISRE and gamma-interferon-activated sites (GAS) elements in their promoters, both being necessary for maximum induction by IFNγ [21,22]. A third possible mechanism was revealed in the regulation of vascular cell adhesion molecule 1 (VCAM1) [19] and CCL19 [23] expression, where a combination of interferon-stimulated response element (ISRE) and NFκB elements appeared to be responsible for optimal transcription. In case of CXCL10, an even more complex mechanism seems to be involved, since synergistic induction by IFNγ and TNFα relied on ISRE and NFκB elements [11], whereas IFNγ-TLR4 cross-talk was suggested to depend on STAT1 [24].
Until now, the existence of these different cross-talk mechanisms has been revealed predominantly by in vitro experiments conducted in the individual cell types involved in atherosclerotic plaque formation. To date there is limited information available on the role of IFNγ and TLR4 signaling cross-talk in the regulation of pathophysiological processes underlying atheroma development. Here, by applying an in silico approach, we analyzed gene expression profiles in combination with gene ontology (GO) classification and promoter analysis of human coronary and carotid lesions (extracted from GEO: GSE40231 [25] and GSE21545 [26]) for potential evidence that STAT1-dependent inflammatory signal integration may be involved in plaque development. Indeed, our analysis highly suggests that STAT1-NFκB and STAT1-IRF regulatory modules are over-represented in promoters of inflammatory genes up-regulated in human coronary and carotid plaques and points to a possible involvement of IFNγ and TLR4 cross-talk.
Moreover, based on GO classification of these up-regulated genes, we detected high similarity in molecular processes and cellular interactions underlying plaque development in both vessel types, predicting overlap in pathophysiology. Finally, this comparative gene expression analysis revealed the presence of a common subset of inflammatory chemokine, cytokine and matrix remodeling genes, encoding for secreted proteins. These could serve as a basis of a non-invasive diagnostic assay for early detection and monitoring of the atherosclerotic process.

2. Results

2.1. STAT1 Target Genes Are Profoundly Present in Coronary Plaques

First, we analyzed a microarray dataset obtained from human coronary plaques. The dataset is available in the GEO NCBI database (acc. no. GSE40231 [25]). We evaluated 254 genes up-regulated at least two-fold as compared to the control, healthy arterial tissue. The full list of differentially regulated genes is available in supplementary data. When we examined the top 20 up-regulated genes, we immediately recognized a distinctive inflammatory transcriptional program, including proteins involved in cell-cell adhesion and trans-endothelial migration (DSC3, CDH2), cell-matrix adhesion (HAPLN1) and cytokine signaling (TNFRSF11, CYTL1, CARTPT, IL13RA2). Secreted phosphoprotein 1 (SPP1) was shown to be up-regulated in calcified lesions [27] and also to up-regulate IFNγ [28], and IL12 and SCG2 expression is stimulated by oxidized low-density lipoprotein (oxLDL) in macrophages and dendritic cells [29]. We then inspected the whole list in a more global way by conducting functional analysis and gene ontology (GO) enrichment studies, which showed significant over-representation of GO terms related to: adhesion, migration, immune response, etc. Table 1 shows the 20 terms with the strongest enrichment, a full list of these terms is available in the supplementary data file. These processes are hallmarks of atherosclerotic plaque formation. We then analyzed the presence of potential STAT1 binding sites in the promoters of the up-regulated genes and found that GAS and ISRE elements were both over-represented with z-scores of 5.14 and 3.63, respectively (Table 2). Among 62 genes possessing one of these elements we found multiple chemokines (CCL2, CCL5, CCL19, CXCL10, CXCL9), cytokines (VEGFC), adhesion molecules (integrins), and proteins involved in matrix remodeling (MMP19) (also consult the supplementary data).
Table 1. Twenty most enriched gene ontology (GO) terms in the coronary dataset.
Table 1. Twenty most enriched gene ontology (GO) terms in the coronary dataset.
GO TermGO IDlog10 p-Value
response to external stimulusGO:0009605−12.1409
biological adhesionGO:0022610−11.2716
single-multicellular organism processGO:0044707−11.2716
locomotionGO:0040011−11.0788
multicellular organismal processGO:0032501−10.5817
single-organism processGO:0044699−10.0362
developmental processGO:0032502−9.0752
anatomical structure developmentGO:0048856−8.7645
single-organism cellular processGO:0044763−8.4034
response to stimulusGO:0050896−8.2916
signallingGO:0023052−5.8633
single organism signallingGO:0044700−5.8633
localization of cellGO:0051674−5.4306
response to chemical stimulusGO:0042221−4.3116
immune system processGO:0002376−3.6253
regulation of biological qualityGO:0065008−3.4145
behaviourGO:0007610−3.2549
response to abiotic stimulusGO:0009628−2.6289
rhythmic processGO:0048511−2.3675
Table 2. Over-representation of regulatory elements in promoters of genes from all datasets. nuclear factor kappa-light-chain-enhancer of activated B cells (NFκB) and interferon-regulated factor (IRF) refer to signal transducer and activator of transcription-1 (STAT1)-NFκB and STAT1-IRF modules. The numbers shown are z-scores.
Table 2. Over-representation of regulatory elements in promoters of genes from all datasets. nuclear factor kappa-light-chain-enhancer of activated B cells (NFκB) and interferon-regulated factor (IRF) refer to signal transducer and activator of transcription-1 (STAT1)-NFκB and STAT1-IRF modules. The numbers shown are z-scores.
STAT1NFκBIRF
GASISRE
coronary5.143.635.053.66
carotid4.926.395.236.29
common3.242.062.874.86
As explained in the introduction, a distinct feature allowing for signaling cross-talk between IFNγ and TLR4 is the presence of regulatory modules in promoters of synergistically induced genes. We found that STAT1-NFκB and STAT1-IRF modules were over-represented with z-scores of 5.05 and 3.66, respectively (Table 2). Genes with putative STAT1-IRF modules again reflected mechanisms underlying plaque formation: cell adhesion (NCAM1, VCAM1, THBS1) and migration (CCL19, CCL2, CCL4 and CCRL1), matrix remodeling and calcification (ADAMTS9, SPP1, MMP19) and inflammatory signaling (IL7R, IL13RA2, TLR4, TNFRSF10B, TNFRSF11B). Putative STAT1-NFκB modules were present in a similar set of genes (also consult the supplementary data file).

2.2. Carotid Plaques Share with Coronary Plaques Similar Pro-Inflammatory Gene Expression Patterns Governed by STAT1

In addition to the coronary dataset (GSE40231) we analyzed microarrays from dataset GSE21545 [26], containing carotid plaque samples. A full list of genes with significantly changed expression is available in the supplementary data. We analyzed 845 two-fold up-regulated genes. Similar to coronary plaques we found that the top 20 genes represented inflammatory response and pathophysiology of plaque formation: matrix remodeling and calcification (SPP1, MMP12, MMP9, ACP5), markers of macrophage activation (CHI3L) disturbed composition of lipoproteins (APOC1, PLA2G7), molecules involved in leukocyte adhesion and transmigration (S100A8, S100A9), cell apoptosis (BCL2A1) and antigen processing (IFI30). We took a similar approach as in the case of coronary plaque analysis and looked at global functional changes in the dataset by analyzing enrichment in specific GO terms. We found that genes up-regulated in carotid plaques had functions involved in immune response, adhesion and migration, and response to stress (Table 3).
Table 3. Twenty most enriched GO terms in the carotid dataset.
Table 3. Twenty most enriched GO terms in the carotid dataset.
GO TermGO IDlog10 p-Value
immune system processGO:0002376−133.4342
immune responseGO:0006955−120.0516
response to stimulusGO:0050896−77.3635
response to biotic stimulusGO:0009607−55.3516
response to stressGO:0006950−54.1255
response to chemical stimulusGO:0042221−53.7212
multi-organism processGO:0051704−47.5129
signallingGO:0023052−28.9431
single organism signallingGO:0044700−28.9431
single-organism cellular processGO:0044763−25.752
single-organism processGO:0044699−25.5622
biological adhesionGO:0022610−23.9788
regulation of biological processGO:0050789−23.9706
immune effector processGO:0002252−23.7825
biological regulationGO:0065007−22.2381
interspecies interaction between organismsGO:0044419−22.1319
locomotionGO:0040011−22.0991
regulation of molecular functionGO:0065009−20.7825
regulation of biological qualityGO:0065008−18.4841
The next step was to assess possible over-representation of putative STAT1 binding sites (Table 2). We found that both GAS and ISRE elements were over-represented with z-scores of 4.92 (GAS) and 6.39 (ISRE). Genes possibly regulated by STAT1 in carotid plaques were mostly chemokines (CCL2, CCL19, CCL5, CCRL2, CCL13). Next, we checked for the presence of STAT1-NFκB and STAT1-IRF regulatory modules and found z-scores of 5.23 (STAT1-NFκB) and 6.29 (STAT1-IRF) (Table 2). STAT1-NFκB and STAT1-IRF modules were present in a similar set of genes: adhesion molecules (ITGAM, PECAM1, VCAM1, ITGB2), chemokines (CCL2, CCL5, CCL18, CCL19, CCL7, CCL8), and matrix remodeling molecules (MMP1, MMP9, MMP12, SPP1).
Together these results highly imply that there are common features between coronary and carotid plaques, possibly regulated by STAT1 and by STAT1-dependent regulatory modules. The similarity of GO terms enriched in both datasets reflects an underlying common biology of the two plaque types.

2.3. 72 Genes Are Expressed in both Coronary and Carotid Plaques and Form “Plaque Signature”

Finally, we compared genes up-regulated in carotid and coronary datasets and found 72 common genes. Among these genes there are multiple chemokines (CCL18, CCL19, CCL2, CCL4, CCL5, CXCL10, CXCL2, CXCL9), adhesion molecules (HAPLN1, THBS1, THBS2, VCAM1, ALCAM, SELE), matrix remodeling molecules (VCAN, MMP9, LAMB1, HPSE) and also genes with proven involvement in atherosclerosis (APOC1, APOE, CD55, CD69, SCG2, SPP1, TLR4, SLC16A3). When we organized these genes according to localization of their products we found that 30 were secreted and 25 were cell membrane-bound, which makes them good potential biomarkers of the plaque (Figure 1 and Table 4).
Figure 1. Localization of products of genes expressed in both plaque types. Gene products are organized according to their localization: extracellular (green), cell membrane (blue), and intracellular (red). Wider ribbons represent a known link to atherosclerosis, as determined by literature mining (see Methods).
Figure 1. Localization of products of genes expressed in both plaque types. Gene products are organized according to their localization: extracellular (green), cell membrane (blue), and intracellular (red). Wider ribbons represent a known link to atherosclerosis, as determined by literature mining (see Methods).
Ijms 15 14313 g001
Table 4. Genes up-regulated in both plaque types, whose products are secreted. The table also shows presence of putative regulatory elements in promoters, if regulation by STAT1 has been proven.
Table 4. Genes up-regulated in both plaque types, whose products are secreted. The table also shows presence of putative regulatory elements in promoters, if regulation by STAT1 has been proven.
Gene SymbolKnown STAT1 RegulationISRESTAT1-NFκBSTAT1-IRF
APOC1no+
APOEno+
CCL18no+++
CCL19no+++
CCL2yes [30,31]+++
CCL4yes [32]+++
CCL5yes [33]+++
CPno++
CXCL10yes [34]++
CXCL2no+
CXCL9yes [34]+
GZMAno++
GZMKno++
HAPLN1no+++
HPSEno
IBSPno+
IGJno++
LAMB1no+++
MMP9yes [35]++
PTX3no++
SCG2no++
SERPINE1yes [36]+++
SPINK1no++
SPP1no+++
THBS1no++
THBS2no+
TNFRSF11Bno+++
TNS3no+
VCANno+++
VEGFCyes [37]++

2.4. Genes Common for Carotid and Coronary Plaques Are Strongly Involved in Plaque Formation Processes and Could Be Regulated by STAT1

Gene ontology classification and enrichment of genes common for both types of plaques showed strong bias towards immune response, cell adhesion and migration, and response to stress (Table 5).
Table 5. Twenty most enriched GO terms obtained from the list of genes common for carotid and coronary plaques.
Table 5. Twenty most enriched GO terms obtained from the list of genes common for carotid and coronary plaques.
GO TermGO IDlog10 p-Value
biological adhesionGO:0022610−7.8861
immune system processGO:0002376−6.2765
locomotionGO:0040011−6.2765
single-organism processGO:0044699−5.5575
single-organism cellular processGO:0044763−5.2832
immune responseGO:0006955−4.8041
response to stressGO:0006950−4.6402
localization of cellGO:0051674−4.5607
response to external stimulusGO:0009605−4.1475
response to chemical stimulusGO:0042221−4.0329
multi-organism processGO:0051704−3.2677
multi-multicellular organism processGO:0044706−2.7118
regulation of biological qualityGO:0065008−2.7082
response to biotic stimulusGO:0009607−2.6917
localizationGO:0051179−2.5934
response to stimulusGO:0050896−2.121
response to abiotic stimulusGO:0009628−2.1199
single-organism developmental processGO:0044767−2.0234
Also, binding sites of STAT1 and STAT1-containing regulatory modules were over-represented with z-scores of 3.24 (GAS), 2.06 (ISRE), 2.87 (STAT1-NFκB), 4.86 (STAT1-IRF) (Table 2). All of the 72 common genes had a putative GAS element. STAT1-NFκB module was present in 31 genes from the group, and in eight it was present exclusively. STAT1-IRF module was detected in 45 genes, and 23 did not have the other module. In fact, some of these genes (CCL5, CCL19, CCL4, CXCL10, CXCL2, CXCL9 and MMP9) were actually shown to be co-regulated by STAT1, NFκB and IRFs. Fourty of the 72 genes have a known link to atherosclerosis, as determined by literature mining (Figure 1).

3. Discussion

IFNγ and TLR4-mediated signaling pathways, activated in different immunomodulatory and vascular cell types, have been implicated in plaque development and progression [8,13]. Both pathways utilize STAT1 [15] to regulate expression of inflammatory and pro-atherosclerotic genes, such as chemokines (CXCL10) and adhesion molecules (ICAM1) [38]. In our previous research we were able to show in vitro that in endothelial cells (ECs) and vascular smooth muscle cells (VSMCs) cross-talk between IFNγ and LPS exists and in vitro facilitates STAT1-dependent increase chemokine expression and monocyte to endothelial cells adhesion, a hallmark of early atherosclerosis [17]. A similar STAT1-dependent mechanism was described by others in immune cells [39]. Thus, IFNγ-TLR4 signaling cross-talk importantly controls behavior and interactions of all cells involved in atherosclerotic plaque formation, although this has not been studied in the context of atherosclerotic plaques.
Here we applied an in silico approach on deposited in GEO gene expression profiles of coronary and carotid atherosclerotic plaques. Using GO classification, we first conducted functional analysis on 254 up-regulated genes selected from the coronary dataset. This disclosed a statistically significant over-representation of genes involved in cell adhesion, migration, response to external stimulus, and immune response (Table 1). Likewise, functional analysis of 845 up-regulated genes from the carotid dataset revealed strong over-representation of GO terms such as: immune response, adhesion, migration, etc. (Table 3), which was similar to our findings in coronary plaques.
These results are in line with the current view of the different processes involved in plaque development and progression [40,41]. First, injured endothelium becomes pro-thrombotic and releases cytokines and chemokines, allowing for leukocyte adhesion and migration into the vessel wall. Then, lymphocytes differentiate into T cells (mostly of the TH1 phenotype) and monocytes into macrophages. Subsequent production of inflammatory mediators by these cells activates proliferation and migration of smooth muscle cells facilitated by secretion of matrix-metalloproteinases (by macrophages and smooth muscle cells) [1]. Together with a disturbed composition of plasma lipoproteins, a known risk factor for atherosclerosis [42], all of these processes are reflected in the over-represented GO terms listed in Table 1 and the supplementary data file.
Interestingly, we have observed that 62 out of 254 of the up-regulated genes in coronary plaques possess potential STAT1 binding sites (GAS or ISRE) in their promoters (see Results, Table 2 and supplementary data) highly implicating STAT1 in the development of atherosclerotic plaques. Genes containing either GAS or ISRE elements included chemokines (CCL2, CCL5, CCL19, CXCL10, CXCL9), cytokines (VEGFC), adhesion molecules (integrins), and proteins involved in matrix remodeling (MMP19). In the carotid plaques 208 of 845 up-regulated genes demonstrated presence of putative STAT1 binding sites and involved chemokines (CCL2, CCL19, CCL5, CCRL2, CCL13), chemokine receptors (CX3CR1, CXCR6), cytokines, and cytokine receptors (IL18, IL2RG, IL2RB). These findings are in correlation with results from in vitro studies, which showed that many of these genes are indeed STAT1 targets (e.g., CCL5 [16], CXCL10 [43]).
Since our aim was also to provide evidence for the potential role of IFNγ-TLR4 cross-talk in plaque development, we searched for the presence of transcriptional modules of STAT1 and NFκB or STAT1 and IRF binding sites (not more than 50 nucleotides apart). Indeed, we uncovered in coronary plaques that 116 genes contained a STAT1-NFκB module and 150 a STAT1-IRF module, pointing to a potential mechanism of STAT1 dependent co-regulation of gene expression in the cell types present in atherosclerotic plaques.
GO classification of up-regulated genes in coronary plaques with putative STAT1-IRF modules in their promoters reflect mechanisms underlying plaque formation: cell adhesion (NCAM1, VCAM1, THBS1) and migration (CCL19, CCL2, CCL4 and CCRL1), matrix remodeling and calcification (ADAMTS9, SPP1, MMP19), and inflammatory signaling (IL7R, IL13RA2, TLR4, TNFRSF10B, TNFRSF11B). In carotid plaques STAT1-IRF and STAT1-NFκB modules were present in a partially similar set of genes: adhesion molecules (ITGAM, PECAM1, VCAM1, ITGB2), chemokines (CCL2, CCL5, CCL18, CCL19, CCL7, CCL8), and matrix remodeling molecules (MMP1, MMP9, MMP12, SPP1).
Collectively, these results predict that key processes involved in plaque development are regulated by STAT1, either alone or in close co-operation with NFκB and IRFs providing a platform for cross-talk between different inflammatory stimuli. It also provides further proof for the crucial role of STAT1 in human atherosclerosis. This notion is further supported by an additional study from our group that was recently submitted elsewhere [Chmielewski S. et al.], in which we uncovered a subset of 30 STAT1-dependent genes that are highly amplified by co-administration of IFNγ and LPS in VSMCs in vitro [Chmielewski S. et al.]. These genes included the chemokines CCL5, CXCL10, CCL8, CXCL9, CCRL2, which were also up-regulated in carotid or coronary human plaques, as shown here, and together reflect pro-atherogenic responses in human atherosclerosis. Based on our findings and other studies from the literature, Figure 2 summarizes the complex role of STAT1 in signal integration between IFNγ and TLRs that we propose to take place in the forming atheroma.
Regulatory modules containing elements for STAT1 and IRF or NFκB have been shown to play key role in signal integration of various signaling pathways. For instance, maximal transcription of ICAM1 is possible only when both NFκB and STAT1 are bound to the promoter. The same mechanism was observed for CXCL10 in human smooth muscle cells [44] and monocytic THP-1 cells [45]. We observed a similar synergy between IFNγ and TLR4 signaling in human endothelial cells and murine smooth muscle cells, which was STAT1 dependent and led to amplified expression of CXCL10 and other genes [17] [Chmielewski S. et al.]. Furthermore, induction of iNOS expression by LPS is managed by NFκB binding, but presence of a GAS element in the promoter and co-treatment with IFNγ enables maximal expression [46]. Co-operative action of binding sites for STAT1 and IRF1 has been shown to manage expression of indoleamine 2,3-dioxygenase 1 (IDO1) (also termed INDO) gene [21], which is involved in sustaining chronic inflammation and is one of the hallmark genes of the M1 inflammatory macrophage phenotype [47]. Involvement of regulatory modules in complex signal integration of IFNγ and TLR4 (also other TLRs), in addition to other pro-inflammatory signals, has long been appreciated by the scientific community and the evidence presented here suggests that it is a functional mechanism in the plaque development contributing to chronic production of inflammatory modulators.
Figure 2. Signal integration between interferon gamma (IFNγ) and Toll-like receptors (TLRs) involves interactions of STAT1, interferon-regulated factor (IRF)1, IRF8 and nuclear factor kappa-light-chain-enhancer of activated B cells (NFκB). Analysis of gene expression profiles of atherosclerotic plaques together with the results of microarray experiments conducted by our group on STAT1 -/- SMC [Chmielewski S. et al.] create this complex signal integration picture (see Discussion). STAT1 regulates genes by itself, but also interacting with IRF1, IRF8 or NFκB.
Figure 2. Signal integration between interferon gamma (IFNγ) and Toll-like receptors (TLRs) involves interactions of STAT1, interferon-regulated factor (IRF)1, IRF8 and nuclear factor kappa-light-chain-enhancer of activated B cells (NFκB). Analysis of gene expression profiles of atherosclerotic plaques together with the results of microarray experiments conducted by our group on STAT1 -/- SMC [Chmielewski S. et al.] create this complex signal integration picture (see Discussion). STAT1 regulates genes by itself, but also interacting with IRF1, IRF8 or NFκB.
Ijms 15 14313 g002
It is known that the pathophysiology of coronary plaques resembles that of carotid plaques [48]. Prevalence of stenosis in one artery type has been associated with prevalence of stenosis in the other. Both plaque types have been shown to respond in a similar way to shear stress and both display an inflammatory background [49]. We were interested if the degree of similarity between coronary and carotid plaques could be recognized at the gene expression level. Indeed, overlap was observed in up-regulated genes from both datasets, together forming a 72-gene “plaque signature”. As expected, also distinct gene sets could be identified. This was in agreement with the fact that atherosclerosis is a multi-factorial disease, and it is highly likely that factors inducing and accelerating the disease will differ between various plaque types and even between different specimens [50,51,52].
When we characterized in more detail the biological functions of these 72 “plaque signature” genes, we could distinguish many GO terms related to plaque development processes: cell adhesion (ALCAM [53,54], VCAM1 [55]), T-cell migration (CCL5 [56], CCR5 [57], and CXCR4 [58]), response to LPS (TLR4), etc. (Table 4 and supplementary data file). In addition, promoter analysis of these genes revealed over-representation of STAT1 binding sites as well as STAT1 containing modules (Table 2). The STAT1-NFκB module was present in 31 genes of this common group while the STAT1-IRF module could be detected in 45 genes. Some of these genes (CCL5, CCL19, CCL4, CXCL10, CXCL2, CXCL9 and MMP9) actually possessed binding sites for STAT1, NFκB as well as IRFs. The idea that STAT1-dependent cross-talk between IFNγ and TLR4 potentially exists in human plaques and the discovery of a “plaque gene signature” present the possibility for development of a novel non-invasive screening assay. Indeed, with respect to the cellular localization of the proteins encoded by the 72 common signature genes many of them are secreted and localize to the extracellular space or cell membrane making them ideal serum markers of atherosclerosis (Figure 1 and Table 5). In addition, the majority has a known connection to atherosclerosis, as determined by literature mining (Figure 1).

4. Materials and Methods

4.1. Microarray Data Normalization and Analysis

GSE40231 [25], GSE21545 [26], and GSE13760 [59] datasets were downloaded from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus repository (GEO). We used a microarray dataset obtained from human coronary plaques (GSE40231) published and deposited in NCBI GEO by Hagg et al. [25]. The authors acquired 40 coronary plaques from patients, collected during coronary artery bypass surgery and analyzed gene expression profiles using Affymetrix Human Genome U133 Plus2 microarrays. This dataset also included samples obtained from inferior mesenteric arteries (IMA), which were atherosclerosis-free. In addition we analyzed a microarray dataset (GSE21545) obtained from carotid plaques collected during carotid endarterectomy and deposited by Folkersen et al. This dataset contained 124 Affymetrix Human Genome U133 Plus2 arrays.
Datasets were normalized in Chipster software using RMA algorithm [60]. Signals were log-transformed and probes were combined to genes using Chipster’s “Combine probes to genes” utility. Healthy arterial tissue microarrays from GSE13760 were used as controls for GSE21545, GSE40231 contained own controls (IMA samples). Batch effects between the combined datasets were removed using ComBat tool, a commonly used method for removing variations between batches of microarrays [61]. Then, log of fold change was calculated for GSE21545 and GSE40231. For further analysis, genes which were statistically significantly up-regulated at least two times were used. The significance was tested with empirical Bayes test adjusted by Bonferroni correction for multiple testing.

4.2. Gene Ontology Enrichment Studies

Gene lists were mapped to gene ontology terms from biological_process category and term enrichment was calculated using GOEAST advance tool: Multi-Batch genes [62]. The default settings were used for the analysis: hypergeometric test and Yekutieli multi-test adjustment method. Terms with p-values of less than 0.01 were considered significantly enriched. Next, GO terms level 1 and 2 were fed to REVIGO online tool to produce a short list of meaningful non-redundant GO terms. Tables presented in the article contain 20 most enriched terms. Lists of all terms are available in the supplementary data file.

4.3. Pathway Analysis and Literature Mining

Pathway analysis and literature mining were done using Genomatix Pathway System [63]. Lists of genes were uploaded to the system and then over-representation of genes related to specific terms based on literature (“literature mining”) or to canonical pathways (“pathway analysis”) was calculated. In case of pathway analysis, the software relies on canonical pathways from databases such as KEGG [64], Biocarta [65] and Reactome [66]. Literature mining uses PubMed database of scientific articles. In both cases, the system calculates a p-value, which indicates how strong over-representation of a certain group of genes is in the analyzed dataset. Based on these analyses a Circos plot was generated using the online tool [67].

4.4. Promoter Analysis

First, the Genomatix Gene2Promoter tool was used to retrieve human promoters corresponding to gene symbols from the ElDorado database. A total of 1.5 kb upstream and 0.5 kb downstream of the transcription start site was retrieved. The list of STAT1 target genes up-regulated by IFNγ and LPS in murine smooth muscle cells (SMCs) [Chmielewski S. et al., personal communication] was used as the starting point for promoter analysis. That list was fed into the pSCAN online promoter analysis tool in order to look for GAS, ISRE and NFκB binding sites. The software was set to analyze 950 bp upstream and 50 bp downstream of the transcription start site. PSCAN produced a list of over-represented transcription factor binding sites together with occurrences of each site and a matrix similarity score. Occurrences having the score of at least 0.8 were fed into MatDefine (Genomatix software package) to create a highly conserved matrix for each transcription factor binding site. The settings were as follows: touple size—8; no. of sequences containing touple—60%; matrix similarity score for sequence inclusion—0.9. Matrices for GAS, ISRE, and NFκB binding sites were then used in pSCAN as “user supplied matrices” to search for occurrences in genes up-regulated two-fold in atherosclerotic plaque datasets. RegionMiner (Genomatix) was used to assess co-appearance of STAT1 binding sites with either NFκB or IRF binding sites as modules. All presented z-score values were calculated against promoter background. Z-scores are a measure of distance of the tested population from the general population. A positive value means that the number of given transcription factor binding sites is found more often in the tested sample than it would be in a random sample. A negative value means that it is found less often. It is generally accepted to consider z-scores of more than two or less than minus two as significant. It corresponds to a p-value of 0.05 [68].

4.5. Supplementary Data

The enclosed supplementary data file is available in MS Excel format. It contains several sheets with: two-fold up-regulated genes in carotid and coronary plaques with indicated putative regulatory elements (GAS, ISRE, NFkB, STAT1-NFkB, STAT1-IRF); GO terms significantly enriched in each plaque dataset; a list of common genes between the datasets with indicated link to atherosclerosis (determined by literature mining) and a list of GO terms significantly enriched in the common genes. In addition, gene lists were compared to putative STAT1 targets obtained from a microarray study, in which STAT1-deficient murine vascular smooth muscle cells were treated with IFNγ and/or LPS [Chmielewski S. et al., personal communication]. Genes that came up as STAT1-targets were marked with a “+”. The lists are organized and formatted as tables, which allows for easy filtering of data by using the filter arrows in the heading of each column.

5. Conclusions

We present further proof that IFNγ-TLR cross-talk, which has been shown in vitro to exist in many immunomodulatory and vascular cell types, potentially contributes to atherosclerotic plaque development in vivo. We are aware of the limitations of our in silico study, and realize that detailed promoter mutagenesis or chromatin immunoprecipitation-polymerase chain reaction (ChIP-PCR) experiments on a selection of newly identified STAT1-target genes is necessary to provide final proof for our proposed role of STAT1 in human atherosclerotic plaque development. Nevertheless, it is tempting to speculate that targeted intervention at STAT1 in the vasculature would be beneficial in the management of vascular disease [69]. We have confirmed by data mining that carotid and coronary plaques are indeed very similar. In addition, we observed that many genes involved in plaque development could be regulated in a concerted way by STAT1, NFκB and IRFs. Further analysis of transcription factors activated by other stimuli, such as TGFβ or Angiotensin II, could provide insight into how these pathways interact with the IFNγ-TLR axis.
With the proven role of STAT1-dependent signal amplification in endothelial and vascular smooth muscle cells and also atheroma-interacting immune cells, STAT1 target genes represent promising markers of endothelial dysfunction and atherosclerosis development. Among the 72 genes that were shown to be commonly expressed in carotid and coronary plaques, the majority encoded secreted or membrane bound proteins that could potentially be detected in serum of patients [70]. In addition, many of these genes have known important roles particularly in the early stages of plaque development. Therefore, we propose that a selection of these markers could reflect a specific profile of early plaque formation, allowing for timely atherosclerosis detection, disease progression monitoring, and treatment outcome assessment by a simple blood test. In the clinic, the availability of such an assay would mean non-invasive diagnosis of atherosclerosis and the ability to assess the efficacy of drugs in individual patients. In addition, early diagnosis could possibly improve outcomes and decrease treatment costs.

Supplementary Files

Supplementary File 1

Acknowledgments

The work presented here was supported by the Polish Ministry of Science and Higher Education by grant Iuventus Plus no. 0493/IP1/2011/71 (to KS) and by the Polish National Science Centre grant no. 2012/07/B/NZ1/02710 (to HB).

Author Contributions

Krzysztof Sikorski participated in the study design, conducted analysis, and prepared the manuscript. Hans Bluyssen participated in the study design, data analysis, and manuscript preparation. Joanna Wesoly participated in the study design and manuscript preparation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Raines, E.W. The extracellular matrix can regulate vascular cell migration, proliferation, and survival: Relationships to vascular disease. Int. J. Exp. Pathol. 2000, 81, 173–182. [Google Scholar] [CrossRef]
  2. Ameriso, S.F.; Fridman, E.A.; Leiguarda, R.C.; Sevlever, G.E. Detection of helicobacter pylori in human carotid atherosclerotic plaques. Stroke J. Cereb. Circ. 2001, 32, 385–391. [Google Scholar]
  3. Sitia, S.; Tomasoni, L.; Atzeni, F.; Ambrosio, G.; Cordiano, C.; Catapano, A.; Tramontana, S.; Perticone, F.; Naccarato, P.; Camici, P.; et al. From endothelial dysfunction to atherosclerosis. Autoimmun. Rev. 2010, 9, 830–834. [Google Scholar] [CrossRef]
  4. Dauphinee, S.M.; Karsan, A. Lipopolysaccharide signaling in endothelial cells. Lab. Investig. J. Tech. Methods Pathol. 2006, 86, 9–22. [Google Scholar]
  5. Robertson, A.K.; Hansson, G.K. T cells in atherogenesis: For better or for worse? Arterioscler. Thromb. Vasc. Biol. 2006, 26, 2421–2432. [Google Scholar] [CrossRef]
  6. Hansson, G.K.; Robertson, A.K.; Soderberg-Naucler, C. Inflammation and atherosclerosis. Annu. Rev. Pathol. 2006, 1, 297–329. [Google Scholar]
  7. Hidalgo, L.G.; Halloran, P.F. Role of ifn-gamma in allograft rejection. Crit. Rev. Immunol. 2002, 22, 317–349. [Google Scholar]
  8. Tellides, G.; Tereb, D.A.; Kirkiles-Smith, N.C.; Kim, R.W.; Wilson, J.H.; Schechner, J.S.; Lorber, M.I.; Pober, J.S. Interferon-gamma elicits arteriosclerosis in the absence of leukocytes. Nature 2000, 403, 207–211. [Google Scholar]
  9. Tamura, T.; Yanai, H.; Savitsky, D.; Taniguchi, T. The irf family transcription factors in immunity and oncogenesis. Annu. Rev. Immunol. 2008, 26, 535–584. [Google Scholar]
  10. Liu, J.; Guan, X.; Ma, X. Interferon regulatory factor 1 is an essential and direct transcriptional activator for interferon {gamma}-induced rantes/ccl5 expression in macrophages. J. Biol. Chem. 2005, 280, 24347–24355. [Google Scholar]
  11. Ohmori, Y.; Hamilton, T.A. The interferon-stimulated response element and a kappa b site mediate synergistic induction of murine ip-10 gene transcription by ifn-gamma and tnf-alpha. J. Immunol. 1995, 154, 5235–5244. [Google Scholar]
  12. Methe, H.; Kim, J.O.; Kofler, S.; Weis, M.; Nabauer, M.; Koglin, J. Expansion of circulating toll-like receptor 4-positive monocytes in patients with acute coronary syndrome. Circulation 2005, 111, 2654–2661. [Google Scholar]
  13. Michelsen, K.S.; Wong, M.H.; Shah, P.K.; Zhang, W.; Yano, J.; Doherty, T.M.; Akira, S.; Rajavashisth, T.B.; Arditi, M. Lack of toll-like receptor 4 or myeloid differentiation factor 88 reduces atherosclerosis and alters plaque phenotype in mice deficient in apolipoprotein e. Proc. Natl. Acad. Sci. USA 2004, 101, 10679–10684. [Google Scholar]
  14. Sakaguchi, S.; Negishi, H.; Asagiri, M.; Nakajima, C.; Mizutani, T.; Takaoka, A.; Honda, K.; Taniguchi, T. Essential role of irf-3 in lipopolysaccharide-induced interferon-beta gene expression and endotoxin shock. Biochem. Biophys. Res. Commun. 2003, 306, 860–866. [Google Scholar]
  15. Sikorski, K.; Chmielewski, S.; Olejnik, A.; Wesoly, J.; Heemann, U.; Baumann, M.; Bluyssen, H.A. Stat1 as a central mediator of ifnγ and tlr4 signal integration in vascular dysfunction. JAK-STAT 2012, 1, 241–249. [Google Scholar] [CrossRef]
  16. Ohmori, Y.; Schreiber, R.D.; Hamilton, T.A. Synergy between interferon-gamma and tumor necrosis factor-alpha in transcriptional activation is mediated by cooperation between signal transducer and activator of transcription 1 and nuclear factor kappab. J. Biol. Chem. 1997, 272, 14899–14907. [Google Scholar]
  17. Sikorski, K.; Chmielewski, S.; Przybyl, L.; Heemann, U.; Wesoly, J.; Baumann, M.; Bluyssen, H.A. Stat1-mediated signal integration between ifnγ and lps leads to increased ec and smc activation and monocyte adhesion. Am. J. Physiol. Cell Physiol. 2011, 300, C1337–C1344. [Google Scholar]
  18. Hiroi, M.; Ohmori, Y. The transcriptional coactivator creb-binding protein cooperates with stat1 and nf-kappa b for synergistic transcriptional activation of the cxc ligand 9/monokine induced by interferon-gamma gene. J. Biol. Chem. 2003, 278, 651–660. [Google Scholar]
  19. Collins, T.; Read, M.A.; Neish, A.S.; Whitley, M.Z.; Thanos, D.; Maniatis, T. Transcriptional regulation of endothelial cell adhesion molecules: Nf-kappa b and cytokine-inducible enhancers. FASEB J. 1995, 9, 899–909. [Google Scholar]
  20. Ganster, R.W.; Guo, Z.; Shao, L.; Geller, D.A. Differential effects of tnf-alpha and ifn-gamma on gene transcription mediated by nf-kappab-stat1 interactions. J. Interferon Cytokine Res. 2005, 25, 707–719. [Google Scholar]
  21. Chon, S.Y.; Hassanain, H.H.; Gupta, S.L. Cooperative role of interferon regulatory factor 1 and p91 (stat1) response elements in interferon-gamma-inducible expression of human indoleamine 2,3-dioxygenase gene. J. Biol. Chem. 1996, 271, 17247–17252. [Google Scholar]
  22. Brucet, M.; Marques, L.; Sebastian, C.; Lloberas, J.; Celada, A. Regulation of murine tap1 and lmp2 genes in macrophages by interferon gamma is mediated by stat1 and irf-1. Genes Immun. 2004, 5, 26–35. [Google Scholar]
  23. Pietila, T.E.; Veckman, V.; Lehtonen, A.; Lin, R.; Hiscott, J.; Julkunen, I. Multiple nf-kappab and ifn regulatory factor family transcription factors regulate ccl19 gene expression in human monocyte-derived dendritic cells. J. Immunol. 2007, 178, 253–261. [Google Scholar]
  24. Tamassia, N.; Calzetti, F.; Ear, T.; Cloutier, A.; Gasperini, S.; Bazzoni, F.; McDonald, P.P.; Cassatella, M.A. Molecular mechanisms underlying the synergistic induction of cxcl10 by lps and ifn-gamma in human neutrophils. Eur. J. Immunol. 2007, 37, 2627–2634. [Google Scholar]
  25. Hagg, S.; Skogsberg, J.; Lundstrom, J.; Noori, P.; Nilsson, R.; Zhong, H.; Maleki, S.; Shang, M.M.; Brinne, B.; Bradshaw, M.; et al. Multi-organ expression profiling uncovers a gene module in coronary artery disease involving transendothelial migration of leukocytes and lim domain binding 2: The stockholm atherosclerosis gene expression (stage) study. PLoS Genet. 2009, 5, e1000754. [Google Scholar] [CrossRef]
  26. Folkersen, L.; Persson, J.; Ekstrand, J.; Agardh, H.E.; Hansson, G.K.; Gabrielsen, A.; Hedin, U.; Paulsson-Berne, G. Prediction of ischemic events on the basis of transcriptomic and genomic profiling in patients undergoing carotid endarterectomy. Mol. Med. 2012, 18, 669–675. [Google Scholar]
  27. Canfield, A.E.; Farrington, C.; Dziobon, M.D.; Boot-Handford, R.P.; Heagerty, A.M.; Kumar, S.N.; Roberts, I.S. The involvement of matrix glycoproteins in vascular calcification and fibrosis: An immunohistochemical study. J. Pathol. 2002, 196, 228–234. [Google Scholar]
  28. Cantor, H.; Shinohara, M.L. Regulation of t-helper-cell lineage development by osteopontin: The inside story. Nat. Rev. Immunol. 2009, 9, 137–141. [Google Scholar]
  29. Chen, T.; Yan, H.; Li, Z.; Jing, T.; Zhu, W.; Ge, J.; Zheng, X.; Pan, X.; Zhu, J. Microrna-155 regulates lipid uptake, adhesion/chemokine marker secretion and scg2 expression in oxldl-stimulated dendritic cells/macrophages. Int. J. Cardiol. 2011, 147, 446–447. [Google Scholar]
  30. Scharl, M.; Hruz, P.; McCole, D.F. Protein tyrosine phosphatase non-receptor type 2 regulates ifn-gamma-induced cytokine signaling in thp-1 monocytes. Inflamm. Bowel Dis. 2010, 16, 2055–2064. [Google Scholar]
  31. Kok, S.H.; Hong, C.Y.; Kuo, M.Y.; Wang, C.C.; Hou, K.L.; Lin, Y.T.; Galson, D.L.; Lin, S.K. Oncostatin m-induced ccl2 transcription in osteoblastic cells is mediated by multiple levels of stat-1 and stat-3 signaling: An implication for the pathogenesis of arthritis. Arthritis Rheumatol. 2009, 60, 1451–1462. [Google Scholar]
  32. Dai, X.; Sayama, K.; Tohyama, M.; Shirakata, Y.; Yang, L.; Hirakawa, S.; Tokumaru, S.; Hashimoto, K. The nf-kappab, p38 mapk and stat1 pathways differentially regulate the dsrna-mediated innate immune responses of epidermal keratinocytes. Int. Immunol. 2008, 20, 901–909. [Google Scholar]
  33. Fleetwood, A.J.; Dinh, H.; Cook, A.D.; Hertzog, P.J.; Hamilton, J.A. Gm-csf- and m-csf-dependent macrophage phenotypes display differential dependence on type i interferon signaling. J. Leukoc. Biol. 2009, 86, 411–421. [Google Scholar]
  34. Hardison, S.E.; Herrera, G.; Young, M.L.; Hole, C.R.; Wozniak, K.L.; Wormley, F.L., Jr. Protective immunity against pulmonary cryptococcosis is associated with stat1-mediated classical macrophage activation. J. Immunol. 2012, 189, 4060–4068. [Google Scholar]
  35. Tyagi, A.; Agarwal, C.; Dwyer-Nield, L.D.; Singh, R.P.; Malkinson, A.M.; Agarwal, R. Silibinin modulates tnf-alpha and ifn-gamma mediated signaling to regulate cox2 and inos expression in tumorigenic mouse lung epithelial lm2 cells. Mol. Carcinog. 2012, 51, 832–842. [Google Scholar]
  36. Kosaka, H.; Yoshimoto, T.; Fujimoto, J.; Nakanishi, K. Interferon-gamma is a therapeutic target molecule for prevention of postoperative adhesion formation. Nat. Med. 2008, 14, 437–441. [Google Scholar]
  37. Eastabrook, G.D.; Hu, Y.; Tan, R.; Dutz, J.P.; Maccalman, C.D.; von Dadelszen, P. Decidual nk cell-derived conditioned medium (dnk-cm) mediates vegf-c secretion in extravillous cytotrophoblasts. Am. J. Reprod. Immunol. 2012, 67, 101–111. [Google Scholar]
  38. Hu, X.; Ivashkiv, L.B. Cross-regulation of signaling pathways by interferon-gamma: Implications for immune responses and autoimmune diseases. Immunity 2009, 31, 539–550. [Google Scholar]
  39. Schroder, K.; Sweet, M.J.; Hume, D.A. Signal integration between ifngamma and tlr signalling pathways in macrophages. Immunobiology 2006, 211, 511–524. [Google Scholar]
  40. Cagnin, S.; Biscuola, M.; Patuzzo, C.; Trabetti, E.; Pasquali, A.; Laveder, P.; Faggian, G.; Iafrancesco, M.; Mazzucco, A.; Pignatti, P.F.; et al. Reconstruction and functional analysis of altered molecular pathways in human atherosclerotic arteries. BMC Genomics 2009, 10, 13. [Google Scholar]
  41. Hansson, G.K. Immune mechanisms in atherosclerosis. Arterioscler. Thromb. Vasc. Biol. 2001, 21, 1876–1890. [Google Scholar]
  42. Poston, R.N.; Davies, D.F. Immunity and inflammation in the pathogenesis of atherosclerosis. A review. Atherosclerosis 1974, 19, 353–367. [Google Scholar] [CrossRef]
  43. Toshchakov, V.; Jones, B.W.; Perera, P.Y.; Thomas, K.; Cody, M.J.; Zhang, S.; Williams, B.R.; Major, J.; Hamilton, T.A.; Fenton, M.J.; et al. Tlr4, but not tlr2, mediates ifn-beta-induced stat1alpha/beta-dependent gene expression in macrophages. Nat. Immunol. 2002, 3, 392–398. [Google Scholar] [CrossRef]
  44. Clarke, D.L.; Clifford, R.L.; Jindarat, S.; Proud, D.; Pang, L.; Belvisi, M.; Knox, A.J. Tnfalpha and ifngamma synergistically enhance transcriptional activation of cxcl10 in human airway smooth muscle cells via stat-1, nf-kappab, and the transcriptional coactivator creb-binding protein. J. Biol. Chem. 2010, 285, 29101–29110. [Google Scholar]
  45. Qi, X.F.; Kim, D.H.; Yoon, Y.S.; Jin, D.; Huang, X.Z.; Li, J.H.; Deung, Y.K.; Lee, K.J. Essential involvement of cross-talk between ifn-gamma and tnf-alpha in cxcl10 production in human thp-1 monocytes. J. Cell. Physiol. 2009, 220, 690–697. [Google Scholar]
  46. Gao, J.; Morrison, D.C.; Parmely, T.J.; Russell, S.W.; Murphy, W.J. An interferon-gamma-activated site (gas) is necessary for full expression of the mouse inos gene in response to interferon-gamma and lipopolysaccharide. J. Biol. Chem. 1997, 272, 1226–1230. [Google Scholar]
  47. Butcher, M.J.; Galkina, E.V. Phenotypic and functional heterogeneity of macrophages and dendritic cell subsets in the healthy and atherosclerosis-prone aorta. Front. Physiol. 2012, 3, 44. [Google Scholar]
  48. Rohde, L.E.; Lee, R.T. Pathophysiology of atherosclerotic plaque development and rupture: An overview. Semin. Vasc. Med. 2003, 3, 347–354. [Google Scholar]
  49. Jashari, F.; Ibrahimi, P.; Nicoll, R.; Bajraktari, G.; Wester, P.; Henein, M.Y. Coronary and carotid atherosclerosis: Similarities and differences. Atherosclerosis 2013, 227, 193–200. [Google Scholar]
  50. Prasad, A.; Zhu, J.; Halcox, J.P.; Waclawiw, M.A.; Epstein, S.E.; Quyyumi, A.A. Predisposition to atherosclerosis by infections: Role of endothelial dysfunction. Circulation 2002, 106, 184–190. [Google Scholar]
  51. Xu, Q. Role of heat shock proteins in atherosclerosis. Arterioscler. Thromb. Vasc. Biol. 2002, 22, 1547–1559. [Google Scholar]
  52. Frantz, S.; Ertl, G.; Bauersachs, J. Mechanisms of disease: Toll-like receptors in cardiovascular disease. Nat. Clin. Pract. Cardiovasc. Med. 2007, 4, 444–454. [Google Scholar]
  53. Rauch, S.J.; Rosenkranz, A.C.; Bohm, A.; Meyer-Kirchrath, J.; Hohlfeld, T.; Schror, K.; Rauch, B.H. Cholesterol induces apoptosis-associated loss of the activated leukocyte cell adhesion molecule (alcam) in human monocytes. Vasc. Pharmacol. 2011, 54, 93–99. [Google Scholar]
  54. Kurata, H.; Matsumoto, A.; Fujiwara, Y.; Kondo, K.; Itakura, H.; Mitchell, A.; Fidge, N. A candidate high density lipoprotein (hdl) receptor, hb2, with possible multiple functions shows sequence homology with adhesion molecules. J. Atheroscler. Thromb. 1998, 4, 112–117. [Google Scholar] [CrossRef]
  55. Barks, J.L.; McQuillan, J.J.; Iademarco, M.F. Tnf-alpha and il-4 synergistically increase vascular cell adhesion molecule-1 expression in cultured vascular smooth muscle cells. J. Immunol. 1997, 159, 4532–4538. [Google Scholar]
  56. Murooka, T.T.; Rahbar, R.; Platanias, L.C.; Fish, E.N. Ccl5-mediated t-cell chemotaxis involves the initiation of mrna translation through mtor/4e-bp1. Blood 2008, 111, 4892–4901. [Google Scholar]
  57. Yurchenko, E.; Tritt, M.; Hay, V.; Shevach, E.M.; Belkaid, Y.; Piccirillo, C.A. Ccr5-dependent homing of naturally occurring cd4+ regulatory t cells to sites of leishmania major infection favors pathogen persistence. J. Exp. Med. 2006, 203, 2451–2460. [Google Scholar]
  58. Ottoson, N.C.; Pribila, J.T.; Chan, A.S.; Shimizu, Y. Cutting edge: T cell migration regulated by cxcr4 chemokine receptor signaling to zap-70 tyrosine kinase. J. Immunol. 2001, 167, 1857–1861. [Google Scholar]
  59. Cangemi, C.; Skov, V.; Poulsen, M.K.; Funder, J.; Twal, W.O.; Gall, M.A.; Hjortdal, V.; Jespersen, M.L.; Kruse, T.A.; Aagard, J.; et al. Fibulin-1 is a marker for arterial extracellular matrix alterations in type 2 diabetes. Clin. Chem. 2011, 57, 1556–1565. [Google Scholar]
  60. Kallio, M.A.; Tuimala, J.T.; Hupponen, T.; Klemela, P.; Gentile, M.; Scheinin, I.; Koski, M.; Kaki, J.; Korpelainen, E.I. Chipster: User-friendly analysis software for microarray and other high-throughput data. BMC Genomics 2011, 12, 507. [Google Scholar]
  61. Johnson, W.E.; Li, C.; Rabinovic, A. Adjusting batch effects in microarray expression data using empirical bayes methods. Biostatistics 2007, 8, 118–127. [Google Scholar]
  62. Zheng, Q.; Wang, X.J. Goeast: A web-based software toolkit for gene ontology enrichment analysis. Nucleic Acids Res. 2008, 36, W358–W363. [Google Scholar]
  63. Data Analysis Suite GENOMATIX. Available online: www.genomatix.de (accessed on 14 August 2014).
  64. Kanehisa, M.; Goto, S. Kegg: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar]
  65. BioCarta. Available online: http://www.biocarta.com (accessed on 20 February 2014).
  66. Croft, D.; O’Kelly, G.; Wu, G.; Haw, R.; Gillespie, M.; Matthews, L.; Caudy, M.; Garapati, P.; Gopinath, G.; Jassal, B.; et al. Reactome: A database of reactions, pathways and biological processes. Nucleic Acids Res. 2011, 39, D691–D697. [Google Scholar]
  67. Krzywinski, M.; Schein, J.; Birol, I.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar]
  68. Genomatix Regionminer Manual. Available online: http://www.genomatix.de/online_help/help_regionminer/overrepresented_TFs.html (accessed on 20 February 2014).
  69. Sikorski, K.; Czerwoniec, A.; Bujnicki, J.M.; Wesoly, J.; Bluyssen, H.A. Stat1 as a novel therapeutical target in pro-atherogenic signal integration of ifngamma, tlr4 and il-6 in vascular disease. Cytokine Growth Factor Rev. 2011, 22, 211–219. [Google Scholar]
  70. Kuvin, J.T.; Karas, R.H. Clinical utility of endothelial function testing: Ready for prime time? Circulation 2003, 107, 3243–3247. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Sikorski, K.; Wesoly, J.; Bluyssen, H.A.R. Data Mining of Atherosclerotic Plaque Transcriptomes Predicts STAT1-Dependent Inflammatory Signal Integration in Vascular Disease. Int. J. Mol. Sci. 2014, 15, 14313-14331. https://doi.org/10.3390/ijms150814313

AMA Style

Sikorski K, Wesoly J, Bluyssen HAR. Data Mining of Atherosclerotic Plaque Transcriptomes Predicts STAT1-Dependent Inflammatory Signal Integration in Vascular Disease. International Journal of Molecular Sciences. 2014; 15(8):14313-14331. https://doi.org/10.3390/ijms150814313

Chicago/Turabian Style

Sikorski, Krzysztof, Joanna Wesoly, and Hans A. R. Bluyssen. 2014. "Data Mining of Atherosclerotic Plaque Transcriptomes Predicts STAT1-Dependent Inflammatory Signal Integration in Vascular Disease" International Journal of Molecular Sciences 15, no. 8: 14313-14331. https://doi.org/10.3390/ijms150814313

APA Style

Sikorski, K., Wesoly, J., & Bluyssen, H. A. R. (2014). Data Mining of Atherosclerotic Plaque Transcriptomes Predicts STAT1-Dependent Inflammatory Signal Integration in Vascular Disease. International Journal of Molecular Sciences, 15(8), 14313-14331. https://doi.org/10.3390/ijms150814313

Article Metrics

Back to TopTop