- Research article
- Open Access
Serum levels of B-cell activating factor of the TNF family (BAFF) correlate with anti-Jo-1 autoantibodies levels and disease activity in patients with anti-Jo-1positive polymyositis and dermatomyositis
Arthritis Research & Therapyvolume 20, Article number: 158 (2018)
B-cell activating factor of the tumour necrosis factor family (BAFF) plays a role in autoantibody production and is elevated in dermatomyositis (DM) and anti-Jo-1-positive polymyositis (PM). We investigated the inter-relationships between serum levels of BAFF, anti-Jo-1 autoantibodies, and disease activity.
Serum levels of BAFF and anti-Jo-1 antibodies measured by enzyme-linked immunosorbent assay (ELISA) were compared to levels of myoglobin, creatine kinase (CK), aminotransferases (alanine (ALT) and aspartate (AST)), C-reactive protein (CRP), and disease activity assessed by the Myositis Disease Activity Assessment Tool in 63 anti-Jo-1 antibody-positive DM/PM patients. Serial serum samples collected at 2 (46 cases) and 3–5 time points (23 cases) were included. Relationships between BAFF, anti-Jo-1, disease activity, CRP, and their longitudinal changes were evaluated using correlation analysis, multiple regression (MR), path analysis (PA), and hierarchical linear models (HLM).
Cross-sectional assessment demonstrated significant correlations between the levels of BAFF and anti-Jo-1 antibodies which were associated with levels of CK, myoglobin, AST, and CRP, as well as multivariate associations between BAFF, anti-Jo-1 antibodies, and CK levels. PA revealed direct effects of anti-Jo-1 antibodies on CK (β = 0.41) and both direct (β = 0.42) and indirect (through anti-Jo-1 antibodies; β = 0.17) effects of BAFF on CK. Changes in levels of both BAFF and anti-Jo-1 between two time points (Δ) were associated with Δmyoglobin and Δaminotransferases and changes of BAFF correlated with ΔCK, Δcutaneous, Δmuscle, Δglobal, and Δskeletal disease activities.
The longitudinal analysis showed a high intra-individual variability of serum levels of BAFF over time (97%) which could predict 79% of the variance in anti-Jo-1 levels. The anti-Jo-1 variability was explained by inter-individual differences (68%). The close longitudinal relationship between levels of BAFF, anti-Jo-1, and disease activity was supported by high proportions of their variance explained with serum levels of CK and CRP or pulmonary and muscle activities.
Our findings of associations between levels of BAFF and anti-Jo-1 antibodies in serum and myositis activity suggest a role of this cytokine in disease-specific autoantibody production as part of disease mechanisms, and support BAFF as a potential target for intervention in anti-Jo-1-positive myositis patients.
Polymyositis (PM) and dermatomyositis (DM) are chronic, inflammatory disorders characterised by muscle weakness and by the presence of inflammatory infiltrates in the skeletal muscle . Other organs such as the skin and lungs are frequently involved. Myositis-specific antibodies (MSA) or myositis-associated antibodies are present in up to 80% of PM/DM patients . The anti-histidyl-tRNA synthetase (anti-Jo-1) autoantibodies are the most frequent MSA (present in 20–30% of DM/PM patients)  and are associated with a distinct clinical phenotype (i.e. anti-synthetase syndrome), characterised by myositis, Raynaud’s phenomenon, interstitial lung disease (ILD), arthritis, and skin changes of the hands . The observation that anti-Jo-1 antibodies could be present before the onset of clinical symptoms may suggest a possible role of the antibody in the pathogenesis of this subset of myositis [5, 6].
B-cell activating factor of the tumour necrosis factor family (BAFF; also known as B lymphocyte stimulator, or BLyS) is crucial for B-cell maturation and survival. BAFF is also believed to play a role in autoantibody production . High serum levels of BAFF have been reported in patients with autoimmune diseases [8,9,10,11,12,13,14,15] and are associated with disease activity and the presence or levels of autoantibodies [9, 10, 12,13,14, 16]. In patients with myositis, elevated serum levels of BAFF were found to correlate with serum levels of creatine kinase (CK). In addition, the linear regression analysis of variance (ANOVA) model confirmed that anti-Jo-1 antibodies and ILD are the main influencing factors for levels of BAFF, particularly in PM patients . Based on this observation and on the previously reported correlations between serum levels of anti-Jo-1 antibodies and disease activity , we aimed to study associations between BAFF and anti-Jo-1 antibody levels in longitudinally collected serum samples and their relation to standardised clinical measures and laboratory markers of disease activity in patients with DM/PM, with a focus on early cases and the subgroup defined by the presence of ILD. C-reactive protein (CRP) is currently not used as a biomarker of disease activity in myositis. However, correlations between serum levels of CRP and lung function tests among anti-Jo-1-positive myositis patients with ILD have been reported , and elevated CRP levels have been described as a risk factor for developing ILD in patients with myositis [20, 21]. Moreover, an association between increased BAFF and CRP levels in the serum was reported in patients with systemic lupus erythematosus (SLE) [13, 22]. Therefore, we also included analysis of associations with serum levels of CRP.
Patients and controls
All patients with myositis who had tested positive for anti-Jo-1 antibodies over a period of 9 years in the Institute of Rheumatology Prague, Czech Republic (n = 53), and ten patients from the Rheumatology Clinic at Karolinska University Hospital, Stockholm, Sweden, with available longitudinally collected clinical data and paired serum samples were included in the study. Of these, 41 patients were diagnosed as probable or definite PM and 22 as DM based on the Bohan and Peter criteria [23, 24] and 48 of them (76%) had myositis-associated ILD defined by pathological chest x-ray or computed tomography (CT), and abnormal pulmonary function tests.
Clinical data and serum samples were collected prospectively at consecutive clinical visits over 9 years using a disease register in both hospitals. The prospectively collected disease activity assessments were available in 58 patients (20 DM and 38 PM) in the form of the Myositis Disease Activity Assessment Tool (MDAAT) according to the International Myositis Assessment and Clinical Studies (IMACS) Group, including extramuscular, muscular, and physician’s global score of disease activity by visual analogue scale (VAS) and MYOACT score . The Swedish and Czech cohorts had comparable demographic and clinical characteristics, but Swedish patients were older and had longer disease duration (data not shown). Sera from 41 age- (mean 51.1 ± 11.7 years) and gender (female:male = 28:13)-comparable Czech healthy individuals without any known inflammatory disease or recent infection were used as controls. Ethics committees of both institutions approved the study and informed consent was obtained from the participants.
Altogether, 143 serum samples were collected within the study. The earliest available serum sample obtained from each patient was used for the cross-sectional analysis. Paired samples obtained during regular clinic appointments (median 1.5 years, range 0.4 – 7.2) with clinical activity assessment by MDAAT at both visits were available in 46 cases (36 Czech and 10 Swedish, 17 DM and 29 PM patients). Sera from three to five time points collected from 23 patients were included to analyse longitudinal associations between BAFF and anti-Jo-1 antibody levels and their relation with disease activity. Intervals between the first and last sample ranged from 0.6 to 8.9 years and the time flow of individual visits is represented in Additional file 1: Figure S1.
Serum levels of CK, myoglobin, alanine aminotransferase (ALT), aspartate aminotransferase (AST), and CRP were routinely measured in local laboratories with comparable reference levels for Czech and Swedish cohorts (Table 1). All other laboratory analyses were performed in the same laboratory at the Institute of Rheumatology in Prague in serum samples stored at −80 °C until analysis.
The levels of BAFF were measured in serum samples by an enzyme-linked immunosorbent assay (ELISA) according to the manufacturer’s instructions (R&D Systems, Inc., Minneapolis, MN, USA). The cut-off was defined as the mean plus two standard deviations of the control subjects.
Anti-Jo-1 levels were measured in the same patient’s serum samples with the use of ELISA according to the manufacturer’s instructions (Orgentec, Mainz, Germany). The declared cut-off range for positivity was < 15 kU/l with borderline values 15–25 kU/l. To confirm anti-Jo-1 positivity, patient samples were tested using myositis Western blot (Anti-Myositis–Antigen EUROLINE-WB, Euroimmun, Lubeck, Germany) and line blot assay (Myositis-LIA, IMTEC, Berlin, Germany).
The ELISA test revealed levels of anti-Jo-1 antibodies above the upper detection limit of the assay in six out of 143 serum samples that were not available for further titration; therefore, they were excluded from parametric statistical analyses. Six patients had initial levels of anti-Jo-1 antibodies below the cut-off 15 kU/l, but samples were found to be confirmed positive using line blot and/or Western blot assays.
Statistical analyses were performed using GraphPad Prism 5 (GraphPad Software, Inc., San Diego, CA, USA) and SPSS 17.0 (SPSS, Inc., Chicago, IL, USA). With respect to non-normal distribution of data, non-parametric tests were used or the logarithmical transformation to normality was applied where appropriate. For analysis of differences between groups, Mann-Whitney U test and for comparison of changes between two time-points Wilcoxon’s signed rank test were performed. Spearman’s rank order test (correlation coefficient rho abbreviated here as r) was used for correlations of parameters in cross-sectional evaluation and changes between two time points (Δ; computed by subtraction of values at the second visit from values obtained at first visit). Contingency tables were evaluated with Fisher’s exact test. The p values below 0.05 were considered as statistically significant. The Bonferroni correction of alpha values was introduced where appropriate. Based on the distribution, data are presented as median (range; minimum–maximum) values or (for age) as mean ± SD.
Multiple regression (MR) of serum levels of BAFF, anti-Jo-1 antibodies, and CRP (as independent variables), and CK (as dependent variable) was performed using SPSS package 17.0 (SPSS, Inc., Chicago, Illinois, USA). The path analysis (PA) of these variables was performed with use of the software package LISREL 8.8 (Scientific Software International, Inc., Lincolnwood, IL, USA)  on a covariance matrix derived from the logarithmically transformed variables with function ln(1 + x), due to their positively skewed distributions . The maximum likelihood (ML) method was selected to estimate correct standard errors and the associated goodness-of-fit statistics [28, 29]. Models and their variants were evaluated by several fit indices: the chi-square statistics, root mean square error of approximation (RMSEA), comparative fit index (CFI), and parsimony normed fit index (PNFI). RMSEA < 0.05 and CFI > 0.95 were considered as acceptable model fit. Lower values of chi-square statistics and higher values of PNFI also indicate a good fit. Only significant standardised path coefficients (β) were accepted.
Analysis of longitudinal data
None of the traditional ANOVA-like models is appropriate for the analysis of repeated measures over time of the unbalanced design of longitudinal collection of samples (as illustrated in Additional file 1: Figure S1). Therefore, the two-level hierarchical linear model (HLM) was designed separately for BAFF, anti-Jo-1, and CK using the HLM 6 package . The principles of HLM are explained in detail in Additional file 2 .
A sequence of three consecutive models was accomplished. Firstly, the unconditional model was used for evaluation of the variance components situated within patients (lower level 1; 3 and more time points) and between patients (higher level 2; 23 cases) . Their ratio, the intra-class correlation coefficient (ICC), described the proportion of the total variance between patients. Next, the unconditional growth model with time as the independent variable was applied as an indicator of intensity of associations over time. Finally, the model with one time-varying independent variable was used to assess the percentage of the variance of either BAFF, anti-Jo-1, or CK explained by a relevant variable measured repeatedly over time along with a given dependent variable (pseudo R2).
Since the limited sample size of longitudinally collected samples entered the HLM, only few variables representing the hypothesis were included to keep it in proportion with the sample size and to let the path-analytic model and HLM to be identified.
Patient group characteristics
Patient demographic, laboratory, and clinical characteristics at initial assessment are summarised in Table 1. DM and PM cohorts had comparable disease activity (except for cutaneous VAS), descriptive characteristics, and frequency of treatment with glucocorticoids (GC) or disease-modifying anti-rheumatic drugs (DMARDs), with comparable daily doses of GC. Patients with ILD had higher pulmonary activity, extramuscular global assessment and global disease activity VAS, and MYOACT score compared with non-ILD patients (median values in ILD: 20, 19, 19 mm and 0.75 index; non-ILD: 0, 4, 6.5 mm and 0.27 index; p < 0.05 for all) but otherwise comparable disease activity, laboratory parameters, and clinical characteristics. A shorter duration of symptoms (p = 0.04) and time from diagnosis (p = 0.006) at inclusion was seen in patients with ILD.
Early cases (up to 6 months after diagnosis of myositis) constituted 48% of the total cohort (n = 30). These had higher disease activity as assessed by MYOACT score (median 0.87 vs 0.31), extramuscular global assessment, global, muscle and pulmonary activity VAS or serum levels of ALT (23.0, 25.0, 16.5, 23.0 mm and 0.57 μkat/l) when compared with patients with longer disease duration (7.5, 8.0, 6.5, 6.5 mm, 0.35 μkat/l; p < 0.05 for all after Bonferroni correction).
Serum levels of BAFF and anti-Jo-1 antibodies and their mutual correlation at initial evaluation
Patients with PM/DM had higher serum levels of BAFF (median 1.8, range 0.3–20.9 ng/ml) compared with healthy individuals (0.8, 0.4–2.0 ng/ml, p < 0.0001), and elevated levels were documented in 54% of patients.
The serum levels of anti-Jo-1 antibodies had a high variability between individual patients with a wide range (Table 1) and the distribution was skewed towards lower values.
Serum levels of both BAFF and anti-Jo-1 were comparable between subgroups of patients with DM and PM (Table 1; BAFF, p = 0.30; anti-Jo-1, p = 0.13), and with or without ILD (median and range of BAFF: ILD = 1.8, 0.5–20.9, and non-ILD = 1.5, 0.3–8.6 ng/ml; p = 0.56; or anti-Jo-1: ILD = 173.6, 0.8–3604, and non-ILD = 151.8, 0.55–2063 kU/l; p = 0.36).
Positive correlations were found in the whole patient cohort between serum levels of BAFF and anti-Jo-1 antibodies (r = 0.42, p = 0.0006, n = 63), as well as in subgroups with DM or PM and with or without ILD (Table 2).
Associations of serum levels of BAFF and anti-Jo-1 antibodies with laboratory and clinical parameters of disease activity and CRP at initial evaluation
Serum levels of BAFF or anti-Jo-1 were positively associated with serum markers of muscle involvement (levels of CK, myoglobin, and AST; Table 2 and Additional file 3), particularly in patients with early disease.
When analysing the disease activity measures, a positive correlation was found between levels of BAFF and cutaneous activity (Spearman’s correlation coefficient r = 0.46, p = 0.04, n = 20) in patients with DM.
A significant correlation of serum levels of CRP with levels of BAFF (p = 0.0004) and anti-Jo-1 antibodies (p = 0.003; Table 2) was recorded. Therefore, we incorporated associations of CRP into further analyses. The baseline levels of CRP in the serum were significantly higher in patients compared with healthy individuals (median 3.25 mg/l, range 0.2–54.8 mg/l vs 1.3 mg/l, 0.4–3.5 mg/l; p < 0.001). Elevation of CRP levels above the reference limit (5 mg/l) occurred in 38% of the entire patient group and in 47% of patients with early disease. CRP levels were associated with muscle involvement (serum levels of CK (r = 0.30; p = 0.02, n = 60) and myoglobin (r = 0.42, p = 0.002, n = 60)0 as analysed in the whole group of patients and in subsets with PM or with ILD (CK: r = 0.42 or 0.30, p = 0.009 or 0.045, n = 38 or 45; myoglobin: r = 0.53 or 0.40; p = 0.008 or 0.01, n = 33 or 40, respectively). No correlation was found between CRP levels and clinical assessments of disease activity, including constitutional and skeletal VAS.
Multivariate analysis of the cross-sectional relationships between serum levels of BAFF, anti-Jo-1 antibodies, CK, and CRP
The associations among bivariate correlated variables were analysed with multiple regression in the entire cohort of patients (n = 63). Serum levels of BAFF, anti-Jo-1 antibodies, and CRP were included as independent variables, and serum levels of CK as a dependent variable. This analysis showed significant associations both for serum levels of CK and BAFF (with path coefficient β = 0.38; p = 0.005) and levels of CK and anti-Jo-1 antibodies (β = 0.28; p = 0.02).
The possible causal relationships among these variables were further evaluated with the path analysis (PA) as shown by the path diagram (Fig. 1a). The modification of the MR model (explained in detail in Additional file 4) resulted in the PA model (Fig. 1b) where the dependence of anti-Jo-1 levels on BAFF levels (β = 0.42) appeared. Also, the positive effect of serum levels of anti-Jo-1 antibodies (β = 0.41) and BAFF on serum CK levels was recorded. The effect of BAFF on CK was both direct (β = 0.29) and indirect through anti-Jo-1 antibodies (following multiplicative rules: β = 0.42 × 0.41 = 0.17). The effect of CRP on CK was only indirectly mediated by BAFF and anti-Jo-1, but was significant (β = (0.52 × 0.29) + (0.52 × 0.42 × 0.41) = 0.24; p = 0.001). Modification of this model for subgroups of patients with DM and PM also supported direct effects of BAFF and anti-Jo-1 antibodies on CK and of CRP to BAFF (data not shown).
Alternative models testing whether CRP levels could result from the muscular or lung activity did not exhibit acceptable fit.
Time changes in serum levels of BAFF, anti-Jo-1 antibodies, and disease activity, their mutual associations, and effect of treatment with glucocorticoids
A time-related decrease in disease activity and BAFF levels was documented by a negative correlation of time elapsed from diagnosis until sampling with laboratory markers of muscle involvement (CK, myoglobin, LDH, ALT, and AST; r = −0.32, −0.32, −0.30, −0.40, and −0.32, respectively; p < 0.05 for all; n = 63, 53, 55, 62, and 61, respectively) and clinical disease activity (cutaneous, skeletal, pulmonary, extramuscular global, muscle, physician’s global disease activity and MYOACT Score; r = −0.34, −0.36, −0.41, −0.62, −0.40, −0.65 and, −0.55 respectively; p < 0.01 and n = 58 for all). The negative association of serum levels of BAFF (r = −0.55, p = 0.002, n = 30) as well as marginally significant negative association of anti-Jo-1 levels (r = −0.35, p = 0.06, n = 30) with disease duration was documented in early patients.
Between the first two samplings, a significant reduction in myoglobin (median from 94 to 62.8 μg/l; p = 0.001) and a borderline significant drop in CK (from 2.62 to 1.22 μcat/l; p = 0.07) in the total group was recorded, as well as a significant decrease in other muscle enzymes (LDH, ALT, and AST; p < 0.05) in the early patient group. The clinical activity decreased significantly (p < 0.05 for all) in cutaneous, pulmonary, global extraskeletal muscle, muscle, and physician’s global disease activity between assessments at the first (VAS medians 1, 12.5, 15, 11, and 17.5 mm) and second (0, 4, 5, 6, and 11 mm) time points.
The significant reduction in BAFF levels between the two time points was documented within the entire patient group (from 1.8 to 1.3 ng/ml; p = 0.04) and was more pronounced in early patients (from 1.6 to 1.05 ng/ml; p = 0.02). In total, serum BAFF levels decreased to below the cut-off in 37.5% of patients with high BAFF, but an increase above the cut-off was recorded in 27% of patients with initially low BAFF levels. This non-linear variability of BAFF levels over time was also seen in a growth model as described below.
No significant change in anti-Jo-1 antibody levels between two time points was found (from 153.9 to 102.8 kU/l, p = 0.93). The serum levels of anti-Jo-1 antibodies detected in the first and second samplings were closely correlated (r = 0.70, p < 0.0001, n = 46), which corresponds to lower variability of anti-Jo-1 with time.
We further analysed correlations between changes in serum levels of BAFF or anti-Jo-1 and changes in laboratory markers and clinical disease activity in individual patients (Δ = value at first visit minus value at second visit). The changes in BAFF were associated with ΔCK, Δmyoglobin, ΔALT, and ΔAST, and showed a trend to association with Δcutaneous, Δmuscle, Δglobal, and Δskeletal disease activities (Table 3 and Additional file 5, left panel). The Δanti-Jo-1 antibody levels correlated only with the change in laboratory markers ΔAST and marginally with ΔALT, Δmyoglobin, and ΔCRP (Table 3 and Additional file 5, centre panel).
Changes in BAFF were associated with the daily dose of GC used at the time of first blood sampling. The change in disease activity assessed by Δmyoglobin, ΔALT, ΔAST, Δpulmonary, and Δglobal VAS also reflected the effect of GC therapy (Table 3, right column). Thus, higher doses of GC were followed by a larger decrease of serum levels of BAFF and with a decrease in some indicators of disease activity, but not with change in the anti-Jo-1 levels (r = 0.10, p = 0.54).
Finally, the changes in levels of anti-Jo-1 and BAFF were not mutually significantly associated, with a trend to improve after the exclusion of a single outlying value of Δanti-Jo-1 (Additional file 5).
The proportion of patients with high serum CRP decreased from an initial 38% to 28% at the second sampling, but no significant quantitative decrease in CRP levels was recorded. However, the changes in CRP were associated with changes in BAFF and anti-Jo-1 (Table 3), Δmyoglobin (r = 0.33, p < 0.05, n = 36) and weakly with Δconstitutional disease activity (r = 0.31, p = 0.07, n = 36).
Longitudinal relations between serum levels of BAFF or anti-Jo-1 antibodies and disease activity markers or CRP
The time variability in serum levels of BAFF, anti-Jo-1, and CK showed comparable trends within time in a majority of patients who were followed from the early period after the diagnosis (illustrated in Fig. 2 with diagrams of 12 patients who had samples from three and more time points and short disease duration at first sampling; 1.0 ± 1.4 months after diagnosis determination). In general, high levels of serum BAFF declined after initial treatment with a high dose of GC, tended to rise in individual cases after GC tapering, and were parallel to CK levels more closely than anti-Jo-1 antibody levels (Fig. 2).
Anti-Jo-1 antibody levels became negative during follow-up in seven out of 40 patients (17.5%) within 0.5–4.3 years, in association with a significant decrease in serum levels of BAFF (median from 1.9 to 0.7 ng/ml; p = 0.02) and a decrease in disease activity measured as CK, myoglobin, ALT, and AST (p < 0.05). These patients were treated with a combination of GC (median daily dose 27.5 mg) and DMARDs (3 × azathioprine, 2 × methotrexate, 1 × cyclosporine A, and 1 × cyclophosphamide).
At a group level, statistical analysis of repeated measures with HLM showed high intra-individual variability in serum levels of BAFF over time (Table 4). This suggests that the majority (97%) of the variability of BAFF levels depends on time-related changes within the individual patient (ICC = 0.03), whereas 68% of the variability of anti-Jo-1 levels was explained by inter-individual differences between patients (ICC = 0.678).
The percentage of explained variance within patients in the unconditional growth model with time from diagnosis as an independent variable (Table 4) showed that linear regression is a good model for anti-Jo-1 (62%) and marginally acceptable for CK (31%), but not appropriate for BAFF (0%) levels in serum. BAFF was a good predictor of within-patient variability of anti-Jo-1 (it explained 79% of anti-Jo-1 variance) but the opposite prediction was not found, which corresponds to the dependence of anti-Jo-1 levels on BAFF detected by the PA model as described above (Fig. 1b).
The close longitudinal associations between serum levels of BAFF or anti-Jo-1 antibodies with markers of muscle involvement and clinical disease activity were supported, with a high proportion of their variance explained by CK, CRP, pulmonary, and muscle VAS (BAFF: 76%, 74%, 61%, and 85%; or anti-Jo-1: 74%, 70%, 60%, and 72%). Moreover, anti-Jo-1 antibody levels could be predicted with global disease activity (71%).
The extended interpretation of the further results of HLM is provided in detail in Additional file 2.
The main finding of this study was a strong correlation between serum levels of BAFF and anti-Jo-1 antibodies and their association with disease activity variables and clinical phenotype of myositis-associated ILD. In addition, we characterised differences in time-related variability of BAFF and anti-Jo-1 serum levels during longitudinal follow-up. Whereas serum levels of BAFF were associated with treatment within the same individual, the levels of anti-Jo-1 were more stable and did not reflect currently used therapy. Finally, we recorded associations of the serum levels of CRP with BAFF, anti-Jo-1, and markers of muscle impairment.
Higher serum levels of BAFF compared with healthy individuals confirm previous reports in unselected cohorts of patients with myositis [17, 33,34,35]. Here, we analysed the whole accessible cohort of anti-Jo-1-positive patients with myositis. The higher proportion of PM patients compared with DM is in accordance with the known higher frequency of anti-Jo-1 in PM , and the high proportion of patients with ILD compared with the overall myositis population  corresponds with the strong association of ILD with anti-Jo-1 positivity [19, 38]. The numerically more frequent ILD in patients with PM compared with those with DM is in agreement with already published data and, again, can be explained by the anti-Jo-1 status .
Contrary to earlier reports in autoantibody heterogeneous  or not defined  cohorts of myositis patients, we did not find a difference in serum levels of BAFF between patients with DM and PM subgroups or higher BAFF levels in patients with ILD within the current autoantibody homogenous, anti-Jo-1-positive cohort. However, the present findings confirm our previous results of multiple regression analysis showing that diagnosis of PM or DM was less discriminatory for BAFF levels than the presence of anti-Jo-1 antibodies or lung involvement, and explaining higher serum BAFF levels in myositis-associated ILD by the presence of anti-Jo-1 antibodies . We also found a close correlation between serum levels of BAFF and anti-Jo-1, particularly in early cases and in patients with PM or in PM/DM patients with ILD. A similar association between serum levels of BAFF and other disease-specific autoantibodies has been reported, for example anti-dsDNA or anti-Smith antibodies in SLE [9, 10, 13, 16, 22, 40, 41], rheumatoid factor (RF) or anti-CCP antibodies in rheumatoid arthritis (RA) [10, 42], and RF or anti-SSA/B antibodies in Sjögren’s syndrome , suggesting a close relationship between BAFF and autoantibody production. A biologic association between BAFF and anti-Jo-1 antibody levels was supported by our path analysis which revealed a dependence of anti-Jo-1 levels on BAFF levels in serum. A possible causal relationship was further corroborated by longitudinal analysis with the HLM, where a high proportion of anti-Jo-1 variance was explained by serum levels of BAFF, but the converse was not true. This analysis, together with the reported quantitative association of BAFF with another MSA anti-MDA5 (melanoma differentiation-associated gene 5) in juvenile DM , further supports involvement of a disease-specific inflammatory process in the production of autoantibodies triggered by BAFF release into the circulation in inflammatory myopathies.
Studies from other systemic diseases (SLE and RA) reported limited variation of serum levels of BAFF over time and no relation to disease flares [9, 13]. However, our findings of changes in serum levels of BAFF over time in association with changes in disease activity was similar to a longitudinal follow-up in systemic sclerosis (SSc), early RA, and SLE [8, 41, 42], indicating that BAFF may have a higher impact in early disease.
Only 54% of our anti-Jo-1-positive patients had elevated baseline serum levels of BAFF, and the levels varied considerably in association with the laboratory markers of muscle impairment, which confirms previous reports [17, 34]. These correlations were expressed particularly early after the diagnosis and in subgroups of patients with PM or ILD, further supporting the hypothesis that BAFF is more involved in early phases of this disease and that other factors may also influence autoantibody production.
Moreover, serum levels of BAFF were associated with cutaneous disease activity in DM patients. An association of serum BAFF levels with severity of skin sclerosis in SSc, with higher numbers of plaque lesions in localised scleroderma and with malar rash in SLE [8, 12, 40], has been reported previously. All these associations point to a possibility that active skin lesions might be a source of BAFF.
The main variability of serum levels of BAFF was seen in time-related changes within individual patients and a substantial number of patients expressed reduction from high to low levels between two time points, but the opposite, an increase in BAFF levels, also occurred. These changes correlated with changes in markers of muscle impairment and also with changes in cutaneous, muscular, skeletal, and global disease activity. The close longitudinal relationship between disease activity and serum BAFF levels was further supported by a high proportion of its variance being explained by serum CK levels and muscle or pulmonary VAS in analysis with HLM. This, together with the position of BAFF in the PA model, may suggest BAFF as a possible direct mediator of disease mechanisms in patients with anti-Jo-1-positive myositis.
The changes between two samplings in serum levels of BAFF and laboratory and clinical disease activity measures showed a trend towards an association with initial daily glucocorticoid doses. Similar associations were also seen in patients with SLE or SSc [9, 12]. We have previously reported an association of an early decrease of BAFF levels with cumulative glucocorticoid doses at the beginning of treatment in early myositis cases , and our current findings provide further support for the sensitivity of BAFF production to glucocorticoid treatment in myositis.
The anti-Jo-1 antibody became negative in patients with a decrease of disease activity (in 18%), which was preceded by intensive treatment and associated with a decrease in serum BAFF levels. This is similar to previous reports of the disappearance of anti Jo-1 autoantibodies in periods of disease inactivity [18, 43, 44]. This could also be an explanation for the low anti-Jo-1 antibody levels in six already treated patients who had positive anti-Jo-1 levels in sera taken before the study start.
Similar to BAFF, the serum levels of anti-Jo-1 antibodies correlated with serum markers of muscle impairment. The results of our unconditional growth model demonstrated associations between variations in anti-Jo-1 antibody levels over time and with muscular, pulmonary, and global disease activity comparable to longitudinal observations analysed with mixed linear regression for repeated measures reported by Stone et al. , and supports the role of anti-Jo-1 in disease mechanisms.
An unexpected finding in our study was the association of serum levels of CRP with levels of BAFF and anti-Jo-1 antibodies. Besides the association with change of constitutional activity, CRP also reflected muscle impairment by correlation with CK and myoglobin in serum both at baseline and longitudinal follow-up. These associations were recorded despite the modest elevations in CRP levels. Although CRP is not considered to be a marker of disease activity in most systemic diseases, some studies indicates its role mainly in patients with lung involvement , including myositis [19,20,21, 46]. We recorded more pronounced associations of muscle impairment markers with CRP in a subset of patients with ILD, but no correlation of CRP levels in serum with pulmonary activity. In our PA model, CRP displayed significant effects on serum levels of BAFF, similar to that reported in SLE (Fig. 1b) [13, 22, 47]. This finding could be explained by the described capability of CRP to trigger the release of soluble, biologically active BAFF from myeloid cells by FcγΙa receptor engagement . These mechanisms, however, need further exploration.
A limitation of our study is the wide heterogeneity of disease activity with low median values (Table 1), which may explain the absence of associations between serum levels of BAFF and clinical activity measures or anti-Jo-1 at baseline, as seen by other authors [18, 34]. The low disease activity in some patients could possibly be explained by the intensive treatment that patients received before inclusion in our study.
In summary, levels of anti-Jo-1 in serum are proposed as a biomarker of disease activity in patients with anti-Jo-1-positive myositis. Our study demonstrated that serum levels of BAFF correlate with serum levels of anti-Jo-1 antibodies. Moreover, the longitudinal analysis showed a high intra-individual variability of serum levels of BAFF over time which could predict 79% of anti-Jo-1 level variance, whereas the anti-Jo-1 variability was explained by inter-individual differences. A relationship between serum levels of BAFF, anti-Jo-1, and disease activity was supported by high proportions of their variance being explained by serum levels of CK and CRP or pulmonary and muscle activities. The finding of an association between serum levels of CRP and levels of BAFF, anti-Jo-1 antibodies, or markers of muscle impairment was surprising, although consistent with some reports in other systemic diseases.
Associations between serum levels of BAFF and anti-Jo-1 antibodies with disease activity measures support a possible role of BAFF in the pathogenesis of anti-Jo-1-positive myositis. The anti-Jo-1 autoantibody was recently described as one of the strongest predictors of response in rituximab-treated myositis patients in the RIM trial . Our finding of an association of anti-Jo-1 levels with BAFF, particularly in patients with PM and ILD, also implies that BAFF-blocking therapy could be an attractive novel treatment for these patients.
Analysis of variance
B-cell activating factor of the tumour necrosis factor family
Comparative fit index
Disease-modifying anti-rheumatic drug
Enzyme-linked immunosorbent assay
Hierarchical linear model
Intra-class correlation coefficient
Interstitial lung disease
International Myositis Assessment and Clinical Studies
Myositis Disease Activity Assessment Tool
Parsimony normed fit index
Root mean square error of approximation
Systemic lupus erythematosus
Visual analogue scale
Plotz PH, Rider LG, Targoff IN, Raben N, O'Hanlon TP, Miller FW. NIH conference. Myositis: immunologic contributions to understanding cause, pathogenesis, and therapy. Ann Intern Med. 1995;122(9):715–24.
Gunawardena H, Betteridge ZE, McHugh NJ. Myositis-specific autoantibodies: their clinical and pathogenic significance in disease expression. Rheumatology (Oxford). 2009;48(6):607–12.
Ronnelid J, Barbasso Helmers S, Storfors H, Grip K, Ronnblom L, Franck-Larsson K, Nordmark G, Lundberg IE. Use of a commercial line blot assay as a screening test for autoantibodies in inflammatory myopathies. Autoimmun Rev. 2009;9(1):58–61.
Love LA, Leff RL, Fraser DD, Targoff IN, Dalakas M, Plotz PH, Miller FW. A new approach to the classification of idiopathic inflammatory myopathy: myositis-specific autoantibodies define useful homogeneous patient groups. Medicine. 1991;70(6):360–74.
Yoshifuji H, Fujii T, Kobayashi S, Imura Y, Fujita Y, Kawabata D, Usui T, Tanaka M, Nagai S, Umehara H, et al. Anti-aminoacyl-tRNA synthetase antibodies in clinical course prediction of interstitial lung disease complicated with idiopathic inflammatory myopathies. Autoimmunity. 2006;39(3):233–41.
Matsushita T, Hasegawa M, Fujimoto M, Hamaguchi Y, Komura K, Hirano T, Horikawa M, Kondo M, Orito H, Kaji K, et al. Clinical evaluation of anti-aminoacyl tRNA synthetase antibodies in Japanese patients with dermatomyositis. J Rheumatol. 2007;34(5):1012–8.
Mackay F, Schneider P, Rennert P, Browning J. BAFF AND APRIL: a tutorial on B cell survival. Annu Rev Immunol. 2003;21:231–64.
Matsushita T, Hasegawa M, Matsushita Y, Echigo T, Wayaku T, Horikawa M, Ogawa F, Takehara K, Sato S. Elevated serum BAFF levels in patients with localized scleroderma in contrast to other organ-specific autoimmune diseases. Exp Dermatol. 2007;16(2):87–93.
Stohl W, Metyas S, Tan SM, Cheema GS, Oamar B, Xu D, Roschke V, Wu Y, Baker KP, Hilbert DM. B lymphocyte stimulator overexpression in patients with systemic lupus erythematosus: longitudinal observations. Arthritis Rheum. 2003;48(12):3475–86.
Pers JO, Daridon C, Devauchelle V, Jousse S, Saraux A, Jamin C, Youinou P. BAFF overexpression is associated with autoantibody production in autoimmune diseases. Ann N Y Acad Sci. 2005;1050:34–9.
Vallerskog T, Heimburger M, Gunnarsson I, Zhou W, Wahren-Herlenius M, Trollmo C, Malmstrom V. Differential effects on BAFF and APRIL levels in rituximab-treated patients with systemic lupus erythematosus and rheumatoid arthritis. Arthritis Res Ther. 2006;8(6):R167.
Matsushita T, Hasegawa M, Yanaba K, Kodera M, Takehara K, Sato S. Elevated serum BAFF levels in patients with systemic sclerosis: enhanced BAFF signaling in systemic sclerosis B lymphocytes. Arthritis Rheum. 2006;54(1):192–201.
Becker-Merok A, Nikolaisen C, Nossent HC. B-lymphocyte activating factor in systemic lupus erythematosus and rheumatoid arthritis in relation to autoantibody levels, disease measures and time. Lupus. 2006;15(9):570–6.
Mariette X, Roux S, Zhang J, Bengoufa D, Lavie F, Zhou T, Kimberly R. The level of BLyS (BAFF) correlates with the titre of autoantibodies in human Sjogren’s syndrome. Ann Rheum Dis. 2003;62(2):168–71.
Szodoray P, Alex P, Jonsson MV, Knowlton N, Dozmorov I, Nakken B, Delaleu N, Jonsson R, Centola M. Distinct profiles of Sjogren’s syndrome patients with ectopic salivary gland germinal centers revealed by serum cytokines and BAFF. Clin Immunol. 2005;117(2):168–76.
Cheema GS, Roschke V, Hilbert DM, Stohl W. Elevated serum B lymphocyte stimulator levels in patients with systemic immune-based rheumatic diseases. Arthritis Rheum. 2001;44(6):1313–9.
Krystufkova O, Vallerskog T, Helmers SB, Mann H, Putova I, Belacek J, Malmstrom V, Trollmo C, Vencovsky J, Lundberg IE. Increased serum levels of B cell activating factor (BAFF) in subsets of patients with idiopathic inflammatory myopathies. Ann Rheum Dis. 2009;68(6):836–43.
Stone KB, Oddis CV, Fertig N, Katsumata Y, Lucas M, Vogt M, Domsic R, Ascherman DP. anti-Jo-1 antibody levels correlate with disease activity in idiopathic inflammatory myopathy. Arthritis Rheum. 2007;56(9):3125–31.
Richards TJ, Eggebeen A, Gibson K, Yousem S, Fuhrman C, Gochuico BR, Fertig N, Oddis CV, Kaminski N, Rosas IO, et al. Characterization and peripheral blood biomarker assessment of anti-Jo-1 antibody-positive interstitial lung disease. Arthritis Rheum. 2009;60(7):2183–92.
Cen X, Zuo C, Yang M, Yin G, Xie Q. A clinical analysis of risk factors for interstitial lung disease in patients with idiopathic inflammatory myopathy. Clin Dev Immunol. 2013;2013:648570.
Marie I, Hachulla E, Cherin P, Dominique S, Hatron PY, Hellot MF, Devulder B, Herson S, Levesque H, Courtois H. Interstitial lung disease in polymyositis and dermatomyositis. Arthritis Rheum. 2002;47(6):614–22.
Eilertsen GO, Van Ghelue M, Strand H, Nossent JC. Increased levels of BAFF in patients with systemic lupus erythematosus are associated with acute-phase reactants, independent of BAFF genetics: a case-control study. Rheumatology (Oxford). 2011;50(12):2197–205.
Bohan A, Peter JB. Polymyositis and dermatomyositis (second of two parts). N Engl J Med. 1975;292(8):403–7.
Bohan A, Peter JB. Polymyositis and dermatomyositis (first of two parts). N Engl J Med. 1975;292(7):344–7.
Isenberg DA, Allen E, Farewell V, Ehrenstein MR, Hanna MG, Lundberg IE, Oddis C, Pilkington C, Plotz P, Scott D, et al. International consensus outcome measures for patients with idiopathic inflammatory myopathies. Development and initial validation of myositis activity and damage indices in patients with adult onset disease. Rheumatology (Oxford). 2004;43(1):49–54.
Jöreskog KG, Sörbom D. LISREL (version 8.8). Lincolnwood: Scientific Software International, Inc; 2006. p. 22–69.
Montfort K, Mooijaart A, Meijerink F. Estimating structural equation models with nonnormal variables by using transformations. Statistica Nederlandica. 2009;63(2):213–26.
Kaplan DW. Structural equation modeling: foundation and extension. Thousand Oaks: Sage; 2008. p. 13–38.
Kline RB. Principles and practice of structural equation modeling. 2nd ed. New York: Guilford Press; 2005. p. 93–161.
Raudenbush SW, Bryk AS, Cheong YF, Congon R, duToit M. HLM 6: hierarchical linear and nonlinear modeling. Lincolnwood: Scientific Software International, Inc; 2004. p. 14–64.
Hedeker D, Gibbons R. Longitudinal Data Analysis. New Jersey: Wiley; 2006. p. 47.
Singer JD, Willett JB. Applied longitudinal data analysis: modeling change and event occurrence. New York: Oxford University Press; 2003. p. 45–74.
Szodoray P, Alex P, Knowlton N, Centola M, Dozmorov I, Csipo I, Nagy AT, Constantin T, Ponyi A, Nakken B, et al. Idiopathic inflammatory myopathies, signified by distinctive peripheral cytokines, chemokines and the TNF family members B-cell activating factor and a proliferation inducing ligand. Rheumatology (Oxford). 2010;49(10):1867–77.
Peng QL, Shu XM, Wang DX, Wang Y, Lu X, Wang GC. B-cell activating factor as a serological biomarker for polymyositis and dermatomyositis. Biomark Med. 2014;8(3):395–403.
Kobayashi N, Kobayashi I, Mori M, Sato S, Iwata N, Shigemura T, Agematsu K, Yokota S, Koike K. Increased serum B cell activating factor and a proliferation-inducing ligand are associated with interstitial lung disease in patients with juvenile dermatomyositis. J Rheumatol. 2015;42(12):2412–8.
Brouwer R, Hengstman GJ, Vree Egberts W, Ehrfeld H, Bozic B, Ghirardello A, Grondal G, Hietarinta M, Isenberg D, Kalden JR, et al. Autoantibody profiles in the sera of European patients with myositis. Ann Rheum Dis. 2001;60(2):116–23.
Lilleker JB, Vencovsky J, Wang G, Wedderburn LR, Diederichsen LP, Schmidt J, Oakley P, Benveniste O, Danieli MG, Danko K, et al. The EuroMyositis registry: an international collaborative tool to facilitate myositis research. Ann Rheum Dis. 2018;77(1):30–9.
Schmidt WA, Wetzel W, Friedlander R, Lange R, Sorensen HF, Lichey HJ, Genth E, Mierau R, Gromnica-Ihle E. Clinical and serological aspects of patients with anti-Jo-1 antibodies—an evolving spectrum of disease manifestations. Clin Rheumatol. 2000;19(5):371–7.
Fathi M, Dastmalchi M, Rasmussen E, Lundberg IE, Tornling G. Interstitial lung disease, a common manifestation of newly diagnosed polymyositis and dermatomyositis. Ann Rheum Dis. 2004;63(3):297–301.
McCarthy EM, Lee RZ, Ni Gabhann J, Smith S, Cunnane G, Doran MF, Howard D, O'Connell P, Kearns G, Jefferies CA. Elevated B lymphocyte stimulator levels are associated with increased damage in an Irish systemic lupus erythematosus cohort. Rheumatology (Oxford). 2013;52(7):1279–84.
Petri M, Stohl W, Chatham W, McCune WJ, Chevrier M, Ryel J, Recta V, Zhong J, Freimuth W. Association of plasma B lymphocyte stimulator levels and disease activity in systemic lupus erythematosus. Arthritis Rheum. 2008;58(8):2453–9.
Bosello S, Youinou P, Daridon C, Tolusso B, Bendaoud B, Pietrapertosa D, Morelli A, Ferraccioli G. Concentrations of BAFF correlate with autoantibody levels, clinical disease activity, and response to treatment in early rheumatoid arthritis. J Rheumatol. 2008;35(7):1256–64.
Yoshida S, Akizuki M, Mimori T, Yamagata H, Inada S, Homma M. The precipitating antibody to an acidic nuclear protein antigen, the Jo-1, in connective tissue diseases. A marker for a subset of polymyositis with interstitial pulmonary fibrosis. Arthritis and rheumatism. 1983;26(5):604–11.
Miller FW, Twitty SA, Biswas T, Plotz PH. Origin and regulation of a disease-specific autoantibody response. Antigenic epitopes, spectrotype stability, and isotype restriction of anti-Jo-1 autoantibodies. J Clin Invest. 1990;85(2):468–75.
Muangchan C, Harding S, Khimdas S, Bonner A, Canadian Scleroderma Research Group, Baron M, Pope J. Association of C-reactive protein with high disease activity in systemic sclerosis: results from the Canadian Scleroderma Research Group. Arthritis Care Res. 2012;64(9):1405–14.
Xu Y, Yang CS, Li YJ, Liu XD, Wang JN, Zhao Q, Xiao WG, Yang PT. Predictive factors of rapidly progressive-interstitial lung disease in patients with clinically amyopathic dermatomyositis. Clin Rheumatol. 2016;35(1):113–6.
Roth DA, Thompson A, Tang Y, Hammer AE, Molta CT, Gordon D. Elevated BLyS levels in patients with systemic lupus erythematosus: associated factors and responses to belimumab. Lupus. 2016;25(4):346–54.
Li X, Su K, Ji C, Szalai AJ, Wu J, Zhang Y, Zhou T, Kimberly RP, Edberg JC. Immune opsonins modulate BLyS/BAFF release in a receptor-specific fashion. J Immunol. 2008;181(2):1012–8.
Aggarwal R, Bandos A, Reed AM, Ascherman DP, Barohn RJ, Feldman BM, et al. Predictors of clinical improvement in rituximab-treated refractory adult and juvenile dermatomyositis and adult polymyositis. Arthritis Care Res. 2014;66(3):740–9.
The authors are thankful to Johan Rönnelid, M.D., prof., Uppsala University, Uppsala, Sweden, for kind consultations on methodological issues of autoantibody quantification.
This study was supported by grant NT/12438–4 of IGA MZCR of the Project for Conceptual Development of Research Organisation 00023728 from the Ministry of Health in the Czech Republic, the Swedish Research Council K2014-52X-14045-14, and through the regional agreement on medical training and clinical research (ALF) between Stockholm County Council and Karolinska Institutet.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The study was approved by the Ethics Committee of the Institute of Rheumatology in Prague and by the local ethics committee at the Karolinska Hospital Nord, Stockholm. Informed consent to participate was obtained from all participants.
Consent for publication
Not applicable as no individual personal data are included.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Design of data collection. Graphical representation of the time flow of individual patient visits. The blue background highlights cases followed for more than two time points. (PDF 122 kb)
The application method of the hierarchical linear model (HLM). Detailed explanation of HLM principles, its use for the unbalanced designs, and interpretation of results. (PDF 199 kb)
The scatter plots of source cross-sectional data for correlational analysis presented in Table 2. The values of levels of BAFF, anti-Jo-1 antibodies, and CRP in serum are plotted in columns against the serum levels of markers of muscle impairment (CK, myoglobin, and AST) and CRP in rows. Based on the non-normal distribution, the logarithmically transformed data are plotted. Statistics are: r = Spearman’s correlation coefficient; p = p value. (PDF 540 kb)
Multivariate analysis of the cross-sectional relationships between serum levels of BAFF, anti-Jo-1 antibodies, CK, and CRP. The detailed explanation of modifications of the multiple regression model and its variants accomplished in the path analysis and principles of interpretation. (PDF 193 kb)
The scatter plots of source data for correlational analysis presented in Table 3. Changes between first two visits (Δ = 1st visit – 2nd visit) of (A) BAFF plotted against Δanti-Jo-1, changes in both BAFF and anti-Jo-1 plotted in columns against parameters of activity in rows. These are: (B) changes in markers of muscle impairment (ΔCK, Δmyoglobin, ΔALT and ΔAST) and (C) changes in clinical disease activity assessments (Δmuscle, Δglobal, Δskeletal within the entire patient group, Δcutaneous within patients with dermatomyositis (DM), and Δpulmonary within patients with lung involvement (ILD)). Statistics are: r = Spearman’s correlation coefficient; p = p value. A single outlying value of Δanti-Jo-1 is highlighted by a red circle. The graphs with exclusion of the outlier are plotted in the right column. The significance of some correlations became even stronger after exclusion of the outlier. (PDF 535 kb)