Introduction

Systemic lupus erythematosus (SLE) is a chronic multisystem autoimmune disease, where treatment is typically long term or even lifelong. With a paucity of advanced therapeutics, most patients are treated with systemic glucocorticoids in addition to immunosuppressants (IS). Despite these, irreversible organ damage and mortality remain unacceptably high. A recent study from the investigators’ center reported that the standardized mortality rate (SMR) of SLE patients in China compared to general population was 3.2 [1]; a meta-analysis of 15 studies comprising 26,101 SLE patients of various ethnicities revealed a similar SMR of 2.7 [2]. Understanding of factors which lead to increased SLE morbidity is therefore paramount.

The principle of treat-to-target (T2T) has been successfully applied to many rheumatological and non-rheumatological diseases. For instance, T2T strategies have dramatically improved the prognosis of patients with rheumatoid arthritis. The concept of T2T in SLE was relatively more recently proposed [3], with several “remission” or “low disease activity” states proposed by an expert task force as treatment targets [4]. The definitions of remission in SLE (DORIS) criteria for remission and lupus low disease activity state (LLDAS) criteria from the Asia–Pacific Lupus Collaboration (APLC) for low disease activity were the most accepted and used definitions [5, 6]. The LLDAS has been both retrospectively and prospectively validated as being associated with reduced damage accrual in several studies [7,8,9], and found also to be more attainable than remission, whilst being no less protective [10]. With only one prospective validation study previously reported [6], the purpose of the present study was to investigate the frequency and determinants of achieving LLDAS, and the influence of LLDAS on disease flare and damage accrual in a single-center prospective longitudinal cohort of Chinese SLE patients.

Materials and methods

Study population

Baseline and follow-up data are prospectively collected from all consecutive patients treated at the Peking University First Hospital in Beijing. Ethics approval was obtained from the ethics committee of the Peking University First Hospital (Project Number: 2017[1284]), and written informed consent obtained from all participants. Subjects fulfilled either the 1997 American College of Rheumatology (ACR) updated classification criteria for SLE [11] or the 2012 Systemic Lupus International Collaborating Clinics (SLICC) classification criteria [12]. No specific treatment algorithm was predefined. Patients were usually followed every 3 months, with 6 months as the maximum allowable interval between consecutive visits. Data from adult SLE patients (≥ 18 years) with at least one follow-up visit between January 2017 and June 2020 were analyzed in this study. Incident patients were defined as the patients who were recruited within 1 year of the onset of their SLE symptoms. Patients were considered lost to follow-up if they had no data recorded for 12 months, and were unable to be contacted after two attempts. The data for these patients collected prior to their loss to follow-up was included in the analysis.

Data collection

The demographics, disease duration at recruitment, follow-up duration, SLE-related manifestations, and organ involvement, as determined by the ACR classification criteria on an “ever present” basis, were collected. SLEDAI-2K [13], PGA (scale 0–3) [14], and disease flare assessed using the SLE flare index (SFI) [15] were collected at each visit. Irreversible disease damage was captured using the Systemic Lupus International Collaborating Clinics Damage Index (SDI) [16] and health-related quality of life was captured using the short form (36) health survey (SF36) [17], both annually. Damage accrual was defined as an increase of ≥ 1 in SDI. All data were recorded in a standardized electronic case report form as part of the APLC longitudinal cohort study [18]. Of note, data in the present study have not been previously included in APLC publications due to the recency of the Peking University Hospital joining this collaboration. Current use and doses of glucocorticoids, hydroxychloroquine (HCQ) and IS, and laboratory results including complete blood count, renal function, serum albumin (ALB), 24-h urine total protein (UTP), complement levels (C3 and C4), and anti-double stranded DNA (dsDNA) antibody titers at baseline and each follow-up visit were also collected.

Definitions

Cutaneous and mucosal involvement, leukopenia, thrombocytopenia, and serositis were defined using the 2012 SLICC criteria [12]. Lupus nephritis (LN) was defined as (i) proteinuria > 0.5 grams per day or > 3+, or (ii) cellular casts that may be red cell, hemoglobin, granular, tubular, or mixed, or (iii) biopsy-proven nephritis compatible with SLE [12, 19]. Neuropsychiatric lupus (NPSLE) included a series of disorders ranging from diffuse central nervous system (CNS) disorders (acute confusional state, psychosis, anxiety and depression, and clinical to subclinical cognitive disorders with variable functional significance) to focal CNS syndromes (seizures, cerebrovascular diseases, chorea, myelopathy, transverse myelitis, demyelinating syndrome, aseptic meningitis, headaches) and peripheral nervous system disorders (polyneuropathies, mononeuropathies, autonomic disorders, plexopathy, acute inflammatory demyelination and polyradiculo-neuropathy) [20, 21]. The diagnosis of autoimmune hemolytic anemia (AIHA) was based on evidence of hemolysis including reticulocytosis, bilirubinemia, increased lactate dehydrogenase (LDH), and a positive direct antiglobulin test [22]. Pulmonary arterial hypertension (PAH) was defined as a mean pulmonary artery pressure ≥ 25 mmHg and a pulmonary arterial wedge pressure ≤ 15 mmHg by right heart catheterization, or a mean estimated pulmonary systolic pressure ≥50 mmHg on echocardiography [23].

