Article Text

Download PDFPDF

Extended report
Comparative effectiveness of biologic monotherapy versus combination therapy for patients with psoriatic arthritis: results from the Corrona registry
  1. Philip J Mease1,
  2. David H Collier2,
  3. Katherine C Saunders3,
  4. Guo Li4,
  5. Joel M Kremer5 and
  6. Jeffrey D Greenberg3,6
  1. 1Department of Rheumatology, Swedish Medical Center and University of Washington, Seattle, Washington, USA
  2. 2Inflammation Global Development, Amgen Inc., Thousand Oaks, California, USA
  3. 3Epidemiology and Outcomes Research, Corrona, LLC, Southborough, Massachusetts, USA
  4. 4Department of Biostatistics, Axio Research LLC, Seattle, Washington, USA
  5. 5Center for Rheumatology, Albany Medical College and the Center for Rheumatology, Albany, New York, USA
  6. 6Division of Rheumatology, New York University School of Medicine, New York, New York, USA
  1. Correspondence to Professor Philip J Mease; pmease{at}philipmease.com

Abstract

Objectives To characterise the comparative effectiveness of combination therapy (a tumour necrosis factor inhibitor (TNFi) and a conventional synthetic disease-modifying antirheumatic drug (csDMARD) such as methotrexate) and monotherapy (TNFi only) for psoriatic arthritis (PsA) from a large US registry.

Methods The analysis included adult patients with PsA who were enrolled in the Corrona database (ClinicalTrials.gov, NCT01402661), had initiated a TNFi, were biologic naïve, and had a follow-up visit ≥90 days after drug initiation. The endpoints of the analysis were TNFi persistence (drug survival) and time to Clinical Disease Activity Index (CDAI) remission. All analyses were performed using propensity scoring, which were estimated using CDAI and patient sex, to control for channelling bias.

Results Of 519 patients meeting the inclusion criteria (318 with combination therapy and 201 with monotherapy), the analysis population was 497 for TNFi persistence and 380 for time to remission. The difference between combination therapy (TNFi+methotrexate, 91% of patients; TNFi+other csDMARD, 9%) and monotherapy was not statistically significant for TNFi persistence (32 and 31 months, p=0.73) and time to remission (21 and 25 months, p=0.56). Predictors of TNFi persistence included Hispanic ethnicity (longer persistence), PsA duration (longer persistence), history of methotrexate use (shorter persistence), body mass index (shorter persistence) and disease activity (shorter persistence).

Conclusions Patients with PsA from a large US registry experienced similar TNFi persistence on combination therapy and monotherapy. Prospective, randomised clinical trials evaluating the efficacy of combination therapy versus monotherapy would provide much-needed clarity on treatment options for patients with PsA.

Trial registration number NCT01402661.

  • Anti-TNF
  • Psoriatic Arthritis
  • DMARDs (biologic)

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.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/4.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.

Key messages

What is already known about this subject?

  • There is conflicting evidence about the clinical effectiveness of tumour necrosis factor inhibitor (TNFi) monotherapy versus combination therapy (ie, with a conventional synthetic disease-modifying antirheumatic drug such as methotrexate) in psoriatic arthritis. We performed an analysis of TNFi persistence (ie, time to therapy change) in a large, US-based, prospective registry.

What does this study add?

  • We performed an analysis of TNFi persistence (ie, time to therapy change) in a large, US-based, prospective registry.

  • There was no statistically significant difference in TNFi persistence between patients using monotherapy and combination therapy (32 and 31 months) or time to remission (21 and 25 months).

How might this impact on clinical practice?

  • Our analysis indicates that TNFi persistence is similar between monotherapy and combination therapy, a result that is consistent with similar analyses of registries from other countries.

  • Clinical trials should be conducted to verify the findings of this and other observational registries.

Introduction

Psoriatic arthritis (PsA) is a chronic inflammatory disease that occurs in approximately 0.3% of the US population, as suggested by previous studies.1 ,2 However, it has recently been shown that up to 30% of patients with psoriasis may have PsA,3 and recent population surveys show that psoriasis occurs in 3.2% of the US population.4 Therefore, previous estimates of PsA general population prevalence may underestimate the true prevalence. Early intervention is important, because the persistent inflammation can cause progressive joint damage leading to physical limitations and severe disability in some patients.5 ,6 In addition, PsA can increase the risk for cardiovascular complications and death.7–9 Therefore, optimal treatment strategies are important to understand.

