Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

IgG4:IgG RNA ratio differentiates active disease from remission in granulomatosis with polyangiitis: a new disease activity marker? A cross-sectional and longitudinal study

Arthritis Research & Therapy volume 21, Article number: 43 (2019) Cite this article

Abstract

Objectives

An important limitation in granulomatosis with polyangiitis (GPA) is the lack of disease activity markers. Immunoglobulin G4-positive (IgG4+) B cells and plasma cells are implicated in the pathogenesis of GPA. We hypothesized that the presence of these cells in peripheral blood could serve as disease activity parameter in GPA.

Methods

We included 35 proteinase 3-antineutrophil cytoplasmic antibodies-positive patients with GPA in a cross-sectional study. Active disease was defined as Birmingham Vasculitis Activity Score (BVAS) ≥ 3 (n = 15), remission as BVAS of 0 (n = 17), and low disease activity (LDA) as BVAS of 1–2 and clinical remission (n = 3). Healthy subjects (n = 10), patients with systemic lupus erythematosus (n = 24), and patients with rheumatoid arthritis (n = 19) functioned as control subjects. An additional longitudinal study was performed in ten patients with GPA. Using a validated qPCR test, we measured the IgG4:IgG RNA ratio in all groups and compared the results with known biomarkers.

Results

The median qPCR score was higher in active GPA (21.4; IQR 12.1–29.6) than in remission/LDA (3.3; IQR 1.6–5.6) (Mann-Whitney U test, p < 0.0001) and outperformed other known disease activity parameters in detecting activity. A cutoff qPCR score of 11.2% differentiated active disease from remission/LDA accurately (AUC 0.993). The qPCR test correlated well with the BVAS (Spearman r = 0.77, p < 0.0001). In the longitudinal study, a decrease in BVAS correlated with qPCR score reduction (paired t test, p < 0.05).

Conclusions

The IgG4:IgG RNA ratio in GPA accurately distinguishes active disease from remission and correlates well with disease activity in these single-center studies. If these results are confirmed in larger longitudinal studies, this test might help to steer treatment decisions in patients with GPA.

Background

Granulomatosis with polyangiitis (GPA), formerly called Wegener’s disease, is a hallmark form of primary small vessel vasculitis. In GPA, the disease activity often fluctuates over time, making titration of therapy challenging [1]. Both disease activity and immunosuppressive therapies cause morbidity and mortality. Although therapeutic options have improved considerably, mortality is still high (15–25% mortality in the first 5 years after diagnosis) [2, 3]. In order to optimize treatment and improve survival, there is a pressing need for reliable biomarkers that reflect disease activity and can help to steer titration of therapy.

Current practice relies on general and unspecific inflammation parameters, such as erythrocyte sedimentation rate (ESR) and C-reactive protein (CRP), that do not correlate well with disease activity in vasculitis. Although GPA is associated with antineutrophil cytoplasmic antibodies (ANCA), these antibodies are more reliable in diagnosing the disease than in monitoring disease activity [4,5,6,7,8,9]. Instruments measuring disease activity, such as the Birmingham Vasculitis Activity Score (BVAS), have provided a solution and have been validated [10, 11]. However, they require that a clinician is convinced that symptoms scored are attributable to vasculitis activity and not to clinical mimics such as infections. This is often challenging in daily clinical practice.

A potential marker might be the presence of immunoglobulin G subtype 4 (IgG4)-positive B cells in peripheral blood. IgG4 itself and IgG4-positive (IgG4+) cells have been implicated in the pathogenesis of GPA. Several studies have reported increased levels of IgG4 antibodies in both serum and tissue in ANCA-associated vasculitis (AAV) [12,13,14,15]. IgG4+ B cells were also detected in inflamed tissue [13, 14, 16, 17].

In IgG4-related disease (IgG4-RD), a multiorgan disease known for high IgG4+ plasma cell infiltration in tissue, we recently demonstrated that the presence of IgG4+ B-cell receptor (BCR) clones in blood can serve as a diagnostic and disease activity biomarker [18]. These clones can be identified not only by next-generation sequencing (NGS) of the BCR repertoire but also by a simplified protocol using a qPCR test. The qPCR test accurately detected the IgG4+ BCR clones by measuring the RNA ratio between IgG4+ BCR and total IgG+ BCR [19].

Given the presence of IgG4+ B-cells in AAV, we hypothesized that the marker mentioned might also be informative in vasculitis as a disease activity marker. To this end, we tested the qPCR test in a cross-sectional and longitudinal cohort of patients with GPA.

Materials and methods

Study subjects

Patients meeting the American College of Rheumatology criteria and the Chapel Hill Consensus Conference definition for GPA were included in the cross-sectional study. In addition to these criteria, elevation in proteinase 3 (PR3)-ANCA had to be present at one point in the patient’s history. The diagnosis and these criteria were confirmed by two independent rheumatologists. Patients with concomitant other inflammatory diseases were excluded.

Patients were defined as having active disease if the BVAS version 3 showed a score of 3 points or higher [20]. Patients with a score of 0 were defined as being in remission. Patients with a BVAS of 1–2 who did not require ongoing therapy to be altered were defined as having low disease activity (LDA). The BVAS results were confirmed by two BVAS-certified clinicians. Both clinicians were blinded for the qPCR result, which was scored independently, and they subsequently reached agreement on the final BVAS. Limited disease, in contrast to systemic disease, was in agreement with present literature and was defined as the presence of manifestations of GPA that posed no immediate threat to either the patient’s life or the function of a vital organ [21]. In addition, activity parameters such as ESR (cutoff 30 mm/h), CRP (cutoff 5 mg/L), PR3-ANCA level (cutoff 7 kIU/L), and physician visual analogue scale (VAS) were determined. Age, sex, medication use, and the number of years since diagnosis were annotated.

Ten patients (seven with active disease) from the cross-sectional cohort were resampled at a later time point in the same manner, thus allowing longitudinal follow-up. One of these patients was resampled twice (Fig. 3, patient A).

