Reprogramming of bone marrow myeloid progenitor cells in patients with severe coronary artery disease
Abstract
Atherosclerosis is the major cause of cardiovascular disease (CVD). Monocyte-derived macrophages are the most abundant immune cells in atherosclerotic plaques. In patients with atherosclerotic CVD, leukocytes have a hyperinflammatory phenotype. We hypothesize that immune cell reprogramming in these patients occurs at the level of myeloid progenitors. We included 13 patients with coronary artery disease due to severe atherosclerosis and 13 subjects without atherosclerosis in an exploratory study. Cytokine production capacity after ex vivo stimulation of peripheral blood mononuclear cells (MNCs) and bone marrow MNCs was higher in patients with atherosclerosis. In BM-MNCs this was associated with increased glycolysis and oxidative phosphorylation. The BM composition was skewed towards myelopoiesis and transcriptome analysis of HSC/GMP cell populations revealed enrichment of neutrophil- and monocyte-related pathways. These results show that in patients with atherosclerosis, activation of innate immune cells occurs at the level of myeloid progenitors, which adds exciting opportunities for novel treatment strategies.
Introduction
Atherosclerotic cardiovascular disease (CVD), including myocardial infarction and stroke, is the leading cause of death worldwide. Major risk factors for CVD include dyslipoproteinemia, smoking, hypertension, obesity, and diabetes. However, despite optimal pharmacological treatment of these risk factors, a considerable residual CVD risk remains. Evidence is rapidly accumulating that chronic low-grade inflammation of the vascular wall is a key pathophysiological component of atherosclerosis, and that treatment with anti-inflammatory drugs, such as canakinumab and colchicine, additionally lowers future CVD risk (Nidorf et al., 2020; Ridker et al., 2017; Tardif et al., 2019). More detailed knowledge of this inflammatory process and the role of individual immune cells would allow the development of more specific and safe anti-inflammatory therapies.
Within atherosclerotic plaques, monocyte-derived macrophages are the most abundant immune cells, which play a key role in initiation, progression, and destabilization of these plaques, although the mechanisms driving the persistent inflammatory activation are poorly understood (Moore et al., 2013). In patients with CVD or risk factors for CVD, circulating leukocytes show differences in composition and individual cell phenotypes, compared to healthy individuals. Isolated monocytes from patients with CVD are characterized by an increased cytokine production capacity and higher glycolytic metabolism (Bekkering et al., 2016; Elsenberg et al., 2013; Shirai et al., 2016). This also holds true for patients with increased CVD risk due to familial hypercholesterolemia (Bekkering et al., 2019). In addition, higher circulating levels of granulocytes are associated with an increased CVD risk in a general population study (Fani et al., 2020).
Experimental studies in animal models of atherosclerosis have revealed reprogramming of hematopoietic stem and progenitor cells (HSPCs) in the bone marrow as a cause of these changes in circulating innate immune cells under conditions that predispose to CVD, including hyperlipidemia (Christ et al., 2018), and stress (Heidt et al., 2014). Western-type diet feeding of low-density lipoprotein receptor deficient (Ldlr-/-) mice induces inflammatory activation of monocytes by epigenetic reprogramming of HSPCs which persists in normocholesterolemic conditions (Christ et al., 2018). This resembles the finding of persistent functional and transcriptional hyperresponsive monocytes in patients with dyslipidemia despite cholesterol-lowering treatment (Bekkering et al., 2019). In a retrospective analysis of bone marrow samples from cancer patients, HSPCs had a higher proliferative potential in patients with known atherosclerotic CVD than in control patients (van der Valk et al., 2016). A detailed assessment, however, of HSPC composition and function in patients with atherosclerosis is lacking.
Here we performed a comprehensive study in patients with and without coronary artery disease (CAD). We obtained circulating leukocytes and HSPCs and assessed composition with flow cytometry, as well as cytokine production capacity, metabolism, and the transcriptional profile. In addition, we performed 2’-deoxy-2’-[18F]fluoro-D-glucose positron-emission-tomography ([18F]FDG PET/CT) scanning to detect bone marrow and splenic activity. We report that bone marrow mononuclear cells in patients with atherosclerosis are characterized by enhanced cytokine production capacity and an increased metabolic rate. Flow cytometry and RNAseq analysis revealed inflammatory transcriptional reprogramming and myeloid skewing. These results show that in patients with atherosclerosis, activation of the innate immune system occurs at the level of bone marrow myeloid progenitors, which adds exciting opportunities for novel treatment strategies.
Results
Group characteristics
Thirteen individuals with CAD due to severe atherosclerosis with a median Total Plaque Score (TPS) of 14[9-15], and thirteen individuals without coronary atherosclerosis (TPS = 0) participated in the study (Table 1). As expected, patients with atherosclerosis more often received lipid-lowering therapy with statins, which was associated with a lower total and LDL-cholesterol concentration. All outcomes were corrected for the age and BMI because these are known modulators of innate immune cell function (Ter Horst et al., 2016).
Circulating inflammatory markers
The distribution of immune cells and monocyte subpopulations was similar between groups (Table 2). Integrin CD11b expression on monocytes tended to be higher in patients with CAD (p=0.09). Circulating endothelial dysfunction marker E-selectin was higher in patients with atherosclerosis (p<0.05).
Cytokine production capacity of circulating PBMCs is higher in atherosclerosis
Previously, we showed in a comparable study cohort that the LPS-induced production of IL-6, TNFα, IL-1β, and IL-8 in PBMCs was higher in patients with CAD compared to controls (Bekkering et al., 2016). We could confirm this in the current study for IL-8 (p<0.01), with a similar pattern for IL-6 (p=0.08) and TNFα (p=0.13) (Figure 1A). There were no significant differences in response to TLR2 agonist Pam3Cys stimulation (Figure 1B).
HSPC composition is changed in patients with CAD
We assessed the HSPC composition using flow cytometry. Although the total percentage of HSPCs was similar between groups, the percentage of multipotent progenitors (HSC/MPP, p<0.05) was higher in patients with CAD (Figure 2A,B). Trends for a higher percentage of common myeloid progenitors (CMP, p=0.06) and common lymphoid progenitors (CLP, p=0.08) were observed in patients with CAD (Figure 2C,D). In addition, the percentage of pre-monocytes was higher (p<0.05), which was associated with lower megakaryocyte erythrocyte progenitor percentages (MEP, p<0.05) (Figure 2F,G). Within the bone marrow, no difference in the percentage of monocytes (3.3 [3.1–4.4]% in patients versus 3.4 [3.0–4.6] in controls), and monocyte subpopulations, that is, classical (87 [82–88] % versus 82 [79-89]), intermediate (7.5 [7.1–10.8]% versus 9.3 [6.1–11.5]) and non-classical monocytes (5.1 [4.2–8.1]% versus 8.1 [3.6–11.4]), was observed. Interestingly, the percentage of circulating multipotent progenitors (HSC/MPP) was higher in patients with CAD (p<0.05) (Figure 2I). Similar patterns were found using absolute cell counts of progenitor populations.
Functional and metabolic reprogramming of bone marrow MNCs in atherosclerosis
The immune response of BM-MNCs was determined after ex vivo 24-hour stimulation with TLR agonists. TNFα production after LPS stimulation was higher in patients with atherosclerosis (p<0.05), with similar patterns for IL-6 (p=0.12), IL-8 (p=0.10), and IL-1Ra (p=0.11) although these differences did not reach statistical significance (Figure 3A). There were no significant differences in response to Pam3Cys stimulation (Figure 3B).
Seahorse respirometry revealed that both the basal and maximal OCR was higher (p<0.01 and p<0.05) in BM-MNCs from patients with CAD compared to control subjects (Figure 4A,C), as well as the basal and maximal ECAR (p<0.01 and p<0.05) (Figure 4B,D). The measurements were repeated 2 hours after IFN-γ/LPS pre-incubation, in order to assess metabolic function during activation. Again, the maximal OCR (p<0.01) and maximal ECAR (p<0.05) were higher in patients with CAD.
Proliferation of BM-MNCs
Proliferation assays of BM-MNCs did not reveal significant differences in erythroid-, myeloid-, or granulocyte-macrophage progenitor proliferation rates between patients with CAD and healthy controls (Figure 5).
Bone marrow progenitors of patients with CAD are primed to differentiate into an inflammatory myeloid lineage
To further understand the functional changes in HSPCs, we explored the transcriptional signature of HSC, MPP, and GMP cell populations. Differential expression analysis was applied to each population separately and a combined analysis of the populations to identify differentially expressed (DE) genes between patients and control samples. This identified 1747 genes that were differentially regulated in at least one of the cell populations (p≤0.05). A PCA analysis based on DE genes of HSC populations revealed a clear separation between patients and controls (Figure 6A). This was less pronounced for the MPPs and GMPs (Figure 6—figure supplement 1). Among those DE genes, we observed that four genes were significantly upregulated in the patients compared to the controls in a combined analysis of HSCs, MPPs, and GMPs (Figure 6B, padj <0.1), including CCR2, EPB42, FNDC3B, and RBMS1. For individual log fold change (FC) and adjusted p-values, please see Figure 6—source data 1 and the combined and separated heatmaps for the top 50 DE genes for each population (Figure 6—figure supplements 2 and 3). Seven genes were differentially downregulated in the patients, including PFKP, CCDC163P, ARMCX4, PTK7, WDR90, ROBO3, and FAM84B. Within the HSC population, only PROK2 showed a significant upregulation (Figure 6—figure supplement 3), whereas in the MPPs and GMPs no genes showed significantly differential upregulation (Figure 6—figure supplement 3).
Subsequently, we performed pathway enrichment analyses on disease-specific DE genes. The differentially upregulated genes of HSCs (p<0.05) showed significant enrichment for neutrophil and monocyte biological processes. These include neutrophil activation pathways, cytokine production pathways, and macrophage activation pathways (FDR <0.05) (Figure 6C). Transcriptional signatures of the genes enriched in the neutrophil-related pathways indicated upregulation of genes involved in the development and activation of neutrophils, such as insulin like growth factor two receptor (IGF2R), S100A11, and TNFRSF1B (Figure 6—source data 2). Additionally, there was an upregulation of signalling genes in the myeloid lineage, such as CCR2, innate immune signal transduction adaptor (MYD88), IL1RN, IL18R1, and Toll-like receptors TLR2 and TLR4. The differentially upregulated genes in the GMPs showed a similar enrichment for neutrophil-related pathways, and pathways related to myeloid cell differentiation and migration (Figure 6C). Upregulated genes in myeloid cell differentiation and regulation of hematopoiesis pathways revealed upregulation of important signalling and transcription factors, such as colony stimulating factor three receptor (CSF3R), NFKB inhibitor alpha (NFKBIA), signal transducer and activator of transcription 1 (STAT1), STAT3, and transforming growth factor beta 1 (TGFB1). Repeating these analyses using p<0.01 showed that the observed enriched terms were robust (FDR <0.05).
RNA-seq data have been deposited in the ArrayExpress database at EMBL-EBI (www.ebi.ac.uk/arrayexpress) under accession number E-MTAB-9399.
Vascular wall inflammation and hematopoietic tissue activation measured with [18F]FDG PET/CT is not higher in patients with CAD
Vascular wall inflammation and hematopoietic tissue activity, as determined by [18F]FDG PET/CT-scan, was similar between patients with CAD and individuals without atherosclerosis (Figure 7). Correction for age, BMI, and glucose concentrations did not influence the results.
Interestingly, although the splenic activity was not significantly higher in patients with CAD, the splenic [18F]FDG-uptake correlated with HSPCs and with circulating immune cells (Figure 7—figure supplement 1). Splenic [18F]FDG-uptake correlated positively with CCR2 expression on HSPCs (rs = 0.604, p<0.01), with GMPs (rs = 0.405, p<0.05), and with circulating leukocyte (rs = 0.393, p<0.05) and monocyte counts (rs = 0.588, p<0.01).
Discussion
The role of activated innate immune cells in the development of atherosclerotic plaques is well-established. In the current paper we add an important dimension to this pathophysiological framework by showing for the first time that inflammatory reprogramming of innate immune cells in humans with coronary artery disease occurs at the level of HSPCs in the bone marrow. In patients with chronic CAD due to severe coronary atherosclerosis, isolated bone marrow mononuclear cells showed a higher cytokine production capacity and increased metabolic rate compared to individuals without coronary atherosclerosis. Transcriptional signatures of HSPCs were enriched for neutrophil and monocyte related pathways, indicating that the HSPCs are primed to differentiate into an inflammatory myeloid lineage. These findings addinclude an additional layer of inflammatory regulation in atherogenesis which potentially offers novel targets for pharmacological strategies to prevent or treat atherosclerotic CVD.
Our findings align with previous studies in animal models in atherosclerosis that show increased myelopoiesis in specific conditions that promote CVD. Accumulation of cholesterol in HSPC increases proliferation and mobilization from the bone marrow and a differentiation bias towards the myeloid lineage, which is associated with accelerated atherosclerosis (Yvan-Charvet et al., 2010). In an animal model of diabetes, higher circulating numbers of inflammatory Ly6-Chi monocytes and neutrophils are present, which is due to increased myelopoiesis. In this model, the increased myelopoiesis is driven by the increased expression of S100A8 and S100A9 in neutrophils (Nagareddy et al., 2013). Interestingly, the expression of one of the S100 proteins (S100A11) was also increased in the HSCs in patients with CAD in our study. Chronic psychological stress and disturbed sleep have also been reported to increase myelopoiesis, which is associated with elevated circulating levels of neutrophils and Ly6Chigh monocytes and augmented atherosclerosis compared with control mice (Heidt et al., 2014; McAlpine et al., 2019).
Since we did not include patients with specific cardiovascular risk factors, but rather with established symptomatic coronary atherosclerosis, we cannot conclude which factors are responsible for the HSPC reprogramming in our patients. BMI and blood pressure were slightly higher in the patients compared to the controls. Cholesterol concentrations were lower at the moment of inclusion in our study, but might well have been higher previously, since the patients also used more statins. We do not have information about stress and sleep patterns. Another mechanism that might contribute to the HSPC functional reprogramming is trained immunity. Trained immunity describes the phenomenon that brief stimulation of innate immune cells leads to the development of a long-term hyperresponsive phenotype. This is mediated by profound intracellular metabolic and epigenetic reprogramming and occurs both at the level of mature circulating monocytes as well as at the level of their bone marrow progenitors (Netea et al., 2020). Christ et al have recently described that a 4-week period of Western-type diet feeding in atherosclerosis prone Ldlr-/-mice induced a functional, transcriptional and epigenetic reprogramming of circulating myeloid cells and their bone marrow progenitors, which persisted for at least for weeks after switching back to a chow diet (Christ et al., 2018). This was associated with increased circulating concentrations of inflammatory monocytes and granulocytes. Functional enrichment analysis of GMPs showed overrepresentation of TNF and Toll-like receptor signaling pathways, which align with our observations. In contrast, although monocytic signature genes were enriched, most of the granulocytic signature genes were downregulated, including S100A8 and S100A9, in GMPs isolated from WD-fed as compared to CD-fed mice (Christ et al., 2018). Patients with hypercholesterolemia also have circulating monocytes with augmented cytokine production capacity, increased glycolytic metabolism, and enrichment of activating histone modifications, which persist for three months despite cholesterol lowering with statin therapy. Circulating monocyte and neutrophil numbers, however, were not increased in these patients (Bekkering et al., 2019). In our current patient cohort, it remains to be established whether trained immunity is present since we did not assess epigenetic markers in the HSPCs.
The observation that absolute numbers of myeloid progenitors and circulating monocytes and neutrophils are not increased, despite functional and transcriptional inflammatory programming, fits with the effects of alternative inducers of trained immunity in humans. Vaccination with the Bacille Calmette-Guérin (BCG) vaccine induces trained immunity in humans in vivo (Arts et al., 2018). This is accompanied by an increased cytokine production capacity of bone marrow mononuclear cells, comparable to the result in patients with CAD in the current study (Cirovic et al., 2020). Also, 90 days post-BCG administration, HSPCs showed myeloid skewing with activation of neutrophil-associated pathways and pathways associated with regulation of immune responses, comparable to our patients. This functional and transcriptional reprogramming was accompanied by an increase in numbers of circulating myeloid cells in large cohorts of infants vaccinated with BCG (Cirovic et al., 2020).
In the patient group with CAD, the increased cytokine production capacity of the bone marrow mononuclear cell fraction was accompanied by an upregulation of both glycolytic metabolism as well as oxygen consumption by oxidative phosphorylation. It is important to realize that the bone marrow mononuclear fraction also contains large amounts of mature monocytes, so we cannot conclude on the metabolic landscape of HSPCs and specific progenitor subtypes. Isolated monocytes and monocyte-derived macrophages from patients with atherosclerosis are characterized by increased glycolysis, which is coupled to the hyperresponsiveness in terms of cytokine release (Bekkering et al., 2016; Shirai et al., 2016). Also, HSPCs of mice after induction of trained immunity with β-glucan have an upregulation of the glycolytic pathway (Mitroulis et al., 2018).
Analysis of the differentially upregulated genes transcriptome of the flow-sorted HSCs and GMPs showed significant enrichment for neutrophil activation pathways. This further supports a role for neutrophils in the pathophysiology of atherosclerotic CVD, which is suggested by various recent findings in the literature. Although neutrophils are the most abundant circulating white blood cell type, neutrophils only recently received attention in the context of cardiovascular inflammation. A large epidemiological study shows that circulating granulocyte counts are strongly associated with the occurrence of CVD (Fani et al., 2020). During atherogenesis, platelet-derived chemokines, such as CC-chemokine ligand 5, promote neutrophil activation and recruitment (Silvestre-Roig et al., 2020). At the luminal side, these activated neutrophils secrete granule proteins, including cathepsin G, which further promotes myeloid cell recruitment. Several actions further fuel atherosclerotic plaque formation, including secretion of reactive oxygen species and myeloperoxidase, which mediates oxidation of LDL, promoting foam cell formation (Silvestre-Roig et al., 2020). Also, neutrophils can form neutrophil extracellular traps (NETs) (Döring et al., 2020) and they can secrete pro-inflammatory microvesicles, which further fuel vascular wall inflammation and atherosclerosis (Gomez et al., 2020). Interestingly, a recent study used plasma proteomics analysis to explore the mechanism by which colchicine lowers cardiovascular risk in patients with CAD and showed mainly a downregulation of neutrophil activation pathways (Opstal et al., 2020).
Analysis of the individual significantly differentially expressed genes in the HSPCs showed that the glycolytic enzyme phosphofructokinase (PFKP) was downregulated in patients with CAD. Quiescent HSCs depend mainly on anaerobic glycolysis to sustain survival, quiescence, and retention within the bone marrow (Takubo et al., 2010). Once stimulated to divide, HSCs start to activate OXPHOS to meet metabolic demands of proliferation and differentiation (Karigane and Takubo, 2017). Therefore, the downregulation of PFKP might be related to the activation of HSCs. This finding is in contrast to the higher extracellular acidification rate of the bone marrow mononuclear cells which points to the activation of glycolysis. We speculate that this finding is dominated by the inflammatory activation of bone marrow monocytes, which constitute a large percentage of the bone marrow mononuclear cells. The chemokine receptor CCR2 showed significant upregulation in the HSPCs of the patients with CAD. CCR2 is critical for monocyte egress from the bone marrow and for the recruitment of circulating monocytes to the arterial wall (Tsou et al., 2007). The upregulation of CCR2 in our study might contribute to the increased circulating HSC percentage in patients with CAD. In mice, increased levels of circulating HSPCs are found after myocardial infarction that migrate to the spleen for extramedullary hematopoiesis (Dutta et al., 2012). In humans, splenic extramedullary hematopoiesis can occur in response to extreme physiologic stress (e.g. sepsis or hematologic malignancies), but has not been described yet during the development of CVD. Indirect evidence shows that splenic activity in humans is increased after myocardial infarction assessed with PET/CT scanning (van der Valk et al., 2016), suggesting increased immune cell activation. In our study, we observed a correlation between splenic activity and CCR2 expression on circulating HSPCs which might suggest that the spleen performs a role in human atherosclerosis.
Our study has several limitations. First, the sample size is relatively small. However, this sample size needs to be considered in the context of the invasive and complex nature of the study, and this sample size usually allows to detect relevant differences (Cirovic et al., 2020; van der Heijden et al., 2020). The sample size might also be an explanation for the lack of a significant difference in FDG uptake in the large arteries, which has been shown previously in patients with CVD. (Tarkin et al., 2014; van Wijk et al., 2014). Another possible explanation for our results is that our patients had severe mature coronary atherosclerotic plaques with calcium deposits, which are known to have less FDG uptake than uncalcified arterial plaques at the early stage of atherosclerosis (Fernández-Friera et al., 2019). In addition, lipid-lowering therapy is effective in reducing vascular wall inflammation (Pirro et al., 2019). Second, we did not explore epigenetic programs of the HSPCs. As such, we cannot conclude whether the mechanism of trained immunity contributes to our findings. This should be the focus of future studies. Third, in a study on symptomatic patients with CVD, it is inevitable to have differences in medication use compared to the healthy control group.
To summarize, our study is the first to prospectively assess HSPC phenotype in patients with CAD and is unique in providing a comprehensive phenotype by the integration of powerful tools such as RNAseq on flow-sorted HSPC populations, functional parameters, Seahorse respiratory, and PET imaging of vascular wall and hematopoietic activity. These results show that in patients with CAD, activation of the innate immune system occurs at the level of bone marrow myeloid progenitors, which adds exciting opportunities for novel treatment strategies.
Materials and methods
Participant selection
Request a detailed protocolParticipants (male, 18–75 years) were recruited among patients who were admitted for the evaluation of chest pain and underwent cardiac imaging at the Cardiology department of the Canisius Wilhelmina hospital or the Radboud University Medical Centre, Nijmegen, The Netherlands after January 1, 2015. Severe coronary atherosclerosis was defined as a calcium score >400 HU on computed tomography (CT), or a total plaque score (TPS) >4 on coronary CT angiography (CCTA), according to previously described standards (Bekkering et al., 2016; Pen et al., 2013). Control participants, with a calcium score and TPS of zero, were matched for age, body mass index (BMI), and smoking.
Criteria for exclusion were previous cardiovascular events, malignancies, auto-immune or auto-inflammatory diseases (including diabetes mellitus), chronic immunomodulatory drug use, chronic kidney disease (MDRD <45 mL/min), liver disease (ALAT > 135 U/L), or thrombocytopenia (<50 × 106/mL). Additionally, participants were excluded if they had an infection (>38,5°C or antibiotic treatment), hospital admission, or vaccination within 3 months before study entry. The study protocol was approved by the Institutional Review Board Arnhem/Nijmegen, the Netherlands, and registered at ClinicalTrials.gov (NCT03172507). All individuals gave written informed consent.
Evaluation of atherosclerotic burden
Request a detailed protocolCCTA total plaque score was calculated as previously described (Pen et al., 2013). Briefly, a 64-slice MDCT scanner (Philips) was used to obtain CAC and MDCT Image acquisition using an ECG-synchronized axial scan protocol and post-processing CT and CAC studies using IntelliSpace Philips software. Before image acquisition, beta-blockers were administered targeting a heart rate of <60 beats per minute, and patients received nitroglycerin 0.8 mg sublingually. Prospective electrocardiographically gated step-and-shoot contrast-enhanced MDCT imaging was performed, initiated from 10 mm above the level of the left main artery to 10 mm below the inferior myocardial apex with scan parameters being 64 × 0.625 mm sections (2.5 mm), collimation tube currents of 350 to 780 mAs and tube voltage of 100 or 120 kV. In the rare event that prospective scanning was not possible, retrospective or helix scanning was used.
Reconstruction of the MDCT scans was performed with reconstructed images obtained, using an ECG-triggered protocol, at 75% from the previous RR-interval, or at 75% and 40% from the previous RR-interval if a helix scan-protocol was used. The TPS was determined by summing the number of evaluable coronary segments with calcific or non-calcific plaque, or mixed plaque, where non-calcified and mixed plaque was assigned with one point and calcified plaque with 0 (maximum score = 16). Two independent experienced operators scored all CT images and both operators were blinded for all clinical information. In case of disagreement, the opinion of a third independent observer was asked.
Study design
Request a detailed protocolParticipants underwent [18F]FDG PET/CT scanning, follow by venepuncture and bone marrow aspiration within 2–14 days. The participants were invited in pairs (1 patient:1 control) for blood sampling and bone marrow aspiration from December 2017 till July 2018.
Cardiovascular risk assessment
Request a detailed protocolMedical history, smoking status, medication use, BMI, and fasting glucose concentrations were obtained from all individuals. Blood pressure was measured three times by a manual sphygmomanometer after 5 min seated rest according to AHA guidelines. Fasting total cholesterol (Tchol), high-density lipoprotein cholesterol (HDLc), and triglycerides (TG) were measured using standardized methods, and low-density lipoprotein cholesterol (LDLc) was calculated with the Friedewald formula.
Blood sampling and bone marrow aspiration
Request a detailed protocolBone marrow was aspirated from the posterior iliac crest according to standard practice by an experienced physician assistant. Blood was sampled through venous puncture. Sample collection was performed at 8.00–10.00 to avoid interference of circadian rhythms of immune parameters, and sample processing occurred within 2 hours. Plasma and serum were stored at −80°C until further use. Total blood cell counts were determined with an automated Sysmex-XN 450 hematology analyzer (Sysmex, Hamburg, Germany).
Mononuclear cell enrichment and stimulation
Request a detailed protocolBefore mononuclear cell (MNC) enrichment, the bone marrow aspirate was filtered and washed with sterile PBS. Thereafter, the same procedures were followed for peripheral blood MNCs (PBMCs) and bone marrow MNCs (BM-MNCs). PBMCs/BM-MNCs were isolated by Ficoll-Paque density gradient centrifugation (GE Healthcare, Chicago, IL). Cell composition was evaluated by Sysmex analyzer (Sysmex) and with flow cytometry (Table 3, see key resource table (KRT) for RRIDs). PBMCs/BM-MNCs were resuspended in Roswell Park Memorial Institute 1640 Dutch-modified culture medium (RPMI) (Life Technologies/Invitrogen, Waltham, USA) supplemented with 2 mmol/L glutamine (Invitrogen), 10 mg/mL gentamicin (Centrafarm, Etten-Leur, The Netherlands) and 1 mmol/L pyruvate (Invitrogen). Per well, 5 × 105 PBMCs/BM-MNCs were stimulated for 24 hours in duplicate in round-bottom 96-well plates (Corning, NY) with the following stimuli: RPMI, 10 ng/mL Escherichia coli lipopolysaccharide (LPS) (serotype 055:B5 Sigma-Aldrich, St. Louis, MO), and 10 μg/mL Pam3CysK4 (P3C) (EMC Microcollections, Tübingen, Germany). After 24-hour incubation, supernatants were stored after plate centrifugation at −80°C until cytokine assessment.
Cytokine measurements
Request a detailed protocolCytokine and chemokine concentrations were determined in supernatants using Human Cytokine Magnetic Magpix 25-plex panel (Invitrogen) on the MAGPIX platform (Luminex, Austin, TX). Circulating IL-1β, IL-1Ra, IL-6, and IL-18 concentrations were measured with the SimplePlex cartridge on the Ella platform (ProteinSimple, San Jose, CA). Additional circulating cytokines/chemokines concentrations were measured using ELISA (see KRT).
Flow cytometry
Request a detailed protocolIn the circulation, monocyte subpopulations and expression markers were determined with flow cytometry. 50 µL EDTA blood was stained after the lysis-no-wash strategy (BD Pharm Lyse lysing buffer, Becton Dickinson) by monoclonal antibodies CD45 Krome Orange ([KO], clone J33; Beckman Coulter, Cat# B36294, RRID:AB_2833027), HLA-DR PE (clone immu-357; Beckman Coulter, Cat# IM1639U, RRID:AB_2876782), CD14 PC7 (clone 61D3e Bioscience, Cat# 25-0149-42, RRID:AB_1582276), CD16 FITC (clone CB16; eBioscience Cat# 11-0168-42, RRID:AB_10805747), CD3 APC-Alexa750 (clone UCTH1; Beckman Coulter, Cat# AA66329, RRID:AB_2876783), CD56 APC (clone N901; Beckman Coulter, Cat# IM2474U, RRID:AB_2876784), CD192 Brilliant Violet421 ([BV421] clone 48607; Becton Dickinson, Cat#564067, RRID:AB_2738573), CD11b BV785 (clone ICRF44; Biolegend, Cat#301346, RRID:AB_2563794), CD41 PC5.5 (clone Hip8; Biolegend, Cat# 303719, RRID:AB_2561731) and measured with CytoFLEX flow cytometer (Beckman Coulter, RRID:SCR_017217). The gating strategy applied is shown in Supplementary file 1, gates were set with the fluorescence-minus-one method (Weber et al., 2016; Ziegler-Heitbrock et al., 2010). In short, monocytes were selected based on CD45+ HLA-DR+ and monocyte scatter properties, then CD3+ T-lymphocytes and CD56+ NK-cells were excluded, and monocyte subsets were identified in the CD14/CD16 plot as percentage of gated. Data was analyzed with Kaluza 3.1 software (Beckman Coulter, RRID:SCR_016182). Characterization of monocytes subsets is according to current recommendations (Weber et al., 2016; Ziegler-Heitbrock et al., 2010).
Bone marrow progenitors were identified with Navios flow cytometer (Beckman Coulter, RRID:SCR_014421). 5 × 106 Bone marrow cells were washed, lysed, and stained for 20 min with the monoclonal antibodies: CD90 FITC (Clone 5E10; Biolegend Cat#328107, RRID:AB_893438), CD123 PE (Clone 9F5; BD Biosciences, Cat#555644, RRID:AB_396001), CD19 ECD (Clone J3.119; Beckman Coulter, Cat# IM2708U, RRID:AB_130854), CD38 PC5.5 (Clone LS198-4-3; Beckman Coulter, Cat# IM2651U, RRID:AB_131166), CD117 PEC7 (Clone 104D2D1; Beckman Coulter, Cat# IM3698, RRID:AB_131184), CD45RA APC (Clone 2H4LDH11LD89; Beckman Coulter, Cat# B14807, RRID:AB_2876787), DRAQ7 (Biostatus), CD34-APC A750 (Clone 581; Beckman Coulter, Cat# A89309, RRID:AB_2876786), CD192 PB (Clone 48607; Becton Dickinson, Cat# 564067, RRID:AB_2738573), CD45 KO (Clone J33; Beckman Coulter, Cat# B36294, RRID:AB_2833027). The target input was set at 1 × 106 measured events, panels were balanced with the fluorescence-minus-one method to determine spill over, spectral overlap, and nonspecific binding.
The gating strategy is displayed in Figure 2—figure supplement 1. In short, HSPCs were defined as CD45+CD34+CD38dim/+ cells, after selecting for singlets and alive cells. Next, the lymphoid lineage was excluded by the gating of CD19-CD117+ cells. In CD45RAdimCD38+ gated cells, CMP, GMP, MEP, and R1-3 progenitor populations were identified using CD123 and CD45RA expression. CD90 expression in CD38-CD45RA- cells determined MPP and HSC populations. Data were analyzed with Kaluza 2.1 software.
Proliferation assays
Request a detailed protocol2.5 × 104 BM-MNCs were cultured in methylcellulose medium containing erythropoietin, G(M)-CSF, IL-3 and IL-6 (MethoCult GF H84435, Stemcell technologies, Vancouver, Canada) and 2% fetal bovine serum (Integro B.V., Zaandam, The Netherlands) in Petri dishes for 14 days in 37°C 5% CO2 incubator. After 14 days, erythroid progenitor (BFU-E), common myeloid progenitor (CFU-GEMM), and granulocyte-macrophage progenitor (CFU-GM) cell colonies were counted in duplo by an experienced hematology operator.
Mitochondrial respiration and glycolysis assays
Request a detailed protocolUsing a modified protocol from Shirai et al., 2016, 2 × 106 BM-MNCs were plated per well in fivefold to overnight-calibrated cartridges in assay medium (DMEM with 1 mM L-Glutamine, 11 mM glucose, and 1 mM pyruvate [pH adjusted to 7.4]) and incubated for 1 hr in a non-CO2-corrected incubator at 37°C. In addition, BM-MNCs were stimulated with 50 ng/mL interferon gamma (IFN-γ) (Immukine, Boehringer Ingelheim BV, Alkmaar, The Netherlands) and 10 ng/mL LPS (Sigma-Aldrich) for 2 hr. Oxygen consumption rate (OCR) and extracellular acidification rate (ECAR) were measured via XFp Analyzer (Seahorse Bioscience, North Billerica, MA), with final concentrations of 1 mM oligomycin, 1 mM FCCP, and the combination of 1.25 mM rotenone and 2.5 mM antimycin A.
RNA sequencing of sorted progenitor populations and circulating monocytes
Request a detailed protocolHSPCs were sorted by BD FACSAria II SORP flow cytometer sorter (Becton Dickinson, RRID:SCR_018091), using the flow cytometric panel as described above (Figure 2—figure supplement 1). BM-MNCs stored in liquid nitrogen were gently thawed in fetal calf serum (in house) containing 1.25 mM MgCl2 and 0.1 mg/mL DNase I (Sigma). HSCs, MPPs, and GMPs populations were directly sorted in 100 µL TRIzol (ThermoFisher, Waltham, MA) before processing for RNA sequencing.
Total RNA was extracted using the standard TRIzol (ThermoFisher) protocol and used for library preparation and sequencing. mRNA was processed as described previously, following an adapted version of the single-cell mRNA seq protocol of CEL-Seq (Hashimshony et al., 2012; Simmini et al., 2014). In brief, samples were barcoded with CEL-seq primers during reverse transcription and pooled after second strand synthesis. The resulting cDNA was amplified with an overnight in vitro transcription reaction. From this amplified RNA, sequencing libraries were prepared with Illumina Truseq small RNA primers (Illumina, San Diego, CA). Paired-end sequencing was performed on the Illumina Nextseq500 platform (RRID:SCR_014983). Read one was used to identify the Illumina library index and CEL-Seq sample barcode. Read two was aligned to the hg19 human RefSeq transcriptome using BWA (Li and Durbin, 2010). Reads that mapped equally well to multiple locations were discarded. RNA input for all samples was normalized, and libraries for each progenitor population were sequenced in a single run.
Bioinformatics, differential gene expression, and pathway analysis
Request a detailed protocolReads were mapped to hg19 human reference genome using BWA (Li and Durbin, 2010). MultiQC (RRID:SCR_014982) was used to quality check all the samples (Ewels et al., 2016). In total, 107,565,838 reads were mapped in 59 progenitor populations, with one GMP population failed to pass quality control. Raw read counts provided the input for differential gene expression analysis with DESeq2 (RRID:SCR_015687) (Love et al., 2014) after excluding genes with less than 50 reads. Internal variance stabilizing transformation was used for regressing out age and to normalize the data (Benjamini–Hochberg correction was used to adjust p-value for multiple testing correction). In the differential expression analyses within HSC populations, one HSC population was marked as an outlier and further filtered out as it located three standard deviations away from mean values in principle component analysis of this study.
Differential gene expression analysis was carried out after regressing out age with DESeq2 v1.22.0 in R (RRID:SCR_001905) (Love et al., 2014), with internal statistical and normalization method (i.e. correction for multiple testing with Benjamini–Hochberg). The average expression between individuals with atherosclerosis (n = 10) and matched individuals without atherosclerosis (n = 10) was analyzed for each progenitor population. Enrichment analysis was performed using R package clusterProfiler (RRID:SCR_016884) (Yu et al., 2012) with Gene Ontology terms (Ashburner et al., 2000). After BH adjustment, a FDR <0.05 in enrichment analysis was considered significant.
[18F]FDG PET and low-dose CT scanning
Request a detailed protocolParticipants underwent [18F]FDG PET with low-dose non-contrast-enhanced CT from skull base to the trochanter major on a dedicated Siemens Biograph 40 mCT scanner (Siemens Healthineers, Erlangen, Germany). After adhering to a 24-hour low-carbohydrate diet and 6 hours of fasting,~2.1 MBq/kg [18F]FDG was administered intravenously, as described previously (Bucerius et al., 2016). Glucose concentrations were obtained (5.3 ± 0.5 mmol/L) after injection. Before scanning, participants rested in the supine position for 2 hours.
Images were reconstructed according to EARL protocols; using a TrueX algorithm with point spread function (PSF) and time-of-flight (TOF) measurements, using three iterations, 21 subsets, matrix size 200 × 200 (pixel spacing of 4.07 mm), full width half maximum (FWHM) of 3 mm and using 2 min of PET data. Postprocessing was performed using a 3D Gaussian filter kernel, 3.0 mm, using the Inveon Research Workspace 4.2 (Preclinical Solutions, Siemens Medical Solutions USA, Knoxville, TN).
[18F]FDG-uptake in the vascular wall was determined in seven regions of interest (ROI) by a single investigator (MPN); the aorta ascendens, aorta descendens, abdominal aorta, the left and right common carotid arteries, and the left and right iliac arteries. [18F]FDG-uptake in hematopoietic tissue was assessed in the spleen, lumbal vertebrae L2 and L3, and in the left and right medullary bone of the femur. These regions of interest were evaluated using the Inveon Research Workspace 4.2. The standardized uptake value (SUV) was extracted from each ROI after correction for [18F]FDG dose (MBq) and BMI using the PyRadiomics toolbox (van Griethuysen et al., 2017). The SUVs of left and right ROIs (e.g. left and right carotid artery) were averaged. Next, the target-to-background ratio (TBR) was calculated as the ratio of the vascular wall SUV and the mean thoracic arterial blood pool SUV. The TBRs of hematopoietic tissue were expressed as the ratio of the mean liver SUV. The primary outcome is the TBR, as recommended by the European guideline (Bucerius et al., 2016).
Statistical analysis
Request a detailed protocolThis study is exploratory, hence no sample size calculation is performed. Normal distribution of the data was checked with the Shapiro-Wilk test, when the p-value reached <0.05 this assumption was violated and non-parametric tests were used. Data are reported as mean ± SD with independent samples T-test according to Levene’s test for equality of variances, as mean (number of participants) with Χ2 test for categorical data, and as median [interquartile range] with Mann-Whitney U test for non-parametrical data. Outliers were removed with a standard deviation >± 2.5 of Z-scores. All outcomes were log(10)-transformed and thereafter corrected with ANCOVA for confounding demographics such as age and BMI. SPSS V25.0 (SPSS Inc, Chicago, IL, RRID:SCR_002865) and Graphpad Prism v6.0 (GraphPad software, La Jolla, CA, RRID:SCR_002798) were used for data analysis and visualization. A two-sided p-value<0.05 was considered statistically significant.
Data availability
RNA-seq data have been deposited in the ArrayExpress database at EMBL-EBI (https://www.ebi.ac.uk/arrayexpress) under accession number E-MTAB-9399.
-
ArrayExpressID E-MTAB-9399. Bone marrow progenitor cell populations of patients with severe coronary artery disease.
References
-
Position paper of the cardiovascular committee of the european association of nuclear medicine (EANM) on PET imaging of atherosclerosisEuropean Journal of Nuclear Medicine and Molecular Imaging 43:780–792.https://doi.org/10.1007/s00259-015-3259-3
-
Vascular inflammation in subclinical atherosclerosis detected by hybrid PET/MRIJournal of the American College of Cardiology 73:1371–1382.https://doi.org/10.1016/j.jacc.2018.12.075
-
Chronic variable stress activates hematopoietic stem cellsNature Medicine 20:754–758.https://doi.org/10.1038/nm.3589
-
Metabolic regulation of hematopoietic and leukemic stem/progenitor cells under homeostatic and stress conditionsInternational Journal of Hematology 106:18–26.https://doi.org/10.1007/s12185-017-2261-x
-
Macrophages in atherosclerosis: a dynamic balanceNature Reviews Immunology 13:709–721.https://doi.org/10.1038/nri3520
-
Defining trained immunity and its role in health and diseaseNature Reviews Immunology 20:375–388.https://doi.org/10.1038/s41577-020-0285-6
-
Colchicine in patients with chronic coronary diseaseNew England Journal of Medicine 31:2021372.https://doi.org/10.1056/NEJMoa2021372
-
Discordance between framingham risk score and atherosclerotic plaque burdenEuropean Heart Journal 34:1075–1082.https://doi.org/10.1093/eurheartj/ehs473
-
Antiinflammatory therapy with canakinumab for atherosclerotic diseaseNew England Journal of Medicine 377:1119–1131.https://doi.org/10.1056/NEJMoa1707914
-
The glycolytic enzyme PKM2 bridges metabolic and inflammatory dysfunction in coronary artery diseaseJournal of Experimental Medicine 213:337–354.https://doi.org/10.1084/jem.20150900
-
Neutrophils as regulators of cardiovascular inflammationNature Reviews Cardiology 17:327–340.https://doi.org/10.1038/s41569-019-0326-7
-
Efficacy and safety of Low-Dose colchicine after myocardial infarctionNew England Journal of Medicine 381:2497–2505.https://doi.org/10.1056/NEJMoa1912388
-
PET imaging of inflammation in atherosclerosisNature Reviews Cardiology 11:443–457.https://doi.org/10.1038/nrcardio.2014.80
-
Critical roles for CCR2 and MCP-3 in monocyte mobilization from bone marrow and recruitment to inflammatory sitesJournal of Clinical Investigation 117:902–909.https://doi.org/10.1172/JCI29919
-
Catecholamines induce trained immunity in monocytes in vitro and in vivoCirculation Research 127:269–283.https://doi.org/10.1161/CIRCRESAHA.119.315800
-
Increased haematopoietic activity in patients with atherosclerosisEuropean Heart Journal 15:425–432.https://doi.org/10.1093/eurheartj/ehw246
-
Computational radiomics system to decode the radiographic phenotypeCancer Research 77:e104–e107.https://doi.org/10.1158/0008-5472.CAN-17-0339
-
Nonpharmacological lipoprotein apheresis reduces arterial inflammation in familial hypercholesterolemiaJournal of the American College of Cardiology 64:1418–1426.https://doi.org/10.1016/j.jacc.2014.01.088
-
clusterProfiler: an R package for comparing biological themes among gene clustersOMICS: A Journal of Integrative Biology 16:284–287.https://doi.org/10.1089/omi.2011.0118
Article and author information
Author details
Funding
Horizon 2020 (667837)
- Leo AB Joosten
- Mihai G Netea
- Niels P Riksen
Netherlands Organisation for Scientific Research (NWO SPI 94-212)
- Mihai G Netea
European Commission (833247)
- Mihai G Netea
ERA-NET (2018T093)
- Niels P Riksen
Netherlands Organisation for Scientific Research (452173113)
- Siroon Bekkering
Hartstichting (2018T028)
- Siroon Bekkering
Hartstichting (CVON2018-27)
- Leo AB Joosten
- Mihai G Netea
- Niels Peter Riksen
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We would like to thank Cor Jacobs for performing the flow cytometric measurements, and Rob Woestenenk for his expertise with flow sorting at the Radboud technology center flow cytometry.
Ethics
Human subjects: Informed consent was obtained for all participants. The study protocol was approved by the Institutional Review Board Arnhem/Nijmegen, the Netherlands and registered at the ClinicalTrials.gov (NCT03172507).
Copyright
© 2020, Noz et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,043
- views
-
- 280
- downloads
-
- 31
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Medicine
Glucagon-like peptide 1 (GLP-1) is a gut-derived hormone secreted by intestinal L cells and vital for postprandial glycemic control. As open-type enteroendocrine cells, whether L cells can sense mechanical stimuli caused by chyme and thus regulate GLP-1 synthesis and secretion is unexplored. Molecular biology techniques revealed the expression of Piezo1 in intestinal L cells. Its level varied in different energy status and correlates with blood glucose and GLP-1 levels. Mice with L cell-specific loss of Piezo1 (Piezo1 IntL-CKO) exhibited impaired glucose tolerance, increased body weight, reduced GLP-1 production and decreased CaMKKβ/CaMKIV-mTORC1 signaling pathway under normal chow diet or high-fat diet. Activation of the intestinal Piezo1 by its agonist Yoda1 or intestinal bead implantation increased the synthesis and secretion of GLP-1, thus alleviated glucose intolerance in diet-induced-diabetic mice. Overexpression of Piezo1, Yoda1 treatment or stretching stimulated GLP-1 production and CaMKKβ/CaMKIV-mTORC1 signaling pathway, which could be abolished by knockdown or blockage of Piezo1 in primary cultured mouse L cells and STC-1 cells. These experimental results suggest a previously unknown regulatory mechanism for GLP-1 production in L cells, which could offer new insights into diabetes treatments.
-
- Medicine
The zinc transporter Slc30a1 plays an essential role in maintaining cellular zinc homeostasis. Despite this, its functional role in macrophages remains largely unknown. Here, we examine the function of Slc30a1 in host defense using mice models infected with an attenuated stain of Salmonella enterica Typhimurium and primary macrophages infected with the attenuated Salmonella. Bulk transcriptome sequencing in primary macrophages identifies Slc30a1 as a candidate in response to Salmonella infection. Whole-mount immunofluorescence and confocal microscopy imaging of primary macrophage and spleen from Salmonella-infected Slc30a1flag-EGFP mice demonstrate Slc30a1 expression is increased in infected macrophages with localization at the plasma membrane and in the cytosol. Lyz2-Cre-driven Slc30a1 conditional knockout mice (Slc30a1fl/fl;Lyz2-Cre) exhibit increased susceptibility to Salmonella infection compared to control littermates. We demonstrate that Slc30a1-deficient macrophages are defective in intracellular killing, which correlated with reduced activation of nuclear factor kappa B and reduction in nitric oxide (NO) production. Notably, the model exhibits intracellular zinc accumulation, demonstrating that Slc30a1 is required for zinc export. We thus conclude that zinc export enables the efficient NO-mediated antibacterial activity of macrophages to control invading Salmonella.