OBJECTIVES

Rhabdomyolysis in children is a highly variable condition with presentations ranging from myalgias to more severe complications like acute renal failure. We sought to explore demographics and incidence of pediatric rhabdomyolysis hospitalizations and rates of associated renal failure, as our current understanding is limited.

METHODS

This was a retrospective analysis using the Healthcare Cost and Utilization Project Kids’ Inpatient Database to identify children hospitalized with a primary diagnosis of rhabdomyolysis. Data were analyzed for demographic characteristics, as well as geographic and temporal trends. Multivariable logistic regression was used to identify characteristics associated with rhabdomyolysis-associated acute renal failure.

RESULTS

From 2006 to 2016, there were 8599 hospitalized children with a primary diagnosis of rhabdomyolysis. Overall, hospitalizations for pediatric rhabdomyolysis are increasing over time, with geographic peaks in the South and Northeast regions, and seasonal peaks in March and August. Though renal morbidity was diagnosed in 8.5% of children requiring hospitalization for rhabdomyolysis, very few of these patients required renal replacement therapy (0.41%), and death was rare (0.03%). Characteristics associated with renal failure included male sex, age greater than 15 years, and non-Hispanic Black race.

CONCLUSIONS

Though renal failure occurs at a significant rate in children hospitalized with rhabdomyolysis, severe complications, including death, are rare. The number of children hospitalized with rhabdomyolysis varies by geographic region and month of the year. Future studies are needed to explore etiologies of rhabdomyolysis and laboratory values that predict higher risk of morbidity and mortality in children with rhabdomyolysis.

Rhabdomyolysis in children is a highly variable condition most commonly associated with infection, particularly viral myositis, followed by trauma, exercise, and surgery.13  Some children have mild presentations with elevation of serum muscle enzymes, myalgias or muscle weakness, and/or myoglobinuria (dark-colored urine). Others may develop more severe complications, including life-threatening electrolyte or acid-base disturbances, disseminated intravascular coagulation, or compartment syndrome. The most feared complication, acute kidney injury, can vary in severity from transient mild elevation in serum creatinine to florid kidney failure requiring renal replacement therapy (RRT).4 

However, there are limited pediatric studies describing the incidence of morbidity and mortality associated with rhabdomyolysis; among these, reports of rhabdomyolysis-induced acute kidney injury or failure ranged from 0% to 50% and mortality ranged from 0% to 33%.13,513  The largest of these case series involved 319 patients.8  These studies are limited by small size, data from a single context (emergency department, inpatient, or ICU), or a limited number of study sites, precluding generalizable conclusions. General expert consensus recommends hospitalization for intravenous fluid administration to prevent renal morbidity and mortality, but there is limited data to identify those who would benefit most from aggressive therapy for prevention of renal morbidity and mortality.14 

The primary aim of this study was to establish the incidence of hospitalization for rhabdomyolysis among children in the United States. Within this population, we aimed to determine rates of renal morbidity and mortality as well as demographic factors associated with increased renal morbidity or mortality. The secondary aim was to explore geographic and temporal variation.

This was a retrospective analysis of the Healthcare Cost and Utilization Project Kids’ Inpatient Database (KID), a publicly available all-payer database that estimates pediatric hospitalizations (age <21 years) across the United States, estimating roughly 7 million inpatient stays (when weighted). In 2006, the KID database included over 3700 nonrehabilitation hospitals in 38 states; this increased to 4200 hospitals in 47 states in 2016. All hospitalized patients in 2006, 2009, 2012, and 2016 with a primary diagnosis code for rhabdomyolysis (International Classification of Diseases, Ninth and Tenth Revision, Clinical Modification [ICD], ICD-9 code 728.88 [2006–2012]; ICD-10 code M62.82 [2016]) were included. To exclude patients with other diagnoses that may independently cause renal injury, like septic shock or medication-induced, we chose not to include patients who only had a secondary diagnosis of rhabdomyolysis. US Census Bureau yearly population estimates by age group and geographic region were used to determine incidence and rates over time.15  This study, utilizing deidentified data, was not considered human subjects research by the official institutional review board.

Demographic characteristics included US census-designated age groups (infant to young child [less than 5 years], older child [5 to 9 years], early adolescence [10 to14 years], late adolescence [15 to 20 years]), gender (male or female), race (non-Hispanic White, non-Hispanic Black, Hispanic, Asian or Pacific Islander, Native American, or other), and payor (Medicare, Medicaid, private insurance, self-pay, or other). Race was included, as a social construct, to evaluate for diverse sociodemographic representation of groups that may impact generalizability of data or suggest greater burden of disease within various groups. The number of complex chronic conditions (CCCs) were counted (0, 1, 2+CCCs) using a prior described methodology, and patient disease severity was quantified using the All Patient Refined Diagnosis Related Groups classification system with relative weights from Hospitalization Resource Intensity Scores for Kids.16,17  Hospital length of stay (days) represented resource utilization. Finally, we grouped hospitals by United States census-designated geographic region (Northeast, Midwest, South, West), and admission month.18 

