Abdominal subcutaneous fat quantification in obese patients from limited field-of-view MRI data.

Different types of adipose tissue can be accurately localized and quantified by tomographic imaging techniques (MRI or CT). One common shortcoming for the abdominal subcutaneous adipose tissue (ASAT) of obese subjects is the technically restricted imaging field of view (FOV). This work derives equations for the conversion between six surrogate measures and fully segmented ASAT volume and discusses the predictive power of these image-based quantities. Clinical (gender, age, anthropometry) and MRI data (1.5 T, two-point Dixon sequence) of 193 overweight and obese patients (116 female, 77 male) from a single research center for obesity were analyzed retrospectively. Six surrogate measures of fully segmented ASAT volume (VASAT) were considered: two simple ASAT lengths, two partial areas (Ap-FH, Ap-ASIS) and two partial volumes (Vp-FH, Vp-ASIS) limited by either the femoral heads (FH) or the anterior superior iliac spine (ASIS). Least-squares regression between each measure and VASAT provided slope and intercept for the computation of estimated ASAT volumes (V~ASAT). Goodness of fit was evaluated by coefficient of determination (R2) and standard deviation of percent differences (sd%) between V~ASAT and VASAT. Best agreement was observed for partial volume Vp-FH (sd% = 14.4% and R2 = 0.78), followed by Vp-ASIS (sd% = 18.1% and R2 = 0.69) and AWFASIS (sd% = 23.9% and R2 = 0.54), with minor gender differences only. Other estimates from simple lengths and partial areas were moderate only (sd% > 23.0% and R2 < 0.50). Gender differences in R2 generally ranged between 0.02 (dven) and 0.29 (Ap-FH). The common FOV restriction for MRI volumetry of ASAT in obese subjects can best be overcome by estimating VASAT from Vp-FH using the equation derived here. The very simple AWFASIS can be used with reservation.

Introduction

Obesity is one of the major healthcare problems of the twenty-first century with a steadily increasing prevalence over the last decades[1]. In the United States, for example, data from the years 2013 to 2016 have shown that more than two-third of the adults were either obese (38.9%, BMI: 30–40 kg/m2) or severely obese (7.6%, BMI: > 40 kg/m2)[2]. Overweight and obesity are closely associated with an increased overall mortality and morbidity often caused by metabolic or cardiovascular diseases. Direct and indirect treatment costs generate an enormous socio-economic burden on society[3-5].

Abdominal subcutaneous (ASAT) and visceral adipose tissue (VAT) contribute differently to metabolic homeostasis. Fat amounts and their distribution are important risk factors in the pathogenesis of cardiometabolic diseases[6]. VAT is an extremely dyslipidemic and atherogenic fat depot due to its endocrine activity[7]. In contrast, accumulation of ASAT is an independent predictor of lower cardiometabolic mortality[8]. Metabolic complications may arise when ASAT fails to expand and store fat leading to ectopic fat deposition with subsequent lipotoxicity and insulin resistance[9,10]. In view of limited health resources, there is a considerable interest in identifying robust predictors for cardiometabolic risk stratification and monitoring of clinical outcome in obese patients. The amount and distribution of ASAT and VAT are currently among the most promising parameters for that purpose. A weight-loss study by Mayo-Smith et al. has found that the relative volume reduction in ASAT was larger than in VAT, suggesting that ASAT might better reflect the nutritional status[11].

Imaging techniques enable a non-invasive identification of individual fat compartments. MRI provides an excellent soft-tissue contrast and is characterized by the absence of ionizing radiation and should therefore be preferred for younger patients or longitudinal studies with repeated assessments[12-14]. Since the beginning of MRI-based fat quantification in the mid-1980s, methods have been validated post mortem against sections of human and animal cadavers[15-18]. A number of studies have argued that MRI should be preferred over CT in terms of safety[19-22].

For subjects with a higher degree of obesity, it is common that the field of view (CT or MRI) is technically not large enough to image the whole cross section. With most patients being examined in supine position, this will typically exclude ventral or lateral subcutaneous adipose tissue and prevent exact quantification. This shortcoming, although referred to in the literature, has not been resolved so far. In our study, for example, the field of view was not large enough for about 82% of our patients.

The objective of this study was therefore to determine the power of various surrogate measures (distances, fat areas and partial fat volumes) for the prediction of total ASAT volume in overweight or obese adults with missing ('out-of-view') fat data. Gender-specific (linear) equations and measures of agreement are provided to identify more (or less) suitable approaches of ASAT volume estimation.

Material and methods

Study design and data

