Article Text

Download PDFPDF

Trajectory of intensive treat-to-target disease modifying drug regimen in an observational study of an early rheumatoid arthritis cohort
  1. Douglas White1,
  2. Helen Pahau1,
  3. Emily Duggan1,
  4. Sanjoy Paul2,
  5. Ranjeny Thomas1
  1. 1Diamantina Institute, University of Queensland, Princess Alexandra Hospital, Ipswich Road, Woolloongabba, Queensland, Australia
  2. 2Queensland Clinical Trials and Biostatistics Centre, School of Population Health, University of Queensland, Translational Research Institute, Princess Alexandra Hospital, Brisbane, Queensland, Australia
  1. Correspondence to Professor Ranjeny Thomas; ranjeny.thomas{at}uq.edu.au

Abstract

Objectives Studies of early rheumatoid arthritis (RA) cohorts have analysed treatment response and prognostic factors at fixed time points. However, in treat-to-target protocols, therapeutic decision-making is dynamic and responsive to disease activity over time. To determine when a minimal residual disease response target should be expected, our primary objective was to identify the time-dependent therapeutic response to combination disease modifying antirheumatic drugs (DMARDs) for 12 months. Our secondary objective determined factors affecting this response trajectory.

Design Observational cohort.

Setting Treat-to-target early RA clinic in Australian tertiary referral hospital.

Participants We enrolled consecutive patients attending an early arthritis clinic with symptom duration less than 12 months, who were diagnosed with RA for the first time between 2004 and 2008. 101 met these eligibility criteria and data were available at baseline through 12 months.

Interventions intensive DMARDs according to a treat-to-target protocol.

Primary and secondary outcome measures We measured disease activity scores (DAS) at each visit, then analysed therapeutic response and associated factors in a time-dependent fashion over 12 months.

Results The median DAS4vESR of 4.46 at baseline decreased 12 weeks later by 24%, while the proportion with DAS4v ≤ 2.6 increased (p<0.01). DAS4v continued to decrease over 52 weeks. DAS4v reduction of at least −0.45 at 4 weeks was predictive of DAS4v at 28 and 52 weeks. Female gender, current smoking, primary education and an interaction between baseline weight and C reactive protein (CRP) negatively impacted DAS4v reduction over 4 and 52 weeks. Time-varying effects of blood pressure, neutrophils, erythrocyte sedimentation rate and CRP also significantly influenced DAS4v over 52 weeks.

Conclusions Time-dependent data suggest that the largest reduction of DAS4v to combination DMARDs occurs in the first month of therapy, and this predicts subsequent response. Variables known to impact long-term treatment response in RA also impacted early DAS4v response to combination DMARDs.

  • RHEUMATOLOGY
  • EPIDEMIOLOGY
  • STATISTICS & RESEARCH METHODS

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 3.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/3.0/

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Article summary

Article focus

  • Best-practice early rheumatoid arthritis (RA) treatment aims to achieve a target response. In clinical settings of many countries, first-line therapies are disease modifying antirheumatic drugs (DMARDs), including combination DMARDs.

  • We followed an observational cohort for 12 months in a treat-to-target early RA clinic to identify the time-dependent therapeutic response to combination DMARDs for 12 months and factors affecting this response trajectory.

Key messages

  • After initiation of combination disease modifying antirheumatic drugs (DMARDs), the largest reduction in disease activity score occurred in the first month, and its magnitude predicted subsequent response.

  • Disease activity score over 12 months was influenced by female gender and current smoking, education level and an interactive effect of weight and either C reactive protein or erythrocyte sedimentation rate.

  • The data suggest clinical response to combination DMARDs may be more rapid than previously appreciated, and treatment response in the first month may have prognostic significance.

  • These hypotheses require further testing in other cohorts.

Strengths and limitations of this study

  • Monthly observation allowed precise determination of time-dependent therapeutic response and demonstrated an unexpectedly rapid response to combination disease modifying antirheumatic drugs (DMARDs).

  • Standardised combination DMARD treat-to-target protocol.

  • Real-life clinical setting with dynamic therapeutic decision-making.

  • Observational cohort study limits conclusions that can be drawn regarding causality, without further testing in a randomised controlled trial.

  • Relatively small cohort derived from a single centre, with treatment regimen determined within Australian prescribing context, and exclusions due to missing data limit generalisability.

  • Number of participants limited by lack of baseline or 12-month follow-up data and may have introduced selection bias.

  • Owing to incomplete radiographic data, factors associated with radiographic outcomes could not be determined.

