{'fraction_non_alphanumeric': 0.041271242551313175, 'fraction_numerical': 0.026925623482674906, 'mean_word_length': 5.0639009702241555, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 2, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Prediction of disease flare by biomarkers after discontinuing biologics in patients with rheumatoid arthritis achieving stringent remission

To elucidate the disease-flare process in rheumatoid arthritis (RA) after discontinuing biological disease-modifying antirheumatic drugs (bDMARDs), we first focused on RA-flare prediction after achieving stringent remission criteria. Patients with RA who maintained a simplified disease activity index ≤ 3.3 for ≥ 3 months during November 2014–January 2018 in our medical centre in Tokyo, Japan, were eligible. The primary endpoint was flare (disease activity score 28—erythrocyte sedimentation rate ≥ 3.2 with increase from baseline > 0.6) within 2 years after bDMARD discontinuation. Comprehensive clinical assessments, ultrasonographic evaluation of 40 joints, and blood sampling for 12 biomarkers were performed every 2–3 months for 2 years unless patients experienced flare. Flare-positive and flare-negative patients were compared using univariate and Kaplan–Meier analyses. Thirty-six patients (80.6% female, median disease duration, 5.2 years; median treatment period with discontinued bDMARD, 2 years; median remission duration, 18 months) were enrolled. Twenty patients (55.6%) experienced RA flare 43–651 (median, 115) days after the first skipped date of bDMARDs. Two patients who withdrew without disease flare were excluded from the comparison. Clinical and ultrasonographic evaluations did not show significant between-group differences; Kaplan–Meier analysis showed that higher baseline soluble tumour necrosis factor receptor 1 (sTNFR1) concentration impacted subsequent disease flare (p = 0.0041); higher baseline interleukin (IL)-2 concentration was exclusively beneficial to patients with lower sTNFR1 (p = 0.0058), resulting in remission maintenance in 83.3% of patients with lower sTNFR1 and higher IL-2. We demonstrated the usefulness of combined biomarker evaluation for predicting sustained remission after bDMARD discontinuation in RA.

activity by re-treatment with the same bDMARD when they experienced disease flare after bDMARD discontinuation in the above studies.
Therefore, this prospective study was conceived to elucidate the process of RA flare by sequential assessment of clinical, imaging, and biomarker parameters. In this first report, we focused on the prediction of RA flare after bDMARD discontinuation by assessing the above-mentioned parameters starting at the point of achieving ACR/EULAR-defined clinical remission.

Results
Demographic and clinical characteristics of enrolled patients. Demographic and clinical characteristics of the 36 patients are summarised in Table 1. Two patients did not complete the 2-year study, one due to sudden relocation on day 43 and the other because of death due to comorbidity (cardiomyopathy) on day 175; neither had experienced disease flare at the time of drop out. Thus, these patients were excluded from the comparative study between flare-positive and flare-negative patients. Twenty-nine (80.6%) patients were female, and the median (interquartile range [IQR]) age at enrolment was 69 (53-75) years. The median (IQR) disease duration was 5.2 years, treatment with the discontinued bDMARD had lasted for 2.1 (1.4-3.5) years, and the median (IQR) remission duration was 20 (9-31) months (< 6 months in only two patients). Most (72.2%) patients had been under treatment with first bDMARD and 76.5% and 77.8% of the patients were respectively positive for anti-cyclic citrullinated peptide antibody (anti-CCP) and rheumatoid factor.
In line with the fact that all enrolled patients fulfilled the SDAI ≤ 3.3 criterion, the medians (IQRs) of TJC28 and SJC were both 0 (0-0) as were those of Health Assessment Questionnaire-disability index (HAQ-DI). The median PtGA (mm VAS) was 3, and the median PhGA (mm VAS) was 0. The median CRP value was 0.04 mg/ dL, and the median disease activity score 28-erythrocyte sedimentation rate (DAS28-ESR) 27 and SDAI were 2.0 and 0.4, respectively. Importantly, 31 (86.1%) of 36 patients had achieved the ACR/EULAR Boolean remission criteria, which have been shown to be the most stringent criteria of clinical RA remission 4 . The ultrasonographic (US) examination of the 40 joints showed that the total power Doppler (PD) score was '0' in 26 patients and '1' in 5 patients. Therefore, imaging remission was confirmed in at least 72.2% of the enrolled patients.
Twenty-six (72.2%) patients had discontinued anti-tumour necrosis factor (TNF) bDMARD, while four and six patients had discontinued the non-TNF abatacept and tocilizumab, respectively. Methotrexate (MTX) was the predominant (69.4%) concomitant csDMARD as an anchor drug for RA in Japan and other countries 1,28 , while concomitant glucocorticoids were rarely administrated to patients (8.3%) who achieved SDAI remission.
Clinical comparison between patients experiencing RA flare and those maintaining remission. Among the 36 enrolled patients, 20 patients (55.6%) experienced RA flare within 2 years, with a range of 43-651 (median 115) days after the first skipped date of bDMARD infusion/injection (Fig. 1A). Notably, 18 (90%) of those patients received the same bDMARD thereafter and again achieved SDAI remission within 2 months. Another bDMARD was successfully introduced for one patient, and the remaining patient achieved SDAI remission after 9 months without treatment adjustment. Thus, the outcome of patients after experiencing disease flare was very favourable. In addition, two patients had not received any csDMARDs at baseline. One patient flared after tocilizumab discontinuation, and another patient did not flare for 2 years after etanercept discontinuation.
The comparison between flare-positive patients (n = 20) and patients who remained in remission (n = 14) is also shown in Table 1. The demographic characteristics were comparable between the two groups. Among the activity-related parameters, the ESR value and DAS28-ESR were numerically greater in flare-positive patients than in flare-negative patients, although the differences were not statistically significant (p = 0.077 and p = 0.058, respectively). US findings were also comparable for the greyscale (GS), PD, and total scores between flare-positive and flare-negative patients (p = 0.54, p = 0.88, and p = 0.53. respectively). Further, the proportion of patients who had a total PD score ≥ 2 was similar between the flare-positive and flare-negative groups (15.0% and 14.3%, respectively; p = 1.00).
As for the RA treatments, disease flare was similarly observed among patients who had discontinued TNF inhibitors versus non-TNF inhibitors, and those who had discontinued the first bDMARD versus the second or third bDMARD (Fig. 1B,C, Table 1). Patients who received MTX showed a trend for maintaining remission compared with those who did not receive MTX ( Fig. 1D and Table 1), although the differences were not statistically significant. Therefore, clinical and imaging prediction of disease flare after bDMARD discontinuation in the patients who achieved the stringent criteria of RA remission was not suggested by this study.
Biomarker comparison between patients experiencing RA flare and those maintaining remission. Next, we compared the baseline plasma concentrations of nine cytokines and three soluble cytokine receptors relevant to the pathogenesis of RA and/or molecular targets of RA treatment (Table 2). Interestingly, although plasma IL-2 concentration was generally very low and undetectable in eight of 36 patients, it was significantly greater in flare-negative patients than in flare-positive patients (p = 0.017). A univariate logistic regression analysis and subsequent receiver operating characteristics (ROC) curve analysis suggested a cut-off value of 0.06287 pg/mL with an area under the curve (AUC) value of 0.74 ( Supplementary Fig. 1A).
Moreover, the Kaplan-Meier analysis showed that lower baseline IL-2 concentration was associated with subsequent disease flare (p = 0.020, Supplementary Fig. 1B).
Although plasma concentrations of IL-6, TNF, and soluble tumour necrosis factor receptor 2 (sTNFR2) may be affected by some bDMARDs, a marginal difference in soluble tumour necrosis factor receptor 1 (sTNFR1) concentration between flare-positive and flare-negative patients (p = 0.064) led us to further analyses. The plasma sTNFR1 concentration was not correlated with CRP (r 2 = 0.0044, p = 0.70; Fig. 2A), but it was significantly     www.nature.com/scientificreports/ of that of IL-2 (p = 0.73), and a higher baseline IL-2 concentration was beneficial for patients with a sTNFR1 concentration below the cut-off value (p = 0.0058, Fig. 3B), but not for those with sTNFR1 concentration above the cut-off value (p = 0.81, Fig. 3C).

Discussion
This is the first study to comprehensively evaluate patients achieving clinical remission defined by the ACR/ EULAR who had discontinued bDMARDs from the clinical, imaging, and biomarker perspectives. Our results indicated the usefulness of combining biomarker, but not clinical and imaging, assessments for the prediction of disease flare referring to the Outcome Measures in Rheumatology (OMERACT) validation study 29 after bDMARD discontinuation in patients fulfilling the stringent criteria of clinical remission. In this study, the median baseline values of enrolled patients were '0' for TJC28, SJC28, PhGA (mm VAS), HAQ-DI, and total PD score of 40 joints. Despite such sufficient suppression of disease activity, 20 of 36 patients (55.6%) experienced disease flare within the 2 years of a bDMARDs-free period. Interestingly, baseline US examination and the confirmation of imaging remission added little to the prediction of subsequent RA flare in this study, which enrolled patients who had achieved the stringent criteria for RA remission. This fact is concordant with a recent idea that does not support the routine use of magnetic resonance imaging or US joint examination in the treat-to-target strategy 30 . However, two patients with US PD scores of 22 and 10, experienced disease flare at 55 and 105 days, respectively, after bDMARD discontinuation. These patients had one swollen joint in 28 joints and elevated MMP-3 (229.4 and 176.2 ng/mL, respectively). Therefore, US examination may be useful for patients achieving the stringent remission criteria despite some signs of residual arthritis activity. Nonetheless, it should be noted that imaging remission by US does not guarantee sustained remission after bDMARD discontinuation.
Several previous studies have suggested that longer disease duration, anti-CCP positivity, baseline radiographic damage, and residual disease activity at baseline were associated with disease flare after bDMARD discontinuation 1,31,32 . In fact, disease duration, anti-CCP positivity, and DAS28-ESR were numerically greater in flare-positive patients than in flare-negative patients, although statistically not significant in this study with a small sample size (Table 1)   www.nature.com/scientificreports/ Remission study 19 found that multi-biomarker disease activity (MBDA) at baseline was significantly higher in patients with RA who relapsed than in those remaining in stable remission and those tapering/stopping DMARDs 33 . An MBDA score is calculated by the expression of 12 biomarkers including sTNFR1, IL-6, vascular endothelial growth factor, CRP, and MMP-3 which were also evaluated in this study 34 . The sTNFR1 is generated as a result of proteolytic shedding of TNFR1 by tumour necrosis factor-α-converting enzyme (TACE), a member of the ADAM (a disintegrin and metalloproteinase) family of zinc metalloproteinases 35 . As shown in Fig. 2, sTNFR1 concentration was positively correlated with DAS28-ESR, and the DAS28-ESR value was significantly higher in patients with sTNFR1 above the cut-off value than in those with sTNFR1 below the cut-off value (medians, 2.3 and 1.6, respectively, p = 0.022). These results suggest that sTNFR1 which is incorporated in the MBDA score, is a valid biomarker of RA disease activity even in patients achieving SDAI ≤ 3.3. In addition to the results shown in Fig. 3, a sensitivity analysis in patients discontinuing TNF inhibitors (n = 26) revealed that RA flare was observed in 83.3% (10 of 12) of patients with sTNFR1 above the cut-off value which was significantly greater than the flare rate of 28.6% (4 of 14) patients with sTNFR1 below the cut-off value (p = 0.0079). Furthermore, combining the sTNFR1 and IL-2 biomarkers is promising for the prediction of RA flare after bDMARD discontinuation (Fig. 2).
Although the blood concentration of IL-2 was also found to be very low (approximately 1 pg/mL) in several other studies of non-rheumatic diseases [36][37][38] and mostly below the lower limit of detection in a very recent study of predicting drug-free RA remission by a composite score incorporating three transcripts, IL-27, and Boolean remission 39 , and thus only detectable using high-sensitivity measurements, IL-2 is crucial for the development and expansion of regulatory T cells 40,41 which may be associated with sustained RA remission. The limitations of this study include the small sample size from a single centre and the lack of histological or cellular analyses. However, the comprehensive and intensive evaluations from multiple aspects including clinical, imaging and biomarker analyses are the strengths of the study.
In conclusion, this was the first study to identify the usefulness of a combined biomarker evaluation of sTNFR1 and IL-2 for the prediction of sustained remission for 2 years after bDMARD discontinuation in patients with RA. Recently, the difference between clinical remission and 'molecular remission' has been proposed by multi-omics monitoring of responses to bDMARDs in patients with RA 42 . Further analyses focusing on the next objectives including earlier detection of RA flare by comprehensive evaluations and the elucidation of relationships among clinical, imaging, and biomarker changes are ongoing in our laboratory.  Study design. This was a 2-year, prospective, observational study. The primary endpoint was flare DAS28-ESR ≥ 3.2 with an increase from the baseline, > 0.6 29 ) within 2 years after bDMARD discontinuation. The study protocol was approved by the Ethics Committee of Toho University Ohashi Medical Center (Ohashi approval number, , and this study was carried out in compliance with the Declaration of Helsinki. Written informed consent was obtained from all participants. The target sample size was set to 40 patients to obtain data of 20 patients experiencing disease flare for intensive analyses of clinical, imaging, and biomarker changes estimating a 50% flare rate. Patient recruitment was started in November 2014 and closed in January 2018 with the enrolment of 36 patients.
Clinical assessments. We obtained the following medical information from medical records: age, sex, height, weight, disease duration from the onset of symptoms to enrolment, duration of bDMARD treatment to be discontinued, remission duration, Steinbrocker radiographic stage, anti-CCP, rheumatoid factor, MMP-3, and the treatment for RA. After enrolment in this study, the patients were examined every 2-3 months for TJC28, SJC28, PtGA, PhGA, HAQ-DI, blood tests for ESR, serum CRP, and MMP-3 as well as plasma storage for future examination of cytokine and soluble cytokine receptor concentrations, and also underwent US examination of the 40 joints.

US examination.
A Xario (Toshiba Medical Systems, Tochigi, Japan) US machine equipped with a multifrequency linear array probe (7-14 MHz) was used. The PD settings were as follows: PD pulse repetition frequency, 16.5 kHz (flow range, 3.8 cm/s); Doppler frequency, 6.1 MHz, and low wall filter. Colour gain was set just below the level at which noise appeared. The US examination was performed according to the EULAR guidelines for musculoskeletal US in rheumatology 44,45 by rheumatologists and ultrasonographers certified by the Japan College of Rheumatology. These examiners evaluated 40 of the joints including the bilateral shoulders, elbows wrists, metacarpophalangeal joints, thumb interphalangeal joints, second through fifth finger proximal interphalangeal joints, knees, ankles, and metatarsophalangeal joints using longitudinal and transverse scans, without knowledge of the other medical information of the patients.
The US examiners performed the final scoring of the recorded US findings according to the OMERACT definitions 46 . Briefly, the GS was graded semiquantitatively from 0 to 3 (0 = absent, 1 = mild, 2 = moderate, and 3 = marked) in a combined measure of synovial hypertrophy and fluid retention of the articular recess. The intraarticular PD signals were also graded on a 0-3 scale, and the total joint US score was determined by summing all GS and PD scores for each patient.

Measurement of plasma concentration of cytokines and soluble cytokine receptors.
Plasma was stored at − 80 °C until measurement. Plasma concentrations of interleukin (IL)-1β, IL-2, IL-6, IL-8, IL-10, interferon gamma, TNF, granulocyte-macrophage colony stimulating factor, vascular endothelial growth factor, sTNFR1, sTNFR2, and soluble IL-6 receptor were measured by electrochemiluminescence assay with the Ultra-Sensitive Kit (Meso Scale Discovery, Rockville, MD, USA) according to the manufacturer's protocol 47 . Statistical analyses. The statistical analyses were performed using JMP Pro (version 14.2, SAS Institute Japan Ltd., Tokyo, Japan). Continuous variables are presented as the median and interquartile range (IQR) and were analysed using the Mann-Whitney U test. Binomial data were compared between two groups using Fisher's exact test. Kaplan-Meier analysis and log-rank tests were used for the comparison of survival curves between the groups. The correlation coefficient (r) was determined by Spearman's rank test. Logistic regression analyses followed by ROC curve analyses were conducted for the prediction of RA flare based on the baseline biomarkers. p-values < 0.05 were considered statistically significant.