Article Text

Risk of malignant lymphomas in patients with inflammatory bowel disease: a population-based cohort study
  1. Jingru Yu1,
  2. Erle Refsum2,3,
  3. Paulina Wieszczy2,3,4,
  4. Lise M Helsingen2,3,
  5. Vera Perrin1,
  6. Amanda Högdén1,
  7. Magnus Løberg2,3,
  8. Johannes Blom5,
  9. Michael Bretthauer2,3,
  10. Hans-Olov Adami1,2,3,
  11. Weimin Ye1,6,
  12. Mette Kalager2,3
  1. 1Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden
  2. 2Clinical Effectiveness Research Group, Department of Transplantation Medicine, Oslo University Hospital, Oslo, Norway
  3. 3Clinical Effectiveness Research Group, Institute of Health and Society, Department of Health Management and Health Economics, University of Oslo, Oslo, Norway
  4. 4Department of Gastroenterology, Hepatology and Clinical Oncology, Centre of Postgraduate Medical Education, Warsaw, Poland
  5. 5Department of Clinical Science and Education, Södersjukhuset, Karolinska Institutet, Stockholm, Sweden
  6. 6Department of Epidemiology and Health Statistics & Key Laboratory of Ministry of Education for Gastrointestinal Cancer, Fujian Medical University, Fuzhou, Fujian, P.R. China
  1. Correspondence to Dr Erle Refsum; erle.refsum{at}medisin.uio.no

Abstract

Objective To estimate the risk of non-Hodgkin’s lymphoma (NHL) and Hodgkin’s lymphoma (HL) in patients with inflammatory bowel disease (IBD).

Design We undertook a two-country population cohort study with all patients diagnosed with IBD in Norway and Sweden from 1987 and 1993 through 2015 and 2016, respectively, and analysed the risk of NHL and HL. In Sweden, we also analysed prescriptions of thiopurines and anti-tumour necrosis factor (TNF)-α therapy from 2005. We calculated standardised incidence ratios (SIRs) with 95% CIs using the general populations as reference.

Results Among 131 492 patients with IBD with a medium follow-up of 9.6 years, we identified 369 cases of NHL and 44 cases of HL. The SIR of NHL was 1.3 (95% CI 1.1 to 1.5) in ulcerative colitis and 1.4 (95% CI 1.2 to 1.7) in Crohn’s disease. We found no compelling heterogeneity in analyses stratified by patient characteristics. We found a similar pattern and magnitude of excess risks for HL. At 10 years, cumulative incidence was 0.26% (95% CI 0.23% to 0.30%) and 0.06% (95% CI 0.04% to 0.08%) for NHL and HL, respectively. Higher excess risks were found among patients with NHL with concomitant primary sclerosing cholangitis (SIR 3.4; 95% CI 2.1 to 5.2) and in those prescribed thiopurines alone (SIR 2.8; 95% CI 1.4 to 5.7) or with anti-TNF-α agents (SIR 5.7; 95% CI 2.7 to 11.9).

Conclusion Patients with IBD have a statistically significant increased risk of malignant lymphomas compared with the general population, but the absolute risk remains low.

  • lymphoma
  • inflammatory bowel disease
  • epidemiology

Data availability statement

Data are not publicly available due to their containing information that could compromise the privacy of research participants. Summary data can be made available from the corresponding author on reasonable request.

https://creativecommons.org/licenses/by/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/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.

WHAT IS ALREADY KNOWN ABOUT THIS TOPIC

  • Several autoimmune diseases are associated with increased risk lymphoma, but the risk in patients with inflammatory bowel disease remains unclear due to inconsistent evidence.

WHAT THIS STUDY ADDS

  • We present a large, population-based study with long-term and near complete follow-up.

  • We demonstrate that the absolute risk of lymphoma in inflammatory bowel disease is low.

HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE, OR POLICY

  • There is a need for even larger studies with detailed data on disease characteristics and treatment to quantify the risk associated with immunosuppressive therapies.

  • Such studies need to consider the effect of confounding by indication, which requires valid data on disease activity.

Introduction

