Study TPX-100-5: Intra-Articular TPX-100 Significantly Delays Pathological Bone Shape Change and Stabilizes Cartilage in Moderate to Severe Bilateral Knee OA

doi:10.21203/rs.3.rs-259770/v1

Download PDF

Research article

Study TPX-100-5: Intra-Articular TPX-100 Significantly Delays Pathological Bone Shape Change and Stabilizes Cartilage in Moderate to Severe Bilateral Knee OA

https://doi.org/10.21203/rs.3.rs-259770/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background: TPX-100, a promotor of osteoblast and chondroblast differentiation, is a potential osteoarthritis (OA) therapy. This retrospective study compared MRI 3D femoral bone shape changes (B-scores) after intra-articular TPX-100 or placebo and analyzed the relationship between cartilage thickness and bone shape change over 12 months.

Methods: 104 participants with bilateral knee moderate to severe (ICRS 2-4) knee cartilage defects were randomized for evaluation of efficacy and safety of 200mg of TPX-100. Each subject’s contralateral placebo-treated knee served as a paired internal control. After MRI quality control, 78/93 subjects (84%; 156 knees) were analyzed for quantitative femoral B-score and cartilage thickness. All analyses were centrally performed, blind to treatment assignment and clinical data.

Results: TPX-100-treated knees (n=78) demonstrated a statistically significant decrease in pathologic bone shape change compared with placebo-treated knees at 6 and 12 months: 0.0298 (95% C.I. -0.037, 0.097) vs 0.1246 (95% C.I. 0.067, 0.182) (P=0.02); and 0.0856 (95% C.I. 0.013, 0.158) vs. 0.1969 (95% C.I. 0.123, 0.271) (P = 0.01), respectively. The correlation between bone shape change and medial and total tibiofemoral cartilage thickness changes at 12 months was statistically significant in TPX-100-treated knees (P<0.01).

Conclusions: This is the first report of a potential therapy demonstrating a significant effect on bone shape measured by B-score in knee OA. These data, in combination with previously reported statistically significant and clinically meaningful improvements in WOMAC function versus placebo, support TPX-100 as a candidate for disease modification in knee OA.

Trial Registration: NIH clinicaltrials.gov, NCT01925261. Registered 15 August 2013, https://clinicaltrials.gov/ct2/show/NCT01925261?term=NCT01925261

Health Economics & Outcomes Research

Osteoarthritis

DMOAD

TPX-100

Bone shape

Machine learning

TPX-100 is a 23-amino acid peptide derived from Matrix Extracellular Phosphoglycoprotein (MEPE), a small integrin-binding ligand N-linked glycoprotein (SIBLING) family member. MEPE is highly expressed by osteocytes, is downregulated in osteoarthritis, and may play a role in osteoarthritic bone remodeling [1]. TPX-100 acts selectively on cells committed to hard tissue lineage, viz. bone, cartilage or dentin. In animal models, the actions of TPX-100 are limited to bone and cartilage formation and repair [2]. Based on preliminary in vitro data, the putative mechanism of action of TPX100 is via integrin binding with formation of TGF-β receptor complexes and activation of its effector, Smad3. Smad3 can regulate cartilage matrix synthesis via induction of type II collagen in chondrocyte lineage cells and play a role in homeostasis of knee joint bone [3][4].

TPX-100 has been shown to induce articular cartilage formation in goats (N = 8 / dose group) after a standardized full-thickness chondral defect and treatment with 4 weekly intra-articular (IA) injections of TPX-100 (25, 125 or 250 mg/injection) vs. vehicle in a dose-dependent manner. After 6 months, histopathological staining in TPX-100-treated joints demonstrated robust articular (hyaline) cartilage formation and increased Type II collagen compared with vehicle-treated controls [2]. In a standardized rat model of ACL transection and partial medial meniscectomy, IA TPX-100 dose dependently reduced joint damage and improved osteoarthritis scores in actively treated animals compared with controls.

