You are here

Article Text

Article menu

Download PDFPDF

XML
Epidemiology
Research
Long-term musculoskeletal morbidity after adult burn injury: a population-based cohort study
  1. Sean M Randall1,
  2. Mark W Fear2,
  3. Fiona M Wood2,3,
  4. Suzanne Rea2,3,
  5. James H Boyd1,
  6. Janine M Duke2
  1. 1Centre for Data Linkage, Curtin University, Perth, Western Australia, Australia
  2. 2Burn Injury Research Unit, School of Surgery, University of Western Australia, Perth, Western Australia, Australia
  3. 3Burns Service of Western Australia, Royal Perth Hospital and Princess Margaret Hospital, Perth, Western Australia, Australia
  1. Correspondence to Professor Janine M Duke; janine.duke{at}uwa.edu.au

Abstract

Objective To investigate if adults who are hospitalised for a burn injury have increased long-term hospital use for musculoskeletal diseases.

Design A population-based retrospective cohort study using linked administrative health data from the Western Australian Data Linkage System.

Subjects Records of 17 753 persons aged at least 20 years when hospitalised for a first burn injury in Western Australia during the period 1980–2012, and 70 758 persons who were age and gender-frequency matched with no injury admissions randomly selected from Western Australia's electoral roll.

Main outcome measures Admission rates and cumulative length of stay for musculoskeletal diseases. Negative binomial and Cox proportional hazards regression modelling were used to generate incidence rate ratios (IRR) and HRs with 95% CIs, respectively.

Results After adjustment for pre-existing health status and demographic characteristics, the burn cohort had almost twice the hospitalisation rate for a musculoskeletal condition (IRR, 95% CI 1.98, 1.86 to 2.10), and spent 3.70 times as long in hospital with a musculoskeletal diagnosis (95% CI 3.10 to 4.42) over the 33-year period, than the uninjured comparison cohort. Adjusted survival analyses of incident post-burn musculoskeletal disease admissions found significant increases for the 15-year post burn discharge period (0–6 months: HR, 95% CI 2.51, 2.04 to 3.11; 6 months–2 years: HR, 95% CI 1.77, 1.53 to 2.05; 2–15 years: HR, 95% CI 1.32, 1.23 to 1.42). Incident admission rates were significantly elevated for 20 years post-burn for minor and severe burn injury for a range of musculoskeletal diseases that included arthropathies, dorsopathies, osteopathies and soft tissue disorders.

Conclusions Minor and severe burn injuries were associated with significantly increased post-burn incident admission rates, long-term hospital use and prolonged length of stay for a range of musculoskeletal diseases. Further research is required that facilitates identification of at-risk patients and appropriate treatment pathways, to reduce the long-term morbidity associated with burns.

  • EPIDEMIOLOGY
  • SURGERY

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/