As control groups, we included healthy control subjects (HC) and patients diagnosed with other autoimmune diseases with important B-cell involvement: systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA). In addition, individuals who were age-matched to the active GPA cohort were included from our IgG4-RD cohort and used as positive controls [19]. In the SLE and RA groups, disease activity was determined according to current clinical standard using the Systemic Lupus Erythematosus Disease Activity Index (> 3 points) and Disease Activity Score in 28 joints (> 2.6 points) scores, respectively. The IgG4-RD control subjects had untreated disease and had active disease at the time of inclusion. Enrollment of all patients took place in the Academic Medical Center in Amsterdam, the Netherlands.

Peripheral blood sampling and preparation

Peripheral blood was collected in PAXGene blood RNA tubes (catalog 762165; PreAnalytiX, Breda, The Netherlands) according to manufacturer’s instructions. RNA isolation was performed according to the manufacturer’s instructions using RNase-free DNase (Qiagen, Hilden, Germany) in the process. Complementary DNA was synthesized from 1000 ng of total RNA input using SuperScript III Reverse Transcriptase (Life Technologies, Carlsbad, CA, USA).

qPCR

For a detailed description of the qPCR test used in this study, see an earlier publication by our group [19] and Additional file 1: Supplementary methods. Briefly, qPCR uses a common forward primer that binds upstream on the constant region of the BCR heavy chain and two specific reverse primers that bind downstream on the constant region of the BCR heavy chain. One reverse primer binds all possible IgG subtypes, and the other reverse primer is specific for the IgG4 subtype only. Subsequently, the IgG4:IgG ratio can be calculated from the qPCR results.

Data analysis

qPCR data were analyzed using a LightCycler 480 instrument (Roche Diagnostics, Indianapolis, IN, USA) and LinRegPCR [22]. The analysis was always done by two blinded scientists separately. Statistical analyses were performed using Prism (GraphPad Software, La Jolla, CA, USA) and SPSS (IBM, Armonk, NY, USA) software. (Paired) t tests and analysis of variance (ANOVA) were used for normal distributions. Comparison of nonnormal distributed data was performed using the Mann-Whitney U test or Kruskal-Wallis test. Corrections for multiple testing were performed according to the Bonferroni or Dunn method. Similarly, correlation tests were performed using the Pearson method in case of normal distribution and the Spearman method in case of nonnormal distribution of data.

Results

Patient inclusion and characteristics

Thirty-five PR3-ANCA-positive patients with GPA were included. The median age was 56 years. Fifty-four percent of patients were female. The cohort was divided into an active disease group and a remission group on the basis of BVAS. A BVAS of 3 points or higher was regarded as active disease (15 patients, 43%). Seventeen patients (49%) had a BVAS of 0 (remission). Three patients (9%) had a BVAS of 1–2 (LDA). Patients in the LDA group had mild myalgia, arthralgia, and transient bloody nasal discharge that did not require change of treatment. The remission and LDA groups were combined for analysis (remission/LDA group, 20 patients [57%]). Sex, age, and the number of years since diagnosis did not differ significantly between both groups (Table 1). The symptoms scored in the BVAS per active disease and LDA subject are shown in Additional file 1: Table S1.

Table 1 Patient characteristics and disease activity parameters in total granulomatosis with polyangiitis cohort, active patients as determined by Birmingham Vasculitis Activity Score ≥ 3, and patients in remission/low disease activity as determined by Birmingham Vasculitis Activity Score < 3
Full size table

The active disease group had a higher BVAS than the remission/LDA group (median 8 points [IQR 5–10] vs 0 [0–0]; p < 0.0001). Other indicators of disease activity (median [IQR]) were significantly higher in the active group, such as physician VAS (6 [2–8] vs 0 [0–0.5] points; p < 0.0001), CRP (30.9 [6.5–70.7] vs 2.5 [1.9–8.7] mg/L; p < 0.01), and ESR (20.5 [5–38] vs 5 [5–10] mm/h; p < 0.05). However, prednisolone use (7.5 [5–15] vs 5.8 [0.3–13.8] mg/d; p = 0.39) and PR3-ANCA titer (124.6 [4.2–409] vs 22 [2.4–80.8] kIU/L; p = 0.17) did not differ significantly between the groups.

qPCR test distinguishes active disease from remission/LDA in GPA

We hypothesized that the qPCR test could be used as a disease activity parameter because IgG4+ B-cells are implicated in the pathogenesis of GPA. To this end, we compared the qPCR scores between patients in the active disease group (BVAS ≥ 3) and those in the remission/LDA group (BVAS 0–2). The BVAS symptoms can be found in Additional file 1: Table S1.

The median qPCR score was higher in the active disease group than in the remission/LDA group (median 21.4% [IQR 12.1–29.6] vs 3.3% [1.6–5.6]; p < 0.0001) (Fig. 1a). The ROC curve allowed accurate distinction between active disease and remission/LDA at a cutoff value of 11.2%, yielding a specificity of 100% and sensitivity of 87% (AUC 0.993) (Additional file 1: Table S2 and Additional file 2: Figure S1). In this case, specificity was preferred over sensitivity because of the expected clinical relevance (certainty of disease activity). When the LDA samples were excluded from the remission/LDA group, the results remained significant (Additional file 3: Figure S2).

Fig. 1
figure 1

Disease activity parameters in active vs remission/low disease activity (LDA). a Scatter dot plot portraying the percentage of immunoglobulin G subtype 4 (IgG4) from total IgG RNA molecules in the active vs remission/LDA groups in patients with granulomatosis with polyangiitis (GPA). The gray dots within the remission/LDA group represent the three LDA patients. b–e Scatter dot plots showing erythrocyte sedimentation rate (ESR), C-reactive protein (CRP), antineutrophil cytoplasmic antibodies (ANCA), and serum IgG4, respectively, in the active vs remission/LDA groups in GPA (median with IQR). **** p < 0.0001, ** p < 0.01, * p < 0.05, ns = not significant

Full size image

To be of clinical relevance, the qPCR test has to outperform currently used disease activity parameters in GPA, such as ESR, CRP, and ANCA. To this end, the qPCR test was compared with the aforementioned markers using the BVAS as a gold standard for disease activity (Fig. 1b–e). We also tested serum IgG4 as a potential marker in patients when serum was available at the same time point.