Treatment strategies for patients with moderate to severe PsA have historically been based on extrapolation of efficacy data from patients with rheumatoid arthritis (RA). Like in RA, tumour necrosis factor (TNF)-α plays a key role in the pathogenesis of PsA, and a number of clinical trials have shown that TNF inhibitors (TNFis) are effective in patients with PsA.10–14 Coadministration of a conventional synthetic disease-modifying antirheumatic drug (csDMARD), such as methotrexate (MTX), improves the efficacy of TNFis in the treatment of RA15 ,16; therefore, combination therapy is commonly used to treat PsA.17 ,18 PsA is, however, distinct from RA with respect to its pathology, clinical characteristics and outcomes,5 ,17 ,19 ,20 and no randomised controlled clinical trials have prospectively compared the effectiveness of combination therapy versus monotherapy in PsA.

The phase 3 studies of TNFis in PsA allowed patients receiving MTX at time of enrolment to continue MTX, and similar clinical responses were observed among patients who received concomitant MTX therapy versus those who received TNFi as monotherapy.10–14 A number of observational studies of real-world registry data found that higher percentages of patients with PsA receive MTX in combination with a TNFi versus a TNFi alone.21–23 An analysis of PsA data from the South Swedish Arthritis Treatment Group (SSATG) register showed that concomitant MTX did not affect treatment responses.22 Likewise, the Norwegian DMARD (NOR-DMARD) study found no benefit of concomitant MTX.21 There was a small beneficial effect of concomitant MTX on the achievement of American College of Rheumatology 20% (ACR20) response among patients enrolled in the Danish Biologics (DANBIO) registry.23 The effect of concomitant MTX on drug survival varied among the three registry studies, with a trend towards longer drug survival in the SSATG study22 and the DANBIO study (significant after the analysis was adjusted for other baseline variables)23 and significantly longer drug survival in the NOR-DMARD study.21

There is variable evidence regarding the effectiveness of combination therapy versus monotherapy in PsA. Moreover, MTX is the most widely used systemic agent for PsA, whether used alone or in combination with a biologic. The data to date, which include two neutral placebo-controlled studies24 ,25 and some positive open-label studies, indicate that MTX modifies symptoms but not the underlying disease.26 Clinicians must make individualised decisions about whether to prescribe MTX with a biologic in PsA, so we performed a study to provide more data on the comparative effectiveness of combination therapy and monotherapy. The objectives of our study were to examine the effects of these treatment regimens on TNFi persistence (termed drug survival in other studies) and measures of effectiveness using a large, US-based registry: the Corrona registry (ClinicalTrials.gov, NCT01402661). We also investigated predictors of TNFi persistence and treatment response.

Methods

Data source

Corrona is a large, US-based, independent, prospective, observational, disease-based registry initiated in 2001 that has always collected data on patients with RA and PsA. Patients are enrolled in the registry by participating rheumatologists, and questionnaires are completed by both patients and physicians during the visits as part of routine clinical care. Corrona is a prospective observational registry where data collection is not mandated at regular intervals. Physicians are encouraged to complete the questionnaires at least every 6 months (mean interval for current study, 8.6±5.4 months). Corrona is reviewed by the New England Institutional Review Board.

The analysis was performed on data collected between 2005 and 1 October 2012. During that time Corrona had enrolled 5408 adult patients with PsA, comprising 25 865 visits to 54 rheumatologists from 84 sites in 33 states. The mean±SD follow-up for these patients was 2.1±2.2 years (range 0–7.72 years; total of 11 550 patient-years).

Study population

Patients were included in the analysis if they had a diagnosis of PsA; were at least 18 years old; initiated a biologic in 2005 or later, either alone or in combination with a csDMARD; were biologic naïve; and had a follow-up visit at least 90 days after drug initiation. The diagnosis of PsA was according to the physician's designation on the form; the registry existed before the ClASsification criteria for Psoriatic ARthritis (CASPAR) criteria were published in 2006.27 For the analysis of TNFi persistence, patients had a follow-up visit of at least 1 year after drug initiation.

