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.

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/m 2 ) or severely obese (7.6%, BMI: > 40 kg/ m 2 ) 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 socioeconomic burden on society [3][4][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  www.nature.com/scientificreports/ 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][13][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][16][17][18] . A number of studies have argued that MRI should be preferred over CT in terms of safety [19][20][21][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/m 2 ). 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. 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 (A p-FH and A p-ASIS ) and partial volumes (V p-FH and V p-ASIS ). Nomenclature and definitions are given in Fig. 2 and Table 1.

Statistical analysis.
A least-squares regression (fit) between surrogate measure S and total ASAT volume was performed according to the equation with resulting fit parameters slope (m s ) and intercept (b s ). For a given surrogate measurement, total ASAT volume estimates, denoted with a tilde, for example, V~A SAT (d ven ), V~A SAT (A p-ASIS ) or V~A SAT (V p-FH ), can then be computed with these conversion parameters. The predictive value of a surrogate measure was quantified with www.nature.com/scientificreports/ the coefficient of determination R 2 and the standard deviation s d% of the percent differences between estimated and fully measured volume (V~A SAT -V ASAT )/V ASAT × 100%. All statistical analyses were performed with SPSS version 24 (IBM, Armonk, NY). 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 (V p-FH ) and simple abdominal wall fat at the level of the anterior superior iliac spine (AWF ASIS ). There is generally good to moderate agreement of the data with the linear fits. The coefficients of determination R 2 are considered as good for V p-FH and moderate for AWF ASIS . For both measures, the mere R 2 values were slightly higher for males. The resulting slope m s and intercept b s 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 V p-FH (vs. AWF ASIS ) and also smaller for males (vs. females) for both surrogates.

Results
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   www.nature.com/scientificreports/ all surrogate measures. At about 14%, V p-FH represents the best surrogate followed by V p-ASIS at around 17.5%. The s d% values of partial areas and simple distances were all above 20% with the next best agreement for AWF ASIS . Gender differences were not systematic and on the order of a few percents except for V p-FH . Table 2 provides a summary of all results and allows for a detailed inspection of both measures of agreement (R 2 and s d% ) and gender differences. V p-FH happens to be the surrogate with both smallest s d% and highest coefficient of determination as well as no substantial difference between females and males. By design, the reported m s and b s 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.
Supplementary Table S3 provides a summary of literature results and recommendations for V ASAT 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 R 2 ) 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, www.nature.com/scientificreports/ 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, singleslice 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][29][30][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 V p-FH (Pearson's r p = 0.88; R 2 = 0.78) followed by V p-ASIS (r p = 0.83; R 2 = 0.69) which is likely explained by the fact that volumetric measures encode more information about fat distributions than planimetric ones. Simpler measures like AWF ASIS , HGF ASIS , A p-FH and A p-SIAS with s d% values over 20% are regarded as less useful. It is worth noting, that simple lengths at the ASIS level-r p /R 2 (AWF) = 0.73/0.54 and r p /R 2 (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 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 V ASAT . A look at five related studies (Table S3) shows that the focus until now has been on practical rather than  www.nature.com/scientificreports/ 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 r p around 0.7 or higher). Equations ignoring the gender are provided for comparison. For practical considerations, parameter s d% appears to be more instructive than R 2 . 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.