The qPCR test outperformed CRP, ESR, and PR3-ANCA level in discriminating active disease from remission/LDA (Fig. 1, Table 2). CRP and ESR levels were also significantly different between the active disease group and the remission/LDA group, but lacked discriminatory capacity (Fig. 1b, c). Serum IgG4 and ANCA titers did not discriminate active disease from remission/LDA (Fig. 1d, e). The IgG4:IgG RNA ratio as determined by the qPCR test did not correlate with serum IgG4 (Additional file 4: Figure S3).

Table 2 Test characteristics of the different available tests detecting disease activity in our granulomatosis with polyangiitis cohort, using Birmingham Vasculitis Activity Score version 3 as gold standard, with Birmingham Vasculitis Activity Score ≥ 3 reflecting active disease
Full size table

Subsequently, the active disease group was further subdivided into a systemic disease group and a limited disease group (for definitions, see the “Materials and methods” section above). Both qPCR score and BVAS did not differ significantly between both groups (Additional file 5: Figure S4). Subdivision of the patients based on the most dominant organ with matching average qPCR score can be found in Additional file 1: Table S3.

Prednisolone use did not differ significantly between both groups (Table 1). More potent agents, such as rituximab (RTX) and cyclophosphamide (CYC) were used by 3 of 15 patients with active disease and by 1 patient in remission/LDA at the time of inclusion. Use of these agents did not result in lower qPCR scores in the active group (data not shown).

In summary, we demonstrated that the qPCR test distinguished active GPA from remission/LDA with very high accuracy using a cutoff at 11.2%. The predictive values of the qPCR test outperformed other frequently used laboratory parameters and the physician VAS score for disease activity.

Elevated qPCR test scores observed in GPA are not observed in other autoantibody-associated diseases

These results prompted us to look in more detail into the specificity of these findings. We wondered whether the high qPCR could also be found in other autoantibody-associated diseases. To this end, we compared active GPA with active SLE and RA (Fig. 2). We added ten healthy control subjects as a negative reference group. We also made a comparison against 15 age-matched patients from the IgG4-RD cohort [19] with active disease as a positive reference. Details on the SLE, RA, and IgG4-RD cohorts can be found in Additional file 1: Table S4.

Fig. 2
figure 2

Granulomatosis with polyangiitis (GPA) vs disease control subjects and healthy control subjects (HC). Scatter dot plot portraying the percentage immunoglobulin G subtype 4 (IgG4) from total IgG RNA molecules in (a) active GPA, systemic lupus erythematosus (SLE), and rheumatoid arthritis (RA) and (b) GPA, SLE, and RA in remission/low disease activity (LDA) and HC. * p < 0.05, ** p < 0.01, *** p < 0.001 (Kruskal-Wallis test with Dunn’s correction)

Full size image

The qPCR score in the active GPA cohort was significantly higher than that in active disease control subjects (median 21.4% [IQR 12.1–29.6] vs 1.3% [0.6–9.7] in SLE [p < 0.001] and 3.9% [2–9] in RA [p < 0.01]) (Fig. 2a). The median qPCR score in the active GPA cohort was similar to that in the active IgG4 RD cohort (Additional file 6: Figure S5).

Interestingly, even when comparing scores of patients in remission/LDA, the qPCR scores were higher in GPA (median 3.8% [IQR 2–5.8] vs SLE 1.1% [0.8–3.1] [p < 0.05] and RA 0.7% [0.5–2.8] [p < 0.05]). The healthy control subjects had the lowest qPCR scores, and their result is shown as a reference value (Fig. 2b).

We observed that, collectively, the qPCR score was significantly higher in active GPA than in other active autoantibody-associated diseases. The same held true when we looked at the qPCR scores in patients with GPA in remission vs patients with control diseases in remission. This suggests that the strong upregulation of cellular IgG4 might be a phenomenon that is specific for GPA. Furthermore, this finding might contain important clues regarding the pathogenesis of GPA.

qPCR test seems to be a promising monitoring tool in GPA

For a disease activity marker to be used in a clinical setting, the result should change when the disease activity changes and should stay stable when the disease activity is stable. To test whether this was true for the qPCR test, we repeated the qPCR test in ten patients at a second time point, creating a longitudinal cohort. Six patients were resampled while in steady state (BVAS change < 3 points). Four patients were in stable remission/LDA, and two (A1 and B) were unsuccessfully treated between the initial and follow-up time points (Fig. 3a). Five patients were sampled after a BVAS decrease of 3 or more points. Of interest, patient A is the same in both panels of Fig. 3. This patient was sampled three times: once after unsuccessful treatment (Fig. 3a) and once after successful treatment (Fig. 3b).

Fig. 3
figure 3

Longitudinal pilot. Longitudinal data in a scatter dot plot portraying the percentage of immunoglobulin G subtype 4 (IgG4) from total IgG RNA molecules in (a) patients who showed no or minor changes in Birmingham Vasculitis Activity Score (BVAS) and (b) patients who showed a decrease of 3 points or more in BVAS after follow-up

Full size image

In the group with a stable BVAS, the qPCR scores did not change (Fig. 3b). We observed a significant decrease in qPCR ratio in the group of patients with a BVAS reduction > 3 points (paired t test, p < 0.05) (Fig. 3b).

qPCR differentiates true vasculitis activity from mimics in six individuals

The unmet need for biomarkers in GPA lies in the ability to distinguish GPA activity from mimics such as infections. To investigate whether the IgG4 qPCR would be able to differentiate between GPA activity and mimics, we performed the IgG4 qPCR test in a real-life clinical situation where clinical symptoms could be caused by either GPA activity or GPA mimics. We tested six patients. Two patients were from the cross-sectional cohort (CS1 and CS2), and four patients were from the longitudinal cohort (E, A2, A1, B) (Table 3). We used additional diagnostic tests and clinical follow-up to determine whether the symptoms were caused by GPA activity or a mimic.

Table 3 Overview of six individuals tested while suspected of having granulomatosis with polyangiitis flare
Full size table