Background

Intervention with early combination disease modifying antirheumatic drug (DMARD) therapy favourably influences progression of rheumatoid arthritis (RA) independent of treatment in later years, suggesting that there is a ‘window of opportunity’ in which the disease process can be altered.1 ,2 Moreover, a good response at 6 months to tight disease control using methotrexate predicted outcome after 5 years of treatment in participants in the CAMERA study.3 The severity of disease varies in RA patients. In those with aggressive disease, damage to articular structures occurs early in the disease process: erosions were detected in 12.8% of patients after a median of 8 weeks in one study.4 Thus, early evidence and determinants of treatment response to a given regimen are critical, in order to channel patients at greatest risk of poor outcome to more intensive induction regimens or more expensive biological therapies within that window.

Studies of prognostic factors by statistical modelling have analysed disease progression outcomes including erosions, disease activity score (DAS28) and disability index as measured by Health Assessment Questionnaire (HAQ) at fixed time points—usually 6 or 12 months, with the earliest being 3 months—to determine treatment response and associated factors influencing this. Factors associated with poor radiological outcome include smoking, rheumatoid factor (RF) positivity, the presence of anticitrullinated peptide autoantibodies (ACPA), HLA-DR genotype, low socioeconomic status and bone oedema on MRI.5–9 However, poor outcome measured by HAQ was associated with high baseline disease activity or HAQ, including RF, DAS28 score, tender and swollen joint counts, ESR and CRP.10 ,11 However, in treat-to-target protocols, such as was used in the TICORA trial and which occur in real-life clinic settings, decision-making about dose and drugs is a dynamic process, responsive to the patient's disease activity over time.12 In many early arthritis protocols, including the current study, patients are treated and monitored intensively during the first 3–6 months, followed by a reduced visit frequency. Longitudinal analysis of all available data, while modelling the trajectories and drawing inferences on the significance of various risk factors, provides higher power and better insight into the dynamic process.

Aims

In the current study, our primary objective was to identify the time-dependent therapeutic response in an observational study of combination DMARDs for 12 months in order to determine when a minimal residual disease response target should be expected. Our secondary objective was to determine factors affecting this response trajectory. We therefore gathered disease activity data at each treatment visit then analysed the disease activity response in a time-dependent fashion. We then determined factors which influenced this time-dependent response to an intensive DMARD regimen.

Methods

We enrolled consecutive patients referred by general practitioners from a relatively socioeconomically disadvantaged catchment (60% referrals of employed individuals working in manual industries) to an early arthritis clinic in a public teaching hospital, with symptom duration less than 2 years, who were diagnosed with RA for the first time between 2004 and 2008. Patients were selected for inclusion in the current study if data were available at baseline through 12 months; however, data were not required at every time point for inclusion. Two hundred and six patients were referred with possible RA and 101 patients met these eligibility criteria; 49 did not have RA, and 54 patients who met all other criteria were excluded as data were unavailable at 12 months. Of these, seven were seen once and diagnosed with RA then treated elsewhere, and the remainder were reviewed at least once but not at 12 months. All study participants met the American College of Rheumatology 1987 revised criteria for the classification of RA.13