An increased risk of Hodgkin’s lymphoma (HL) and non-Hodgkin’s lymphoma (NHL) has been documented in several autoimmune diseases, such as Sjögren’s syndrome, coeliac disease and rheumatoid arthritis.1–4 A link between disease severity and lymphoma risk has been established, and although the exact molecular mechanisms remain to be explained, longstanding chronic inflammation with B cell and T cell activation and/or stimulation are likely predisposing factors. Several mechanisms may be ongoing in parallel and the type of immune dysfunction involved in lymphoma development is thought to be disease specific, as different types of lymphoma are associated with the different disease entities.4 Immunosuppression is generally thought to play less of a role in lymphoma development in most autoimmune diseases, although such drugs may increase the risk of cancer through reduced immunological surveillance.4

Inflammatory bowel disease (IBD) is also an autoimmune disease characterised by chronic inflammation, but, contrasting to other autoimmune diseases, the available evidence on risk of lymphoma is scarce and too inconsistent to allow reliable quantification and guide patient care.5 6 Meta-analyses have concluded with both increased7 8 and no increased9 risk of lymphoma in patients with IBD treated with immunosuppressants. Limited statistical power, short follow-up, conflation of different haematopoietic malignancies, and lack of detailed data on immunosuppressive therapy are limitations of published studies.1 3 10–17 With these limitations, overall the available literature points towards a minimally increased risk in untreated patients with IBD18 and an increased risk among patients with IBD treated with immunomodulators, particularly thiopurine analogues.18–21 To what extent the disease activity itself or the use of different immunosuppressants drive an increased risk of lymphoma in IBD is still unclear.21

To add to the existing literature on risk of lymphoma in patients with IBD, we undertook a nationwide population study of IBD and lymphoma risk in Norway and Sweden. We took advantage of cross-linkage between numerous high-quality registries to empanel a large patient cohort with complete, long-term follow-up.

Methods

Study cohort

Norway and Sweden are countries with similar populations with regard to distribution of age, socioeconomic factors, ethnicity, and access to public healthcare systems. We defined our study population using data from databases of all public hospitals in Norway (from 1 January 1987 to 31 December 2015), and from the Swedish Patient Register (from 1 January 1993 to 31 December 2016).

We used the individually unique national registration numbers in registries in each country to link patients with IBD to the following public registries for data on outcomes and censoring: cancer (Cancer Registry of Norway, the Swedish Cancer Register), mortality (the Norwegian Cause of Death Registry, the Swedish Cause of Death Register), demographics and migration (the total Population Register in both countries). The empanelled registries are nationwide and entail compulsory, automated registration of all patients and are close to 100% complete.

Patients were followed from the date of diagnosis of IBD until diagnosis of lymphoma, diagnosis of any other cancer, emigration, death, or end of follow-up (Norway: 31 December 2015; Sweden: 31 December 2016), whichever occurred first. Patients with a cancer diagnosis prior to IBD diagnosis were excluded. A flow chart summarising the enrolment procedures is shown in online supplemental figure 1.

To ensure high specificity of IBD diagnosis, we required at least two inpatient and/or outpatient diagnosis codes of IBD. This approach yields a positive predictive value of an IBD diagnosis of 90% in our cohort.22 We used the International Classification of Disease (ICD) codes from the first two visits to define IBD subtypes at baseline: patients were classified as having ulcerative colitis (ICD-9: 556; ICD-10: K51) or Crohn’s disease (ICD-9: 555; ICD-10: K50) if their first two codes were for ulcerative colitis or Crohn’s disease, respectively, and unclassified IBD if their first two codes were indeterminate colitis (ICD-10 K52.3) or any combination of codes (online supplemental table 1). In main analyses, the date of the second diagnostic listing was defined as date of IBD onset.

Comparator

The general population of Norway and Sweden served as reference, using anonymised individual data of all cancer diagnoses in Norway and Sweden diagnosed during the study period. The datasets include information on cancer topography and morphology, in addition to date of diagnosis, patient sex and birth year. Similar to the patients with IBD, we required that incident cases of lymphoma in the general population had no history of cancer.

Covariates