In patients in whom the symptoms were caused by mimics, the qPCR scores were normal (Table 3, columns 2–5). In CS1, CS2, and A2, the CRP levels were elevated due to the ongoing infection. ANCA titer was elevated during infection in two cases (CS2, A2). The qPCR scores were very high in the patients in whom the symptoms were caused by GPA activity (A1, B), indicating that the test is a promising way to distinguish true vasculitis activity from (infectious) mimics. These results need to be confirmed in a larger and more diverse group of patients with GPA.

Correlation between the qPCR test and disease activity as measured by BVAS

The qPCR test accurately distinguished active disease from remission/LDA and decreased when BVAS decreased after treatment. To investigate the correlation in more detail, we analyzed the relationship between the BVAS and qPCR test (Fig. 4a).

Fig. 4
figure 4

Correlation between qPCR test and Birmingham Vasculitis Activity Score (BVAS). a Correlation between BVAS (x-axis) and the percentage of immunoglobulin G subtype 4 (IgG4) from total IgG RNA molecules (y-axis) in all individuals. The blue parts of the graph depict the areas where patients are expected in case of remission/low disease activity (BVAS < 3 and qPCR < 11.2%, lower left quadrant) and activity (BVAS ≥ 3 and qPCR ≥ 11.2%, upper right quadrant). b Correlation between ∆BVAS and ∆qPCR in the ten longitudinal cohort subjects

Full size image

We observed that the BVAS and qPCR showed a significant correlation (p < 0.0001, Spearman r = 0.77). In the ten patients from the longitudinal cohort, we calculated a ∆BVAS and ∆qPCR (Fig. 4b). Overall, these values correlated significantly (p < 0.01, Spearman r = 0.81). Correlation was excellent in the lower disease range. However, the correlation slightly decreased, especially in the patients showing high ∆BVAS and ∆qPCR. Further studies are needed to analyze whether this reflects nonlinearity when measuring vasculitis activity by the qPCR, by the BVAS, or both.

Discussion

In these single-center studies, we demonstrated that the qPCR test associates with disease activity in GPA as measured by BVAS-V3 and outperforms current markers such as PR3-ANCA titer, CRP, ESR, and serum IgG4 level. We confirmed the association in a small longitudinal study. As such, we suggest the qPCR test as a potential new disease activity marker in GPA.

Over the last decade, the BVAS has been instrumental in creating uniformity in scoring GPA disease activity. However, an important limitation is the fact that the scoring system assumes that the physician is certain that the scored symptoms are caused by active vasculitis. When in doubt, physicians turn to available biomarkers such as ESR, CRP, ANCA titer, or sometimes even serum IgG4. In line with previous studies, we found that these markers have limitations in specificity (ANCA, CRP) and sensitivity (ESR) (Table 2) [4, 7,8,9]. To the best of our knowledge, no other more sensitive or more specific disease activity markers have been identified until now [4, 5]. Although our findings appear robust, our study is limited by its size and the limited amount of longitudinal data. The value of the qPCR test should be evaluated in more detail in further longitudinal studies in which there is special attention to disease activity mimics such as infections. Importantly, in four cases in this study, we demonstrated that the qPCR remained low in episodes of mimicked vasculitis activity (Table 3). This further suggests that the qPCR test might be a specific marker that can identify true vasculitis activity.

The qPCR test is based on earlier NGS results in IgG4-RD, where we showed that the qPCR test can accurately discriminate the biliary form of IgG4-RD from its mimics and quantitatively correlates with disease activity in the individual patient [19]. The qPCR test was validated against NGS as a gold standard and showed reproducible results in independent cohorts. We do want to point out that experience with the qPCR test is limited to our center.

The fact that similar results are now found in AAV is not unexpected, given the clear overlap in clinical, histological, and even laboratory features between IgG4-RD and GPA [13,14,15, 17, 23,24,25,26,27,28,29,30,31]. IgG4 might also play a role in other forms of (ANCA-associated) vasculitis. For example, IgG4 serum levels have been correlated with disease activity in eosinophilic GPA [12, 32]. Although there are also many differences in regard to disease presentation between vasculitis and IgG4-RD, the involvement of IgG4+ cells, as measured by the qPCR test in this study, raises the hypothesis that there is a shared (patho-)physiological mechanism involved. This mechanism might by exploited as a potential therapeutic target in the future.

In peripheral blood, the qPCR test quantifies RNA from both IgG4- and IgG-producing cells. The test result is subsequently calculated by dividing the amount of IgG4 RNA by that of total IgG RNA. Thus, an increase in this test result does not reflect an increase in the amount of IgG4 protein present in serum, but rather represents either an increase in the number of IgG4-producing B-lineage cells and/or an increase in the number of IgG4 mRNA molecules expressed per cell, such as during differentiation from mature B cells to plasmablasts. Although we extensively studied the relationship between the qPCR test and the presence of such BCR clones in IgG4-RD previously, we once more confirmed this in one of the patients with GPA (Additional file 7: Figure S6).

It is puzzling that the presence of IgG4+ B cells in the blood correlates with disease activity in GPA. We propose two hypotheses that may explain our findings. The first is that IgG4+ clones are pathogenic in GPA. This may seem unlikely owing to the theoretical inability of IgG4 to initiate pathogenic responses [33,34,35], but in diseases such as pemphigus and myasthenia gravis, IgG4 antibodies have been demonstrated to have pathogenic properties [18, 36,37,38,39,40,41]. The second hypothesis is that the IgG4 response represents a counterregulatory mechanism following vasculitis activity to dampen the inflammation. This would be more compatible with the theory of IgG4 being a suppressive immunomodulator [34]. The lack of correlation between protein IgG4 in serum and disease activity is probably a reflection of continued stable production of IgG4 by plasma cells located in bone marrow or other locations that is unaffected by disease activity elsewhere. We demonstrated that our IgG4/IgG mRNA ratio as measured by the qPCR test, which correlates with disease activity, reflects the real-time differentiation of IgG4-producing B cells to plasmablasts and migration of these cells through the blood to other locations, such as from site of inflammation to bone marrow [18, 19, 42]. It is conceivable that this process begins before the onset of disease activity, thus predicting the onset of a flare. This concept will need further exploration in longitudinal studies.