Referrals from local general practitioners were triaged within 1 week, and patients were generally diagnosed within the next 4 weeks. Since full clinical and laboratory evaluation was available at the first visit to the early arthritis clinic, patients received combination methotrexate (MTX), sulfasalazine (SSZ) and hydroxychloroquine (HCQ),14 unless contraindicated, immediately after diagnosis and confirmation of RA by the treating rheumatologist. Treatment was intensified according to a response-driven step-up algorithm, as previously described,15 with remission as the target.16 ,17 Briefly, criteria for dose escalation were either >2 swollen joints and abnormal ESR or CRP, or at least two of the following four criteria: morning stiffness >30 mins, pain or fatigue visual analogue scale (VAS) >30 mm, or >2 tender joints. The following medications were prescribed at baseline: MTX 10 mg/week, folic acid 5 mg/week, SSZ 500 mg daily increasing by 500 mg at weekly intervals to 1000 mg twice daily, HCQ 200 mg daily for 1 week then 400 mg daily thereafter. Patients were seen at 4-weekly intervals and the MTX dose was escalated according to treatment response at a conservative rate by 5 mg at each visit to a maximum of 25 mg weekly. If disease remained active on this combination, SSZ was stopped, MTX reduced to 10 mg weekly and leflunomide started at a dose of 20 mg daily. MTX dose was titrated back to 25 mg weekly, and if this combination failed and the Australian Pharmaceutical Benefit Scheme criteria were met, the patient started biologic therapy. Based on these criteria, 10% of patients in this setting started biologics per year. In general, the use of non-steroidal anti-inflammatory drugs and oral corticosteroids was minimised, but intra-articular or oral steroids could be administered at the discretion of the treating physician. Large joints were injected with 40–80 mg DepoMedrol and smaller joints with 1 mL (5.7 mg) Celestone. Oral and intra-articular dosage of corticosteroids was recorded monthly.

Response to therapy was measured as follows: the four variable DAS28ESR (hereafter referred to as DAS4v) was used as an index of inflammatory control,18 and the mHAQ as an index of disability.19 Each index was calculated at each visit: baseline, and weeks 4, 8, 12, 16, 22, 28, 36, 44 and 52.

Demographic details were ascertained by questionnaire and included: age at presentation, symptom duration, level of education, gender, current, ex-smokers and non-smokers and mHAQ. Patients completed VAS for pain, fatigue and their global assessment of disease. The 28 tender and swollen joint counts, height, weight and blood pressure (BP) were recorded by the clinical research nurse. Blood was collected at baseline for the analysis of ESR, CRP, LFT, FBC, RF and ACPA. ACPA were measured at Queensland Health Pathology using the anti-CCP2 ELISA (Axis-Shield) test, with the cut-off of six for a positive test.

Basic statistics were presented by number (%) or mean (SD) or median (IQR), as appropriate. Five imputations for missing data on clinical, biochemical and score data were performed using Bayesian-Markov chain Monte Carlo multiple-imputation technique. Of those who met eligibility criteria for study inclusion, not all patients attended for all visits, however the patterns of missingness were random for all the study parameters. The consistency in the distributions of the five imputed data was checked for all study parameters. Given the skewed DAS and mHAQ scores, the medians and their 95% CI are presented. The changes in these scores over the study period are presented by median and 95% CI. Significance levels (p values) are based on the appropriate non-parametric test.

Generalised multivariate linear regression models with γ distribution and identity link were used to identify the statistically significant (p≤0.10) risk factors and their possible interaction effects on disease activity scores at week 4 of the study. The possible consistency in the effect sizes of the statistically significant risk factors (at week 4) were also assessed on the disease activity scores at week 12 of the study. Combining the 10 longitudinal measurements obtained over 1 year of the study, the time varying effects of individual risk factors on the disease activity scores were explored using generalised estimating equation regression approach with γ distribution and identity link function under the assumption of ‘unstructured’ correlation structure.

Results

Time-dependent therapeutic response to combination DMARDs for 12 months

One hundred and one patients were included in the study and 54 (of whom 7 were only seen once) were excluded due to missing 12 month follow-up data. The baseline characteristics of the included and the excluded patients are shown in table 1. Except for a lower systolic BP in excluded subjects, there were no significant differences between included and excluded subjects. All except four patients took at least two and up to three of the following DMARDs in combination during the 12 months study: methotrexate, sulfasalazine, hydroxychloroquine and leflunomide. These four patients took methotrexate monotherapy.

Table 1

Baseline characteristics of the patients

The median disease activity score at baseline was 4.46 for DAS4v (table 2). Four of the 12 patients with baseline DAS≤2.6 (minimal disease activity) were taking steroids prior to referral. There was a highly significant (p<0.001) DAS reduction of 24% at 28 weeks (table 2). We also observed a significant increase in the proportion of patients with minimal residual, (DAS28≤2.6) and low disease (DAS28 s≤3.2) over the treatment period (p<0.01) (table 2). Consistent with this, the patients’ pain scores improved highly significantly by 31% and 56% at the end of 6-month and 1-year of treatment, respectively. The improvement in mHAQ from baseline to 6 months, but not between 6 and 12 months of treatment, was significant. The average annual change was 0.30 units (table 2).