In Norway, outpatient diagnoses became gradually available in hospital databases by geographic roll-out during the study period, while in Sweden, outpatient data from the whole country were available from 2001; patients diagnosed up until 2003 could represent a mix of patients with incident and prevalent IBD. We therefore analysed the data by calendar year periods of enrolment (1987/1993 to 2002 and 2003 to 2015/2016).

We used the Montreal classification to define extent of disease during follow-up, based on ICD-9 and ICD-10 codes.23 The most severe disease extension during follow-up was used to classify patients (online supplemental table 2). Indicators of disease severity such as laboratory markers, endoscopic or clinical scores were not available in our databases. We therefore used bowel surgery during follow-up as an indicator of disease severity (online supplemental table 3). Primary sclerosing cholangitis (PSC) was defined using ICD-9 (Norway: 5761; Sweden: 576B) and ICD-10 codes (Norway: K83.0; Sweden: K83.0). Bowel surgery and PSC were analysed as time-varying covariates.

Pharmacotherapy

Data on pharmacotherapy were available for all Swedish patients from the Swedish Drug Register from 2005 onwards. We focused analyses on thiopurines and anti-tumour necrosis factor (TNF)-α therapy, as these are frequently used and have been associated with lymphoma risk in previous studies.24 The anatomical therapeutic chemical (ATC) codes are listed in online supplemental table 4. During each day of follow-up, patients were categorised as being exposed to one or two classes of drugs. Patients were considered unexposed until the day of their first prescription. For each drug, any overlapping prescriptions within a 90-day period were bridged, and duration of exposure estimated to last until end of the last prescription. Duration of prescriptions was assumed to last 3 months. To minimise prevalent user bias, we included a 1-year wash-out window and excluded patients with their first prescription in the year prior to cohort entry (their second IBD diagnosis). To eliminate any reverse causality due to initiation of treatment due to cancer symptoms, we introduced a 1-year lag between initiation of treatment and the follow-up window. Additionally, the discontinuation time between prescriptions (where the interval was longer than >3 months) was removed, since patients in the discontinuation period were different from patients without any treatment or with each class of the treatments. See online supplemental figure 2 for design diagram. Individual pharmacotherapy data were not available in Norway.

Outcome

The cancer registries have different coding practices: in Norway, all cancer diagnoses are retrospectively coded according to ICD-O-3. In Sweden, ICD-O-2 and Systematized Nomenclature of Medicine (SNOMED) codes were introduced in the Swedish Cancer Register from 1 January 1993. To allow for accurate subgrouping of lymphomas, we restricted the IBD cohort in Sweden to patients diagnosed from that date (see online supplemental table 5 for outcome definitions in each country). Lymphomas were divided into NHL and HL, as the codes used in the Swedish Cancer Register do not allow adequate resolution to accurately define subcategories such as B-cell, T-cell or NK-cell lymphomas.

Statistical analyses

We calculated crude incidence rates (IRs) for each type of lymphoma, and present IRs standardised for age and sex to the Norwegian population of 1987 graphically. Cumulative incidence was calculated by Kaplan-Meier method, using the cumulative hazard function as a direct estimate. Patients were considered at risk from 1 year after IBD onset to reduce bias due to reverse causality.

We calculated standardised incidence ratios (SIRs, the ratios of observed cancer cases divided by expected number of cancer cases in the population) with 95% CIs, overall and stratified by covariates, to estimate the relative risk of lymphomas in patients with IBD compared with the general population. The expected number was derived by multiplying the person-years at risk with country-specific, age-specific (5-year strata), sex-specific, and calendar year-specific IRs in the general population.

The Cochran-Armitage trend test was used to evaluate cancer risk across age groups and follow-up duration. SIRs using data from Norway and Sweden were analysed separately and combined.25

Sensitivity analyses

In sensitivity analyses, patients were considered at risk from their second diagnostic listing of IBD (inpatient or outpatient), that is, the 1-year lag was removed.

Patient and public involvement statement

Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

Results

Patient characteristics

