- Research article
- Open Access
- Open Peer Review
Multi-parametric muscle and fat correlation of computed tomography parameters to outcomes in a total hip arthroplasty population
BMC Musculoskeletal Disordersvolume 19, Article number: 4 (2018)
Cross-sectional imaging is not currently used in planning Total Hip Arthroplasty (THA). The aim of our study is to determine correlations between CT parameters and outcomes following THA.
A prospective registry of patients who underwent total joint arthroplasty was reviewed for patients who: (1) underwent THA, (2) had a CT between 1 year before and 6 months after surgery, and (3) completed perioperative WOMAC and Harris Hip Score (HHS) questionnaires. Two readers measured CT parameters, yielding mean Hounsfield Units, area, average diameter, and perimeter of the psoas major, gluteus medius and minimus muscles. A segmentation algorithm determined visceral and subcutaneous fat area, and waist circumference. ICC was calculated for each measurement to examine inter-reader agreement. Regression analyses were performed to select measurements with most impact on outcome scores.
Twenty-eight patients met inclusion criteria (17 female, 11 male), having mean (+/− standard deviation) age of 54.4 +/− 14.8 years and BMI 29.0 +/− 6.3 kg/m2. Correlations were found between HHS and age (0.650, p = 0.018), height (−1.263, p = 0.009), visceral-to-subcutaneous fat area ratio at the psoas level (0.511, p = 0.018), and waist circumference at the psoas level (1.759, p = 0.002). Inter-reader analysis showed ICC > 0.850 for all measurements.
Age and height, as well as CT-derived visceral-to-subcutaneous fat area ratio and waist circumference significantly correlate with postsurgical HHS scores following THA. Our study suggests that parameters derived from cross-sectional CT imaging can be useful additional preoperative planning tool for THA.
Almost 700,000 Total Joint Arthroplasties (TJA) are performed in the United States each year, which include over 230,000 Total Hip Arthroplasties (THA) . It is estimated that 15.1% of the THAs performed each year are revision procedures rather than primary operations, the economic impact of which is profound . From 1997 to 2003, 19% of all Medicare spending on THA was consumed by revision operations . THA is a widely used procedure, and the prevalence is projected to increase to as many as 668,700 procedures in 2030, with 96,700 (14.4%) of those estimated to be revision THA . Methods to identify patients who will most likely benefit and likewise identify patients who will fail THA are immensely important because they can potentially aid in optimizing healthcare spending.
Cross-sectional imaging modalities, such as Magnetic Resonance Imaging (MRI) and Computed Tomography (CT), are not commonly used in preoperative screening of prospective THA patients. In 2015, Wiater et al. retrospectively examined a registry of patients who underwent reverse total shoulder arthroplasty (RTSA) and found significant correlations between radiographic parameters from MRI and postoperative outcome assessments .
CT measurement of the psoas major muscle has been used as a surrogate measure of sarcopenia . It has been shown in some surgical and trauma populations that increased rates of overall morbidity [4,5,6], mortality , and decreased quality-of-life  significantly correlate with decreased psoas major area, volume, and density. The gluteal (abductor) muscles are considered highly important for overall hip function, especially for function after THA . Due to superior soft tissue contrast and spatial resolution, CT examination of the pelvis can be used to evaluate both the abductor and psoas major muscle groups. In addition, evaluation of muscle size, muscle density, and fat segmentation with an analysis of fat fractions can be performed using independent software in a reproducible fashion.
Overall, there is a knowledge gap as to how these CT parameters ultimately can be related to outcomes of THA. No previous study has evaluated multi-compartment tissue analysis in the same setting and its relationship to outcomes of THA. The aim of our exploratory study was therefore to determine correlations between CT parameters and outcomes following THA.
Patients were selected from a prospectively constructed registry of TJA patients at our institution. A retrospective review of this registry was performed under a HIPPA waiver and as part of an Institutional Review Board-approved study. This review yielded 802 TJA patients, 318 of whom underwent THA between 2006 and 2014.
At their initial visit, 6 months postoperatively, and 12 months postoperatively, THA patients completed the Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) Hip Score and Harris Hip Score (HHS) questionnaires in addition to normal follow-up.
We included patients who : underwent THA , had a pelvic CT study performed for any reason in the time 1 year before or 6 months after the surgery, and  completed at least one set of postoperative questionnaires in addition to the initial visit questionnaires.
Image analysis and reader training
All CT images were reconstructed as 3 mm slices in a soft tissue window for evaluation. Images were processed using an independent software program, Aquarius intuition (TeraRecon, Foster City, CA).
Two second-year medical students (MH, RB) served as readers, and independently measured parameters at two levels. Each reader was trained by a board-certified, fellowship trained radiologist (AC) to identify and analyze defined slices and measure selected regions of interest (ROIs) (See Slice and ROI Selection). This training consisted of two sessions, each two-hours in length, where the readers were taught to identify the defined slices and measure the selected ROIs on CT images from non-study patients. Readers then practiced this workflow under the supervision of the radiologist.
Additionally, 10% of the study sample (3 CTs) were used as a training set. These images were measured by the two readers using separate but identical workstations, working simultaneously and cooperatively to ensure a standardized workflow. Subsequently, the readers collected measurements using identical workstations, independent of one another on the remaining 90%, such that each CT was eventually measured twice, once by each reader.
Slice and ROI selection
Two slices of interest were defined, and on each slice ROIs were selected. The psoas slice was defined as the first axial slice superior to the sacroiliac joint, but which does not visualize the sacroiliac joint (Fig. 1a and b). The ROIs on the psoas slice included the psoas major muscles, both ipsilateral and contralateral to surgery. These muscles were selected as ROIs because cross-sectional psoas area measured from axial CT has been used as a marker for core muscle mass and sarcopenia, and has been shown to be correlated with postoperative functional outcomes after THA [4, 10]. The gluteal slice was defined as the axial slice demonstrating fully open S1 anterior foramina (Fig. 1c and d). The ROIs on the gluteal slice included the gluteus medius-minimus complex, again ipsilateral and contralateral to surgery. These muscles were identified as ROIs because they contribute to the overall musculature surrounding the hip and gluteus medius atrophy has been found to predict limping after THA at 6 months postoperatively [9, 11].
Measurement of ROIs
On the psoas slice, three functions were performed: measurement of the psoas major ROI ipsilateral and contralateral to surgery and an automated fat-analysis segmentation algorithm. On the gluteal slice, three functions were performed: measurement of the gluteus medius-minimus ROI ipsilateral and contralateral to surgery and the same fat-analysis segmentation algorithm.
Measurement of each muscle was accomplished by manually constructing a polygon around the indicated muscle, yielding four parameters: mean Hounsfield units (Mean), cross-sectional area (Area), mean diameter (Diameter), and Circumference (Fig. 1a and c).
The fat-analysis segmentation algorithm also yielded four parameters: visceral fat area, subcutaneous fat area, ratio of visceral fat to total fat areas (Area Ratio, i.e. visceral fat area divided by the sum of visceral fat and subcutaneous fat areas), and the waist circumference of the patient (Outer Circumference) (Fig. 1b and d).
A sensitivity analysis was performed between the study group and the overall THA group to determine whether the study sample differed significantly from the excluded registry population in demographic or clinical factors. Inter-reader agreement was assessed using the intraclass correlation coefficient (ICC) value for each measurement. Average measurements between the two readers were used for this analysis.
Additionally, univariate and multiple regression analyses were performed to determine the association between measurements and postsurgical outcome scores. The CT measurements and demographic data were used as independent variables and change in WOMAC and HHS were response variables. Change in WOMAC and HHS was defined as the difference between 12-month and initial scores.
In the multiple regression analysis, all CT measurements together with age, BMI, height, weight, and sex were initially included. A stepwise selection algorithm was then used to select the independent variables that most affected the response variables based on the F statistic. Using the selected independent variables, a multivariate model was constructed which was used to predict change in WOMAC and HHS.
All p-values of less than 0.05 were considered to be statistically significant. All statistical analysis was performed using Statistical Analysis Software (SAS Institute, Cary, NC).
Of the 318 THA patients in the registry, 13.2% (42/318) had a CT study performed in the indicated time frame. Of those 42 patients who underwent CT imaging, 66.7% (28/42, 17 females, 11 males) had a complete set of postoperative questionnaires (Fig. 2). Mean (+/− standard deviation) age was 54.4 +/− 14.8 years, height was 165.9 +/− 9.0 cm, weight was 80.0 +/− 19.3 kg, and mean BMI was 29.0 +/− 6.3 kg/m2. Mean (+/− standard deviation) initial WOMAC score was 62.7 +/− 18.4 and initial HHS was 38.4 +/− 17.8.
The study population was found to be similar to the overall THA group, with height being slightly different (Table 1). Analysis of agreement between the two readers showed an ICC > 0.85 for all measurements (Table 2). No significant correlation was found between any demographic or CT parameters and change in WOMAC Hip Score.
Significant univariate correlations were found between HHS and age (0.650, p = 0.018), height (−1.263, p = 0.009), Area Ratio at the psoas slice (Psoas Area Ratio, 0.511, p = 0.018), and Outer Circumference at the psoas slice (Psoas Outer Circumference, 1.759, p = 0.002) (Table 3, Fig. 3). Additionally, a non-significant trend was found between HHS and the mean Hounsfield units of the gluteus minimus-medius ipsilateral to surgery (Ipsilateral Gluteal Mean, 0.359, p = 0.108).
Correlations are reported such that a one unit increase in a given CT measurement is associated with the reported increase in HHS. For example, an increase in height of 1 cm is associated with a 1.263 decrease in HHS, and a one Hounsfield unit increase in Ipsilateral Gluteal Mean is associated with a 0.359 increase in HHS.
The resulting multivariate model that was produced was:
where PAR is Psoas Area Ratio, POC is Psoas Outer Circumference, and IGM is Ipsilateral Gluteal mean. When used to predict change in HHS for each patient, the model performed with R2 = 0.5858 (Fig. 4).
To the best of the authors’ knowledge, this is the first attempt at evaluation of pelvis CT muscle and fat parameters and characterization of their bearing on THA outcomes. The software proved to be reliable in deriving muscle area, density, and circumference from CT images as confirmed by the inter-observer performance. The fat ratios are easily evaluated, yielding reproducible data regarding subcutaneous and visceral fat proportions. Outer waist circumference was also reproducibly measured.
While no significant correlation was found between any CT measurement and change in WOMAC Hip Score, it is apparent that some CT parameters correlated with the postoperative HHS assessment tool. However, some of the correlations we observed are not what would be expected. For instance, larger waist circumferences and higher visceral fat fractions, surrogate markers of obesity, have been shown to be associated with negative surgical outcomes [12, 13], but in this study, they were found to be associated with improved outcome after THA. We believe this anomaly may be due to the limitations of this study, such as its retrospective nature and relatively small sample size. However, a sensitivity analysis showed that the data adequately represented the overall THA study group.
Our data reflects a preliminary exploration of this topic and future studies should include prospective enrollment and a standardized multivariate statistical approach including intra-reader performance. Quantitative information obtained from preoperative imaging not only holds diagnostic value, it may also aid in prognostication. This information may prove beneficial in preoperative patient counseling and might aid preoperative and postoperative decision-making by identifying subpopulations of patients who may benefit by therapy aimed at improving muscle properties before undergoing THA. A follow-up study should involve prospective recruitment of patients with focus on uniform imaging procedures and consistent collection of outcome data in order to ensure a large sample size. While MRI is also a potential technique for the above measurements, in this retrospective evaluation CTs were available for these patients who uniformly had hip scoring performed prospectively. In addition, the Aquarius software used in this study is more useful for segmentation using CT data. Although MRI can provide fat fraction analysis with recent advanced Dixon techniques, it happens to be a more expensive study and the imaging may be limited by metal artifacts due to hip prosthesis. There has been recent encouraging work on muscle and fat composition using MRI with qualitative and quantitative muscle and fat fraction analyses yielding meaningful correlations of muscle volume and intramuscular steatosis with muscle power and post-surgical outcomes [14,15,16,17]. Similar work could be potentially performed in THA population using metal reduction techniques. Lastly, volumetric multi-slice evaluation is now available with recent software upgrades. In a future study, it could be employed using a similar dataset to determine whether correlations are markedly different from those presented in this study.
To conclude, multi-compartment tissue parameters derived from cross-sectional CT imaging can be useful as an additional preoperative planning tool for THA and establish a foundation for future prospective studies in this domain to expand validity of our preliminary results.
Harris Hip Score
- ICC :
Intra-class Correlation Coefficient
Magnetic Resonance Imaging
Total Hip Arthroplasty
- TJA :
Total Joint Arthroplasty
Western Ontario and McMaster Universities Osteoarthritis Index
Kurtz S, Ong K, Lau E, Mowat F, Halpern M. Projections of primary and revision hip and knee arthroplasty in the United States from 2005 to 2030. J Bone Joint Surg Am. 2007;89:780–5.
Ong KL, Mowat FS, Chan N, Lau E, Halpern MT, Kurtz SM. Economic burden of revision hip and knee arthroplasty in Medicare enrollees. Clin Orthop Relat Res. 2006;446:22–8.
Wiater BP, Koueiter DM, Maerz T, et al. Preoperative deltoid size and fatty infiltration of the deltoid and rotator cuff correlate to outcomes after reverse total shoulder arthroplasty. Clin Orthop Relat Res. 2015;473:663–73.
Jones KI, Doleman B, Scott S, Lund JN, Williams JP. Simple psoas cross-sectional area measurement is a quick and easy method to assess sarcopenia and predicts major surgical complications. Color Dis. 2015;17:O20–6.
Joglekar S, Asghar A, Mott SL, et al. Sarcopenia is an independent predictor of complications following pancreatectomy for adenocarcinoma. J Surg Oncol. 2015;111:771–5.
Cornet M, Lim C, Salloum C, et al. Prognostic value of sarcopenia in liver surgery. J Visc Surg. 2015;152:297–304.
Englesbe MJ, Patel SP, He K, et al. Sarcopenia and mortality after liver transplantation. J Am Coll Surg. 2010;211:271–8.
Fairchild B, Webb TP, Xiang Q, Tarima S, Brasel KJ. Sarcopenia and frailty in elderly trauma patients. World J Surg. 2015;39:373–9.
Lachiewicz PF. Abductor tendon tears of the hip: evaluation and management. J Am Acad Orthop Surg. 2011;19:385–91.
Wakabayashi H, Watanabe N, Anraku M, Oritsu H, Shimizu Y. Pre-operative psoas muscle mass and post-operative gait speed following total hip arthroplasty for osteoarthritis. J Cachexia Sarcopenia Muscle. 2016;7:95–6.
Nankaku MPP, Tsuboyama TPM, Aoyama TPM, Kuroda YPM, Ikeguchi RPM, Matsuda SPM. Preoperative gluteus medius muscle atrophy as a predictor of walking ability after total hip arthroplasty. Phys Ther Res. 2016;19:8–12.
Kovesdy CP, Czira ME, Rudas A, et al. Body mass index, waist circumference and mortality in kidney transplant recipients. Am J Transplant. 2010;10:2644–51.
Aquina CT, Rickles AS, Probst CP, et al. Visceral obesity, not elevated BMI, is strongly associated with incisional hernia after colorectal surgery. Dis Colon rectum 2015;58:220-227.
Marcon M, Berger N, Manoliu A, et al. Normative values for volume and fat content of the hip abductor muscles and their dependence on side, age and gender in a healthy population. Skelet Radiol. 2016;45:465–74.
Lang T, Cauley JA, Tylavsky F, et al. Computed tomographic measurements of thigh muscle cross-sectional area and attenuation coefficient predict hip fracture: the health, aging, and body composition study. J Bone Miner Res. 2010;25:513–9.
Inacio M, Ryan AS, Bair WN, Prettyman M, Beamer BA, Rogers MW. Gluteal muscle composition differentiates fallers from non-fallers in community dwelling older adults. BMC Geriatr. 2014;14:37.
Lee S, Lucas RM, Lansdown DA, et al. Magnetic resonance rotator cuff fat fraction and its relationship with tendon tear severity and subject characteristics. J Shoulder Elb Surg. 2015;24:1442–51.
The authors acknowledge and thank Robert W. Bucholz, MD, who assembled the TJA registry over many years.
Only institutional sources of funding were used for this project.
Availability of data and materials
The datasets used and analysed in this study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
This retrospective study was approved by the Institutional Review Board of UT Southwestern Medical center with reference number STU 042014–073, and informed consent was waived by the IRB due to the retrospective nature of the study. The TJA registry examined in the present study was was created in concordance with and after approval by the Institutional Review Board of UT Southwestern Medical Center with reference number STU 122010–071, and written, informed consent was obtained prior to inclusion in the registry.
Consent for publication
AC is a consultant for Siemens regarding MSK CAD in areas not related to this study. DH and AC also serve on the editorial board of this journal, BMC Musculoskeletal Disorders. The authors declare no other competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.