The primary outcome was incidence of hospitalizations for pediatric rhabdomyolysis and associated renal morbidity or death. Overall, incidence was calculated using the number of pediatric hospitalizations for rhabdomyolysis compared with the general population of children reported by the United States census from 2006 to 2016. To generate national estimates, we used the discharge-level weight available in the database to estimate patients treated at community hospitals (excluding rehabilitation and long-term acute care facilities) in the United States. These rates were analyzed for trends over time, by geographic region, and by month of the year. To estimate the disease burden among hospitalized children, the number of hospitalizations for rhabdomyolysis was also compared with the overall number of children hospitalized over the same period. Renal morbidity was defined as acute renal insufficiency and/or acute renal failure (including need for RRT). ICD-9 and ICD-10 diagnosis codes were used to identify patients with concurrent diagnoses of rhabdomyolysis and acute renal insufficiency, injury and/or failure; ICD-9 and ICD-10 procedure codes for dialysis or hemodialysis were used to identify patients who received RRT (Table 3). Patients with the ICD-10 diagnosis code for “disorder of kidney and ureter, unspecified” (N28.9) were not included, as this was less specific for acute renal failure. As a secondary outcome, patient characteristics were analyzed for independent predictors of renal morbidity diagnoses.

All analyses made use of the sampling frames and weights available in the Healthcare Cost and Utilization Project KID database. Bivariate analyses describing renal failure by baseline demographics was conducted using χ2 and Wilcoxon Rank Sum tests as appropriate. Because of the non-normal distribution of continuous data, geometric means and standard deviation (SD) were used. Rates of rhabdomyolysis were analyzed over time with a generalized linear model with a Poisson distribution and log link using a fixed effect for year. The log transformation of the general population of children for each year (derived from census data) was used as an offset in the model. Comparisons across geographic regions were performed by adding an interaction term to this model between census region and year. Multivariable logistic regression model with backward selection was used to estimate adjusted odds ratio. Variables that were statistically significant in the bivariate analysis were considered for the multivariable logistic regression. All variables were evaluated in the full model and retained at 0.05 by backward elimination. Analysis was performed with SAS Enterprise Guide v7.1 (SAS Institute, Cary, NC) at a significant level of 0.05. Figures were produced in GraphPad Prism v.9 (GraphPad Software, San Diego, CA).

Among the 8599 included patients with rhabdomyolysis diagnosis, the majority were male (80.9%), in late adolescence or early adulthood (15 to 20 years; 75.1%), and predominantly non-Hispanic White and non-Hispanic Black (41.7% and 36.5%, respectively). Few patients had a complex chronic condition (17.8%), and the majority had private insurance (51.2%). The geometric mean length of stay was 3.46 days (SD 0.06). Most patients presented to hospitals in the Southern region (43.7%) (Table 1).

TABLE 1

Descriptive Characteristics of Hospitalized Pediatric Patients Diagnosed With Rhabdomyolysis, With or Without Acute Renal Failure

CharacteristicOverall Cohort (N = 8599)Renal Failure = No (N = 7866)Renal Failure = Yes (N = 733)
N (%)N (%)N (%)
Male 6911 (80.9) 6235 (79.8) 676 (92.2) 
Age, y    
 <5 274 (3.2) 266 (3.4) 9 (1.2) 
 5–9 881 (10.3) 873 (11.2) 7 (1) 
 10–14 963 (11.3) 922 (11.8) 41 (5.7) 
 15–20 6394 (75.1) 5729 (73.5) 664 (92.1) 
Race    
 Non-Hispanic White 3186 (41.7) 2943 (42.2) 243 (36.5) 
 Non-Hispanic Black 2451 (32.1) 2174 (31.2) 277 (41.6) 
 Hispanic 1356 (17.8) 1243 (17.8) 112 (16.9) 
 Othera 643 (7.5) 610 (8.8) 33 (4.9) 
Payor    
 Medicaid 2808 (32.7) 2582 (32.8) 226 (30.8) 
 Private insurance 4393 (51.1) 4030 (51.2) 363 (49.5) 
 Self-pay 671 (7.8) 599 (7.6) 73 (10.0) 
 Otherb 714 (8.3) 648 (8.2) 66 (9.0) 
Region    
 Northeast 1526 (17.8) 1444 (18.4) 83 (11.3) 
 Midwest 1693 (19.7) 1526 (19.4) 167 (22.9) 
 South 3760 (43.7) 3395 (43.2) 365 (49.8) 
 West 1619 (18.8) 1501 (19.1) 117 (16) 