Study objectives

The objectives of the analysis were to estimate effectiveness measured by the length of time that patients with PsA maintained their initial biologic therapy (TNFi persistence), either as monotherapy or in combination with MTX or other csDMARDs (an a priori comparison); time to remission; and effects of patient risk factors on these outcomes.

Change in monotherapy was defined as any change in the initial biologic, including the addition of a csDMARD or changing or stopping the biologic. Change in combination therapy was defined as changing or stopping the biologic. Remission was defined as a Clinical Disease Activity Index (CDAI) score <2.8.

Baseline disease activity

Measurements of baseline disease activity included CDAI, physician global assessment of disease activity and skin, patient global assessment of disease activity and skin, modified Health Assessment Questionnaire (mHAQ), tender and swollen joint counts, erythrocyte sedimentation rate levels, and levels of C reactive protein (CRP).

Statistical analyses

The data sets used in the analyses are shown in figure 1. Patients initiating a non-TNFi biologic were excluded from the analysis owing to the small number of patients (n=12).

Figure 1

Data sets used in the analyses (CDAI, Clinical Disease Activity Index; PsA, psoriatic arthritis).

TNFi persistence (measured as the time from the date of initiation of biologic therapy to the date of change in therapy) was estimated using the Kaplan–Meier method. Kaplan-Meier curves were estimated for the monotherapy and combination therapy groups, as well as the individual TNFis. Log-rank tests were used to compare the monotherapy and combination therapy groups. The Kaplan-Meier method was also used to estimate the median time to therapy change and the proportion of patients who remained on their initial biologic therapy at different follow-up time points.

To control for channelling bias, all patients receiving monotherapy or combination therapy were matched using a propensity score matching method.28 ,29 The propensity score was estimated using CDAI and patient sex. These variables were chosen following an analysis of the data using other possible risk factors (swollen and tender joint counts, CDAI, physician global assessment of disease activity, history of csDMARD use and history of MTX use). CDAI has correlated with many other disease severity measures. Therefore, we felt CDAI would be a good indicator to represent overall disease activity at baseline, which we wanted to balance between the monotherapy and combination therapy groups. We also included patient sex in the model, because we have observed more women initiating combination therapy compared with men. In order to minimise the propensity score missingness due to missing covariates, we limited the covariates to CDAI and gender in the propensity score calculation.

We derived propensity scores for matching patients in the monotherapy and combination therapy groups based on the conditional probability of receiving either treatment. Weights were also constructed using propensity scores from logistic models and were normalised using inverse probability treatment weighting, rather than directly matching the monotherapy and combination therapy groups. Weights information was taken into account in proportional hazards regression models. A separate propensity score was calculated for each analysis of TNFi persistence by TNFi.

Patients who remained on their initial biologic therapy for at least 90 days were included in the analysis of time to remission.

Both univariate and multivariate Cox regression models were used to identify predictors of TNFi persistence and time to remission. In order to control the residual confounding effects after inverse probability treatment weighting, we adjusted covariates that were significant in the univariate models. Age and patient sex were not always significant in univariate analysis, but these variables were always included in multivariate models to adjust potential confounding or risk modification effects. Body mass index (BMI) was not always significant in univariate regression, but it was included in the multivariate model because it is an important clinical risk factor. Race and ethnicity variables were sometimes significant in the univariate model but were no longer significant after adjusting for other variables; therefore, these were not included in the multivariate model. We did not use a strict stay criterion for identifying predictors in the multivariate model; predictors were included in the multivariate model if they were marginally significant in the univariate model (p<0.10), clinically important, or strongly prognostic.

All analyses were done using SAS V.9.1 (SAS Institute Inc, Cary, North Carolina, USA) and Stata V.10.1 (Stata Corp, College Station, Texas, USA). p Values <0.05 were considered statistically significant.

Results

Baseline demographics, disease characteristics and medical history