In humans, a Phase II randomized, double-blind, placebo-controlled, 12-month trial of TPX-100 (NCT01925261) has been completed to evaluate safety, tolerability and efficacy of IA TPX-100 in subjects with bilateral, mild-moderate (ICRS grades 2–3) patellofemoral cartilage defects, with or without tibiofemoral cartilage defects [2][5]. Each subject’s contralateral placebo-treated knee served as a paired internal control, intended to control for effects of age, sex, weight, genetic factors, and activity levels on outcome measures. The pre-selected primary efficacy outcome measure in this trial was the 6-month change in patellar cartilage thickness as measured using standardized magnetic resonance imaging (MRI) in TPX-100-treated knees compared with placebo-exposed knees (See Fig. 1; CONSORT diagram). Subjects were screened clinically and with MRI. Synovitis and meniscal damage were among the exclusion criteria for enrollment. Briefly, Part A of the study evaluated safety of four, once-weekly IA doses of TPX-100 in sequential dose cohorts (25, 50, 100 and 200 mg/injection) of 6–9 subjects, with progression to the next dose following Safety Review Committee approval. All doses were reasonably safe and well tolerated, and the highest dose, 200 mg/injection, was selected for Part B of the study for evaluation of efficacy and safety. Of the 118 subjects enrolled for Parts A and B, 93 subjects received 4 injections of 200 mg TPX-100 and had baseline MRIs with at least one follow-up scan. Per the statistical analysis plan, these subjects made up the primary analysis population. Efficacy outcome measures included MRI cartilage measures (6 and 12 months) and patient-reported outcomes (WOMAC, KOOS, NRS for pain at 3, 6 and 12 months). Semi-quantitative MRI analysis by central readers blind to clinical data and treatment assignment demonstrated that only 14% of knees had measurable patellar cartilage thickness changes, with no significant treatment differences at 6 or 12 months. MRIbased MOAKS evaluation, including cartilage defects, meniscal pathology and Hoffa's synovitis, did not show differences at baseline or follow up between TPX-100-treated knees and placebo-exposed knees, and there were no significant within-knee changes in bone marrow lesions at 6 or 12 months. In contrast, patient benefit, measured by WOMAC and KOOS scores, was statistically significant and clinically meaningful in favor of TPX-100-treated knees at 6 months compared with placebo, with robust functional benefits sustained through the end of the study at 12 months [2]. Posthoc analyses revealed that 68 (73%) of subjects had, in addition to bilateral patellofemoral cartilage defects, moderate to severe (ICRS 2–4) bilateral tibiofemoral cartilage defects. In these subjects, sustained, statistically significant and clinically meaningful clinical benefits in favor of TPX-100-treated knees were observed, nearly identical to those in the whole population [5].

The present study (TPX-100-5) was designed as a retrospective MRI study to investigate femoral boneshape change at 6 and 12 months after TPX100 or placebo administration and to analyze relationships between cartilage thickness and femoral bone shape change at 6 and 12 months after TPX100 or placebo administration.

Bone shape change, measured by MRI, has been shown to predict radiographic onset of OA [6], is associated with radiographic structural progression [7], discriminates people with knee OA from those without knee OA [8] and is more responsive to change over time than is radiographic assessment [9]. In each of these studies, the femur (defined as the whole of the lateral and medial femoral condyles) had greater discrimination and responsiveness to change than did the tibia or patella. The femoral bone shape (“B-score”) metric is a form of statistical z-score that represents the position of a femoral bone shape along a shape vector from a non-OA knee shape (origin) towards an OA knee shape (positive direction). Non-OA and OA knees used to define this 3-dimensional (3D) shape vector were categorized using centrally-read and adjudicated Kellgren-Lawrence grading [10].

These data provide construct validity for this new, continuous scalar measurement. In addition, Bowes et al showed that bone shape is directly associated with the risk of clinical outcome measures such as knee pain, functional deficit and joint failure, as indicated by total joint replacement, with only small effect sizes from adjusting models for potential covariates such as age, sex, ethnicity, body mass index (BMI), alignment, previous knee surgery, non-steroidal anti-inflammatory drugs (NSAIDs) use and smoking status [10]. These findings support femoral B-score as a structural endpoint in clinical trials of disease-modifying osteoarthritis drugs (DMOADs).

Subjects between the ages of 25–75 years with bilateral knee OA were enrolled at 15 sites across the United States.

All subjects were screened by MRI, centrally read, to confirm structural inclusion criteria: bilateral PF cartilage defects (ICRS Grade 1–3) with intact menisci. One knee was randomly assigned to receive 4 weekly injections of TPX100, while the contralateral knee received placebo (saline) injections that were identical in appearance and viscosity. Investigators, subjects, sites and sponsor were blinded to treatment assignment.