Subdivision of the active GPA group into limited (five individuals) and systemic (ten individuals) disease groups did not show any differences in qPCR score, suggesting that the IgG4+ B-cell upregulation is not linked to disease manifestation (Additional file 5: Figure S4). However, it has to be noted that these two sample groups were small and that further subdivision based on organ involvement was not possible, because sample groups would be too small to record any significant differences.

Furthermore, one could argue that although not significantly different between the groups, the remission/LDA group had a longer duration since diagnosis and therefore possibly had more cumulative treatment, which could explain these results. However, years since diagnosis and qPCR score did not correlate, indicating that a longer disease duration is unlikely to yield lower qPCR scores (Additional file 8: Figure S7). Instead, we believe that disease activity and shorter duration since diagnosis are linked because it has extensively been demonstrated that disease flares tend to occur within the first few years after GPA diagnosis.

Finally, there were several interesting observations with regard to immunosuppressive drugs in our cohort. Use of prednisone was similar in both the active and remission/LDA groups. Prednisone is a potent agent frequently used in induction of remission. However, only after reaching remission (while prednisone is often continued) did the qPCR scores decrease. Second, most patients were included at first presentation of active disease (de novo, flare), when more potent agents such as RTX and CYC are yet to be initiated. Only four subjects were included at a moment during induction therapy when they had received one or more cycles of RTX or CYC. Our findings do not indicate that the IgG4 qPCR score is directly altered by these agents, implying that the IgG4 qPCR reflects disease activity. Nevertheless, we want to emphasize that this study is too small to make definite conclusions. A larger longitudinal study is needed to study the relationship between drugs and IgG4 qPCR scores. Such a study would need enough responders and nonresponders receiving each drug to differentiate the effect of a drug from a change in disease activity as a cause of the change in IgG4:IgG mRNA ratio. Moreover, the long-term effects of B-cell-depleting therapies such as RTX on the IgG4:IgG mRNA ratio are unknown and would require a study that addresses the relationship between the IgG4 qPCR score and circulating B-cell subsets, including plasmablasts. In IgG4-RD, prednisone was also shown to decrease plasmablast counts, which in turn was correlated with disease activity [18, 19]. Dissecting these relationships might help elucidate the pathophysiology of IgG4+ cells and their responses in GPA.

Conclusions

We report that the IgG4:IgG RNA ratio as measured by a qPCR test might be a new and promising disease activity marker in GPA. The qPCR test is able to accurately distinguish disease activity from remission/LDA and outperforms frequently used disease activity parameters in GPA. Furthermore, it seems feasible to use the qPCR test as a monitoring tool and discriminate disease activity from mimics such as infection. A major limitation of this study is the lack of independent validation of the results. Therefore, the value of the test should be further explored in longitudinal studies. If the present results are confirmed, it will help to optimize treatment, thus preventing morbidity and mortality from under- and overtreatment.

Abbreviations

AAV:

Antineutrophil cytoplasmic antibodies-associated vasculitis

ANCA:

Antineutrophil cytoplasmic antibodies

BCR:

B-cell receptor

BVAS:

Birmingham Vasculitis Activity Score

CRP:

C-reactive protein

CYC:

Cyclophosphamide

ESR:

Erythrocyte sedimentation rate

GPA:

Granulomatosis with polyangiitis

IgG:

Immunoglobulin G

IgG4:

Immunoglobulin G subtype 4

IgG4-RD:

Immunoglobulin G subtype 4-related disease

LDA:

Low disease activity

NGS:

Next-generation sequencing

PR3:

Proteinase 3

RA:

Rheumatoid arthritis

RTX:

Rituximab

SLE:

Systemic lupus erythematosus

VAS:

Visual analogue scale

