Objective To evaluate exposure–response between 1,3-butadiene, styrene and lymphohaematopoietic cancers in an updated cohort of workers at six North American plants that made synthetic rubber polymers.
Methods Employees were followed from 1943 through 2009 to determine mortality outcomes. Cox regression analyses estimated rate ratios (RRs) and 95% CIs by quartile of cumulative exposure to butadiene or styrene, measured in parts per million-years (ppm-years), and exposure–response trends for all leukaemia, lymphoid leukaemia, myeloid leukaemia, acute myeloid leukaemia, non-Hodgkin’s lymphoma (NHL), multiple myeloma and all B-cell malignancies.
Results Among 21 087 workers, adjusted RRs for butadiene and all leukaemia (132 deaths) rose with increasing exposure, with an RR of 2.53 (95% CI 1.37 to 4.67) in the highest exposure quartile (≥363.64 ppm-years), and the exposure–response trend was statistically significant for all leukaemia (p=0.014) and for lymphoid leukaemia (52 deaths, p=0.007). Styrene exposure–response trends for all leukaemia and lymphoid leukaemia were less consistent than those for butadiene. Cumulative exposures to butadiene and styrene were not associated consistently with myeloid leukaemias or the B-cell malignancies, NHL and multiple myeloma.
Conclusions We confirmed a positive exposure–response relationship between butadiene and all leukaemia among workers, most of whom had coexposure to styrene. Results supported an association between butadiene and lymphoid leukaemia, but not myeloid leukaemia, and provided little evidence of any association of butadiene or styrene exposures with major subtypes of B-cell malignancies other than lymphoid leukaemia, including NHL and multiple myeloma.
- occupational health
Data availability statement
All data relevant to the study are included in the article or uploaded as supplementary information.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
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?
Workers in the synthetic rubber industry are exposed to butadiene and styrene. Our previous research on the largest cohort of synthetic rubber industry workers indicated that male workers had an excess of leukaemia that was likely to have been due to butadiene or butadiene plus styrene and other chemicals. The International Agency for Research on Cancer (IARC) and other agencies have classified butadiene as a human carcinogen causing lymphohaematopoietic cancer, especially leukaemia, and classified styrene as probably carcinogenic to humans.
What are the new findings?
The study confirmed a positive exposure–response relationship between butadiene and all leukaemia among workers with coexposure to styrene, supporting the IARC classification of butadiene as a known human carcinogen. Results supported an association between butadiene and lymphoid leukaemia, but not myeloid leukaemia. There was less support for an independent causal association between styrene and leukaemia.
How might this impact on policy or clinical practice in the foreseeable future?
In 2019, The US Environmental Protection Agency included butadiene in its list of 20 substances designated as high-priority for risk re-evaluation. Thus, butadiene will move through the process required by Toxic Substances Control Act (TSCA) to evaluate any unreasonable risk it may present to human health or the environment. The results will inform the management of this chemical under TSCA.
Synthetic rubber polymer manufacturing began in the USA and Canada in the early 1940s. Operations at the plants included in the present study have been described in detail previously.1 Styrene-butadiene rubber (SBR), a copolymer of 1,3-butadiene and styrene, initially was the main type of synthetic rubber produced. At all of the plants, SBR production areas included polymerisation, coagulation and finishing, with tank farm, laboratory, maintenance, warehouse and utilities support operations. Workers at the plants were exposed potentially to butadiene and styrene monomers and other chemicals. Monomer exposure potential varied by work area, type of job and time period. Exposure levels were relatively high in the 1940s and 1950s and declined thereafter due to changes in production, work practices and engineering controls. The International Agency for Research on Cancer (IARC),2 3 the US National Toxicology Programme4 and other agencies have classified butadiene as a human carcinogen causing lymphohaematopoietic cancer (LHC), especially leukaemia. In 2019, The US Environmental Protection Agency listed butadiene as a high-priority for risk evaluation.5 IARC has classified styrene as probably carcinogenic to humans (Group 2A), based on limited evidence of carcinogenicity in humans and sufficient evidence in experimental animals.6
This paper describes internal analyses of the exposure–response relation between cumulative exposure to butadiene and styrene and mortality from LHCs in the largest cohort of synthetic rubber polymer workers, adding to two previous papers that presented preliminary exposure–response analyses of all leukaemia, non-Hodgkin’s lymphoma (NHL) and multiple myeloma.7 8 New aspects of this paper include analyses by quartile of cumulative exposure to monomers for women and men combined, with beta coefficients (β) for the slope of the exposure–response curve and corresponding 95% CIs; results pertaining to B-cell malignancies; results for acute myeloid leukaemia (AML), reported to be associated with styrene exposure5; results for each monomer and leukaemia, stratified by exposure to the other monomer and results of analyses using lagged monomer exposure data. The present analyses provide quantitative data on exposure–response trends that could be useful for risk assessments by regulatory agencies.
Overview of study design and cohort data
Sathiakumar et al 8 described in detail the methods used to update the cohort study of mortality among workers employed at eight North American synthetic rubber polymer plants, with follow-up extended through 2009. The update included 17 924 men employed for at least 1 year and 4861 women employed for at least 1 day before 1 January 1992. Turnover among male employees in the first year of employment at the plants was large. The 1 year duration of employment criterion was imposed on the large male cohort to maximise informativeness. A similar restriction was deemed unnecessary for the smaller female cohort.
Quantitative butadiene and styrene monomer exposure estimates previously were developed for six of the eight plants by investigators who were blinded to disease outcomes.9 All analyses were restricted to 21 087 employees who had worked only at these six plants.
Work histories and monomer exposure estimates were available through the end of 1991. We did not obtain post-1991 job histories for 4079 workers who were actively employed at the end of 1991. In their last 1991 job, 46% were exposed to monomers, but their exposures were relatively low (median values, 1.1 parts per million (ppm) for butadiene and 0.4 ppm for styrene). Exposure estimation entailed identifying for each plant-specific work area/job combination its component tasks and documenting historical changes in those tasks; calculating plant-specific, work area/job-specific and time-specific average exposure indices (8-hour time-weighted average concentration in ppm) and compiling these into job-exposure matrices; and linking the time-specific and work area/job-specific exposure estimates in the job-exposure matrices with each employee’s work history to obtain cumulative exposure estimates as of each day of follow-up.
Updated vital status information through 2009 was available for 99% of the cohort.8 Data on underlying and contributing causes of death came from death certificates, the US National Death Index and the Canadian Mortality Data Base of Statistics Canada. In a prior update, we attempted to obtain medical records of men whose death certificate mentioned any type of LHC.10 Medical records were retrieved for 86% of leukaemia, 84% of NHL and 81% of multiple myeloma cases and confirmed the diagnoses of 100% of leukaemia and multiple myeloma and 96% of NHL cases with records. We did not obtain any additional medical records for the current update. Causes of death were determined without knowledge of monomer exposure status.
Monomer exposure variables
No new exposure estimates were developed for the present analysis. Exposure variables were time-dependent butadiene and styrene cumulative ppm-years. Analyses evaluated unlagged cumulative ppm-years and cumulative ppm-years lagged by 10 or 20 years to allow for potential latency. We did not analyse other exposure variables previously studied,11 including number of high-intensity tasks (HITs), average ppm and ppm-years below and above a threshold, all of which were highly correlated with ppm-years (eg, Spearman rank correlation coefficient=0.86 for butadiene ppm-years and HITs).
Outcomes were all leukaemia, lymphoid leukaemia, myeloid leukaemia, AML, NHL, multiple myeloma and all B-cell malignancies, including lymphoid leukaemia, NHL and multiple myeloma.12 The B-cell malignancy cases may have included a few T-cell neoplasms, which are thought to comprise 10%–15% of NHL and lymphoid leukaemia. For each LHC category, events included any decedent with the condition as the underlying or a contributing cause of death or with a medical record indicating that the condition was present.
Follow-up began for male workers on their date of accruing 1 year of employment or on the earliest date when complete plant records were available, whichever was later, and for female workers on their hire date or on the earliest date when complete plant records were available, whichever was later.8 For all workers, follow-up ended on the earliest of their death date, their loss-to-follow-up date or 31 December 2009.
Analyses of the relation between monomer exposure and LHC mortality used multivariable Cox regression methods to estimate LHC HRs and exposure–response trends within the cohort of workers, without reference to an external comparison population. These ‘internal’ Cox regression analyses provided maximum partial likelihood estimates of disease-specific HRs, interpreted as rate ratios (RRs) and 95% CIs for each quartile of ppm-years of monomer exposure compared with no exposure, with quartiles specified according to the exposure distribution of cases with each form of LHC. Additional Cox regression analyses estimated beta coefficients (β) and 95% CIs for trends in exposure–response using butadiene or styrene ppm-years. For all leukaemia, we further described the exposure–response curve with restricted cubic spline (RCS) Cox regression models, fitted to all exposure data and to the trimmed data, with five knots corresponding to the 5, 27.5, 50, 72.5 and 95 percentile boundaries among the exposed. Also, for all leukaemia, we analysed butadiene exposure–response separately in the following two strata of styrene exposure: lower styrene exposure, defined as below the median value of 27 styrene ppm-years among leukaemia decedents; and higher styrene exposure, defined as exposure at or above the median value. We also analysed styrene exposure–response, stratified by lower (below the median value, 121 ppm-years) versus higher (≥121 ppm-years) butadiene exposure.
All analyses used age as of each person-day of follow-up as the time scale and treated monomer ppm-years as time-dependent. Our main analyses assessed exposure–response trends using all person-day records, and models included as covariates, age at hire, calendar year of hire, sex, race, plant and payroll status (ever hourly paid or always salaried). Salaried employees mainly held supervisory or managerial positions, while hourly employees typically worked directly and regularly with manufacturing operations and had higher potential exposure to monomers. Payroll status also was associated with socioeconomic factors.
We conducted several series of sensitivity analyses. One of these excluded person-day records having zero cumulative exposure, in order to eliminate the possibility that any observed exposure–response trend was due to differences in uncontrolled factors between unexposed and exposed person-time.13
To investigate the influence of data at extreme exposure values, additional sensitivity analyses estimated exposure–response trends using ‘trimmed’ data that excluded all unexposed person-time and all person-time with ppm-years values above the 95th percentile of the exposure distribution of leukaemia decedents (for analyses of all leukaemias, lymphoid leukaemia, myeloid leukaemia and AML) or of B-cell malignancy decedents (for analyses of NHL, multiple myeloma and all B-cell malignancies). Several considerations prompted these analyses. Cohort studies in other industry settings have reported that exposure–response curves tend to diminish at higher exposure levels.14 Two of our earlier studies of male synthetic rubber polymer workers found stronger exposure–response trends for butadiene and leukaemia in analyses that excluded exposures above the 95th percentile11 or categorised butadiene into deciles.7 Both of the latter procedures can reduce the impact of exposure outliers. In addition, an investigation at the largest study plant, performed to validate our butadiene exposure estimates, found greater misclassification for jobs entailing higher exposures than for jobs with lower exposures.15
Additional sensitivity analyses used ‘reduced’ models that contained fewer covariates. The goal of the reduced models was to preserve the control of confounding, while providing more precise results. Thus, we examined models containing all possible subsets of covariates and selected a reduced model on the basis of having: (1) a monomer exposure parameter estimate within 5% of that obtained in the corresponding full model and (2) the fewest covariates.
Preliminary models evaluated possible differences in monomer exposure-LHC associations between women and men. Those analyses found no statistically significant sex-monomer exposure interaction and are not described further.
We used the SAS V.9.4 Cox proportional hazard model procedure Proportional Hazards REGression (PHREG) for the Cox regression analyses. We used the Akaike information criterion (AIC) to compare the statistical fit of reduced versus full models.16
Of the 21 087 workers in the cohort, 14 004 (66%) were classified as ever exposed to butadiene, and 15 422 (73%) were ever exposed to styrene (table 1). Cumulative exposure to monomers was right-skewed: overall, exposed workers had median and mean values of 48 and 187, respectively, for butadiene ppm-years and of 11 and 38, respectively, for styrene ppm-years. Monomer exposure was higher among men than among women and also varied by plant, payroll status, race and period of hire. The cohort had a median of 8.3 years of employment at the end of 1991. At the end of the follow-up, the median time since hire was 40 years, the median age was 69 years, and 46% of the cohort was deceased.
Total numbers of decedents with an outcome event were 132 for all leukaemia, 52 for lymphoid leukaemia, 67 for myeloid leukaemia, 41 for AML, 110 for NHL, 60 for multiple myeloma and 213 for B-cell malignancy. Butadiene and styrene ppm-years were strongly correlated, with Spearman correlation coefficients ranging from 0.81 for NHL decedents to 0.87 for all leukaemia and lymphoid leukaemia decedents. Online supplemental table S1 shows decedents with each outcome, cross-classified by quartiles of butadiene and styrene ppm-years.
For butadiene, RCS analyses (figure 1A,B and online supplemental figure S1A) indicated that the adjusted RR for all leukaemia increased in an approximately linear fashion at exposures below about 1000 ppm-years, with attenuation of the curve at higher exposures. Analyses by quartile of ppm-years also indicated that adjusted RRs rose with increasing exposure, with an RR of 2.53 (95% CI 1.37 to 4.67) in the highest exposure category (≥363.64 ppm-years) (table 2). Using untrimmed butadiene ppm-years, the exposure–response trend was statistically significant, regardless of the inclusion of unexposed person-time (including unexposed: β=2.55×10−4, (95% CI 0.52 to 4.57)×10−4, trend p=0.014) (excluding unexposed: β=2.50×10−4, (95% CI 0.27 to 4.73)×10−4, trend p=0.028). Trimming to restrict data to ppm-years >0 and ≤95th percentile (1144 ppm-years) of all leukaemia decedents yielded a somewhat stronger exposure–response trend for butadiene (β=9.94×10−4, (95% CI 1.88 to 18.00)×10−4, trend p=0.016). This result was consistent with RRs from the analysis by exposure quartile. For example, the beta coefficient obtained from the trimmed data, evaluated at the mean (1156.83 ppm-years) or median (725.08 ppm-years) value of exposure in quartile 4, yielded RRs of 3.16 and 2.06, respectively, as compared the RR of 2.53 from the analysis by quartile. In contrast, the beta coefficient from the analysis of untrimmed butadiene data, evaluated at the same mean and median ppm-years, yielded RRs of 1.34 and 1.20, respectively.
For styrene, RCS results (figure 1C,D and online supplemental figure S1B) indicated a positive exposure–response curve, with attenuation at higher exposures. Other analyses (table 2) yielded results similar to those for butadiene ppm-years, but those results were not statistically significant.
Analyses of each monomer, stratified by lower vs higher exposure to the other monomer, found that the butadiene exposure–response association was evident only in the category of higher styrene exposure (β=2.89×10−4, (95% CI 0.41 to 5.36)×10−4, trend p=0.022) (table 2). Styrene exposure–response was not observed in either category of butadiene exposure.
Lymphoid and myeloid leukemia
For lymphoid leukaemia (table 3), the pattern of RRs by quartile of butadiene ppm-years was irregular, with RRs of 2.61 (95% CI 1.02 to 6.67) and 1.95 (95% CI 0.76 to 5.03) in the third (213.43–<376.31 ppm-years) and fourth (≥376.31 ppm-years) quartiles, respectively; exposure–response trends were statistically significant (including unexposed: β=3.81×10−4, (95% CI 1.05 to 6.58)×10−4, trend p=0.007; excluding unexposed: β=4.78×10−4, (95% CI 1.34 to 8.23)×10−4, trend p=0.007). For all myeloid leukaemia and AML, the RR for each butadiene exposure quartile was elevated but statistically imprecise, and none of the exposure–response trends was significant. Analysis of trimmed exposure data yielded a statistically significant exposure–response trend for butadiene and lymphoid leukaemia (β=15.40×10−4, (95% CI 4.19 to 26.53)×10−4, trend p=0.007) but not for all myeloid leukaemia or AML.
For styrene, analyses by quartile of ppm-years yielded results similar to those for butadiene ppm-years for lymphoid leukaemia, all myeloid leukaemia and AML (table 3), but the only statistically significant RR was that for AML in the third exposure quartile. The only statistically significant exposure–response trend was that for lymphoid leukaemia, with the inclusion of the unexposed (trend p=0.046).
NHL, multiple myeloma and B-cell malignancy
For butadiene and NHL (table 3), no exposure–response was detected in analyses of exposure quartile or trends using untrimmed butadiene ppm-years. However, trimming to restrict data to ppm-years >0 and ≤95th percentile (1083 ppm-years) yielded a trend p value of 0.002. No statistically significant butadiene exposure–response was found for multiple myeloma. For butadiene and all B-cell malignancies combined, the trend was statistically significant (p=0.002) only in analyses of trimmed exposure data. Styrene did not appear to be associated with NHL, multiple myeloma or all B-cell malignancy.
Analyses that lagged cumulative exposure by 10 or 20 years (table 4) found that for all leukaemia and lymphoid leukaemia, butadiene exposure–response beta coefficients were similar regardless of the exposure lag applied. For all leukaemia, trends were not statistically significant for styrene exposure lagged 10 or 20 years. For lymphoid leukaemia, only the trend for styrene lagged 10 years approached statistical significance. No trends were detected for lagged butadiene or styrene exposure and other LHCs.
Seventy-eight per cent of reduced models yielded exposure parameter estimates that were within 5% of those in full models, and all were within 10% online supplemental tables S2–S5. Use of the reduced models did not identify any additional statistically significant results.
The present cohort comprises the largest butadiene-exposed group studied to date and includes the only large group of butadiene-exposed women. The update8 on which the current analyses were based substantially augmented information from our previous studies of male1 17 and female18 synthetic rubber polymer workers.
Internal Cox regression analyses of all leukaemia found a statistically significant exposure–response trend for butadiene but not for styrene. Analyses of major subtypes of leukaemia indicated that for lymphoid leukaemia a trend of increasing risk with increasing cumulative exposure to butadiene was present, with evidence of a similar pattern for styrene that, however, was not statistically significant in most analyses. For all myeloid leukaemia and AML a positive exposure–response trend was not evident for either monomer. Restricting data to exposed person-time did not substantively alter these results. For all leukaemias and lymphoid leukaemia, the butadiene trends were three to four times stronger when data were restricted to ppm-years >0 and ≤95th percentile (1144 ppm-years). Alternative explanations for this pattern include biological irrelevance of higher exposures, exposure misclassification, depletion of the pool of susceptible people and other possible reasons.14 Lagging exposure had little effect on results for leukaemias or other outcomes, as exposure diminished over calendar time,9 more than 90% of the LHC cases died 20 years or more after hire, and no exposure was recorded after 1991.
With regard to butadiene and leukaemia, the results of the present study are consistent with those reported earlier by other investigators19–22 and us1 7 8 10 11 17 18 23 24 in supporting a positive association in the synthetic rubber polymer industry. In contrast, investigations of relatively small cohorts of butadiene monomer production workers, who were exposed to butadiene but not to styrene, have reported results for leukaemia that were null or weakly positive.25–28
With regard to styrene, increased risks of LHC, particularly leukaemia and lymphoma, have been reported among styrene-exposed workers in both synthetic rubber polymer and reinforced plastic industries.6 Reinforced plastics industry workers were exposed to styrene concentrations at higher levels than were typically found in the synthetic rubber industry and were not exposed to butadiene. Positive exposure–response relationships between styrene exposure and leukaemia were reported in an older multinational European study of reinforced plastics workers29 but not in a recent reanalysis of that study.30 A study of mortality among workers in the reinforced plastics boatbuilding industry reported an association between several indices of exposure to styrene and leukaemia.31–33 However, studies of several other cohorts of reinforced plastics workers found little evidence of an association between styrene and overall leukaemia.34–37
With regard to exposure to butadiene or styrene and cell-type specific leukaemias, data from other studies are limited.6 Christensen et al 36 reported a positive exposure–response trend for a styrene cumulative exposure score and AML. Collins et al 35 found no statistically significant exposure–response trend among workers employed in the US reinforced plastics and composite industry.
NHL, multiple myeloma and B-cell malignancy
We found little evidence of any association between cumulative exposure to butadiene or styrene and NHL or multiple myeloma. These results are consistent with our previous analyses.7 8 The positive exposure–response for butadiene and all B-cell malignancies combined, seen only in analyses using trimmed exposure data, may have reflected the association between butadiene and lymphoid leukaemia.
Several studies have assessed mortality from NHL or subtypes of NHL in three cohorts of butadiene monomer production workers exposed to butadiene but not to styrene.25–28 The observed number of NHL deaths exceeded the expected number in each study, but all results were based on small numbers.
An investigation that included many of the subjects in our study reported a positive association between styrene and NHL.21 The reanalysis of the European reinforced plastics cohort reported an association between mean level of exposure to styrene and NHL but found no association with cumulative exposure to styrene.30 In contrast, updated studies of the reinforced plastics industry did not find a clear excess of NHL, overall or in subgroups with higher styrene exposure.31 32 34–36 38
Matanoski et al 21 reported that multiple myeloma was associated with butadiene in a study that included most of the subjects in our male cohort, in contrast to our analyses by cumulative exposure, which provided no support for an association with butadiene or styrene. Divine and Hartman25 reported slightly more than expected deaths from multiple myeloma among butadiene production workers, but data on multiple myeloma were sparse and internally inconsistent.
Strengths and limitations
Study strengths included the long follow-up period, the use of objective procedures to classify workers according to monomer exposure and cause of death, the inclusion of female employees and the use of sensitivity analyses to facilitate interpretation of results. Potential misclassification of monomer exposure remains a concern. Previous investigations have validated the butadiene exposure estimates15 and assessed the impact of exposure misclassification on the exposure–response relation between butadiene and leukaemia.1 39 However, sensitivity analyses to assess the potential impact of monomer exposure uncertainties would strengthen the data if used for risk assessment.
Other limitations include lack of information on lifestyle factors and the use of mortality, rather than incidence, data to ascertain LHCs. Mortality data are not optimal for cancers with relatively long survival including certain leukaemias and B-cell malignancies. Although our results provided little support for the hypothesis that styrene causes LHC, the styrene concentrations experienced by this cohort were relatively low, as compared with the exposures of reinforced plastics industry workers.35 We performed multiple comparisons, and it is possible that the observed associations occurred by chance.
This study confirmed a positive exposure–response relationship between butadiene and all leukaemia and supports the classification of butadiene as a human carcinogen. Results supported an association between butadiene and lymphoid leukaemia, but not myeloid leukaemia. Evidence of an independent causal association between styrene and leukaemia remained less convincing. The study found little evidence that butadiene or styrene exposures were associated with major subtypes of B-cell malignancy other than lymphoid leukaemia, including NHL and multiple myeloma.
Data availability statement
All data relevant to the study are included in the article or uploaded as supplementary information.
Patient consent for publication
The Institutional Review Board of the University of Alabama at Birmingham reviewed and approved the study protocol.
We are grateful to the late James McGraw of the International Institute of Synthetic Rubber Producers for his effective coordination of many sponsor–investigator interactions since the inception of the original study in 1991. Maurizio Macaluso and Rodney Larson carried out exposure estimation in previous investigations.
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.
Contributors NS and ED contributed substantially to the conception and design of the work. IB, LC, ML and TA analysed the data. NS, ED, IB, LC, ML, MT, BEB and TA interpreted the data. NS, ED, IB, MT and BEB drafted the work and NS, ED, IB, LC, ML, MT and BEB revised it critically for important intellectual content. NS, ED, IB, LC, ML, MT, BEB and TA approved of the final version published. NS, ED, IB, LC, ML, MT, BEB and TA agreed to be accountable for all aspects of the work to ensure questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Funding International Institute of Synthetic Rubber Producers, American Chemistry Council (Olefins Panel) and Styrene Information and Research Center. The sponsors were given an opportunity to provide comments on a draft of this paper. However, the contract between the University of Alabama at Birmingham and the sponsors stipulated that the academic investigators should independently carry out the design, conduct and reporting of the study. Accordingly, the authors made all decisions about the contents of this paper.
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.