Article Text

Download PDFPDF

Do smokers have greater risk of delayed and non-union after fracture, osteotomy and arthrodesis? A systematic review with meta-analysis
  1. R G Pearson,
  2. R G E Clement,
  3. K L Edwards,
  4. B E Scammell
  1. Division of Rheumatology, Orthopaedics and Dermatology, School of Medicine, University of Nottingham, Queen's Medical Centre, Nottingham, UK
  1. Correspondence to Dr R G Pearson; richard.pearson{at}nottingham.ac.uk

Abstract

Objectives Systematic review and meta-analysis of published observational cohort studies. To quantify the increased risk smokers have of experiencing a delayed and/or non-union in fractures, spinal fusion, osteotomy, arthrodesis or established non-unions.

Setting Medical Literature Analysis and Retrieval System Online (MEDLINE), Excerpta Medica database (EMBASE), Allied and Complementary Medicine Database (AMED) and Web of Science Core Collection from 1966 to 2015.

Study eligibility criteria, participants and interventions Observational cohort studies that reported adult smokers and non-smokers with delayed and/or non-union or time to union of the fracture, spinal fusion, osteotomy, arthrodesis or established non-union were eligible.

Data extraction and outcome measures 2 authors screen titles, abstracts and full papers. Data were extracted by 1 author and checked independently by a second. The relative risk ratios of smoking versus non-smoking and the mean difference in time to union patients developing a delayed and/or non-union were calculated.

Results The search identified 3013 articles; of which, 40 studies were included. The meta-analysis of 7516 procedures revealed that smoking is linked to an increased risk of delayed and/or non-union. When considered collectively, smokers have 2.2 (1.9 to 2.6) times the risk of experiencing delayed and/or non-union. In all the subgroups, the increased risk was always ≥1.6 times that of non-smokers. In the patients where union did occur, it was a longer process in the smokers. The data from 923 procedures were included and revealed an increase in time to union of 27.7 days (14.2 to 41.3).

Conclusions Smokers have twice the risk of experiencing a non-union after fracture, spinal fusion, osteotomy, arthrodesis or treatment of non-union. Time to union following fracture, osteotomy, arthrodesis or treatment of an established non-union is longer in smokers. Smokers should be encouraged to abstain from smoking to improve the outcome of these orthopaedic treatments.

  • Non-union
  • Delayed-uion
  • Fracture
  • Smoking
  • meta-analysis

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 and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/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.

Strengths and limitations of this study

  • This is a comprehensive systematic review of the published literature, including 40 studies achieving a large sample size of over 8000 adults.

  • This meta-analysis identifies that smokers take 27.7 days (14.2 to 41.3) longer for union to occur and smokers have double the risk of non-union 2.2 (1.9 to 2.6).

  • Included data originated from observational cohort studies as randomised controlled trials (RCTs) are not available addressing this topic. The feasibility of conducting RCTs with patients being randomly allocated to continuing smoking or smoking cessation is faced with many ethical issues and potential biasing. The quality of included studies as defined by the Meta-Analyses and Systematic Reviews of Observational Studies (MOOSE) guidelines was variable and notably often did not adjust for baseline differences. Similarly, included studies often did not thoroughly investigate potential confounders such as alcohol consumption.

Introduction

The incidence of fractures in England is reported to be as high as 3.6% with over a third of the population sustaining at least one fracture in their lifetime.1 The majority of fractures heal uneventfully, but those that ultimately fail to unite have wide-ranging implications for the patient and are a huge financial burden on the healthcare service. The cost of non-union in a tibia shaft fracture can be upwards of £10 000.2

A number of clinical factors are known to be deleterious to fracture healing, including the fracture classification severity, advancing age, infection, the use of non-steroidal anti-inflammatory drugs and corticosteroids and smoking.3 Despite the considerable morbidity and mortality associated with smoking, it is estimated that 19% of the adult UK population continues to smoke.4 Smoking is therefore an important modifiable factor associated with delayed healing of fracture and allied orthopaedic procedures.