Our cohort consisted of 131 492 patients with IBD with a median follow-up of 9.6 years (IQR 4.3–14.9 years), and a total of 1 322 831 person-years of follow-up (table 1). In Sweden, 473 (1%) patients were included in analyses of anti-TNF-α therapy, 10 271 (25%) included in analyses of thiopurines and 3909 (10%) patients analysed as treated with a combination of these drugs since 2005 (after first year of follow-up).

Table 1

General characteristics of patients with IBD in Norway and Sweden, 1987–2016

Non-Hodgkin’s lymphoma

During follow-up, 369 patients were diagnosed with NHL, 231 patients with ulcerative colitis and 117 patients with Crohn’s disease (table 2).

Table 2

Crude incidence rates per 100 000 person-years, number of cancer cases and SIRs of non-Hodgkin's lymphoma and Hodgkin’s lymphoma in patients with IBD in Norway and Sweden, 1987–2016*

We observed a crude IR of NHL of 30.8 (95% CI 27.8 to 34.2) per 100 000 person-years, corresponding to three patients with lymphoma if 10 000 patients with IBD were followed for 1 year. We found an overall 30% excess risk as compared with the background population (SIR 1.3; 95% CI 1.2 to 1.5). Cumulative incidence of NHL was 0.26% after 10 years of follow-up (95% CI 0.23% to 0.30%) (figure 1).

Figure 1

Cumulative incidence of non-Hodgkin’s lymphoma (NHL) and Hodgkin’s lymphoma (HL) by year since onset of inflammatory bowel disease (IBD). Patients were at risk after 1 year of follow-up. Cumulative incidence was calculated directly from Nelson-Aalen cumulative hazard estimates.

We found a higher relative risk, but low IRs with younger age of onset of IBD, while the IRs were highest in those diagnosed at a later age, with a lower relative risk compared to the background. While the combined results showed no difference between the two study periods (before and after 2002), in Norway, the excess risk was most pronounced in the second half of the study period, while in Sweden the results were opposite (online supplemental table 6). Compared with the general population of Norway, patients with IBD in Norway had similar age and sex SIRs of NHL until 2004, when an excess risk became evident. In Sweden, the excess risk was more pronounced in the earlier years of the study period and remained increased (figure 2, panel A).

Figure 2

Age-standardised and sex-standardised incidence rates of non-Hodgkin's lymphoma (panel A) and Hodgkin’s lymphoma (panel B). The Norwegian background population of 1987 was used as reference. Panel A: the peak in 1997 in Sweden was due to cases of lymphoma in strata with few patients with inflammatory bowel disease (IBD). Moreover, the number of patients with incident IBD was artifactually low before 2001, when outpatient records became available in the Swedish Patient Register (online supplemental figure 2). Panel B: the figure is based on 14 cases in Norway and 29 cases in Sweden. The increased rates in 2002 and 2003 in Norway were due to two cases in strata with very few patients with IBD.

Among patients with ulcerative colitis, the excess risk was about 30% regardless of extent of disease. Likewise, among patients with Crohn’s disease we found no apparent relation between location or extent of disease and the excess risk of NHL, except that those with perianal disease encountered a markedly increased risk (SIR 2.4; 95% CI 1.7 to 3.3). Of all categories in table 2, patients diagnosed with PSC experienced the highest excess risk (SIR 3.4; 95% CI 2.1 to 5.2). For patients who underwent bowel surgery, the risk was only slightly increased (SIR 1.8; 95% CI 1.4 to 2.3). Although we excluded the first year of follow-up to reduce selection bias, the overall excess risk was higher 1–2 years after diagnosis compared with subsequent years with follow-up beyond 10 years.

Hodgkin’s lymphoma

We identified 44 incident cases of HL in patients with IBD; 27 in patients with ulcerative colitis and 16 in patients with Crohn’s disease. There was no clear pattern in IRs over time (figure 2, panel B). We found a crude IR of 3.7 (95% CI 2.7 to 4.9) per 100 000 person-years, corresponding to <1 case of lymphoma if 10 000 patients with IBD were followed for 1 year. Compared with the background, we found a 50% excess risk with an overall SIR of 1.5 (95% CI 1.1 to 2.1). The 10-year cumulative incidence of HL was 0.06% (95% CI 0.04% to 0.08%).

