Online ISSN: 1945-4589
Research Paper Volume 13, Issue 4 pp 5650—5673

Synergistic inflammatory signaling by cGAS may be involved in the development of atherosclerosis

Guan-Feng Lu1, *, , Sheng-Cai Chen1, *, , Yuan-Peng Xia1, *, , Zi-Ming Ye1,2, , Fei Cao1, , Bo Hu1, ,

  • 1 Department of Neurology, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan 430022, China
  • 2 Department of Neurology, The First Affiliated Hospital, Guangxi Medical University, Nanning 530021, Guangxi, China
* Equal contribution

Received: May 26, 2020       Accepted: October 23, 2020       Published: February 11, 2021      

https://doi.org/10.18632/aging.202491
How to Cite

Copyright: © 2021 Lu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Inappropriate activation or overactivation of cyclic GMP-AMP synthase (cGAS) by double-stranded deoxyribonucleic acid (dsDNA) initiates a regulatory signaling cascade triggering a variety of inflammatory responses, which are a great threat to human health. This study focused on identifying the role of cGAS in atherosclerosis and its potential mechanisms. The relationship between cGAS and atherosclerosis was identified in an ApoE -/- mouse model. Meanwhile, RNA sequencing (RNA-seq) analysis of the underlying mechanisms of atherosclerosis in RAW264.7 macrophages treated with cGAS inhibition was conducted. Results showed that cGAS was positively correlated with atherosclerotic plaque area, and was mainly distributed in macrophages. RNA-seq analysis revealed that inflammatory response, immune response and cytokine–cytokine receptor interaction may play important roles in the development of atherosclerosis. Real-time quantitative polymerase chain reaction (RT-qPCR) results showed that the expression of the pro-inflammatory factors, signal transducer and activator of transcription (Stat), interferon regulatory factor (Irf), toll-like receptors (Tlrs), and type I interferons (Ifns) were synergistically reduced when cGAS was inhibited. Furthermore, cGAS inhibition significantly inhibited RAW264.7 macrophage M1 polarization. These results demonstrate that cGAS may contribute to the development of atherosclerosis through synergistic inflammatory signaling of TLRs, STAT/IRF as well as IFNs, leading to macrophage M1 polarization.

Introduction

The identification of exogenous deoxyribonucleic acid (DNA) is the basic function of host defense. Intracellular signaling induced by innate DNA sensing can lead to strong anti-infective immune responses [1]. Cyclic GMP–AMP synthase (cGAS) located in the cytoplasm can bind to double-stranded DNA (dsDNA). Activated cGAS can catalyze the synthesis of cyclic dinucleotide GMP–AMP (cGAMP), which mediates the canonical stimulator of interferon response cGAMP interactor (STING), and the noncanonical signaling pathways [2, 3]. Although cGAS was initially identified as a critical part of anti-infective immune defense [4, 5], recent research has revealed that mis-localized self-DNA triggers improper activation or overactivation of cGAS and thus provokes endogenous inflammation, ultimately, leading to disease [3, 6–9].

Atherosclerosis, the primary etiological factor of cardio-cerebrovascular diseases which are the leading cause of death worldwide, is a chronic inflammatory process with characteristic modified low-density lipoprotein (mLDL) deposition which drives the recruitment of circulating immune cells that subsequently trigger inflammatory cascades [10]. Great efforts have been made to elucidate the mechanisms of atherosclerosis, but nevertheless, it has not yet been fully clarified. Accumulating evidence has shown that mis-localized self-DNA and DNA damage occur in atherosclerosis. The Health 2000 Survey showed that in women without hormone replacement therapy (HRT), higher circulating cell-free DNA levels indicate worse arterial elasticity [11]. Mitochondrial DNA(mtDNA), a small double-stranded circular form of DNA, is increasingly recognized as a potent stimulus in innate immune responses and inflammation [12]. A study carried out by Zhang et al. demonstrated that human atherosclerotic plasma contains higher concentrations of mtDNA than those in healthy controls [13]. Furthermore, mtDNA is released into the cytoplasm under oxidized LDL(oxLDL) stimulation [14]. In addition, it is reported that increased plaque necrosis core and induction by mtDNA damage are attributable to apoptosis of smooth muscle cells and monocytes rather than reactive oxygen species [15]. Furthermore, mtDNA can directly provoke inflammatory response by engaging the cGAS-STING pathway. Besides, extracellular release of oxidized mtDNA is also demonstrated to be pro-inflammatory in vitro, and type-I interferon (IFN) signaling was stimulated in response to the DNA sensor when oxidized mtDNA is injected into mice [16]. Hypothetically therefore, cGAS may play a pivotal role in the development of atherosclerosis. Exploring the role of cGAS in atherosclerosis and the possible underlying mechanisms may provide a new potential therapeutic target for atherosclerosis.

Bioinformatics technology is increasingly applied to unearth the potential targets of diseases, thereby enabling researchers to identify the underlying mechanisms. Macrophages are critical in immune inflammation and play a central role in the development of atherosclerosis, and a study carried out by Bai et al. found that the cGAS-cGAMP-STING pathway is activated in macrophages from high-fat diet (HFD)-induced obese mice [17]. These imply that cGAS may regulate atherosclerosis through macrophages.

This study focused on exploring the relationship between cGAS and atherosclerosis using in vivo and in vitro atherosclerosis model experiments and identifying the underlying differentially-expressed genes (DEGs) and signaling of cGAS inhibition in macrophages using RNA-sequencing (RNA-seq) analysis, which were associated with atherosclerosis.

Results

Mitochondrial DNA damage in atherosclerosis

It is well known that mtDNA damage correlates with atherosclerosis progression [15, 18, 19]. DNA with oxidative damage is resistant to degradation by cytosolic nuclease three prime repair exonuclease 1(TREX1), and enhances STING-dependent immune sensing [20]. Released oxidized mtDNA acts as a powerful inflammatory stimulus [16]. In accordance with previous findings, immunofluorescence revealed that 8-Oxo-2′-deoxyguanosine (8-OH-dG), a marker for oxidative damage to DNA, can be detected in the plaque of ApoE -/- mice and strong 8-OH-dG staining was discovered mainly in mitochondria, irrespective of diet. In addition, part of the 8-OH-dG staining overlapped with neither nucleus nor mitochondria (Figure 1A), suggesting that oxidized DNA can be released into the cytoplasm or outside the cell. Additionally, human atherosclerotic plasma contained higher concentrations of dsDNA than those in healthy controls (Figure 1B), consistent with results demonstrated by previous studies [11, 13]. This released DNA could be the actuator for the development of atherosclerosis.

Mitochondrial DNA damage in atherogenesis. ApoE -/- mice were fed a western-type diet or a chow diet for 16 weeks. Immunofluorescence was used to analyze aortic root plaques. (A) Oxidative damage to DNA was measured by immunostaining of 8-OH-dG (green) (n = 6, Scale bars: 100μm for 20× images and 10μm for 80× images). Quantitative analysis of 8-OH-dG+ cells in selected areas showed that strong 8-OH-dG staining was discovered mainly in mitochondria (TOMM20, red) compared with nuclei (DAPI, blue) (P = 0.037 in the chow diet group, and P = 0.018 in the western-type diet group). (B) The plasma dsDNA levels of atherosclerosis patients and volunteers were analyzed using a PicoGreen ® dsDNA quantitative kit. ***P

Figure 1. Mitochondrial DNA damage in atherogenesis. ApoE -/- mice were fed a western-type diet or a chow diet for 16 weeks. Immunofluorescence was used to analyze aortic root plaques. (A) Oxidative damage to DNA was measured by immunostaining of 8-OH-dG (green) (n = 6, Scale bars: 100μm for 20× images and 10μm for 80× images). Quantitative analysis of 8-OH-dG+ cells in selected areas showed that strong 8-OH-dG staining was discovered mainly in mitochondria (TOMM20, red) compared with nuclei (DAPI, blue) (P = 0.037 in the chow diet group, and P = 0.018 in the western-type diet group). (B) The plasma dsDNA levels of atherosclerosis patients and volunteers were analyzed using a PicoGreen ® dsDNA quantitative kit. ***P< 0.001. Data are means ± SD. 8-OH-dG, 8-Oxo-2′-deoxyguanosine; SD, standard deviation.

Differences in cGAS expression identified by database analysis based on the Gene Expression Omnibus

Although cGAS acts mainly as a DNA sensor and the role in broadly eliciting immune responses and inflammation has also been reported, whether cGAS contributes to atherosclerosis is yet to be uncovered. In the present study, the dataset GSE40156 from the Gene Expression Omnibus (GEO) was re-analyzed by bioinformatics methods and it was found that the expression of cGAS was increased in the total aorta of ApoE -/- mice. Furthermore, cGAS expression in ApoE -/- mice was significantly increased with age (Table 1). It is well known that macrophages are of great importance in the initiation, progression and regression of atherosclerosis and different subtypes of macrophages affect the outcome of atherosclerosis [21]. To explore the difference of cGAS in macrophage polarization, GSE57614 from the GEO was assayed and it was found that cGAS expression remarkably increased in M1 macrophages when compared with M2 macrophages (Table 2). Collectively, the results of these bioinformatics analyses results mentioned above imply that cGAS may play a pivotal role in atherosclerosis.

Table 1. cGAS expression in aorta.

Comparisonlog2FCP value
ApoE-Aorta-6-weeks vs wt-Aorta-6-weeks1.2040.947
ApoE-Aorta-32-weeks vs wt-Aorta-32-weeks0.2940.599
ApoE-Aorta-78-weeks vs wt-Aorta-78-weeks1.3730.067
ApoE-Aorta-32-weeks vs ApoE-Aorta-6-weeks1.1170.073
ApoE-Aorta-78-weeks vs ApoE-Aorta-6-weeks2.6710.000
# 0.000 means < 0.0001.

Table 2. cGAS expression in macrophage polarization.