Month    
 January 603 (7.3) 559 (7.4) 44 (6.2) 
 February 712 (8.6) 669 (8.8) 44 (6.3) 
 March 862 (10.4) 820 (10.8) 43 (6.1) 
 April 698 (8.4) 664 (8.8) 33 (4.8) 
 May 647 (7.8) 614 (8.1) 33 (4.7) 
 June 722 (8.7) 628 (8.3) 94 (13.4) 
 July 714 (8.6) 636 (8.4) 78 (11.2) 
 August 1072 (13) 924 (12.2) 147 (21.1) 
 September 709 (8.6) 629 (8.3) 80 (11.5) 
 October 597 (7.2) 553 (7.3) 44 (6.3) 
 November 494 (6) 461 (6.1) 33 (4.7) 
 December 438 (5.3) 412 (5.4) 26 (3.7) 
Year    
 2006 1241 (14.4) 1159 (11.9) 82 (11.2) 
 2009 1672 (19.4) 1527 (14.7) 145 (19.8) 
 2012 2316 (26.9) 2157 (27.4) 159 (21.7) 
 2016 3370 (39.2) 3024 (38.4) 347 (47.3) 
CCCc, n    
7067 (82.2) 6530 (83) 537 (73.3) 
1240 (14.4) 1094 (13.9) 146 (19.9) 
2+ 291 (3.4) 242 (3.1) 50 (6.8) 
Severityd 1.43 (0.02) 1.24 (0.01) 3.46 (0.11) 
Length of stay, dayse 3.46 (0.06) 3.32 (0.06) 4.99 (0.42) 
CharacteristicOverall Cohort (N = 8599)Renal Failure = No (N = 7866)Renal Failure = Yes (N = 733)
N (%)N (%)N (%)
Male 6911 (80.9) 6235 (79.8) 676 (92.2) 
Age, y    
 <5 274 (3.2) 266 (3.4) 9 (1.2) 
 5–9 881 (10.3) 873 (11.2) 7 (1) 
 10–14 963 (11.3) 922 (11.8) 41 (5.7) 
 15–20 6394 (75.1) 5729 (73.5) 664 (92.1) 
Race    
 Non-Hispanic White 3186 (41.7) 2943 (42.2) 243 (36.5) 
 Non-Hispanic Black 2451 (32.1) 2174 (31.2) 277 (41.6) 
 Hispanic 1356 (17.8) 1243 (17.8) 112 (16.9) 
 Othera 643 (7.5) 610 (8.8) 33 (4.9) 
Payor    
 Medicaid 2808 (32.7) 2582 (32.8) 226 (30.8) 
 Private insurance 4393 (51.1) 4030 (51.2) 363 (49.5) 
 Self-pay 671 (7.8) 599 (7.6) 73 (10.0) 
 Otherb 714 (8.3) 648 (8.2) 66 (9.0) 
Region    
 Northeast 1526 (17.8) 1444 (18.4) 83 (11.3) 
 Midwest 1693 (19.7) 1526 (19.4) 167 (22.9) 
 South 3760 (43.7) 3395 (43.2) 365 (49.8) 
 West 1619 (18.8) 1501 (19.1) 117 (16) 
Month    
 January 603 (7.3) 559 (7.4) 44 (6.2) 
 February 712 (8.6) 669 (8.8) 44 (6.3) 
 March 862 (10.4) 820 (10.8) 43 (6.1) 
 April 698 (8.4) 664 (8.8) 33 (4.8) 
 May 647 (7.8) 614 (8.1) 33 (4.7) 
 June 722 (8.7) 628 (8.3) 94 (13.4) 
 July 714 (8.6) 636 (8.4) 78 (11.2) 
 August 1072 (13) 924 (12.2) 147 (21.1) 
 September 709 (8.6) 629 (8.3) 80 (11.5) 
 October 597 (7.2) 553 (7.3) 44 (6.3) 
 November 494 (6) 461 (6.1) 33 (4.7) 
 December 438 (5.3) 412 (5.4) 26 (3.7) 
Year    
 2006 1241 (14.4) 1159 (11.9) 82 (11.2) 
 2009 1672 (19.4) 1527 (14.7) 145 (19.8) 
 2012 2316 (26.9) 2157 (27.4) 159 (21.7) 
 2016 3370 (39.2) 3024 (38.4) 347 (47.3) 
CCCc, n    
7067 (82.2) 6530 (83) 537 (73.3) 
1240 (14.4) 1094 (13.9) 146 (19.9) 
2+ 291 (3.4) 242 (3.1) 50 (6.8) 
Severityd 1.43 (0.02) 1.24 (0.01) 3.46 (0.11) 
Length of stay, dayse 3.46 (0.06) 3.32 (0.06) 4.99 (0.42) 

All comparisons across groups were statistically significant with P < .001. Data may not sum to total because of a small number of missing values.

a

Asian or Pacific Islander, Native American, and other.

b

Medicare, self-pay, no charge, and other.

c

Complex chronic condition.

d

Hospitalization Resource Intensity Scores for Kids mean (SD).

e

Geometric mean (SD).

TABLE 2

Multivariable Logistic Regression Modeling Demographic Characteristics Analyzed for Associations With Outcome of Acute Renal Failure in Hospitalized Pediatric Patients With Rhabdomyolysis