LLDAS was defined as an SLE disease activity index (SLEDAI)-2 K of ≤ 4, no activity in any major organ system and no features of new disease activity, a physician global assessment (PGA, 0–3) ≤ 1, prednisone dose ≤ 7.5 mg/day, and allowance for maintenance IS and anti-malarials [6].

Data analysis

Data are reported as mean (standard deviation (SD)) for normally distributed continuous variables, median (interquartile range (IQR)) for skewed continuous data, and percentages or proportions for categorical variables. The Student’s t-test was used for comparisons of normally distributed continuous variables, Kruskal–Wallis and Mann–Whitney U tests for comparisons of non-normally distributed continuous variables, and the chi-squared test for comparisons of categorical data. Univariate and multivariable logistic regression was used to identify factors associated with achieving LLDAS or damage accrual. Variables with a p value ≤ 0.05 in simple logistic regression analyses were included in stepwise multivariable regression analysis. Variables which relate directly to LLDAS assessment such as SLEDAI, PGA, and prednisone dose were not included in the logistic regression model for achieving LLDAS. Flare rate analysis during the follow-up was performed using the Kaplan–Meier method with comparisons performed using the log-rank test. The primary endpoint was disease flare or data censoring. The duration of follow-up was defined as the time from recruitment until the first flare or last follow-up.

All analyses were performed with STATA version 13.1 (StataCorp, College Station, TX, USA) for Windows and a p-value of < 0.05 was considered statistically significant.

Results

Subject characteristics

Between January 2017 and June 2020, 200 SLE patients were recruited in the prospective cohort. After excluding 15 patients who did not return for follow-up after the first visit, 185 patients including 54 incident patients (29.2%) with total 1203 visits were included in this study. The mean age at disease onset was 33.5 ± 14.9 years with a female predominance (88.1%). The median (IQR) duration at recruitment was 2.3 (0.8–7.7) years, median follow-up duration was 2.2 (1.0–2.9) years, and median number of visits was 7 (4–9). Among the 185 patients, 19 patients were lost to follow-up.

In terms of clinical characteristics, cutaneous and mucosal involvement was most common (61.6% of the patients), followed by LN (55.7%), leukopenia (44.3%), arthritis (39.5%), thrombocytopenia (25.4%), serositis (15.3%), AIHA (9.7%), and NPSLE (8.7%). Other less common disease manifestations included myositis (6.5%), gastrointestinal tract involvement (6.0%), and PAH (4.3%). Renal biopsy was conducted in 62 (33.5%) patients. The median (IQR) SLEDAI-2 K, PGA, and SDI at recruitment were 2 (2–6), 1 (1–2), and 0 (0–1), respectively. One hundred and eighty-two patients (98.4%) received glucocorticoid treatment, 171 (92.4%) received HCQ, and 170 (91.9%) received IS for at least 3 months. During follow-up, 111 (82.1%) patients received IS for at least 6 months and 109 (73.2%) for at least 12 months.

Frequency of LLDAS attainment

In the 185 patients, 58 (31.4%) patients fulfilled LLDAS at recruitment, 81 (43.8%) patients achieved LLDAS during follow-up, and 45 (24.3%) never achieved LLDAS. Therefore, there were 139 (75.1%) patients who achieved LLDAS at least once, including 29 (15.7%) patients who were in LLDAS for 100% of observations, 53 (28.7%) for 50 to 100% of observations, and 58 (31.4%) in LLDAS for < 50% of observations (Fig. 1A). In a subgroup of 144 patients who completed at least 12 months of follow-up, the attainment of LLDAS during the follow-up were shown in Fig. 1B.

Fig. 1
figure 1

The distribution of percentage observations in LLDAS during follow-up. A All patients; B patients with at least 12 months of follow-up

Characteristics of patients achieving LLDAS attainment for ≥ 50% of observations

