Inflammatory expression profiles in monocyte-to-macrophage differentiation in patients with systemic lupus erythematosus and relationship with atherosclerosis
- Benjamin D Korman1Email author,
- Chiang-Ching Huang1,
- Carly Skamra1,
- Peggy Wu1,
- Renee Koessler1,
- David Yao1,
- Qi Quan Huang1,
- William Pearce1,
- Kim Sutton-Tyrrell^2,
- George Kondos3,
- Daniel Edmundowicz4,
- Richard Pope†1 and
- Rosalind Ramsey-Goldman†1
© Korman et al.; licensee BioMed Central Ltd. 2014
Received: 7 March 2014
Accepted: 10 June 2014
Published: 10 July 2014
Our objectives were to examine mononuclear cell gene expression profiles in patients with systemic lupus erythematosus (SLE) and healthy controls and to compare subsets with and without atherosclerosis to determine which genes’ expression is related to atherosclerosis in SLE.
Monocytes were obtained from 20 patients with SLE and 16 healthy controls and were in vitro-differentiated into macrophages. Subjects also underwent laboratory and imaging studies to evaluate for subclinical atherosclerosis. Whole-genome RNA expression microarray was performed, and gene expression was examined.
Gene expression profiling was used to identify gene signatures that differentiated patients from controls and individuals with and without atherosclerosis. In monocytes, 9 out of 20 patients with SLE had an interferon-inducible signature compared with 2 out of 16 controls. By looking at gene expression during monocyte-to-macrophage differentiation, we identified pathways which were differentially regulated between SLE and controls and identified signatures based on relevant intracellular signaling molecules which could differentiate SLE patients with atherosclerosis from controls. Among patients with SLE, we used a previously defined 344-gene atherosclerosis signature in monocyte-to-macrophage differentiation to identify patient subgroups with and without atherosclerosis. Interestingly, this signature further classified patients on the basis of the presence of SLE disease activity and cardiovascular risk factors.
Many genes were differentially regulated during monocyte-to-macrophage differentiation in SLE patients compared with controls. The expression of these genes in mononuclear cells is important in the pathogenesis of SLE, and molecular profiling using gene expression can help stratify SLE patients who may be at risk for development of atherosclerosis.
Systemic lupus erythematosus (SLE) is a chronic multisystem autoimmune disease with a diverse array of clinical manifestations. More recently, it has become evident that SLE is associated with both accelerated atherosclerosis and an increased risk of cardiovascular complications. Women with SLE have a high prevalence of coronary heart disease, and in some cohorts this figure is as high as 6% to 10% [1–4]. The estimated incidence of new cardiovascular events in patients with SLE is approximately 1.2% to 1.5% per year, which represents a five- to six-fold age-adjusted increased risk compared with women without SLE, which is even more striking since it occurs in young pre-menopausal women . Although chronic inflammation is thought to be an important risk factor for these manifestations and a number of cytokines (including interferons) have been proposed as important in this process, the exact pathogenesis of atherosclerosis in SLE remains to be elucidated .
Monocytes recruited into tissues from peripheral blood differentiate into macrophages, which are critical in the pathogenesis of many diseases, including both atherosclerosis and SLE . Atherosclerosis has long been associated with chronic inflammation, and previous work has shown that monocytes and macrophages play an important role in the development of atherosclerosis [8, 9]. Recent studies have shown that patients with SLE have monocyte/macrophage defects involving surface protein expression, cytokine production, and phagocytic capacity, suggesting that monocytes and macrophages are key players in the pathobiology of SLE .
RNA expression microarray has been previously used to analyze mononuclear cells from patients with SLE and has identified an interferon-inducible gene expression profile which is seen in a significant proportion of patients with SLE and tends to correlate with more severe disease manifestations [11–14]. Similar analyses have been done in patients with coronary artery disease and atherosclerosis and have found strong differential expression of genes related to innate and adaptive immunity, particularly those involved in leukocyte transendothelial migration [15–18].
There are dramatic changes in gene expression profiles during monocyte-to-macrophage differentiation, including genes important in lipid metabolism and foam cell function [19, 20]. Furthermore, since monocytes differentiate into macrophages as they egress from the circulation into the inflammatory or atherosclerotic milieu, we reasoned that comparing this process in patients with SLE and controls might reveal differences that would provide insights into the increased risk for atherosclerosis seen in patients with SLE. To test this hypothesis, we examined gene expression profiles in monocytes and in vitro-differentiated macrophages in SLE patients and healthy controls with and without subclinical atherosclerosis to look for differential expression patterns that might help explain the increased coincidence of these diseases.
Materials and methods
Peripheral blood was obtained from 36 females: 20 patients with SLE and 16 healthy controls. The institutional review boards of Northwestern University and the University of Illinois at Chicago approved the study, and all participants gave written informed consent prior to their involvement.
SLE patients and control subjects were selected from the SOLVABLE study (Study of SLE Vascular and Bone Long-term Endpoints, a prospective study assessing risk for subclinical and clinical cardiovascular disease) on the basis of having either evidence of atherosclerosis or no evidence of atherosclerosis. All patients had a physician-verified SLE diagnosis and met at least four American College of Rheumatology (ACR) criteria for the classification of SLE. Controls were selected from the general population and represent healthy women who were matched to patients with SLE by age (±5 years), race/ethnicity, and zip code of residence . Controls did not have evidence of SLE or other autoimmune disease.
At study visits, subjects filled out a questionnaire in which they provided demographic, personal, and family history and information about cardiovascular risk factors and underwent routine physical examination and laboratory tests. Patient subphenotypes were determined on the basis of whether patients met each of the 11 ACR classification criteria. Information on age, self-reported race/ethnicity, smoking history, current estrogen use, current aspirin use, and menopause status was obtained from a questionnaire. (If indeterminate by questionnaire, menopause was determined by follicle-stimulating hormone status). Blood pressure was measured twice, and the mean of the two measurements was used for analysis. Height, weight, and waist/hip measurements were obtained by using the protocol in the Multi-ethnic Study of Atherosclerosis .
Laboratory tests were run on all patients and controls. These included fasting lipids (total cholesterol, high-density lipoprotein cholesterol (HDL-c), and triglyceride) and fasting glucose measured in the Lipid Laboratory at the University of Pittsburgh Graduate School of Public Health. The Friedewald equation was used to estimate low-density lipoprotein cholesterol (LDL-c) unless the triglyceride level was more than 400 mg/dL, in which case LDL-c was measured directly. Plasma glucose levels were determined by enzymatic assay, and plasma insulin levels were measured by radioimmunoassay. C-reactive protein (CRP) was measured by using an immunonephelometric assay at the Laboratory for Clinical Biochemistry Research at the University of Vermont. C3 and C4 were measured by nephelometry, and double-stranded DNA (dsDNA) antibodies were measured by using the Crithidia luciliae method. dsDNA antibodies were considered positive if the titer was at least 1:10.
All patients and controls were examined and validated measures of SLE disease activity (Systemic Lupus Erythematosus Disease Activity Index 2000, or SLEDAI-2K) as well as disease damage (Systemic Lupus International Collaborating Clinics/American College of Rheumatology Damage Index, or SLICC/ACR-DI), which were completed by trained physicians. The SLEDAI-2K and SLICC/ACR-DI instruments were used to calculate disease activity and damage, respectively [23, 24]. Disease duration was calculated by using the date the subject fulfilled the 4th ACR classification criteria for SLE for disease onset and the study visit date as the end date. Participants provided information on corticosteroid treatment (current use and duration of treatment) as well as current use of hydroxychloroquine and immunosuppressants (which included cyclophosphamide, azathioprine, methotrexate, mycophenolate mofetil, cyclosporine, and tacrolimus). Renal disease was defined as present if the subject had fulfilled ACR classification criteria for SLE renal involvement (greater than 0.5 gm/day or 3 or more proteinuria or the presence of cellular casts or a combination of these) or had a renal biopsy with evidence of World Health Organization class IIb, III, IV, or V SLE nephritis. International Society of Nephrology/Renal Pathology Society classification was not used, because most renal biopsies were obtained prior to its widespread use.
All 36 individuals underwent imaging evaluations of their carotid arteries, coronary arteries, and aorta. For the purpose of this study, we defined an atherosclerosis phenotype by the presence of at least three of the following four abnormalities on carotid ultrasound or electron beam computed tomography (EBCT): a score of more than 0 on the carotid plaque index, an intima-media thickness (IMT) of greater than mean of the study group (0.67 mm), higher coronary artery calcium (CAC) score (>10), or higher aorta calcium (AC) score (>100) [25, 26].
In the coronary arteries and aorta, EBCT scanning was performed to measure vascular calcium by using an Imatron C-150 Ultrafast CT scanner (General Electric Medical Systems, South San Francisco, CA, USA). Calcium scores were calculated with a densitometric program available on the Imatron C-150 scanner by using the Agatston method. The outcome measures used for analyses were dichotomized with a low CAC score of not more than 10 and an AC score of not more than 100 and a high CAC score of greater than 10 and an AC score of greater than 100 signifying the presence of atherosclerosis at these arterial beds [25, 26].
Subclinical cardiovascular disease was measured in the carotid arteries by B-mode ultrasound performed by trained sonographers. Carotid ultrasound was performed, and B-mode images of the right and left carotid artery, carotid bifurcation, and the first centimeter of the internal carotid were obtained in multiple planes and assessed for plaque. Carotid ultrasound was completed in a single center. Ultrasonographic measurements were performed by using a Siemens Sequioa model C256 (Siemens, Munich, Germany) equipped with a transducer 8 L5. Sonographers scanned the right and left common carotid artery, carotid bulb, and the first 1.5 cm of the internal and external carotid arteries. For each location, the sonographer imaged the vessel in multiple planes. The sonographers scored the ultrasound images for plaque in the distal common, carotid bulb, and internal arteries. Plaque was defined as a distinct area protruding into the vessel lumen with at least 50% greater thickness than that found in surrounding areas. IMT was measured by using specialized reading software across 1-cm segments of both the right and left sides of the near and far walls of the distal common carotid artery and the far wall of the carotid bulb and internal carotid artery. The mean of all average IMT readings across the eight sites was used for analysis. Reproducibility of IMT and plaque was assessed in five women who underwent two ultrasound examinations within 1 week. Each time, the scanning was performed by a different sonographer, and each scan was scored by two readers. Sonographers and readers were trained as part of a reproducibility protocol study in carotid duplex scanning .
Mononuclear cells were isolated from peripheral blood of SLE patients and healthy controls, and then monocytes were isolated and in vitro-differentiated into macrophages . Specifically, peripheral blood was drawn in tubes containing citrate-phosphate-dextrose anti-coagulant solution (Sigma-Aldrich, St. Louis, MO, USA). Mononuclear cells were isolated by Histopaque-1077 density gradient centrifugation. Two different approaches were applied to isolate monocytes. The first approach employed a Human CD14 Positive Selection Kit (StemCell Technologies, Vancouver, BC, Canada), which isolates monocytes by directly binding to anti-CD14 antibody-coupled magnetic beads. The cells attached to the beads are immediately frozen in Trizol at -70°C for future RNA isolation. Monocytes obtained from this approach were more than 95% pure. To isolate monocytes for in vitro differentiation into macrophages, negative selection was employed (Human Monocyte Enrichment Kit Without CD16 Deletion) and resulted in improved survival during differentiation. The mononuclear cells were incubated with an antibody cocktail including antibodies to CD2, CD3, CD20, CD56, CD66b, CD123, and glycophorin attached to magnetic beads, which were employed because they resulted in higher purity. Cells positive for these markers were depleted, leaving CD14+ monocytes from peripheral blood. The isolated CD14+ cells were allowed to attach to plastic for 1 hour and then were in vitro-differentiated in RPMI-1640 culture medium supplemented with 20% fetal bovine serum, 10 ng/mL macrophage colony-stimulating factor, 1 μg/mL polymyxin-B, 100 U penicillin, and 100 mg/mL streptomycin for 7 days at 37°C in 5% CO2.
RNA was extracted from monocytes and macrophages by Trizol reagent (Invitrogen, Carlsbad, CA, USA). The quality of RNA was determined by 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA), evaluated by the RNA integrity number of at least 8.0 and the ratio of 28 s and 18 s ribosomal RNA. RNA that passed quality control was then used to create cDNA and biotinylated cRNA by using an Illumina TotalPrep RNA Amplification Kit (Illumina, San Diego, CA, USA). RNA was hybridized to a Sentrix® Human-6 Expression BeadChip (Illumina) to measure expression of more than 47,000 transcripts in each sample.
Microarray data were normalized and corrected for batch effect; we then performed significance analysis of microarrays (SAM), multidimensional scaling (MDS) analysis, cluster analysis, and gene ontology (GO) analysis. These analyses were performed first in monocytes and macrophages separately and then during monocyte-to-macrophage differentiation (the difference in expression between an individual’s monocytes and macrophages). For each of these groups, the comparisons included were (a) SLE versus controls, (b) all individuals with versus without an atherosclerosis phenotype, and (c) SLE patients with versus without an atherosclerosis phenotype.
The microarray data were normalized by a quantile normalization procedure by using the bioconductor package affy, and batch effect was corrected by using the ComBat algorithm . Cluster analysis was performed to identify gene expression patterns by using CLUSTER and TREEVIEW software. GO analysis was performed by using the online application PANTHER (Protein Analysis Through Evolutionary Relationships) . The significance of enrichment of differentially expressed genes for each GO term was determined by the Bonferroni-corrected chi-squared P value as the default setting in PANTHER, whereas the significance of differential expression was determined by nominal P value of less than 0.05 from the t test and a fold change of greater than 2. MDS analysis was performed to reveal the sample relation by using global gene expression profiles. For monocyte-to-macrophage differentiation analysis, we calculated fold change for each probe as the ratio of macrophage expression over monocyte expression from each pair of monocyte-macrophage samples (from the same subject). With these fold-change data, we excluded transcripts that have average fold changes of less than 2 in both SLE patient and control groups. This filtering process resulted in 3,044 (out of 47,231) transcripts for data analysis. To determine the fold-change difference between SLE patients and controls for each of these 3,044 transcripts, we calculated the P value by using the two-sample t test. The microarray data are uploaded to Gene Expression Omnibus with accession number GSE37356.
Demographic, clinical, laboratory, and imaging results
Demographic, clinical, laboratory, medication, and imaging data from 20 systemic lupus erythematosus patients and 16 controls
SLE mean ± SD
Control mean ± SD
46.9 ± 8.9
52.8 ± 11.0
Race, % who were Caucasian
Body mass index, kg/m2
27.7 ± 7.7
28.3 ± 4.9
0.85 ± 0.07
0.86 ± 0.08
Systolic blood pressure, mm Hg
117.4 ± 13.2
116.8 ± 14.4
Diastolic blood pressure, mm Hg
72.3 ± 9.0
71.8 ± 11.5
Current smoking, %
Family history of cardiovascular disease, %
Total cholesterol, mg/dL
188.9 ± 51.9
218.9 ± 41.4
101.5 ± 41.5
139.0 ± 8.7
125.6 ± 99.9
101.0 ± 48.9
98.4 ± 10.7
103.7 ± 13.9
Glomerular filtration rate, mL/min
80.7 ± 17.8
77.9 ± 18.8
C-reactive protein, mg/L
3.5 ± 3.6
3.8 ± 6.3
101.5 ± 21.5
20.3 ± 8.4
dsDNA (crithidia) level
Presence of carotid plaque
Higher CAC score (>10), %
Higher AC score (>100), %
Intima-media thickness (mean ± SD)
0.65 ± 0.14
0.71 ± 0.21
Atherosclerosis phenotype, %
4.2 ± 4.3
2.2 ± 1.8
Disease duration, years
16.3 ± 8.2
Total ACR SLE classification criteria (median)
Gene expression profiles in monocytes and monocyte-derived macrophages
Biologic processes in which genes associated with differential expression during monocyte-to-macrophage differentiation
Genes in PANTHER database
Immune system process
9.98 × 10-22
5.75 × 10-13
2.18 × 10-7
B cell-mediated immunity
2.49 × 10-6
Lipid metabolic process
4.59 × 10-6
6.08 × 10-6
1.97 × 10-5
Response to interferon-gamma
6.83 × 10-5
Induction of apoptosis
1.58 × 10-4
5.33 × 10-4
5.43 × 10-4
Of the biologically relevant pathways observed, the signal transduction pathway contained 147 differentially regulated genes, the most of any pathway (Table 2). We looked at the expression pattern of differentially expressed genes in this pathway and found that there were clear differences between the SLE patients and the controls for genes that were either upregulated or downregulated during macrophage differentiation (Figure 1C, D). When cluster analysis was performed with both upregulated and downregulated genes in this pathway, SLE patients and controls were distinguishable. In the case of downregulated genes, the clusters were able to differentiate 15 of the 20 patients from all of the 16 controls and five SLE patients of whom none had an atherosclerosis phenotype (Figure 1C). This pattern was not seen when genes from other relevant pathways were analyzed.
Correlations with atherosclerosis
Unsupervised expression profiling in monocytes and macrophages by using the 1,000 most differentially expressed genes—those with the highest coefficient of variation (Cv)—was unable to clearly distinguish individuals with and without atherosclerosis. However, in monocytes from patients with SLE, we demonstrated that patients with subclinical atherosclerosis demonstrated substantial enrichment of genes identified as associated with atherosclerosis based on a previously described 344-gene signature used to differentiate individuals with and without a substantial burden of atherosclerosis . A complete gene list for this signature is detailed in Additional file 6. When used to look at all individuals, this signature did not effectively delineate disease status or atherosclerosis status (Additional file 7). We found that only 140 (4.6%) of 3,044 genes demonstrated a P value of less than 0.05 when comparing the 19 participants with atherosclerosis with the 17 participants without atherosclerosis (Additional file 3A). Among the 20 patients with SLE, we found 163 significant genes (5.3%) with a P value of less than 0.05 between those with or without atherosclerosis (Additional file 3B). Among the 16 controls, we found only 74 (2.4%) genes with a P value of less than 0.05 between those with or without atherosclerosis (Additional file 3C). Of the genes identified in each of these analyses, 11 (7.8%) out of 140, 22 (13.5%) out of 163 genes, and three (1.3%) out of 74 overlapped the previously identified 344-gene atherosclerosis signature (Additional file 3D). Upon hierarchical clustering analysis, these differentially expressed genes did not differentiate any of the groups by atherosclerosis status.
Clinical characteristics of patients stratified by microarray clusters using 344-gene atherosclerosis signature
45.8 ± 9.3
46.9 ± 5.4
49.0 ± 7.4
45.0 ± 13.0
13.6 ± 9.6
21.2 ± 3.31
17.3 ± 8.8
13.1 ± 9.2
SLE disease activity/damage
8.7 ± 6.4
7.0 ± 3.8
2.0 ± 3.1
2.7 ± 2.4
2.3 ± 2.3
4.0 ± 0.8
2.3 ± 1.6
0.7 ± 1.2
Carotid intima-media thickness
Coronary artery calcium score
124 ± 85
349 ± 786
103 ± 252
Aortic calcium score
2,638 ± 3,261
828 ± 995
89 ± 149
Traditional cardiovascular risk factors
Average number of risk factors
Current medication use
When we compared global gene expression of patients with SLE and controls, we confirmed the importance of type-1 interferon signature in monocytes. In addition, we found that signal transduction, macrophage activation, response to interferon-gamma, and apoptosis were differentially regulated during the process of monocyte-to-macrophage differentiation.
About half of patients possess the interferon signature and it appears to be associated with more severe disease, including renal disease, which is a leading cause of morbidity and mortality in SLE. Although previous work has shown a mechanistic link between interferon-alpha and atherosclerosis in SLE with interferon promoting macrophage lipid uptake and foam cell formation , interestingly, we did not see an association between this signature and atherosclerosis and instead saw that atherosclerosis was more common in individuals who did not have an interferon signature with seven out of 11 without the interferon signature having atherosclerosis versus two out of nine among individuals with the interferon signature. Although the number of patients who expressed a cytokine/chemokine expression signature is small (four out of 20), the genes seen differentially expressed in this subgroup suggest that alteration in innate rather than adaptive immunity may play an important role in the pathogenesis of SLE in a subset of patients. In addition to having high disease activity, these individuals tended to have slightly shorter disease duration, and none had an atherosclerosis phenotype (Table 3). Although larger studies will be needed to confirm these types of gene expression signature associations with clinical data, if robust, they could potentially be used help stratify patients’ risk for various disease manifestations or risk of atherosclerosis.
Interestingly, we also identified lipid and carbohydrate metabolism genes as differentially expressed in SLE patients compared with controls. One could hypothesize that these genes’ expression may play a role in how these cells participate in the processes of energy metabolism and mitochondrial function during cellular differentiation or during atheroma formation. These observations suggest that monocytes from patients with SLE upon entering sites of inflammation or vascular injury may be programmed as they differentiate into macrophages with an enhanced pathogenic potential to promote an atherosclerosis phenotype.
Signal transduction genes during monocyte-to-macrophage differentiation were predictive of the presence or absence of SLE; differentially downregulated genes were able to identify 15 of 20 SLE patients from the controls and five SLE patients who did not have atherosclerosis (Figure 1B). Many of the genes differentially regulated during monocyte-to-macrophage differentiation (downregulated genes include JAK2, STAT6, TLR8, and TLR 2, and upregulated genes include VEGFB, TGFB1, FN1, IL-1R2, SCARB1, MSR1, and CD163) have been previously reported in the pathogenesis of atherosclerosis [16, 18, 31, 32].
Foamy macrophages in atherosclerotic plaques demonstrate both M1 and M2 markers, and similar trends were observed during monocyte-to-macrophage differentiation in SLE. Mechanistically, we found that JAK2 and STAT6 were suppressed to a greater degree during monocyte-to-macrophage differentiation in patients with SLE compared with controls, especially among patients with atherosclerosis. JAK2 and STAT6 are both important in mediating M2 macrophage differentiation . Expression of SOCS3, a JAK/STAT suppressor whose deficiency leads to M1 differentiation, was also reduced to a greater degree in patients with SLE. Furthermore, IL15, which mediates its signal through STATs including STAT6, was reduced more during monocyte-to-macrophage differentiation in patients with SLE, especially in those with atherosclerosis. The JAK2 fold-change expression was highly correlated with IL15 and IFITM1, which supports the internal consistency of the data. These results indicate that SLE patients’ monocytes favor M1 rather than M2 differentiation compared with controls. Furthermore, we found that, compared with controls, patients with SLE had significantly increased levels of vascular endothelial growth factor (VEGFB) which promotes ischemic myocardial revascularization . CD163, an M2 marker, is highly expressed in atherosclerotic plaques, whereas MSR1 is important in plaque formation [32, 34]. Patients with SLE therefore appear to have an increased predilection for creating inflammatory M1 macrophages but also exhibit features associated with M2 differentiation which may contribute to the pathogenesis of atherosclerosis.
One inherent limitation of the present study is the small sample size. Although clinical heterogeneity in SLE is also a concern, the SLE patients analyzed in this study were representative of the larger SOLVABLE population of 165 individuals who did not participate in this study in terms of clinical parameters and medication use. Because of potential confounding from the relatively small sample size, we looked at other variables, including medication usage, given the higher use of immunomodulatory medications and statin medications in the patients with SLE. Although these differences may contribute to some of the differences observed, we did not find that individuals’ expression clustered by medication use, and clustering based on medication status was not able to identify patients with or without atherosclerosis (data not shown). Given that patients on medication still had more severe disease, the medications may have dampened these patients’ inflammatory signatures.
Monocyte gene expression patterns of patients with SLE and subclinical atherosclerosis were associated with a 344-gene signature previously reported as associated with atherosclerosis (Figure 2) . This signature includes primarily immune and inflammatory genes, including six Toll-like receptor genes, IL1B, IRAK3, and MAPK14, and a number of genes involved in apoptosis and mobilization of calcium. Strikingly, hierarchical clustering using this gene signature divided patients into groups based on their disease activity, and within these clusters, those with and without an atherosclerosis phenotype. Although the presence/absence of atherosclerosis was not absolute in any group, the differences between groups were readily apparent. Differences in patient demographics such as age and disease duration were not statistically significant and cannot fully account for the differences we saw between clusters in terms of atherosclerosis. We did note a trend (p = 0.14) for clusters 2 and 3 to have longer standing disease, and therefore without longitudinal follow-up of clusters 1 and 4, it is difficult to say whether these individuals will develop atherosclerosis at similar rates once they have had longer disease duration. Cluster 2 (high disease activity, 100% renal disease, and abundant atherosclerosis), like cluster 3 (low disease activity and abundant atherosclerosis), had substantial traditional cardiovascular risk factors (diabetes, smoking, dyslipidemia, and hypertension). The data suggest that both disease activity and traditional risk factors appear to contribute to cardiovascular risk in patients with SLE and raise the question of whether these risks are additive. One potential confounder for this gene expression pattern could be medication use, especially because immune modulating medications may affect the expression of genes relevant to the signatures being examined. For example, patients in cluster 2 were most likely to be on steroids (P = 0.03) and immunosuppressants (p = 0.004), which could affect their expression profile. Also, the use of statins in patients with SLE may have had an effect on gene expression and was likely to have had an effect on progression (and possibly development) of atherosclerosis . These data are intriguing and hypothesis-generating but will need to be validated in a larger cohort of patients than was possible in this pilot study.
Our data support the importance of the interferon-inducible signature in SLE and suggest that gene expression profiles in monocytes from patients with SLE are different between individuals with high and low disease activity and with and without atherosclerosis. The changes in expression seen during monocyte-to-macrophage differentiation suggest that this process of differentiation may contribute to disease pathogenesis, possibly by polarizing macrophages toward classic M1 activation. Further studies may be directed to determine whether factors present in SLE serum, genetic predisposition, or stimuli such as the stress response of oxidized lipid uptake might more clearly identify factors that promote the differentiation of pathogenic macrophages in SLE.
American College of Rheumatology
coronary artery calcium
electron beam computed tomography
intima-media thickness (carotid)
low-density lipoprotein cholesterol
Protein Analysis Through Evolutionary Relationships
significance analysis of microarrays
systemic lupus erythematosus
systemic lupus erythematosus disease activity index
Systemic Lupus International Collaborating Clinics/American College of Rheumatology Damage Index
Study of Lupus Vascular and Bone Long-term Endpoints
- VEGFB :
vascular endothelial growth factor.
Funding was provided by National Institutes of Health grants T32 AR 07611, AR055240, P60 2 AR30692, P60 AR30492, K24 AR 002138, and UL1RR025741 and The Mary Kirkland Center for SLE Research. Expression analysis services were performed at the Northwestern University Genomics Core Facility.
- Esdaile JM, Abrahamowicz M, Grodzicky T, Li Y, Panaritis C, du Berger R, Cote R, Grover SA, Fortin PR, Clarke AE, Senécal JL: Traditional Framingham risk factors fail to fully account for accelerated atherosclerosis in systemic lupus erythematosus. Arthritis Rheum. 2001, 44: 2331-2337.View ArticlePubMedGoogle Scholar
- Manzi S, Meilahn EN, Rairie JE, Conte CG, Medsger TA, Jansen-McWilliams L, D’Agostino RB, Kuller LH: Age-specific incidence rates of myocardial infarction and angina in women with systemic lupus erythematosus: comparison with the Framingham Study. Am J Epidemiol. 1997, 145: 408-415.View ArticlePubMedGoogle Scholar
- Petri M, Perez-Gutthann S, Spence D, Hochberg MC: Risk factors for coronary artery disease in patients with systemic lupus erythematosus. Am J Med. 1992, 93: 513-519.View ArticlePubMedGoogle Scholar
- Petri M, Spence D, Bone LR, Hochberg MC: Coronary artery disease risk factors in the Johns Hopkins Lupus Cohort: prevalence, recognition by patients, and preventive practices. Medicine. 1992, 71: 291-302.View ArticlePubMedGoogle Scholar
- Roman MJ, Shanker BA, Davis A, Lockshin MD, Sammaritano L, Simantov R, Crow MK, Schwartz JE, Paget SA, Devereux RB, Salmon JE: Prevalence and correlates of accelerated atherosclerosis in systemic lupus erythematosus. N Engl J Med. 2003, 349: 2399-2406.View ArticlePubMedGoogle Scholar
- Li J, Fu Q, Cui H, Qu B, Pan W, Shen N, Bao C: Interferon-alpha priming promotes lipid uptake and macrophage-derived foam cell formation: a novel link between interferon-alpha and atherosclerosis in lupus. Arthritis Rheum. 2011, 63: 492-502.View ArticlePubMedGoogle Scholar
- McMahon M, Hahn BH: Atherosclerosis and systemic lupus erythematosus: mechanistic basis of the association. Curr Opin Immunol. 2007, 19: 633-639.PubMed CentralView ArticlePubMedGoogle Scholar
- Libby P: Inflammation in atherosclerosis. Nature. 2002, 420: 868-874.View ArticlePubMedGoogle Scholar
- Woollard KJ, Geissmann F: Monocytes in atherosclerosis: subsets and functions. Nat Rev Cardiol. 2010, 7: 77-86.PubMed CentralView ArticlePubMedGoogle Scholar
- Li Y, Lee PY, Reeves WH: Monocyte and macrophage abnormalities in systemic lupus erythematosus. Arch Immunol Ther Exp. 2010, 58: 355-364.View ArticleGoogle Scholar
- Crow MK: Type I interferon in systemic lupus erythematosus. Curr Top Microbiol Immunol. 2007, 316: 359-386.PubMedGoogle Scholar
- Bennett L, Palucka AK, Arce E, Cantrell V, Borvak J, Banchereau J, Pascual V: Interferon and granulopoiesis signatures in systemic lupus erythematosus blood. J Exp Med. 2003, 197: 711-723.PubMed CentralView ArticlePubMedGoogle Scholar
- Pascual V, Farkas L, Banchereau J: Systemic lupus erythematosus: all roads lead to type I interferons. Curr Opin Immunol. 2006, 18: 676-682.View ArticlePubMedGoogle Scholar
- Baechler EC, Batliwalla FM, Karypis G, Gaffney PM, Ortmann WA, Espe KJ, Shark KB, Grande WJ, Hughes KM, Kapur V, Gregersen PK, Behrens TW: Interferon-inducible gene expression signature in peripheral blood cells of patients with severe lupus. Proc Natl Acad Sci U S A. 2003, 100: 2610-2615.PubMed CentralView ArticlePubMedGoogle Scholar
- Sinnaeve PR, Donahue MP, Grass P, Seo D, Vonderscher J, Chibout SD, Kraus WE, Sketch M, Nelson C, Ginsburg GS, Goldschmidt-Clermont PJ, Granger CB: Gene expression patterns in peripheral blood correlate with the extent of coronary artery disease. Plos One. 2009, 4: e7037-PubMed CentralView ArticlePubMedGoogle Scholar
- Wingrove JA, Daniels SE, Sehnert AJ, Tingley W, Elashoff MR, Rosenberg S, Buellesfeld L, Grube E, Newby LK, Ginsburg GS, Kraus WE: Correlation of peripheral-blood gene expression with the extent of coronary artery stenosis. Circ Cardiovasc Genet. 2008, 1: 31-38.View ArticlePubMedGoogle Scholar
- Charchar FJ, Bloomer LD, Barnes TA, Cowley MJ, Nelson CP, Wang Y, Denniff M, Debiec R, Christofidou P, Nankervis S, Dominiczak AF, Bani-Mustafa A, Balmforth AJ, Hall AS, Erdmann J, Cambien F, Deloukas P, Hengstenberg C, Packard C, Schunkert H, Ouwehand WH, Ford I, Goodall AH, Jobling MA, Samani NJ, Tomaszewski M: Inheritance of coronary artery disease in men: an analysis of the role of the Y chromosome. Lancet. 2012, 379: 915-922.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang CC, Liu K, Pope RM, Du P, Lin S, Rajamannan NM, Huang QQ, Jafari N, Burke GL, Post W, Watson KE, Johnson C, Daviglus ML, Lloyd-Jones DM: Activated TLR signaling in atherosclerosis among women with lower Framingham risk score: the multi-ethnic study of atherosclerosis. Plos One. 2011, 6: e21067-PubMed CentralView ArticlePubMedGoogle Scholar
- Liu H, Shi B, Huang CC, Eksarko P, Pope RM: Transcriptional diversity during monocyte-to-macrophage differentiation. Immunol Lett. 2008, 117: 70-80.PubMed CentralView ArticlePubMedGoogle Scholar
- Martinez FO, Gordon S, Locati M, Mantovani A: Transcriptional profiling of the human monocyte-to-macrophage differentiation and polarization: new molecules and patterns of gene expression. J Immunol. 2006, 177: 7303-7311.View ArticlePubMedGoogle Scholar
- Rhew EY, Lee C, Eksarko P, Dyer AR, Tily H, Spies S, Pope RM, Ramsey-Goldman R: Homocysteine, bone mineral density, and fracture risk over 2 years of followup in women with and without systemic lupus erythematosus. J Rheumatol. 2008, 35: 230-236.PubMedGoogle Scholar
- Bild DE, Bluemke DA, Burke GL, Detrano R, Diez Roux AV, Folsom AR, Greenland P, Jacob DR, Kronmal R, Liu K, Nelson JC, O’Leary D, Saad MF, Shea S, Szklo M, Tracy RP: Multi-ethnic study of atherosclerosis: objectives and design. Am J Epidemiol. 2002, 156: 871-881.View ArticlePubMedGoogle Scholar
- Gladman D, Ginzler E, Goldsmith C, Fortin P, Liang M, Urowitz M, Bacon P, Bombardieri S, Hanly J, Hay E, Isenberg D, Jones J, Nived O, Petri M, Richter M, Sanchez-Guerrero J, Snaith M, Sturfelt G, Symmons D: Systemic lupus international collaborative clinics: development of a damage index in systemic lupus erythematosus. J Rheumatol. 1992, 19: 1820-1821.PubMedGoogle Scholar
- Gladman DD, Ibanez D, Urowitz MB: Systemic lupus erythematosus disease activity index 2000. J Rheumatol. 2002, 29: 288-291.PubMedGoogle Scholar
- Budoff MJ, McClelland RL, Nasir K, Greenland P, Kronmal RA, Kondos GT, Shea S, Lima JA, Blumenthal RS: Cardiovascular events with absent or minimal coronary calcification: the Multi-Ethnic Study of Atherosclerosis (MESA). Am Heart J. 2009, 158: 554-561.PubMed CentralView ArticlePubMedGoogle Scholar
- Detrano RC, Anderson M, Nelson J, Wong ND, Carr JJ, McNitt-Gray M, Bild DE: Coronary calcium measurements: effect of CT scanner type and calcium measure on rescan reproducibility–MESA study. Radiology. 2005, 236: 477-484.View ArticlePubMedGoogle Scholar
- Thompson TS-TK, Wildman R: Continuous quality assessment programs can improve carotid duplex scan quality. J Vas Tech. 2001, 25: 33-39.Google Scholar
- Johnson WE, Li C, Rabinovic A: Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007, 8: 118-127.View ArticlePubMedGoogle Scholar
- The PANTHER (Protein ANalysis THrough Evolutionary Relationships) Classification System. [http://www.pantherdb.org/]
- National Cholesterol Education Program Expert Panel on Detection E, Treatment of High Blood Cholesterol in A: Third report of the national cholesterol education program (NCEP) expert panel on detection, evaluation, and treatment of high blood cholesterol in adults (adult treatment panel III) final report. Circulation. 2002, 106: 3143-3421.Google Scholar
- Li X, Tjwa M, Van Hove I, Enholm B, Neven E, Paavonen K, Jeltsch M, Juan TD, Sievers RE, Chorianopoulos E, Wada H, van Wildemeersch M, Noel A, Foidart JM, Springer ML, von Degenfeld G, Dewerchin M, Blau HM, Alitalo K, Eriksson U, Carmeliet P, Moons L: Reevaluation of the role of VEGF-B suggests a restricted role in the revascularization of the ischemic myocardium. Arterioscler Thromb Vasc Biol. 2008, 28: 1614-1620.PubMed CentralView ArticlePubMedGoogle Scholar
- Neyen C, Pluddemann A, Roversi P, Thomas B, Cai L, van der Westhuyzen DR, Sim RB, Gordon S: Macrophage scavenger receptor A mediates adhesion to apolipoproteins A-I and E. Biochemistry. 2009, 48: 11858-11871.PubMed CentralView ArticlePubMedGoogle Scholar
- O’Sullivan LA, Liongue C, Lewis RS, Stephenson SE, Ward AC: Cytokine receptor signaling through the Jak-Stat-Socs pathway in disease. Mol Immunol. 2007, 44: 2497-2506.View ArticlePubMedGoogle Scholar
- Stoger JL, Gijbels MJ, van der Velden S, Manca M, van der Loos CM, Biessen EA, Daemen MJ, Lutgens E, de Winther MP: Distribution of macrophage polarization markers in human atherosclerosis. Atherosclerosis. 2012, 225: 461-468.View ArticlePubMedGoogle Scholar
- Kleemann R, Bureeva S, Perlina A, Kaput J, Verschuren L, Wielinga PY, Hurt-Camejo E, Nikolsky Y, van Ommen B, Kooistra T: A systems biology strategy for predicting similarities and differences of drug effects: evidence for drug-specific modulation of inflammation in atherosclerosis. BMC Syst Biol. 2011, 5: 125-PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.