MR images were acquired at baseline, 6 and 12 months with 1.5T clinical MRI scanners using a sagittal T1-weighted 3D SPGR, FLASH or FFE sequence with fat saturation or water excitation. An identical scanner and knee coil were used for baseline and follow-up measurements of each participant. Acquisition parameters were as follows: contiguous slice thickness 1.5 mm, in-plane resolution of 0.31 mm, repetition time: 17 msec, echo time: 7 msec, flip angle: 15°. Identical acquisition parameters were used at baseline and follow-up.

In Study TPX-100-5, MRIs from all subjects (n = 93) who received, as randomly assigned, four weekly injections of TPX-100 (200 mg/injection) in one knee (Index) and placebo in the contralateral knee (Control) were eligible for inclusion. Image quality was assessed centrally and blind to clinical data and treatment assignment. Quality was sufficient to include bilateral knee images from 78 of the 93 subjects (84%) analyzed in study TPX100-1 (See Fig. 1; CONSORT diagram).

Cartilage measurement was performed centrally, blind to treatment assignment and clinical data, at a single center (Chondrometrics GmbH, Ainring, Germany). The subchondral bone and cartilage surface areas of the medial and lateral tibia and medial and lateral weight-bearing central femoral condyles were traced manually, excluding cartilage cover of osteophytes [11]. Quantitative measures of cartilage, including mean cartilage thickness averaged over the total area of subchondral bone, were computed using Chondrometrics software [12]. Mean cartilage thickness over the total area of subchondral bone for the medial tibiofemoral compartment was computed by summing the values of the medial tibia and those of the weight-bearing medial femoral condyle; and for the lateral tibiofemoral compartment by summing the values of the lateral tibia and those of the lateral weight-bearing femoral condyle. Definitions of the tibiofemoral subregions as well as measurement reliability (test–retest reproducibility with repositioning of the joint between acquisitions) have been published previously [12].

Femur bone surfaces were automatically segmented using active appearance models (AAMs) at a single center (Imorphics Ltd, Manchester, UK). An AAM is a type of statistical shape model trained using machine learning to search images. AAMs have proven to be a successful supervised machine learning method that can produce a segmented knee bone surface with sub-millimeter accuracy [8][13]. AAMs were constructed using a training set of 96 knee DESS sequence MRIs that were selected to provide examples of OA across KLG grades from subjects enrolled in the OAI [14]. The performance of the AAM has been tested with various other MRI sequences and slice thicknesses. Mean automated segmentation accuracy using MRI FLASH sequences is -0.0009 mm, with ± 95th percentiles of error of + 0.34/-0.43 mm, where a positive error represents the automated surface being outside the reference surface [15]. The construction of an AAM produces a “shape space,” spanned by the set of principal components used to describe the training set of examples. Within this shape space, an “OA vector” is defined as the line which passes through the mean shape of two populations: a population with OA (OA Group, defined as all knees with KLG ≥ 2 over 4 years of follow-up) and a population without OA (Non-OA Group, defined as knees with KLG of 0 in the same period). Distances along the femur OA vector are termed ‘B-score’, with the origin (B-score = 0) defined as the mean shape of the Non-OA Group for each sex, and 1 unit defined as 1 standard deviation of the Non-OA Group along the OA vector (positive values toward the OA Group). Each parameterized femur bone shape was projected orthogonally onto the OA vector to specify the corresponding B-score value [10]. Representative examples of differences in femur bone shape at various B-scores, and a heat map of the areas which change most with increasing B-score are shown in Fig. 2.

B-score changes over 6 and 12 months were graphically compared with the 12-month trajectories of “non-progressor” and “progressor” B-score groups from OAI data. The nonprogressor group was defined as all of those who do not change by more than 95% of the smallest detectable difference (SDD) using the slope of change measured over 4 years. The average slope over this 4-year period is about 0.04 B-scores per annum. The progressor group was defined as all of those who do change by more than 95% of the SDD using the slope of change measured over slope over 4 years. The average slope over this 4-year period is about 0.24 B-scores per annum. SDD was calculated as the 95% limit of agreement between the first and second image measurements from test-retest data, using the Bland-Altman method [10].