This research focuses on the systematic review and meta-analysis of observational cohort studies. These were considered collectively but also as subgroups based on the injury or the treatment. The subgroups for the delayed and/or non-union in smokers and non-smokers were long bone fracture,5–18 spinal fusions,19–24 osteotomy,25–29 arthrodesis30–37 and established non-union.23 ,38–42 The mean difference in the time to union between smokers and non-smokers was also subjected to meta-analysis.5 ,7 ,13 ,14 ,26 ,27 ,43 ,44

There have been a large number of observational studies investigating the effect of smoking on bone healing, and these have been summarised in a number of recent review articles.45–50 The objective of this review was to identify if there is an increase in delayed and/or non-union in smokers compared with non-smokers and if this is the case to quantify this as risk ratios. In addition, we quantified the mean difference in the time to union between smokers and non-smokers. Hence, we are able to provide evidence-based information on the difference in risk of developing a non-union or the time taken for union to occur in smokers.

Methods

We followed the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines in conducting a systematic review of the literature (research checklist 1). In reporting our meta-analysis, we followed the guidance provided by MOOSE; proposals for reporting meta-analysis of observational studies (research checklist 2).51

Data sources and search strategy

A literature search was conducted through MEDLINE, EMBASE, AMED, Web of Science Core Collection from 1966 to May 2015 using the keywords ‘smoking’, ‘nicotine’, ‘fracture’, ‘arthrodesis’ (and the alternative spelling ‘arthrodeses’) and osteotomy as subject headings and abstract, using the search string (smoking OR nicotine) AND (fracture OR arthrodesis OR arthrodeses OR osteotomy). Non-English language articles were included. An example of a detailed search strategy for Web of Science is provided (see online supplementary file 1). The search results were exported to an Endnote database. Duplicates were deleted leaving us with the potentially relevant studies. We narrowed down the papers through a staged process. In stage 1, we reviewed the titles and if the title indicated likelihood that the article might meet the eligibility criteria, we retrieved the abstract. In stage 2, we reviewed the abstracts and selected potentially eligible studies for retrieval of the full text. In stage 3, we checked the full-text publication against the eligibility criteria. Finally, we reviewed the bibliography of each study that met our inclusion criteria for additional studies. Conference abstracts were excluded.

Study selection

Inclusion criteria

We identified studies presenting data which enabled the calculation of the risk ratio of delayed and/or non-union between smokers and non-smokers. Similarly, through systematic review, we identified observational cohort studies reporting or providing raw data which enabled the calculation of the mean difference in time to union between smokers and non-smokers. Only peer-reviewed published literature of adults who had fractured a bone or who had undergone a spinal fusion, osteotomy or arthrodesis or treatment of an established non-union was included. Occasionally, instances arose where a small amount of additional data was required to include the study and the corresponding author was contacted. Where data were forthcoming in a personal communication, the study was included.

Exclusion criteria

Papers that did not make a direct comparison between smokers and non-smokers were excluded. Limb distraction studies were not included unless this was conducted as a combined procedure with arthrodesis.

Data extraction and quality assessment

Data extraction was independently performed by two authors (RGP, RGEC). The following data were extracted from each paper: author's name, publication year, study design, patient demographics, definition of smoker, fracture/procedure type, definition of union and non-union, length of follow-up, loss to follow-up, time to fracture union, number of patients included in the study and number experiencing delayed and/or non-union. The subgroup meta-analysis was conducted on fractures, spinal fusion, osteotomy, arthrodesis and established non-union were recorded.

The quality of each selected study was assessed with the Newcastle–Ottawa Scale (NOS) for assessing the quality of non-randomised studies in meta-analyses.52 This checklist was applied to make an assessment covering three broad perspectives: the selection of the study groups; the comparability of the groups; and the ascertainment of either the exposure or outcome of interest in each observational cohort study. The studies are graded 0–9 with higher scores indicating higher quality studies. κ analysis was performed to assess agreement between the quality raters.

Data synthesis and analysis