References

  1. Anwar S, Karim MY. Update on systemic vasculitides. J Clin Pathol. 2017;70(6):476–82 https://doi.org/10.1136/jclinpath-2016-203875.

    Article  CAS  Google Scholar 

  2. Phillip R, Luqmani R. Mortality in systemic vasculitis: a systematic review. Clin Exp Rheumatol. 2008;26(5 Suppl 51):S94–104.

    CAS  PubMed  Google Scholar 

  3. Wallace ZS, Lu N, Unizony S, Stone JH, Choi HK. Improved survival in granulomatosis with polyangiitis: a general population-based study. Semin Arthritis Rheum. 2016;45(4):483–9 https://doi.org/10.1016/j.semarthrit.2015.07.009.

    Article  Google Scholar 

  4. Draibe JB, Fulladosa X, Cruzado JM, Torras J, Salama AD. Current and novel biomarkers in anti-neutrophil cytoplasm-associated vasculitis. Clin Kidney J. 2016;9(4):547–51 https://doi.org/10.1093/ckj/sfw056.

    Article  CAS  Google Scholar 

  5. Monach PA. Biomarkers in vasculitis. Curr Opin Rheumatol. 2014;26(1):24–30 https://doi.org/10.1097/bor.0000000000000009.

    Article  CAS  Google Scholar 

  6. Kyndt X, Reumaux D, Bridoux F, Tribout B, Bataille P, Hachulla E, et al. Serial measurements of antineutrophil cytoplasmic autoantibodies in patients with systemic vasculitis. Am J Med. 1999;106(5):527–33.

    Article  CAS  Google Scholar 

  7. Sanders JS, Huitma MG, Kallenberg CG, Stegeman CA. Prediction of relapses in PR3-ANCA-associated vasculitis by assessing responses of ANCA titres to treatment. Rheumatology (Oxford). 2006;45(6):724–9 https://doi.org/10.1093/rheumatology/kei272.

    Article  CAS  Google Scholar 

  8. Tomasson G, Grayson PC, Mahr AD, Lavalley M, Merkel PA. Value of ANCA measurements during remission to predict a relapse of ANCA-associated vasculitis—a meta-analysis. Rheumatology (Oxford). 2012;51(1):100–9 https://doi.org/10.1093/rheumatology/ker280.

    Article  CAS  Google Scholar 

  9. Thai LH, Charles P, Resche-Rigon M, Desseaux K, Guillevin L. Are anti-proteinase-3 ANCA a useful marker of granulomatosis with polyangiitis (Wegener’s) relapses? Results of a retrospective study on 126 patients. Autoimmun Rev. 2014;13(3):313–8 https://doi.org/10.1016/j.autrev.2013.11.003.

    Article  CAS  Google Scholar 

  10. Luqmani RA, Bacon PA, Moots RJ, Janssen BA, Pall A, Emery P, et al. Birmingham Vasculitis Activity Score (BVAS) in systemic necrotizing vasculitis. QJM. 1994;87(11):671–8.

    CAS  PubMed  Google Scholar 

  11. Stone JH, Hoffman GS, Merkel PA, Min YI, Uhlfelder ML, Hellmann DB, et al. A disease-specific activity index for Wegener’s granulomatosis: modification of the Birmingham Vasculitis Activity Score. International Network for the Study of the Systemic Vasculitides (INSSYS). Arthritis Rheum. 2001;44(4):912–20.

    Article  CAS  Google Scholar 

  12. Vaglio A, Strehl JD, Manger B, Maritati F, Alberici F, Beyer C, et al. IgG4 immune response in Churg-Strauss syndrome. Ann Rheum Dis. 2012;71(3):390–3 https://doi.org/10.1136/ard.2011.155382.

    Article  CAS  Google Scholar 

  13. Chang SY, Keogh KA, Lewis JE, Ryu JH, Cornell LD, Garrity JA, et al. IgG4-positive plasma cells in granulomatosis with polyangiitis (Wegener’s): a clinicopathologic and immunohistochemical study on 43 granulomatosis with polyangiitis and 20 control cases. Hum Pathol. 2013;44(11):2432–7 https://doi.org/10.1016/j.humpath.2013.05.023.

    Article  CAS  Google Scholar 

  14. Carruthers M, Shinagare S, Khosroshahi A, Deshpande V, Stone JH. IgG4 plasma cell infiltration in granulomatosis with polyangiitis (formerly Wegener’s) lung biopsies [abstract 1534]. Presented at the 2012 ACR/ARHP Annual Meeting, November 9–14, 2012, Washington, DC.

  15. Perez Alamino R, Martinez C, Espinoza LR. IgG4-associated vasculitis. Curr Rheumatol Rep. 2013;15(8):348 https://doi.org/10.1007/s11926-013-0348-9.

    Article  Google Scholar 

  16. Alba MA, Milisenda J, Fernandez S, Garcia-Herrera A, Hernandez-Rodriguez J, Grau JM, et al. Small-vessel vasculitis with prominent IgG4 positive plasma cell infiltrates as potential part of the spectrum of IgG4-related disease: a case report. Clin Exp Rheumatol. 2015;33(2 Suppl 89):S-138–41.

    PubMed  Google Scholar 

  17. Chang SY, Keogh K, Lewis JE, Ryu JH, Yi ES. Increased IgG4-positive plasma cells in granulomatosis with polyangiitis: a diagnostic pitfall of IgG4-related disease. Int J Rheumtol. 2012;2012:121702 https://doi.org/10.1155/2012/121702.

    Google Scholar 

  18. Maillette de Buy Wenniger LJ, Doorenspleet ME, Klarenbeek PL, Verheij J, Baas F, Elferink RP, et al. Immunoglobulin G4+ clones identified by next-generation sequencing dominate the B cell receptor repertoire in immunoglobulin G4 associated cholangitis. Hepatology. 2013;57(6):2390–8 https://doi.org/10.1002/hep.26232.

    Article  CAS  Google Scholar 

  19. Doorenspleet ME, Hubers LM, Culver EL, Maillette de Buy Wenniger LJ, Klarenbeek PL, Chapman RW, et al. Immunoglobulin G4+ B-cell receptor clones distinguish immunoglobulin G4-related disease from primary sclerosing cholangitis and biliary/pancreatic malignancies. Hepatology. 2016;64(2):501–7 https://doi.org/10.1002/hep.28568.

    Article  CAS  Google Scholar 

  20. Geetha D, Specks U, Stone JH, Merkel PA, Seo P, Spiera R, et al. Rituximab versus cyclophosphamide for ANCA-associated vasculitis with renal involvement. J Am Soc Nephrol. 2015;26(4):976–85 https://doi.org/10.1681/asn.2014010046.

    Article  CAS  Google Scholar 

  21. Stone JH, Wegener’s Granulomatosis Etanercept Trial Research Group. Limited versus severe Wegener’s granulomatosis: baseline data on patients in the Wegener’s granulomatosis etanercept trial. Arthritis Rheum. 2003;48(8):2299–309 https://doi.org/10.1002/art.11075.

    Article  Google Scholar 

  22. Ramakers C, Ruijter JM, Deprez RH, Moorman AF. Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci Lett. 2003;339(1):62–6.

    Article  CAS  Google Scholar 

  23. Ayuzawa N, Ubara Y, Keiichi S, Masayuki Y, Hasegawa E, Hiramatsu E, et al. Churg-Strauss syndrome with a clinical condition similar to IgG4-related kidney disease: a case report. Intern Med. 2012;51(10):1233–8.

    Article  Google Scholar 

  24. Bando H, Iguchi G, Fukuoka H, Taniguchi M, Kawano S, Saitoh M, et al. A diagnostic pitfall in IgG4-related hypophysitis: infiltration of IgG4-positive cells in the pituitary of granulomatosis with polyangiitis. Pituitary. 2015;18(5):722–30 https://doi.org/10.1007/s11102-015-0650-9.

    Article  CAS  Google Scholar 

  25. Hanioka Y, Yamagami K, Yoshioka K, Nakamura T, Kishida M, Nakamura T, et al. Churg-Strauss syndrome concomitant with chronic symmetrical dacryoadenitis suggesting Mikulicz’s disease. Intern Med. 2012;51(17):2457–61.

    Article  Google Scholar 

  26. Holland M, Hewins P, Goodall M, Adu D, Jefferis R, Savage CO. Anti-neutrophil cytoplasm antibody IgG subclasses in Wegener’s granulomatosis: a possible pathogenic role for the IgG4 subclass. Clin Exp Immunol. 2004;138(1):183–92 https://doi.org/10.1111/j.1365-2249.2004.02566.x.

    Article  CAS  Google Scholar 

  27. Kronbichler A, Gut N, Zwerina J, Neuwirt H, Rudnicki M, Mayer G. Extending the spectrum of a chameleon: IgG4-related disease appearing as interstitial nephritis and mimicking anti-neutrophil cytoplasmic antibody-associated vasculitis. Rheumatology (Oxford). 2015;54(10):1936–8 https://doi.org/10.1093/rheumatology/kev248.

    Article  Google Scholar 

  28. Lomborg N, Marcussen N, Junker P. Retroperitoneal fibrosis as the presenting manifestation of IgG4-related disease in a patient previously diagnosed with granulomatosis with polyangiitis (Wegener’s): one disease or two? Scand J Rheumatol. 2014;43(4):345–7 https://doi.org/10.3109/03009742.2014.905628.

    Article  CAS  Google Scholar 

  29. Touge H, Tomita K, Yamasaki A, Shimizu E. A case of proteinase 3 anti-neutrophil cytoplasmic antibody (PR3-ANCA) positive/IgG4-related lung disease. Respir Med Case Rep. 2017;20:92–4 https://doi.org/10.1016/j.rmcr.2017.01.001.

    PubMed  PubMed Central  Google Scholar 

  30. Della-Torre E, Lanzillotta M, Campochiaro C, Bozzalla E, Bozzolo E, Bandiera A, et al. Antineutrophil cytoplasmic antibody positivity in IgG4-related disease: a case report and review of the literature. Medicine. 2016;95(34):e4633 https://doi.org/10.1097/md.0000000000004633.

    Article  Google Scholar 

  31. Massey J. IgG4-related hypertrophic pachymeningitis coexpressing antineutrophil cytoplasmic antibodies. Neurol Neuroimmunol Neuroinflamm. 2017;4(3):e341 https://doi.org/10.1212/nxi.0000000000000341.

    Article  Google Scholar 

  32. Yamamoto M, Takahashi H, Suzuki C, Tabeya T, Ohara M, Naishiro Y, et al. Analysis of serum IgG subclasses in Churg-Strauss syndrome – the meaning of elevated serum levels of IgG4. Intern Med. 2010;49(14):1365–70.

    Article  Google Scholar 

  33. Aalberse RC, Platts-Mills TA, Rispens T. The developmental history of IgE and IgG4 antibodies in relation to atopy, eosinophilic esophagitis, and the modified TH2 response. Curr Allergy Asthma Rep. 2016;16(6):45 https://doi.org/10.1007/s11882-016-0621-x.

    Article  Google Scholar 

  34. Aalberse RC, Schuurman J. IgG4 breaking the rules. Immunology. 2002;105(1):9–19.

    Article  CAS  Google Scholar 

  35. Rispens T, Ooijevaar-de Heer P, Bende O, Aalberse RC. Mechanism of immunoglobulin G4 Fab-arm exchange. J Am Chem Soc. 2011;133(26):10302–11 https://doi.org/10.1021/ja203638y.

    Article  CAS  Google Scholar 

  36. Huijbers MG, Zhang W, Klooster R, Niks EH, Friese MB, Straasheijm KR, et al. MuSK IgG4 autoantibodies cause myasthenia gravis by inhibiting binding between MuSK and Lrp4. Proc Natl Acad Sci U S A. 2013;110(51):20783–8 https://doi.org/10.1073/pnas.1313944110.

    Article  CAS  Google Scholar 

  37. Klooster R, Plomp JJ, Huijbers MG, Niks EH, Straasheijm KR, Detmers FJ, et al. Muscle-specific kinase myasthenia gravis IgG4 autoantibodies cause severe neuromuscular junction dysfunction in mice. Brain. 2012;135(Pt 4):1081–101 https://doi.org/10.1093/brain/aws025.

    Article  Google Scholar 

  38. Lo AS, Mao X, Mukherjee EM, Ellebrecht CT, Yu X, Posner MR, et al. Pathogenicity and epitope characteristics do not differ in IgG subclass-switched anti-desmoglein 3 IgG1 and IgG4 autoantibodies in pemphigus vulgaris. PLoS One. 2016;11(6):e0156800 https://doi.org/10.1371/journal.pone.0156800.

    Article  Google Scholar 

  39. Qian Y, Jeong JS, Ye J, Dang B, Abdeladhim M, Aoki V, et al. Overlapping IgG4 responses to self- and environmental antigens in endemic pemphigus foliaceus. J Immunol. 2016;196(5):2041–50 https://doi.org/10.4049/jimmunol.1502233.

    Article  CAS  Google Scholar 

  40. Mellbye OJ, Mollnes TE, Steen LS. IgG subclass distribution and complement activation ability of autoantibodies to neutrophil cytoplasmic antigens (ANCA). Clin Immunol Immunopathol. 1994;70(1):32–9.

    Article  CAS  Google Scholar 

  41. Brouwer E, Tervaert JW, Horst G, Huitema MG, van der Giessen M, Limburg PC, et al. Predominance of IgG1 and IgG4 subclasses of anti-neutrophil cytoplasmic autoantibodies (ANCA) in patients with Wegener’s granulomatosis and clinically related disorders. Clin Exp Immunol. 1991;83(3):379–86.

    Article  CAS  Google Scholar 

  42. von Borstel A, Abdulahad WH, Rutgers A, Land J, Stegeman CA, Heeringa P, Sanders JSF. Increased CD38hiCD27+ Plasmablast Frequency in Remission Predicts Relapsing Disease in Granulomatosis with Polyangiitis Patients [abstract]. Arthritis Rheumatol. 2017;69(suppl 10). https://acrabstracts.org/abstract/increased-cd38hicd27-plasmablast-frequency-in-remission-predicts-relapsing-disease-in-granulomatosis-with-polyangiitis-patients/. Accessed 10 Jan 2019.

  43. Klarenbeek PL, de Hair MJ, Doorenspleet ME, van Schaik BD, Esveldt RE, van de Sande MG, et al. Inflamed target tissue provides a specific niche for highly expanded T-cell clones in early hum autoimmune disease. Annals of the rheumatic diseases. 2012;71(6):1088-93. Epub 2012/02/02. https://doi.org/10.1136/annrheumdis-2011-200612. PubMed PMID: 22294635.

  44. Klarenbeek PL, Tak PP, van Schaik BD, Zwinderman AH, Jakobs ME, Zhang Z, et al. Human T-cell memory consists mainly of unexpanded clones. Immunol Lett. 2010;133(1):42-8. https://doi.org/10.1016/j.imlet.2010.06.011. PubMed PMID: 20621124.

