Abstract
Background:
The aim of this study was to evaluate disease-activity-guided stepwise tapering or discontinuation of rhTNFR:Fc, an etanercept biosimilar, in patients with ankylosing spondylitis (AS) in a prospective, randomized, open-label, multicentric study.
Methods:
Active AS patients with AS disease activity score (ASDAS) ⩾2.1 recruited from 10 hospitals were treated with rhTNFR:Fc 50 mg weekly for 12 weeks, and further randomized into different tapering or discontinuation groups according to ASDAS at week 12. Patients who achieved clinical remission (ASDAS < 1.3) were assigned randomly to stepwise tapering group or discontinuation group. Patients who achieved low disease activity (LDA, 1.3⩽ASDAS < 2.1) were assigned randomly to stepwise tapering, delayed tapering, or discontinuation group. All patients were evaluated every 12 weeks until week 48. The primary endpoint was cumulative flare rates in different groups at week 48.
Results:
A total of 311 patients were enrolled with an average ASDAS of 3.6 ± 1.0, and 259 completed 12 weeks of rhTNFR:Fc induction therapy, with 148 patients (57.1%) achieved clinical remission, 100 (38.6%) achieved LDA, and 11 (4.3%) remained as high disease activity (ASDAS⩾2.1). In patients who achieved clinical remission at week 12, stepwise tapering of rhTNFR:Fc demonstrated significantly lower flare rates at each evaluation compared with discontinuation. In patients who achieved LDA, there was no significant difference of flare rates between stepwise tapering, delayed tapering, and discontinuation. With stepwise tapering of rhTNFR:Fc, flare rates were comparable in AS patients, irrespective of initial ASDAS before tapering.
Conclusion:
Stepwise tapering of rhTNFR:Fc when patients achieved clinical remission was able to maintain favorable low flare rates in 48 weeks. LDA was an alternative therapeutic target, as well as an viable timing for initiation of rhTNFR:Fc tapering. rhTNFR:Fc 25 mg monthly maintained flare-free status in a considerable number of patients. However, abrupt discontinuation of rhTNFR:Fc even if patients achieved clinical remission should be avoided.
Trial registration:
ClinicalTrials.gov: NCT03880968,
Introduction
Ankylosing spondylitis (AS), a subset of axial spondyloarthritis (axSpA), is a chronic inflammatory disorder characterized by inflammatory back pain and predominant involvement of sacroiliac joints and spine, leading to bony fusion of vertebrae and, eventually, disability in some patients. 1 Nonsteroidal anti-inflammatory drugs (NSAIDs) are recognized as a first-line therapy for AS, 2 but the overall response rates to NSAIDs are considerably unsatisfactory. 3 With the advent of biologics, the outcomes of AS patients have been greatly improved. Biologics, including tumor necrosis factor α (TNFα) inhibitors (TNFi), have been included in many recommendations for the treatment of AS.1,2,4,5 Etanercept, a recombinant human TNFα receptor, is capable of binding to TNFα and blocking its biological activities. 6 It is effective in relieving symptoms, improving physical function, and reducing disease activity in patients with AS, and generally no severe adverse effects have been reported.6,7 However, the high expense of biologics restricts their long-term use, which urges a viable strategy to reduce the dosage of biologics while maintaining an optimal therapeutic efficacy.
Several studies have investigated the tapering and discontinuation strategies for biologics in AS.8–20 However, when to start tapering or discontinuing biologics, and what is the optimal tapering strategy in AS, remain undetermined. In some studies, dose reduction was initiated on the basis of the time schedule after a certain period of standard etanercept therapy,8,17 but, in most studies, tapering was considered when remission or low disease activity (LDA) was achieved. Nevertheless, the definitions of remission and LDA have been discrepant across different studies. For instance, some defined remission as Bath ankylosing spondylitis disease activity index (BASDAI) <4,11,18 whereas some defined remission as BASDAI < 2, 13 and in still other studies, LDA was defined as BASDAI < 4.10,14,15,19 Multiple strategies of tapering have also been suggested in different studies, including prolongation of administration intervals, reduction in dosages, or a combination of both. In most studies, the reduced dosage was predetermined and unchanged across the whole observation. Few studies adopted the patient-tailored tapering strategy, and the dosage of etanercept was reduced gradually. 15
To investigate the stepwise tapering and discontinuation of TNFi based on disease activity in patients with AS, a 48-week, prospective, randomized, multicentric study was conducted. An etanercept biosimilar, rhTNFR:Fc (recombinant TNF receptor: Fc fusion protein, Yisaipu), which is one of the most widely used biosimilars in China, was used in this study.
Methods
Study design
This 48-week prospective, randomized, open-label study was conducted in 10 hospitals in southeast China, including the Second Affiliated Hospital of Zhejiang University School of Medicine, Shanghai Guanghua Integrative Medicine Hospital Affiliated to Shanghai University of Tradition Chinese Medicine, the First Affiliated Hospital of Wenzhou Medical University, the First Hospital of Jiaxing, Shaoxing Second Hospital, Shanghai Jiaotong University Affiliated Sixth People’s Hospital, Ningbo Medical Center Lihuili Hospital, Wenzhou Central Hospital, Zhejiang Provincial People’s Hospital, and Shaoxing People’s Hospital. The trial was designed to evaluate disease activity guided tapering and discontinuation strategies of rhTNFR:Fc. Patients with AS were treated with standard dosage of rhTNFR:Fc for 12 weeks as the induction therapy in the first stage of this study. Afterwards, based on whether patients achieved clinical remission or LDA at the end of 12 weeks, they were further assigned randomly to different tapering or discontinuation groups.
The study was approved by local ethics committee (SAHZU2012-13), and was performed in accordance with Good Clinical Practice guidelines of the International Conference on Harmonization and the principles of the Declaration of Helsinki. This trial was registered [ClinicalTrials.gov identifier: NCT03880968].
Patients
Patients with AS were recruited between 1 March 2012 and 30 September 2013. Eligible patients were aged between 18 years old and 65 years old, and were diagnosed with AS according to 1984 revised New York classification criteria. 21 Inclusion criteria included an active disease of AS disease activity score (ASDAS) using C reactive protein (ASDAS-CRP) ⩾2.1, a disease duration of 6 months to 30 years, and no exposure to biologics in the 6 months before recruitment. Concomitant medications including NSAIDs, conventional disease modifying anti-rheumatic drugs (cDMARDs), or prednisone or a prednisone equivalent (⩽10 mg/day) were retained if they had been maintained at a stable dose for at least 4 weeks before recruitment. Patients with late-stage spinal fusion, severe cardiac, hepatic, renal, hematologic or endocrine diseases, multiple sclerosis, current or past malignancies, active or recurrent infections, use of oral antibiotics 2 weeks or intravenous antibiotics 4 weeks before screening, or current or past or potential tuberculosis were excluded. Patients who were pregnant or breastfeeding, or planning to become pregnant were also excluded. Written informed consents were obtained for all patients.
Randomization and masking
Randomization was implemented using a random number table generated by computer, and cards were placed in sealed opaque envelopes (HXW). Randomization lists were kept at each participating center. In patients who achieved clinical remission at the end of induction therapy, they were allocated in a ratio of 3:1 to rhTNFR:Fc stepwise tapering or discontinuation group. In patients who achieved LDA at the end of induction therapy, they were allocated in a ratio of 3:1:1 to rhTNFR:Fc stepwise tapering, delayed tapering, or discontinuation group. Patients and investigators were unblinded once the interventions were assigned. Sample size was calculated based on the risk ratio, and noninferiority tests in the trial arms were performed.
Procedures
In the induction phase of this study (Figure 1), all AS patients were treated with 50 mg rhTNFR:Fc subcutaneously each week for the first 12 weeks. ASDAS (refers to ASDAS-CRP below unless specified otherwise) was then calculated at the end of week 12. Patients who achieved clinical remission with ASDAS < 1.3 were designated as group A, and those achieved LDA with 1.3⩽ASDAS < 2.1 as group B (previously 1.3⩽ASDAS < 2.1 was recognized as moderate disease activity, but recently LDA was used to indicate status within this range).22,23