For statistical analysis and the preparation of figures, several software packages were used; IBM SPSS Statistics V.20, StatsDirect statistical software V.2.8 and Review Manager (RevMan) software, V.5.3 (Copenhagen: The Nordic Cochrane Centre, The Cochrane Collaboration, 2012) were used. Raw data for non-union events and time to union from each study were entered and used to calculate a relative risk (RR) or mean difference for each study. The subtotal and total pooled risk ratios were obtained using a random effects model based on estimating the amount of between-study variation by comparing each study's result with a Mantel-Haenszel fixed-effect meta-analysis result. The subgroups of fractures, spinal fusion, osteotomy, arthrodesis and established non-union were identified. The analysis also attributed a weight to each study defined by inverse variance. Forest plots were used to express risk ratios or mean differences with 95% CI for each study. Hence, the main outcome measures were the risk ratio between smokers and non-smokers in the development of delayed union and/or non-union. The mean difference in the time to union between smokers and non-smokers was also determined.

Heterogeneity between studies and assessment of meta-analysis bias

The Cochrane Q and the I2 statistics were used to assess study heterogeneity. Bias was initially assessed by funnel plot asymmetry. The data were also subjected to Egger analysis. Where bias was identified through this methodology, it was investigated further, considering parameters such as size of study, bone site or confounders such as diabetes, excessive alcohol consumption, drug abuse or socioeconomic group. Calculation of Rosenthal's fail-safe number was considered as suggested in the MOOSE guidelines but was discounted due to the reasons reported in the Cochrane guidance on meta-analysis.53 Essentially, it has been reported that available calculation methods lead to diverse estimates of the number of additional studies required to bring the p value above 0.05. The method also runs against the principle that in medical systematic reviews, the focus is on effect size and CI of the estimated observed effect as opposed to the p value reaching an arbitrary value.

Results

Identification of studies

The search identified 1164 articles in MEDLINE, 2613 in EMBASE, 35 in AMED and 1991 in Web of Science Core Collection at stage 1. At stage 2, duplicates were deleted leaving 3013 potentially relevant publications. At stage 3, selection was made from the 3013 based on the titles, which identified 396 publications. At stage 4, selection was based on reading the abstracts, giving a total of 101 publications. At stage 5, from studying the 101 full-text publications, and 46 study citation arising from these manuscripts (15 of which were included), a total of 40 studies were included and used for the meta-analysis (figure 1).

Figure 1

PRISMA flow diagram describing numbers of included and excluded studies at the progressive stages of searching and reviewing the literature.

Characteristics of included studies

Of the 3013 publications, we extracted data from 40 studies that met our eligibility criteria (see online supplementary file 2). These studies had been conducted in a broad range of clinical settings, ranging from provincial clinics to level 1 trauma centres within university teaching hospitals, with some studies being multicentre. The range of nations where the research originated was diverse, including the Netherlands, Finland, Australia, Sweden, USA, Germany, Denmark, West Indies, China, Belgium, Greece, Turkey, Italy, Canada and the UK. Classification of the extent a smoker smoked was broad across of the range of studies which meant that it was only possible to classify participants into smokers and non-smokers for the meta-analysis. Exclusion of studies from the meta-analysis covered both types of case–control studies, each with their own complexities. Those where the case was smoking54 ,55 and those where the case was non-union.56–59 Within the included studies, there were retrospective and prospective observational cohort studies, with the majority being retrospective. This often gave rise to confusion regarding the reporting of exclusion criteria and loss to follow-up. Many studies chose to exclude patients due to a combination of reasons with those being excluded on the basis of a lack of suitable follow-up not being specifically identified. Inconsistent definition of delayed and non-union was prevalent in the studies with several grouping their defined definition of delayed with their definition of non-union. Hence, we reached the decision not to discriminate between delayed and non-union in the meta-analysis. In general, the follow-up was for a sufficient length of time and numbers lost <20%; however, on occasion, this information was not reported. The definition of union tended to be a combination of radiological and clinical factors, often based on weight-bearing capacity of a lower limb. Radiological union criteria tended to be described as bridging of at least three of the four bone cortices. Comorbidities were mentioned in a number of studies; however, it was rare to identify an approach to combat such confounders using multivariate analysis. The quality of the eligible studies was broad and often lacking in detail, which made assessment of the quality using the NOS difficult to apply (see online supplementary file 2). This was highlighted by the κ analysis indicating relatively poor agreement, κ=0.18 (0.05 to 0.31) and weighted κ 0.24 (0.09 to 0.39).

Delayed and/or non-union and time to union