Comparisonlog2FCP value
M2a macrophage 6h vs M1 macrophage 6h-3.2000.000
M2c macrophage 6h vs M1 macrophage 6h-2.5600.001
M2a macrophage 12h vs M1 macrophage 12h-2.6100.000
M2c macrophage 12h vs M1 macrophage 12h-2.2500.001
M2a macrophage 24h vs M1 macrophage 24h-2.700.001
M2c macrophage 24h vs M1 macrophage 24h-2.4200.002
# Human macrophage M1 polarization was obtained by Lipopolysaccharide (LPS) and interferon gamma (IFNγ), M2a by IL4, whereas IL10 induced a “deactivated” state (M2c). 0.000 means < 0.0001.

cGAS as a novel factor for modulating atherosclerosis

To explore the effects of cGAS on atherosclerosis, ApoE -/- mice were kept on a western-type or chow diet for 16 weeks. Results demonstrated that cGAS was expressed in the plaque of ApoE -/- mice regardless of diet, and was positively correlated with plaque area, and mainly distributed in macrophages (Figure 2A, 2B), whereas, the protein expression of cGAS in the aorta of ApoE -/- mice fed a western-type diet was higher than that of ApoE -/- mice fed a chow diet (Figure 2C, 2D). Importantly, cGAS inhibition reduced lipid deposition and foam cell formation in RAW264.7 macrophages (Figure 3A, 3B). These data demonstrated that cGAS contributed to atherosclerosis, and suggested that exploring the underlying mechanism could be of great significance.

cGAS was positively correlated with atherosclerosis. ApoE -/- mice were fed either a western-type diet or a chow diet for 16 weeks. (A) Immunostaining of cGAS (green) and the macrophage marker MOMA-2 (red), the smooth muscle cell marker α-SMA (red), the endothelial cell marker CD31 (red) and their co-localization in atherosclerotic plaques. (n = 6, Scale bar: 100 μm). (B) Histological analysis of the aortic root stained with hematoxylin and eosin. (C) Oil Red-O analysis of lesion area as a percentage of in total aortic area in ApoE -/- mice. (D) The protein expression of cGAS in aorta was detected by western blotting. * P P P P

Figure 2. cGAS was positively correlated with atherosclerosis. ApoE -/- mice were fed either a western-type diet or a chow diet for 16 weeks. (A) Immunostaining of cGAS (green) and the macrophage marker MOMA-2 (red), the smooth muscle cell marker α-SMA (red), the endothelial cell marker CD31 (red) and their co-localization in atherosclerotic plaques. (n = 6, Scale bar: 100 μm). (B) Histological analysis of the aortic root stained with hematoxylin and eosin. (C) Oil Red-O analysis of lesion area as a percentage of in total aortic area in ApoE -/- mice. (D) The protein expression of cGAS in aorta was detected by western blotting. * P < 0.05, ** P < 0.01, *** P < 0.001, **** P < 0.0001. Data are means ± SD.

cGAS inhibition reduced lipid deposition and foam cell formation in RAW264.7 macrophages. RAW264.7 cells were treated for 24 hours with oxLDL(100 μg/mL) and with or without RU.521(cGAS inhibitor 2 μg/mL). (A) Oil Red-O staining or (B) BODIPY staining were used to identify lipid deposition and foam cell formation. RAW264.7 macrophages were treated with RU.521 (2 μg/mL) for 24 h before analysis. (C) Cholesterol uptake was determined using DiI-oxLDL (Scale bar: 20 μm). (D) HDL-mediated cholesterol efflux was assayed using NBD 485/535 (Scale bar: 20 μm). *** P P

Figure 3. cGAS inhibition reduced lipid deposition and foam cell formation in RAW264.7 macrophages. RAW264.7 cells were treated for 24 hours with oxLDL(100 μg/mL) and with or without RU.521(cGAS inhibitor 2 μg/mL). (A) Oil Red-O staining or (B) BODIPY staining were used to identify lipid deposition and foam cell formation. RAW264.7 macrophages were treated with RU.521 (2 μg/mL) for 24 h before analysis. (C) Cholesterol uptake was determined using DiI-oxLDL (Scale bar: 20 μm). (D) HDL-mediated cholesterol efflux was assayed using NBD 485/535 (Scale bar: 20 μm). *** P < 0.001, **** P < 0.0001. Data are means ± SD.

Transcriptome profiling of RU.521-treated RAW264.7 macrophages

In order to have a better view of the mechanisms of atherosclerosis mediated by cGAS, we performed RNA-seq transcriptomic analysis of the underlying DEGs and the mechanisms in RAW264.7 macrophages treated for 12 h with or without RU.521, a cGAS inhibitor. The genes expressed were sequenced using the Illumia platform. DEGs between the cGAS inhibition group and the control group were detected. First, the intra-group samples showed a strong correlation (Figure 4A). Further, the repeatability of intra-group data based on principal component analysis (PCA) was acceptable, where the distances of intra-group samples were close. (Figure 4B).

Transcriptome profiling of RU.521-treated RAW264.7 macrophages. (A) Pearson’s correlation analysis of samples. The color reflects the intensity of the correlation. The higher the correlation coefficient between samples, the closer the expression pattern. (B) PCA analysis of samples. Principal component 1 (PC1) and principal component 2 (PC2) are indicated on the X-axis and Y-axis, respectively. Cos2 stands for variables. The closer the two samples were to each other, the smaller the difference was between the two samples in terms of gene expression patterns. (C) Volcano plot of the DEGs. The X-axis represents the difference in log 2 conversion, the Y-axis represents the significant difference after log 10 conversion, the blue represents the down-regulated DEGs, the red represents the up-regulated DEGs, and the gray represents the non-DEGs. The DEGs were defined according to the criteria of log2FC ≥2 and Padj D) Bar graph showing DEGs. The X-axis represents the type, the Y-axis represents gene counts, the green represents the downregulated DEGs, the blue represents the upregulated DEGs, and the red is the non-DEGs. (E) Hierarchical clustering heatmap. Red and green represent high and low relative expression, respectively. Rows represent DEGs, and columns represent samples. DEGs, differentially-expressed genes; PCA, principal component analysis.

Figure 4. Transcriptome profiling of RU.521-treated RAW264.7 macrophages. (A) Pearson’s correlation analysis of samples. The color reflects the intensity of the correlation. The higher the correlation coefficient between samples, the closer the expression pattern. (B) PCA analysis of samples. Principal component 1 (PC1) and principal component 2 (PC2) are indicated on the X-axis and Y-axis, respectively. Cos2 stands for variables. The closer the two samples were to each other, the smaller the difference was between the two samples in terms of gene expression patterns. (C) Volcano plot of the DEGs. The X-axis represents the difference in log 2 conversion, the Y-axis represents the significant difference after log 10 conversion, the blue represents the down-regulated DEGs, the red represents the up-regulated DEGs, and the gray represents the non-DEGs. The DEGs were defined according to the criteria of log2FC ≥2 and Padj <0.005. (D) Bar graph showing DEGs. The X-axis represents the type, the Y-axis represents gene counts, the green represents the downregulated DEGs, the blue represents the upregulated DEGs, and the red is the non-DEGs. (E) Hierarchical clustering heatmap. Red and green represent high and low relative expression, respectively. Rows represent DEGs, and columns represent samples. DEGs, differentially-expressed genes; PCA, principal component analysis.

A total of 1,462 DEGs were identified with the criteria of log2 (Fold Change, log2FC) ≥ 2 and Padj < 0.05, including 845 upregulated genes and 617 down-regulated genes (Figure 4C, 4D and Supplementary Table 2). Good distinction of DEGs between different groups was presented in the hierarchical clustering heatmap (Figure 4E).

DEGs and functional enrichment analysis associated with atherosclerosis

In order to understand whether cGAS inhibition changes expression of the genes involved in modulating atherosclerosis, we identified DEGs in cGAS inhibition probably linked to atherosclerosis according to research published in the PubMed database. A total of 275 DEGs associated with atherosclerosis were found when cGAS was inhibited (Supplementary Table 3). These DEGs included colony-stimulating factors (Csfs such as Csf2, Csf2rb, Csf1r, and Lif.), interleukins (Ils, including Il1b ,Il7, Il18, Il10.), tumor necrosis factor family members (such as TNF superfamily members–Tnfsf–Tnfsf10, Tnfsf8, and Tnfsf11a.), toll-like receptors (Tlrs including Tlr3, Tlr5, and Tlr7), interferon regulatory factors (Irfs, Irf5, Irf7, and Irf9), genes related to cholesterol metabolism (Msr1, Acat2, Apoa1, Abca1 and Abcg5), signal transducers and activators of transcription (Stats, such as Stat1, and Stat2), IFN and IFN-stimulated genes (Isgs, Ifna, Ifnb, Ifit1, Ifit2, and Isg15.) and so on. In addition, expression of Il6, Tlr9, Irf3, Cd36, Acat1 and Sting were also decreased despite absolute log2FC < 2.

Gene Ontology (GO) analysis showed that the main functional enrichments occurred in biological processes (BP), including inflammatory response, immune response, cholesterol biosynthetic process and so on, which are of great importance in the formation and progression of atherosclerotic plaque. In addition, GO terms linked to molecular function (MF) and cellular compartments (CC) were also discovered, such as extracellular space, cell surface, endoplasmic reticulum, and cytokine activity (Figure 5A, 5B and Table 3). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed comparable results, with the most significant KEGG terms associated with atherosclerosis being cytokine–cytokine receptor interaction, Jak–STAT signaling pathway, steroid biosynthesis, chemokine signaling pathway, TNF signaling pathway, fat digestion and absorption, PPAR signaling pathway and TLR signaling pathway (Figure 5C, 5D and Table 4).

Table 3. GO enrichment analysis.