Rates and SIRs were higher among men and those diagnosed with IBD at or after 60 years of age. There were no significant differences for patients with perianal disease, PSC or bowel surgery. In Norway, the SIR was higher in the latter study period, while it was similar in the two periods in Sweden (online supplemental table 7).

Pharmacotherapy

Among patients with IBD in Sweden, 14 653 (36%) were treated with thiopurines or anti-TNF-α therapy from 2005. The IR of NHL was 10.5 (95% CI 1.5 to 74.3) per 100 000 person-years for patients not treated with thiopurines or anti-TNF-α therapy (SIR 0.9; 95% CI 0.6 to 1.3) (table 3). For patients treated with thiopurines or anti-TNF-α therapy, SIRs for NHL were significantly increased among patients on combination therapy with both drug classes (SIR 5.7; 95% CI 2.7 to 11.9), while monotherapy with thiopurines also increased the risk of NHL with SIR 2.8 (95% CI 1.4 to 5.7). For HL, the overall IR was 5.0 (95% CI 2.0 to 10.4) per 100 000 person-years with an SIR of 1.9 (95% CI 0.9 to 4.0).

Table 3

Crude incidence rates per 100 000 person-years, number of cancer cases and SIRs of non-Hodgkin's lymphoma and Hodgkin’s lymphoma in patients with IBD in Sweden, 2005–2016*

Sensitivity analysis

Analyses where the first year of follow-up was included showed similar results (online supplemental table 8 (NHL) and online supplemental table 9 (HL)).

Discussion

Our nationwide study in Norway and Sweden showed a 30% increased risk of NHL in patients with ulcerative colitis and Crohn’s disease. This excess risk was of similar magnitude and in most instances statistically significant across numerous subgroups defined by patient and disease characteristics. The risk pattern was largely similar for HL, but was statistically insignificant, possibly due to smaller event rates and thus limited statistical power. The only compelling deviations from this pattern were more than threefold increased risk of lymphomas in patients with IBD with concomitant PSC and in those receiving combination immunosuppressive therapy with anti-TNF-α therapy and thiopurines.

Strengths of our study include its considerable size, prospective population-based design, complete long-term follow-up, and access to drug prescription data. We also tried to accommodate the substantial methodological challenges in studies of cancer risk in patients with IBD.26 Several limitations of our study are notable: both NHL and HL are complex, heterogenous disease entities with many subgroups; these subgroups likely have at least partly different causal patterns.5 6 Analyses of specific subgroups, particularly of NHL, may therefore be required for causal inference. Unfortunately, this would limit our follow-up time, since detailed classification is only available with the introduction of ICD-O-3 in the cancer registries, which was introduced in 2005 in Sweden. Furthermore, because the baseline incidence of HL is so low, even larger studies are needed to achieve informative statistical power. One important drawback of a registry-based approach is the lack of valid information on disease severity in register data. We used bowel surgery as a proxy for disease severity, but this may have low sensitivity and misclassify subgroups of patients with severe disease who did not undergo surgery. We also used the Montreal classification for disease severity, but the positive predictive value for extensive disease (pancolitis) seemed low, at only 60.8% in the Norwegian part of the cohort (data not shown). More importantly, disease severity could not be properly adjusted for in analyses on pharmacotherapy, which entails that we cannot exclude confounding by indication. Data on prescriptions were not available in the Norwegian part of the cohort, and as such the combined analyses with data from both countries represent a mix of treated and untreated patients with IBD. Finally, our analyses on anti-TNF-α therapy may be limited in power due to low sensitivity for certain agents in the Swedish Prescribed Drug Register and the Swedish National Patient Register, which could dilute any harmful effects due to misclassification of exposure.27