The characteristics of patients who achieved LLDAS for ≥ 50% of observations were compared to those who did not (Table 1). Compared with subjects who achieved LLDAS for at least 50% of observations (referred to as the LLDAS ≥ 50% group), those who achieved LLDAS for less than 50% of observations (referred to as LLDAS < 50% group) had a higher UTP (390 [50–2300] vs 50 [20–150] mg, p < 0.001), lower albumin (39.2 ± 5.3 vs 42.4 ± 5.2 g/L, p < 0.001), higher Scr (86.7 ± 34.2 vs 75.6 ± 17.2, p = 0.027), lower C3 (714 ± 209 vs 857 ± 200 mg/L, p < 0.001), higher SLEDAI-2 K (4 [2,3,4,5,6] vs 2 [0–2], p < 0.001), higher PGA (1 [1, 2] vs 0 [0–1], p < 0.001), and higher prednisone dose (11.9 ± 3.9 vs 5.4 ± 2.6 mg, p < 0.001) at recruitment. The LLDAS < 50% group had a significantly higher proportion of patients who experienced a disease flare during the follow up period (58.6% vs 27.5%, p < 0.001) and were more likely to have damage accrual after 2 years of follow-up (48.8% vs 21.4%, p = 0.004) than the LLDAS ≥ 50% group. There were no significant differences in individual organ involvement or auto-antibodies between the two groups. Gender, education level, disease duration and SDI at recruitment, and frequency of HCQ and IS use at recruitment and during follow-up were also comparable between the two groups (Table 1).

Table 1 Characteristics of SLE patients, categorized by proportion of follow-up observations with achievement of LLDAS

Determinants of achieving LLDAS for ≥ 50% of observations

Univariate logistic regression analyses showed that higher C3 (per 100 mg/L) (OR = 1.43, 95%CI [1.13–1.74], p = 0.000) and Alb (per g/L) (OR = 1.14, 95% CI [1.05–1.24], p = 0.001) level at recruitment were positively associated with achieving LLDAS for ≥ 50% of observations. Similarly, higher UTP (per g) (OR = 0.426, 95%CI [0.225–0.805], p = 0.009) and Scr (per 10 μmol/L) (OR = 0.79, 95%CI [0.64–0.97], p = 0.026) level at recruitment were negatively associated with achieving LLDAS for ≥ 50% of observations. In the multivariable logistic regression model, UTP (per g) (OR = 0.447, 95%CI [0.207–0.968], p = 0.041), Scr (per 10 μmol/L) (OR = 0.72, 95%CI [0.52–0.99], p = 0.040), and C3 (mg/L) (OR = 1.60, 95%CI [1.18–2.17], p = 0.003) at recruitment remained significantly associated with achieving LLDAS for ≥ 50% of observations (Table 2).

Table 2 Determinants of achieving LLDAS in at least 50% of observations by univariate and multivariable logistic regression analysis

Protective effect of LLDAS against flare during follow-up

By the end of the study, 73 (39.4%) patients had experienced a disease flare, a mild to moderate flare in 53 (28.6%) patients and severe flare in 20 (10.8%) patients. Fifty-four (23.8%) patients flared within the first 12 months of follow-up and 19 (10.3%) patients after 12 months.

Kaplan–Meier analyses showed that the overall cumulative flare rate was 12% at 6 months, 23% at 12 months, 45% at 24 months, and 52% at 36 months. The flare rate showed a significant difference among patients who never achieved LLDAS, LLDAS < 50%, and LLDAS ≥ 50%, with the lowest flare rate in the LLDAS ≥ 50% group (Log-rank p = 0.000, Fig. 2). Meanwhile, the attainment of LLDAS achievement during the first follow-up year was associated with a reduced the risk of flare after the 1st year of follow-up (Log-rank p = 0.049, Fig. 3).

Fig. 2
figure 2

Occurrence of flare curves during follow-up (Kaplan–Meier analysis). The comparison of flare rate was among patients who achieved LLDAS for more than 50% of observations, those who achieved LLDAS less than 50% of observations, and those who never achieved LLDAS (log rank p = 0.000)

Fig. 3
figure 3

Flare in the second year of follow-up (Kaplan–Meier analysis) in patients categorized by time spent in LLDAS (LLDAS more than 50% of observations, LLDAS less than 50% of observations, and never achieved LLDAS) during the first year of follow-up (log rank p = 0.049)

Protective effect of LLDAS on damage accrual during follow-up