Table 2

Change in median pain VAS scores, DAS and HAQ scores (95% CI) over 1 year

Analysis of the change in DAS4v over time showed a progressive reduction over 52 weeks, with the steepest drop between baseline and 4 weeks (figure 1). The median (95% CI) of changes in DAS4v at 4, 28 and 52 weeks were −0.45 (−0.84, −0.07), −0.86 (−1.30, −0.41) and −1.35 (−1.67, −1.03), respectively (p<0.01 at week 52). The changing patterns of the distribution of DAS4v over time are evident from the density plots in figure 2. Although a significant shift in the distribution of DAS4v at 4 weeks from baseline is evident from the density plot, the distributions overlap at 4, 28 and 52 weeks.

Figure 1

Distribution of DAS4v over the study period. (A) The median and 95% CI are plotted for each visit over the 52 week study period. (B) Changes in DAS4v over 4, 28 and 52 weeks are indicated. (C) The changing distribution in DAS4v in the sample is plotted at baseline, 4, 28 and 52 weeks.

Figure 2

Variation in the disease activity parameters over the study period. The median and 95% CI are plotted for each visit over the 52-week study period for erythrocyte sedimentation rate, C reactive protein, tender joint count, swollen joint count, fatigue, morning stiffness, patient global and physician global scores.

Analysis of the individual components of the DAS over this period demonstrated that patient global score, swollen and tender joint counts all fell most steeply between baseline and 4 weeks (figure 2). While this was not the case for the fall in either ESR or CRP, similar steep falls in fatigue score, morning stiffness and physician global scores occurred between baseline and 4 weeks. Thus most measures of disease activity fell most rapidly in the first 4 weeks after DMARD initiation. In contrast, ESR fell for 3 months before reaching a plateau, while CRP fell progressively for 6 months.

Factors affecting the response trajectory in early RA patients treated with combination DMARDs

To determine whether the fall in DAS4v at 4 weeks predicted the DAS at 28 and 52 weeks, we first calculated that the median level of change in DAS4v at 4 weeks was −0.45. This was clinically discriminatory: at 4 weeks, 52% had no change or an increase in DAS4v while 48% improved from baseline DAS4v. The number and proportion of patients receiving steroids is indicated in table 3. While baseline steroids impacted the likelihood of improvement at 4 weeks, this was not statistically significant (69% of patients receiving steroids improved and 53% not receiving steroids improved; OR for improvement with steroids 1.95, p=0.12). Patients with reduction in DAS4v at 4 weeks of at least −0.45 were three times more likely (OR (95% CI) 3.10 (1.2 to 8.0)) at 28 weeks and 17 times more likely (OR (95% CI) 17.14 (4.52 to 64.94)) at 52 weeks to maintain the same or reduced DAS4v as achieved after 4 weeks of treatment. Univariate modelling of factors affecting outcome showed that female sex, smoking status and increasing alanine aminotransferase (ALT) at baseline negatively affected DAS4v at 4 weeks, but these effects became less significant by 12 weeks (table 4). An interaction between baseline weight and CRP negatively affected DAS at both week 4 and 12. Patients taking steroid did not have a significantly different disease score, and symptom duration before RA diagnosis, anti-CCP or RF titre did not impact 4-week DAS. The reduction in DAS4v at 4 weeks was significantly greater in patients with tertiary than with primary education.

Table 3

Frequency of steroid use over the study

Table 4

Variables influencing DAS scores at 4 and 12 weeks of study—univariate regression

Combining baseline characteristics and the longitudinal measurements obtained over 1 year, we explored the time-varying effects of individual risk factors on DAS4v in a univariate model (table 5). DAS4v over 52 weeks was again influenced by female gender and current smoking, and an interactive effect of weight and either CRP or ESR. Time varying effects of systolic and diastolic blood pressure, neutrophil counts, ESR and CRP also significantly influenced DAS4v observed over 52 weeks (table 4). At week 52, the largest reduction in DAS4v was observed in patients with tertiary education (3.57), compared with that observed among patients with secondary (2.56) or primary education (1.33). Symptom duration prior to diagnosis did not significantly influence DAS4v over 52 weeks. Over the course of the study, DAS4v was increased by 0.66 in those patients taking steroids (p<0.01). These data are, in keeping with the use of steroid in this study at the clinician's discretion, to provide additional control for disease activity that was not controlled by the DMARD protocol.