Download references

Acknowledgements

Not applicable.

Funding

PLK and NdV were supported by grants from ZonMw, the Netherlands Organisation for Health Research and Development, in the 2Treat program (grant 436001001), as well as by the Veni program (PLK) and the Innovative Medicines Initiative-funded Anti-Biopharmaceutical Immunization: prediction and analysis of clinical relevance to minimize the RISK (ABIRISK) project (www.abirisk.eu) (PLK and NdV).

Availability of data and materials

Patient information, qPCR and NGS data, and the vast majority of original RNA material are available. This data is available from the corresponding author on reasonable request.

Author information

Author notes

  1. P. L. Klarenbeek and N. de Vries contributed equally to this work.

Authors and Affiliations

  1. Department of Rheumatology & Clinical Immunology and Amsterdam Rheumatology and Immunology Center (ARC), Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands

    A. Al-Soudi, M. E. Doorenspleet, R. E. Esveldt, L. T. Burgemeister, A. E. Hak, S. W. Tas, R. F. van Vollenhoven, P. L. Klarenbeek & N. de Vries

  2. Department of Genome Analysis, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands

    A. Al-Soudi & R. E. Esveldt

  3. Department of Experimental Immunology, Amsterdam Infection & Immunity Institute, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands

    A. Al-Soudi, M. E. Doorenspleet, R. E. Esveldt, S. W. Tas, P. L. Klarenbeek & N. de Vries

  4. Department of Vascular Medicine, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands

    B. J. H. van den Born