This retrospective IRB-approved study was performed at a single center–the Integrated Research and Treatment Center for Adiposity Diseases, University Medicine Leipzig, Leipzig, Germany. Data collection, analysis and publication were approved by the Institutional Review Board (IRB) of the Leipzig University Faculty of Medicine, Leipzig, Germany (reference numbers 283/11-ff, 284/10-ff, 363/10-ff, 363/11-ff) and informed consent was obtained from all subjects. All methods were carried out in accordance with relevant guidelines and regulations (Declaration of Helsinki).

A total of 193 patients (116 females and 77 males) were selected, primarily with respect to age (at least 18 years) and BMI (> 24.5 kg/m2). Patients with visibly missing parts of the abdominal subcutaneous tissue on any of the image slices were excluded. Figure 1 shows a histogram of the BMI distribution for all 193 patients along with a Gaussian fit to the data. Further details of the study population are given in Supplementary Table S1.

Figure 1

Histogram of BMI distribution for all 193 patients (range 24.8–41.2 kg/m2). Triangles at integer BMI values reflect number of patients (frequency) with actual BMI ± 0.5 kg/m2 around that value. Dotted curve indicates Gaussian fit to BMI data (32.1 ± 3.7 kg/m2).

MR imaging and analysis

The patient was examined in supine position using a standard 1.5 T system (Achieva XR, Philips Healthcare, Best, Netherlands) and the integrated whole-body coil for signal reception. Fat-sensitive imaging was based on an axial two-point Dixon sequence with two stacks of 25 10-mm-thick slices (0.5 mm interslice gap) covering the whole abdominal cavity. The main imaging parameters are listed in Supplementary Table S2. Image analysis was performed with a custom made software[23-25] that was developed using a predefined framework for radiological image analysis[26]. Surrogate measures of full ASAT volume were thickness of the abdominal wall fat or hip girdle fat (AWF and HGF), partial ASAT areas (Ap-FH and Ap-ASIS) and partial volumes (Vp-FH and Vp-ASIS). Nomenclature and definitions are given in Fig. 2 and Table 1.

Figure 2

Overview of surrogate parameters. left: length (thickness) of abdominal wall fat (AWFASIS) and hip girdle fat (HGFASIS) at the level of the anterior superior iliac spine (ASIS). middle: partial ASAT area within rectangles bound by the femoral heads or ASIS (Ap -FH, Ap-SIAS). right: partial ASAT volume (Vp-FH, Vp-ASIS) is the sum over all local ASAT volumes (partial ASAT areas Ap × interslice spacing of 10.5 mm) between diaphragm and pelvic floor.

Table 1

Nomenclature of surrogate parameters.

Quantity	Symbol	Definition
Reference	VASAT	ASAT volume measured from segmented ASAT areas of all slices between diaphragm and pelvic floor (typically 40–50 slices with effective spacing of 10.5 mm)
Length (thickness)	V~ASAT(AWFASIS)	ASAT volume estimated by length (thickness) of abdominal wall fat (AWFASIS) or hip girdle fat (HGFASIS) on single slice (at the level of the anterior superior iliac spine, ASIS)
	V~ASAT(HGFASIS)
Partial area	V~ASAT(Ap-ASIS)	ASAT volume estimated from partial ASAT area on single slice within rectangular box bound by either centers of femoral heads (Ap-FH) or ASIS (Ap-ASIS)
	V~ASAT(Ap-FH)
Partial volume	V~ASAT(Vp-ASIS)	ASAT volume estimated from partial ASAT areas of all slices within rectangular box bound by either centers of femoral heads (Vp-FH) or anterior superior iliac spines (Vp-ASIS)
	V~ASAT(Vp-FH)

AWF maximum paramedian distance between outer rectus sheath and skin surface, HGF maximum perpendicular distance between tangent to iliac plate and skin surface.

Statistical analysis

A least-squares regression (fit) between surrogate measure S and total ASAT volume was performed according to the equationwith resulting fit parameters slope (ms) and intercept (bs). For a given surrogate measurement, total ASAT volume estimates, denoted with a tilde, for example, V~ASAT(dven), V~ASAT(Ap-ASIS) or V~ASAT(Vp-FH), can then be computed with these conversion parameters. The predictive value of a surrogate measure was quantified with the coefficient of determination R2 and the standard deviation sd% of the percent differences between estimated and fully measured volume (V~ASAT – VASAT)/VASAT × 100%. All statistical analyses were performed with SPSS version 24 (IBM, Armonk, NY).

Results

Figure 3 shows gender-specific scatter plots and linear fits between fully measured ASAT volume and two selected surrogate measures: partial ASAT volume bound by femoral heads (Vp-FH) and simple abdominal wall fat at the level of the anterior superior iliac spine (AWFASIS). There is generally good to moderate agreement of the data with the linear fits. The coefficients of determination R2 are considered as good for Vp-FH and moderate for AWFASIS. For both measures, the mere R2 values were slightly higher for males.