Study design.
In the second phase of this study, patients in group A with clinical remission were randomly assigned to stepwise tapering group A1 (treated with rhTNFR:Fc 25 mg weekly for 12 weeks, 25 mg every 2 weeks for 12 weeks, and 25 mg monthly for the last 12 weeks), or discontinuation group A2 (rhTNFR:Fc was discontinued whereas previous concomitant NSAIDs and/or cDMARDs were retained). Patients in group B with LDA were assigned randomly to stepwise tapering group B1 (same tapering strategy as group A1), delayed tapering group B2 (treated with rhTNFR:Fc 50 mg weekly for another 12 weeks, 25 mg weekly for 12 weeks, and 25 mg every 2 weeks for 12 weeks), and discontinuation group B3 (same discontinuation strategy as group A2) (Figure 1).
Outcomes
The primary endpoint was cumulative flare rates based on ASDAS at week 48 with different tapering or discontinuation strategies. All patients were evaluated at week 0, 12, 24, 36, and 48 for symptoms, physical signs, cotherapies as well as adverse events. ASDAS, BASDAI, Bath AS metrology index (BASMI), and Bath AS functional index (BASFI) were documented, and complete blood count, chemistry panel, CRP, erythrocyte sedimentation rate (ESR), and urinalysis were assayed at each visit. Antinuclear antibody (ANA) and chest radiographs were examined at week 0 and 48. After week 12, flares were defined as ASDAS > 2.1, 24 and flared patients restarted rhTNFR:Fc 50 mg weekly with subsequent tapering.
Statistical analysis
Differences in treatment strategies were evaluated by unpaired two-sided t test, Wilcox test, ANOVA, or Kruskall–Wallis test (p < 0.05). Kaplan–Meier analysis was used to estimate the cumulative probability of the relapse rates. Survival curves were generated for different rhTNFR:Fc tapering or discontinuation strategies from week 12 to endpoint. χ 2 test was used to compare cumulative relapse rates, and log-rank test was applied to compare survival curves. Cox proportional hazard model was used to determine baseline predictors that were associated with risks of disease flares. Crude hazard ratios (HRs) were calculated for each candidate characteristic, whereas adjusted HRs were calculated by including redefined statistically significant characteristics (p < 0.10) in one model simultaneously, or adjusted for statistically significant variables at baseline when comparing relapse rates in a disease activity guided manner. All analyses were done using R Statistical Software, version 3.5.0 (R Foundation for Statistical Computing, Vienna, Austria).
Results
Patients characteristics
A total of 311 patients with active AS were recruited, with 223 men and 88 women and an average age of 31.0 ± 11.3 years old. ASDAS was 3.6 ± 1.0 (range 2.1–6.4) at inclusion. All patients were managed with rhTNFR:Fc 50mg subcutaneously each week, and 259 patients completed induction therapy with rhTNFR:Fc for 12 weeks. At the end of week 12, 148 (57.1%) patients achieved clinical remission (ASDAS < 1.3) with the absolute change of ASDAS from week 0 (|∆ASDAS|) ranging from 0.9 to 4.5, 100 (38.6%) achieved LDA (1.3⩽ASDAS < 2.1) with |∆ASDAS| from 0.1 to 3.9, and 11 (4.3%) remained high disease activity (ASDAS > 2.1). Subsequently, 142 patients with clinical remission and 93 patients with LDA were randomized to different tapering or discontinuation groups. Patients with high disease activity despite rhTNFR:Fc for 12 weeks withdrew from the study and switched to other medications (Figure 2). The characteristics of patients in subgroups were comparable at the end of week 12 before they were further randomized with different therapeutic strategies (Table 1).