CharacteristicReference GroupOdds Ratio (OR)95% Confidence LimitsP
Age 15–20 y Age < 15 y 3.55 2.53–4.98 <.0001 
Males Females 3.07 2.12–4.46 <.0001 
Non-Hispanic Black Non-Hispanic White 1.35 1.06–1.71 .0146 
Northeast region Southern region 0.53 0.38–0.73 .0001 
1 CCC No CCC 1.84 1.41–2.39 <.0001 
2+CCC No CCC 3.97 2.56–6.16 <.0001 
CharacteristicReference GroupOdds Ratio (OR)95% Confidence LimitsP
Age 15–20 y Age < 15 y 3.55 2.53–4.98 <.0001 
Males Females 3.07 2.12–4.46 <.0001 
Non-Hispanic Black Non-Hispanic White 1.35 1.06–1.71 .0146 
Northeast region Southern region 0.53 0.38–0.73 .0001 
1 CCC No CCC 1.84 1.41–2.39 <.0001 
2+CCC No CCC 3.97 2.56–6.16 <.0001 

In 2006, the rate of pediatric hospitalizations for rhabdomyolysis was 0.8 per 100 000 children, increasing to 2.8 per 100 000 children by 2016, a 218% increase (Fig 1). In 2006 and 2009, there were no significant differences in the rates of rhabdomyolysis across the 4 geographic regions. In 2012, the South and Midwest had significantly higher rates than the West (29% difference, P = .009 and 33% difference, P = .021, respectively). In 2016, rates in the South were 47% higher than the West (P < .001) and 36% higher than the Midwest (P < .001), and rates in the Northeast were 27% higher than the Midwest (P = .033) and 47% higher than the West (P < .001). There was notable seasonal variation in the incidence of rhabdomyolysis with peaks in March and August (Fig 2). The peak in March was dominated by the Northeast and Midwest, whereas the peak in August was dominated by the South and Northeast regions.

FIGURE 1

10-Year incidence of pediatric rhabdomyolysis by geographic region. Each geographic region’s rate of hospitalization for pediatric rhabdomyolysis is plotted by month and year from 2006 to 2016. The national rate is shown in red. There is a significant increase in the overall rate of hospitalizations for pediatric rhabdomyolysis, along with increasing numbers in the Northeast, Midwest, and Southern regions.

FIGURE 1

10-Year incidence of pediatric rhabdomyolysis by geographic region. Each geographic region’s rate of hospitalization for pediatric rhabdomyolysis is plotted by month and year from 2006 to 2016. The national rate is shown in red. There is a significant increase in the overall rate of hospitalizations for pediatric rhabdomyolysis, along with increasing numbers in the Northeast, Midwest, and Southern regions.

Close modal
FIGURE 2

Seasonal variation in incidence of pediatric rhabdomyolysis by geographic region. For each geographic region, the average rate of hospitalization for rhabdomyolysis per 100 000 children from 2006 to 2016 is plotted by month of the year. The red line represents the national average monthly rate of hospitalization for rhabdomyolysis per 100 000 children.

FIGURE 2

Seasonal variation in incidence of pediatric rhabdomyolysis by geographic region. For each geographic region, the average rate of hospitalization for rhabdomyolysis per 100 000 children from 2006 to 2016 is plotted by month of the year. The red line represents the national average monthly rate of hospitalization for rhabdomyolysis per 100 000 children.

Close modal

Diagnosis of acute renal failure occurred at a rate of 85.2 per 1000 hospitalized patients with rhabdomyolysis (95% confidence interval [CI] 79.3–91.1). Over time, rates increased from 6.6% in 2006 to 10.3% in 2016. Males had significantly higher rates of renal failure compared with females (9.8% versus 3.5%, respectively); similarly, the older patients compared with the younger cohort (<15 years old) had higher rates of renal failure (10.4% versus 2.7%, respectively). With concurrent renal failure diagnosis, non-Hispanic Black patients surpassed non-Hispanic White patients to represent the dominant group (41.6% versus 36.5%, respectively). Similarly, a majority of patients diagnosed with renal failure (57.2% total) were hospitalized during the summer and early fall months (June–September). The largest peak of renal failure diagnoses was in August (21.1%), with almost twice as many cases than other months. Patients with 2 or more CCCs also had a higher proportion of concurrent renal failure diagnoses than patients without any CCCs (17% versus 7.6%, respectively).

In the overall cohort of hospitalized patients with rhabdomyolysis, only 35 required reported RRT (4.1 per 1000 rhabdomyolysis patients [95% CI 2.7–5.4]); all of these patients had an associated renal failure diagnosis. Among the 733 children with rhabdomyolysis and acute renal failure diagnoses, 4.8% required RRT. Most of these patients were male and older adolescents or young adults (77.1% and 85.7%, respectively), consistent with the rates in the overall rhabdomyolysis cohort and in those without a renal failure diagnosis. Those requiring RRT had significantly higher disease severity overall (geometric mean 5.37 [SD 0.32] versus 1.43 [SD 0.02], P < .001), longer hospital length of stay (geometric mean 18.5 days [SD 2.14] versus 3.46 days [SD 0.06]), and at least 1 CCC. Though diagnoses of renal failure have increased over time, associated RRT has decreased from 12 cases in 2006 to none reported in 2016. There was no significant pattern of seasonality or regional predominance in the group of patients with RRT, and none were reported to have died.