CategoryTermP valueFold enrichment
Biologyical processGO:0006954~inflammatory response0.0006.942
Biologyical processGO:0002376~immune system process0.0006.414
Biologyical processGO:0006955~immune response0.0007.275
Biologyical processGO:0045087~innate immune response0.0005.629
Biologyical processGO:0006695~cholesterol biosynthetic process0.00027.720
Biologyical processGO:0009615~response to virus0.00012.185
Biologyical processGO:0048661~positive regulation of smooth muscle cell proliferation0.00011.941
Biologyical processGO:0008203~cholesterol metabolic process0.00010.733
Biologyical processGO:0007584~response to nutrient0.00012.151
Biologyical processGO:0051607~defense response to virus0.0006.946
Celluar componentGO:0005615~extracellular space0.0003.196
Celluar componentGO:0005576~extracellular region0.0002.576
Celluar componentGO:0009986~cell surface0.0003.705
Celluar componentGO:0009897~external side of plasma membrane0.0004.824
Celluar componentGO:0016020~membrane0.0001.446
Celluar componentGO:0005783~endoplasmic reticulum0.0002.147
Celluar componentGO:0043231~intracellular membrane-bounded organelle0.0002.618
Celluar componentGO:0045121~membrane raft0.0004.169
Celluar componentGO:0005887~integral component of plasma membrane0.0002.134
Celluar componentGO:0005578~proteinaceous extracellular matrix0.0003.457
Molecular functionGO:0005125~cytokine activity0.0007.129
Molecular functionGO:0005515~protein binding0.0001.572
Molecular functionGO:0004713~protein tyrosine kinase activity0.0006.030
Molecular functionGO:0016491~oxidoreductase activity0.0002.746
Molecular functionGO:0004714~transmembrane receptor protein tyrosine kinase activity0.0009.827
Molecular functionGO:0005102~receptor binding0.0003.220
Molecular functionGO:0042803~protein homodimerization activity0.0002.411
Molecular functionGO:0003725~double-stranded RNA binding0.0007.581
Molecular functionGO:0005164~tumor necrosis factor receptor binding0.00012.438
Molecular functionGO:0004871~signal transducer activity0.0002.457
# 0.000 means < 0.0001.

Table 4. KEGG pathway analysis.

TermP valueFold enrichment
mmu04060: Cytokine-cytokine receptor interaction0.0005.189
mmu05164: Influenza A0.0004.360
mmu00100: Steroid biosynthesis0.00014.457
mmu05134: Legionellosis0.0006.884
mmu05321: Inflammatory bowel disease (IBD)0.0006.651
mmu04380: Osteoclast differentiation0.0004.360
mmu04668: TNF signaling pathway0.0004.680
mmu05162: Measles0.0004.039
mmu04630: Jak-STAT signaling pathway0.0003.789
mmu05160: Hepatitis C0.0003.751
mmu00900: Terpenoid backbone biosynthesis0.00010.236
mmu04975: Fat digestion and absorption0.0007.228
mmu04062: Chemokine signaling pathway0.0003.003
mmu05132: Salmonella infection0.0014.528
mmu05168: Herpes simplex infection0.0012.830
mmu03320: PPAR signaling pathway0.0014.414
mmu04620: Toll-like receptor signaling pathway0.0013.885
mmu05140: Leishmaniasis0.0014.905
mmu01130: Biosynthesis of antibiotics0.0032.567
mmu05161: Hepatitis B0.0042.956
# 0.000 means <0.0001.
Functional enrichment analysis associated with atherosclerosis. (A, B) The GO function classification map and bubble chart of DEGs associated with atherosclerosis. (C, D) KEGG pathway classification map and bubble chart of DEGs. The X-axis represents the functional classification, and Y-axis represents the significant differences after the log 10 conversion or gene ratio. GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes.

Figure 5. Functional enrichment analysis associated with atherosclerosis. (A, B) The GO function classification map and bubble chart of DEGs associated with atherosclerosis. (C, D) KEGG pathway classification map and bubble chart of DEGs. The X-axis represents the functional classification, and Y-axis represents the significant differences after the log 10 conversion or gene ratio. GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes.

Moreover, Metascape analysis showed remarkable enrichment in regulation of the cytokine-mediated signaling pathway (Figure 6A–6C). There were 22 overlapped DEGs in both GO enriched terms “Inflammatory response” and “Immune system process”/ “Immune response”/ “Innate immune response” and seven overlapped DEGs in the KEGG analysis terms “Cytokine–cytokine receptor interaction”, “TNF signaling pathway”, “Jak–STAT signaling pathway”, “ Chemokine signaling pathway” and “Toll-like receptor signaling pathway”(Figure 6D, 6E and Tables 5, 6).

Table 5. GO enrichment analysis of overlapped DEGs.

GeneGO term
Olr1GO:0006954~inflammatory response
GO:0002376~immune system process
Cxcl1GO:0006954~inflammatory response
GO:0006955~immune response
Cxcl2GO:0006954~inflammatory response
GO:0006955~immune response
Il10GO:0006954~inflammatory response
GO:0006955~immune response
Cxcl3GO:0006954~inflammatory response
GO:0006955~immune response
Ccl7GO:0006954~inflammatory response
GO:0006955~immune response
Tnfrsf8GO:0006954~inflammatory response
GO:0002376~immune system process
GO:0006955~immune response
GO:0045087~innate immune response
CleceGO:0006954~inflammatory response
GO:0006955~immune response
Ccr1GO:0006954~inflammatory response
GO:0006955~immune response
Tnfsf4GO:0006954~inflammatory response
GO:0006955~immune response
Tnfrsf1bGO:0006954~inflammatory response
GO:0006955~immune response
Tlr7GO:0006954~inflammatory response
GO:0002376~immune system process
GO:0006955~immune response
GO:0045087~innate immune response
Nlrc4GO:0006954~inflammatory response
GO:0045087~innate immune response
Samhd1GO:0006954~inflammatory response
GO:0006955~immune response
Tnfrsf25GO:0006954~inflammatory response
GO:0006955~immune response
Il1bGO:0006954~inflammatory response
GO:0006955~immune response
Cd5lGO:0006954~inflammatory response
GO:0002376~immune system process
Il18GO:0006954~inflammatory response
GO:0006955~immune response
Il27GO:0006954~inflammatory response
GO:0002376~immune system process
GO:0045087~innate immune response
AxlGO:0006954~inflammatory response
GO:0006955~immune response
GO:0045087~innate immune response
Tlr5GO:0006954~inflammatory response
GO:0006955~immune response
GO:0045087~innate immune response
Tlr3GO:0006954~inflammatory response
GO:0006955~immune response
GO:0045087~innate immune response

Table 6. KEGG pathway analysis of overlapped DEGs.

GeneKEGG pathway
Cxcl1mmu04060: Cytokine-cytokine receptor interaction
mmu04668: TNF signaling pathway
mmu04062: Chemokine signaling pathway
Cxcl2mmu04060: Cytokine-cytokine receptor interaction
mmu04668: TNF signaling pathway
mmu04062: Chemokine signaling pathway
Lifmmu04060: Cytokine-cytokine receptor interaction
mmu04668: TNF signaling pathway
mmu04630: Jak-STAT signaling pathway
Il1bmmu04060: Cytokine-cytokine receptor interaction
mmu04668: TNF signaling pathway
mmu04620: Toll-like receptor signaling pathway
Csf2mmu04060: Cytokine-cytokine receptor interaction
mmu04668: TNF signaling pathway
Stat1mmu04630: Jak-STAT signaling pathway
mmu04062: Chemokine signaling pathway
mmu04620: Toll-like receptor signaling pathway
Ifnb1mmu04060: Cytokine-cytokine receptor interaction
mmu04630: Jak-STAT signaling pathway
mmu04620: Toll-like receptor signaling pathway
Enrichment analysis of Metascape and overlapped DEGs. (A) Bar graph of enriched terms, colored according to P value. Networks of enriched terms: (B) colored according to cluster ID and (C) colored according to P value. (D, E) Chord plot of the relationship between the overlapped DEGs and their corresponding GO and KEGG terms, together with the log2FC of the genes. Left half of Chord indicates whether the gene expression was up- or down-regulated. The right half represents different terms with different colors.

Figure 6. Enrichment analysis of Metascape and overlapped DEGs. (A) Bar graph of enriched terms, colored according to P value. Networks of enriched terms: (B) colored according to cluster ID and (C) colored according to P value. (D, E) Chord plot of the relationship between the overlapped DEGs and their corresponding GO and KEGG terms, together with the log2FC of the genes. Left half of Chord indicates whether the gene expression was up- or down-regulated. The right half represents different terms with different colors.

Protein–protein interaction network of DEGs associated with atherosclerosis

Next, the protein–protein interaction (PPI) network of DEGs was built as shown below (Figure 7A). The hub genes identified by the MCC method were identified as follows: Stat1, Irf7, Ifih1, Isg15, Ifit1, Rsad2, Stat2, Irf9, Ifi44, and Dhx58, exhibiting lower expression in cGAS inhibition (Figure 7B). A functional subnet module containing hub genes was selected from the PPI network and functional enrichment analysis was performed using ClueGO. The nodes in this module related to atherosclerosis were mainly enriched in the “response to interferon” (Figure 7C–7E).

Results of analysis of the PPI network associated with atherosclerosis. (A) The PPI network of DEGs associated with atherosclerosis. The green circles represent down-regulated DEG-encoded proteins and the red circles represent upregulated DEG-encoded proteins. The size of the circle indicates the abundance of gene expression. (B) The top 10 hub genes in the PPI network by MCC method. (C–E) A function subnet module analysis performed by ClueGO. Terms are presented in different colors. PPI, Protein–Protein Interaction.

Figure 7. Results of analysis of the PPI network associated with atherosclerosis. (A) The PPI network of DEGs associated with atherosclerosis. The green circles represent down-regulated DEG-encoded proteins and the red circles represent upregulated DEG-encoded proteins. The size of the circle indicates the abundance of gene expression. (B) The top 10 hub genes in the PPI network by MCC method. (C–E) A function subnet module analysis performed by ClueGO. Terms are presented in different colors. PPI, Protein–Protein Interaction.

Verification of the genes expressed in synergistic inflammatory signaling