The SDI distribution at recruitment and during follow-up is shown in Supplementary Table 1. During follow-up, damage accrual was observed in 21/141 (14.9%) patients within 1 year, 37/104 (35.6%) patients within 2 years, and 11/21 (52.4%) patients within 3 years. In a multivariable logistic regression model, age at disease onset and duration of follow-up (per year) were independent risk factors for damage accrual (OR = 1.05, 95%CI [1.01–1.11], p = 0.047; and OR = 6.94, 95%CI [1.89–25.44], p = 0.003; respectively), and achieving LLDAS at least 50% of observations during follow-up had a significant protective effect on damage accrual (OR = 0.19, 95%CI [0.04–0.99], p = 0.049) (Supplementary Table 2).

Discussion

LLDAS, an outcome measure developed by the multinational Asia–Pacific Lupus Collaboration, has been validated retrospectively and prospectively to be associated with a reduced risk of organ damage in studies from Latin America, North America, Europe, and the Asia–Pacific region [6, 9, 24]. Here, we present a prospective single-center longitudinal study of LLDAS attainment and its association with disease flare and damage accrual in Chinese patients.

In our cohort, 75.1% of the patients achieved LLDAS at least once and 44.3% were in LLDAS at least 50% of visits, which were comparable with other reports [7, 25] and confirmed that LLDAS is an achievable goal, even in newly diagnosed patients who are more likely to have active disease [26].

Our analyses showed that there were significant differences in baseline characteristics of patients who achieved LLDAS for at least 50% of follow-up observations and those who did not, in terms of baseline UTP, Scr, Alb, and C3. As expected, the variables which are directly involved in the LLDAS assessment criteria such as SLEDAI-2 K, PGA, and prednisone dose were also significantly different between the two groups. Compared with other serological markers such as C4 and anti-dsDNA, baseline C3 was shown in logistic regression analyses to be the most significant marker for subsequent LLDAS attainment. Two previous studies showed that renal involvement was negatively associated with achieving LLDAS [7, 27], but the present study identified that two more specific indicators (UTP and Scr) of renal involvement were independent negative predictive factors for achieving LLDAS ≥ 50% of observations. Of course, we need to clarify that the patients with LN accounted for more than 50% of all patients in this cohort. In case of without these LN patients, the predicting effect of UTP and Scr would not be statistically significant. A series of epidemiological studies had consistent results with us that nephropathy is one of the most common major organ involvements in SLE patients in China and Asian-Pacific region [1, 7, 28]. The present study further emphasizes the critical role of renal involvement in affecting therapeutic target achievement.

The present study revealed that percentage of time in LLDAS was negatively associated with disease flare. Disease flare is often managed with the use of glucocorticoids and IS, which in turn may lead to further damage accrual or infection risk. A study from the Netherlands showed that having at least one major flare was associated with future damage accrual [25]. The recommendations from the SLE-T2T taskforce underlined that prevention of flares should be a therapeutic goal in SLE and the present study identifies LLDAS as an achievable and important target for preventing flare. These data are encouraging for the future application of the LLDAS in clinical practice as a treatment target.

In the present cohort, we found 14.9% of patients had damage accrual within 1 year, and we observed a significant protective association of LLDAS with damage accrual. Achieving LLDAS ≥ 50% of observations was an independent protective factor for damage accrual. Early damage accrual is a particularly important outcome indicator for SLE patients as it is predictive of further damage accrual [29] and lower survival. For example, one study showed that patients with initial damage had a fourfold higher mortality rate compared to those with no early damage [30]. Another study showed that damage 1 year after diagnosis was a significant predictor of death within 10 years of diagnosis [31].

Interestingly, age of disease onset was found to have different associations with disease flare and damage accrual, wherein older age at disease onset was an independent protective factor for disease flare, while it was an independent risk factor for damage accrual. One study found that patients with younger age at disease onset were more likely to have active disease than those with older age at disease onset [32]. This implies that more vigilant care is needed for these patients. However, consistent with another report [8], older patients were more likely to develop permanent organ damage which need closer monitor and prophylactic treatment.

There are some limitations to this study. First, our sample size was relatively small, the follow-up duration was relatively short especially in those who did not achieve LLDAS, and the loss to follow-up ratio was relatively high, which limits the evaluation for organ damage development and the correlation between LLDAS achievement and damage accrual. Secondly, our assumption that patients are in LLDAS for the entire duration between two visits at which they are found to be in LLDAS may overestimate the duration of being in LLDAS.

In conclusion, LLDAS was an attainable goal in clinical practice in Chinese SLE patients and therefore a potentially applicable goal in a T2T strategy. Higher urine protein and serum creatinine level and lower C3 level at baseline were negative influence factors for achieving LLDAS subsequently. LLDAS is a promising treatment target in SLE, negatively associated with disease flare and damage. Future prospective studies, using large sample sizes, are expected to validate these findings.