Paired Students’ T-tests were used to compare B-score change from baseline between the knee receiving IA placebo (control) or TPX-100 (index) at 6-months and 12-months. For each knee (index, control) at 6months and 12-months, the Pearson coefficients of the changes from baseline between the femur B score and cartilage thickness of tibiofemoral variables were estimated. The two-sided p-value and 95% confidence intervals from the test of the null hypothesis (that the true correlation coefficient is equal to zero) were also computed. Statistical significance was set at P < 0.05. The data analysis for this paper was generated using SAS software version 9.

Characteristics of the per-protocol (TPX-100-1) and B-score analysis (TPX-100-5) cohorts were as follows:

All Per Protocol (PPT) subjects: TPX-100-1, N = 93

Mean Age: 58.4 (95% C.I.: 56.4, 60.4)

Sex: 38 Males, 55 Females (59.1% Females)

Mean BMI: 30.4 (95% C.I.: 29.1, 31.7)
Evaluable B-Score Analyzed Cohort: TPX-100-5, N = 78

Mean Age: 58.4 (95% C.I.: 56.2, 60.6)

Sex: 30 Males, 48 Females (61.5% Females)

Mean BMI: 30.9 (95% C.I.: 29.5, 32.3)

Relative proportions of ICRS grades of Index and Control knees of the subjects analyzed in TPX-100-1 and TPX-100-5 studies are provided in Table 1.

Table 1. ICRS Grades of Index and Control Knees of the Subjects Analysed in TPX-100-1 and TPX-100-5 Studies, Respectively at Baseline. ICRS grade in each knee was determined as the maximum degree of cartilage defect observed in tibiofemoral knee compartment.

Comparison of the 78 evaluable TPX-100-treated knees demonstrated a statistically significant decrease in B-score change in the femur compared with placebo-exposed knees at 6 months: 0.0298 (95% C.I. -0.037, 0.097) vs 0.1246 (95% C.I. 0.067, 0.182) (P = 0.02) and at 12 months: 0.0856 (95% C.I. 0.013, 0.158) vs 0.1969 (95% C.I. 0.123, 0.271) (P = 0.01). Graphical comparisons of Bscore change with progressor and non-progressor groups from OAI data showed B-score trajectories in TPX-100-treated knees that were strikingly similar to those of OA non-progressors, and, among placebo-exposed knees, B-score trajectories that were similar to OA progressors (Fig. 3).

Analysis of quartiles of index knee B-score change at 12 months showed a statistically significant association between cartilage thickness loss, particularly in the lower medial femoral condyle, and lower positive B-score change (Fig. 4).

In TPX-100-treated knees, correlation analyses revealed a statistically significant association between B-score change and medial femoral condyle cartilage thickness change at 6 months; and a statistically significant association between B-score change and cartilage thickness change in entire femoral condyle and entire TF compartments at 12 months (P < .05; Table 2). There was also a significant association between B-score change and cartilage thickness change in the medial femoral condyle and medial TF compartments at 12 months (P < 0.005; Table 2).

Table 2: Pearson Correlations between Femur Bone Shape Stabilization and Tibiofemoral Cartilage Thickening/Stabilization in TPX-100-treated Knee.

Cartilage thickness change in the lateral femoral condyle and lateral TF compartments did not show an association with B-score changes at either 6 or 12 months. There was little or no association between B-score change and cartilage thickness change in control knees (Table 3).

Table 3

Pearson Correlations between Femur Bone Shape Stabilization and Tibiofemoral Cartilage Thickening/Stabilization in Control Knee.
Change Period	Knee Compartment	n	Correlation Coefficient	p-value
Baseline to 6 months	Entire Femoral Condyle	78	-0.068	0.55
	Medial Femoral Condyle	78	-0.053	0.65
	Entire TF	78	0.043	0.71
	Medial TF	78	0.011	0.93
Baseline to 12 months	Entire Femoral Condyle	78	-0.174	0.13
	Medial Femoral Condyle	78	-0.246	0.03
	Entire TF	78	0.011	0.92
	Medial TF	78	-0.085	0.46