The analysis of the data reveals that smokers have a greater risk of incidence of delayed and/or non-union. This is presented as risk ratios in the forest plot (figure 2). Thirty-eight studies were included. The patients underwent a broad range of orthopaedic treatments and procedures from fractures, spinal fusions, osteotomies, arthrodesis to treatment of established non-union. The meta-analysis of 7500 procedures revealed that smokers have an increased risk of delayed and/or non-union. When considered collectively, smokers have 2.2 (1.9 to 2.6) times the risk of experiencing delayed and/or non-union. In all the subgroups, the increased risk was always ≥1.6 times that of non-smokers.

Figure 2

Delayed and non-union in smokers and non-smokers—fractures, spinal fusions, osteotomy, arthrodesis and treatment of established non-union.

Time to union was a longer process in smokers as identified in the eight studies included in the mean difference to union meta-analysis. Six of these studies also provided data for the RR meta-analysis. Figure 3 shows that the bones of smokers take longer to unite following tibia fractures, including when treated with an external-fixator; or tibial osteotomies, or established non-union of the tibia. This was also identified when treating hallux valgus with an osteotomy or arthrodesis. The data from 923 procedures revealed it took an additional 27.7 days (14.2 to 41.3) for union to occur in smokers.

Figure 3

Time to union in smokers and non-smokers—fractures, spinal fusions, osteotomy, arthrodesis and treatment of established non-union.

Heterogeneity and bias

Considering all the included studies as one group with df=37 in the delayed and/or non-union data, there is moderate heterogeneity, with an I2 value of 47%. This heterogeneity was further investigated, the majority of the diversity originating from studies where the number of events, that is, delayed and or non-union, were small and in particular where there was a zero value in one of the exposures, usually non-smokers. Cochrane RevMan software assigns a very small value to enable the risk ratio to be calculated. It is worth noting that this then results in the large risk ratio estimates being made, with large CIs. When there are a large number of studies in the subgroup, heterogeneity is not observed within that subgroup, with I2 values of 0% as seen in the fractures subgroup 1.7 (1.4 to 2.0). There were fewer studies in each of the remaining subgroups, however in spinal fusions 1.6 (1.3 to 2.1), osteotomy 2.9 (1.9 to 4.4) and arthrodesis 3.0 (2.1 to 4.3) no heterogeneity was identified, I2 equal to 0%. In the treatment of the established non-union subgroup 11.8 (0.8, 178.8), I2 was equal to 23%. This heterogeneity can be attributed to the Hak et al study,38 a relatively small study that did not identify any non-unions after surgical treatment of the non-smoking patients. All of these three procedures are linked to a higher risk ratio. It could be considered that arthrodesis and treatment of non-union are likely to be linked to comorbidities which impact on bone biology and hence the increased risk of incidence of non-union in smokers. Alternatively, these procedures are less likely to result in a non-union in non-smokers which could increase the risk ratio with smokers. The heterogeneity between the subgroups is considerable and is supported by an I2 of 84%, clearly being attributable to the difference between the fractures/spinal fusion subgroups and osteotomy, arthrodesis and treatment of non-union.

When considering the bias across the studies within the non-union analysis, there was little identified, exemplified in the SE—log RR funnel plot, the majority of studies lying within the demarcated funnel (figure 4). However, when considering the bias among the time to union studies, it is evident that those studies with larger SE tended to be the studies which identified the largest effect size between the smokers and non-smokers.

Figure 4

Bias funnel plots; (A) non-union log relative risk, (B) time to union effect size.

Discussion

Smokers have poor bone healing. The meta-analysis reveals that smokers have a greater risk of incidence of delayed and/or non-union. This is the case for a broad range of orthopaedic treatments and procedures from fractures, spinal fusions, osteotomies, arthrodesis to established non-union. This meta-analysis assessed a total of 7516 procedures. When considered collectively, smokers have 2.2 (1.9 to 2.6) times the risk of experiencing delayed and/or non-union. In all the subgroups, the increased risk was always ≥1.6 times that of non-smokers.

