Objectives Osteoarthritis (OA) patient stratification is an important challenge to design tailored treatments and drive drug development. Biochemical markers reflecting joint tissue turnover were measured in the IMI-APPROACH cohort at baseline and analysed using a machine learning approach in order to study OA-dominant phenotypes driven by the endotype-related clusters and discover the driving features and their disease-context meaning.
Method Data quality assessment was performed to design appropriate data preprocessing techniques. The k-means clustering algorithm was used to find dominant subgroups of patients based on the biochemical markers data. Classification models were trained to predict cluster membership, and Explainable AI techniques were used to interpret these to reveal the driving factors behind each cluster and identify phenotypes. Statistical analysis was performed to compare differences between clusters with respect to other markers in the IMI-APPROACH cohort and the longitudinal disease progression.
Results Three dominant endotypes were found, associated with three phenotypes: C1) low tissue turnover (low repair and articular cartilage/subchondral bone turnover), C2) structural damage (high bone formation/resorption, cartilage degradation) and C3) systemic inflammation (joint tissue degradation, inflammation, cartilage degradation). The method achieved consistent results in the FNIH/OAI cohort. C1 had the highest proportion of non-progressors. C2 was mostly linked to longitudinal structural progression, and C3 was linked to sustained or progressive pain.
Conclusions This work supports the existence of differential phenotypes in OA. The biomarker approach could potentially drive stratification for OA clinical trials and contribute to precision medicine strategies for OA progression in the future.
Trial registration number NCT03883568.
Data availability statement
Data are available on reasonable request. Data are available on request to the APPROACH Steering Committee.
Statistics from Altmetric.com
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.
What is already known about this subject?
There is an unmet need for new therapies that target the underlying pathology in osteoarthritis (OA).
Computational methods based on unsupervised machine learning have the potential to stratify OA cohorts into subsets that correspond to distinct molecular endotypes.
What does this study add?
By applying these methods to the IMI-APPROACH cohort, we identified three dominant clusters and characterised them as inflammatory, low-repair and subchondral bone/articular cartilage-driven phenotypes.
Patients in the discovered clusters had statistically significant differences in clinical characteristics.
How might this impact on clinical practice or future developments?
The biomarker-based endotype discovery approach could potentially drive stratification for OA clinical trials and contribute in the future to precision medicine strategies for OA care.
Osteoarthritis (OA) is the most common form of arthritis among older people, affecting more than 500 million people (7% of the global population).1 It is one of the most frequent causes of physical disability among older individuals and a major contributor to healthcare and societal costs globally.2 The risk factors for the development of OA include age, sex, obesity, previous joint injuries, repeated stress on the joint, malalignment, genetics, bone shape (including deformities) and certain metabolic diseases.3 According to studies on the global burden of disease, knee OA represents the greatest burden.4 5 However, despite the ever-increasing rise in the incidence and burden of OA, there is an unmet need for new therapies that target the underlying pathophysiologies.6 The currently available pharmacological treatments are only able to target the symptoms of OA, and they have adverse side effects, especially in older adults with common comorbidities.
The development of effective treatments and disease-modifying OA drugs (DMOADs) for this debilitating condition is extremely challenging.7 Many of the approaches that have been tried thus far have either failed or produced unsatisfactory outcomes. One of the greatest challenges in OA drug development is the heterogeneity of the disease.8 9 However, despite being a multifaceted and heterogeneous syndrome, there is an opportunity to target different treatments to patients according to their disease drivers characterised by molecular endotypes (a description of a subset of patients with common molecular characteristics) and clinical phenotypes (an observable characteristic or trait of a disease).9 10 OA may be amenable to tailored treatments that target specific phenotypes, including inflammatory, low repair, subchondral bone, metabolic or articular cartilage-driven phenotypes.11–17
Therefore, development of computational tools that includes objectively measured markers, such as biochemical markers, may facilitate OA drug development through patient subgrouping based on endotypic characteristics.9 11 An example of an OA endotype could be a group of patients with elevated bone biochemical markers, as compared with the remaining of the OA population. Then, based on the link with clinical data, this subgroup could be annotated as having a bone-driven disease (ie, a OA disease phenotype), and hypothetically, this group of patients should be enriched for in clinical trials testing the efficacy of a bone-modulating drug.18
At present, defining the appropriate outcome measures that are needed for OA clinical trials and the objective assessment of new therapies is challenging.19 Therefore, new computational methods based on machine learning (ML) and big data analytics can help advance this field of research by enabling protocols for patient classification into subtypes, using a combination of clinical, biochemical and/or imaging data.20–22
The aim of this study was to develop a methodology based on unsupervised ML (specifically, clustering) to identify/discover OA endotypes in the IMI-APPROACH cohort of patients with knee OA from a panel of 16 biochemical markers related to different joint tissue processes (eg, degradation, formation or inflammation), measured at the baseline of the study. The properties of the discovered clusters were thoroughly analysed using a combination of statistical and ML techniques, and the consistency of the discovered clusters was validated using data from an external cohort.
Cohort description and data collection
Applied Public-Private Research enabling OsteoArthritis Clinical Headway funded by the Innovative Medicines Initiative (IMI-APPROACH, trial registration number: NCT03883568) is a prospective cohort study including 297 patients with tibiofemoral OA according to the American College of Rheumatology classification criteria. Patients were (pre)selected from existing cohorts using ML models, developed on data from the CHECK cohort, to display a high likelihood of radiographic joint space width (JSW) loss and/or knee pain progression.23 24 The ultimate objective of APPROACH is to use real-world data to develop analysis methodologies to define disease subtypes and identify different knee OA clusters/phenotypes, to allow targeted treatment.
The IMI-APPROACH cohort screened 433 patients with OA (at five centres: Utrecht and Leiden, The Netherlands; A Coruña, Spain; Paris, France; Oslo, Norway) and enrolled 297 patients most likely to be pain and/or structural progressors at 2-year follow-up.24 Enrolled patients were predominantly women (n=230), predominantly Caucasian/white (n=283), aged 44–82 years (median age: 67.5, IQR 62–71 years) and mostly overweight (median body mass index (BMI): 27 kg/m2, IQR 24.4–31.6). At baseline, serum (S) and urine (U) samples were collected for analyses of 16 biochemical markers (table 1). The biomarkers were measured in International Organization for Standardization-certified laboratories at Nordic Bioscience (S_RE_C1M, S_C2M, S_C3M, S_C10C, S_CRPM, S_PRO_C2, U_CTXII, S_CTXI, U_CTXI_ALPHA, S_NMID, S_HA, S_COMP and S_hsCRP), Artialis (S_COLL2_1 and S_COLL2_1NO2) and Lund University (S_ARGS). The list of biomarkers was selected based on present knowledge of joint tissue turnover and OA.
In addition to the biochemical markers data ( B ), extra information ( E ) was collected as part of the IMI-APPROACH cohort.23 These included assessment of radiographs of knees and hands, MRIs and CT scans of the knees, and outcomes of physical examinations and questionnaires: Function Index of Hand OA (FIHOA), Hip Disability and Osteoarthritis Outcome Score, Intermittent and Constant Osteoarthritis Pain Score, Knee-Injury and Osteoarthrosis Outcome Score (KOOS) and the 36-Item Short Form Health Survey. See online supplemental table B1. All data used in this paper were collected at the baseline visit of the study, except for the data on progression (Relation of clusters to progression section).
The biochemical markers data ( B ) were log transformed to account for long-tailed distributions. Missing data in B (<0.01% of values) were estimated (imputed) using either random Forest (RF) or k-nearest neighbour (KNN) regression models (see online supplemental appendix A, section 1.1).
As not all patients fasted before the sample collection, the fasting sensitivity of the biomarkers had to be assessed. The Spearman rank correlation with the patient’s fasting status was found to be weak, except for U_CTXI (r=0.41). The values for this biomarker were corrected with an imputation approach (see online supplemental appendix A, section 1.2). We opted for a model-agnostic correction (ie, correcting the data rather than altering the analysis model) because it is more suitable for the downstream ML analysis we performed. identifies the processed biomarkers data.
The extremes values of were trimmed with a combination of Tukey and Winsor methods25 to reduce the effect of outliers. Afterwards, principal component analysis was used to eliminate correlated biomarkers (see online supplemental appendix A, section 1.4). This resulted in 13 principal components which were found to explain 95% of data variance. These components were clustered using the k -means algorithm.26 The optimum value for k (number of clusters) was identified from the consensus of silhouette score, the j-score and adjusted mutual information score. To obtain a robust estimate of these metrics, for each the k-means algorithm was run 10 times with different random seeds (see online supplemental appendix A, section 1.5). The clustering with the highest quality was found for . We chose for the rest of the analysis in this paper as we aimed to investigate the highest number of meaningful clusters. The final cluster membership was taken from the algorithm run with the highest silhouette score for .
Using data in , we trained three RF classification models predicting membership to each cluster (one cluster vs the rest) and then interpreted the model decisions using the SHAP (SHapley Additive exPlanations) TreeExplainer method,27 to understand which variables determine the cluster membership. RF hyperparameters were tuned through a nested cross-validation procedure with recursive feature elimination (RFE-CV). See online supplemental appendix A, section 1.6 for more details.
Statistical analysis of cluster differences
To further describe the clusters, statistical tests were conducted for each feature in and E , to assess whether the clusters had statistically different distributions for individual markers. The Mann-Whitney U test was used for continuous and ordinal features, and the χ2 test for categorical ones. The clusters were compared pairwise, and the null hypothesis was rejected following the Benjamini-Hochberg correction procedure for multiple comparisons applied across features.28 The features in were inverse log transformed to operate on actual biomarker concentrations (see online supplemental appendix A, section 1.3, for normality tests).
Figure 1 shows an overview of the entire data analysis pipeline described in this section, including data preprocessing, clustering, cluster’s interpretation and the statistical analysis.
Validation on an external cohort
The proposed clustering pipeline was also applied to FNIH/OAI. The FNIH/OAI is the largest available OA cohort that was similar to IMI-APPROACH in terms of biomarkers.29 The two cohorts had 11 biomarkers in common. Incurrent sample remeasurement for the adjustment of technical batch effects could not be performed, as no samples were left and available from the FNIH/OAI cohort for this purpose. Therefore precise data alignment on the absolute mean concentrations and variance between the two cohorts was not possible to conduct.30 31 As a result, the only possible type of external validation consisted in replicating the clustering pipeline for the two cohorts restricted to the common set of 11 biomarkers and evaluating the consistency of the identified clusters across cohorts.
Potential age and sex-based bias
To investigate the potential bias of age and gender in the clustering process, we statistically analysed the differences in age and sex across clusters, and we applied our clustering pipeline separately to the male and female subcohorts for both IMI-APPROACH and FNIH/OAI, to assess the consistency across clusters.
Relation of clusters to progression
To verify a relation between the clusters and disease progression, we used 2-year follow-up data to decide for each patient whether and how they have progressed, available only for a subset of 221 IMI-APPROACH participants. We defined one non-progressive category and three progressive categories related to pain, structure, and combined pain and structure.23 24 Then we analysed the distribution of progressors in each cluster. See online supplemental appendix E for more details.
Our clustering pipeline identified three clusters. These are shown in figure 2 as a two-dimensional projection obtained with UMAP (Uniform Manifold Approximation and Projection).32 UMAP hyperparameters were optimised via grid search to maximise the two-dimensional silhouette score. The projection preserves the local neighbourhood structure and gives an idea of the strength of the global separation between the clusters in the original multidimensional space.
The classification models trained to predict patient’s cluster membership achieved high F1 scores (C1 vs rest: 0.85, C2 vs rest: 0.91, C3 vs rest: 0.88). As a result, the subsequently performed model interpretation was expected to be meaningful. Figure 3 shows which biomarkers were predominantly used by each model to decide the cluster membership. Figure 4 compares the median biomarker concentrations for each cluster in a radar plot. Figure 5 shows the differences in biomarker value distributions across clusters. Bringing all these results together, the three clusters were interpreted as follows:
Cluster 1 represents a low tissue turnover phenotype: patients have all the inflammation and structural damage related biomarkers in the mid/low ranges.
Cluster 2 represents a structural damage phenotype: patients have high values of the bone and cartilage markers: S_CTXI, U_CTXIALPHA, S_NMID and U_CTXII.
Cluster 3 represents a systemic inflammation phenotype: patients have high values of the inflammatory and MMP-driven markers: S_hsCRP, S_RE_C1M, S_CRPM and S_C3M. In contrast, these patients show low values of bone and cartilage related markers: U_CTXIALPHA, S_NMID, and S_CTXI.
The clustering stability was investigated by comparing the results obtained for with those obtained for and . We found that clusters and interpretation were reasonably preserved at least until . This demonstrates that the three clusters analysed in this work are well-defined in the data space and robust with respect to finer clustering (see online supplemental appendix A, section 1.7).
Statistical analysis of differences between clusters
Several statistically significant differences in clinical scores were found. Full results are provided in online supplemental appendix B, and here we only present highlights of those findings. All figures cited in this section are provided in online supplemental appendix B.
Clusters 2 and 3 had a higher percentage of women than cluster 1, and cluster 3 had a higher mean BMI (online supplemental figure B15).
There was no difference in median age and range, smoking status, comorbidities and use of OA medication (online supplemental figure B14) across the clusters.
Cluster 3 had statistically more patients experiencing substantial pain when standing (KOOS_P09, online supplemental figure B9), burning sensation (pain detect 09, online supplemental appendix B, B14) and more pain now and on average over the past 4 weeks (pain detect 01 and 03, online supplemental figure B14) than clusters 1 and 2. Patients in cluster 2 also experienced more pain in the past week than those in cluster 1 (pain detect 03, B14). Maximum Numeric Rating Scale (NRS) pain for hands were higher in cluster 3 (online supplemental figure B15), as well as having worse overall health self-assessment (SF36_11d, online supplemental figure B17).
Cluster 1 has higher knee JSW (mean) than cluster 2 and less severe carpometacarpal Kellgren-Lawrence scores compared with cluster 3 (online supplemental figure B17).
External validation using FNIH/OAI data
We reduced the set of data features to the common subset of 11 biomarkers across the IMI-APPROACH and FNIH/OAI cohorts and applied the same clustering pipeline to both datasets. Figure 6 shows the comparison of obtained clusters. Despite the removal of five biomarkers, the IMI-APPROACH clusters still corresponded to structural damage, inflammatory and low tissue turnover endotypes. The FNIH/OAI clusters were found to consistently exhibit the same patterns, demonstrating cross-cohort robustness of our approach (see online supplemental appendix C).
Analysis of age and gender-based bias
We found no statistical difference between clusters in terms of age, as well as no statistical difference between male and female subcohorts in terms of age (see online supplemental figure D1). However, the male and female subcohorts had statistically different distributions for the following eight biomarkers: S_ARGS, S_C10C, S_COLL2_1, S_COLL2_1NO2, S_CTXI, S_NMID, U_CTXII and U_CTXI_ALPHA. Moreover, the clusters were significantly different in terms of gender, suggesting that it plays an important role in driving the clustering results (online supplemental figure D4). Similar patterns could be found for the FNIH/OAI cohort (see online supplemental appendix D).
Relation of clusters to progression
Table 2 summarises the progression status of the clusters. While we found progressors in all clusters, they were not distributed uniformly by progression type. There was more pain-related progressors and combined pain and structure progressors in the inflammation cluster (C3). Similarly, there were more structure-related progressors in the structural damage cluster (C2). The highest relative number of non-progressive patients was found in the low tissue turnover cluster (C1) and the lowest in the inflammation cluster (C3).
The aim of this work was to test if ML techniques can be used to identify biologically meaningful subgroups of patients with OA in the APPROACH cohort based on selected biochemical markers. By using clustering, that is, an unsupervised ML approach that does not exploit domain knowledge, we were able to identify molecular endotypes from 16 well-defined biochemical markers reflecting different molecular pathways and ongoing pathophysiological processes. We discovered three distinct OA phenotypes associated with the clusters (endotypes): C1—a low tissue turnover phenotype, C2—a structural damage phenotype and C3—systemic inflammation phenotype. The clustering reflects well the current biological and mechanistic understanding of the respective biomarkers, in that distinct patterns could be identified for the subtypes. In particular, the combination of different markers describes the underlying biology in the clusters. This result is in line with published results from the FNIH/OAI biomarker initiative,29 33 and the progression status of the members of each cluster is consistent with the cluster interpretation provided above: C1 has the highest proportion of non-progressors, C2 has the highest proportion of structural progressors and C3 has the highest proportion of pain-related progressors, and those progressing both in pain and in structure. However, although the proportions varied (ranging from 43% to 56%) progressive patients were found in all clusters. This means that the clusters represent different disease subtypes, within which the progression may occur.
Putting this in context of the work conducted on markers in clinical interventional trials, a few things can be learnt. Oral salmon calcitonin was tested as an antiresorptive treatment for OA. The phase III clinical trials failed to meet their clinical endpoints. Interestingly, calcitonin did significantly modulate CTX-I and CTX-II.34 There are likely several reasons why this study failed, however, it begs to wonder what would the outcome have been if the study was enriched for C2 patients? Another failed trial was testing the efficacy of the IL-1 monoclonal antibody in OA and found markers from C3 modulated by treatment.35 Would it still fail if it was enriched for C3 patients?
Despite a large and growing disease burden in OA, many pharmaceutical companies have de-emphasized or even abandoned OA drug development due to perceived hurdles. Crucial in this is the lack of appropriate predictive and outcome measures that can robustly identify patients early in the disease, which may benefit from a specific therapy. The lack of specific and sensitive baseline characteristics and subsequent endpoints to differentiate between responders and non-responders, both at the level of pain and tissue structure modification (ie, DMOAD), has led to trials that included hundreds of patients in each arm with at least 3-year follow-up. Despite these enormous trials, European Medicines Agency and Food and Drug Administration have not approved any DMOAD yet.36 There is a general lack of understanding of OA pathogenesis which appears rather variable and likely reflects different phenotypes with fundamental differences in disease aetiology, tissue alterations, clinical manifestations (pain/mobility) and disease progression. Although the current mindset for drug treatment in the field is moving to a more personalised medicine and patient stratification approach, there are no accepted methods or guidelines to classify patients with OA, for example, to predict the underlying pathophysiology, to select patients according to their prognosis or to differentiate between patients in terms of diagnosis methodology and treatment plan. However, several initiatives have been initiated to generate more focus on the development of projects for identifying endotypes. For example, a framework for conducting and reporting phenotyping research was provided37—this may very well be the first step toward integrating the concept of phenotyping in research.
A better understanding of disease stratification and acceptance of a guideline to classify patients with OA will provide clear phenotype-directed protocols for DMOAD trials that enable us to target subgroups with OA that have uniform disease characteristics, thereby increasing the chances of success. We propose that the biomarker clustering analysis performed herein can be used to stratify patients with OA into groups with distinct molecular endotypes. This approach could potentially drive OA clinical trials stratification and serve as the basis for precision medicine strategies for OA progression in the future. Although there are limited data publicly available, there have been a few attempts to identify multimarker endotypes in OA. Sonh et al showed that several cytokines were elevated in synovial fluid and serum of patients with OA compared with normal samples when looking at an average level; however, it was also obvious that the pattern was very heterogeneous.38 Werdyani et al identified three distinct endotypes using metabolomics.39 One of those clusters showed some association with muscle weakness. These data suggest that a subset of patients could belong to an inflammatory endotype.
Moreover, we focused on biochemical markers measured at the baseline of the study, and not their longitudinal changes, as this analysis would be more useful to inform future clinical trials. Longitudinal monitoring of biomarkers can give insight in the pharmacodynamic effects or provide early proof of effectiveness of a compound in interventional clinical trials, however often fail to predict progression in the study population in these trials.34 40–42 Therefore, although longitudinal monitoring of individual biomarkers are only modestly predictive (if at all) of knee OA progression, they might have some utility as patient stratification like described herein for enriching OA trials for progressors.29
As more longitudinal data of the IMI-APPROACH cohort becomes available (currently an ongoing process), future investigations could explore the longitudinal data on biomarkers, imaging and other markers in IMI-APPROACH to further refine the description of the phenotypes and possibly explore more detailed stratifications. This analysis could take many different directions, for example, analyse cluster membership differences between visits or on comparison of the entire patient trajectories over 2 years of the study.
The main limitations of this work were the small numbers of patients in the IMI-APPROACH cohort and being able to perform only a partial validation with an external cohort, limited to a common subset of biomarkers. It would be beneficial for the field if future biomarker studies use a superset of the FNIH/OAI and IMI-APPROACH biomarkers, to allow for a complete validation of the discovered clusters. The use of predefined set of biochemical markers limits the discovery potential to certain molecular mechanisms. This could be avoided if clustering was performed on data generated by an untargeted platform (eg, RNA-seq); however, the analysis of such high-dimensional data is often much less robust, especially on small sample sizes. Finally, more research should be conducted on more abundant cohorts to fully evaluate the gender bias in clustering analysis of OA-related biochemical markers. From our analysis in the IMI-APPROACH and FNIH/OAI cohorts, we believe it is advisable for future studies to consider male and female patients separately and possibly draw conclusions that are gender based, if sample sizes are large enough.
Data availability statement
Data are available on reasonable request. Data are available on request to the APPROACH Steering Committee.
Patient consent for publication
This study involves human participants. The study is being conducted in compliance with the protocol, Good Clinical Practice (GCP), the Declaration of Helsinki, and applicable ethical and legal regulatory requirements (for all countries involved), and is registered under ClinicalTrials.gov: NCT03883568. All participants have received oral and written information and provided written informed consent. Participants gave informed consent to participate in the study before taking part.
We thank Agnes Lalande (Institut de Recherches Internationales Servier, Suresnes, France) for clinical study coordination, Thomas Lohr (GSK) for sample handling logistics and Janneke Boere and colleagues at Lygature for coordinating the research activity.
Handling editor Josef S Smolen
Twitter @larhumato, @jaumebp
Contributors FA, PW, JBa contributed to the conception and design of the study, AS, MU, YH, ACM, MK, FJB, IKH, FB and ACBJ contributed to the acquisition of data. AM, JBl, CL, JL and ACB-J contributed to the analysis and interpretation of data. JBa is the paper's guarantor. All authors contributed to drafting the article or revising it critically for important intellectual content. All authors gave final approval of the version to be submitted.
Funding The research leading to these results has received support from the Innovative Medicines Initiative Joint Undertaking under Grant Agreement no 115770, resources of which are composed of financial contributions from the European Union’s Seventh Framework Programme (FP7/2007–2013) and EFPIA companies’ in-kind contribution. See http://www.imi.europa.eu/ and http://wwwapproachprojecteu/.
Competing interests ACB-J is a full-time employee and shareholder of Nordic Bioscience, a privately owned company involved in the development and commercialisation of biomarkers for fibroinflammatory disorders. YH is the founder and president of Artialis, and MU is a full-time employéé of Artialis, a spin-off company of the University of Liège. YH has also received fees from Tilman, Genequine, Seikagaku, Expanscience, Nestlé, Immubio, Biose and Labhra. CL was an employee of Merck at project start. IKH consults for Abbvie and Novartis and has received funding from Pfizer. JL is employed by and shareholder in GlaxoSmithKline. FB reports personal fees from AstraZeneca, Boehringer, Bone Therapeutics, CellProthera, Expanscience, Galapagos, Gilead, Grunenthal, GSK, Eli Lilly, Merck Sereno, MSD, Nordic, Nordic Bioscience, Novartis, Pfizer, Roche, Sandoz, Sanofi, Servier, UCB, Peptinov, 4P Pharma, 4Moving Biotech and grants from TRB Chemedica, outside the submitted work. FJB reports funding from Gedeon Richter, Bristol-Myers Squibb, Sun Pharma Global FZE, Celgene, Janssen Cilag, Janssen Research & Development, Viela Bio, Astrazeneca, UCB BIOSCIENCES, UCB BIOPHARMA SPRL, AbbVie Deutschland, Merck, Amgen, Novartis Farmacéutica, Boehringer Ingelheim España, CSL Behring, Glaxosmithkline Research & Development, Pfizer, Lilly, Corbus Pharmaceuticals, Biohope Scientific Solutions for Human Health, Centrexion Therapeutics, Sanofi, TEDEC-MEIJI FARMA, Kiniksa Pharmaceuticals, Fundación para la Investigación Biomédica Del Hospital Clínico San Carlos, Grünenthal and Galapagos. MK receives consulting fees from Abbvie, Pfizer, Levicept, GlaxoSmithKline, Merck-Serono, Kiniksa, Flexion, Galapagos, Jansen, CHDR, Novartis, UCB. AM receives fees/funding from Merck, Kolon TissueGene, Pfizer, Galapagos-Servier, Image Analysis Group (IAG), Artialis, Aché Laboratórios Farmacêuticos, AbbVie, Guidepoint Global, Alphasights, Science Branding Communications, GSK, Flexion Therapeutics, Pacira Biosciences, Sterifarma, Bioiberica, SANOFI, Genacol, Kolon Life Science, BRASIT/BRASOS, GEOS, MCI Group, Alcimed, Abbot, Laboratoires Expansciences, SPRIM Communications, Frontiers Media and University Health Network Toronto.
Provenance and peer review Not commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.