We sought to verify and extend the main findings suggested by RNA-seq. RAW264.7 cells were treated with or without the cGAS inhibitor RU.521(2 μg/mL) for 12 h and then real-time quantitative polymerase chain reaction (RT-qPCR) was performed. As shown in Figure 8A, a robust anti-inflammatory effect can be induced by cGAS inhibition via increasing the expression of anti-inflammatory factors (Il10, Arg1, Lif, Il1rn) and reducing inflammatory factors (Il1b, Il7, Tnfsf10, Il18). Additionally, Il6 expression was also significantly decreased despite absolute log2FC < 2 in the RNA-seq results. In general, macrophages are classified as classical M1 or alternative M2. M1 macrophages are pro-inflammatory with high levels of cytokines such as IL 1B, IL6, and TNFα. Meanwhile M2 macrophages show anti-inflammatory properties, secreting anti-inflammatory factors, such as IL10 and ARG1. When stimulated by TLR ligands and IFNs, macrophages can transform into the M1 phenotype [22]. Furthermore, it is widely known that abnormal TLR signaling contributes to chronic inflammation [23]. Our results demonstrated that cGAS inhibition reduced the mRNA expression of Tlrs and, Ifns as well as the IFN-responsive genes (Figure 8B, 8C). These indicated that cGAS inhibition can restrain M1 polarization, promote M2 polarization, and that TLR or IFN signaling mediated by cGAS may participate in this process. Evidence suggests that M1 polarization can also be regulated by STAT1, IRF3, and IRF5 [24, 25], and inflammatory genes can be regulated by STAT1–IFN–IRF through over-expression in their promoters [26]. Previous studies have reported that STAT and IRF in the atherosclerotic tissues are displayed at significantly higher levels when compared with the matched normal tissues [27–31]. Here, Stat (Stat1, Stat2) and Irf (Irf3, Irf5, Irf7, Irf9) mRNA were significantly down regulated under cGAS inhibition (Figure 8B), consistent with the down regulation of inflammatory factors. RAW264.7 macrophages were pretreated with RU.521 for 12 h and then incubated for another 6 h with lipopolysaccharide (LPS, 10 ng/mL) plus interferon gamma (IFNγ, 20 ng/mL) that polarizes macrophages to the M1 phenotype. The results showed that cGAS significantly inhibited the M1 phenotype by decreasing Il1b, Il6, Tnfa, and Cd86 while increasing Il10, Arg1, and Nos2 mRNA expression which are characteristic of the M2 phenotype. Consistently, cGAS inhibition also inhibited the up-regulation of Stat and Irf in response to LPS plus IFNγ stimulation (Figure 8F). These results indicated that STAT and IRF signaling mediated by cGAS may be the key point of macrophage polarization through inflammatory regulation.

Inflammatory signaling factors TLRs, STAT/IRF as well as IFNs were verified by RT-qPCR. (A–E) RAW264.7 cells were treated for 12 h with RU.521 (2 μg/mL) or with DMSO as vehicle control. (F) RAW264.7 macrophages were pretreated with RU.521 for 12 h and then incubated for another 6 h with LPS (10 ng/mL) plus IFNγ (20 ng/mL) that polarizes macrophages to M1. Relative expression by RT-qPCR analysis. * P P P P

Figure 8. Inflammatory signaling factors TLRs, STAT/IRF as well as IFNs were verified by RT-qPCR. (A–E) RAW264.7 cells were treated for 12 h with RU.521 (2 μg/mL) or with DMSO as vehicle control. (F) RAW264.7 macrophages were pretreated with RU.521 for 12 h and then incubated for another 6 h with LPS (10 ng/mL) plus IFNγ (20 ng/mL) that polarizes macrophages to M1. Relative expression by RT-qPCR analysis. * P < 0.05, ** P < 0.01, *** P < 0.001, **** P < 0.0001, compared with control. Data are means ± SD. LPS, Lipopolysaccharide; IFNγ, interferon gamma; RT-qPCR, real time fluorescence quantitative PCR.

Given that foam cell formation acts as a hallmark of atherosclerosis, the expression of genes contributing to cholesterol uptake, cholesterol esterification and hydrolysis as well as cholesterol efflux was evaluated via RT-qPCR. cGAS inhibition blocked expression of Cd36, Msr1, Acat1, and Acat2, which mediate cholesterol uptake, cholesterol esterification and hydrolysis. In contrast, cGAS inhibition significantly increased the mRNA expression of Apoa1 and Abcg5, which contribute to cholesterol efflux, consistent with the RNA-Seq results. However, cGAS failed to increase the expression of Abca1 that promotes cholesterol efflux (Figure 8D). These results indicated that cGAS inhibition could suppress foam cell formation by hindering cholesterol uptake. Obviously, this assumption was confirmed in the RAW264.7 in vitro model (Figure 3C, 3D). In addition, the mRNA expression of colony-stimulating factors as well as chemokines was also verified (Figure 8E).

Discussion

As a chronic inflammatory disease, atherosclerosis shows the following characteristics: recruitment of circulating monocytes and their migration into the vascular intima, lipid accumulation, vascular local inflammation, smooth muscle cell (SMCs) proliferation, apoptosis, necrosis, and fibrosis [32–34]. Even in asymptomatic atherosclerosis, patients can suffer from coronary and carotid artery disease [35]. Therefore, identifying the molecular targets of atherosclerosis is imperative.

Following analysis of the microarray data of GSE40156 and GSE57614, cGAS was found to be increased in the total aorta of ApoE -/- mice and M1 macrophages. Consistently, our experiment first demonstrated that cGAS was expressed in the plaque of ApoE -/- mice regardless of diet, and was positively correlated with plaque area, mainly distributed in macrophages. Intriguingly, cGAS can also be found located in DAPI negative regions, and western blot results found that cGAS was clearly detected in the culture supernatant of RAW264.7 cells when treated with deoxyribonucleic acid sodium salt from herring testes (HT-DNA) or LPS plus IFNγ (Supplementary Figure 1), suggesting that cGAS may be secreted extracellularly. Further research is required to confirm this phenomenon and elaborate on how cGAS is secreted extracellularly. cGAS inhibition reduced foam cell formation in RAW264.7 macrophages and restrained M1 polarization as well as expression of inflammatory factors. Accordingly, cGAS may act as a critical molecule for the development of atherosclerosis via an inflammatory response mediated by macrophages. In order to reveal the possible mechanisms, RNA-seq was conducted in this research. A total of 1,462 DEGs were identified, 275 of which could be associated with atherosclerosis. Simultaneously, functional enrichment analysis and PPI network construction were performed on the 275 DEGs contributing to atherosclerosis, and some intriguing results were found. The DEGs were primarily concentrated in immune response, inflammatory response, and cytokine–cytokine receptor interaction signaling pathways, which have been proven to be critical in the development of atherosclerosis. We also found that the overlapped genes were mainly concentrated in the GO terms of “Immune response, Inflammatory response”, and “Cytokine–cytokine receptor interaction” pathways. The research results herein indicate that cGAS is probably a novel clue for the development of atherosclerosis through the inflammatory signaling pathway.

It has been reported that atherosclerosis can be regulated by abnormal activation of the STAT and IRF signaling pathways [27, 28, 30, 36, 37]. It is universally acknowledged that inflammatory cytokines are mediators of atherosclerosis, and previous studies have demonstrated that inflammatory genes can be regulated by STAT1–NF-κB or STAT1–IRF through over-expression in their promoters [26]. Pro-inflammatory transcription can be enhanced by signal integration of IFN with TLR, which is involved in sequential recruitment of STAT1-complexes and NF-kB [38]. This suggests that type I IFN and TLR cross-talk could promote the process of atherogenesis. As the important hub genes, Stat and Irf mRNA expression were demonstrated to be reduced, accompanied by decreased expression of genes related to pro-inflammatory and immunoregulatory effects, and increased expression of genes related to anti-inflammatory effects. Meanwhile, cGAS inhibition also reduced expression of Tlrs and Ifns. Type I IFNs exert pro-inflammatory immune effects, subsequently promoting atherosclerosis [37, 39, 40]. Up-regulated ISGs are also detected in macrophages of mouse plaques, indicating a type-I IFN responsive subset [40]. Extensive study has shown that IRF3/7 activates the IFNα/β transcriptional promoter [40–42]. The cGAS–STING signaling pathway triggers the downstream TANK binding kinase (TBK)1, followed by phosphorylation of IRF3 as well as IRF7, which form homodimer and enter the nucleus from cytoplasm, accompany with or without other transcription factors, such as STAT and NF-kB, thereby allowing initiation of the subsequent production of CXCL10 and additional type I IFN [41, 43, 44]. On the other hand, the binding of type I IFN to the IFNα/β receptor (IFNAR)2 recruits IFNAR1. This complex enables activation of the receptor-associated JAK1 and tyrosine kinase (TYK)2, followed by STAT1 and STAT2 phosphorylation, which bind to IRF9, forming IFN-stimulated gene factor (ISGF)3 complex. The ISGF3 complex translocates into nucleus and promote the production of ISGs and IRF7 by binding to IFN-stimulated regulatory elements (ISRE) as well as IRF7 elements in DNA. Phosphorylated STAT1 can also form a homodimer, which binds to a comparable γ-activated sequence (GAS) in DNA, inducing the expression of IRF1 and pro-inflammatory genes [45, 46]. One study also found a positive correlation between lipid quantification by magnetic resonance imaging and the upregulation of genes of the IFN/STAT1 pathway [47]. TLRs, a crucial class of pattern recognition receptors contributing to inflammatory responses, can induce a subset of ISGs via activation of IRF [48]. Microbial membrane components such as lipoproteins can be recognized by cell surface TLRs. By contrast, the intracellular TLRs such as TLR3, TLR7, and TLR9 recognize pathogenic RNA or DNA [49]. In combination with our findings, these imply that cGAS inhibition can further downregulate IFNβ production and the subsequent inflammatory cascade by inhibiting TLRs. In addition, it has been reported that there is cross-talk between TLR and STAT signaling. STAT1 S727 phosphorylation can be induced by multiple TLRs dependent on MyD88/TRIF signaling instead of IRF and IFN signaling. Phosphorylated STAT1 transfers into the nucleus and augments TLR–NF-kB activation, promoting the expression of pro-inflammatory genes [50] (Figure 9).