The eight eligible time to union studies identified that union was a longer process in smokers and is reported as the mean difference in time to union in the meta-analysis. A total of 923 procedures comprised the data set and revealed time to union was 27.7 days (14.2 to 41.3) longer in smokers. Bones of smokers take longer to unite; following tibia fractures, tibial osteotomies or established non-union. This was also the case for osteotomy or arthrodesis used in the treatment of hallux valgus.

Contribution of smoking and injury to healing

Cigarette smoke contains in excess of 4800 constituents with at least 200 known to be toxic. It is well documented that smoking has an adverse effect on bone physiology, leading to decreased bone mineral density and increased incidence of osteoporotic fractures.60–62 The effect of smoking on fracture healing is considered to be multifactorial; physiological changes at the fracture site include hypoxia and modification of cellular metabolic activity. Nicotine is a powerful vasoconstrictor causing reduction in peripheral blood flow; carbon monoxide reduces the oxygen-carrying capacity of the blood through the formation of carboxyhaemoglobin and hydrogen cyanide inhibits cytochrome c oxidase preventing aerobic metabolism.

We hypothesise that a significant contribution to the deleterious effects of smoking is made through changes to the blood supply inducing hypoxia, implying the effect could be greater in severe fracture subtypes where oxygen delivery to the fracture site is most extensively compromised. The accelerated failure time model constructed by Kyro et al5 implies that as age and injury severity increases, the effect of cigarette smoking becomes more important. This is supported by evidence from Schmitz et al6 who found healing times in smokers were 1.5 times longer than non-smokers following simple fractures but at least 2.5 times longer after complex fractures where the damage to the blood supply is likely to be more significant. Adams et al7 showed a trend for increasing times to union in open tibia fractures from Gustilo grade I to IIIA, but the trend did not continue into IIIB subtypes. The effect of smoking is likely to be linked to smoking-related comorbidities, for example, chronic obstructive pulmonary disease.63–65 There are also considerations regarding secondhand smoke; research has established that there is a link with biomarkers of cardiovascular disease risk in passive smokers.66

An interesting question is whether smokers develop a larger zone of injury than non-smokers which could be mediated by hypoxia and hyperoxidant stress. We are aware that smoking increases free radical production and that a high yield of free radicals is produced when a bone fractures.67 Hyperoxidant stress may play a larger part in smokers following fracture where their natural antioxidant mechanisms are already depleted and could contribute to why smokers’ bones healed more slowly.5

Strengths and weaknesses

This is a comprehensive systematic review of the published literature presenting data synthesised through meta-analysis. This review adds to the literature of smoking and non-union. To the best of our knowledge, Scolaro et al48 is the only published systematic review on this topic which has implemented meta-analysis which addressed non-union in long bones OR 2.32 (1.76 to 3.06) with subgroup analysis of tibial fractures and open fractures as subgroup analysis; whereas here we report a broad grouping of orthopaedic procedures with subgroups of fractures, spinal fusions, osteotomy, arthrodesis and non-union. Hence, our research indicates that smoking impacts a broad range of elected orthopaedic surgical procedures in addition to primary long bone fracture treatment.

However, there was a wide range in quality of studies included in the analysis from which we have drawn our conclusions. There was on the whole very little information regarding how many cigarettes were being smoked and it is also of note that it was extremely rare for any biochemical confirmation of non-smoking status, that is, cotinine or carbon monoxide.68 Study quality, severity of fracture or surgical procedure, age of patients, bone involved and non-union definition were identified as potential sources of heterogeneity. We have noticed that in all the studies where surgery for non-union is discussed, the percentage of smokers requiring surgical intervention is always higher than the non-smokers, but none of these studies have ascertained a statistical significance for this association.5 ,7 ,26 In the light of these reports, we should consider this as a potential source of bias which merits further investigation. An additional topic which warrants further investigation is the covariate in smokers of alcohol consumption.

Comparison of our univariate analysis with published RR or OR

Some studies reported risk ratio, whereas others reported ORs along with the raw data we extracted for meta-analysis. In four publications, OR and RR had been generated using multivariate regression analysis as opposed to univariate methods With an aim to facilitate comparison between our univariate analysis, employed for our meta-analysis, and published ORs, we calculated OR in addition to RR in these instances (table 1).

Table 1

Comparison of published odds and risk ratios for included studies in comparison with the values generated through univariate meta-analysis