Overall, 3 patients died (0.4 deaths per 1000 rhabdomyolysis patients [95% CI 0.0, 0.7]; none of these patients reported RRT while hospitalized, though 1 had a concurrent renal failure diagnosis. Compared with the patients who required RRT and survived, the patients who died had higher disease severity (geometric mean 8.7 [SD 0.89] versus 5.4 [SD 0.32], P < .001). Of the 3 patients who died, associated diagnoses included seizures, disseminated intravascular coagulation, cardiac arrest, shock, acute respiratory distress or failure, and nonspecific viral infection.

In multivariable modeling, the characteristics most significantly associated with renal failure diagnosis included male sex (odds ratio [OR] 3.07; 95% CI 2.12–4.98) and older patients compared with those less than 15 years old (OR 3.55; 95% CI 2.53–4.98). Presence of CCCs was also associated with diagnosis of renal failure, with 2 or more CCCs nearly double the likelihood of only 1 CCC (OR 3.97 [95% CI 2.56–6.16] versus OR 1.84 [95% CI 1.41–2.39]). Finally, patients presenting to hospitals in the Northeast region were less likely to have a diagnosis of renal failure, compared with those presenting to hospitals in the Southern region (OR 0.53; 95% CI 0.38–0.73) (Table 3).

TABLE 3

ICD-9 and ICD-10 Diagnosis and Procedure Codes Used for Definitions of Acute Renal Failure and Renal Replacement Therapy (RRT), Respectively

ICD-9ICD-10
Code Code Description Code Code Description 
Diagnosis codes 
 584.9 Acute kidney injury, acute kidney failure N17.0 Acute kidney failure with tubular necrosis 
N17.1 Acute kidney failure with acute cortical necrosis 
N17.2 Acute kidney failure with medullary necrosis 
N17.8 Other acute kidney failure 
N17.9 Acute kidney failure, unspecified 
5A1D70Z Performance of urinary filtration, intermittent, less than 6hrs per day 
5A1D80Z Performance of urinary filtration, prolonged intermittent, 6–18 h per day 
5A1D90Z Performance of urinary filtration, continuous, greater than 18 h per day 
Procedure Codes 
 38.95 Venous catheterization for renal dialysis 
 39.95 Hemodialysis 
ICD-9ICD-10
Code Code Description Code Code Description 
Diagnosis codes 
 584.9 Acute kidney injury, acute kidney failure N17.0 Acute kidney failure with tubular necrosis 
N17.1 Acute kidney failure with acute cortical necrosis 
N17.2 Acute kidney failure with medullary necrosis 
N17.8 Other acute kidney failure 
N17.9 Acute kidney failure, unspecified 
5A1D70Z Performance of urinary filtration, intermittent, less than 6hrs per day 
5A1D80Z Performance of urinary filtration, prolonged intermittent, 6–18 h per day 
5A1D90Z Performance of urinary filtration, continuous, greater than 18 h per day 
Procedure Codes 
 38.95 Venous catheterization for renal dialysis 
 39.95 Hemodialysis 

Yearly hospitalizations for rhabdomyolysis in children are relatively infrequent but have increased from 0.8 to 2.8 per 100 000 children over the past decade. The majority of rhabdomyolysis hospitalizations occur in the Southern region with higher peaks in March and August. Though renal morbidity was diagnosed in 8.5% of children requiring hospitalization for rhabdomyolysis, very few of these patients required RRT (0.41%), and death was rare (0.03%). Demographic characteristics associated with concurrent diagnosis of renal failure included male sex, adolescent, or young adult age (15 to 20 years old), non-Hispanic Black race, Southern region, and presence of 1 or more CCC.

Prior studies were mostly single or dual institution case series and unable to estimate incidence, rates over time, or regional variability. By using a large national database, we were able to establish the number of hospitalizations for rhabdomyolysis over time, noting geographic and temporal peaks; this seasonal variation remains consistent across all geographic regions. Overall, there are increased rates during March and August; prior studies were unable to observe any seasonal variability.1  Additionally, the majority of patients with concurrent diagnosis of renal failure presented in the summer months (57.2% in June to September). We were unable to make conclusions about the rates of specific etiologies of rhabdomyolysis in this cohort of patients because of the limitations of the administrative database used; however, this seasonal pattern introduces the hypothesis of possible environmental contributors. Prior studies in children report exercise or exertion and viral myositis, particularly influenza, as the most common etiologies for pediatric rhabdomyolysis.13,7,8,11,12  Our study’s reports of higher rates of rhabdomyolysis in males and adolescents and young adults is consistent with military literature on exertional rhabdomyolysis in active duty military.1921  Higher temperatures in the Southern region, paired with increased intensity of school athletic training outdoors in the late summer or early fall (August) and during spring break (March), may contribute to increased rates of heat-related illnesses, including exertional rhabdomyolysis.22  However, heat-related illnesses, including dehydration, may independently cause elevations in serum creatinine. Additionally, there is persistence of influenza and other respiratory viruses through May, which may contribute to the increase in rhabdomyolysis seen in March, though it not fully consistent with peak winter respiratory viral season.13,10,11,13 

Most pediatric patients hospitalized with rhabdomyolysis did not have CCCs, suggesting pediatric patients with rhabdomyolysis are generally healthy. Presence of CCCs did not increase risk for rhabdomyolysis hospitalization, but greater number of CCCs did confer greater likelihood of diagnosis of renal failure with rhabdomyolysis. Pediatric patients with rhabdomyolysis also have generally low disease severity, which suggests lower acuity of illness while hospitalized. Though rhabdomyolysis and associated renal failure diagnoses are increasing over time, need for renal replacement has decreased over time, suggesting either decreasing acuity of hospitalized patients versus early diagnosis and intervention preventing complications. Increased use of diagnostic testing may contribute to increased diagnoses without true increased disease incidence.

To our knowledge, this is the largest study to date to report incidence of renal morbidity in children hospitalized with rhabdomyolysis. We report a much lower rate of associated renal morbidity, at 8.5% out of 8599 hospitalized patients, compared with previously reported rates up to 50%.13,5,713,23  Our reported rate is also lower than the 14% to 48% reported for rhabdomyolysis-associated renal failure in adults.6,2428  It is important to note that definitions of renal morbidity (including acute kidney injury, acute renal insufficiency, and acute renal failure) vary across these studies, thus direct comparison is challenging. Despite a clinically significant rate of renal failure, few patients require RRT (0.41%) and death was rare (0.03%), consistent with most prior studies, with the exception of Watanabe et al who reported rates of 33% for both renal failure and death (total 18 patients).2,3,8,10,11  Further studies are needed to investigate whether this reported difference is because of the natural history of underlying etiology or because of regional differences in management.

Non-Hispanic Black patients, along with non-Hispanic White patients, made up the majority of pediatric patients hospitalized with rhabdomyolysis. A greater percentage of non-Hispanic Black children are hospitalized with rhabdomyolysis (32%) compared with the percentage of non-Hispanic Black children in the general population (13.5% in 2016).29  Previous case reports have noted an association between Black race, exertional rhabdomyolysis, and sudden death with exertion.3033  Thus, presence of sickle cell trait was hypothesized as a risk factor, since there is a greater incidence of sickle cell trait in Black newborns.34  A large scale military study of over 47 000 Black soldiers determined that sickle cell trait conferred greater adjusted risk of exertional rhabdomyolysis.35  In a similar study using an administrative database, Black race was found to be an independent risk factor for rhabdomyolysis and subsequent acute kidney injury in trauma patients, despite only having <0.1% of Black patients reporting sickle cell trait. Sickle cell trait did not independently confer greater risk rhabdomyolysis or acute kidney injury, though it may not have been powered enough to show any difference. Our study shows the majority of those with renal failure identified as non-Hispanic Black, though they were not the predominant group in the overall cohort with rhabdomyolysis. Prior studies have shown worse outcomes for Black patients presenting with trauma or sepsis caused by systemic under-recognition and under-resuscitation in this group, rather than biologic differences. Similarly, when looking at these differences in outcomes for pediatric rhabdomyolysis, the possible contribution of systemic racism and inherent biases cannot be ignored; however, further investigation is needed, as we are unable to comment on etiology for rhabdomyolysis, severity at presentation, or progression of disease while hospitalized.

This study should be viewed with several limitations in mind. Prior studies, in both pediatric and adult populations, have used varying definitions of rhabdomyolysis and acute renal failure. By using administrative billing data, patients were identified as having rhabdomyolysis or renal failure based on discharge diagnosis codes assigned by individual clinicians, which may limit the consistency in the diagnostic definition across patients, institutions, and regions in this study. The KID itself is dependent on the completeness of data reported by each participating hospital; subsequently, there are missing values reported in both race and diagnosis coding. The decision to not include patients with a secondary diagnosis of rhabdomyolysis may have limited inclusion of higher acuity patients with greater rates of renal failure and longer hospital stays. From ICD-9 to ICD-10, the specific diagnosis code for acute renal insufficiency was replaced with a more generic code for “disorder of the kidney and ureter, unspecified,” which we chose not to include. Although this kept our scope narrow to higher acuity renal failure, it may have excluded lower acuity patients with mild renal insufficiency. Whereas there was a change from ICD-9 to ICD-10, the code definition and scope for rhabdomyolysis did not change, though provider billing practices might have changed. Additionally, we were unable to account for changes in diagnostic testing patterns or thresholds for diagnosis that may influence apparent trends in incidence, since the KID does not provide specific information about diagnostic testing or individual laboratory values.36 

Rhabdomyolysis is an uncommon disease process in children with infrequent rate of hospitalizations. However, among patients who are hospitalized, a significant portion of patients can have associated renal failure, though severe complications continue to be rare. Clinicians should have a high index of suspicion for rhabdomyolysis in children during peak months (March and August), especially in geographic areas with higher incidence (South and Northwest). Further studies are needed to correlate etiology with outcomes and investigate other risk factors for morbidity and mortality.

FUNDING: No external funding.

CONFLICT OF INTEREST DISCLOSURES: The authors have indicated they have no conflicts of interest to disclose.

Dr Agharokh conceptualized and designed the study, interpreted the data, drafted the initial manuscript, and reviewed and revised the manuscript; Dr Zaniletti conceptualized and designed the study, collected data, conducted the initial analyses, interpreted the data, and reviewed and revised the manuscript; Drs Yu, Lee, Hall, Williams, and Wilson conceptualized and designed the study, interpreted the data, and critically reviewed the manuscript for important intellectual content; and all authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.

RRT

renal replacement therapy

HCUP

Healthcare Cost and Utilization Project

KID

Kids’ Inpatient Database

ICD

International Classification of Diseases, Clinical Modification

CCC

complex chronic conditions

APR-DRG

All Patient Refined Diagnosis Related Groups

1.
Mannix
R
,
Tan
ML
,
Wright
R
,
Baskin
M
.
Acute pediatric rhabdomyolysis: causes and rates of renal failure
.
Pediatrics
.
2006
;
118
(
5
):
2119
2125
2.
Chen
CY
,
Lin
YR
,
Zhao
LL
, et al
.
Clinical spectrum of rhabdomyolysis presented to pediatric emergency department
.
BMC Pediatr
.
2013
;
13
:
134
3.
Perreault
S
,
Birca
A
,
Piper
D
,
Nadeau
A
,
Gauvin
F
,
Vanasse
M
.
Transient creatine phosphokinase elevations in children: a single-center experience
.
J Pediatr
.
2011
;
159
(
4
):
682
685
4.
Luck
RP
,
Verbin
S
.
Rhabdomyolysis: a review of clinical presentation, etiology, diagnosis, and management
.
Pediatr Emerg Care
.
2008
;
24
(
4
):
262
268
5.
Agyeman
P
,
Duppenthaler
A
,
Heininger
U
,
Aebi
C
.
Influenza-associated myositis in children
.
Infection
.
2004
;
32
(
4
):
199
203
6.
Chen
CY
,
Lin
YR
,
Zhao
LL
,
Yang
WC
,
Chang
YJ
,
Wu
HP
.
Clinical factors in predicting acute renal failure caused by rhabdomyolysis in the ED
.
Am J Emerg Med
.
2013
;
31
(
7
):
1062
1066
7.
Gelbart
B
,
DeMarco
R
,
David Hussey
A
, et al
.
Rhabdomyolysis in a tertiary PICU: a 10-year study
.
Pediatr Crit Care Med
.
2018
;
19
(
1
):
e51
e57
8.
Gardner
HM
,
Askenazi
DJ
,
Hoefert
JA
,
Helton
A
,
Wu
CL
.
Acute kidney injury among children admitted with viral rhabdomyolysis
.
Hosp Pediatr
.
2021
;
11
(
8
):
878
885
9.
Hummel
K
,
Gregory
A
,
Desai
N
,
Diamond
A
.
Rhabdomyolysis in adolescent athletes: review of cases
.
Phys Sportsmed
.
2016
;
44
(
2
):
195
199
10.
Rosenberg
T
,
Heitner
S
,
Scolnik
D
,
Levin Ben-Adiva
E
,
Rimon
A
,
Glatstein
M
.
Outcome of benign acute childhood myositis: the experience of 2 large tertiary care pediatric hospitals
.
Pediatr Emerg Care
.
2018
;
34
(
6
):
400
402
11.
Watanabe
T
.
Rhabdomyolysis and acute renal failure in children
.
Pediatr Nephrol
.
2001
;
16
(
12
):
1072
1075
12.
Watemberg
N
,
Leshner
RL
,
Armstrong
BA
,
Lerman-Sagie
T
.
Acute pediatric rhabdomyolysis
.
J Child Neurol
.
2000
;
15
(
4
):
222
227
13.
Zepeda-Orozco
D
,
Ault
BH
,
Jones
DP
.
Factors associated with acute renal failure in children with rhabdomyolysis
.
Pediatr Nephrol
.
2008
;
23
(
12
):
2281
2284
14.
Yang
Y
,
Carter
LP
,
Cook
RE
,
Paul
E
,
Schwartz
KR
.
A case of exertional rhabdomyolysis: a cheer for standardizing inpatient management and prevention
.
Hosp Pediatr
.
2016
;
6
(
12
):
753
756
15.
U.S. Census Bureau
.
Population and housing unit estimates
.
Available at: https://www.census.gov/programs- surveys/popest.html. Accessed September 1, 2021
16.
Feudtner
C
,
Christakis
DA
,
Connell
FA
.
Pediatric deaths attributable to complex chronic conditions: a population-based study of Washington State, 1980-1997
.
Pediatrics
.
2000
;
106
(
1 Pt 2
):
205
209
17.
Feudtner
C
,
Feinstein
JA
,
Zhong
W
,
Hall
M
,
Dai
D
.
Pediatric complex chronic conditions classification system version 2: updated for ICD-10 and complex medical technology dependence and transplantation
.
BMC Pediatr
.
2014
;
14
:
199
18.
U.S. Census Bureau
.
2010 census regions and divisions of the United States
.
19.
Armed Forces Health Surveillance Branch
.
Update: exertional rhabdomyolysis, active component, U.S. Army, Navy, Air Force, and Marine Corps, 2011-2015
.
MSMR
.
2016
;
23
(
3
):
21
24
20.
Armed Forces Health Surveillance Bureau
.
Update: exertional rhabdomyolysis, active component, U.S. Armed Forces, 2012-2016
.
MSMR
.
2017
;
24
(
3
):
14
18
21.
Armed Forces Health Surveillance Bureau
.
Update: exertional rhabdomyolysis, active component, U.S. Armed Forces, 2014-2018
.
MSMR
.
2019
;
26
(
4
):
21
25
22.
Choudhary
E
,
Vaidyanathan
A
;
Centers for Disease Control and Prevention
.
Heat stress illness hospitalizations — environmental public health tracking program, 20 states, 2001–2010
.
Available at: https://www.cdc.gov/mmwr/preview/mmwrhtml/ss6313a1.htm. Accessed August 17, 2021
23.
Chamberlain
MC
.
Rhabdomyolysis in children: a 3-year retrospective study
.
Pediatr Neurol
.
1991
;
7
(
3
):
226
228
24.
McKenna
MC
,
Kelly
M
,
Boran
G
,
Lavin
P
.
Spectrum of rhabdomyolysis in an acute hospital
.
Ir J Med Sci
.
2019
;
188
(
4
):
1423
1426
25.
McMahon
GM
,
Zeng
X
,
Waikar
SS
.
A risk prediction score for kidney failure or mortality in rhabdomyolysis
.
JAMA Intern Med
.
2013
;
173
(
19
):
1821
1828
26.
Melli
G
,
Chaudhry
V
,
Cornblath
DR
.
Rhabdomyolysis: an evaluation of 475 hospitalized patients
.
Medicine (Baltimore)
.
2005
;
84
(
6
):
377
385
27.
Premru
V
,
Kovač
J
,
Ponikvar
R
.
Use of myoglobin as a marker and predictor in myoglobinuric acute kidney injury
.
Ther Apher Dial
.
2013
;
17
(
4
):
391
395
28.
Simpson
JP
,
Taylor
A
,
Sudhan
N
,
Menon
DK
,
Lavinio
A
.
Rhabdomyolysis and acute kidney injury: creatine kinase as a prognostic marker and validation of the McMahon Score in a 10-year cohort: a retrospective observational evaluation
.
Eur J Anaesthesiol
.
2016
;
33
(
12
):
906
912
30.
Anzalone
ML
,
Green
VS
,
Buja
M
,
Sanchez
LA
,
Harrykissoon
RI
,
Eichner
ER
.
Sickle cell trait and fatal rhabdomyolysis in football training: a case study
.
Med Sci Sports Exerc
.
2010
;
42
(
1
):
3
7
31.
Harrelson
GL
,
Fincher
AL
,
Robinson
JB
.
Acute exertional rhabdomyolysis and its relationship to sickle cell trait
.
J Athl Train
.
1995
;
30
(
4
):
309
312
32.
Kark
JA
,
Posey
DM
,
Schumacher
HR
,
Ruehle
CJ
.
Sickle-cell trait as a risk factor for sudden death in physical training
.
N Engl J Med
.
1987
;
317
(
13
):
781
787
33.
Makaryus
JN
,
Catanzaro
JN
,
Katona
KC
.
Exertional rhabdomyolysis and renal failure in patients with sickle cell trait: is it time to change our approach?
Hematology
.
2007
;
12
(
4
):
349
352
34.
Ojodu
J
,
Hulihan
MM
,
Pope
SN
,
Grant
AM
.
Incidence of sickle cell trait — United States, 2010
.
35.
Nelson
DA
,
Deuster
PA
,
Carter
R
III
,
Hill
OT
,
Wolcott
VL
,
Kurina
LM
.
Sickle cell trait, thabdomyolysis, and mortality among U.S. Army soldiers
.
N Engl J Med
.
2016
;
375
(
5
):
435
442
36.
Müskens
JLJM
,
Kool
RB
,
van Dulmen
SA
,
Westert
GP
.
Overuse of diagnostic testing in healthcare: a systematic review
.
BMJ Qual Saf
.
2022
;
31
(
1
):
54
63