Signal integrations of TLRs, STAT/IRF as well as type-I IFN exacerbate synergistic amplification of gene expression that leads to an inflammatory cascade and pro-atherogenic responses. cGAS is activated by dsDNA and triggers the downstream TANK binding kinase (TBK)1, followed by phosphorylation of IRF3 as well as IRF7, which form homodimer and enter the nucleus from cytoplasm, accompany with or without other transcription factors such as STAT and NF-kB, thereby allowing initiation of the subsequent production of type I IFNα/β. On the other hand, the binding of type I IFN to the IFNα/β receptor (IFNAR)2 recruits IFNAR1. This complex enables activation of the receptor-associated JAK1 and tyrosine kinase (TYK)2, followed by STAT1 and STAT2 phosphorylation, which bind to IRF9, forming IFN-stimulated gene factor (ISGF)3 complex. The ISGF3 complex translocates into nucleus and promote the production of ISGs and IRF7 by binding to IFN-stimulated regulatory elements (ISRE) as well as IRF7 elements in DNA. Phosphorylated STAT1 can also form a homodimer, which binds to a comparable γ-activated sequence (GAS) in DNA, inducing the expression of IRF1 and pro-inflammatory genes. Moreover, cGAS can also result in up-regulated expression of TLRs (TLR3, TLR7, TLR9) and STAT (STAT1, STAT2). STAT1 phosphorylation can be induced by multiple TLRs dependent on MyD88 and TRIF signaling. Phosphorylated STAT1 translocates into the nucleus and augments TLR-NF-kB activation, promoting the expression of pro-inflammatory genes. Multiple TLRs also activate IRF5 and IRF7 as well as IRF3 via MyD88 and TRIF signaling.

Figure 9. Signal integrations of TLRs, STAT/IRF as well as type-I IFN exacerbate synergistic amplification of gene expression that leads to an inflammatory cascade and pro-atherogenic responses. cGAS is activated by dsDNA and triggers the downstream TANK binding kinase (TBK)1, followed by phosphorylation of IRF3 as well as IRF7, which form homodimer and enter the nucleus from cytoplasm, accompany with or without other transcription factors such as STAT and NF-kB, thereby allowing initiation of the subsequent production of type I IFNα/β. On the other hand, the binding of type I IFN to the IFNα/β receptor (IFNAR)2 recruits IFNAR1. This complex enables activation of the receptor-associated JAK1 and tyrosine kinase (TYK)2, followed by STAT1 and STAT2 phosphorylation, which bind to IRF9, forming IFN-stimulated gene factor (ISGF)3 complex. The ISGF3 complex translocates into nucleus and promote the production of ISGs and IRF7 by binding to IFN-stimulated regulatory elements (ISRE) as well as IRF7 elements in DNA. Phosphorylated STAT1 can also form a homodimer, which binds to a comparable γ-activated sequence (GAS) in DNA, inducing the expression of IRF1 and pro-inflammatory genes. Moreover, cGAS can also result in up-regulated expression of TLRs (TLR3, TLR7, TLR9) and STAT (STAT1, STAT2). STAT1 phosphorylation can be induced by multiple TLRs dependent on MyD88 and TRIF signaling. Phosphorylated STAT1 translocates into the nucleus and augments TLR-NF-kB activation, promoting the expression of pro-inflammatory genes. Multiple TLRs also activate IRF5 and IRF7 as well as IRF3 via MyD88 and TRIF signaling.

Overall, signal integration of TLRs, STAT/IRF as well as IFN exacerbates synergistic amplification of gene expression, leading to inflammatory cascade and pro-atherogenic responses. It is well known that macrophages are of great importance in the initiation, progression and regression of atherosclerosis, and that different subtypes of macrophages affect the outcome of atherosclerosis. Classically, the pro-inflammatory M1 macrophage phenotype is responsible for atherosclerotic plaque vulnerability, whereas the anti-inflammatory M2 macrophage phenotype can increase atherosclerotic plaque stability [21]. Here, our result demonstrated that cGAS inhibition reduced expression of cytokines related to M1 polarization such as Il1b and Il6, while it increased expression of cytokines related to M2 polarization such as Arg1 and Il10. Accumulating evidence suggests that STAT1, IRF3, and IRF5 contribute to M1 polarization. Conversely, STAT3, STAT6, and PPARγ contribute to M2 polarization [24, 25]. Moreover, the conversion from M2 to M1 can be regulated by activated IRF7 signaling [51]. IL10 inhibits production of IL1B and IL6 and promotes STAT3 phosphorylation [52–55], which causes the macrophage phenotype to revert from M1 to M2. Of note, LIF, a pleiotropic cytokine, has been proven to be anti-inflammatory, and increases IL10 expression, which can ameliorate atherosclerosis [56, 57]. Data have demonstrated that LIF is effective in not only inhibiting the formation of plaque but also retarding the progression of pre-existing atherosclerotic plaques in atherosclerosis. atherosclerosis [58].

APOA1, the main protein of high-density lipoprotein (HDL) which shows protective effects on atherosclerosis [59], was increased when cGAS was inhibited. Raising APOA1/HDL levels dampens the inflammatory activities of leukocytes and macrophages [60]. A previous study demonstrated that after infection with Salmonella (a Gram-negative bacterium that expresses LPS), the plasma levels of IFNβ significantly increase in ApoA1 deficient mice compared with wild-type mice [61]. RNA-seq transcriptomic analysis also showed that GO analysis for the rHDL-repressed genes revealed the following categories: cytokine biosynthetic process, cellular response to IFNβ, defense response to virus, and regulation of IL1B production [59]. Likewise, cGAS inhibition reduced Cd36 and Msr1, by which the modified LDL is taken up. A previous study suggested STAT1 acetylation and its interaction with PPARγ induced by cholesterol crystals increase CD36 expression, which facilitates oxLDL uptake and foam cell formation [62]. Similarly, MSR1 upregulation is dependent on activation of the canonical Jak-STAT signaling pathway [63]. Furthermore, Acat1 and Acat2, which mediate cholesterol esterification, were reduced when cGAS was inhibited, and their expression as well as overall activity in macrophages can be induced by IFNβ stimulation [64].

In conclusion, cGAS may contribute to the development of atherosclerosis. Gene expression profiling of cGAS inhibition in macrophages successfully identified important biological processes and pathways specific to atherosclerosis. The cGAS exacerbates the inflammatory cascade through synergistic signaling of TLRs, STAT, and IRF as well as IFN, triggering modulation of macrophage phenotypes to M1, thus increasing lipid deposition by up-regulating molecules related to cholesterol uptake. As the results are limited to sequencing analysis, a large number of experiments are still warranted. Despite this, these findings offer the theoretical foundation for future research into the underlying mechanisms of atherosclerosis and may provide potential therapeutic targets.

Materials and Methods

Animals

Eight-week-old ApoE -/- male mice on a C57BL/6 background from Hua Fukang Experimental Animal Center (Beijing, China) were kept under standard conditions with constant humidity (approximately 40%) and temperature (approximately 22° C) as well as a 12-h light/dark cycle. A standard chow diet was offered for 1 week prior to all procedures, after which a western-type diet containing 0.21% cholesterol and 41 kCal% fat (Special Diet Services, Huafukang, China) or continuation on the standard chow diet was offered for 16weeks. Hearts and aortas were harvested for assessment of atherosclerotic burden.

Cell culture and treatment

RAW 264.7 macrophages purchased from the Cell Resource Center of Shanghai Institute for Biological Sciences were grown in Dulbecco’s Modified Eagle’s Medium (DMEM, HyClone, Logan, UT, USA) containing 10% heat-inactivated fetal bovine serum (FBS, ScienCell Research Laboratories, Carlsbad, CA, USA), 1%penicillin and 1% streptomycin (HyClone). Cells were incubated under standard conditions with constant temperature (37° C) and an atmosphere of 5% CO2. Cells were grown to 60–80% confluence before use. In order to perform RNA-seq, cells were treated for 12 h with RU.521 (cGAS inhibitor 2 μg/mL, InvivoGen, San Diego, CA, USA) or with DMSO as vehicle control. In addition, RAW264.7 macrophages were pretreated with RU.521 (2 μg/mL) for 12 h and then incubated for another 6 h with LPS (10 ng/mL, HY-D1056, MedChemExpress, Monmouth Junction, NJ, USA) plus IFNγ (20 ng/mL, HY-P7071, MedChemExpress) that polarize macrophages to M1.

Histologic and morphometric analysis

As previously described, mouse hearts were embedded in OCT and frozen at −80° C after perfusion, fixation, and dehydration. Then serial sections were cut at 10μm thickness using a cryostat (CryoStar NX50, Thermo Scientific, Waltham, USA). Hematoxylin and eosin (HE) staining of sections was performed for quantification of lesion areas.

Immunofluorescent staining

Immunofluorescent staining was performed as previously described [65]. The primary antibodies used were as follows: cGAS (Mouse, 1:100, sc-515777, Santa Cruz Biotechnology, Santa Cruz, CA, USA), MOMA-2 (Rat, 1:100, ab33451, Abcam, Cambridge, UK), CD31 (1:100, AF3628, R&D Systems, Minneapolis, MN, USA) and α-SMA (Rabbit, 1:100, GB13044, Servicebio, Wuhan, China). 4′,6-diamidino-2-phenylindole (DAPI, Invitrogen, Carlsbad, CA, USA) was applied to stain nuclei. Fluorescent images were captured with a Nikon A1Si confocal microscope (Nikon, Tokyo, Japan), and analyzed with NIS Elements AR Imaging Software 4.10 (Nikon) and ImageJ 1.41 software (NIH, Bethesda, MD, USA).

Dil-oxLDL uptake assay

Cells were treated with or without RU.521(2 μg/mL) for 24 h, followed by incubation with 20 μg/mL Dil-LDL (YB-0010, Yiyuan Corporation, Guangzhou, China) for 6 h or 8 h at 37° C. Then the cells were fixed and analyzed by immunofluorescent staining.

Cholesterol efflux

Cholesterol efflux assays were performed as previously described [66]. Briefly, RAW264.7 macrophages were labeled with NBD (N1148, Life Technologies, Carlsbad, CA, USA) in the presence of RU.521 (2 μg/mL). After the cholesterol pools were equilibrated, cells were incubated with HDL (50 μg/mL, YB-003, Yiyuan Corporation) for 6 h. The fluorescence intensity of NBD was measured at 485nm excitation and 535nm emission. Efflux was measured as a percentage of fluorescence intensity in medium/(fluorescence intensity in medium + fluorescence intensity in cells) ×100%.

Foam cell formation assays

