- Research article
- Open Access
Monocyte transcriptomes from patients with axial spondyloarthritis reveal dysregulated monocytopoiesis and a distinct inflammatory imprint
Arthritis Research & Therapy volume 23, Article number: 246 (2021)
In patients with axial spondyloarthritis (axSpA), monocytes show a pre-activated phenotype. Gut inflammation is a trigger of monocyte activation and may also affect their development in the bone marrow (BM). As gut inflammation is commonly observed in axSpA patients, we performed a detailed analysis of monocyte transcriptomes of axSpA patients in two cohorts and searched for signs of activation and developmental adaptations as putative imprints of gut inflammation.
Transcriptomes of blood CD14+ monocytes of HLA-B27+ axSpA patients and healthy controls (HC) were generated by microarrays from cohort 1 and by RNA-sequencing from cohort 2. Differentially expressed genes from both analyses were subjected to gene set enrichment analysis (GSEA) and to co-expression analysis in reference transcriptomes from BM cells, blood cells and activated monocytes. As serological markers of translocation, 1,3 beta-glycan, intestinal fatty acid binding protein, and lipopolysaccharide binding protein (LBP) were determined by LAL and ELISA.
Transcriptome analysis identified axSpA-specific monocyte signatures showing an imprint of LPS/cytokine-activated monocytes, late granulopoietic BM cells, blood neutrophils, and G-CSF-mobilized blood cells, which suggests LPS/TNF activation and more prominent BM adaptation promoting a neutrophil-like phenotype. GSEA mapped axSpA upregulated genes to inflammatory responses and TNFα signaling and downregulated probe-sets to metabolic pathways. Among translocation markers, LBP levels were significantly increased in axSpA patients vs. HC (p < 0.001).
Stratified analysis by disease activity and stage identified an “active disease signature” (BASDAI ≥ 4) with an imprint of LPS/cytokine-activated monocytes and CD16+ monocyte subsets. The “AS signature” (vs. non-radiographic axSpA) showed a reinforced neutrophil-like phenotype due to deprivation of dendritic cell transcripts.
The neutrophil-like phenotype of axSpA monocytes points towards a biased monocytopoiesis from granulocyte-monocyte progenitors. This shift in monocytopoiesis and the LPS/cytokine imprint as well as the elevated LBP levels are indicators of systemic inflammation, which may result from bacterial translocation. The BM adaptation is most prominent in AS patients while disease activity appears to be linked to activation and trafficking of monocytes.
Axial spondyloarthritis (axSpA), comprising non-radiographic axial SpA (nr-axSpA) and ankylosing spondylitis (AS), is characterized by inflammation within the sacroiliac joints and the spine.
Although disease development is strongly linked to the presence of the MHC class I molecule HLA-B27 , innate immune cells appear particularly important for disease development. Thus, genome wide association studies identified multiple risk genes for AS that are related to innate immune functions such as autophagy, NFkB regulation, and myeloid cell function [2,3,4]. Moreover, TNFα, which is successfully targeted in axSpA, is predominantly produced by innate immune cells, specifically monocytes .
First, Wright et al. reported inflammatory changes in monocytes of AS patients . They found an increase of proinflammatory proteins and an activation of the ubiquitin proteasome pathway in ex vivo isolated AS monocytes by mass spectrometry. In our own previous study in axSpA patients, we observed signs of preactivation and changes in the monocyte subsets . Thus, we found a higher percentage of classical, i.e., CD14+CD16− monocytes as compared to healthy controls. The classical monocytes comprise more than 90% of blood monocytes  and give rise to the less abundant subsets of intermediate (CD14+CD16+) and non-classical monocytes (CD14−CD16+) . As the percentage of individual populations is dynamically changed upon inflammatory triggers such as endotoxemia [9, 10] the change in monocyte subsets in axSpA may be related to monocyte request and activation [9, 10].
In this line, monocytes of axSpA patients showed a pre-activated phenotype as indicated by an increased spontaneous and induced production of proinflammatory cytokines, including TNFα .
A recent study showed that circulating microbial antigens does not only activate monocytes directly but also affect their progenitors in the bone marrow (BM). This results in a shift in progenitor use during monocytopoiesis, in addition to enhanced myeloid cell output . Thus, under homeostatic conditions, classical monocytes comprise a mixture of cells that are derived either from granulocyte-monocyte progenitors (GMP) or monocyte-dendritic cell progenitors (MDP) . GMP-derived monocytes share transcriptional pattern with neutrophils, while MDP-derived monocytes share transcriptional signatures with DCs. It was shown that circulating LPS promotes GMP-dependent generation of neutrophil-like monocytes while CpG favors MDP-dependent monocyte production and by that shifts the ratio between neutrophil-like monocytes and DC-like monocytes .
To search for signs of BM adaptation and transcriptional preactivation of monocytes in axSpA patients, we assessed monocyte transcriptomes of axSpA patients. In addition, we determined levels of intestinal fatty acid binding protein (I-FABP), (1→3)-β-D-Glucan (βDG) and lipopolysaccharide binding protein (LBP) as putative serological markers of bacterial translocation.
Monocyte transcriptomes were generated from two independent cohorts by either GeneChip microarrays (cohort 1) or RNA-sequencing (RNA-seq — cohort 2). Cohort 1 included 25 axSpA patients and 10 healthy controls (HC). Cohort 2 included 32 axSpA patients and 22 HC. All axSpA patients were HLA-B27+ and fulfilled the ASAS criteria of axSpA . AS was defined according to the modified New York criteria . Disease activity was determined by BASDAI  and ASDAS . Patients in cohort 1 were treated with non-steroidal anti-inflammatory drugs (NSAIDs) on demand or continuously but not by biological DMARDs; in cohort 2, 34% (11/32 patients) were treated by biological DMARDs. Patient characteristics are given in Table 1.
Analysis of serological translocation markers was performed in cohort 1 and in 59 axSpA patients from a third cohort, the Gent Inflammatory Arthritis and spoNdylitis cohort (GIANT), in which microscopic evaluation of gut biopsies taken upon colonoscopy was performed . None of GIANT patients were treated with biological DMARDs. Gut biopsies were classified according to microscopic evaluation as “no colonic inflammation”, “acute” (neutrophilic) and “chronic” (lymphocytic) inflammation.
All participants gave written consent to the study, which was approved by the local ethical committees of the university hospitals Charité (Berlin, Germany), Ile-de-France XI (Saint-Germaine-en-Laye, France), and Ghent University Hospital (Ghent, Belgium).
Measurement of translocation markers in the serum
Serum was collected and stored until measurement at least at − 20 to – 80 °C. Commercial ELISA or limulus based detection kits were used for measurement of I-FABP (Hycult Biotech), (1→3)-β-D-Glucan (Fujifilm Wako), and LBP (USCN).
Monocyte purification, transcriptomic profiling, and statistical analysis
Extended methodological details and details on statistical analysis are given in a supplementary file (Supplementary data 1). In brief, CD14+ monocytes were isolated from heparinized whole blood by positive magnetic selection using CD14 microbeads. In cohort 1, RNA isolation and gene chip hybridization were performed as previously described . In cohort 2, monocyte transcriptomes were generated by RNA-sequencing using the Illumina HiSeq4000 sequencer. For GeneChip microarray analysis, significant differences in expression were identified using high performance chip data analysis (HPCDA) (for details see Supplementary data 1). A HPCDA score of > 100 was considered significant. For RNA-sequencing, analysis of differential gene expression was carried out with edgeR using a linear model including the disease status, the gender and the sequencing experiment. Scores reaching a p-value < 0.05 were considered significant.
Hierarchical clustering of transcriptomes was performed with Genes@Work software . Principal component analysis (PCA) was performed using Qlucore (Lund, Sweden). For functional allocation, gene set enrichment analysis (GSEA) was performed using the Molecular Signature Database (MSigDB) v7.2 [23, 24]. Differentially expressed probe-sets and genes were analyzed for co-expression in 70 reference transcriptomes (all Affymetrix HG-U133 Plus 2.0 transcriptomes) generated by our own or retrieved from Gene expression omnibus (GEO) data repository as previously described . In brief, data sets were harmonized by quantile normalization before application in co-expression analysis. Pearson correlation coefficients were calculated between the differentially expressed probe-sets of individual comparisons on the basis of the 70 reference transcriptomes. Euclidian distance and average linkage as an agglomeration rule were applied for hierarchical clustering of this gene-to-gene correlation. For analysis of cohort 2, differentially expressed genes were matched to GeneChip probe-sets before analysis.
The axSpA-specific monocyte signature shows alterations related to inflammation and metabolism
For transcriptome analysis, blood monocytes were isolated by positive magnetic selection using CD14 microbeads. This selection method resulted in ≥ 95% purity and induces minor transcriptional changes . First, transcriptomes from monocytes of 25 HLA-B27+ axSpA patients and 10 HC from cohort 1 were generated by microarray. Using pair-wise comparisons and a HPCDA score > 100 as cut-off, 957 probe-sets referring to 805 genes were found differentially expressed between axSpA patients and HC (Fig. 1A, Supplementary Table 1). Of these 957 probe-sets, 661 probe-sets were up- and 296 probe-sets were downregulated. Although differences in expression of individual probe-sets were low (mean FC 1.20), both PCA and hierarchical clustering demonstrated clear separation of monocyte transcriptomes of axSpA patients from HC (Fig. 1B, C). For functional classification of differentially expressed genes, the MSigDB as integral part of GSEA was applied. This analysis mapped axSpA upregulated probe-sets with functional groups such as “protein secretion”, “TNFα signaling via NFkB”, “inflammatory responses”, “apoptosis”, and “interferon gamma responses” suggesting a proinflammatory activation of the axSpA monocytes. For axSpA-downregulated genes GSEA showed overlap with pathways such as “oxidative phosphorylation”, “hypoxia”, and “G2M checkpoint” mostly indicative of metabolic changes (Fig. 1D). A few of the downregulated genes also mapped to the “TNF signaling via NFkB” functional group. Most of these genes including IFNGR2, ETS2, and FOS rather downregulate proinflammatory responses in monocytes/ macrophages [27, 28] suggesting profound alteration of TNF signaling pathway in monocytes in axSpA patients.
In cohort 2, monocyte transcriptomes were generated by RNA-seq. Group-wise comparison of transcriptomes of HLA-B27+ axSpA patients and HLA-B27- HC identified 48 genes as up- and 47 as downregulated (corresponding to 109 and 129 GeneChip probe-sets, respectively,) in axSpA patients compared to HC (p-value < 0.05; Supplementary Table 1). The lower number of differentially expressed genes is most likely related to more stringent analysis (HPCDA scoring compared to edgeR computed p-values). Furthermore, the overall overlap in significant genes between cohort 1 and 2 was distressingly small (only one gene of the up- and the downregulated gene sets was shared). This is partly explained by limited coverage of RNA-seq genes by the microarrays, which was about 60%. Apart from different methodology, low similarities in individual genes in independent studies analyzing primary tissue are a known problem . However, GSEA can reveal shared pathways between independent studies. Indeed, GSEA of the upregulated genes in cohort 2 overlapped with functional pathways “TNF signaling via NFkB”, “inflammatory responses”, and “IL-2-STAT5 signaling”, which showed strong similarity with cohort 1 functional gene sets. Downregulated genes were enriched in functional groups such as “notch signaling”, “epithelial-mesenchymal transition”, “G2M checkpoint”, and “hypoxia function” according to GSEA which overlapped with cohort 1 in “G2M checkpoint” and “hypoxia” functional groups (Fig. 1D).
Thus, GSEA of both cohorts suggests an imprint of inflammation in axSpA monocytes.
Co-expression analysis of the axSpA signature genes points towards distinct LPS/TNFα activation and reveals a neutrophil-like phenotype of monocytes
To decipher transcriptional alterations of axSpA monocytes in the context of BM adaptation and LPS/cytokine activation, we applied mapping of monocyte profiles with previously published profiles from GEO data repositories. Up- and downregulated probe-sets are analyzed separately; for analysis of cohort 2 the matched GeneChip probe-sets were used.
The reference profiles comprised (1) 34 transcriptomes of 11 different cell types of early and late myelopoiesis from bone marrow (BM; GSE42519) ; (2) 15 transcriptomes from blood cells comprising BDCA1+ (n = 3) and BDCA3+ DCs (n = 3), CD15+ PMN (n = 3) and blood monocyte subsets including classical CD14+CD16− monocytes (Mo-CD16−; n = 3), and non-classical monocytes CD14−CD16+ (Mo-CD16+; n = 3)  (GSE18565); (3) six transcriptomes of blood leukocytes from healthy donors before and after treatment with G-CSF (GSE7400) ; and (4) 12 transcriptomes of blood purified monocytes incubated without stimulus or stimulated with TNFα, LPS, IFNγ, or IFNα (GSE38351) [21, 31].
By applying the differentially expressed probe-sets to these transcriptomes, their regulation during myelopoiesis, in subsets of blood cells and upon in vitro activation of monocytes can be studied. The list of differentially expressed probe-sets is retrieved from the reference transcriptomes; quantile normalized and co-regulated probe-sets are identified by calculating Pearson correlation coefficients. Hierarchical clustering of this gene-to-gene correlation (correlation matrix) provides the order for display of the reference transcriptomes showing relative expression values. This allows visual identification of co-expression clusters among the differentially expressed probe-sets in the reference transcriptomes.
Among the probe-sets upregulated in axSpA patients of cohort 1, we found five clusters with distinct co-expression in BM and blood cells (C1-C5, Fig. 2A). Most strikingly, genes of cluster 1–3 and cluster 5 were overrepresented in late stage granulocytopoiesis and PMNs but weakly expressed in DCs. Moreover, cluster 1 and 5 probe-sets were also enriched after G-CSF treatment. These co-expression cluster included genes such as TLR6, CCR1, NFkB, STAT3, and CSF3R which are related to inflammatory and PMN effector functions such as bacterial recognition, ROS production and phagocytosis .
In addition, a co-expression cluster comprising fewer genes was found in in vitro activated monocytes (cluster 4). GSEA of genes in this cluster indicated overlap with TNF signaling via NFkB, IL-6-JAK-STAT3-signaling and inflammatory responses among others.
Co-expression analysis of the 296 downregulated probe-sets of the axSpA signature of cohort 1 revealed three clusters (Fig. 2B). Most strikingly, genes of cluster 1 and 3, comprising > 50% of the downregulated probe-sets, were highly expressed in DCs but strongly downregulated in PMNs and late granulocytopoiesis. GSEA of these cluster genes revealed overlap with metabolic pathways such as oxidative phosphorylation, cholesterol homeostasis, and adipogenesis, cell cycle control, and mTORC1 signaling. These metabolic pathways are constitutively active in resting DCs while PMNs relay on glycolytic pathways. Genes of cluster 2, comprising less than 10% of the downregulated genes in axSpA, were upregulated in late granulopoiesis and related to TNFα signaling and hypoxia.
Co-expression profiling of significant genes of cohort 2 revealed similar changes as in cohort 1. Thus, co-expression of the upregulated genes was found in late granulocytopoiesis and blood PMNs as well as in LPS/cytokine stimulated monocyte signatures. Among the downregulated genes a co-expression cluster was found in DCs (Fig. 2C, D). Genes of this cluster were in opposite downregulated in PMNs and late granulocytopoiesis.
Altogether, upregulated genes of the axSpA signature are found in late PMNs of the BM, in blood PMNs, and in response to G-CSF while downregulated genes are DC associated. Thus, axSpA monocytes exhibit a neutrophil-like phenotype but deprivation of DC features which is compatible with a shift towards GMP-dependent generation of monocytes. Together with the LPS/TNF imprint, this may reflect response and adaptation to continuous microbial stimulation possibly originating from the gut.
LBP serum levels are increased in axSpA patients, predominantly in those with chronic gut inflammation
To assess putative microbial translocation we measured serum concentration of I-FABP, βDG and LBP in patients of cohort 1, in which none of the patient suffered from clinical IBD. Pathological elevations of these parameters were reported before in clinical conditions of intestinal inflammation such as coeliac disease, HIV enteropathy, and Crohns disease [33,34,35].
As a result, we found no difference in I-FABP serum levels between axSpA and HC (Fig. 3A) and βDG levels were undetectable in both axSpA and HC samples (while pathological elevation was found in patients with metabolic syndrome in our hands; unpublished observations B.S.). However, LBP serum levels were significantly higher in axSpA patients compared to HC (p < 0.001; Fig. 3B). LBP levels tended to be higher in patients with active disease (BASDAI > 4) although this did not reach statistical difference (Fig. 3C) and did not differ between nr-axSpA and AS patients (Fig. 3D).
To determine whether LBP levels are related to intestinal inflammation in axSpA patients, we analyzed LBP levels in axSpA patients of the GIANT cohort in whom presence of intestinal inflammation was assessed by colonoscopy and microscopic evaluation. Interestingly, LBP levels were higher in patients with microscopic gut inflammation than in patients without gut inflammation. The highest LBP levels were found in patients with chronic gut inflammation suggesting a relation between LBP levels and gut inflammation in axSpA patients (Fig. 3E).
“Active-axSpA” monocyte signature reveals association with inflammatory response, gene pattern of terminally differentiated CD16+ blood monocytes, and suppression of DC-associated genes
To determine, if disease activity is associated with distinct transcriptional changes in monocytes, we compared transcriptomes of axSpA patients from cohort 1 separated into active patients (BASDAI ≥ 4 = BASDAIhigh, n = 12) and inactive patients (BASDAI < 4 = BASDAIlow, n = 13). The group-wise comparison of these two groups of patients identified 1367 probe-sets differentially expressed (HPCDA score > 100). Hierarchical clustering showed the heterogeneity between patients but very distinctive separation of nine BASDAIhigh patients from eight BASDAIlow patients, while three BASDAIhigh patients clustered outside these groups and five BASDAIlow patients clustered in-between these groups (Fig. 4A).
Out of 1367 probe-sets of this “active axSpA signature”, 535 were up- and 832 were downregulated in monocytes of BASDAIhigh patients. GSEA of upregulated genes identified associations with inflammatory pathways such as “IL-2-STAT5”, “TNFα”, “TGFβ”, and “Notch signaling”, as well as with “IFNγ response” among other pathways. The downregulated genes showed an overlap with functional groups such as “myc target”, “unfolded protein response”, and “adipogenesis” (Fig. 4B).
Up- and downregulated genes of the “active axSpA signature” were analyzed for co-expression in 70 reference transcriptomes as described above. For the upregulated genes, in total four gene clusters were identified, where cluster 1 and 2 genes showed co-expression in early and late granulocytopoiesis, cluster 3 in LPS/cytokine-activated monocytes and cluster 4 showed distinct co-expression in the non-classical CD16+ monocyte subset in the blood.
GSEA of cluster 3 und 4 genes indicated overlap with functional pathways “interferon gamma response”, “IL2-STAT5 signaling”, “inflammation response”, “TNFα signaling via NFkB”, “TGFβ signaling”, “mTORC1”, and “Notch signaling”.
In contrast, among the 832 probe-sets downregulated in BASDAIhigh patients (and hence upregulated in BASDAIlow patients), high expression was found in DCs (cluster 1a–d; Fig. 4D). Genes of this cluster were underrepresented in PMNs and progressively downregulated in late granulocytopoieses. GSEA revealed overlap with metabolic pathways such as “mTORC1” (cluster 1a), “Myc targets”, “G2M checkpoint”, “E2F targets”, “DNA repair”, and “adipogenesis” (cluster 1b). A small co-expression cluster was found in PMNs (cluster 2).
Thus, the “active SpA signature” points towards LPS/cytokine activation, changes in monocyte subset mobilization (CD16+-associated gene enhancement of non-classical monocytes) and suppression of DC-associated genes in monocytes of active axSpA patients.
Monocyte transcriptomes of nr-axSpA and AS patients differ in pathways controlling metabolism in DCs
Nr-axSpA and AS are considered consecutive stages of axSpA differing in presence of structural damage. However, progression to the radiographic stage is slow and limited to a small group of patients  suggesting that genetic heterogeneity among patients may determine progression to AS rather than time.
Since monocytes give rise to osteoclast which promote bone destruction in AS , we performed pair-wise comparisons of transcriptomes of nr-axSpA and AS patients of cohort 1 to search for an AS-specific monocyte transcriptomic profile. In total, 562 probe-sets were identified to be differentially expressed, again with low magnitudes of changes in their expression. Hierarchical clustering demonstrated that all 10 nr-axSpA clustered separately from the AS group, but 3 AS clustered with the nr-axSpA patients (Fig. 5A). PCA with these 562 significant probe-sets confirmed distinctive differences in monocyte transcriptomes of 10 nr-axSpA and 15 AS patients (Fig. 5B).
GSEA of probe-sets upregulated in AS showed overlap with the functional group “apoptosis”, while probe-sets downregulated in AS patients mapped to metabolic pathways such as oxidative phosphorylation, myc and mTORC1 (Fig. 5C).
Co-expression analysis of the 145 upregulated AS signature probe-sets revealed co-expression of cluster 1 and 2 in early and late myelopoiesis (Fig. 5D). In contrast, probe-sets downregulated in AS (Fig. 5E) where clearly overrepresented in DCs (cluster 1, 2, and 4) and repressed in PMNs and late granulocytopoiesis. Only genes of cluster 3, comprising less than 25% of the downregulated probe-sets, were higher expressed in late granulocytopoiesis and downregulated in DCs (cluster 3). Cluster 1, 2, and 4 genes overlapped with metabolic pathways such as “myc pathway”, “unfolded protein response”, and “oxidative phosphorylation”.
Thus, monocytes of AS patients show a stronger shift towards a neutrophil-like phenotype as compared to nr-axSpA patients which express more DC-like features suggesting a stronger adaptation of monocytopoiesis in the BM of AS patients.
Gut inflammation is a common comorbidity in axSpA; however, the impact on systemic immunity is poorly understood. Monocytes as sensors of microbes might be particularly affected and previous studies reported a proinflammatory, i.e., activated phenotype of monocytes in axSpA patients [6, 7]. To study this monocytic activation in more detail, we examined transcriptional pattern from monocytes of axSpA patients for signs of activation and developmental adaptations.
As a result of pathway and co-expression analysis of differentially expressed genes from two independent cohorts, we found an “axSpA monocyte signature” showing a small, but distinct, LPS/TNF imprint and more prominent signs of BM adaptation. Specifically, we found an enrichment of inflammatory pathways such as “TNF signaling”, “inflammatory response”, “allograft rejection”, and “interferon gamma” in the axSpA signature. However, co-expression analysis using reference transcriptomes of myelopoietic cells, blood cell subsets, and in vitro activated monocytes placed these transcriptional changes into the context of a shifted monocytopoiesis. In fact, only a small part of the axSpA signature genes was found to be regulated by LPS/cytokine treatment in monocytes. In contrast, there was co-expression of upregulated axSpA signature genes in late granulocytopoietic BM cells, BM PMNs, and CD15+ blood PMNs while downregulated axSpA genes were DC-associated genes. Thus axSpA monocytes have a neutrophil-like phenotype with deprivation of DC features compatible with a skewing of monocytopoiesis towards GMP-driven monocytopoiesis . Similarly, an expansion of GMP-driven hematopoiesis was recently described in experimental SpA in SKG mice . The severe inflammation of the small intestine that accompanies SpA induction in this model is considered the stimulus for the changes in hematopoiesis.
In this study here, none of the axSpA patients suffered from overt IBD and of translocation markers we only found increased LBP levels, but no increase of I-FABP and βDG levels. However, in axSpA cohort 3 with histological assessment of gut inflammation, we found highest LBP levels in patients with signs of chronic gut inflammation suggesting a relation of LBP levels to gut pathology in axSpA. As LBP levels are controlled even under homeostatic conditions by gut microbiota , they appear more sensitive to changes in gut permeability than for instance I-FABP. LBP is also an acute phase protein like CRP ; however, in axSpA patients, we found elevated LBP levels without elevation of CRP in most patients. Both, the elevated LBP levels but also the transcriptional axSpA signature of monocytes with signs of BM adaptation and LPS/TNF activation are compatible with low level but chronic bacterial translocation in these patients. Alternatively, a direct effect of HLA-B27 on monocyte development and hematopoiesis could be discussed and needs experimental verification.
A methodological issue in our study is the low overlap of individual genes in transcriptomic analysis of cohort 1 and 2. This may be related to differences in molecular techniques (geneChip vs. RNA-seq) , data analysis (HPCDA scoring vs. edgeR based scoring), and cohort population (age, male/female matching, use of biologics). A low overlap of individual genes in transcriptomes from primary tissue from independent studies, even if standardized by technology, is a known hurdle and pathway analysis using GSEA or co-expression analysis can be used to detect functional similarities between signatures . In this line, GSEA found great overlap in functional groups in particular of upregulated genes between both cohorts and co-expression analysis revealed the abovementioned LPS/TNF imprint and signs of BM adaptation in both cohorts. However, even though GSEA and co-expression analysis link both studies further confirmation of the transcriptional axSpA signature for instance by protein expression or single cell profiling of monocytes is required and would significantly enhance the relevance of our findings. Furthermore, additional studies should compare monocyte transcriptomes between different rheumatic and inflammatory disease to better define the specificity of the signature.
To determine if disease activity in axSpA is associated with transcriptional changes in monocytes we compared transcriptomes from active (BASDAI ≥ 4) vs. inactive patients (BASDAI < 4) in cohort 1, in which none of the patients was treated with biologics. The “active axSpA signature” resulted in only reasonable clustering of active and inactive patients, which is not unexpected given that BASDAI is patient-reported without molecular parameters of inflammation. However, “active axSpA signature” genes were enriched in inflammatory pathways such as IL-2-STAT5 signaling, TNFα signaling, and interferon gamma response while downregulated pathways comprised metabolic pathways. Distinct co-expression of the upregulated “active axSpA” genes was found in LPS/cytokine-activated monocytes, G-CSF mobilized cells and in CD16+ non-classical monocytes. This CD16+ non-classical monocyte pattern may even be underrepresented in our study since we enriched for CD14 which is lower expressed on CD16+ cells. Non-classical CD16+ monocytes usually recirculate in blood while classical CD16− monocytes migrate into peripheral tissues . In view of the co-signature with G-CSF mobilized cells, the enrichment of the CD16+ non-classical monocyte pattern may reflect enhanced monocyte mobilization and enhanced extravasation of classical monocytes, for instance into the gut, in active axSpA.
We also compared transcriptomes of patients classified as nr-axSpA and AS which are considered successive stages of axSp A. Interestingly, we found an “AS monocyte signature” which clearly separated most AS from nr-axSpA patients. Upregulated genes were coexpressed in myelopoiesis and involved apoptosis pathways while downregulated genes were highly expressed in DCs but strongly suppressed in PMNs and late granulocytopoietic cells. These genes are enriched in metabolic pathways such as fatty acid metabolism, mTORC1 signaling and oxidative phosphorylation used by resting DCs but not PMNs for energy consumption [43, 44]. This deprivation of resting DC features in AS patients suggests a stronger BM adaptation in AS patients towards GMP-driven monopoiesis. Interestingly and in line with similar clinical disease activity in AS and nr-axSpA patients, no signs of enhanced acute inflammatory responses (i.e., monocyte mobilization and LPS/TNF response) are found in the AS signature.
Changes in monocytopoiesis may also affect development of structural damage as osteoclasts derived from monocytes. In AS patients, overactivity of RANK-mediated osteoclastogenesis was found related to sacroiliac joint ankylosis . Favored oxidative phosphorylation in nr-axSpA monocytes may result in lower bone destruction as oxidative phosphorylation is linked to reduced osteoclastic degradation and glycolysis to osteoclastic destruction .
In conclusion, we here describe changes in monocyte transcriptomes of axSpA patients which are compatible with LPS/TNF activation of monocytes and a more prominent adaptation of monocytopoiesis resulting in a shift from MDP towards GMP-driven monopoiesis. Together with the LBP increase in axSpA patients, this may result from low-grade translocation from the gut. These findings provide a mechanistic link between gut inflammation and bone marrow, the site of skeletal inflammation and destruction.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Lipopolysaccharide binding protein
Non-radiographic axial SpA
Inflammatory bowel disease
Monocyte-dendritic cell progenitors
Intestinal fatty acid binding protein
lipopolysaccharide binding protein
Non-steroidal anti-inflammatory drugs
Gent Inflammatory Arthritis and spoNdylitis cohort
High performance chip data analysis
Principal component analysis
Gene set enrichment analysis
Molecular Signature Database
Gene expression omnibus
Brewerton DA, Hart FD, Nicholls A, Caffrey M, James DC, Sturrock RD. Ankylosing spondylitis and HL-A 27. Lancet. 1973;1(7809):904–7. https://doi.org/10.1016/S0140-6736(73)91360-3.
Cortes A, Hadler J, Pointon JP, Robinson PC, Karaderi T, Leo P, et al. Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci. Nat Genet. 2013;45(7):730–8. https://doi.org/10.1038/ng.2667.
Ellinghaus D, Jostins L, Spain SL, Cortes A, Bethune J, Han B, et al. Analysis of five chronic inflammatory diseases identifies 27 new associations and highlights disease-specific patterns at shared loci. Nat Genet. 2016;48(5):510–8. https://doi.org/10.1038/ng.3528.
Gracey E, Vereecke L, McGovern D, Frohling M, Schett G, Danese S, et al. Revisiting the gut-joint axis: links between gut inflammation and spondyloarthritis. Nat Rev Rheumatol. 2020;16(8):415–33. https://doi.org/10.1038/s41584-020-0454-9.
Menegatti S, Guillemot V, Latis E, Yahia-Cherbal H, Mittermuller D, Rouilly V, et al. Immune response profiling of patients with spondyloarthritis reveals signalling networks mediating TNF-blocker function in vivo. Ann Rheum Dis. 2020;80(4):475–86. https://doi.org/10.1136/annrheumdis-2020-218304.
Wright C. Edelmann M, diGleria K, Kollnberger S, Kramer H, McGowan S, McHugh K, Taylor S, Kessler B, Bowness P: Ankylosing spondylitis monocytes show upregulation of proteins involved in inflammation and the ubiquitin proteasome pathway. Ann Rheum Dis. 2009;68(10):1626–32. https://doi.org/10.1136/ard.2008.097204.
Conrad K, Wu P, Sieper J, Syrbe U. In vivo pre-activation of monocytes in patients with axial spondyloarthritis. Arthritis Res Ther. 2015;17(1):179. https://doi.org/10.1186/s13075-015-0694-2.
Ziegler-Heitbrock L, Ancuta P, Crowe S, Dalod M, Grau V, Hart DN, et al. Nomenclature of monocytes and dendritic cells in blood. Blood. 2010;116(16):e74–80. https://doi.org/10.1182/blood-2010-02-258558.
Patel AA, Zhang Y, Fullerton JN, Boelen L, Rongvaux A, Maini AA, et al. The fate and lifespan of human monocyte subsets in steady state and systemic inflammation. J Exp Med. 2017;214(7):1913–23. https://doi.org/10.1084/jem.20170355.
Tak T, van Groenendael R, Pickkers P, Koenderman L. Monocyte subsets are differentially lost from the circulation during acute inflammation induced by human experimental endotoxemia. J Innate Immun. 2017;9(5):464–74. https://doi.org/10.1159/000475665.
Stolwijk C, van Tubergen A, Castillo-Ortiz JD, Boonen A. Prevalence of extra-articular manifestations in patients with ankylosing spondylitis: a systematic review and meta-analysis. Ann Rheum Dis. 2013;72(Suppl 3):A520.2–A521. https://doi.org/10.1136/annrheumdis-2013-eular.1558.
Mielants H, Veys EM, Cuvelier C, de Vos M. Ileocolonoscopic findings in seronegative spondylarthropathies. Br J Rheumatol. 1988;27(Suppl 2):95–105.
Van Praet L, Van den Bosch FE, Jacques P, Carron P, Jans L, Colman R, et al. Microscopic gut inflammation in axial spondyloarthritis: a multiparametric predictive model. Ann Rheum Dis. 2013;72(3):414–7. https://doi.org/10.1136/annrheumdis-2012-202135.
Chavakis T, Mitroulis I, Hajishengallis G. Hematopoietic progenitor cells as integrative hubs for adaptation to and fine-tuning of inflammation. Nat Immunol. 2019;20(7):802–11. https://doi.org/10.1038/s41590-019-0402-5.
Yanez A, Coetzee SG, Olsson A, Muench DE, Berman BP, Hazelett DJ, et al. Granulocyte-monocyte progenitors and monocyte-dendritic cell progenitors ndependently produce eunctionally distinct monocytes. Immunity. 2017;47(5):890–902 e894. https://doi.org/10.1016/j.immuni.2017.10.021.
Rudwaleit M, van der Heijde D, Landewe R, Listing J, Akkoc N, Brandt J, et al. The development of Assessment of SpondyloArthritis International Society classification criteria for axial spondyloarthritis (part II): validation and final selection. Ann Rheum Dis. 2009;68(6):777–83. https://doi.org/10.1136/ard.2009.108233.
van der Linden S, Valkenburg HA, Cats A. Evaluation of diagnostic criteria for ankylosing spondylitis. A proposal for modification of the New York criteria. Arthritis Rheum. 1984;27(4):361–8. https://doi.org/10.1002/art.1780270401.
Garrett S, Jenkinson T, Kennedy LG, Whitelock H, Gaisford P, Calin A. A new approach to defining disease status in ankylosing spondylitis: the Bath Ankylosing Spondylitis Disease Activity Index. J Rheumatol. 1994;21(12):2286–91.
Lukas C, Landewe R, Sieper J, Dougados M, Davis J, Braun J, et al. Assessment of SpondyloArthritis International S: Development of an ASAS-endorsed disease activity score (ASDAS) in patients with ankylosing spondylitis. Ann Rheum Dis. 2009;68(1):18–24. https://doi.org/10.1136/ard.2008.094870.
Cypers H, Varkas G, Beeckman S, Debusschere K, Vogl T, Roth J, et al. Elevated calprotectin levels reveal bowel inflammation in spondyloarthritis. Ann Rheum Dis. 2016;75(7):1357–62. https://doi.org/10.1136/annrheumdis-2015-208025.
Smiljanovic B, Grun JR, Steinbrich-Zollner M, Stuhlmuller B, Haupl T, Burmester GR, et al. Defining TNF-alpha- and LPS-induced gene signatures in monocytes to unravel the complexity of peripheral blood transcriptomes in health and disease. J Mol Med (Berl). 2010;88(10):1065–79. https://doi.org/10.1007/s00109-010-0648-8.
Lepre J, Rice JJ, Tu Y, Stolovitzky G. Genes@Work: an efficient algorithm for pattern discovery and multivariate feature selection in gene expression data. Bioinformatics. 2004;20(7):1033–44. https://doi.org/10.1093/bioinformatics/bth035.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50. https://doi.org/10.1073/pnas.0506580102.
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34(3):267–73. https://doi.org/10.1038/ng1180.
Smiljanovic B, Grutzkau A, Sorensen T, Grun JR, Vogl T, Bonin M, et al. Synovial tissue transcriptomes of long-standing rheumatoid arthritis are dominated by activated macrophages that reflect microbial stimulation. Sci Rep. 2020;10(1):7907. https://doi.org/10.1038/s41598-020-64431-4.
Grützkau A, Radbruch A. Separation of whole blood cells and its impact on gene expression. In: Bosio A, Gerstmayer B, editors. Birkhäuser Basel Microarrays in inflammation progress in inflammation research; 2008.
Wesemann DR, Benveniste EN. STAT-1 alpha and IFN-gamma as modulators of TNF-alpha signaling in macrophages: regulation and functional implications of the TNF receptor 1:STAT-1 alpha complex. J Immunol. 2003;171(10):5313–9. https://doi.org/10.4049/jimmunol.171.10.5313.
Ma X, Jiang Z, Li N, Jiang W, Gao P, Yang M, et al. Ets2 suppresses inflammatory cytokines through MAPK/NF-kappaB signaling and directly binds to the IL-6 promoter in macrophages. Aging (Albany NY). 2019;11(22):10610–25. https://doi.org/10.18632/aging.102480.
Rapin N, Bagger FO, Jendholm J, Mora-Jensen H, Krogh A, Kohlmann A, et al. Comparing cancer vs normal gene expression profiles identifies new disease entities and common transcriptional programs in AML patients. Blood. 2014;123(6):894–904. https://doi.org/10.1182/blood-2013-02-485771.
Buzzeo MP, Yang J, Casella G, Reddy V. Hematopoietic stem cell mobilization with G-CSF induces innate inflammation yet suppresses adaptive immune gene expression as revealed by microarray analysis. Exp Hematol. 2007;35(9):1456–65. https://doi.org/10.1016/j.exphem.2007.06.001.
Smiljanovic B, Grun JR, Biesen R, Schulte-Wrede U, Baumgrass R, Stuhlmuller B, et al. The multifaceted balance of TNF-alpha and type I/II interferon responses in SLE and RA: how monocytes manage the impact of cytokines. J Mol Med (Berl). 2012;90(11):1295–309. https://doi.org/10.1007/s00109-012-0907-y.
Evrard M, Kwok IWH, Chong SZ, Teng KWW, Becht E, Chen J, et al. Developmental analysis of bone marrow neutrophils reveals populations specialized in expansion, trafficking, and effector functions. Immunity. 2018;48(2):364–79 e368. https://doi.org/10.1016/j.immuni.2018.02.002.
Adriaanse MP, Leffler DA, Kelly CP, Schuppan D, Najarian RM, Goldsmith JD, et al. Serum I-FABP detects gluten responsiveness in adult celiac disease patients on a short-term gluten challenge. Am J Gastroenterol. 2016;111(7):1014–22. https://doi.org/10.1038/ajg.2016.162.
Lakatos PL, Kiss LS, Palatka K, Altorjay I, Antal-Szalmas P, Palyu E, et al. Serum lipopolysaccharide-binding protein and soluble CD14 are markers of disease activity in patients with Crohn's disease. Inflamm Bowel Dis. 2011;17(3):767–77. https://doi.org/10.1002/ibd.21402.
Ramendra R, Isnard S, Mehraj V, Chen J, Zhang Y, Finkelman M, et al. Circulating LPS and (1-->3)-beta-D-Glucan: a folie a deux contributing to HIV-associated immune activation. Front Immunol. 2019;10:465. https://doi.org/10.3389/fimmu.2019.00465.
Ramiro S, van der Heijde D, van Tubergen A, Stolwijk C, Dougados M, van den Bosch F, et al. Higher disease activity leads to more structural damage in the spine in ankylosing spondylitis: 12-year longitudinal data from the OASIS cohort. Ann Rheum Dis. 2014;73(8):1455–61. https://doi.org/10.1136/annrheumdis-2014-205178.
Bleil J, Maier R, Hempfing A, Schlichting U, Appel H, Sieper J, et al. Histomorphological and histomorphometric characteristics of zygapophyseal joint remodelling in ankylosing spondylitis. Arthritis Rheum. 2014;66(7):1745–54. https://doi.org/10.1002/art.38404.
Regan-Komito D, Swann JW, Demetriou P, Cohen ES, Horwood NJ, Sansom SN, et al. GM-CSF drives dysregulated hematopoietic stem cell activity and pathogenic extramedullary myelopoiesis in experimental spondyloarthritis. Nat Commun. 2020;11(1):155. https://doi.org/10.1038/s41467-019-13853-4.
Molinaro A, Koh A, Wu H, Schoeler M, Faggi MI, Carreras A, et al. Hepatic expression of lipopolysaccharide-binding protein (Lbp) is induced by the gut microbiota through Myd88 and impairs glucose tolerance in mice independent of obesity. Mol Metab. 2020;37:100997. https://doi.org/10.1016/j.molmet.2020.100997.
Blairon L, Wittebole X, Laterre PF. Lipopolysaccharide-binding protein serum levels in patients with severe sepsis due to gram-positive and fungal infections. J Infect Dis. 2003;187(2):287–91. https://doi.org/10.1086/346046.
Zhao S, Fung-Leung WP, Bittner A, Ngo K, Liu X. Comparison of RNA-Seq and microarray in transcriptome profiling of activated T cells. PLoS One. 2014;9(1):e78644. https://doi.org/10.1371/journal.pone.0078644.
Deodhar A, Strand V, Kay J, Braun J. The term 'non-radiographic axial spondyloarthritis' is much more important to classify than to diagnose patients with axial spondyloarthritis. Ann Rheum Dis. 2016;75(5):791–4. https://doi.org/10.1136/annrheumdis-2015-208852.
Pearce EJ, Everts B. Dendritic cell metabolism. Nat Rev Immunol. 2015;15(1):18–29. https://doi.org/10.1038/nri3771.
Injarabian L, Devin A, Ransac S, Marteyn BS. Neutrophil metabolic shift during their lifecycle: impact on their survival and activation. Int J Mol Sci. 2019;21(1). https://doi.org/10.3390/ijms21010287.
Im CH, Kang EH, Ki JY, Shin DW, Choi HJ, Chang EJ, et al. Receptor activator of nuclear factor kappa B ligand-mediated osteoclastogenesis is elevated in ankylosing spondylitis. Clin Exp Rheumatol. 2009;27(4):620–5.
Lemma S, Sboarina M, Porporato PE, Zini N, Sonveaux P, Di Pompo G, et al. Energy metabolism in osteoclast formation and activity. Int J Biochem Cell Biol. 2016;79:168–80. https://doi.org/10.1016/j.biocel.2016.08.034.
Funding: This work was supported by the FOREUM — Foundation for Research in Rheumatology. B.S. was supported by the German Federal Ministry of Education and Research (BMBF) with the network project ArthroMark (01EC1009A). Open Access funding enabled and organized by Projekt DEAL.
Ethics approval and consent to participate
All participants gave written consent to the study, which was approved by the local ethical committees of the university hospitals Charité (Berlin, Germany), Ile-de-France XI (Saint-Germaine-en-Laye, France), and Ghent University Hospital (Ghent, Belgium).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Karow, F., Smiljanovic, B., Grün, J.R. et al. Monocyte transcriptomes from patients with axial spondyloarthritis reveal dysregulated monocytopoiesis and a distinct inflammatory imprint. Arthritis Res Ther 23, 246 (2021). https://doi.org/10.1186/s13075-021-02623-7
- Gut translocation