In this study, we measured femoral bone shape using the B-score metric that describes changes along a vector tracing a linear path from a normal, non-osteoarthritic state to the pathological shape of the osteoarthritic knee. Pathological shape change is characterized by the broadening and flattening of the femoral condyles and tibial plateaus and the growth of osteophytes. The vector is defined by KLG=0 at its origin and KLG≥2, where a higher B-score is more radiographically pathologic, and, as demonstrated in previous work, is associated with progressive knee OA. Knees treated with IA TPX-100 had a bone shape change at 6 and 12 months similar to that of knees classified as KLG=0. This was substantially different from the trajectory of control knees, which was similar to that of knees classified as KLG≥2. There were no significant baseline group differences in any clinical or MRI measures.

These data represent, to our knowledge, the first report of an investigational agent demonstrating a significant treatment difference compared with placebo on pathologic bone shape (B-score) change in the knee. A previous study of the novel Cathepsin K Inhibitor MIV-711 demonstrated a reduction of pathological bone surface area change in subjects with knee OA at 6 months, with no group differences in clinical outcomes [16]. While clinically meaningful, treatment-related differences in B score changes have not been defined, the difference in the rate of B-score change between TPX-100-treated knees and placebo‑exposed knees were notable, suggesting a slowing of the OA-related bone pathology associated with OA progression and predictive of joint failure. In addition, treatment-associated reductions in B-score correlated significantly with reduced cartilage thickness loss in medial and total tibiofemoral cartilage in TPX-100-treated knees.

There were limitations in this study. TPX-100-5 was a retrospective analysis of the subset of Study TPX-100-1 subjects including those with tibiofemoral OA, mild to severe, in addition to PFOA, whose MRI images were qualified for bone shape analyses. However, the statistical analysis plan for the study was finalized prior to image qualification and analysis, and all analyses were performed blind to timepoint, treatment assignment and all clinical information. As bone-shape analysis was not part of the original (TPX-100-1) study plan, the MR image sequence used was designed for manual cartilage segmentation and therefore was not optimal for automated bone segmentation and shape analysis, for which thinner, 0.75mm slices would have been preferred. Nevertheless, all the image data that could be utilized were included in the analysis without any specific subject selection. The severity of OA in the subject population was defined by using MRI-based cartilage defect (ICRS) scoring after clinical inclusion/exclusion criteria were met. ICRS grading was performed by central readers blind to clinical information in order to establish a baseline of structural severity, and to exclude significant meniscal pathology and synovitis, per study exclusion criteria.

These results of structural imaging in TPX-100-treated knees compared with placebo-exposed knees, in combination with previously reported robust improvements in established clinical outcomes (as measured by WOMAC and KOOS scores), support the potential of TPX-100 as a candidate disease modifying drug in knee OA.

Ethics approval and consent to participate

Written informed consent was obtained prior to study enrollment. Institutional Review Board approval was obtained from Western IRB with approval ID #20131394.

Consent for publication

Not applicable.

Availability of data and materials

Osteoarthritis Initiative (OAI) source data are available from https://data-archive.nimh.nih.gov/oai/. The datasets generated and/or analyzed within this publication are available from the corresponding author on reasonable request.

Competing Interests

AB and MB are employees of Imorphics, Manchester, UK. DM, MM, DR, and YK are employees of OrthoTrophix, Inc., Oakland, CA. NS reports no competing interests.

Funding Sources

This study is sponsored by Orthotrophix, Inc., CA, USA. The study sponsor was involved in the study design, collection, analysis, and interpretation of data and in the writing of the manuscript and decision to submit for publication.

Author Contributions

DM was responsible for the conception and design of the study. AB was responsible for initial drafting the article. DM, YK, MM, DR, AB, MB and NS provided critical revision of the article for important intellectual content. All authors were responsible for analysis and interpretation of the data and for final approval of the article and take responsibility for the integrity of the work.

Acknowledgments

The authors would like to thank Felix Eckstein and the team at Chondrometrics for producing the cartilage measures. The authors would also like to thank Charles S. Davis PhD for his invaluable contribution to the statistical analysis.