https://doi.org/10.1136/bmjopen-2015-009395

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

    • Population-based linked administrative health data provide a cost-effective means to examine long-term health impacts of burn injury.

    • Population-based linked administrative health data minimise issues of selection and reporting bias, and loss to follow-up.

    • The retrospective longitudinal study design included a comparison group.

    • Lack of individual-based risk factor data.

    • Small proportion of patient population with severe burns limited examination of long-term effects.

    Introduction

    Despite advances in surgical and medical treatment, burn injuries continue to present significant challenges to clinicians, and to burn survivors, often leading to long-term psychological and physical impairments.1 It is well documented that severe burn injury induces acute inflammatory and hypermetabolic responses that persist for at least 2 years after the initial injury.2 ,3 The subsequent metabolic demands and energy requirements are profound and induce mobilisation of proteins and amino acids, resulting in an associated increase in protein turnover, degradation and negative nitrogen balance, characteristics of serious illness.4 The physiological impacts of hypermetabolism include protein catabolism, losses of body weight, lean body mass and bone mineral density.5 ,6

    In addition to the hypermetabolic response and muscle wasting there is extensive and sustained suppression of bone formation as a result of the systemic inflammatory and endocrine responses to severe burn injury.7–9 Muscle wasting and immobility after burn injury can also directly alter the balance of bone synthesis and degradation, leading to bone loss.7 ,9–11 Burn patients are also at risk of vitamin D deficiency, contributing to bone loss.12–15 Vitamin D deficiency may develop progressively as a result of a number of factors including prolonged sun avoidance during treatments,16 hypoparathyroidism17 ,18 and low serum levels of cholesterol evident after burn injury,19 ,20 preventing the synthesis and activation of vitamin D. In adults and children, loss of bone mass and changes in bone density predispose burn patients to an increased incidence of fractures and potentially to lifelong issues associated with osteoporosis.7 ,11 ,14 ,15 ,21 ,22

    Minor burn injuries represent the majority of burn related hospital admissions in developed countries,23 ,24 and there is a growing interest in understanding the potential for systemic responses after non-severe burns. Recent research has demonstrated long-term impacts after minor burn injury on bone marrow,25 muscle, innervation and bone,26–28 with population-based findings of increased cancer incidence29 ,30 and long-term mortality.31 ,32 Non-severe burn injury has also been found to have a sustained impact on reducing trabecular bone density long after resolution of inflammation.27 However, the clinical relevance of these bone alterations is not yet clear. Investigation of the long-term effects of severe and minor burn injury is important to provide evidence for improvements in acute burn care.

    To date, limited data have been available to examine the long-term health impacts of burns. Population-based linked health administrative data provide a cost-effective means to examine long-term morbidity trends expressed in the number of hospital admissions and length of stay for specific disease classifications.33 Given the profound growth and musculoskeletal changes experienced during childhood and adolescence with the potential for different post-burn impacts, this study was limited to data of adult burn survivors 20 years and above. The aim of this study was to use population-based linked health administrative data to determine if adults hospitalised for a burn injury have increased long-term hospital use for musculoskeletal diseases, after adjustment for sociodemographic factors and pre-existing comorbidities.

    Methods

    Our study formed part of the Western Australian Population-based Burn Injury Project—a retrospective cohort investigation—that uses administrative health data from the Western Australian Data Linkage System (WADLS), a validated linkage system that links several core datasets for the entire population of Western Australia.34 The project was approved by the human research ethics committees of the University of Western Australia and the Western Australian Department of Health.

    Analyses were performed on a de-identified extraction of hospital morbidity records for all individuals who were aged at least 20 years when admitted to a hospital in Western Australia with a first burn injury between 1 January 1980 and 30 June 2012, undertaken by the WADLS. A first (index) burn injury was defined as the first hospital admission in a patient's medical record in which a burn injury was given as the principal diagnosis or an additional diagnosis, defined by International Classification of Diseases and Related Health (ICD) 9 CM 940–949 or ICD10 AM T20–T31. A population-based comparison cohort was randomly selected from Western Australia's electoral roll. Any person with an injury hospitalisation during the study period was excluded from this cohort by WADLS staff. The resultant comparison cohort was frequency matched (4:1) on birth year and sex of each burn injury case for each year from 1980 to 2012. Cohort selection and analytical methods have been reported previously.31

    Data from Western Australia's Hospital Morbidity Data System and Death Register were linked to the burn and non-injured cohorts for the period 1980–2012. Hospital admissions data included principal and additional diagnoses, age at admission, sex, and Aboriginal status, date of admission, date of discharge or other separation and mode of separation. Data supplied for the burn and non-injured cohorts also included geocoded place of residence and geocoded indices of geographical remoteness35 and social disadvantage.36 Geographical remoteness was classified into five categories: major cities, inner regional, outer regional, remote and very remote. The social disadvantage index was reclassified into quintiles (most to least disadvantaged). The mortality data included date of death and cause of death.

    An individual listed as Aboriginal or Torres Strait Islander on any admission record was categorised as Aboriginal. Supplementary codes ICD9-CM 948 or ICD10-AM T31 were used to classify the patients into those with minor burns (<20% of total body surface area (TBSA)), severe burns (≥20% TBSA) and burns of unspecified TBSA. Comorbidity (baseline) was assessed, with a 5-year look-back period, using the Charlson Comorbidity Index (CCI) and the principal and additional diagnoses included in the hospital morbidity records (0 CCI=0; 1 CCI≥1).37 ,38 A record of an existing congenital anomaly was identified using principal and additional diagnosis data (ICD9 740–759; ICD10 Q00-G99). The final discharge date for the index burn admission was used as the study start for follow-up for the burn cases and the respective frequency matched non-injury controls.

    Categorical and non-parametric continuous variables were compared using χ2 and Kruskal Wallis tests, respectively. A p value of 0.05 or lower was considered statistically significant. The total number of admissions for musculoskeletal disorders after burn injury discharge and the cumulative length of stay for principal diagnosis musculoskeletal disorders classified by subchapter headings were used as outcome measures (ICD10 AM: M00-M25—arthropathies (including infectious and inflammatory arthropathies, arthrosis and other joint disorders); M30-M36—systemic connective tissue disorders; M40-M54—dorsopathies (including deforming spinal conditions, spondylopathies and other disorders of the spine); M60-M79—soft tissue disorders (including disorders of muscles, synovium and tendons, other soft tissue disorders); M80-M94—osteopathies and chrondropathies (including disorders of bone density, osteopathies and chrondropathies). ICD-10 codes were mapped to ICD-9 codes.39 The hospitalisation of the first burn injury was not included in these outcomes. Crude yearly admission rates were calculated for these variables. Adjusted rate ratios (incidence rate ratio (IRR) and 95% CI) between the burn injury and no injury cohorts and the outcome measures were generated using negative binomial regression. Sociodemographic (gender, Aboriginality, 5-year age group, social disadvantage, remoteness of residence and year of admission) and health status information (comorbidity at baseline, history of musculoskeletal disease, congenital anomaly) were included as covariates in the models to adjust for potential confounding.

    Survival analyses of incident hospital use for arthropathies, systemic connective tissue disorders, dorsopathies, soft tissue disorders, osteopathies and chrondropathies (ie, ICD subchapter headings) as well as specific subsets of musculoskeletal diseases, infectious arthropathy, inflammatory polyarthropathy and disorders of bone density and structure, were conducted using multivariate Cox proportional hazards models. Analyses were conducted on the burn and uninjured cohorts excluding those with prior admission for musculoskeletal disease, and additionally, excluding those with a record of an injury admission in the burn cohort. Cox models were adjusted for sociodemographic and health status variables (as above). Analyses were undertaken of subgroups by burn TBSA severity and gender-specific analyses were undertaken to explore potential gender dimorphic differences in disease incidence. The proportional hazard assumption for the burn injured versus non-injured was tested using scaled and unscaled Schoenfeld residuals and by adding a group-by-time interaction term.40 Where preliminary analyses showed non-proportionality, adjusted HR and 95% CI for first record of musculoskeletal admission for burn versus non-injury cohorts were modelled for time periods guided by Aalen's linear hazard models and plots.41

    Attributable risk percentages (AR%) were calculated as the adjusted rate ratio (IRR, HR) minus one, divided by the adjusted rate (IRR, HR) ratio, multiplied by 100.42 AR% was used to estimate the proportion of long-term and incident hospital use for musculoskeletal diseases, where burn injury was a component cause.43 Statistical analyses were performed using Stata V.12 (StataCorp. LP, College Station, Texas, USA).

    Results

    The burn injury cohort included 17 753 persons aged at least 20 years when hospitalised for a first burn injury during the period January 1980–June 2012. The median age was 36 (IQR: 26–51) and 71.1% were male. Four per cent had severe burns of 20% TBSA or greater, 49% had burns of less than 20% TBSA and for 47%, the TBSA was unspecified. Eighteen per cent of the burn cohort had sustained full thickness burns, 37% partial thickness, 19% erythema and 29% had experienced burns for which the burn depth was unspecified; an individual may have had multiple burns sites and depths recorded. Among the burn cohort, 34% had a record of a non-burn injury admission (before or after index burn), and 13.8% had a previous musculoskeletal hospital admission. The burn injury cohort had a median follow-up time of 13.6 years (IQR; minimum (min)–maximum (max): 5.5–22.7; >0–32.5) for a total of 250 670 person years.

    The comparison non-injury cohort comprised 70 758 persons with median age 41 (IQR: 27–51) and males accounting for 71.0%. The uninjured cohort had a median follow-up time of 14.5 years (IQR; min–max: 6.9–23.1; 0.01 to 32.5) for a total of 1 067 568 person years. Refer to table 1 for other baseline sociodemographic and health status variables for the burn and non-injury cohorts. The burn injury cohort comprised significantly higher proportions of Aboriginal people, people who were socially disadvantaged, people living in regions outside of major cities and people who had pre-existing comorbidity when compared with the non-injured comparison cohort.

    View this table:
    Table 1

    Baseline demographic and pre-existing health status factors for the burn injury (total) and the frequency matched non-injury (total) cohorts

    Admissions for diseases of musculoskeletal system—rates and cumulative length of stay

    There were 10 761 hospital admissions occurring after burn hospitalisation discharge with a primary diagnosis of a musculoskeletal disease. Arthropathies were the most common cause of musculoskeletal admission, followed by dorsopathies (see table 2). A total of 55 810 days were spent in hospital for musculoskeletal diagnosis after a burn hospitalisation. The median length of hospital stay was 1 day (IQR: 0–5 days). The total number of days spent in hospital with a musculoskeletal primary diagnosis by the uninjured cohort was 68 946, where the median length of stay of musculoskeletal admissions was 1 day (IQR: 0–3 days). The length of stay of 0 days represents admission and discharge occurring on the same day.

    View this table:
    Table 2

    Number of admissions (%) for musculoskeletal diseases classified by ICD subchapter codes in the non-injury and burn injury cohorts

    Unadjusted incidence rates for musculoskeletal admissions and lengths of stay are shown in figure 1. These graphs show a higher rate of admissions for musculoskeletal diseases, and more time spent in hospital with musculoskeletal diseases for burn patients compared with uninjured patients over the entire 33-year period.

    Figure 1
    Figure 1

    Unadjusted rates of hospital admissions and cumulative length of stay (per 100 person-years (PYs)) for musculoskeletal diseases (total) among adults with burn injury versus no injury.

    After adjustment for pre-existing health and sociodemographic characteristics, the burn cohort had almost twice as many hospitalisations for a musculoskeletal condition (IRR, 95% CI 1.98, 1.86 to 2.10), and spent 3.70 times as long in hospital with a musculoskeletal diagnosis (95% CI 3.10 to 4.42) over the 33-year period, than the uninjured comparison cohort. The adjusted AR% suggested that the burn cohort experienced an excess of 49.5% of admissions for musculoskeletal conditions (n=5329) and 73.0% of all days spent in hospital for musculoskeletal conditions (n=40 727 days), when compared to the uninjured cohort.

    Unadjusted admission rates for the musculoskeletal system sub-categories are shown in figure 2. After controlling for demographic factors and previous health status, those who had a burn had higher rates of admissions for arthropathies (IRR, 95% CI 1.64, 1.54 to 1.75), dorsopathies (IRR, 95% CI 2.16, 1.93 to 2.41), osteopathies and chondropathies (IRR, 95% CI 5.64, 4.56 to 6.97) and soft tissue disorders (IRR, 95% CI 1.91, 1.76 to 2.06). There was no difference in the rate of admissions for connective tissue disorders between the burn and the non-injured cohorts (IRR, 95% CI 0.88, 0.48 to 1.62). The adjusted AR% suggested that the burn cohort experienced excess post-burn admissions of 39.0% for arthropathies (n=1702), 53.7% for dorsopathies (n=1925), 82.3% for osteopathies and chrondropathies (n=629), and 47.6% for soft tissue disorders (n=832), when compared with the uninjured cohort.

    Figure 2
    Figure 2

    Unadjusted rates (per 100 person-years (PYs)) of hospital admissions by musculoskeletal disease subgroup for adults with burn injury versus no injury.

    Examination of more specific subconditions showed a large and significant increase in the rate of admissions for disorders of bone density (IRR, 95% CI 13.87, 9.89 to 19.44), along with increases in infectious arthropathies (IRR, 95% CI 1.98, 1.46 to 2.69) and inflammatory polyarthropathies (IRR, 95% CI 3.82, 2.61 to 5.59) for those hospitalised with a burn injury.

    A similar picture was shown for the length of time spent in hospital, with longer cumulative length of stay after a burn hospitalisation with arthropathies (IRR, 95% CI 1.64, 1.54 to 1.75), dorsopathies (IRR, 95% CI 2.16, 1.93 to 2.41), osteopathies and chondropathies (IRR, 95% CI: 5.64, 4.56 to 6.97), and soft tissue disorders (IRR, 95% CI 1.91, 1.76 to 2.06). Those with a burn injury spent almost 50 times longer in hospital with disorders of bone density compared to the uninjured cohort (IRR, 95% CI 49.24, 28.72 to 84.42), with increases in cumulative time spent in hospital also found for infectious arthropathies (IRR, 95% CI 13.1, 5.62 to 30.56) and inflammatory polyarthropathies (IRR, 95% CI 1.98, 1.46 to 2.69). No significant difference was found between the burn and uninjured cohorts for length of time in hospital with connective tissue disorders (IRR, 95% CI 0.65, 0.23 to 1.85).

    Analysis by TBSA severity (see table 3) showed that increased hospitalisation rates were associated with severe and minor burns, with those with minor burns having higher admission rates for arthropathies, dorsopathies, osteopathies and chondropathies, and soft tissue disorders when compared with the uninjured cohort.

    View this table:
    Table 3

    Adjusted IRR and 95% CIs for recurrent admissions musculoskeletal subconditions, by burn severity compared with the uninjured cohort

    Incidence—survival analysis

    Analysis of time until first musculoskeletal admission (post-burn) was performed on the burn and uninjured cohorts who did not have a prior record of musculoskeletal hospitalisation, with the additional exclusion from the burn cohort of those with a record of a principal diagnosis injury admission. There were 10 440 individuals in this reduced burn cohort, of which 1779 had a first time (incident) musculoskeletal admission within the study period. The uninjured comparison cohort comprised of 66 869 controls, of which 9689 had an incident musculoskeletal admission in this study period.

    Results of multivariate Cox regression modelling for an incident musculoskeletal (combined diseases) admission revealed evidence of non-proportionality. First time musculoskeletal (combined diseases) admissions were most frequent immediately after burn and while remaining significantly elevated, decreased over the study period. Adjusted analyses conducted on partitioned time windows found significant HRs for the first 6 months after burn (HR, 95% CI 2.51, 2.04 to 3.11), from 6 months to 2 years after burn (HR, 95% CI 1.77, 1.53 to 2.05), and from 2 to 15 years after burn (HR, 95% CI 1.32, 1.23 to 1.42). Differences were found for severe (HR, 95% CI 4.37, 2.32 to 8.25), minor (HR, 95% CI 2.25, 1.68 to 3.01) and unknown TBSA (HR, 95% CI 2.42, 1.74 to 3.37) for the first 6 months; only minor burns (HR, 95% CI 6 months–2 years; 1.68, 1.39 to 2.04, 2–15 years; 1.44, 1.30 to 1.59) and unknown TBSA burns (HR, 95% CI 6 months–2 years; 1.83, 1.54 to 2.17, 2–15 years; 1.28, 1.17 to 1.41) showed increased admissions over the first 15 years.

    After adjustment for confounders, both males (HR, 95% CI 0–6 months; 2.62, 2.03 to 3.39, 6 months–2 years; 1.83, 1.54 to 2.17, 2–15 years; 1.31, 1.21 to 1.42) and females (HR, 95% CI 0–6 months; 2.09, 1.35 to 3.23, 6 months–2 years; 1.68, 1.27 to 2.23, 2–15 years; 1.34, 1.17 to 1.53) with burns showed significant increases over the first 15 years post-burn, compared with uninjured males and females, respectively.

    In the burn cohort, there were 968 first time admissions for arthropathies, 468 for dorsopathies, 155 for osteopathies and chrondropathies (n=70 bone density disorders), 541 for soft tissue disorders and 8 first time admissions for connective tissue disorders. First time admissions for infectious arthropathy and inflammatory polyarthropathy accounted for 24 and 67 of the total incident arthrropathy admissions, with median (IQR) times to incident admission of 8.4 years (3.4–16.3) and 9.4 years (3.8–16.9), respectively. Results for adjusted Cox regression models for incident musculoskeletal subconditions are shown in table 4. Evidence of non-proportionality was common for numerous subconditions, with higher rates of admissions in the burn cohort more immediately after the burn injury. All subconditions showed significantly higher rates of first admissions in the first 20 years after burn. The AR% calculated suggested that 567 first admissions to hospital for musculoskeletal diseases could be attributed to burn injury. The largest proportion of these incident admissions attributable to the burn injury were for arthropathies (35.4%), followed by diseases of soft tissue (30.0%), dorsopathies (20.0%) and finally, osteopathies and chrondropathies (14.5%).

    View this table:
    Table 4

    Adjusted HRs for first time post-burn admissions for musculoskeletal subconditions, comparing the burn cohort with the uninjured cohort

    Discussion

    This study quantifies the increased population-based long-term hospital use for musculoskeletal disorders experienced by adults hospitalised with burn injury, after controlling for demographic and pre-existing comorbidities. Members of the burn injury cohort had 1.98 times the rate of hospitalisations and 3.70 times the length of stay in hospital for combined musculoskeletal disorders post-burn when compared with the uninjured cohort. Arthropathies and dorospathies combined accounted for 74% of all readmissions (prevalent and incident disease admissions) for musculoskeletal disorders. After controlling for sociodemographic and pre-existing health status, the burn cohort experienced significantly higher rates of hospitalisations post-burn for arthropathies, dorsopathies, osteopathies and chrondropathies and soft tissue disorders. No difference was found for admission rates for diseases of the connective tissue when compared with the uninjured cohort.

    Assessment of the impact of burn severity on recurrent admissions for musculoskeletal diseases identified significantly elevated rates of similar magnitude for severe and minor burn injury for arthropathies, dorsopathies and soft tissue disorders. However, while admission rates for osteopathies and chrondropathies were significantly elevated for severe burns (3.32 times higher), minor burns were associated with an admisson rate of 8.44 times that observed for the uninjured cohort over the study period. This lack of a dose–response relationship may in part be due to the small numbers of burn survivors with severe burns in this patient population, the TBSA classification used and/ or a ‘healthy’ survivor effect. An additional possibility is that patients with severe burns have supportive treatments to mitigate systemic inflammatory and endocrine responses and are more likely to have higher levels of continued post-burn care that may lead to earlier diagnoses and management of secondary pathologies, resulting in reduced hospitalisations for more serious presentations. The impact of minor burn injury on musculoskeletal disease admissions is interesting in light of recent evidence of less severe burn causing depletion of trabecular bone volume.27

    To exclude any potential additive systemic impacts associated with other non-burn injury among members of the burn cohort,43 ,44 first time admissions for musculoskeletal diseases post-burn were examined excluding those in the burn cohort with a principal diagnosis injury admission. After adjustment for sociodemographic and health factors, significantly elevated rates of first time post-burn hospitalisations for musculoskeletal disorders were identified. The adjusted AR% for incident post-burn musculoskeletal disease admissions suggested that 567 first time admissions experienced in the burn cohort after discharge, could be attributed to burn injury. Gender-specific analyses found males and females with burn injury to have significantly elevated admission rates.

    Long-term admissions for infectious and inflammatory arthropathies were assessed as a means to investigate the potential for chronic impacts of the initial systemic inflammatory and immune response to burn injury.45–47 Significantly increased incident admission rates for infectious and inflammatory arthropathies post-burn were identified with effects persisting for at least 20 years, suggesting that burn patients are at increased long-term susceptibility to infectious and inflammatory diseases. The increases in inflammatory and infectious diseases may be related to the immune suppression induced by burn injury;46 ,48 ,49 however, recent evidence also suggests long-term changes in the immune response post-burn, including reduced dendritic cell activation and inability of dendritic cells to activate T cells.44 While the findings of considerably elevated incident admission rates for bone density disorders during the first 2 years post-burn were consistent with other research,7 ,11 ,14 ,15 ,21 ,22 the burn cohort continued to experience significantly elevated incident admissions over the 33-year study period. Burn injury affects multiple systems of the body including the immune, metabolic and endocrine systems, and the pathophysiological mechanisms that underlie the different musculoskeletal conditions discussed in this paper are most likely variable. Future research will be important to identify how burn injury triggers different mechanisms in order to inform early intervention and prevent these long-term pathologies.

    Strengths and limitations

    Through the use of linked administrative health data, we performed a large scale population-based study with a long follow-up time, accurate pre- and post burn injury measures,50 a non-injured comparison group and sample size sufficient for quantitative analyses. The analytical strategy of the study was based on the assumption that after controlling for confounders, any excess in hospital use for musculoskeletal diseases in the burn cohort (compared to the non-injured), was an outcome principally associated with the incident burn. We were able to include variables of health, social disadvantage and geographic access to services, to determine the individual burden of hospital readmissions for musculoskeletal disorders and quantify at the population level. While health administrative datasets do not include variables of risk-taking behaviour, analyses undertaken of the burn cohort excluding members with a record of a non-burn injury admission and potentially those of high risk-taking behaviour, also found significantly elevated incident musculoskeletal admission rates. Likewise, models were adjusted for social disadvantage, a factor previously correlated with lifestyle risk factors (eg, nutrition, smoking, alcohol)51 ,52 that could be associated with the burn exposure and the outcomes measured. Changes in ICD coding and incomplete TBSA% data may have limited a complete understanding of burn severity on long-term hospital use; however, significantly increased admission rates were observed for severe and minor burns and burns of unspecified TBSA. Any differential effect of burn injury on the incidence of fractures could not be evaluated in this study due to the exclusion of any injury admissions (including fractures) in the comparator uninjured cohort. Further research is planned that will include a non-burn injury cohort where differences in long-term musculoskeletal morbidity will be examined, including fractures.

    The out-of-state migration for Western Australia is consistently low at levels below 3%; such losses were not anticipated to bias the results.53 This study represents burn injury and musculoskeletal diseases serious enough to require hospitalisation and the results may underestimate the impact of burn injury experienced in the community. Research of functional outcomes and quality-of-life of burn patients in Western Australia have provided valuable data on patients during the first year after discharge;54 ,55 however, this study has generated new information on the longer-term musculoskeletal morbidity experienced long after healing of the burn wound and cessation of attendance at outpatient burn clinics. Further work that links individuals’ pharmaceutical and or primary care data with hospital data will provide a clearer picture of the time of diagnosis of secondary musculoskeletal pathologies and treatment pathways post-burn. We expect that these results would be generalisable to other countries of similar demographics and healthcare systems.

    After adjustment of confounders, increased rates of first time admissions and readmissions for musculoskeletal diseases were identified for those with severe and minor burns when compared with a non-injured cohort. These findings highlight long-term effects of burn injury, including minor burn injury, on musculoskeletal morbidity. Further research is required that facilitates identification of at-risk patients, the mechanisms that may be responsible for these long-term effects and appropriate treatment pathways, to reduce the long-term morbidity associated with burns.

    Acknowledgments

    The authors thank the staff of the Health Information Linkage Branch for access to the Western Australian Data Linkage System and for their assistance in obtaining the data, the WA Health Data Custodians for access to the core health datasets and the Western Australian Department of Health. The project was supported by the Raine Medical Research Foundation (Priming grant JMD). JMD is supported by Woodside Corporation. MWF is supported by Chevron Australia.

    References

    1. ↵
      1. Jeschke MG,
      2. Chinkes DL,
      3. Finnerty CC, et al
      . Pathophysiologic response to severe burn injury. Ann Surg 2008;248:387–401. doi:10.1097/SLA.0b013e318176c4b3
      OpenUrlPubMedWeb of Science
    2. ↵
      1. Przkora R,
      2. Barrow RE,
      3. Jeschke MG, et al
      . Body composition changes with time in pediatric burn patients. J Trauma 2006;60:968–71. doi:10.1097/01.ta.0000214580.27501.19
      OpenUrlCrossRefPubMed
    3. ↵
      1. Przkora R,
      2. Jeschke MG,
      3. Barrow RE, et al
      . Metabolic and hormonal changes of severely burned children receiving long-term oxandrolone treatment. Ann Surg 2005;242:384–9, discussion 390-1.
      OpenUrlPubMedWeb of Science
    4. ↵
      1. Chang DW,
      2. DeSanti L,
      3. Demling RH
      . Anticatabolic and anabolic strategies in critical illness: a review of current treatment modalities. Shock 1998;103:155–60. doi:10.1097/00024382-199809000-00001
      OpenUrl
    5. ↵
      1. Rennie MJ
      . Muscle wasting in the muscular dystrophies. Dev Med Child Neurol 1985;27:524–7. doi:10.1111/j.1469-8749.1985.tb04579.x
      OpenUrlPubMed
    6. ↵
      1. Rennie MJ
      . Muscle protein turnover and the wasting due to injury and disease. Br Med Bull 1985;41:257–64.
      OpenUrlAbstract/FREE Full Text
    7. ↵
      1. Klein GL
      . Burn-induced bone loss: importance, mechanisms, and management. J Burns Wounds 2006;5:e5.
      OpenUrlPubMed
    8. ↵
      1. Klein GL,
      2. Wolf SE,
      3. Goodman WG, et al
      . The management of acute bone loss in severe catabolism due to burn injury. Horm Res 1997;48(Suppl 5):83–7. doi:10.1159/000191334
      OpenUrlPubMedWeb of Science
    9. ↵
      1. Leblebici B,
      2. Sezgin N,
      3. Ulusan SN, et al
      . Bone loss during the acute stage following burn injury. J Burn Care Res 2008;29:763–7. doi:10.1097/BCR.0b013e31818480f4
      OpenUrlCrossRefPubMed
    10. ↵
      1. Klein GL,
      2. Bi LX,
      3. Sherrard DJ, et al
      . Evidence supporting a role of glucocorticoids in short-term bone loss in burned children. Osteoporos Int 2004;15:468–74.
      OpenUrlPubMed
    11. ↵
      1. Klein GL,
      2. Herndon DN,
      3. Goodman WG, et al
      . Histomorphometric and biochemical characterization of bone following acute severe burns in children. Bone 1995;17:455–60. doi:10.1016/8756-3282(95)00279-1
      OpenUrlCrossRefPubMed
    12. ↵
      1. Garrel D
      . Burn scars: a new cause of vitamin D deficiency? Lancet 2004;363:259–60. doi:10.1016/S0140-6736(03)15417-7
      OpenUrlPubMed
    13. ↵
      1. Klein GL
      . The interaction between burn injury and vitamin D metabolism and consequences for the patient. Curr Clin Pharmacol 2008;3:204–10. doi:10.2174/157488408785747647
      OpenUrlCrossRefPubMed
    14. ↵
      1. Klein GL
      . Burns: where has all the calcium (and vitamin D) gone? Adv Nutrition 2011;2:457–62. doi:10.3945/an.111.000745
      OpenUrl
    15. ↵
      1. Klein GL
      . Does vitamin D deficiency contribute to post-burn bone loss? F1000Res 2012;1:57. doi:10.12688/f1000research.1-57.v1
      OpenUrl
    16. ↵
      1. Norman AW
      . Sunlight, season, skin pigmentation, vitamin D, and 25-hydroxyvitamin D: integral components of the vitamin D endocrine system. Am J Clin Nutr 1998;67:1108–10.
      OpenUrlPubMedWeb of Science
    17. ↵
      1. Klein GL,
      2. Langman CB,
      3. Herndon DN
      . Persistent hypoparathyroidism following magnesium repletion in burn-injured children. Pediatr Nephrol 2000;14:301–4. doi:10.1007/s004670050763
      OpenUrlCrossRefPubMed
    18. ↵
      1. Murphey ED,
      2. Chattopadhyay N,
      3. Bai M, et al
      . Up-regulation of the parathyroid calcium-sensing receptor after burn injury in sheep: a potential contributory factor to postburn hypocalcemia. Crit Care Med 2000;28:3885–90. doi:10.1097/00003246-200012000-00024
      OpenUrlCrossRefPubMed
    19. ↵
      1. Gottschlich MM,
      2. Alexander JW
      . Fat kinetics and recommended dietary intake in burns. JPEN J Parenter Enteral Nutr 1987;11:80–5. doi:10.1177/014860718701100180
      OpenUrlAbstract/FREE Full Text
    20. ↵
      1. Gottschlich MM,
      2. Jenkins M,
      3. Warden GD, et al
      . Differential effects of three enteral dietary regimens on selected outcome variables in burn patients. JPEN J Parenter Enteral Nutr 1990;14:225–36. doi:10.1177/0148607190014003225
      OpenUrlAbstract/FREE Full Text
    21. ↵
      1. Mayes T,
      2. Gottschlich M,
      3. Scanlon J, et al
      . Four-year review of burns as an etiologic factor in the development of long bone fractures in pediatric patients. J Burn Care Rehabil 2003;24:279–84. doi:10.1097/01.BCR.0000085844.84144.E0
      OpenUrlCrossRefPubMed
    22. ↵
      1. Edelman LS,
      2. McNaught T,
      3. Chan GM, et al
      . Sustained bone mineral density changes after burn injury. J Surg Res 2003;114:172–8. doi:10.1016/S0022-4804(03)00275-0
      OpenUrlCrossRefPubMed
    23. ↵
      1. Akerlund E,
      2. Huss FR,
      3. Sjoberg F
      . Burns in Sweden: an analysis of 24,538 cases during the period 1987–2004. Burns 2007;33:31–6. doi:10.1016/j.burns.2006.10.002
      OpenUrlCrossRefPubMedWeb of Science
    24. ↵
      1. Duke J,
      2. Wood F,
      3. Semmens J, et al
      . A 26-year population-based study of burn injury hospital admissions in Western Australia. J Burn Care Res 2011;32:379–86. doi:10.1097/BCR.0b013e318219d16c
      OpenUrlCrossRefPubMedWeb of Science
    25. ↵
      1. Rea S,
      2. Giles NL,
      3. Webb S, et al
      . Bone marrow-derived cells in the healing burn wound—more than just inflammation. Burns 2009;35:356–64. doi:10.1016/j.burns.2008.07.011
      OpenUrlCrossRefPubMedWeb of Science
    26. ↵
      1. Anderson JR,
      2. Zorbas JS,
      3. Phillips JK, et al
      . Systemic decreases in cutaneous innervation after burn injury. J Invest Dermatol 2010;130:1948–51. doi:10.1038/jid.2010.47
      OpenUrlCrossRefPubMedWeb of Science
    27. ↵
      1. O'Halloran E,
      2. Kular J,
      3. Xu J, et al
      . Non-severe burn injury leads to depletion of bone volume that can be ameliorated by inhibiting TNF-alpha. Burns 2015;41:558–64. doi:10.1016/j.burns.2014.09.004
      OpenUrl
    28. ↵
      1. Morellini NM,
      2. Fear MW,
      3. Rea S, et al
      . Burn injury has a systemic effect on reinnervation of skin and restoration of nociceptive function. Wound Repair Regen 2012;20:367–77. doi:10.1111/j.1524-475X.2012.00787.x
      OpenUrlPubMed
    29. ↵
      1. Duke J,
      2. Rea S,
      3. Semmens J, et al
      . Burn Injury and cancer risk: a state-wide longitudinal study. Burns 2011;38:340–7. doi:10.1016/j.burns.2011.10.003
      OpenUrl
    30. ↵
      1. Duke JM,
      2. Bauer J,
      3. Fear MW, et al
      . Burn injury, gender and cancer risk: population-based cohort study using data from Scotland and Western Australia. BMJ Open 2014;4:e003845. doi:10.1136/bmjopen-2013-003845
      OpenUrlAbstract/FREE Full Text
    31. ↵
      1. Duke JM,
      2. Boyd JH,
      3. Rea S, et al
      . Long-term mortality among older adults with burn injury: a population-based study in Australia. Bull World Health Organ 2015;93:400–6. doi:10.2471/BLT.14.149146
      OpenUrl
    32. ↵
      1. Duke JM,
      2. Rea S,
      3. Boyd JH, et al
      . Mortality after burn injury in children: a 33-year population-based study. Pediatrics 2015;135:e903–10. doi:10.1542/peds.2014-3140
      OpenUrlAbstract/FREE Full Text
    33. ↵
      1. Hendrie D,
      2. Miller TR
      . Assessing the burden of injuries: competing measures. Inj Control Saf Promot 2004;11:193–9. doi:10.1080/156609704/233/289689
      OpenUrlCrossRefPubMed
    34. ↵
      1. Holman CD,
      2. Bass AJ,
      3. Rouse IL, et al
      . Population-based linkage of health records in Western Australia: development of a health service research linked database. Aust NZ J Public Health 1999;23:453–9. doi:10.1111/j.1467-842X.1999.tb01297.x
      OpenUrlCrossRefPubMedWeb of Science
    35. ↵
      1. Glover J,
      2. Tennant S
      . Remote areas statistical geography in Australia: notes on the Accessibility/Remoteness Index for Australia (ARIA+ version). Working Papers Series No. 9. Adelaide: Public Health Information Development Unit, Adelaide, The University of Adelaide, 2003.
    36. ↵
      1. Trewin D
      . Socio-economic indexes for areas (information paper, census of population and housing). Canberra: Australian Bureau of Statistics, 2003.
    37. ↵
      1. Charlson ME,
      2. Pompei P,
      3. Ales KL, et al
      . A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chronic Dis 1987;40:373–83. doi:10.1016/0021-9681(87)90171-8
      OpenUrlCrossRefPubMedWeb of Science
    38. ↵
      1. Preen DB,
      2. Holman CDAJ,
      3. Spilsbury K, et al
      . Length of comorbidity lookback period affected regression model performance of administrative health data. J Clin Epidemiol 2006;59:940–6. doi:10.1016/j.jclinepi.2005.12.013
      OpenUrlCrossRefPubMedWeb of Science
    39. ↵
      Ministry of Health—Manatu Haurora NZ. Mapping between ICD-10 and ICD-9. Secondary Mapping between ICD-10 and ICD-9. 30 May 2015. http://www.health.govt.nz/nz-health-statistics/data-references/mapping-tools/mapping-between-icd-10-and-icd-9
    40. ↵
      1. Hosmer DW,
      2. Lemeshow S
      . Applied survival analysis: regression modeling of time to event data. New York: Wiley, 1999.
    41. ↵
      1. Hosmer DW,
      2. Royston P
      . Using Aalen's linear hazards model to investigate time-varying effects in the proportional hazards regression model. Stata J 2002;2:331–50.
      OpenUrl
    42. ↵
      1. Gordis L
      . Epidemiology. 2nd edn. Philadelphia: W.B. Saunders Company, 2000.
    43. ↵
      1. Cameron CM,
      2. Purdie DM,
      3. Kliewer EV, et al
      . Ten-year health service use outcomes in a population-based cohort of 21,000 injured adults: the Manitoba injury outcome study. Bull World Health Organ 2006;84:802–10. doi:10.2471/BLT.06.030833
      OpenUrlCrossRefPubMedWeb of Science
    44. ↵
      1. Valvis SM,
      2. Waithman J,
      3. Wood FM, et al
      . The Immune Response to Skin Trauma is Dependent on the Etilology of Injury in a Mouse Model of Burn and Excision. J Invest Dermatol 2015;135:2119–28. doi:10.1038/jid.2015.123
      OpenUrl
    45. ↵
      1. Mace JE,
      2. Park MS,
      3. Mora AG, et al
      . Differential expression of the immunoinflammatory response in trauma patients: burn vs. non-burn. Burns 2012;38:599–606. doi:10.1016/j.burns.2011.10.013
      OpenUrlCrossRefPubMedWeb of Science
    46. ↵
      1. O'Sullivan ST,
      2. O'Connor TP
      . Immunosuppression following thermal injury: the pathogenesis of immunodysfunction. Br J Plast Surg 1997;50:615–23. doi:10.1016/S0007-1226(97)90507-5
      OpenUrlCrossRefPubMedWeb of Science
    47. ↵
      1. Rock KL,
      2. Lai JJ,
      3. Kono H
      . Innate and adaptive immune responses to cell death. Immunol Rev 2011;243:191–205. doi:10.1111/j.1600-065X.2011.01040.x
      OpenUrlCrossRefPubMed
    48. ↵
      1. D'Arpa N,
      2. Accardo-Palumbo A,
      3. Amato G, et al
      . Circulating dendritic cells following burn. Burns 2009;35:513–18. doi:10.1016/j.burns.2008.05.027
      OpenUrlCrossRefPubMed
    49. ↵
      1. Patenaude J,
      2. D'Elia M,
      3. Hamelin C, et al
      . Selective effect of burn injury on splenic CD11c(+) dendritic cells and CD8alpha(+)CD4(-)CD11c(+) dendritic cell subsets. Cell Mol Life Sci 2010;67:1315–29. doi:10.1007/s00018-009-0251-8
      OpenUrlPubMed
    50. ↵
      Department of Health Western Australia. Clinical Information Audit Program—Hospital Activity Report. Operational Directive OD 0201/09. Perth Department of Health WA, 2009.
    51. ↵
      1. Mishra G,
      2. Ball K,
      3. Patterson A, et al
      . Socio-demographic inequalities in the diets of mid-aged Australian women. Eur J Clin Nutr 2005;59:185–95. doi:10.1038/sj.ejcn.1602057
      OpenUrlCrossRefPubMedWeb of Science
    52. ↵
      1. Siahpush M,
      2. Borland R,
      3. Scollo M
      . Prevalence and socio-economic correlates of smoking among lone mothers in Australia. Aust NZ J Public Health 2002;26:132–5. doi:10.1111/j.1467-842X.2002.tb00905.x
      OpenUrlPubMed
    53. ↵
      1. Clark A,
      2. Preen DB,
      3. Ng JQ, et al
      . Is Western Australia representative of other Australian States and Territories in terms of key socio-demographic and health economic indicators? Aust Health Rev 2010;34:210–15. doi:10.1071/AH09805
      OpenUrlCrossRefPubMedWeb of Science
    54. ↵
      1. Grisbrook TL,
      2. Elliott CM,
      3. Edgar DW, et al
      . Burn-injured adults with long term functional impairments demonstrate the same response to resistance training as uninjured controls. Burns 2013;39:680–6. doi:10.1016/j.burns.2012.09.005
      OpenUrl
    55. ↵
      1. Grisbrook TL,
      2. Reid SL,
      3. Edgar DW, et al
      . Exercise training to improve health related quality of life in long term survivors of major burn injury: a matched controlled study. Burns 2012;38:1165–73. doi:10.1016/j.burns.2012.03.007
      OpenUrlPubMed

    Footnotes

    • Contributors The manuscript, including related data, figures and tables has not been previously published and that the manuscript is not under consideration elsewhere. All authors have made contributions to the paper and authorised the submission: JMD designed the study, supported data analyses and interpretation and drafted manuscript. SMR provided data management, analyses and data interpretation. JHB provided statistical advice. MWF, SR and FMW contributed clinical interpretation. All authors contributed to manuscript preparation and critical revision.

    • Funding Project data costs were supported by a Raine Medical Research Foundation Priming grant (JMD). JMD is supported by Woodside corporate sponsorship and the Fiona Wood Foundation. MWF is supported by Chevron Australia and the Fiona Wood Foundation.

    • Competing interests None declared.

    • Ethics approval This study has received ethics approval from the University of Western Australia and the Western Australian Department of Health. Being a large population cohort study using de-identified linked data, approval included a waiver of informed consent.

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

    • Data sharing statement No additional data are available.