Study profile.
At the end of week 48, a total of 183 AS patients completed the study, with 33 patients flared and 150 remaining flare-free. From week 12 to week 48, 52 patients withdrew from the study, mostly because of loss to follow up. rhTNFR:Fc was generally safe without severe adverse events. Only five mild adverse events were reported, including two cases with leukopenia, two cases with abnormal liver function, and one case with injection site reaction; all recovered shortly after discontinuing rhTNFR:Fc.
Flare rates with different tapering or discontinuation strategies
The cumulative flare rates were calculated every 12 weeks until the end of week 48. In patients who achieved clinical remission at week 12 with subsequent stepwise tapering of rhTNFR:Fc (A1), the flare rate increased over time when rhTNFR:Fc was gradually tapered every 12 weeks, from 1.0% at week 24, to 4.3% at week 36, and 9.0% at week 48. When compared with patients who discontinued rhTNFR:Fc (A2), patients in the A1 groups with stepwise tapering had significantly lower flare rates at each time point of evaluation (supplemental Table S1). The average time to flare from the initiation of tapering or discontinuation was 28.5 ± 9.0 weeks in A1, and 18.86 ± 9.44 weeks in A2 (p = 0.06) (Figure 3). In patients who achieved LDA at week 12 (B1, B2, B3), the cumulative flare rates all increased over time, but there was no significant difference among three groups at each time point (supplemental Table S1). No significant difference of average time to flare was identified among B subgroups (Figure 3). Likewise, when survival curves were generated using Kaplan–Meier approach and evaluated by log-rank test, significant difference between A1 and A2 was noticed, whereas no significance among B subgroups was indicated (Figure 4).
Demographics, disease characteristics and concomitant medications at the end of week 12.
Median (first quartile, third quartile).
Normal distributed data was demonstrated as mean ± SD.
cDMARDs, conventional disease modifying anti-rheumatic drugs; CRP, C reactive protein; ESR, erythrocyte sedimentation rate; NSAIDs, nonsteroidal anti-inflammatory drugs; SD, standard deviation.