Table 5

Effects of time-varying risk factors individually on DAS28 scores over 1 year of study—univariate regression

We found that the relationship between mHAQ and DAS4v for the cohort was significantly correlated at baseline, 4, 28 and 52 weeks (p<0.001), with this correlation becoming progressively tighter over time as DAS and mHAQ fell. Thus, functional outcome after 1 year of early RA treatment is highly dependent on achievement of low disease activity.

Discussion

Our study describes the response of a group of patients with early RA to intensive conventional DMARD therapy in a time-dependent fashion over the first year. Baseline characteristics were in line with previous cohorts of patients with early RA. However, our baseline DAS scores were relatively low, reflecting our rapid triage and clinical and laboratory assessment of early arthritis referrals. Surprisingly, the time-dependent analysis of DAS response showed that the majority of disease activity measures fall most rapidly in the first 4 weeks after starting intensive DMARD treatment in this population. There was a subsequent slow and progressive reduction in DAS until week 52. This fall in DAS4v at 4 weeks appeared to be clinically meaningful, as it predicted the DAS at 28 and 52 weeks. This observation suggests that for patients who failed to respond within 4 weeks to combination DMARD treatment, few gains were made by continuing to apply the same DMARD treat-to-target algorithm for 6–12 months. This was reflected in the similar proportion of patients with minimal or low disease activity between 6 and 12 months. By this stage patients had progressed through the combination DMARD algorithm, for which the next step would be biologics. However, because their disease activity is minimal or low, they failed to qualify for biologics based on Australia's Pharmaceutical Benefit Scheme requirements (http://www.medicareaustralia.gov. au/provider/pbs/drugs2/rheumatoid.jsp).20 On the other hand, our data suggest the hypothesis that continued effort in applying a treat-to-target combination DMARD algorithm is likely to be effective over the ensuing months in patients who make a moderate response by week 4. Our data further suggest that combination DMARDs act unexpectedly rapidly in this early RA population, as patients’ use of steroids did not influence the reduction in DAS. In support of this conclusion regarding steroids, in a study of 61 patients with early RA treated according to a similar response-driven step-up combination DMARD algorithm, Proudman et al17 obtained an almost identical 6 month minimal disease activity rate (DAS28<2.6 in 29%), despite infrequent use of corticosteroids.

The current study has a number of limitations. First, our interpretation that the magnitude of the fall in DAS4v after 1 month predicts 1 year outcome is limited by the observational study design. However, the question of whether outcome could be improved in patients with a minimal treatment response within 1 month could be tested in a randomised controlled trial comparing switch with biologic therapy with continued combination DMARDs. Second, this is a relatively small cohort derived from a single centre with referrals derived from a relatively socioeconomically disadvantaged catchment, with the treatment regimen determined within the Australian prescribing context. At the time of recruitment, 1987 ACR criteria were used to diagnose RA, which would have limited capacity to diagnose less severe patients. The number of participants was limited by lack of baseline or 12-month follow-up data and this may have introduced selection bias towards a more compliant group. The small sample size and number of exclusions due to incomplete data limit generalisability to other prescribing environments or clinical settings, and further studies are needed to test the generalisability of our findings. For example, it is possible that those excluded had a different disease trajectory due to differences risk for poor outcome or differences in adverse events. A subanalysis of the trajectory excluded patients was not possible because of the low number of paired baseline and 4-week DAS4v measurements in this group. On the other hand, there were no differences in the baseline characteristics of the excluded patients (except systolic BP). Furthermore, almost all factors associated with 4-week DAS4v response have been previously demonstrated to affect disease outcome in longer-term and larger studies. The strengths of this study are that it analyses real-world data, monthly observations allowed precise determination of time-dependent response, and patients received a standardised combination DMARD treat-to-target protocol, reducing the confounding effect of treatment decisions based on individual clinician preference.