Hopwood B, Tsykin A, Findlay DM, et al. Microarray gene expression profiling of osteoarthritic bone suggests altered bone remodelling, WNT and transforming growth factor-β/bone morphogenic protein signalling. Arthritis Res Ther. 2007;9:R100.
McGuire D, Lane N, Segal N, et al. TPX-100 leads to marked, sustained improvements in subjects with knee osteoarthritis: pre-clinical rationale and results of a controlled clinical trial. Osteoarthr Cartil. 2018;26:243.
Chen C, Thuillier D, Chin E. et.al. Chondrocyte-intrinsic Smad3 represses Runx2-inducible MMP-13 expression to maintain articular cartilage and prevent osteoarthritis. Arthritis Rheuml. 2012;64:3278–89.
Dai G, Xiao H, Liao J, et al. Osteocyte TGFβ1-Smad2/3 is positively associated with bone turnover parameters in subchondral bone of advanced osteoarthritis. Int'l J of Mol Medicine. 2020;46:167–8.
McGuire D, Segal N, Metyas S, et al. Intra-Articular TPX-100 in Knee Osteoarthritis: Robust Functional Response at 6 and 12 Months Is Associated with Increased Tibiofemoral Cartilage Thickness [abstract], Arthritis Rheumatol. 2018; 70 (suppl 10).
Neogi T, Bowes MA, Niu J, et al. Magnetic Resonance Imaging-Based Three-Dimensional Bone Shape of the Knee Predicts Onset of Knee Osteoarthritis: Data From the Osteoarthritis Initiative. Arthritis Rheum. 2018;65:2048–58.
Barr AJ, Dube B, Hensor E, et al. The relationship between three-dimensional knee MRI bone shape and total knee replacement—a case control study: data from the Osteoarthritis Initiative. Rheumatology. 2016;55:1585–93.
Bowes MA, Vincent GR, Wolstenholme CB, et al. A novel method for bone area measurement provides new insights into osteoarthritis and its progression. Ann Rheum Dis. 2013;74:519–25.
Hunter D, Nevitt M, Lynch J, et al. Longitudinal validation of periarticular bone area and 3D shape as biomarkers for knee OA progression? Data from the FNIH OA Biomarkers Consortium. Ann Rheum Dis. 2016;75:1607–14.
Bowes MA, Kacena K, Alabas OA, et al. Machine-learning, MRI bone shape and important clinical outcomes in osteoarthritis: data from the Osteoarthritis Initiative. Ann Rheum Dis annrheumdis-2020-217160, 2020.
Eckstein F, Ateshian G, Burgkart R, et al. Proposal for a nomenclature for Magnetic Resonance Imaging based measures of articular cartilage in osteoarthritis. Osteoarthr Cartil. 2006;14:974–83.
Wirth W, Eckstein F. A Technique for Regional Analysis of Femorotibial Cartilage Thickness Based on Quantitative Magnetic Resonance Imaging. IEEE Trans Med Imaging. 2008;27:737–44.
Williams TG, Holmes A, Waterton J, et al. Anatomically Corresponded Regional Analysis of Cartilage in Asymptomatic and Osteoarthritic Knees by Statistical Shape Modelling of the Bone. IEEE Trans Med Imaging. 2010;29:1541–59.
Hunter DJ, Nevitt M, Lynch J, et al. Preliminary assessment of predictive validity of periarticular bone area and shape markers in knee OA. Osteoarthr Cartil 2013 2013; 21: S175.
Bowes MA, Lohmander LS, Wolstenholme CBH, et al. Marked and rapid change of bone shape in acutely ACL injured knees – an exploratory analysis of the Kanon trial. Osteoarthr Cartil. 2019;January:1–8.
Conaghan PG, Bowes MA, Kingsbury S, et al. Disease-Modifying Effects of a Novel Cathepsin K Inhibitor in Osteoarthritis. Ann Intern Med 2019; Dec.:1–7.

Download PDF

Review #1 received at journal
20 Jun, 2021
Reviewer #2 agreed at journal
07 Jun, 2021
Reviews received at journal
03 Jun, 2021
Reviewer #1 agreed at journal
02 Jun, 2021
Reviewers invited by journal
28 Feb, 2021
Editor assigned by journal
21 Feb, 2021
Submission checks completed at journal
21 Feb, 2021
Editor invited by journal
21 Feb, 2021
First submitted to journal
19 Feb, 2021

You are reading this latest preprint version

Study TPX-100-5: Intra-Articular TPX-100 Significantly Delays Pathological Bone Shape Change and Stabilizes Cartilage in Moderate to Severe Bilateral Knee OA

Status:

Version 1

Abstract

Figures

Background

Methods

Results

Discussion

Conclusions

Declarations

Availability of data and materials

References

Status:

Version 1