Figure 3

Gender-specific scatter plots and linear fits (solid lines) between reference ASAT volume (VASAT) and two surrogate measures. (a, b) ASAT volume estimated from partial ASAT areas of all axial (10-mm thick) abdominal slices within rectangular box bound by centers of femoral heads (Vp-FH). (c, d) simple length (thickness) of abdominal wall fat (AWFASIS) on single slice at the level of the anterior superior iliac spine (ASIS). Plots are annotated with corresponding coefficient of determination R2.

The resulting slope ms and intercept bs were used to compute estimated total ASAT volumes from the respective surrogate measures. Supplementary Fig. S1 shows the corresponding Bland–Altman plots between estimated and fully measured total ASAT volumes. Correspondingly, the standard deviation (and limits of agreement) are smaller for Vp-FH (vs. AWFASIS) and also smaller for males (vs. females) for both surrogates.

Each difference was then divided by the respective reference ASAT volume to account for the different volume scales for females and males. The resulting standard deviations of the percent differences are plotted in Fig. 4 for all surrogate measures. At about 14%, Vp-FH represents the best surrogate followed by Vp-ASIS at around 17.5%. The sd% values of partial areas and simple distances were all above 20% with the next best agreement for AWFASIS. Gender differences were not systematic and on the order of a few percents except for Vp-FH.

Figure 4

Standard deviations [%] sd% of the mean differences between ASAT volume estimated from surrogate parameters (V~ASAT) and fully measured one (VASAT) for both genders individually: simple lengths (AWFASIS, HGFASIS), partial areas (Ap-FH, Ap-ASIS) and partial volumes (Vp-FH, Vp-ASIS). ASIS anterior superior illiac spine, FH femoral heads.

Table 2 provides a summary of all results and allows for a detailed inspection of both measures of agreement (R2 and sd%) and gender differences. Vp-FH happens to be the surrogate with both smallest sd% and highest coefficient of determination as well as no substantial difference between females and males. By design, the reported ms and bs values represent the best fit to the respective data. In practice, however, they should only be used with the knowledge of the corresponding measures of agreement given here.

Table 2

Gender-specific measures of agreement and regression parameters between reference ASAT volume and surrogate ASAT measures.

Surrogate	Males (77)				Females (116)				Total (193)
	rp (R2)	sd% [%]	mS	bs [L]	rp (R2)	sd% [%]	ms	bs [L]	rp (R2)	sd% [%]	ms	bs [L]
AWFASIS	0.69 (0.48)	21.2	0.174 L/mm	6.36	0.71 (0.50)	23.1	0.178 L/mm	4.36	0.73 (0.54)	23.9	0.198 L/mm	4.74
HGFASIS	0.55 (0.30)	23.2	0.127 L/mm	4.58	0.62 (0.39)	26.5	0.162 L/mm	2.86	0.69 (0.47)	24.6	0.137 L/mm	3.95
Ap-FH	0.50 (0.25)	24.1	0.127 L/cm2	6.05	0.73 (0.54)	21.0	0.198 L/cm2	2.05	0.69 (0.47)	23.4	0.164 L/cm2	3.70
Ap-ASIS	0.60 (0.36)	22.4	0.097 L/cm2	5.13	0.57 (0.32)	26.4	0.090 L/cm2	4.49	0.69 (0.47)	24.5	0.107 L/cm2	3.97
Vp-FH	0.84 (0.70)	14.3	5.303	2.19	0.88 (0.77)	14.0	6.297	0.06	0.88 (0.78)	14.4	5.774	1.07
Vp-ASIS	0.79 (0.63)	16.9	3.430	2.98	0.77 (0.60)	18.4	3.850	1.04	0.83 (0.69)	18.1	3.794	1.59

rp, Pearson correlation coefficient; R2, coefficient of determination; sd%, standard deviation of the percent differences; ms and bs, slope and intercept of linear regression.

Supplementary Table S3 provides a summary of literature results and recommendations for VASAT quantification from five MRI and CT studies in comparison with the present work. None of them, however, has explored the use of partial areas and volumes and proper measures of agreement (like R2) are only reported for one rather specific study.

Discussion

The current guidelines for assessing body fat distribution in obesity are mainly based on measurements of BMI and waist circumference. These quantitative measures continue to provide researchers and clinicians with some basic information but cannot be used to distinguish a healthy metabolic phenotype from one with an elevated cardiometabolic risk, in particular at an early stage.

Blüher et al., for example, have described metabolically healthy obesity (MHO) as a subentity in which excessive body fat accumulation does not lead to adverse metabolic effects[32]. Such subjects are characterized by an increased expandability of SAT that may be accompanied by a lower fat deposition in visceral or ectopic regions. A personalized therapy therefore calls for a reliable method to quantify the individual abdominal fat compartments.