RAW264.7 macrophages were incubated in DMEM containing oxLDL (100 μg/mL, YB-002, Yiyuan Corporation) with or without RU.521 (2 μg/mL) for 24 h. After fixation for 20 min, cells were blocked with 1% BSA, followed by staining with BODIPY 493/503 (D3922, Invitrogen) for 30 min. Next, nuclei visualized by staining DAPI for 10 min. Finally, antifade mounting medium (P0126, Beyotime Institute of Biotechnology, Jiangsu, China) was applied to prevent fluorescence quenching. Images were acquired using a Nikon A1Si confocal microscope. Analysis of different images was performed using ImageJ and Adobe Photoshop CC. Oil Red O staining was carried out as previously described [66, 67]. RAW264.7 macrophages were incubated in DMEM containing oxLDL (100 μg/mL) with or without RU.521 (2 μg/mL) for 24 h. After fixation for 30 min, cells were stained with oil red O (G1015, Servicebio) for 30 min, and then visualized by light microscopy (Olympus, Tokyo, Japan).

Extraction and measurement of plasma dsDNA

Individuals (healthy volunteers and patients with atherosclerosis) were enrolled from Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China. Atherosclerosis was confirmed through ultrasonography or MRA/CTA/DSA imaging of the carotid or cerebral artery. Individuals with autoimmune diseases were excluded. Blood was collected and the plasma was isolated. Extracellular dsDNA in plasma was extracted using a DNA Extractor SP Kit (ZWK-296-60501, Wako Pure Chemical Industries Ltd., Osaka, Japan) and quantified using the Quant-iT™ PicoGreen ® dsDNA Reagent Kit (P7589, Invitrogen) following the manufacturer’s instructions.

RNA extraction and real-time quantitative PCR analysis

Total cellular RNA was isolated and extracted using the TRIzol reagent (Vazyme, Nanjing, China). Quantitative RT-qPCR was carried out using the SYBR Green detection chemistry (Vazyme, Nanjing, China) in a 10μL reaction volume with an ABI 7500 Fast real-time PCR system. Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) expression was used as the internal normalized reference gene. All samples were measured in triplicate, and the ΔΔCt method was applied to relative quantitative measurements. Primers are shown in Supplementary Table 1.

Western blotting analysis

Western blotting analysis was conducted as previously described [65]. Primary antibodies against cGAS (Mouse, 1:1000, sc-515777, Santa Cruz Biotechnology) and, GAPDH (GAPDH, A01020, Abbkine, Redlands, CA, USA) were used. ImageJ software (NIH) was used to analyze protein expression levels which were normalized to GAPDH.

Microarray data and data processing of DEGs

The microarray data of GSE40156 and GSE57614 derived from Gene Expression Omnibus (GEO) databases(https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE40156 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE57614) which were based on the GPL1261 [Mouse430_2] Affymetrix Mouse Genome 430 2.0 Array and GPL6480 Agilent-014850 Whole Human Genome Microarray 4x44K G4112F. DEGs were picked out using GEO2R online tools. A P value <0.05 and | log2FC | >1 were chosen as cut-off standards.

RNA sequencing

RNA sequencing was performed following a previously published method [67]. Briefly, for each sample, a total amount of 1μg RNA was used as input material. The NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (New England Biolabs (NEB), Ipswich, MA, USA) was applied to prepare sequencing libraries, followed by generation of clusters with the TruSeq PE Cluster Kit v3-cBot-HS (Illumina Inc., San Diego, CA, USA). The Illumina Novaseq platform was utilized to sequence libraries and then about 150 bp paired-end reads were obtained. Analysis of differential expression was performed using the edgeR R package. The Benjamini and Hochberg method was used to adjust P value. A Padj < 0.05 and | log2FC | ≥ 2 were set as the thresholds for significant differential expression.

Functional enrichment analysis of DEGs