A total of 519 patients (combination therapy, n=318; monotherapy, n=201) met the inclusion criteria (table 1). Most patients were white; approximately half were women; and the mean age was 51.6 years. Patient demographics and treatment history were similar between the two groups, except that almost twice as many patients in the combination therapy group had initiated treatment with infliximab (table 1). Patients in the combination therapy group had significantly higher baseline tender and swollen joint counts, CDAI values, and physician's assessment of skin values than those in the monotherapy group (table 1). Combination therapy included MTX for 91% of patients and another csDMARD, such as sulfasalazine or leflunomide, for 9%. The number of patients with valid propensity scores and weights in each analysis was 497 for TNFi persistence and 380 for time to remission (figure 1).

Table 1

Characteristics of patients with PsA in the longitudinal cohort of patients registered in the Corrona database

At baseline, 218 patients were taking adalimumab, 158 etanercept and 114 infliximab. Most baseline characteristics within each individual TNFi subgroup were balanced between the combination therapy and monotherapy groups. There was a greater proportion of patients with a history of MTX use in the combination therapy group in each individual TNFi subgroup. Compared to monotherapy, combination therapy was associated with a smaller proportion of women (etanercept only) and patients with alcohol use (adalimumab and etanercept), longer PsA duration (infliximab only), and higher swollen joint count and erythrocyte sedimentation rate (infliximab only).

TNFi persistence

The difference in TNFi persistence between the combination therapy and monotherapy groups of all TNFis was not statistically significant (32.4 vs 30.8 months, p=0.73; table 2 and figure 2A). Considering individual TNFis, TNFi persistence between combination therapy and monotherapy favoured monotherapy with etanercept, was similar between therapies with adalimumab, and favoured combination therapy with infliximab (table 2 and figure 2B–D).

Table 2

TNFi persistence, time to remission and time to loss of remission

Figure 2

Time to change in initial tumour necrosis factor inhibitor (TNFi) therapy (ie, TNFi persistence) among (A) all patients and those taking (B) adalimumab, (C) etanercept and (D) infliximab. Patients were censored at their last follow-up visit if no change in therapy was observed.

A number of predictors of TNFi persistence were identified for patients receiving combination therapy versus monotherapy (table 3). Hispanic ethnicity (not recorded for a substantial number of patients) and longer PsA duration were significant predictors of longer TNFi persistence in the univariate analysis. A history of prior treatment with MTX, history of coronary artery disease, higher BMI, and most measures of baseline disability and disease activity were significant predictors of shorter TNFi persistence in the univariate analysis. Higher mHAQ and CDAI scores were the only significant predictors (shorter TNFi persistence) in the multivariate analysis (table 3).

Table 3

Univariate and multivariate analyses of predictors for TNFi persistence and time to remission

Achievement of remission

There was no statistically significant difference between the combination and monotherapy groups in median time to achieve remission (20.7 vs 25.1 months; p=0.56; table 2 and figure 3). A number of factors were significant predictors of a longer time to achieving remission in the univariate analysis, including female sex, higher BMI, higher baseline disability and disease activity, and history of hypertension or diabetes (table 3). Higher BMI and mHAQ score were the only significant predictors (longer time to achieving remission) in the multivariate analysis (table 3).

Figure 3

Time to remission. Patients were censored if remission was not achieved, a measure of remission was missing, or there was no further follow-up (TNFi, tumour necrosis factor inhibitor).

Discussion

In a propensity score-matched analysis of patients with PsA who were enrolled in a large US-based registry, combination therapy of TNFi and csDMARD was not significantly different compared with TNFi monotherapy in TNFi persistence or time to remission. When the TNFi persistence analysis was stratified by TNFi, TNFi persistence was significantly longer with combination therapy in patients treated with infliximab, whereas it was either the same with both therapeutic approaches (adalimumab) or longer with monotherapy (etanercept) with the other TNFis. There was no difference between monotherapy and combination therapy in time to remission. Further analysis found significant predictors for shorter/longer TNFi persistence and time to achieving remission.