Authors
  1. A. Al-Soudi
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. M. E. Doorenspleet
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. R. E. Esveldt
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. L. T. Burgemeister
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. A. E. Hak
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. B. J. H. van den Born
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. S. W. Tas
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. R. F. van Vollenhoven
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. P. L. Klarenbeek
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. N. de Vries
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

AAS, PLK, and NdV were responsible for study concept and design. AAS, LTB, AEH, BJHvdB, SWT, RFvV, PLK, and NdV were responsible for patient inclusion. AAS, MED, and REE were responsible for laboratory work. AAS, MED, LTB, REE, AEH, BJHvdB, SWT, RFvV, PLK, and NdV were responsible for analysis and interpretation of data. All authors contributed to drafting the manuscript, critically revised the manuscript, and approved the final version.

Corresponding author

Correspondence to A. Al-Soudi.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the local medical ethics committee of the Amsterdam University Medical Centers, location AMC (MEC_AMC_NL44031.018.13). All study participants gave written informed consent for study participation and publication of any related results. The study was performed according to the principles of the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Supplementary methods. Table S1. BVAS division per active (a) and low disease activity (LDA) subject with the most involved (dominant) organ and qPCR score. The final column shows the number the subject has in Fig. 3. Table S2. Coordinates of the curve: the critical cutoff values for the qPCR test, producing different values for sensitivity and specificity. Table S3. Average qPCR score per dominant organ. Table S4. Patient characteristics for SLE (n = 24) and RA (n = 19) cohorts. IgG4-RD (n = 15) and HC (n = 10) are shown as reference groups. na = not annotated, ns = not significant. (DOCX 25 kb)

Additional file 2:

Figure S1. ROC curve portraying sensitivity and specificity of the qPCR test to distinguish active GPA from remission. The blue arrow points toward the cutoff that yields the highest specificity with the least loss of sensitivity. (JPG 40 kb)

Additional file 3:

Figure S2. Active GPA vs remission GPA (without LDA). Scatter dot plot portraying the percentage of IgG4 from total IgG RNA molecules in the active vs remission groups in GPA without LDA. (BMP 3716 kb)

Additional file 4:

Figure S3. Serum IgG4 vs qPCR test. Scatter dot plot portraying the correlation between serum IgG4 (x-axis) and the percentage of IgG4 from total IgG RNA molecules (y-axis) (Spearman’s r correlation). (BMP 3413 kb)

Additional file 5:

Figure S4. Systemic GPA vs limited GPA. Scatter dot plot with in (a) the qPCR result within the active GPA group divided by systemic and limited disease and in (b) the matching BVAS for this group. (BMP 13494 kb)

Additional file 6:

Figure S5. Active GPA vs active IgG4-RD. Scatter dot plot portraying the percentage of IgG4 from total IgG RNA molecules in the active GPA vs active IgG4-RD control groups. ns = not significant. (BMP 3861 kb)

Additional file 7:

Figure S6. NGS data on IgG4+ BCR clones. NGS data portraying the frequency of IgG4+ BCR clones (red) and IgG clones (black) in a GPA patient (x-axis, left) and a representative healthy control subject (x-axis, right). See an earlier study for a complete description of the amplification protocol [43]. Briefly, a linear amplification is performed using six primers specific for the V segment of the BCR heavy chain. Each primer has a nonbinding common sequence on the 5′ end for the second amplification step. The primers cover all known V-segment alleles. After purification, amplification is performed using the common sequence introduced by the V primers and a common primer on the boundary of the J segment and C segment. Primers are available upon request. Further processing was performed according to the protocol for the MiSeq platform (Illumina, San Diego, CA, USA), which was used for NGS. NGS results were analyzed using a custom pipeline [44]. (PDF 60 kb)

Additional file 8:

Figure S7. Years since diagnosis vs qPCR score. The r and p values demonstrate no correlation between disease duration and qPCR score. (BMP 3328 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Al-Soudi, A., Doorenspleet, M.E., Esveldt, R.E. et al. IgG4:IgG RNA ratio differentiates active disease from remission in granulomatosis with polyangiitis: a new disease activity marker? A cross-sectional and longitudinal study. Arthritis Res Ther 21, 43 (2019). https://doi.org/10.1186/s13075-018-1806-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13075-018-1806-6

Keywords

Arthritis Research & Therapy

ISSN: 1478-6362

Contact us