GO enrichment analysis and KEGG pathway analysis of DEGs was performed using the online tool DAVID (version: 6.8, https://david.ncifcrf.gov/) with thresholds of count ≥ 2 and P value < 0.05. Furthermore, Metascape (https://metascape.org/gp/index.html#/main/step1) was also used to obtain functional enrichment results.

Protein–protein interaction network analysis of DEGs

The PPI network analysis of DEGs was based on the STRING database (https://string-db.org/), which functions as a tool for predicting protein–protein interactions. A threshold of PPI score (medium confidence) ≥ 0.4 was used to predict the DEG-encoded proteins. Then, Cytoscape software (version: 3.6.0) was used to visualize the PPI network. Module and hub gene analyses were carried out with MCODE, cytoHubba plugin as well as ClueGO of Cytoscape software.

Statistical analyses

Results are expressed as mean ± standard deviation (SD). The analysis was completed using GraphPad Prism 6. In addition, the Student's t-test was used for comparison between the two groups. P values < 0.05 or Padj < 0.05 were considered statistically significant. R (3.6.2) software was used to display the distribution of DEGs as well as functional enrichment analysis.

Ethics statement

This study was approved by the Medical Ethics Committee of Tongji Medical College and the Institutional Animal Care and Use Committee, Huazhong University of Science and Technology (HUST), Wuhan, China (S947, S2342). Written informed consents were obtained from Individuals (healthy volunteers and patients with atherosclerosis).

Author Contributions

H.B. and C.F. conceived and designed the experiments. L.G.F and C.S.C performed the experiments, downloaded and analyzed the data. L.G.F., C.S.C. and. X.Y.P. wrote the manuscript and performed the visualization. Y.Z.M helped to revise the manuscript. H.B. and C.F. were responsible for project management and funding acquisition, and all authors read and approved the manuscript.

Acknowledgments

We gratefully acknowledge discussion and suggestions by Shuai Yang and Xingjie Hao that helped us to improve the paper.

Conflicts of Interest

The authors deny any conflicts of interest.

Funding

This work was supported by the National Natural Science Foundation of China (NSFC) (grant numbers 81801172) and the Temple Kallikon Research Foundation (grant numbers KF20120204).

References

  • 1. Tan X, Sun L, Chen J, Chen ZJ. Detection of microbial infections through innate immune sensing of nucleic acids. Annu Rev Microbiol. 2018; 72:447–78. https://doi.org/10.1146/annurev-micro-102215-095605 [PubMed]
  • 2. Sun L, Wu J, Du F, Chen X, Chen ZJ. Cyclic GMP-AMP synthase is a cytosolic DNA sensor that activates the type I interferon pathway. Science. 2013; 339:786–91. https://doi.org/10.1126/science.1232458 [PubMed]
  • 3. Kerur N, Fukuda S, Banerjee D, Kim Y, Fu D, Apicella I, Varshney A, Yasuma R, Fowler BJ, Baghdasaryan E, Marion KM, Huang X, Yasuma T, et al. cGAS drives noncanonical-inflammasome activation in age-related macular degeneration. Nat Med. 2018; 24:50–61. https://doi.org/10.1038/nm.4450 [PubMed]
  • 4. Ghosh A, Shao L, Sampath P, Zhao B, Patel NV, Zhu J, Behl B, Parise RA, Beumer JH, O’Sullivan RJ, DeLuca NA, Thorne SH, Rathinam VA, et al. Oligoadenylate-synthetase-family protein OASL inhibits activity of the DNA sensor cGAS during DNA virus infection to limit interferon production. Immunity. 2019; 50:51–63.e5. https://doi.org/10.1016/j.immuni.2018.12.013 [PubMed]
  • 5. Cohen D, Melamed S, Millman A, Shulman G, Oppenheimer-Shaanan Y, Kacen A, Doron S, Amitai G, Sorek R. Cyclic GMP-AMP signalling protects bacteria against viral infection. Nature. 2019; 574:691–95. https://doi.org/10.1038/s41586-019-1605-5 [PubMed]
  • 6. Benmerzoug S, Ryffel B, Togbe D, Quesniaux VF. self-DNA sensing in lung inflammatory diseases. Trends Immunol. 2019; 40:719–34. https://doi.org/10.1016/j.it.2019.06.001 [PubMed]
  • 7. Ablasser A, Chen ZJ. cGAS in action: expanding roles in immunity and inflammation. Science. 2019; 363:eaat8657. https://doi.org/10.1126/science.aat8657 [PubMed]
  • 8. Chung KW, Dhillon P, Huang S, Sheng X, Shrestha R, Qiu C, Kaufman BA, Park J, Pei L, Baur J, Palmer M, Susztak K. Mitochondrial damage and activation of the STING pathway lead to renal inflammation and fibrosis. Cell Metab. 2019; 30:784–99.e5. https://doi.org/10.1016/j.cmet.2019.08.003 [PubMed]
  • 9. Kato Y, Park J, Takamatsu H, Konaka H, Aoki W, Aburaya S, Ueda M, Nishide M, Koyama S, Hayama Y, Kinehara Y, Hirano T, Shima Y, et al. Apoptosis-derived membrane vesicles drive the cGAS-STING pathway and enhance type I IFN production in systemic lupus erythematosus. Ann Rheum Dis. 2018; 77:1507–15. https://doi.org/10.1136/annrheumdis-2018-212988 [PubMed]
  • 10. Libby P, Buring JE, Badimon L, Hansson GK, Deanfield J, Bittencourt MS, Tokgözoğlu L, Lewis EF. Atherosclerosis. Nat Rev Dis Primers. 2019; 5:56. https://doi.org/10.1038/s41572-019-0106-z [PubMed]
  • 11. Jylhävä J, Lehtimäki T, Jula A, Moilanen L, Kesäniemi YA, Nieminen MS, Kähönen M, Hurme M. Circulating cell-free DNA is associated with cardiometabolic risk factors: the health 2000 survey. Atherosclerosis. 2014; 233:268–71. https://doi.org/10.1016/j.atherosclerosis.2013.12.022 [PubMed]
  • 12. West AP, Shadel GS. Mitochondrial DNA in innate immune responses and inflammatory pathology. Nat Rev Immunol. 2017; 17:363–75. https://doi.org/10.1038/nri.2017.21 [PubMed]
  • 13. Zhang Z, Meng P, Han Y, Shen C, Li B, Hakim MA, Zhang X, Lu Q, Rong M, Lai R. Mitochondrial DNA-LL-37 complex promotes atherosclerosis by escaping from autophagic recognition. Immunity. 2015; 43:1137–47. https://doi.org/10.1016/j.immuni.2015.10.018 [PubMed]
  • 14. Jin H, Ko YS, Park SW, Kim HJ. P2Y2R activation by ATP induces oxLDL-mediated inflammasome activation through modulation of mitochondrial damage in human endothelial cells. Free Radic Biol Med. 2019; 136:109–17. https://doi.org/10.1016/j.freeradbiomed.2019.04.004 [PubMed]
  • 15. Yu E, Calvert PA, Mercer JR, Harrison J, Baker L, Figg NL, Kumar S, Wang JC, Hurst LA, Obaid DR, Logan A, West NE, Clarke MC, et al. Mitochondrial DNA damage can promote atherosclerosis independently of reactive oxygen species through effects on smooth muscle cells and monocytes and correlates with higher-risk plaques in humans. Circulation. 2013; 128:702–12. https://doi.org/10.1161/CIRCULATIONAHA.113.002271 [PubMed]
  • 16. Lood C, Blanco LP, Purmalek MM, Carmona-Rivera C, De Ravin SS, Smith CK, Malech HL, Ledbetter JA, Elkon KB, Kaplan MJ. Neutrophil extracellular traps enriched in oxidized mitochondrial DNA are interferogenic and contribute to lupus-like disease. Nat Med. 2016; 22:146–53. https://doi.org/10.1038/nm.4027 [PubMed]
  • 17. Bai J, Cervantes C, Liu J, He S, Zhou H, Zhang B, Cai H, Yin D, Hu D, Li Z, Chen H, Gao X, Wang F, et al. DsbA-L prevents obesity-induced inflammation and insulin resistance by suppressing the mtDNA release-activated cGAS-cGAMP-STING pathway. Proc Natl Acad Sci USA. 2017; 114:12196–201. https://doi.org/10.1073/pnas.1708744114 [PubMed]
  • 18. Martinet W, Knaapen MW, De Meyer GR, Herman AG, Kockx MM. Elevated levels of oxidative DNA damage and DNA repair enzymes in human atherosclerotic plaques. Circulation. 2002; 106:927–32. https://doi.org/10.1161/01.cir.0000026393.47805.21 [PubMed]
  • 19. Tumurkhuu G, Shimada K, Dagvadorj J, Crother TR, Zhang W, Luthringer D, Gottlieb RA, Chen S, Arditi M. Ogg1-dependent DNA repair regulates NLRP3 inflammasome and prevents atherosclerosis. Circ Res. 2016; 119:e76–90. https://doi.org/10.1161/CIRCRESAHA.116.308362 [PubMed]
  • 20. Gehrke N, Mertens C, Zillinger T, Wenzel J, Bald T, Zahn S, Tüting T, Hartmann G, Barchet W. Oxidative damage of DNA confers resistance to cytosolic nuclease TREX1 degradation and potentiates STING-dependent immune sensing. Immunity. 2013; 39:482–95. https://doi.org/10.1016/j.immuni.2013.08.004 [PubMed]
  • 21. Yang S, Yuan HQ, Hao YM, Ren Z, Qu SL, Liu LS, Wei DH, Tang ZH, Zhang JF, Jiang ZS. Macrophage polarization in atherosclerosis. Clin Chim Acta. 2020; 501:142–46. https://doi.org/10.1016/j.cca.2019.10.034 [PubMed]
  • 22. Barrett TJ. Macrophages in Atherosclerosis Regression. Arterioscler Thromb Vasc Biol. 2020; 40:20–33. https://doi.org/10.1161/ATVBAHA.119.312802 [PubMed]
  • 23. Wang Y, Zhang S, Li H, Wang H, Zhang T, Hutchinson MR, Yin H, Wang X. Small-molecule modulators of toll-like receptors. Acc Chem Res. 2020; 53:1046–55. https://doi.org/10.1021/acs.accounts.9b00631 [PubMed]
  • 24. Vinué Á, Navarro J, Herrero-Cervera A, García-Cubas M, Andrés-Blasco I, Martínez-Hervás S, Real JT, Ascaso JF, González-Navarro H. The GLP-1 analogue lixisenatide decreases atherosclerosis in insulin-resistant mice by modulating macrophage phenotype. Diabetologia. 2017; 60:1801–12. https://doi.org/10.1007/s00125-017-4330-3 [PubMed]
  • 25. Zhao Y, Hao C, Bao L, Wang D, Li Y, Qu Y, Ding M, Zhao A, Yao W. Silica particles disorganize the polarization of pulmonary macrophages in mice. Ecotoxicol Environ Saf. 2020; 193:110364. https://doi.org/10.1016/j.ecoenv.2020.110364 [PubMed]
  • 26. Sikorski K, Wesoly J, Bluyssen HA. Data mining of atherosclerotic plaque transcriptomes predicts STAT1-dependent inflammatory signal integration in vascular disease. Int J Mol Sci. 2014; 15:14313–31. https://doi.org/10.3390/ijms150814313 [PubMed]
  • 27. Szelag M, Piaszyk-Borychowska A, Plens-Galaska M, Wesoly J, Bluyssen HA. Targeted inhibition of STATs and IRFs as a potential treatment strategy in cardiovascular disease. Oncotarget. 2016; 7:48788–812. https://doi.org/10.18632/oncotarget.9195 [PubMed]
  • 28. Sikorski K, Czerwoniec A, Bujnicki JM, Wesoly J, Bluyssen HA. STAT1 as a novel therapeutical target in pro-atherogenic signal integration of IFNγ, TLR4 and IL-6 in vascular disease. Cytokine Growth Factor Rev. 2011; 22:211–19. https://doi.org/10.1016/j.cytogfr.2011.06.003 [PubMed]
  • 29. Wang F, Liu Z, Park SH, Gwag T, Lu W, Ma M, Sui Y, Zhou C. Myeloid β-catenin deficiency exacerbates atherosclerosis in low-density lipoprotein receptor-deficient mice. Arterioscler Thromb Vasc Biol. 2018; 38:1468–78. https://doi.org/10.1161/ATVBAHA.118.311059 [PubMed]
  • 30. Recio C, Oguiza A, Mallavia B, Lazaro I, Ortiz-Muñoz G, Lopez-Franco O, Egido J, Gomez-Guerrero C. Gene delivery of suppressors of cytokine signaling (SOCS) inhibits inflammation and atherosclerosis development in mice. Basic Res Cardiol. 2015; 110:8. https://doi.org/10.1007/s00395-014-0458-1 [PubMed]
  • 31. Seneviratne AN, Edsfeldt A, Cole JE, Kassiteridi C, Swart M, Park I, Green P, Khoyratty T, Saliba D, Goddard ME, Sansom SN, Goncalves I, Krams R, et al. Interferon regulatory factor 5 controls necrotic core formation in atherosclerotic lesions by impairing efferocytosis. Circulation. 2017; 136:1140–54. https://doi.org/10.1161/CIRCULATIONAHA.117.027844 [PubMed]
  • 32. Kobiyama K, Ley K. Atherosclerosis. Circ Res. 2018; 123:1118–20. https://doi.org/10.1161/CIRCRESAHA.118.313816 [PubMed]
  • 33. Libby P, Ridker PM, Hansson GK. Progress and challenges in translating the biology of atherosclerosis. Nature. 2011; 473:317–25. https://doi.org/10.1038/nature10146 [PubMed]
  • 34. Moore KJ, Sheedy FJ, Fisher EA. Macrophages in atherosclerosis: a dynamic balance. Nat Rev Immunol. 2013; 13:709–21. https://doi.org/10.1038/nri3520 [PubMed]
  • 35. López-Melgar B, Fernández-Friera L, Oliva B, García-Ruiz JM, Sánchez-Cabo F, Bueno H, Mendiguren JM, Lara-Pezzi E, Andrés V, Ibáñez B, Fernández-Ortiz A, Sanz J, Fuster V. Short-term progression of multiterritorial subclinical atherosclerosis. J Am Coll Cardiol. 2020; 75:1617–27. https://doi.org/10.1016/j.jacc.2020.02.026 [PubMed]
  • 36. Liu H, Cheng WL, Jiang X, Wang PX, Fang C, Zhu XY, Huang Z, She ZG, Li H. Ablation of interferon regulatory factor 3 protects against atherosclerosis in apolipoprotein E-deficient mice. Hypertension. 2017; 69:510–20. https://doi.org/10.1161/HYPERTENSIONAHA.116.08395 [PubMed]
  • 37. Goossens P, Gijbels MJ, Zernecke A, Eijgelaar W, Vergouwe MN, van der Made I, Vanderlocht J, Beckers L, Buurman WA, Daemen MJ, Kalinke U, Weber C, Lutgens E, de Winther MP. Myeloid type I interferon signaling promotes atherosclerosis by stimulating macrophage recruitment to lesions. Cell Metab. 2010; 12:142–53. https://doi.org/10.1016/j.cmet.2010.06.008 [PubMed]
  • 38. Piaszyk-Borychowska A, Széles L, Csermely A, Chiang HC, Wesoły J, Lee CK, Nagy L, Bluyssen HA. Signal integration of IFN-I and IFN-II with TLR4 involves sequential recruitment of STAT1-complexes and NFκB to enhance pro-inflammatory transcription. Front Immunol. 2019; 10:1253. https://doi.org/10.3389/fimmu.2019.01253 [PubMed]
  • 39. Levy Z, Rachmani R, Trestman S, Dvir A, Shaish A, Ravid M, Harats D. Low-dose interferon-alpha accelerates atherosclerosis in an LDL receptor-deficient mouse model. Eur J Intern Med. 2003; 14:479–83. https://doi.org/10.1016/j.ejim.2003.08.010 [PubMed]
  • 40. Chen HJ, Tas SW, de Winther MP. Type-I interferons in atherosclerosis. J Exp Med. 2020; 217:e20190459. https://doi.org/10.1084/jem.20190459 [PubMed]
  • 41. Watanabe T, Asano N, Fichtner-Feigl S, Gorelick PL, Tsuji Y, Matsumoto Y, Chiba T, Fuss IJ, Kitani A, Strober W. NOD1 contributes to mouse host defense against helicobacter pylori via induction of type I IFN and activation of the ISGF3 signaling pathway. J Clin Invest. 2010; 120:1645–62. https://doi.org/10.1172/JCI39481 [PubMed]
  • 42. Jefferies CA. Regulating IRFs in IFN driven disease. Front Immunol. 2019; 10:325. https://doi.org/10.3389/fimmu.2019.00325 [PubMed]
  • 43. Ning X, Wang Y, Jing M, Sha M, Lv M, Gao P, Zhang R, Huang X, Feng JM, Jiang Z. Apoptotic caspases suppress type I interferon production via the cleavage of cGAS, MAVS, and IRF3. Mol Cell. 2019; 74:19–31.e7. https://doi.org/10.1016/j.molcel.2019.02.013 [PubMed]
  • 44. Dai J, Huang YJ, He X, Zhao M, Wang X, Liu ZS, Xue W, Cai H, Zhan XY, Huang SY, He K, Wang H, Wang N, et al. Acetylation blocks cGAS activity and inhibits self-DNA-induced autoimmunity. Cell. 2019; 176:1447–60.e14. https://doi.org/10.1016/j.cell.2019.01.016 [PubMed]
  • 45. Mogensen TH. IRF and STAT transcription factors - from basic biology to roles in infection, protective immunity, and primary immunodeficiencies. Front Immunol. 2019; 9:3047. https://doi.org/10.3389/fimmu.2018.03047 [PubMed]
  • 46. Kim N, Kukkonen S, Martinez-Viedma Mdel P, Gupta S, Aldovini A. Tat engagement of p38 MAP kinase and IRF7 pathways leads to activation of interferon-stimulated genes in antigen-presenting cells. Blood. 2013; 121:4090–100. https://doi.org/10.1182/blood-2012-10-461566 [PubMed]
  • 47. Chai JT, Ruparelia N, Goel A, Kyriakou T, Biasiolli L, Edgar L, Handa A, Farrall M, Watkins H, Choudhury RP. Differential gene expression in macrophages from human atherosclerotic plaques shows convergence on pathways implicated by genome-wide association study risk variants. Arterioscler Thromb Vasc Biol. 2018; 38:2718–30. https://doi.org/10.1161/ATVBAHA.118.311209 [PubMed]
  • 48. Verweyen E, Holzinger D, Weinhage T, Hinze C, Wittkowski H, Pickkers P, Albeituni S, Verbist K, Nichols KE, Schulert G, Grom A, Foell D, Kessel C. Synergistic signaling of TLR and IFNα/β facilitates escape of IL-18 expression from endotoxin tolerance. Am J Respir Crit Care Med. 2020; 201:526–39. https://doi.org/10.1164/rccm.201903-0659OC [PubMed]
  • 49. Meng S, Chanda P, Thandavarayan RA, Cooke JP. Transflammation: innate immune signaling in nuclear reprogramming. Adv Drug Deliv Rev. 2017; 120:133–41. https://doi.org/10.1016/j.addr.2017.09.010 [PubMed]
  • 50. Luu K, Greenhill CJ, Majoros A, Decker T, Jenkins BJ, Mansell A. STAT1 plays a role in TLR signal transduction and inflammatory responses. Immunol Cell Biol. 2014; 92:761–69. https://doi.org/10.1038/icb.2014.51 [PubMed]
  • 51. Tanaka T, Murakami K, Bando Y, Yoshida S. Interferon regulatory factor 7 participates in the M1-like microglial polarization switch. Glia. 2015; 63:595–610. https://doi.org/10.1002/glia.22770 [PubMed]
  • 52. Sun Y, Ma J, Li D, Li P, Zhou X, Li Y, He Z, Qin L, Liang L, Luo X. Interleukin-10 inhibits interleukin-1β production and inflammasome activation of microglia in epileptic seizures. J Neuroinflammation. 2019; 16:66. https://doi.org/10.1186/s12974-019-1452-1 [PubMed]
  • 53. Rossato M, Curtale G, Tamassia N, Castellucci M, Mori L, Gasperini S, Mariotti B, De Luca M, Mirolo M, Cassatella MA, Locati M, Bazzoni F. IL-10-induced microRNA-187 negatively regulates TNF-α, IL-6, and IL-12p40 production in TLR4-stimulated monocytes. Proc Natl Acad Sci USA. 2012; 109:E3101–10. https://doi.org/10.1073/pnas.1209100109 [PubMed]
  • 54. Marshall JS, Leal-Berumen I, Nielsen L, Glibetic M, Jordana M. Interleukin (IL)-10 inhibits long-term IL-6 production but not preformed mediator release from rat peritoneal mast cells. J Clin Invest. 1996; 97:1122–28. https://doi.org/10.1172/JCI118506 [PubMed]
  • 55. Kondo H, Kira S, Oniki T, Gotoh K, Fukui A, Abe I, Ikebe Y, Kawano K, Saito S, Aoki K, Okada N, Nagano Y, Akioka H, et al. Interleukin-10 treatment attenuates sinus node dysfunction caused by streptozotocin-induced hyperglycaemia in mice. Cardiovasc Res. 2019; 115:57–70. https://doi.org/10.1093/cvr/cvy162 [PubMed]
  • 56. Nicola NA, Babon JJ. Leukemia inhibitory factor (LIF). Cytokine Growth Factor Rev. 2015; 26:533–44. https://doi.org/10.1016/j.cytogfr.2015.07.001 [PubMed]
  • 57. Xue B, Liu D, Song M, Zhao G, Cao Y, Yan G, Dai J, Hu Y. Leukemia inhibitory factor promotes the regeneration of rat uterine horns with full-thickness injury. Wound Repair Regen. 2019; 27:477–87. https://doi.org/10.1111/wrr.12729 [PubMed]
  • 58. Rolfe BE, Stamatiou S, World CJ, Brown L, Thomas AC, Bingley JA, Worth NF, Campbell JH. Leukaemia inhibitory factor retards the progression of atherosclerosis. Cardiovasc Res. 2003; 58:222–30. https://doi.org/10.1016/s0008-6363(02)00832-5 [PubMed]
  • 59. Fotakis P, Kothari V, Thomas DG, Westerterp M, Molusky MM, Altin E, Abramowicz S, Wang N, He Y, Heinecke JW, Bornfeldt KE, Tall AR. Anti-Inflammatory Effects of HDL (High-Density Lipoprotein) in Macrophages Predominate Over Proinflammatory Effects in Atherosclerotic Plaques. Arterioscler Thromb Vasc Biol. 2019; 39:e253–e272. https://doi.org/10.1161/ATVBAHA.119.313253 [PubMed]
  • 60. Barrett TJ, Distel E, Murphy AJ, Hu J, Garshick MS, Ogando Y, Liu J, Vaisar T, Heinecke JW, Berger JS, Goldberg IJ, Fisher EA. Apolipoprotein AI) promotes atherosclerosis regression in diabetic mice by suppressing myelopoiesis and plaque inflammation. Circulation. 2019; 140:1170–84. https://doi.org/10.1161/CIRCULATIONAHA.119.039476 [PubMed]
  • 61. Suzuki M, Pritchard DK, Becker L, Hoofnagle AN, Tanimura N, Bammler TK, Beyer RP, Bumgarner R, Vaisar T, de Beer MC, de Beer FC, Miyake K, Oram JF, Heinecke JW. High-density lipoprotein suppresses the type I interferon response, a family of potent antiviral immunoregulators, in macrophages challenged with lipopolysaccharide. Circulation. 2010; 122:1919–27. https://doi.org/10.1161/CIRCULATIONAHA.110.961193 [PubMed]
  • 62. Kotla S, Singh NK, Rao GN. ROS via BTK-p300-STAT1-PPARγ signaling activation mediates cholesterol crystals-induced CD36 expression and foam cell formation. Redox Biol. 2017; 11:350–64. https://doi.org/10.1016/j.redox.2016.12.005 [PubMed]
  • 63. Pasquin S, Laplante V, Kouadri S, Milasan A, Mayer G, Tormo AJ, Savin V, Sharma M, Martel C, Gauchat JF. Cardiotrophin-like cytokine increases macrophage-foam cell transition. J Immunol. 2018; 201:2462–71. https://doi.org/10.4049/jimmunol.1800733 [PubMed]
  • 64. Gordts PL, Foley EM, Lawrence R, Sinha R, Lameda-Diaz C, Deng L, Nock R, Glass CK, Erbilgin A, Lusis AJ, Witztum JL, Esko JD. Reducing macrophage proteoglycan sulfation increases atherosclerosis and obesity through enhanced type I interferon signaling. Cell Metab. 2014; 20:813–26. https://doi.org/10.1016/j.cmet.2014.09.016 [PubMed]
  • 65. Wu JH, Li YN, Chen AQ, Hong CD, Zhang CL, Wang HL, Zhou YF, Li PC, Wang Y, Mao L, Xia YP, He QW, Jin HJ, et al. Inhibition of Sema4D/PlexinB1 signaling alleviates vascular dysfunction in diabetic retinopathy. EMBO Mol Med. 2020; 12:e10154. https://doi.org/10.15252/emmm.201810154 [PubMed]
  • 66. Yang S, Xia YP, Luo XY, Chen SL, Li BW, Ye ZM, Chen SC, Mao L, Jin HJ, Li YN, Hu B. Exosomal CagA derived from helicobacter pylori-infected gastric epithelial cells induces macrophage foam cell formation and promotes atherosclerosis. J Mol Cell Cardiol. 2019; 135:40–51. https://doi.org/10.1016/j.yjmcc.2019.07.011 [PubMed]
  • 67. Hamczyk MR, Villa-Bellosta R, Quesada V, Gonzalo P, Vidak S, Nevado RM, Andrés-Manzano MJ, Misteli T, López-Otín C, Andrés V. Progerin accelerates atherosclerosis by inducing endoplasmic reticulum stress in vascular smooth muscle cells. EMBO Mol Med. 2019; 11:e9736. https://doi.org/10.15252/emmm.201809736 [PubMed]
Download PDF
Corresponding Authors

Fei Cao
2000xh0644@hust.edu.cn

Bo Hu
hubo@mail.hust.edu.cn

Keywords
cGAS atherosclerosis RNA-seq M1 polarization synergistic inflammatory signaling
Table of Contents