Our OR and RR values tended to be similar or identical to published univariate values with the exception of Chang et al39 OR 5.4 (1.3 to 22.13) where there was a small discrepancy compared with our OR of 4.7 (1.2 to 18.4) (table 1). On investigation, the reason for this could not be identified.

The published multivariate models are likely to have been adjusted for baseline differences such as age and sex, and in some studies for comorbidities, for example, severity of fracture. Considering the published multivariate analysis, W-Dahl and Toksvig-Larsen26 produced a similar RR using multiple regression analysis to our univariate value. However, in two publications where multivariate models had been constructed, the OR reported was approximately double our univariate value;15 ,17 in a third comparison, the published OR was approaching 50% of the univariate value.42

Two studies published ORs without raw data and therefore could not be included in the meta-analysis. Metsemakers et al69 reported on 486 tibial fractures using multiple linear regression smokers OR 1.74 (0.87 to 3.49) for delayed union and 0.96 (0.48 to 1.95) for non-union; however, in neither delayed or non-union was smoking considered to effect the odds statistically, p=0.120 and 0.915, respectively. We postulate that this could be due to the relatively young study population (mean 30.2 years). Similarly, Egol et al70 presented an OR without suitable data for inclusion in the meta-analysis for long bone non-union in smokers 1.3 (0.27 to 6.40) also stating a lack of significance (p>0.05).

Time to union including limitations arising due to primary data

The meta-analysis for the time to union was based on the mean time to union. This is not ideal when dealing with censorship when the union has not occurred during the study due to the length of follow-up. This was reported to have occurred with two patients in one study, where a value of 52 weeks was attributed at the final follow-up, despite healing not having occurred.44 This can be considered to have had a minimal influence in the total of 923 patients in the meta-analysis as censoring was not reported in any other study. An alternative approach to dealing with censoring in non-union data is survivorship modelling such as Kaplan-Meier estimates which provides an estimated median time to union.6 The Cox HRs have also been considered for discriminating between time to union in smokers and non-smokers.10

Smoking cessation and advice to patients

Perioperative smoking cessation is considered to be beneficial for fracture healing, although there are limited data to draw this conclusion.20 ,21 ,71 As a large number of orthopaedic procedures are not elective, postoperative smoking is likely to play a more significant role in orthopaedic clinical practice. Andersen et al showed that fusion rate was higher in patients who gave up smoking postoperatively and Glassman et al observed a trend towards significance in those who gave up for longer than 6 months.20 ,21 It has been suggested that the greater the number of cigarettes smoked, the longer the healing time.5 ,68 The question of how long a patient should abstain from smoking following injury is difficult to answer and may depend on how long a patient has been a smoker. In an animal model, smoking was observed to damage tibial endothelial cell function through nitric oxide production and marked vasospasm. Long-term exposure to smoke damaged the endothelium more severely than short-term exposure, and although the effects were reversible, cessation duration was linked to benefits.72

There should be clear advice to patients on smoking cessation at all clinic visits. Patients should be offered the opportunity to enrol into smoking cessation interventions. Smokers should be encouraged to abstain following a fracture, osteotomy or arthrodesis to maximise the chances of effective healing. With elective procedures, we suggest patients abstain for as long as possible, as much as 6 months preoperatively. We also suggest enrolling patients with high-risk injuries into cessation programmes that have been shown to achieve higher cessation rates.73

Conclusions

This study shows that smokers take 27.7 days (14.2 to 41.3) longer for union to occur for fractures, osteotomy, arthrodesis and established non-union. Smokers have double the risk of non-union 2.2 (1.9 to 2.6) for fractures, osteotomy, arthrodesis and established non-union. It is important that this is explained at the same time as smoking cessation advice is given to orthopaedic patients preferably in good time, prior to elective orthopaedic surgery.

References

Footnotes

  • Contributors RGP, RGEC and BES conceived and designed the study and produced the search strategy. RGP and RGEC conducted the searches. RGP, RGEC, KLE and BES screened abstracts and titles, assessed study inclusion and exclusion, extracted data and assessed study quality. RGP, RGEC, KLE and BES were responsible for manuscript preparation and final approval of the manuscript.

  • Competing interests None declared.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement No additional data are available.