Fat content and distribution are increasingly considered as predictive biomarkers in longitudinal studies of obesity. Monitoring abdominal fat contents during different types of weight-loss interventions (behavior, diet, drug, surgery) is also likely to contribute to the understanding of obesity despite the complexity of the disease. Image segmentation, however, is usually time consuming and may not be appropriate for larger cohorts or clinical settings. A rapid yet reliable quantification method is highly desirable and various works have already investigated the power of selected single-slice data for the prediction of whole-abdominal fat volume.

In both a clinical drug trial over 24 weeks[27] as well as a 6-month exercise intervention[29,33], however, single-slice MRI fat estimates were found to be inappropriate to assess the changes in intraabdominal fat volume. For such trials, especially those involving weight loss, multi-slice MRI measurements seem to be inevitable at the moment[28-31]. Maurovich-Horvat et al. attributed this shortcoming to the inter-individual heterogeneity of the abdominal fat distribution and a different axial dimension of the abdomen. A planimetric assessment of the SAT/VAT ratio was found to be substantially different from a volumetric one[30]. Another study by Kanaley et al. also showed a poor prediction of total abdominal fat volume and related changes during weight loss[31].

The present work provides a comprehensive quantitative evaluation of surrogate markers for the prediction of total ASAT volume in overweight and obese adults. Unlike previous studies in the literature, it specifically addresses the common challenge of image truncations which primarily occur in ventrolateral areas-given that patients are imaged in supine position. The measures of interest were therefore located dorsally except for abdominal wall fat thickness which was ventral but midsagittal.

The best correlation was observed for partial volume Vp-FH (Pearson's rp = 0.88; R2 = 0.78) followed by Vp-ASIS (rp = 0.83; R2 = 0.69) which is likely explained by the fact that volumetric measures encode more information about fat distributions than planimetric ones. Simpler measures like AWFASIS, HGFASIS, Ap-FH and Ap-SIAS with sd% values over 20% are regarded as less useful. It is worth noting, that simple lengths at the ASIS level—rp/R2 (AWF) = 0.73/0.54 and rp/R2 (HGF) = 0.69/0.47—agreed better than any partial area (FH or ASIS). These lengths can be determined without much effort making them favorable for routine application. There is still a considerable improvement when using the best volumetric measure Vp-FH (from rp = 0.73 to 0.88). Pausch et al., House et al., Li et al. and Raman et al. have already found ventral and dorsolateral lengths[27,34-36] to be convenient and reproducible but did not correlate or validate them against a reference.

Total ASAT volume is a likely candidate as a tissue marker for metabolic, epidemiological or pharmaceutical studies. The present study is among the first ones to quantify the reliability of various surrogate measures for VASAT. A look at five related studies (Table S3) shows that the focus until now has been on practical rather than methodological aspects. Most importantly, none of them have considered partial areas (or volumes) indicating that truncated datasets had been absent or excluded. A detailed comparison is therefore difficult.

The gender-specific conversion equations (12 in total) allow any researcher to compute ASAT volumes in a standardized manner, also retrospectively. The corresponding measures of agreement in Table 2 also indicate which combinations (gender and measure) tend to be more reliable (for example, seven cases with rp around 0.7 or higher). Equations ignoring the gender are provided for comparison. For practical considerations, parameter sd% appears to be more instructive than R2. All equations were derived with 10-mm thick slices, similar to other studies (Table S3).

Slice thickness will generally affect the tissue contours in the images (partial volume averaging) and consequently the exact equation. The overall effect for ASAT, however, should be moderate because the contours generally vary rather smoothly between transverse sections. Therefore, the provided equations should also hold approximately for slice thicknesses other than 10 mm which might be useful for other study designs.

This study is generally limited by its retrospective, single-center and mono-ethnic design. The results for men need to be interpreted with more care due to the lower number of participants. Manual segmentation is increasingly accepted as a reference standard despite the known variability between readers. One selection bias comes up because datasets were taken from existing subpopulations. An information bias towards lower degrees of obesity or overweight is inevitably introduced by requiring datasets with fully contained ASAT for validation.

In conclusion, partial subcutaneous fat volumes bound by the femoral heads seem to provide an acceptable workaround for whole-abdominal ASAT volumetry in subjects with truncated datasets frequently encountered for higher degree of obesity. The abdominal wall fat thickness at the ASIS level may be used with some reservation in cases where the amount of time for analysis is of importance.

Supplementary Information.

Supplementary Figure S1.

Supplementary Table S1.

Supplementary Table S2.

Supplementary Table S3.