Average time to flare with different tapering or discontinuation strategies.

Flare-free survival rates.
For patients who shared the same tapering strategy, those who achieved clinical remission before reducing the dosage (A1) demonstrated lower flare rate as well as prolonged average time to flare than those who achieved LDA when tapering was initiated (B1) (supplemental Table S1, Figure 3). However, when patients with same strategies were compared using Cox regression analysis after adjusting baseline ASDAS and ESR, there was no significant difference of flare rates between A1 and B1, or between A2 and B3 (Table 2).
Disease activity guided comparison of flare rates between groups with the same tapering strategy.
Adjusted by ASDAS and ESR level.
ASDAS, ankylosing spondylitis disease activity score; CI, confidence interval; ESR, erythrocyte sedimentation rate; HR hazard ratio.
Risk factors of disease flares during rhTNFR:Fc tapering or discontinuation
To identify the risk factors of flares during tapering or discontinuation, Cox regression analysis was applied to evaluate the parameters before reduction or discontinuation of rhTNFR:Fc. In patients who achieved clinical remission (A), discontinuation of rhTNFR:Fc, high levels of CRP, and severe stage of sacroiliitis were risk factors for future relapse. However, no risk factors were identified in group B (Table 3).
Risk factors of flare before rhTNFR:Fc tapering or discontinuation.
Crude HRs were calculated in Cox regression models that included only one independent variable at a time. Adjusted HRs were calculated in a multivariable Cox regression model that included all redefined statistically significant variables (p < 0.10) at once.
ASDAS-C, ankylosing spondylitis disease activity score using CRP; CI, confidence interval; CRP, C reactive protein; ESR, erythrocyte sedimentation rate; HR hazard ratio.
Sensitivity analysis
Considering that patients withdrew in the rhTNFR:Fc discontinuation groups due mainly to poorer compliance, sensitivity analysis was performed to evaluate the robustness of our results. By assuming that all the patients who withdrew from the study experienced a disease flare, sensitivity analysis demonstrated roughly similar results (supplemental Tables S2 and S3).
Discussion
Several studies have investigated the strategies of dosage tapering or discontinuation for biologics in AS patients. However, no extensively accepted agreement has been established, and the optimal timing to start tapering as well as favorable tapering strategies are yet to be investigated. At least 13 prospective or retrospective studies have focused on the tapering or discontinuation of etanercept in AS, with follow-up duration varying from 8 weeks to 36 months and the maintenance rate of remission or LDA as 47–86.3% (Table 4).8–20 According to a systemic review and a meta-analysis on etanercept tapering in AS patients, etanercept reduction was effectively comparable with nondecreasing standard dosages in most studies.25,26 In our study, the flare-free maintenance rates with stepwise tapering of rhTNFR:Fc at week 48 were 91.0% and 83.3% in patients with clinical remission and LDA, respectively. When to start tapering etanercept in AS patients has been inconsistent among recommendations or studies. According to the 2016 updates of management recommendations for axSpA, tapering of biologics could be considered if the patient was in sustained remission, which was defined as inactive disease based on ASDAS for 6 months or longer. 2 But in the 2018 update of French society for rheumatology (SFR) recommendations, other than patients with disease remission for at least 6 months, those with low level of activity were also qualified for dose reduction. 27 To our knowledge, our study is the first to investigate the tapering timing in the same research by dichotomizing patients into subgroups according to disease activity. Tapering at clinical remission or LDA were both evaluated. When the same tapering strategy was applied to patients with different disease status, patients with initial clinical remission demonstrated lower cumulative flare rates and prolonged time to relapse when compared with those with initial LDA. After adjustment with ASDAS and ESR, there was no significant difference between these two groups of patients, suggesting that LDA can also be considered as an appropriate timing for the initiation of rhTNFR:Fc dosage reduction. Besides, it further confirms that low/minimal disease activity may be an alternative treatment target. 28
Studies on ETN tapering or discontinuation in AS patients.
AS, ankylosing spondylitis; ASAS, assessment in ankylosing spondylitis; BASDAI, Bath AS disease activity index; BASFI, Bath AS functional index; CRP, C reactive protein; ESR, erythrocyte sedimentation rate; ETN, etanercept; LDA, low disease activity; HAQ, health assessment questionnaire; VAS, visual analogue scale.
Other than the timing of, or the disease activity status before, tapering, the specific procedures of reducing etanercept dosage was critically essential in clinical practice. The maximal interval of etanercept injection in the above mentioned 13 studies was 3 weeks at a dosage of 25 mg, 14 while, in others, 25 mg weekly or every 2 weeks was applied. The minimal dosage in our study, 25 mg each month, has not been documented in any previous research. When rhTNFR:Fc was reduced from 25 mg every 2 weeks to 25 mg each month for 12 weeks, the cumulative flare rate increased from 4.3% to 9.0% in patients who were in clinical remission before tapering, and from 14.0% to 16.7% in those with LDA before tapering. Whether rhTNFR:Fc 25 mg each month can sustainably maintain LDA or clinical remission of longer duration should be investigated in further studies.
One previous study evaluated the discontinuation of etanercept in AS patients who achieved ASAS20, and the cumulative probabilities of relapse at 1, 2, and 3 years were 45.7, 57.1, and 60.0%, respectively. 20 In our study, when rhTNFR:Fc was discontinued, the flare rate was 31.8% after 36 weeks in patients with clinical remission, and 42.9% in patients with LDA. Irrespective of the disease status patients had achieved before rhTNFR:Fc tapering or discontinuation, if the rhTNFR:Fc was discontinued, the relapse rates increased over time, even if the NSAIDs and/or cDMARDs were retained. Discontinuation of rhTNFR:Fc was identified as one of the risk factors of disease flare in patients with clinical remission, indicating that abrupt discontinuation of rhTNFR:Fc should be avoided. Longer duration of remission before tapering or gradual dosage reduction before discontinuation may be investigated.
Multiple parameters have been developed to evaluate the disease activity of AS, including the BASDAI, BASFI, BASMI, ASAS, and ASDAS. Most of the above mentioned studies on etanercept tapering used BASDAI or ASAS to define disease activity status, and only one study also included ASDAS as an assessment. 15 According to the 2017 update of recommendations on treating spondyloarthritis to target, ASDAS is a preferred measure in axSpA as it is more associated with various biomarkers of inflammation than BASDAI and correlates better with magnetic resonance imaging (MRI) changes. 28 ASDAS has also been considered as highly discriminatory in the detection of therapeutic differences and improvement from baseline. 29 Our study used ASDAS, and it was correlated highly with BASDAI (data not shown). We did not measure the biomarkers of inflammation or MRI changes, so whether ASDAS outweighs BASDAI on evaluations during rhTNFR:Fc tapering needs further investigation.
The cutoff values of ASDAS were updated in 2018 by ASAS. ASDAS < 1.3 remained as inactive, and 1.3⩽ASDAS < 2.1, previously indicating moderate disease activity, was considered as LDA. 23 Flare was defined as ASDAS⩾2.1 in our study, 24 but we also evaluated a compound definition of flares, with ASDAS⩾2.1 for patients with LDA, and ASDAS⩾1.3 as well as ∆ASDAS⩾0.6 for patients with clinical remission. It was more complicated, seemingly more reasonable, and tailored to individual subgroups when compared with the single criteria of ASDAS⩾2.1 in all patients, yet this combined definition for flares did not change the results (data not shown). Another definition of clinical important worsening in axSpA was an increase in ASDAS of at least 0.9 points, 30 which was also applied to evaluate the flares in our patients, and the conclusions remained unchanged (data not shown).
It is noteworthy that, at the end of induction therapy for 12 weeks, only 4.3% patients remained with high disease activity, whereas the majority achieved either remission or LDA. A possible explanation is that 50 patients withdrew consent before completion of induction therapy at week 12 (Figure 2), with some complaining of inconvenience, lack of efficacy, and other reasons. Presumably, some of them should remain in high disease activity status at week 12 even if they continued with the induction therapy. As they did not return at week 12 for evaluation, the actual disease activity they had was unknown. The proportions we observed were based on the number of patients who completed 12-week induction therapy.
There are some limitations of our study. This trial was completed by the end of 2014, and it took a long time to publish the results, due mainly to the time-consuming data extraction from all case report forms, the unexpected resignation of the initial statistician for personal reasons, and multiple revisions before submission. Of the 259 patients who finished the initial 12-week rhTNFR:Fc induction, only 183 (70.7%) completed the study at week 48. The compliance of patients should be more carefully monitored to improve the follow-up rate. Besides, the flare rates in the delayed tapering group (B2) were higher than those in stepwise tapering group (B1), even though there was no significant difference between these two groups on the flare rates at each evaluation. We currently do not have a reasonable explanation for this phenomenon, and more patients in group B2 should probably be included for further observation. The tapering of rhTNFR:Fc was observed for only 36 weeks, and whether longer duration of low-dose rhTNFR:Fc was effective warrants future investigation. Lastly, comparisons between groups would also be clearer and more direct if we had included an arm with standard dose of rhTNFR:Fc throughout the trial. However, the current arms in this study represent several commonly used methods of tapering and discontinuation in patients who have to reduce or discontinue rhTNFR:Fc for various reasons, and our results provide evidence of optimal strategies in these clinical scenarios.
To conclude, our study suggested that disease-activity-guided stepwise tapering strategy of rhTNFR:Fc in AS patients is applicable. rhTNFR:Fc stepwise tapering when patients achieved clinical remission was able to maintain favorable flare-free rates across week 48. LDA is an alternative therapeutic target, as well as an appropriate timing for initiation of rhTNFR:Fc dosage reduction. rhTNFR:Fc 25mg each month was capable of maintaining flare-free status in a considerable number of patients and warrants future investigation. Abrupt discontinuation of rhTNFR:Fc was associated with high flare rates and should be avoided regardless of the disease activity status that patients have reached.
Supplemental Material
Supplementary_data_2020-5 – Supplemental material for Disease activity guided stepwise tapering or discontinuation of rhTNFR:Fc, an etanercept biosimilar, in patients with ankylosing spondylitis: a prospective, randomized, open-label, multicentric study
Supplemental material, Supplementary_data_2020-5 for Disease activity guided stepwise tapering or discontinuation of rhTNFR:Fc, an etanercept biosimilar, in patients with ankylosing spondylitis: a prospective, randomized, open-label, multicentric study by Ting Zhang, Jianing Zhu, Dongyi He, Xiaowei Chen, Hongzhi Wang, Ying Zhang, Qin Xue, Weili Liu, Guangbo Xiang, Yasong Li, Zhongming Yu and Huaxiang Wu in Therapeutic Advances in Musculoskeletal Disease
Footnotes
Acknowledgements
We thank all participants and employees of participating centers for their time and effort.
Contributors
TZ, DYH, XWC, HZW, YZ, QX, WLL, GBX, YSL, ZMY, and HXW contributed to study design and data collection. TZ and JNZ performed data analyses. All authors contributed to data interpretation and writing of the manuscript. All authors approved the final manuscript for submission.
Data sharing
Individual deidentified participants’ data that underlie the results reported, the study protocol, statistical analysis plan, and analytical code are available with publication to researchers who provide a methodologically sound proposal to achieve aims in the approved proposal or for individual participant data meta-analysis. Proposals should be directed to wuhx8855@zju.edu.cn; to gain access, data requestors will need to sign a data access agreement.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by 3SBIO Inc., which provided some rhTNFR:Fc used in this study. The funders had no role in study design, data collection, data analysis, data interpretation, or writing of the manuscript. The corresponding author had full access to all data in the study and had final responsibility for the decision to submit for publication.
Conflict of interest statement
The authors declare that there is no conflict of interest.
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