An excess of intestinal, mainly B-cell NHL, is widely accepted in Crohn’s disease,28 29 while only limited data exist on overall lymphoma risk in patients with IBD. Conclusions are often hampered by low statistical power, short follow-up, conflation of different haematopoietic malignancies and based on less than a dozen observed lymphomas among patients with IBD. Since the incidence of lymphoma is low, much larger studies are needed if the true excess risk is the order of 30% as suggested by our study. Comprising a large enough dataset with sufficient data on disease severity and pharmacotherapy poses a challenge. Based on our results, we believe that patients with IBD do experience a modestly increased risk of lymphomas compared with the background. High inflammatory activity has been suggested to be the main cause of increased lymphoma risk in autoimmune diseases.4 30 Our finding of higher excess risks in patients with concomitant PSC and perianal disease is in line with this theory. But we also acknowledge that causal inference requires more detailed and validated clinical data on disease severity, lymphoma subtype and location and lifestyle factors such as smoking and body mass index. With these limitations, we could not investigate the direct effect of immunosuppressive therapy on cancer risk.

The possible impact of immunosuppressive treatment on risk of lymphomas among patients with IBD has been extensively investigated, although published studies often suffer from low statistical power due to few cases.8 21 31–33 A larger cohort study based on 336 lymphoma cases comparing exposed and unexposed patients with IBD found a 2.6-fold increased lymphoma risk following thiopurine monotherapy, 2.4-fold increased risk following anti-TNF-α therapy and about 6-fold increased risk after combination therapy.34 We detected an increased risk associated with combination pharmacotherapy and use of thiopurines, in line with previous studies.7 21 28 We could not identify any cancers in patients exposed to anti-TNF-α monotherapy; the null finding is in line with another population-based study,35 which showed no association between anti-TNF-α therapy and ‘haematopoietic and lymphoid’ malignancies (adjusted rate ratio of 0.90 (95% CI 0.41 to 1.91)), although this finding was based on only eight cases.

In conclusion, we found consistent and statistically significant associations between IBD and malignant lymphomas with little evidence of heterogeneity across different patient and disease characteristics. The absolute risk is too low to warrant screening for lymphoma in patients with IBD and should prompt caretakers to emphasise the low absolute risk to patients in any shared decision-making process concerning choice of treatment for IBD. To strengthen the base for causal inference, future studies need to assess the associations separately for at least the major subtypes of both NHL and HL. They also need to investigate in detail the impact of immunosuppressive therapies, and the potential mutual confounding effect between disease severity and drug treatment; such studies require large numbers of observed cases and long-term follow-up, and more detailed data on disease severity than what is currently available in patient registers. A final challenge is to assess whether an increased risk of malignant lymphomas in IBD is indeed causal or arises due to shared risk factors, perhaps acting through immunological mechanisms or genetic predispositions.

Data availability statement

Data are not publicly available due to their containing information that could compromise the privacy of research participants. Summary data can be made available from the corresponding author on reasonable request.

Ethics statements

Patient consent for publication

Ethics approval

This study was approved by the Regional Ethics Committee of South-East Norway (approval no. 2016/359) and the Ethics Review Board in Stockholm (approval no. 2018/1626-31/1).

References

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

  • JY and ER contributed equally.

  • Contributors JY: methodology, investigation, formal analysis, software, writing—review and editing. ER: methodology, investigation, formal analysis, software, data curation, visualisation, writing—original draft, writing—review and editing. PW: methodology, software, writing—review and editing. LMH: data curation, writing—review and editing. VP: writing—original draft, witing—review and editing. AH: writing—original draft, writing—review and editing. ML: supervision, conceptualisation, methodology, writing—review and editing. JB: supervision, conceptualisation, writing—review and editing. MB: supervision, conceptualisation, writing—review and editing. H-OA: supervision, conceptualisation, writing—original draft, writing—review and editing, funding acquisition. WY: conceptualisation, data curation, supervision, writing—review and editing, funding acquisition. MK: conceptualisation, methodology, data curation, supervision, project administration, writing—review and editing, funding acquisition, guarantor.

  • Funding This work was supported by The Swedish Cancer Society (Cancerfonden) (grant numbers 2018/703 to WY and 2018/735 to HOA); Swedish Research Council (Vetenskapsradet) (grant number 2018-02614 to H-OA); Norwegian Cancer Society (Kreftforeningen) (grant number 190345 to MK) and South-Eastern Norway Regional Health Authority (grant numbers 2015051 and 2018061 to MK). JY was partly supported by a scholarship from the China Scholarship Council.

  • Competing interests None declared.

  • 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.