1. Introduction
B lymphocytes represent about 10-20% of total circulating lymphocytes. They can be divided in four major subpopulations, based on the presence of IgD and CD27 surface molecules [
1]. IgD, which is expressed by naïve B cells, progressively diminishes and eliminates after isotype-switching. Instead, CD27, a receptor molecule belonging to tumor necrosis factor (TNF) superfamily, is not expressed by naïve B lymphocytes, yet it appears at later differentiation stages, and its presence signifies memory B cells [
2,
3]. Therefore, IgD+CD27- cells are described as naïve, IgD+CD27+ cells as non-switched memory and IgD-CD27+ as switched memory B cells. Although the above 3 subtypes are clearly defined in terms of origin, differentiation and function, the fourth one, Double Negative (DN) B cells, IgD-CD27- is still an enigmatic subpopulation. The absence of IgD indicates isotype switching, while absence of CD27 specifies less differentiated cells [
4,
5]. For many years these cells were characterized as a naïve B cell population and they were under-recognized, yet only recently, their origin has caused major interest. Hypothetical models of their derivation have been described, based on their phenotypic and functional characteristics. Therefore, based on the increased expression of B cell lymphoma protein, similar to that of CD27+ cells, DN B cells are anticipated to originate after a defective reactions in germinal centers or derive from memory B cells following displacement or down-regulation of the CD27 molecule, during immune-senescent or immune-exhaustion process [
6,
7,
8].
DN B cells, described as DN1 and DN2, according to the presence or absence of CD24, CD38 and CXCR5 proteins, respectively, are increased during ageing and also, in chronic inflammatory or systemic diseases, as a result of chronic stimulation [
4] Systemic Lupus Erythematosus is characterized by increased expression of DN2 B cells, probably as a result of a hyper-responsiveness to Toll-like receptor-7 signaling [
9]. Recent studies have shown that DN2 B cells were mainly increased in active SLE and in the presence of lupus nephritis and they were positively correlated with anti-RNA, anti-Sm and anti-RNP autoantibodies [
9,
10,
11]. Although many efforts have been made to assess the correlation of DN B cells with markers of SLE activity, only poor literature exists in the estimation of their origin [
5,
6]. B cells are activated by pathogens, antigen presenting cells and, also, by T lymphocytes. We have previously described the predominance of DN B cells in SLE patients, followed by relatively increased expression of switched memory B cells (CD19+IgD-CD27+), even within patients being on remission. Moreover, DN B cells in lupus nephritis patients were significantly increased compared to end stage renal disease patients on hemodialysis [
11]. Although patients on hemodialysis were selected as a representative group of patients with senescent adaptive immunity, our results proved a clear predominance of more differentiated B cell subtypes in lupus nephritis patients compared to hemodialysis [
8,
11]. The same was not confirmed with T lymphocyte subpopulations, in which, senescent phenotypic changes were prominent in dialysis patients [
12,
13].
Based on these previous findings, namely, significant increase of DN B lymphocytes in SLE patients, and important alterations in the T lymphocytes and their subsets, we decided to evaluated possible correlations between the two lymphocyte compartments in patients with lupus nephritis. T lymphocytes were divided in subpopulations, according to their differential status, and the potential correlation with DN B cells was estimated.
2. Materials and Methods
In this is a cross sectional study, we included patients with SLE in remission. After analyzing the phenotype of peripheral B and T lymphocytes, we estimated frequencies of DN B lymphocytes, and assessed potential correlation between them and T lymphocyte subsets.
2.1. Patients
The study population involved of 30 Caucasian patients with SLE, and 31 healthy controls (HC), with similar age, sex and ethnicity. Diagnosis of SLE was based on SLICC/ACR 2012 criteria[
14]. Diagnosis of Lupus Nephritis in 20/30 patients was based on kidney biopsy findings.
Inclusion criteria: All participants were adults, aged 18-67years. Patients should be under regular follow-up, with available clinical, biochemical and immunological data. SLE should be stable with no flare up for at least 2 years, and their immunosuppression, at time of evaluation, could include prednisolone, hydroxychloroquine, calcineurin inhibitors and azathioprine or mycophenolate mofetil.
Exclusion criteria: Comorbidities, such as diabetes mellitus, malignancy, hematological disorder, impaired renal function, defined as eGFR<60ml/min/m2, or presence of chronic active or recent (< 6months) infection. Additionally, patients who had received monoclonal antibodies in the past, or cyclophosphamide during the last 24 months, were also excluded.
Information regarding previous historical data, medication, comorbid conditions, and disease flares were retrieved from patients’ records. At time of assessment, patients’ biochemical and immune profile were recorded and the SLEDAI score was used to estimate activity of SLE [
15].
The study was approved by the Institutional Review Board of the Medical School of Aristotle University of Thessaloniki, Greece. All participants were informed prior to enrollment, and signed the written consent.
2.2. Flow cytometry
Heparinized venous blood from SLE patients and healthy control individuals was drown under sterile conditions and collected in EDTA tubes. Within 12hours following collection, samples proceeded for the evaluation, of total white cell count, total lymphocytes, B and T lymphocytes, and further analysis for CD4 and CD8 subtypes.
The expression of IgD and CD27 on B lymphocytes and of CD45RA, CCR7, CD28, CD31, CD57, on both CD4 and CD8 T lymphocytes was assessed using a cytometer (Navios Flow Cytometer, Beckman Coulter) as described before [
8]. The expression of the above certain surface molecules examined, determined divergent subpopulations. Proportions and counts of all lymphocyte subsets were estimated.
2.3. Flow cytometry reagents
The monoclonal antibodies used in flow cytometry to identify lymphocyte surface receptors, were CD45 PC7, CD19 PC5, IgD FITC and CD27 ECD for B lymphocytes and CD45 PC7, CD3 FITC, CD4 Pacific Blue, CD8 PC5, CD45RA APC, CCR7 PE, CD28 ECD, CD57 FITC, CD31 ECD for CD4 and CD8 lymphocytes.
We used three different tubes (for B, CD4 and CD8 lymphocytes) and combination of estimated receptors resulted in the definition of the following subpopulations.
2.4. Definition of lymphocyte subpopulations
B lymphocyte subpopulation definition was based on the expression of CD27 and IgD surface molecules, and described as naïve (CD19+IgD+CD27-), non-switched and switched memory (CD19+IgD+CD27+ and CD19+IgD+CD27-, respectively), and double negative (DN) B cells (CD19+IgD-CD27-).
Classification of CD4 and CD lymphocytes was based on the presence of the above described surface receptors and described as early differentiated, memory and senescent/advanced differentiated cells..
A. Early differentiated CD4 or CD8 lymphocytes which expressed
1. CD31, CD4+CD31+ and CD8+CD31+, defined as Recent Thymic Emigrants (RTE),
2. CD45RA together with CCR7, CD4+CD45RA+CCR7+ and CD8+CD45RA+CCR7+, defined as Naïve CD4 or CD8 lymphocytes,
3. CD28 and not CD57, CD4CD28+CD57- and CD8CD28+CD57-
4. CD45RA and not CD57, CD4CD45RA+CD57- and CD8CD45RA+CD57-
5. neither CD45RA nor CD57, CD4CD45RA-CD57- and CD8CD45RA-CD57-
B. In the group of memory cells, CD4 or CD8 lymphocytes which expressed
1. CCR7 and not CD45RA, CD4+CD45RA-CCR7+ and CD8+CD45RA-CCR7+, defined as Central memory (CM),
2. neither CCR7 nor CD45RA, CD4+CD45RA-CCR7- and CD8+CD45RA-CCR7-, defined as Effector memory (EM),
C. In the group of Senescent/Advanced differentiated cells,
1. T cells expressing CD45RA and not CCR7, CD4+CD45RA+CCR7- and CD8+CD45RA+CCR7-, defined as EMRA,
2. T cells lacking CD28 molecule, CD4CD28- and CD8CD28-, and their subtypes according to the presence of CD57, CD4+CD28-CD57-, CD4+CD28-CD57+ and CD8+CD28-CD57-, CD8+CD28-CD57+
3. T cells expressing both CD45RA and CD57, CD4CD45RA+CD57+ and CD8CD45RA+CD57+
Gating strategies for B and T lymphocyte subpopulations are described on the supplement file,
Figure 1 and
Figure 2.
2.5. Statistical analysis
Statistical analysis was performed using the Statistical Package for Social Sciences (SPSS Inc., Chicago, IL, USA) version 27.0, for windows. All continues variables were initially tested for normality of distribution, by using the Shapiro-Wilk and/or the Kolmogorov-Smirnov test.
Normally distributed variables were expressed as Mean±Standard Deviation, and Student’s t test for independent samples was performed to compare differences between two groups. Data from Non–parametric variables were expressed as Median and Range, or interquartile range, and Mann Whitney U test was used to assess differences between two groups. Pearson’s and Spearman’s correlation tests were applied to compare variables with normal or non-normal distribution, respectively. Multiple regression analysis was performed to estimate independent factors correlated with DN B lymphocytes, as the dependent variable.
Values of p<0.05 (two-tailed) were considered as statistically significant.
3. Results
3.1. Characteristics of SLE patients
Thirty patients with SLE (29 female) and 31 healthy individuals (17 female), of similar age, sex and ethnicity, were included in the study.
According to inclusion criteria, all SLE patients were in complete or partial remission, under stable treatment for the last 24 months, which included hydroxychloroquine in 30/30 (100%), steroids in 24/30 (80%), MMF in 18/30 (60%) and CNIs in 9/30 (30%) or combination of MMF+CNIs in 5/30 (16.7%). Time since patients were on the above regimes were 84(45-125)months for hydroxychloroquin, 78(45-92)months for steroids, 53(37-65)months for MMF and 42(28-56)months for CNIs.
Twenty patients (67%) had been given cyclophosphamide in the past, more than 2 years ago, while no-one had received monoclonal antibodies in the past.
Clinical and laboratory findings of patients are depicted on
Table 1.
Median time since the initial diagnosis of SLE was 84(45-125)months. SLEDAI score was estimated to 10(2-18) at time of diagnosis and 2(1-5) at time of evaluation. Twenty patients had LN, diagnosed with renal biopsy, performed not more than 5 years prior to evaluation. Mean SCr at time of assessment was 0.89±0.17mg/dl, eGFR 81±19ml/min/1.73m2, Uprot 680±5.600mg/24hr, C3 and C4 levels at time of evaluation were 74.65(29.9-127)mg/dl and 15.7(3.86-27.2)mg/dl, respectively.
3.2. Phenotypic analysis of B lymphocytes in patients with SLE and in HC
Differences in B lymphocytes and their subpopulations are shown on
Table 2. SLE patients had significantly reduced percentage of B cells, and also, reduced proportions of CD19+IgD+CD27+ and CD19+IgD-CD27+ cells, with a clear shift to CD19+IgD-CD27- (DN) B cells. Population of B lymphocytes was significantly reduced compared to controls, almost 1/3 of the HC group, and this reduction affected CD19+IgD+CD27- and CD19+IgD+CD27+ cells (
Table 2).
As the DN B cell population was the only one with an increased percentage within B lymphocytes, we calculated the ratio of (IgD-CD27-)/[(IgD+CD27-)+(IgD-CD27+)+ (IgD+CD27+)] in SLE patients and controls, and found a significant difference 0.14(0.02-2.9) vs. 0.08(0.02-0.54), respectively, p=0.04
The distribution of B lymphocyte subtypes and the differences between SLE patients and controls are depicted in
Figure 1.
Figure 1.
Distribution of B lymphocyte subtypes, IgD+CD27-, IgD+CD27+, IgD-CD27+, IgD-CD27-, in SLE patients (A) and HC (B).
Figure 1.
Distribution of B lymphocyte subtypes, IgD+CD27-, IgD+CD27+, IgD-CD27+, IgD-CD27-, in SLE patients (A) and HC (B).
3.3. Phenotypic analysis of T lymphocytes in SLE patients and HC
CD4 and CD8 lymphocyte subpopulations were classified as early differentiated, memory and advanced differentiated/senescent cells. Classification was based on their phenotypic characteristics, namely the presence of specific surface receptors.
In SLE patients, there was a significant reduction of the whole cohort CD4 but not of the CD8 lymphocytes. More specifically, almost all early differentiated and memory CD4 subtypes were significantly reduced in SLE patients, as were the CD4 EMRA cells. The expression of CD28 molecule on CD4 lymphocytes was reduced but not statistically significant (
Table 3).
Regarding CD8 lymphocytes, CD8EMRA (CD8CD45RA+CCR7-) and CD8EMRACD28- (CD8CD45RA+CCR7-CD28-) cells showed a statistically significant reduction in SLE patients, compared to controls, which proved a lessening in the senescent phenotype, together with a shift to CM and early differentiated CD8 subtypes (
Table 4).
3.4. Distribution of CD4 and CD8 subtypes according to their differentiation status
In
Figure 2 and
Figure 3 the differences in CD4 and CD8 subpopulations are depicted. According to the presence or absence of CD45RA, CD28 and CD57 molecule, T lymphocytes where divided in subgroups, and their percentages were calculated in SLE patients and HC.
Figure 2.
Distribution of CD4 subtypes: the combination of CD45RA expression and CD57 in SLE patients (A) and controls (B), CD28 and CD57 in SLE patients (C) and HC (D).
Figure 2.
Distribution of CD4 subtypes: the combination of CD45RA expression and CD57 in SLE patients (A) and controls (B), CD28 and CD57 in SLE patients (C) and HC (D).
Figure 3.
Distribution of CD8 subtypes: the combination of CD45RA expression and CD57 in SLE patients (A) and controls (B), CD28 and CD57 in SLE patients (C) and HC (D).
Figure 3.
Distribution of CD8 subtypes: the combination of CD45RA expression and CD57 in SLE patients (A) and controls (B), CD28 and CD57 in SLE patients (C) and HC (D).
3.5. Correlation of DN B cells with B lymphocyte subpopulations
Table 5 describes correlation of DN B cells with the rest B cell subpopulations, in SLE and HC. Multiple regression analysis, with DN cells as dependent and the rest of B cell populations as independent factors showed that in both patients and controls CD19+IgD-CD27+ was the only independent parameter correlated to DN B cells, R
2=0.68, p=0.001 for patients and R
2=0.7, p<0.0001, for HC.
3.6. Correlation of DN B cells with T lymphocyte subpopulations
3.6.1. Correlation with CD4 cells
Both in patients and controls, DN B cells showed a significant positive correlation with the number of CD4 cells. In HC DN B cells showed a positive correlation with few of the early differentiated CD4 subsets, namely CD4CD45RA+CD57- and CD4CD28+CD57-, while in SLE patients, almost all early differentiated CD4 subtypes had a significant positive correlation with DN B cells, RTEs, and all subtypes lucking CD57 and expressing CD28 on their membrane. CM CD4 cells also correlated to DN B cells.
3.6.2. Correlation with CD8 cells
Regarding CD8 lymphocytes, DN B cells in HC had no significant correlation neither with CD8 cells, nor with their subpopulations. On the other hand, in SLE patients DN B cells showed major association with RTEs, Naïve CD8 cells, CD8CD45RA-CD57-, CD8CD28+CD57- and CD8CD28+CD57+. No correlation was evident with memory or senescent CD8 cells.
Correlations of DN B cells with CD4 and CD8 lymphocyte counts and their susets are depicted on
Table 6, demonstrating the significant association with early differentiated T lymphocytes.
Multiple Regression analysis, including DN B cells as dependent and CD4, CD8 as independent variables showed that, for SLE patients, CD8 cells was the independent parameter associated with DN B cells, R2=0.277, p=0.01, b coefficient=0.032, CI=0.01-0.05, p=0.006
Therefore, CD8 subpopulations were included in a new multiple regression analysis, which revealed a model including CD8CD45RA-CD57- and CD8CD28+CD57- as the only independent factors related to DN B cells, with adjusted R2=0.464, p<0.0001.
The above model remained stable after adding the rest of CD8 and CD4 subtypes, but was benefited by adding CD4CD31+ cells. Therefore, the model including CD8CD45RA-CD57-, CD8CD28+CD57- and CD4CD31+ cells had an adjusted R2=0.534, p<0.0001.
4. Discussion
The present study was performed in an attempt to estimate possible correlation between T cell immunity with DN B lymphocytes in patients with SLE. DN B cells consist certain subpopulation, which is significantly increased in elderly, but also in several chronic inflammatory diseases, including SLE [
7,
11,
16,
17,
18]. Their origin, function and differential status is still unclear, as they carry characteristics of both naïve, switched memory and immune-exhausted cells [
6,
7]. They may originate from early termination of the GC reaction, from switched memory B cells or from naïve resting B cells [
6].
We have previously described the significant elevation, up to three-fold, of DN B cells in SLE patients compared to healthy controls [
11]. Here we describe a close association of this subtype with certain, less differentiated CD4 and CD8 lymphocyte subpopulations, namely, CD4CD31+, CD4CD45RA+CD28+, CD4CD45RA+CD57-, CD4CD45RA-CD57-, CD4CD28+CD57-, CD4CD28+CD57+, CD8CD31+, Naïve CD8, CD8CD45RA+CD28+, CD8CD45RA-CD57-, CD8CD28+CD57-, and CD8CD28+CD57+ and also, CM CD4 cells. Very interestingly, almost all naïve or at early differentiated T lymphocyte subsets studied showed positive correlation with DN B cells, while no such association was evident for advanced differentiated T cells. Similar correlations were not evident in the healthy control population.
SLE, as well as other systemic autoimmune diseases, is associated with significant increase of DN B lymphocytes, both DN1 and DN2. From those DN2 B cells, are most important, as they act as antigen presenting cells and stimulate T cell response, differentiate into antibody secreting cells, secrete cytokines and modulate immune response, sustained autoimmunity and inflammation [
16,
17,
18,
19]. DN2 B lymphocytes share similar characteristics with naïve B lymphocytes, such as the absence of CXCR5, CD24 and CD38 surface molecules, which are common in the rest of peripheral B lymphocyte subsets, but also retain a unique phenotype, expressing CD69, HLA-DR, CD86 and the regulatory receptors CD32b and CD22, and lacking the lymph node homing receptor L-selectin (CD62L) [
7,
9,
20]
Both DN2 and naïve B lymphocytes express a T-bet transcriptional network, with upregulation of Toll-like receptor-7 (TLR7) and simultaneous reduction of its negative regulators, such as TRAF5 [
9]. As a result, DN2 and naïve B cells are hyper-responsive to TLR7. Stimulation with TLR7 induces phosphorylation of ERK and MAPKp38 pathways, resulting in their continuing activation, production of autoantibodies and generation of peripheral plasmablasts [
21]. Most studies described the above changes in active SLE and active lupus nephritis [
22]. Our patients displayed increased DN B lymphocytes while being in remission, indicating that despite clinical repression, B cell compartment still demonstrates abnormal distribution. However, in the present study, we did not estimate expression of CXRC5, the receptor which distinguishes DN1 (CXCR5+) from DN2 (CXCR5-) B lymphocytes. Previous studies have shown that DN2 cells are increased during disease flares and their levels correlate with autoantibody production [
9,
22].
Although there are several studies confirming predominance and expansion of DN B lymphocytes in SLE patients, and significant correlation with disease pathogenesis, there is only scarcely evidence regarding factors that may stimulate their production [
23,
24]. Correlation between DN B lymphocytes and recent thymic emigrant cells, almost all naïve CD4 and CD8 and central memory CD4 lymphocytes, described in our study, is a completely novel finding which indicates, but cannot prove, causality. Previous studies have shown increased expression of DN B lymphocytes in the cerebrovascular fluid of patients with Multiple Sclerosis, and revealed a T cell stimulatory and proinflamatory function [
25]. Similarly, in Rheumatoid arthritis, DN B lymphocytes seemed to actively influence disease activity and were dowregulated by anti-TNFa treatment [
26]. However, in other studies, the significantly increased levels of DN B lymphocytes were characterized by low expression of HLA-DR, CD40 and CD80 surface molecules, indicating their incapacity to act as antigen presenting cells or interact with T lymphocytes [
26]. In elderly patients DN B cells seemed to behave more as senescent or exhausted cells rather than naïve B cells [
27,
28]. Senescent phenomena were not found to predominate in SLE [
11,
12] although they are prominent findings in chronic inflammatory diseases [
29,
30].
Activation of T lymphocytes and upregulation of DN B lymphocytes seem to happen simultaneously, and it remains unclear whether a common pathogenic mechanism orchestrates the whole reaction, or DN B lymphocytes are themselves implicated in the stimulation of naïve, early differentiated and memory peripheral T lymphocytes. Our findings of close association between DN B cells and early differentiated T lymphocytes in patients with SLE, raise the possibility that DN B lymphocytes are characterized by unique and important characteristics and play a major role in the alterations of adaptive immunity.
Author Contributions
Conceptualization, E.M. and M.S.; methodology, E.M.; software, D-I.B. and S.B.; validation, A.X., G.L. and M.C.; formal analysis, D-I.B.; investigation, M.C. and S.S.; resources, A.I. and M.I.; data curation, D.I.; writing—original draft preparation, E.M.; writing—review and editing, A.F. and M.S.; supervision, A.F. and M.S.; funding acquisition, M.S. All authors have read and agreed to the published version of the manuscript.