The exploratory nature of the study in a relatively small sample could introduce false-positive associations. Although it is possible that the rapid 4-week response to the combination of DMARDs represents regression to the mean, the continued good response of these patients argues against this. Our data also are consistent with recent studies demonstrating that early good response to combination therapy (in the TEAR and RAPID 1 trials) is associated with a continued good response.21 ,22 In these studies, rapid response was ascertained 12 weeks after initiation of combination therapy. By regression analysis, we identified female gender, current smoking, education level, ALT and an interaction between weight and CRP as significant determinants of disease activity over 4 and 52 weeks. Women, current smokers and low levels of education were found in several studies, including those of early RA, to achieve lower reductions in disease activity or remission.23–26 However, no study has previously determined that the impact of these variables may occur within weeks of starting treatment. The interaction between weight and inflammation in RA is intriguing and has been noted previously in insulin resistant states.27 In patients with active RA, those with high body mass index responded less well to infliximab.28 We also identified significant time-varying effects of blood pressure, gender, age, weight and inflammatory markers on disease activity. The interaction between disease activity and cardiovascular risk is well documented in RA, including early RA, and traditional cardiovascular risk factors may also impact the activity of inflammatory disease over time.15 ,29–31 However, it is unknown whether control of cardiovascular risk factors can in turn impact inflammatory disease control.

In this study we were limited to analysis of disease and functional score, as radiographic data were not sufficiently complete to allow measurement of structural damage. However, this issue has been addressed by others, where biomarkers such as ACPA antibodies, RF, CRP and cartilage oligomeric matrix protein can add power to predictive models of bone erosion in early RA.32 In contrast, we found no impact of ACPA or RF on DAS. Our data confirm a strong relationship between disease activity and functional score that appears to strengthen over time, a finding that is supported by data from the BeST cohort.16 We would anticipate that functional disability would be minimised by early treatment with combination DMARDs as shown previously.33 ,34

Since they are traditionally thought to be slow acting, previous studies of DMARD monotherapy in early RA have not analysed time-dependent data from 4 weeks. Although it remains possible that a similar response might be observed in some patients starting DMARD monotherapy, we suggest this rapid response may be a unique feature of intensive combination DMARDs (with multiple mechanisms of action) initiation in early RA, which is the RA population most responsive to therapeutic intervention.35 ,36 The risks and benefits of intensive DMARD therapy (combinations allowing switching to achieve tight control) versus monotherapy in early RA deserve further study, considering inconsistent evidence to support combination DMARD therapy in RA.36 ,37 The need to identify patients with more aggressive disease prompted one group to undertake a trial of a stratified treatment plan based on the likelihood of persistent arthritis, with the aim of minimising over treatment and under treatment in early RA.38 Our data suggest the hypothesis that very early response to an intensive DMARD strategy that minimises under treatment predicts response for the first year.

Data from the ERAN study show that patients with moderate disease activity at 1 year are unlikely to achieve better control of their disease if the same protocol is continued, and a good response at 6 months in the CAMERA study predicted outcome at 5 years.3 ,39 Our data, collected in a cohort of early RA patients with relatively low baseline DAS, likely reflect the trajectory of patients meeting criteria for RA early in disease, and which would be captured in organised clinical settings using the recently published new classification criteria.40

Conclusions

With the availability of increasing numbers of treatment options, application of strategies that identify early non-responders to intensive DMARD combinations, has clear implications for treatment stratification within the window of opportunity. Time-dependent data suggest clinical response to combination DMARDs may be more rapid than previously appreciated, and treatment response in the first month may have prognostic significance. Confirmation in other cohorts will be required to determine the generalisability of this notion.

References

Footnotes

  • SP and RT contributed equally.

  • Contributors DW, SP, RT were involved in conception and design. DW, HP, ED contributed by acquiring, analysing and interpreting the data. DW, HP, ED, SP, RT were involved in drafting the article or revising it critically for important intellectual content, and final approval before submission.

  • Funding The study was supported by NHMRC grants 351439 and 569938. RT is supported by Arthritis Queensland and an ARC Future Fellowship.

  • Competing interests None.

  • Ethics approval Ethical approval for retrospective data analysis was obtained from the Metro South Human Research Ethics Committee.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement No additional data are available.