Abstract
Objective. To characterize the effect of anti-tumor necrosis factor (TNF) therapy compared to conventional synthetic disease-modifying antirheumatic drugs (csDMARD) in children with enthesitis-related arthritis (ERA) over the first year after diagnosis.
Methods. We conducted a multicenter retrospective comparative effectiveness study of children diagnosed with ERA. We estimated the effect of anti-TNF therapy on clinical variables (active joint count, tender entheses count) and patient-reported pain and global assessment of disease activity over the first year after diagnosis using state-of-the-art comparative effectiveness analytic methods.
Results. During the study period, 217 patients newly diagnosed with ERA had a total of 965 clinic visits the first year after disease diagnosis. Children [median age 11.6 yrs, interquartile range 10–14] were treated with anti-TNF monotherapy (n = 33, 15.2%), csDMARD monotherapy (n = 73, 33.6%), or both (n = 52, 23.9%) in the first year after disease diagnosis. There was a statistically significant improvement in the primary outcome, active joint count, over time in children who received an anti-TNF drug versus those who did not (p = 0.03). Additionally, use of anti-TNF therapy versus no anti-TNF therapy was associated with less patient-reported pain (p < 0.01) and improved disease activity over time as assessed by the clinical Juvenile Arthritis Disease Activity Score (p < 0.01). The magnitude of estimated effect on clinical outcomes was uniformly greater, with the exception of tender entheses count, in children treated with an anti-TNF drug versus a csDMARD.
Conclusion. During the first year after diagnosis, anti-TNF exposure was associated with benefits for several clinically meaningful outcomes in children with enthesitis-related arthritis.
- ANTI–TUMOR NECROSIS FACTOR AGENTS
- ENTHESITIS-RELATED ARTHRITIS
- CONVENTIONAL SYNTHETIC DISEASE-MODIFYING ANTIRHEUMATIC DRUGS
The comparative effectiveness of different treatment algorithms for children with enthesitis-related arthritis (ERA) remains unclear and without consensus. ERA is a category of arthritis characterized by arthritis, enthesitis, axial arthritis, symptomatic uveitis, and HLA-B27 positivity. Treatment regimens for ERA include monotherapy or combination therapy with any of the following: nonsteroidal antiinflammatory drugs (NSAID); conventional synthetic disease-modifying antirheumatic drugs (csDMARD) such as methotrexate (MTX), sulfasalazine (SSZ), and leflunomide (LEF); or biological DMARD (bDMARD) agents such as etanercept (ETN), adalimumab (ADA), and infliximab (IFX). There are only 2 randomized clinical trials focused on children with ERA1,2; the majority of trials include ERA as one of several juvenile arthritis categories. Both of the published ERA trials included only children who had established disease and failed at least 1 NSAID and 1 csDMARD and had at least 3 active joints1,2. In both trials, treatment with a bDMARD resulted in sustained clinical improvement. In another study that included a subset of children with prevalent ERA disease (average disease duration 2 yrs) and at least 2 active joints, ETN resulted in improvement in the pediatric American College of Rheumatology (ACR) core set response criteria, tender entheses count, back pain, and back mobility3. There are no published trials of therapy for children with a new diagnosis of ERA. The choice of induction treatment algorithms for children with ERA, according to the ACR treatment recommendations, is based solely on the number of active joints4. Using these algorithms, the earliest a child with ERA might be treated with a bDMARD is after 3 months of therapy with a csDMARD. Sacroiliitis is considered separately in the ACR recommendations and earlier exposure to an anti-tumor necrosis factor (TNF) drug is encouraged. Observational studies have shown that ERA is associated with a lower likelihood of good response to TNF inhibitors compared to other categories of juvenile idiopathic arthritis (JIA)5,6. The comparative efficacy of csDMARD versus bDMARD therapy in children with a new diagnosis of ERA remains unclear.
In this retrospective study, we used a repeated measures design to evaluate the effect of bDMARD therapy compared to csDMARD therapy on relevant clinical and patient-reported outcomes in children from 5 centers with a new diagnosis of ERA. Treatments were based upon provider and family preferences, as per routine clinical practice. The use of state-of-the-art comparative effectiveness analytic methods enabled assessment of bDMARD effects in this cohort with appropriate adjustment of time-invariant and time-variant confounders.
MATERIALS AND METHODS
This study was approved by the committees for the protection of human subjects at each of the participating institutions. Children’s Hospital of Philadelphia served as the coordinating center (IRB 12-009267).
Study sites and participants
The source population for this study was all children who fulfilled the International League of Associations for Rheumatology (ILAR) criteria for ERA7 and had at least 6 months of documented followup at a rheumatology clinic at one of the following academic tertiary care referral centers: Children’s Hospital of Philadelphia (Philadelphia, Pennsylvania), Children’s of Alabama (Birmingham, Alabama), Cincinnati Children’s Hospital Medical Center (Cincinnati, Ohio), Texas Scottish Rite Hospital for Children (Dallas, Texas), and Meyer Children’s Hospital (Florence, Italy). Children who met ERA criteria, but had a first-degree relative with psoriasis were not excluded (n = 7). Children and adolescents were excluded if they transferred care from another institution or were already receiving systemic therapy (csDMARD or anti-TNF) at the time of initial evaluation.
Each institution searched its respective clinical databases for all children diagnosed with ERA in the outpatient health record at the initial or subsequent followup visit. The range of diagnosis dates included from each institution varied depending upon availability of searchable medical records and are as follows: Children’s Hospital of Philadelphia 2001–2012, Children’s of Alabama 2007–2012, Cincinnati Children’s Hospital Medical Center 2007–2012, Texas Scottish Rite Hospital for Children 1993–2011, and Meyer Children’s Hospital 1995–2012. All inclusion and exclusion criteria were verified by the coordinating center using the JIA Calculator8. The JIA Calculator is a Web-based tool to help algorithmically classify children according to ILAR criteria7; 39 children were excluded after this process. The disease characteristics and treatment approaches for an expanded selection of this cohort (including additional patients with an ERA diagnosis, but limited followup) have been previously described9.
Clinical characteristics
Baseline visit was defined as the first rheumatology visit at which the child presented with clinical signs of juvenile arthritis (enthesitis, arthritis, acute uveitis, or inflammatory back pain). The following clinical data were abstracted from the medical record: demographics, family history of HLA-B27–associated disease, clinical features [including development of new sacroiliitis demonstrated on magnetic resonance imaging (MRI)], patient-reported outcomes (disease activity assessment and pain), and medication use. Medications evaluated included NSAID, intraarticular glucocorticoid injections, oral glucocorticoids, csDMARD (MTX, LEF, SSZ), and bDMARD therapy with TNF-α blocking agents (ADA, ETN, IFX). Disease activity at each visit was measured using the Juvenile Spondyloarthritis Disease Activity Index (JSpADA)10, and the clinical Juvenile Arthritis Disease Activity Score (cJADAS10)11. The JSpADA is a validated composite measure (composed of arthritis, enthesitis, patient pain assessment, inflammatory markers, morning stiffness, clinical sacroiliitis, uveitis, back mobility) developed specifically for children with juvenile SpA that ranges from 0 to 8, with higher scores indicating more disease activity. The cJADAS10 (comprising the active joint count, physician’s global, and parent global disease activity evaluation) is also a validated composite disease measure, but it was specifically developed for JIA and ranges from 0 (inactive disease) to 30 (highest disease activity).
Statistical analysis
To determine which clinical factors influence the decision to treat with anti-TNF therapy, we fit a multilevel mixed-effects logistic regression model with adjustment for clustering by patient and site. Covariates tested included study day (baseline visit = Day 1), age, sex, race, HLA-B27 status, csDMARD use, glucocorticoid use (yes/no), and presence of hip arthritis, wrist arthritis, sacroiliitis, or uveitis. The model included all visits up to and including the first visit in which a bDMARD was prescribed. We used locally weighted scatterplot smoothing (Lowess) to also visually evaluate how disease activity scores influenced the probability of being prescribed anti-TNF medication. JSpADA scores were plotted against whether a bDMARD was prescribed at that visit.
Marginal structural models (MSM) were used to estimate the causal effect of anti-TNF treatment12. This approach appropriately estimates time-dependent treatment effects in the presence of time-dependent confounders that are themselves affected by previous treatment and also predict the subsequent treatment. The model is fitted in a 2-stage process. First, at each timepoint each subject’s probability of receiving their own treatment and probability of being censored are derived as inverse-probability-of-treatment weights (IPTW) and inverse-probability-of-censoring weights (IPCW), respectively, using pooled logistic regressions. Second, the association between the treatment and outcome measured repeatedly is evaluated in a generalized estimating equation (GEE) model that is weighted using the product of IPTW and IPCW. Weighting by the IPTW in effect creates a pseudo-population in which no confounding exists. Therefore, the estimated treatment effect from the subsequent regression models based on this population can be interpreted as the true causal effect of the treatment on the outcome. Weighting by the IPCW further accounts for the bias due to any loss of followup, which is common in longitudinal studies.
The pooled logistic regression was used to obtain the IPTW at each visit as the conditional probability of receiving the anti-TNF therapy given the past treatment history, the baseline covariates, and time-varying clinical variables that might influence the receipt of treatment. The baseline variables included the demographics sex, race, and age; time-invariant clinical variables included HLA-B27 status and year of diagnosis; and the time-varying clinical variables included days since diagnosis, glucocorticoid use, prior anti-TNF prescription, hip arthritis, wrist arthritis, and acute anterior uveitis. Similarly, the IPCW were obtained as the conditional probability of not receiving treatment (being censored) at each visit. Stability of the weights was assessed graphically at intervals of 60 study days.
For the second step of weighted MSM, the primary outcome was the active joint count measured over time. Secondary outcomes assessed included repeated measures of the tender entheses count, JSpADA index, cJADAS10, patient assessment of disease activity, and patient-reported pain. An exchangeable correlation structure was used in the GEE model, equivalent to assuming random intercept among individuals to allow for potentially different baseline values. A negative binomial distribution with log link was assumed to account for the overdispersion for active joint count and tender entheses count because of the number of zero counts. The trajectories of the outcome variables over time were visualized by making spaghetti plots, confirming the linear assumption of our model.
Normal distribution with identity link was used for the remainder of outcomes. The models for active joint count and tender entheses count included the following variables: anti-TNF use, csDMARD use, age, sex, study day, an interaction between study day and bDMARD use, and HLA-B27 status, and accounted for clustering within site and weighted inverse probability of bDMARD use (from Step 1). The models for the JSpADA index, cJADAS10, patient assessment of disease activity, and patient-reported pain included the same variables with the exclusion of HLA-B27 status.
Data regarding development of new sacroiliitis was only available for 4 sites. The association of development of new sacroiliitis and anti-TNF use at these 4 sites was tested using chi-square test.
Given the retrospective design of the study, patient-reported outcomes and laboratory tests were collected at the discretion of the site. Missing data were not imputed. One site did not have reliable patient-reported outcome collection prior to March 2010, for which the MSM analysis was restricted to children diagnosed after 2010. Another site did not collect patient-reported disease activity, and this site was excluded from that particular MSM analysis.
All analyses were performed using Stata 14 (StataCorp. 2015, Stata Statistical Software, StataCorp LP) and SAS software 9.4 (2011, SAS Institute Inc.).
RESULTS
Subjects
During the study period, 217 newly diagnosed patients with ERA had a total of 965 outpatient visits during the first year after disease diagnosis. Followup for 1 year ± 2 months after baseline was done for 143 children (65.9%). Median followup for all patients was 335 days [interquartile range (IQR) 280–365 days]. The median age for the cohort was 11.6 years and the median symptom duration at the time of baseline visit was 6 months (IQR 3–12 mos). Demographics and clinical characteristics of the cohort at the visit in which ILAR criteria were met are shown in Table 1.
During the first 6 months following the baseline visit, children presented with arthritis and enthesitis (n = 176, 81.1%), arthritis plus 2 or more additional ILAR criteria (n = 22, 10.1%), or enthesitis plus 2 or more additional ILAR criteria (n = 19, 8.8%). There were 198 (91.2%) children who met the ILAR criteria for ERA at the baseline visit; the remainder fulfilled the ILAR criteria by 6 months. The population was predominantly male (71.9%) and 59.7% were HLA-B27–positive. Sixty-four (29.5%) had a polyarticular course. Thirty-two of 142 (23%) and 49 of 128 (38%) with at least 1 calculable score achieved a JSpADA of 0 or a cJADAS10 of ≤ 1 at some point during the first year after diagnosis, indicating inactive disease. The median time to attain a JSpADA of 0 or a cJADAS10 of ≤ 1 was 232 days (IQR 128–339) and 161 days (IQR 72–287), respectively. A total of 127 children (58.5%) achieved a simultaneous active joint and tender entheses count of 0 a median of 145 days after the baseline visit.
The treating provider performed imaging for sacroiliitis based on clinical suspicion (MRI = 65; radiograph = 25; MRI and radiograph within 90 days of each other = 16; 2 radiograph first, 1 MRI first, and 13 simultaneous). Twenty-one of the children (32%) who had an MRI performed had evidence of sacroiliitis on imaging at diagnosis. An additional 14 children developed sacroiliitis (defined by MRI) over the first year of followup. Of the 35 children with MRI-defined sacroiliitis, 10 also had a radiograph performed within 90 days of the MRI, 1 of which was abnormal. Of the 35 with sacroiliitis at some point during the first year, 21 (60%) were HLA-B27+.
Medication use
Children were treated with anti-TNF monotherapy (n = 33, 15.2%), csDMARD monotherapy (n = 73, 33.6%), or simultaneous csDMARD and anti-TNF therapy (n = 52, 23.9%) the first year after diagnosis. Two patients (1%) were switched from a csDMARD to an anti-TNF without any overlap in medication use. Children who received anti-TNF therapy received ADA (n = 17, 19.5%), ETN (n = 63, 72.4%), or IFX (n = 7, 8.1%) as their primary medication. There were no bDMARD drugs other than anti-TNF drugs prescribed for children with newly diagnosed ERA during the study interval. Median time to first bDMARD prescription was 35 days (IQR 0–99 days). Seventy-five children (86%) who started an anti-TNF agent continued that therapy for the duration of followup. Eight (9%) anti-TNF users switched anti-TNF drugs during the course of therapy. The median time to anti-TNF switch was 155 days (IQR 62–163). Of the 8 children who stopped the anti-TNF drug before the end of followup, 1 had previously tried a different anti-TNF drug.
The severity of disease and certain disease manifestations were associated with a higher probability of physicians to prescribe an anti-TNF drug. As disease activity increased from low to high as measured by the JSpADA index, the probability of receiving an anti-TNF drug increased from < 10% to > 20%, a finding consistent with confounding by indication (Figure 1)13. In a multilevel mixed-effects logistic analysis, the presence of hip, wrist, or sacroiliac arthritis were all associated with increased odds of anti-TNF exposure (Table 2).
The csDMARD prescribed for children with ERA included MTX (n = 100, 78.7%), SSZ (n = 26, 20.5%), and LEF (n = 1, 0.8%). Median time to first csDMARD was 0 days (IQR 0–63 days). Of the 127 children treated with any csDMARD, 27 (20%) discontinued csDMARD use before the end of followup. One hundred and eighty-eight (86.6%) and 58 (26.7%) were treated with nonsteroidal antiinflammatory drugs (NSAID) or systemic glucocorticoids, respectively. Sixty-two (28.6%) received at least 1 joint injection. The median number of joint injections in patients who underwent the procedure was 1 (IQR 1–2).
Outcomes
There were missing values on 38.3% and 34.2% of visits across all sites for patient disease activity assessment and patient-reported pain, respectively. No difference in age, sex, or median active joint count was observed between visits with and without patient-reported outcomes, or between visits with and without missing inflammatory markers.
Results of the MSM model, which adjusted for the confounding and censoring through weighting, are shown in Table 3. There was a statistically significant difference in the primary outcome, active joint count over time, in children who received an anti-TNF drug versus those who did not, when holding other covariates constant (p = 0.03). Over the first year after disease, patients who received an anti-TNF drug also reported less pain (p < 0.01) and had improved disease activity over time as assessed by the cJADAS10 (p < 0.01). Use of an anti-TNF drug was associated with improvement, albeit statistically insignificant, in all remaining clinical, disease activity, and patient assessments.
Use of a csDMARD, holding all other covariates constant, was associated with a significantly lower tender entheses count (p = 0.02). The use of csDMARD, similar to anti-TNF drug use, was associated with improvement, albeit statistically insignificant in all other outcomes. The magnitude of the estimate, however, was dampened for all outcomes except tender enthesis count in comparison to the estimate for anti-TNF use.
Fourteen children were diagnosed with sacroiliitis by imaging over the course of followup. Twelve (86%) of them were not being treated with an anti-TNF drug at the time of sacroiliitis diagnosis (p < 0.01). Three of these children were subsequently given bDMARD therapy and all children were treated with some form of medication. When stratified by HLA-B27 status, lack of anti-TNF exposure remained significantly associated with development of new sacroiliitis (both p < 0.01).
The change over time in all outcomes in children treated with anti-TNF drug plus csDMARD, anti-TNF monotherapy, csDMARD monotherapy, and supportive care only (NSAID, glucocorticoids, intraarticular joint injections) is shown in Figure 2 and Figure 3. For all outcomes, the rate of improvement was greatest for those children treated with both an anti-TNF drug and a csDMARD, followed by anti-TNF monotherapy.
DISCUSSION
This large multicenter comparative effectiveness study of clinical and patient-reported outcomes in children with newly diagnosed ERA revealed that anti-TNF exposure is associated with statistically significant improvements in active joint count, the cJADAS10, and patient-reported pain over the first year after disease diagnosis. Further, the direction of our estimates was consistent across all outcomes measures. We found that csDMARD therapy, as expected, also improved outcome measures. The magnitude of estimated effect, however, was uniformly greater in children treated with an anti-TNF drug versus a csDMARD.
Several findings warrant additional discussion. First, as with any observational study of therapeutic intervention, the possibility of confounding by indication bias must be considered. This bias arises when children with more severe disease manifestations are more likely to receive the exposure of interest and experience poorer outcomes14. In our study, we did find evidence that children with higher disease activity and more severe disease manifestations (hip, wrist, axial arthritis) were more likely to receive bDMARD therapy within the first 3 months, as demonstrated in Table 2 and Figure 1. Our MSM model, which adjusted for confounding and censoring through weighting, likely minimized but did not completely remove this bias. Because we demonstrated that the children who received bDMARD had a greater magnitude of beneficial effect than children who received csDMARD and that the rate of improvement over time was greatest in children who received an anti-TNF drug with or without a csDMARD, the possibility exists that anti-TNF agents have an even greater positive effect on clinical and patient-reported outcomes than we were able to demonstrate.
Second, in this multicenter cohort, more than one-third of children attained a cJADAS10 indicating inactive disease activity a median of 161 days after diagnosis. This proportion of responders is in accordance with values previously reported in an observational study of children with ERA (38% after 15 months of therapy)6. In the aforementioned study, the Wallace criteria for inactive disease15 were used, which were developed for use in other categories of JIA. In another observational study, 43% of children with ERA attained inactive disease according to the Wallace criteria during a 1-year followup period5.
Third, our study was not designed to systematically evaluate for the presence of axial arthritis. Imaging for suspicion of axial disease was performed as per the treating physician. Interestingly, 11% of children in the cohort had axial involvement recognized on MRI evaluation at some point during the first year after disease diagnosis. Of these, 86% were not being treated with an anti-TNF drug. Whether early bDMARD use was protective against development of axial arthritis in those treated with an anti-TNF drug was unclear, as was whether early use of an anti-TNF drug suppressed axial disease symptoms and, therefore, the need for subsequent imaging. Prior studies have shown that in children with newly diagnosed juvenile SpA and MRI evidence of sacroiliitis (both active and chronic lesions), up to two-thirds may not have back pain16. Without the use of universal screening to detect subclinical sacroiliitis, the true efficacy of anti-TNF drugs and other bDMARD agents for this disease manifestation will remain unknown. The role of early bDMARD use in juvenile SpA, including ERA, remains unclear and has not been systematically evaluated.
This study supports the effectiveness of anti-TNF drugs within the first year after disease diagnosis in routine clinical practice. Children treated with bDMARD had improvement in all clinical features and patient-reported outcomes, although some results were statistically insignificant. Efficacy trials are critically needed of early bDMARD use versus traditional csDMARD to assess the effect on time to inactive disease, risk and treatment of sacroiliitis, patient-reported outcomes, and cost implications.
Footnotes
Dr. Weiss’s work was supported by the US National Institutes of Health/National Institute of Arthritis and Musculoskeletal and Skin Diseases (grant 1-K23-AR-059749-01A1).
- Accepted for publication July 6, 2017.
REFERENCES
- 1.
- 2.
- 3.
- 4.
- 5.
- 6.
- 7.
- 8.
- 9.
- 10.
- 11.
- 12.
- 13.
- 14.
- 15.
- 16.