Our results are similar to those of previous analyses of SSATG22 and DANBIO,23 although in the latter, longer TNFi persistence was noted with combination therapy after adjustment for baseline variables. A NOR-DMARD analysis likewise found significantly longer TNFi persistence with combination therapy at years 1 and 2.21 A possible explanation for these differing results between registry studies may be due to differing contributions of individual TNFis, which have different efficacy profiles in PsA with or without a concomitant DMARD. In our analysis, by initial TNFi, we found that TNFi persistence was longer in monotherapy with etanercept, equivalent between therapeutic approaches with adalimumab, and longer in combination therapy with infliximab. These individualised results are mostly consistent with the NOR-DMARD study, which found equivalent TNFi persistence between monotherapy and combination therapy with etanercept and adalimumab and a difference favouring combination therapy with infliximab.21 The results by TNFi in our analysis and those of NOR-DMARD suggest that combination therapy in PsA has little effect on TNFi persistence compared with monotherapy with etanercept and adalimumab, whereas patients taking infliximab may experience significant benefit by also taking a csDMARD such as MTX. Phase 3 studies of infliximab, etanercept and adalimumab for PsA stratified patients according to baseline MTX use in the treatment and placebo groups.10–12 The presence or absence of background MTX had no effect on therapeutic response, but because the enrolled patients were not having adequate response to MTX, these trials do not provide adequate evidence for assessing the potential role of combination MTX in PsA for any of the TNFis. The superiority of combination therapy to monotherapy for TNFis in RA, however, is well established.15 ,16 ,30 ,31 The RA indication for infliximab recommends concomitant MTX use to prevent loss of efficacy.30 The exact mechanism of action for increased efficacy of combination therapy for all TNFis in RA and for infliximab combination therapy in PsA is not clear. It is known that MTX increases infliximab and adalimumab drug levels.30 ,32 Analyses of TNFis in RA have indicated that immunogenicity may play a role in reduced drug levels and efficacy,33–35 and some evidence is emerging that a similar phenomenon may occur in PsA.36–38 A recent study has suggested that antidrug antibodies are inversely associated with red blood cell MTX polyglutamates and that the latter are positively correlated with the trough serum level of infliximab.39 Finally, a recent analysis of two nationwide registries found that infliximab monotherapy was associated with shorter TNFi persistence than combination therapy with MTX.40

We investigated whether there were predictors of longer or shorter TNFi persistence. We found that ethnicity, PsA duration (time since diagnosis), BMI, baseline disability and disease activity, history of coronary artery disease, and history of MTX treatment were significant predictors in this analysis. A longer PsA duration was associated with longer TNFi persistence. This characteristic was not a significant predictor of TNFi persistence in an analysis of SSATG22 and was not significantly different between patients with and without minimal disease activity (MDA) in a prospective, clinical cohort of 146 patients.41 Higher BMI was predictive of shorter TNFi persistence, which is consistent with a report that found obesity to be negatively correlated with achievement of MDA.42 Greater severity of baseline disability and disease activity and a history of prior treatment with MTX were also predictive of shorter TNFi persistence. In NOR-DMARD, patient global assessment was significant in a univariate but not multivariate analysis of biologic survival.21 In DANBIO, an unadjusted analysis found that the number of tender joints and HAQ score were predictive of treatment discontinuation, but none of these were included in the final multivariate model.23 Finally, the SSATG study found that baseline CRP but no other marker of disease activity (or disability) predicted longer TNFi drug survival.22

We also investigated the time to achieving remission and predictors for longer/shorter time to remission. Our analysis showed that the median time to achieve remission was not significantly different between the combination therapy and monotherapy groups (21 vs 25 months). Notably, female sex, higher BMI, and higher baseline disability and disease activity measures were significantly associated with longer time to remission. Previous analyses have found results mostly similar to ours. Female sex was negatively predictive of European League Against Rheumatism (EULAR) good response but not associated with ACR20/50/70 response in the DANBIO study.23 In our study and in previously published analyses, obesity has been negatively associated with achieving and maintaining remission and MDA, respectively.42 Finally, the DANBIO study showed that a number of baseline disability and disease activity measures (HAQ, tender joint counts, swollen joint counts and low levels of CRP) were predictive of poorer clinical response (achievement of ACR improvement criteria or EULAR good response).23

Several limitations to our study need to be mentioned. Observational studies are subject to channelling bias, that is, patients will be directed to therapies based on factors, such as TNFi persistence and remission, that may substantially affect follow-up results. Non-random factors may have led to the choice of biologic as well as the choice of combination therapy or monotherapy. We attempted to control for channelling bias between combination therapy and monotherapy by matching patients using propensity scores. Although propensity score matching reduces the effect of baseline characteristics, the resulting analysis is still not randomised and may be influenced by unmeasured or unanalysed factors. Propensity scoring also reduces the statistical power of an analysis, as shown in a Corrona study comparing csDMARDs and biologics.29 Also, we used the inverse probability treatment weighting method to create weights based on the propensity score. Since we did not use propensity score matching to create two comparable monotherapy and combination therapy groups, we cannot provide a comparison of matched monotherapy vs combination therapy cohorts in order to evaluate how well other covariates are balanced. The inverse probability treatment weighting method allowed us, however, to keep a maximal sample size, which would not have been possible with propensity score matching. Similar to propensity score matching, moreover, the inverse probability treatment weight method removed systematic differences between the monotherapy and combination therapy groups.

We were not able to assess MDA owing to the lack of two of seven required elements (body surface area and enthesitis score)43 in the Corrona records we queried. Therefore, our analysis of remission did not include all the relevant disease characteristics of PsA. Another possible limitation of the remission analyses is that the CDAI has been identified as a conservative estimate of disease activity, particularly in its remission criteria.44 Moreover, neither CDAI nor remission has been validated as disease activity measures in PsA but both are considered reasonable surrogates. Another limitation is that we did not assess reasons for TNFi discontinuation and excluded patients with less than 1 year of follow-up, possibly biasing our study population towards longer TNFi persistence. Finally, our analysis of TNFi persistence by individual TNFi may have been underpowered to avoid a type 2 error. While the results of our analysis are interesting and seem to confirm those of NOR-DMARD,21 prospective studies are needed to assess the comparative effectiveness of monotherapy and combination therapy in PsA.

There is variable evidence regarding the efficacy of combination therapy versus monotherapy in PsA.18 ,20 Our analysis of the Corrona registry provides insight from a US-based clinical practice perspective. Even when adjusting for channelling bias with propensity scores, TNFi persistence and time to remission were not significantly different between the combination therapy and monotherapy groups. TNFi persistence was significantly lower with infliximab monotherapy compared with combination therapy, favoured monotherapy with etanercept, and was not significantly different between therapies with adalimumab. The difference between infliximab combination therapy and monotherapy may be related to background disease severity, immunogenicity or other factors. Greater degree of disability and disease activity use were associated with shorter duration of initial therapy. Prospective, randomised clinical trials evaluating the efficacy of combination therapy versus monotherapy would provide much-needed clarity on treatment options for patients with PsA.

Acknowledgments

Janice Carlson, PhD; Kathryn Boorer, PhD; and Tim Peoples, MA, from Amgen Inc provided medical writing assistance.

References

View Abstract

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors The study was designed and the article drafted by PJM, DHC, KCS and JDG. GL conducted the statistical analysis. All authors contributed to the interpretation of data, revised the draft critically for intellectual content and approved the final version of the article. All authors agree to be accountable for all aspects of the work.

  • Funding This study is sponsored by Corrona, LLC. The Corrona rheumatoid arthritis registry has been supported through contracted subscriptions in the past 2 years by AbbVie, Amgen, AstraZeneca, Genentech, Horizon Pharma, Eli Lilly, Janssen Biotech, Novartis, Pfizer, Vertex, and UCB.

  • Competing interests PJM: Grant, research support, consultant and/or speakers bureau from/for AbbVie, Amgen, Biogen Idec, Bristol Myers Squibb, Genentech, Janssen, Lilly, Pfizer, and UCB; grant, research support, and consultant from/for Celgene, Merck, and Novartis. DHC: employee and shareholder of Amgen Inc. KCS: employee of Corrona, LLC. GL: employee of Axio, LLC (contractor of Corrona, LLC). JMK: consulting for Amgen, AbbVie, and Janssen; employee and shareholder of Corrona, LLC. JDG: employee and shareholder of Corrona, LLC; consultant for AstraZeneca, Novartis, and Pfizer.

  • Patient consent Obtained.

  • Ethics approval New England IRB.

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

  • Data sharing statement No additional data are available.