Skip to main content
  • Research article
  • Open access
  • Published:

The relationship between DXA-based and anthropometric measures of visceral fat and morbidity in women

Abstract

Background

Excess accumulation of visceral fat is a prominent risk factor for cardiovascular and metabolic morbidity. While computed tomography (CT) is the gold standard to measure visceral adiposity, this is often not possible for large studies - thus valid, but less expensive and intrusive proxy measures of visceral fat are required such as dual-energy X-ray absorptiometry (DXA). Study aims were to a) identify a valid DXA-based measure of visceral adipose tissue (VAT), b) estimate VAT heritability and c) assess visceral fat association with morbidity in relation to body fat distribution.

Methods

A validation sample of 54 females measured for detailed body fat composition - assessed using CT, DXA and anthropometry – was used to evaluate previously published predictive models of CT-measured visceral fat. Based upon a validated model, we realised an out-of-sample estimate of abdominal VAT area for a study sample of 3457 female volunteer twins and estimated VAT area heritability using a classical twin study design. Regression and residuals analyses were used to assess the relationship between adiposity and morbidity.

Results

Published models applied to the validation sample explained >80% of the variance in CT-measured visceral fat. While CT visceral fat was best estimated using a linear regression for waist circumference, CT body cavity area and total abdominal fat (R 2 = 0.91), anthropometric measures alone predicted VAT almost equally well (CT body cavity area and waist circumference, R 2 = 0.86). Narrow sense VAT area heritability for the study sample was estimated to be 58% (95% CI: 51-66%) with a shared familial component of 24% (17-30%). VAT area is strongly associated with type 2 diabetes (T2D), hypertension (HT), subclinical atherosclerosis and liver function tests. In particular, VAT area is associated with T2D, HT and liver function (alanine transaminase) independent of DXA total abdominal fat and body mass index (BMI).

Conclusions

DXA and anthropometric measures can be utilised to derive estimates of visceral fat as a reliable alternative to CT. Visceral fat is heritable and appears to mediate the association between body adiposity and morbidity. This observation is consistent with hypotheses that suggest excess visceral adiposity is causally related to cardiovascular and metabolic disease.

Peer Review reports

Background

Body fat distribution, particularly abdominal adiposity, is strongly associated with a range of chronic metabolic and cardiovascular morbidities. The abdomen includes depots of both subcutaneous and visceral adipose tissue, which are distinct tissues that present different gene expression [1], endocrinological profiles [2] and pathogenicity [3]. The detrimental effect of excess fat accumulation is contingent on body fat distribution. Kim et al.[3] have shown that the diabetic phenotype of obese leptin knock out mice can be rescued by shifting fat storage away from the viscera and into the subcutaneous depot. The high ratio of subcutaneous to visceral fat in this mouse model is in contrast to the familial partial forms of lipodystrophy observed in humans, a pathology characterised by a redistribution of adipose tissue to the intra-abdominal region, giving rise to severe metabolic abnormalities as a consequence. These extreme cases illustrate the potent effect that adipose distribution can have upon health [4, 5]. The accumulation of visceral fat may also play a role in the aetiology of cardiovascular and metabolic disease.

Adiposity can be measured with a variety of techniques such as dual-energy X ray absorptiometry (DXA), computed tomography (CT) and magnetic resonance imaging (MRI). In whole body composition analysis, CT is regarded as the gold standard and as such has been used extensively as a benchmark for the quantitative assessment of subcutaneous and visceral adiposity. However, CT and MRI are limited in large research studies due to cost [6] and radiation exposure [7]. As such, there is a continued effort to identify informative adiposity measures that are safe, easily acquired and show robust correlation with CT-measured visceral fat content [811].

Family and twin studies can be used to estimate the heritability (the proportion of trait variance explained by genetic factors) and the relative importance of genetic and shared environmental factors in influencing phenotypic variance, by contrasting trait covariance between relatives with degree of relatedness. Establishing the heritable basis of a trait is also an important consideration prior to gene mapping studies [12].

The aims of this study are to 1) identify and estimate a valid DXA and anthropometric based measures of visceral fat that can be used in this and future epidemiological studies, 2) estimate the heritability and the degree of shared familial environment that influences visceral fat accumulation and 3) to assess whether association between body composition measures of adiposity and morbidities can be fully explained by visceral fat. We hypothesise that if visceral fat is causally related to morbidity, we would expect visceral fat to remain associated with morbidity conditional upon all other measures of body fat distribution.

Methods

Subjects and data collection

St. Thomas’ Research Ethics Committee approved the study (EC96/439 Twins UK) and all participants provided informed written consent. Participants were female volunteers enrolled in the TwinsUK resource, originally established in 1992 to study the heritability of osteoarthritis and osteoporosis in women, but subsequently healthy male and female twins have been recruited (mainly via national advertisements) without respect to any particular disease [13]. While the twins are largely representative of the age-matched UK female population [14], our findings relating to adiposity are restricted to Caucasian women over 40 years old with a BMI ≤ 35.

The validation sample comprised 54 female individuals over the age of 40 years with both CT and DXA data, but was otherwise unselected. Previously reported linear regression models [10, 11, 15] and adiposity indices [8] were assessed for this sample to identify the most predictive combination of DXA abdominal fat and anthropometry measures required to estimate visceral fat. A validated best-fit model was then used to realise an estimate of abdominal visceral fat for a study sample of 3,457 female individuals (533 monozygotic twin pairs, 1102 dizygotic pairs, 187 singletons), where date-matched DXA abdominal fat and anthropometric data were available.

CT

Subjects were placed in the supine position and a truncal CT scan was performed using a CT helical scanner (Brilliance CT component of the Philips Precedence SPECT/CT system, Philips, Eindhoven, Netherlands) with slice thickness set at 5 mm. The total cross sectional area (CSA) was intersected at the level of the intervertebral disc between lumbar vertebrae four and five (L4:L5). CT scan were analysed using Osirix X imaging software, version 3.7.1 (Pixmeo, Switzerland) to provide high-resolution 3-dimensional images. The body cavity CSA was obtained by tracing the outer contour of the abdominal wall and adipose tissue within this area was identified as having an attenuation value between −190 and −30 Hounsfield units [16]. For the 54 subjects with CT and DXA scan data (“validation sample” from here on in), visceral fat was directly measured as a pixel count using CT single slice area (VAT) (area). VAT area was defined as (visceral pixel count)/(total body cavity pixel count) × BC CSA. Single slice visceral fat area and volume are highly correlated [17, 18].

DXA

Baseline whole body DXA scans were performed for 3,457 female twins aged 40–80 years (“study sample” from here on in), in the supine position, using a fan beam X-ray bone densitometer (QDR-4500 W, Hologic, Massachusetts) and analysed with QDR System Software Version 12.6 (Hologic, Inc., Massachusetts). The DXA scanner was calibrated daily with a spine phantom and weekly using the step phantom as instructed by the manufacturer’s guidelines. The region of interest for body composition analysis was manually defined in a similar manner to Bertin et al.[8]. The abdominal region was delineated by an upper horizontal border located at half of the distance between acromions and external end of iliac crests and a lower border determined by the external end of iliac crests, but the lateral borders were delineated by the inner body wall rather than laterally past any trunk soft tissue. DXA fat mass from this abdominal region was recorded and abdominal transverse internal and external diameters were measured as defined by Bertin et al.[8].

Anthropometry

Anthropometrical measurements for both the validation and study samples were height, weight, BMI and waist circumference [14]. Waist circumference for the validation sample was based upon the transverse circumference of CT body scan image at the waist, while the study sample was a tape measurement taken at the same time as the DXA scan. CT images were used to measure for sagittal depth, body cavity, subcutaneous fat and total abdominal fat cross sectional areas, transverse internal and external diameters (TID and TED respectively) and subcutaneous fat width (SFW). Subcutaneous fat width (SFW) was estimated as SFW = (TED – TID)/2, which was also used as an estimate for skin fold as a calliper skin fold measure was not taken for the TwinsUK sample. The 54 subjects were also measured for TED and TID (with SFW derived) using DXA software to draw a horizontal line on the DXA image at the upper most point of the iliac crest, which was converted from pixel to mm length by multiplying by a constant of 2.048 (DXA software support, Vertec Scientific Limited).

Type 2 diabetes

Subjects were identified as type 2 diabetic (T2D) if they were classified as hyperglycaemic for one or more of the following diagnostics: fasting glucose serum concentration > 7 mmol/L, fasting 2 hour oral glucose tolerance test > 11.1 mmol/L or glycosylated haemoglobin > 48 mmol/mol. Questionnaire data on physician diagnosis and medication was also used to identify self-reported cases of T2D (2010–2011). Fasting plasma glucose was measured by enzymatic colorimetric slide assay (Johnson and Johnson Clinical Diagnostic Systems, Amersham UK).

Hypertension

Blood pressure was measured twice using an automatic blood pressure monitor (Omron Healthcare, Inc., Bannockburn, IL) using data from the same visit as the DXA scan. Individuals were classified as having hypertension (HT), if currently receiving anti-hypertensive medication and/or the repeated systolic/diastolic blood pressure was greater than 140/90 mm Hg.

Carotid intima-media thickness

Carotid intima-media thickness (cIMT) was measured as previously described by Cecelja et al.[19] and the quantitative trait was dichotomised at the 75th percentile as a surrogate for subclinical atherosclerosis [20, 21]. In brief, the left and right carotid and femoral arteries were visualized with B-mode ultrasound (Siemens CV70, Siemens, Erlangen, Germany, with 13-MHz vascular probe). Common carotid IMT was measured in the near and far walls, 1 to 2 cm proximal to the carotid bifurcation with automated wall-tracking software (Medical Imaging Applications, Coralville, Iowa) during diastole in an area free of overt plaque. Mean values of IMT in the near and far walls of both arteries were used for analysis.

Liver function tests

Data on liver function test (LFT) proteins were available for 87% of the study sample (for details, see Rahmioglu et al.[22]). Abnormal liver function as assessed by LFTs can be used to assist diagnosis of non-alcoholic fatty liver disease (NAFLD) [23]. Serum concentrations of alanine transaminase (ALT), alkaline phosphatase (ALK), total bilirubin (BIL) and γ-glutamyl transpeptidase (GGT) were used in this study as markers of liver function. Of the four tests, ALT is the most sensitive marker for liver cell damage (as a consequence of disease or drug use), while ALK and GGT are more indicative of cholestatic injury and BIL, haem catabolism [24], [25]. The upper limit of normal threshold used for each protein were: ALT: 39.4 IU/L [26], ALK: 81 IU/L (taken as the 75% percentile of the sample distribution), BIL: 17.1 μmol/L [27], GGT: 33 IU/L [28].

Statistical analyses

Validation sample

Multiple regression analyses [10, 11, 15] and indices of adiposity were assessed [8] to address two questions relating to the estimation of total abdominal visceral fat using DXA adiposity and a range of anthropometric measures: 1) which of these predictive models for the “gold standard” CT measure of visceral fat, best fit our validation sample of 54 females; and 2) in relation to previous discussions [8], can visceral fat be reliably estimated based upon anthropometry alone. CT and DXA scans for the same individuals were date-matched to between 0.23 - 2.3 years of one another. The difference in scan date for the validation sample was included in all visceral fat regression models as a nuisance factor. A Bland-Altman analysis was conducted to assess if the predicted VAT error term was constant or varied across the range of CT-measured VAT area.

Heritability analysis (study sample)

The classical twin model was used to estimate the relative influence of genetic and environmental factors upon individual variation about the visceral fat sample mean [29]. Using variance components analysis, the total phenotypic variance was partitioned into estimates of the additive polygenic (A), dominant (D) genetic and shared familial environment (common; C) and unique to the individual, environmental and measurement error (E) components using genetically informative data. The model assumes no epistasis, gene-environment correlation or interactions and that shared environmental factors are not confounded by zygosity (the equal environments assumption). Provided these assumptions hold, on the basis that MZ twins share identical genes and DZ twins share, on average, half their segregating genes, twin data can be used to infer heritability and shared familial estimates. Maximum likelihood models were implemented using Mx [30].

In addition to the univariate heritability analysis, we also conducted a bivariate variance components analysis between the realised estimate of VAT area and total abdominal fat as an indirect means to further validate the proxy measure of VAT. Bivariate analysis facilitates a test of genetic correlation between two variables, but also whether each variable has a specific genetic component that is not shared by the two traits. We used the specific test to assess if there was evidence for a genetic component that is unique to VAT area and not shared with DXA total abdominal fat [30].

Morbidity association with adiposity (study sample)

Type 2 diabetes, hypertension and liver function were associated with the different measures of adiposity variables and age using logistic regression using a cross-sectional study design. By contrast, prospective incident subclinical atherosclerosis (cIMT) was modelled using Cox regression and baseline measures of visceral fat taken on average 10 years (range 5–16 years) previously. Since there was no baseline examination of cIMT, we could not exclude the small proportion of individuals who may have already developed subclinical atherosclerosis at baseline. To the extent this was true, the study design for these individuals would have been also cross-sectional rather prospective.

LFTs were performed using blood samples collected at the same visit as the DXA scan. The mean, standard deviation and median were all examined to identify potential batch effects across years, but no obvious trend was identified. Therefore the year of visit was categorised as quintiles and included in all analyses as a categorical confounding variable.

Co-linearity between the explanatory measures of VAT area, DXA total abdominal fat, BMI and age was assessed using pairwise correlations and residuals analysis [31, 32]. A variance inflation factor (VIF) was calculated for these data, with a score of ≥ 10 sometimes used as a (conservative) threshold indicator for potentially problematic co-linearity between model explanatory variables [32]. To assess consistency of results, the multiple regression analyses were repeated using random subgroups of data and using quantitative traits for all subjects, where morbidity had been defined using a quantitative trait threshold (i.e. systolic and diastolic blood pressure and LFTs).

Parsimonious best-fit multiple regression models for morbidity were identified by a) using a likelihood ratio test (LRT) to assess the contribution of each explanatory variable to the full model; b) assessing model fit using pseudo-R 2 for T2D, HT and LFTs and the Wald statistic for incident cIMT, and c) the implementation of a linear residuals analyses [31] to assess whether adiposity measures are independent of morbidity, when conditioned upon on visceral fat. For example, to assess if VAT completely mediates the association between BMI and T2D, secondary residuals analysis involved taking the residuals for T2D on the logit scale and the ordinary least squares residuals for BMI, both with respect to VAT area and age. Ordinary least squares (OLS) residuals for DXA total abdominal fat were also taken with respect to VAT area and age. The secondary residuals analysis was implemented using linear regression of T2D residuals upon BMI and DXA residuals. No evidence of association between residuals would imply that BMI/DXA association with morbidity is primarily mediated via VAT area, while a significant residual association demonstrates either that BMI/DXA are associated with morbidity through a path not mediated by VAT area, or that the parametric assumptions of the association models are false.

To facilitate assessment of the relative importance of adiposity measures associated with morbidity, unadjusted and adjusted odd ratios (OR) results are presented for VAT area (cm2), DXA total abdominal fat (kg) and BMI (kg/m2). To facilitate comparison between the different adiposity units, all adiposities were standardised for the regression analyses. Robust standard errors were estimated by grouping twin pairs using the cluster option in Stata to account for intra-family relatedness. All statistical analyses were performed using Stata version 11.1 (StatCorp, Texas).

Results

Estimation of visceral fat (validation sample)

Descriptive statistics for anthropometric and body composition for the validation and study samples are presented in Table 1. All subjects were female and over 40 years old at examination. The mean height, weight, BMI, DXA total abdominal fat, blood pressure, cIMT, liver function tests and visceral fat do not statistically differ between the two samples, although the mean age for the validation sample was six years older (p = 6.5 × 10-11) with a larger waist circumference (p = 2.05 × 10-6) than the study sample. The study sample prevalence and 95% confidence interval for morbidity related to the quantitative traits presented in Table 1 were as follows: T2D 0.046 (0.04-0.05), HT 0.08 (0.07-0.09), cIMT 0.27 (0.24-0.30), ALT 0.22 (0.21-0.24), ALK 0.27 (0.26-0.29), BIL 0.02 (0.02-0.03) and GGT 0.20 (0.18-0.21).

Table 1 Validation and study sample characteristics

The Pearson product moment correlation coefficients (r) between the different adiposity measures and CT measured VAT area for the validation sample are presented in Table 2. CT-VAT area was most strongly correlated with CT-measured body cavity cross sectional area (r = 0.85), sagittal depth (r = 0.84) and tape-measured waist circumference (r = 0.86) and DXA total abdominal fat (r = 0.79). Consistent with these data, we observed that reported models of visceral fat in the literature, whether linear regressions or derived anthropometric indices, all attempt to capture information about the body cavity volume (or area) in relation to the subcutaneous volume [8, 10, 11, 15]. We used this insight to guide our choice of linear regression to estimate CT-measured visceral fat.

Table 2 Validation sample (n = 54) correlation coefficients between CT visceral adipose fat (VAT) area, anthropometric and abdominal fat measures

Although the DXA scans were collected between 0.23 – 2.3 years after the CT scans for the validation sample, no evidence was observed for significant change in weight in these individuals nor was change in weight correlated with time lapse between scan dates (data not shown).

Table 3 presents the results for three previously published DXA-based regression models and anthropometric indices for estimating visceral fat, applied to the TwinsUK CT validation sample. The best-replicated regression models included DXA trunk fat and sagittal depth (R 2 ≈ 0.8), while a combination of DXA and skin fold was less predictive of visceral fat. The best individual indices were functions of sagittal depth (SD), SFW, TID and TED (r ≥ 0.85, equivalent to r2 ≥ 0.72). We note that the most reproducible indices of visceral fat all relate to body cavity CSA. By assuming body cavity CSA takes the form of an ellipse, we estimated this as body cavity CSA = π x (SD – 2SFW) × TID for our validation sample.

Table 3 Visceral adipose tissue area (VAT area) linear model estimates and correlational indices

In modelling CT-measured VAT area, the best fit and most interpretable model included a combination of measures for DXA abdominal fat, body cavity cross sectional area (estimated using the ellipse formula above) and waist circumference (Table 4, Model 0), which together explained 91% of the variance in CT-measured VAT area (R 2 = 0.91). However, since sagittal depth was not available for our study sample for which we wished to estimate VAT, we also assessed a model including only DXA total abdominal fat, tape-measured waist circumference and age. For this (Table 4, Model 1) we obtained a model with an R 2 of 0.83 with the following linear regression equation: VAT area = 10.1(DXA abdominal fat mass) + 40.8(waist circumference) + 1.4(age) using standardised explanatory variables and no intercept term. For this estimate, a Bland-Altman analysis showed no evidence of heteroscedascity across the full range of CT-measured VAT area, with only 2 out of 54 (3.7%) values with a difference outside the 95 limits of agreement (the mean difference ± twice the standard deviation of the difference between the two measures).

Table 4 Linear regression models for CT visceral adipose fat (VAT) area using the validation sample (n = 54)

In relation to efforts attempting to estimate visceral fat using only anthropometric measures [8, 33], we also obtained a highly explanatory model (R 2 = 0.86) with a linear equation using only two CT measures of body cavity CSA (estimated as an ellipse) and waist circumference (Table 4, Model 2). This figure rose to R 2 = 0.89 using body cavity components SD, SFW and WC as explanatory variables for CT VAT area (data not shown). Body cavity CSA results are presented, as this model is more interpretable, while the model including SD, SFW and WC yields a negative beta coefficient for SFW due to co-linearity. Again, these simple anthropometric models could not be used for the study sample however, as sagittal depth was not recorded along with DXA scan for these subjects.

In addition to these validation models, we indirectly assessed the validity of our study sample estimates of visceral fat by making two observations:

  1. 1.

    Realised estimates of VAT area (VATModel1 and VATModel2) for the validation sample were equally or more strongly correlated with VATCT (r = 0.89 and r = 0.93) than DXA total abdominal fat (r = 0.88 and r = 0.70, respectively);

  2. 2.

    Bivariate variance component analysis between VATModel1 and DXA total abdominal fat provided strong statistical evidence (Δχ2 1 = 43.7, p = 4x10-11) for a specific heritable component that was unique to VATModel1 and not shared with DXA total abdominal fat (data not shown).

Heritability (study sample)

The data provided a narrow sense heritability estimate of 0.58 (95% CI: 0.51-0.66) and a familial environmental effect of 0.24 (95% CI: 0.17-0.30) for VAT area (Table 5), with the best fit to the data being the ACE model, including additive polygenic (A), common familial (C) and unique environmental (E) components.

Table 5 Visceral adipose tissue (VAT) area estimate of heritability ( h 2 ) and model fit statistics (n = 3457)

Visceral fat as a risk factor of morbidity (study sample)

For the study sample we estimated visceral fat as a linear function of DXA total abdominal fat, waist circumference and age (Model 1, Table 4). Although a model including DXA total abdominal fat, waist circumference and estimated body cavity area was a better fit (Model 0, Table 4), we used Model 1 as we did not have a measure of sagittal depth for the study sample subjects required to estimate body cavity area from the DXA images.

The correlations between study sample explanatory variables are presented in Additional file 1: Table S1. As a function of these variables, VAT area is most strongly correlated with waist circumference and DXA total abdominal fat. While the univariate odd ratios for VAT area, DXA abdominal fat and BMI were all significantly associated with each morbidity (Tables 6, 7, 8 and 9), VAT area was most consistently and strongly associated with four morbidity traits - T2D, HT, cIMT and ALT. Residual analyses for these four morbidities were consistent with visceral fat entirely mediating the association with other measures of adiposity (BMI and total abdominal fat), but this was not true for the liver function tests ALK, BIL and GGT (details provided in supplementary Tables online).

Table 6 Type 2 diabetes (T2D) and adiposity
Table 7 Hypertension and adiposity
Table 8 Sub-clinical atherosclerosis and adiposity
Table 9 Liver function tests (LFTs) and adiposity

For type 2 diabetes, while the univariate ORs for the 3 adiposity measures were all associated with T2D (Table 6A), visceral fat and age provided the best-fit multiple regression model (pseudo-R 2 = 0.08, Table 6B and LRT, Additional file 1: Table S2), with an adjusted OR of 2.08 (95% CI 1.76 – 2.47) per standard deviation increment in VAT area including age. Removing VAT area from the full model resulted in a significant decline in model fit (LRT χ2 1 = 14.4, p = 1E-04), while the removal of DXA total abdominal fat and BMI, either individually or together (χ2 2 = 1.9, p = 0.39) did not reduce the model fit (Additional file 1: Table S2).

Hypertension was equally strongly associated with VAT area, DXA abdominal fat and BMI for univariate analyses (Table 7A), but visceral fat and age provided the best-fit multiple regression model (pseudo-R 2 = 0.10, Table 7B and LRT, Additional file 1: Table S4), with an adjusted OR of 1.90 (95% CI 1.60 – 2.25) per standard deviation increment in VAT area including age. Removing VAT area from the full model resulted in a nominal decline in model fit (LRT χ2 1 = 7.1, p = 0.01), whilst the removal of DXA total abdominal fat and BMI, either individually or together (χ2 2 = 1.99, p = 0.37) did not reduce the model fit (Additional file 1: Table S4).

The prospective analysis of subclinical atherosclerosis had a median follow up time of 9.7 years, during which a total of 221 (27%) individuals were classified as sub-clinically atherosclerotic. Univariate Cox proportional hazard models showed all three measures of adiposity and age at baseline to be associated with incident cIMT (Table 8A), with VAT area (χ2 1 = 43.8) and age (χ2 1 = 40) providing the best-fit parsimonious model (Table 8B). DXA and BMI could be dropped from the full model with no nominal (χ2 2 = 5.7, p = 0.06) deterioration in model fit (Additional file 1: Table S6).

All LFT protein serum levels were positively associated with measures of adiposity (Table 9), except for bilirubin, which was negatively associated. VAT area remained associated with alanine transaminase (ALT) when conditioned upon DXA total abdominal fat and BMI, while VAT area and DXA were still associated with ALK, BIL and GGT conditional upon BMI (Table 9 and related Additional file 1: Tables S8-S15). Removing VAT area from the full model for ALT resulted in a significant decline in model fit (LRT χ2 1 = 19.6, p = 1E-05), while the removal of DXA total abdominal fat and BMI, either individually or together (χ2 2 = 1.56, p = 0.46) did not (Additional file 1: Table S8).

The estimated variance inflation factor (VIF) between VAT area and DXA abdominal fat, BMI and age was 8.84. Analyses repeated using subgroups of data and analyses using underlying quantitative traits for HT, ALT, ALK, BIL and GGT all provided qualitatively the same association results (data not shown).

Discussion

This study shows that DXA and anthropometric measures can be used to derive reliable estimates of visceral fat and provides evidence that visceral fat - independent of BMI and total abdominal body fat - is a risk factor for type 2 diabetes, hypertension, subclinical atherosclerosis and ALT liver function. Validation sample explanatory models were able to explain >80% of the variance in CT-measured visceral fat. A combination of DXA and anthropometric measures (R 2 = 0.91) or models including only anthropometric measures (R 2 = 0.86 – 0.89), both provided equally good estimates of visceral fat. We obtain a heritability of 58% for visceral fat, which is consistent with familial studies that report within the region of 48-57% [3437] and using a twin design, confirms for the first time that the observed familial component is also partly due to shared familial environment (24%).

Our results confirm that visceral fat is the single most important measure of adiposity for risk of type 2 diabetes. In large DXA resources, such as the one used by Leslie et al. (2010) [38], it would be interesting to ascertain whether the risk of developing diabetes is better predicted if a DXA-based estimate of visceral fat was used rather than using DXA total abdominal fat.

Visceral fat mediates all associations between hypertension and adiposity variables for these cross sectional data. Similarly, in a longitudinal study, Hayashi et al. (2004) [39] showed that only visceral fat is significantly associated with hypertension, when visceral fat is adjusted for other adiposity measures (e.g. BMI, waist circumference or abdominal subcutaneous fat). Visceral and subcutaneous adiposity both contribute to the prevalence of hypertension, but when adjusted for BMI or waist circumference, subcutaneous fat is no longer associated with hypertension [37]. This further illustrates that it is the distribution between these two fat depots that appears to determine the risk of morbidity.

The analysis of incident subclinical atherosclerosis is consistent with cross-sectional studies that show association between atherosclerosis, total abdominal fat [40] and visceral fat [41, 42]. Here we provide additional evidence that total abdominal and visceral fat accumulation precedes atherosclerosis.

The LFT protein results suggest that liver function deteriorates with increasing visceral adiposity with VAT area positively associated with ALT, ALK and GGT. By contrast, the inverse relationship between visceral fat and bilirubin is consistent with the inverse association between bilirubin, insulin resistance [43] and GGT activity [44]. In particular, of the four tests, the linear association between ALT and adiposity is specifically mediated via visceral fat, while for ALK, BIL and GGT there is only evidence to suggest this is the case for abdominal fat per se. The visceral adipose depot is thought to be the major source of elevated fatty acid delivered to the liver via the portal vein and it is possible that visceral fat acts as a marker of hepatic fat content [45]. Hepatic fat accumulation is known to impair insulin signalling in hepatocytes [46, 47] and “pathway selective insulin resistance” [48] maybe pivotal in the transition from normal to impaired fasting glucose states. It has been hypothesised that non-alcoholic fatty liver (NALFD) may play a mediatory role in the pathology of CVD [49].

The relative strengths and limitations of this study were: 1) Adiposity inferences are restricted to middle-aged Caucasian women with a BMI ≤ 35. The VAT area estimates for the study sample included the 4% of females with a BMI > 35, which was outside the validation sample data range. Although these estimates strictly represent an invalid data extrapolation, their exclusion did not qualitatively alter the results for the heritability or morbidity association (data not shown) and so these individuals were retained. 2) We have shown that sagittal depth is an extremely useful yet simple anthropometric measurement to estimate the body cavity area. When used in combination with the subcutaneous fat width, this measure contributes to a very good estimation of visceral fat (R 2 approximately 90%). Where DXA is not available, sagittal depth, waist circumference and subcutaneous fat width could therefore be measured to quickly and inexpensively estimate a person’s visceral fat content with just a tape measure and skin callipers. Not having sagittal depth for the study sample limited our ability to estimate visceral fat more accurately. This is indicated by the best anthropometric validation sample estimate (Model 2) predicting VAT area as well or better than a combination of DXA total abdominal fat, waist circumference and age (Model 1); 3) Arising from the second limitation, although our prediction of visceral fat was still valid [32], co-linearity between model explanatory variables – in particular VAT area and DXA total abdominal fat – meant that adjusted Odds Ratios for risk factors were not interpretable. Instead, conditional independence structures were assessed using residuals analysis. This allowed us to eliminate the possibility of causal relationships between adiposity variables and morbidity, except those mediated by visceral fat; 4) The twin-based heritability estimate presented here is likely to be an under-estimate, since random measurement error is included in the denominator of the heritability ratio estimate and the study sample estimate of visceral fat is known to account for approximately only 83% of CT-measured VAT area; 5) The study design for incident subclinical atherosclerosis was not strictly prospective, since no baseline measure of cIMT was taken to facilitate the exclusion of potential baseline cases. Where undetected baseline cases exist, the study design is cross sectional; 6) Inclusion criteria for the validation sample were that study subjects were female, over 40 and had both DXA and CT scans within 2 years of one another. Although not measured the same day, time between visits was included in analyses as a confounding variable and Bland-Altman plots revealed no bias in the relationship between predicted and CT-measured visceral fat values. More importantly, no significant weight gain was observed between the collection of CT and DXA measurements and small fluctuations in weight were not correlated with the time lapse between scan dates.

This study shows that, where a direct measure of visceral fat is unavailable, using an indirect estimate of visceral fat (even without a measure of subcutaneous fat) is more predictive of morbidity than total abdominal fat. Our study suggests that both DXA body fat composition and anthropometric measures of body cavity volume alone can provide reliable estimates of abdominal visceral fat.

In conclusion, we demonstrate that DXA-based measures of visceral fat can be both reliable and valid and are a cheaper and less intrusive alternative to other methods such as MRI and CT to estimate adiposity. Visceral fat is not only strongly associated with morbidity, but also appears to mediate the association between morbidity and other measures of adiposity for type 2 diabetes, hypertension and elevated alanine transaminase serum levels. This is consistent with hypotheses that suggest excess visceral fat accumulation is causally related to cardiovascular and metabolic disease.

Abbreviations

ALK:

Alkaline phosphatase

ALT:

Alanine transaminase

BIL:

Total bilirubin

cIMT:

Carotid intima-media thickness (subclinical atherosclerosis)

CT:

Computed tomography

CSA:

Cross sectional area

DXA:

Dual-energy X-ray absorptiometry

GGT:

Gamma-glutamyl transpeptidase

HT:

Hypertension

SFW:

Subcutaneous fat width

TID:

Transverse internal diameter

TED:

Transverse external diameter

T2D:

Type 2 diabetes

VAT:

Visceral adipose tissue.

References

  1. Palou A, Palou M, Priego T, Sanchez J, Rodriguez AM, Pico C: Gene Expression Patterns in Visceral and Subcutaneous Adipose Depots in Rats are Linked to Their Morphologic Features. Cell Physiol Biochem. 2009, 24 (5–6): 547-556.

    Article  CAS  PubMed  Google Scholar 

  2. Bastard JP, Antuna-Puente B, Feve B, Fellahi S: Adipokines: The missing link between insulin resistance and obesity. Diabetes Metab. 2008, 34 (1): 2-11. 10.1016/j.diabet.2007.09.004.

    Article  PubMed  Google Scholar 

  3. Kim JY, De Wall EV, Laplante M, Azzara A, Trujillo ME, Hofmann SM, Schraw T, Durand JL, Li H, Li G: Obesity-associated improvements in metabolic profile through expansion of adipose tissue. J Clin Invest. 2007, 117 (9): 2621-2637. 10.1172/JCI31021.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Mantzoros CS, Fiorenza CG, Chou SH: Lipodystrophy: pathophysiology and advances in treatment. Nat Rev Endocrinol. 2011, 7 (3): 137-150. 10.1038/nrendo.2010.199.

    Article  PubMed  Google Scholar 

  5. Savage DB, Huang-Doran I, Sleigh A, Rochford JJ, O’Rahilly S: Lipodystrophy: metabolic insights from a rare disorder. J Endocrinol. 2010, 207 (3): 245-255. 10.1677/JOE-10-0272.

    Article  PubMed  Google Scholar 

  6. Silver HJ, Welch EB, Avison MJ, Niswender KD: Imaging body composition in obesity and weight loss: challenges and opportunities. Diabetes Metab Syndr Obes. 2010, 3: 337-347.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Brenner DJ, Hall EJ: Computed tomography–an increasing source of radiation exposure. N Engl J Med. 2007, 357 (22): 2277-2284. 10.1056/NEJMra072149.

    Article  CAS  PubMed  Google Scholar 

  8. Bertin E, Marcus C, Ruiz JC, Eschard JP, Leutenegger M: Measurement of visceral adipose tissue by DXA combined with anthropometry in obese humans. Int J Obesity. 2000, 24 (3): 263-270. 10.1038/sj.ijo.0801121.

    Article  CAS  Google Scholar 

  9. De Lucia RE, Norris SA, Sleigh A, Brage S, Dunger DB, Stolk RP, Ong KK: Validation of ultrasound estimates of visceral fat in black South African adolescents. Obesity (Silver Spring). 2011, 19 (9): 1892-1897. 10.1038/oby.2011.213.

    Article  Google Scholar 

  10. Hill AM, LaForgia J, Coates AM, Buckley JD, Howe PRC: Estimating abdominal adipose tissue with DXA and anthropometry. Obesity. 2007, 15 (2): 504-510. 10.1038/oby.2007.629.

    Article  PubMed  Google Scholar 

  11. Snijder MB, Visser M, Dekker JM, Seidell JC, Fuerst T, Tylavsky F, Cauley J, Lang T, Nevitt M, Harris TB: The prediction of visceral fat by dual-energy X-ray absorptiometry in the elderly: a comparison with computed tomography and anthropometry. Int J Obesity. 2002, 26 (7): 984-993.

    Article  CAS  Google Scholar 

  12. Andrew T, Aviv A, Falchi M, Surdulescu GL, Gardner JP, Lu X, Kimura M, Kato BS, Valdes AM, Spector TD: Mapping genetic loci that determine leukocyte telomere length in a large sample of unselected female sibling pairs. Am J Hum Genet. 2006, 78 (3): 480-486. 10.1086/500052.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Moayyeri A, Hammond CJ, Valdes AM, Spector TD: Cohort Profile: TwinsUK and Healthy Ageing Twin Study. Int J Epidemiol. 2012, 42: 76-85.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Andrew T, Hart DJ, Snieder H, De Lange M, Spector TD, MacGregor AJ: Are twins and singletons comparable? A study of disease-related and lifestyle characteristics in adult women. Twin Res. 2001, 4 (6): 464-477.

    Article  CAS  PubMed  Google Scholar 

  15. Treuth MS, Hunter GR, Kekesszabo T: Estimating Intraabdominal Adipose-Tissue in Women by Dual-Energy X-Ray Absorptiometry. Am J Clin Nutr. 1995, 62 (3): 527-532.

    CAS  PubMed  Google Scholar 

  16. Yoshizumi T, Nakamura T, Yamane M, Islam AHMW, Menju M, Yamasaki K, Arai T, Kotani K, Funahashi T, Yamashita S: Abdominal fat: Standardized technique for measurement at CT. Radiology. 1999, 211 (1): 283-286.

    Article  CAS  PubMed  Google Scholar 

  17. Shen W, Punyanitya M, Wang ZM, Gallagher D, St-Onge MP, Albu J, Heymsfield SB, Heshka S: Visceral adipose tissue: relations between single-slice areas and total volume. Am J Clin Nutr. 2004, 80 (2): 271-278.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. Irlbeck T, Massaro JM, Bamberg F, O’Donnell CJ, Hoffmann U, Fox CS: Association between single-slice measurements of visceral and abdominal subcutaneous adipose tissue with volumetric measurements: the Framingham Heart Study. Int J Obesity. 2010, 34 (4): 781-787. 10.1038/ijo.2009.279.

    Article  CAS  Google Scholar 

  19. Cecelja M, Hussain M, Greil G, Spector TD, Chowienczyk P: Arterial stiffening relates to arterial calcification but not to non-calcified atheroma as determined by multimodality imaging of the aorta. J Human Hypertens. 2010, 24 (10): 693-693.

    Google Scholar 

  20. Stein JH, Korcarz CE, Hurst RT, Lonn E, Kendall CB, Mohler ER, Najjar SS, Rembold CM, Post WS: Use of carotid ultrasound to identify subclinical vascular disease and evaluate cardiovascular disease risk: A consensus statement from the American society of echocardiography carotid intima-media thickness task force endorsed by the society for vascular medicine. J Am Soc Echocardio. 2008, 21 (2): 93-111. 10.1016/j.echo.2007.11.011.

    Article  Google Scholar 

  21. Benzaquen BS, Nguyen-Thanh HT: Screening for Subclinical Coronary Artery Disease Measuring Carotid Intima Media Thickness. Am J Cardio. 2009, 104 (10): 1383-1388. 10.1016/j.amjcard.2009.07.005.

    Article  Google Scholar 

  22. Rahmioglu N, Andrew T, Cherkas L, Surdulescu G, Swaminathan R, Spector T, Ahmadi KR: Epidemiology and Genetic Epidemiology of the Liver Function Test Proteins. PLoS One. 2009, 4 (2): e4435-10.1371/journal.pone.0004435.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Armstrong MJ, Houlihan DD, Bentham L, Shaw JC, Cramb R, Olliff S, Gill PS, Neuberger JM, Lilford RJ, Newsome PN: Presence and severity of non-alcoholic fatty liver disease in a large prospective primary care cohort. J Hepatol. 2012, 56 (1): 234-240. 10.1016/j.jhep.2011.03.020.

    Article  PubMed  Google Scholar 

  24. Aragon G, Younossi ZM: When and how to evaluate mildly elevated liver enzymes in apparently healthy patients. Cleve Clin J Med. 2010, 77 (3): 195-204. 10.3949/ccjm.77a.09064.

    Article  PubMed  Google Scholar 

  25. Burke MD: Liver function: test selection and interpretation of results. Clin Lab Med. 2002, 22 (2): 377-390. 10.1016/S0272-2712(01)00002-6.

    Article  PubMed  Google Scholar 

  26. Kariv R, Leshno M, Beth-Or A, Strul H, Blendis L, Kokia E, Noff D, Zelber-Sagie S, Sheinberg B, Oren R: Re-evaluation of serum alanine aminotransferase upper normal limit and its modulating factors in a large-scale population study. Liver Int. 2006, 26 (4): 445-450. 10.1111/j.1478-3231.2006.01197.x.

    Article  CAS  PubMed  Google Scholar 

  27. Lazo M, Selvin E, Clark JM: Brief communication: Clinical implications of short-term variability in liver function test results. Ann Int Med. 2008, 148 (5): 348-W376. 10.7326/0003-4819-148-5-200803040-00005.

    Article  PubMed  Google Scholar 

  28. Sabanayagam C, Shankar A, Li JL, Pollard C, Ducatman A: Serum gamma-glutamyl transferase level and diabetes mellitus among US adults. Eur J Epidemiol. 2009, 24 (7): 369-373. 10.1007/s10654-009-9346-7.

    Article  CAS  PubMed  Google Scholar 

  29. Martin N, Boomsma D, Machin G: A twin-pronged attack on complex traits. Nat Genet. 1997, 17 (4): 387-392. 10.1038/ng1297-387.

    Article  CAS  PubMed  Google Scholar 

  30. Neale MC, Cardon LR: North Atlantic Treaty Organization. Scientific Affairs Division.: Methodology for genetic studies of twins and families. Dordrecht. 1992, Boston: Kluwer Academic Publishers

    Book  Google Scholar 

  31. David Clayton MH: Statistical Models in. 1993, Epidemiology: Oxford University Press

    Google Scholar 

  32. O’Brien RM: A caution regarding rules of thumb for variance inflation factors. Quality & Quantity. 2007, 41 (5): 673-690. 10.1007/s11135-006-9018-6.

    Article  Google Scholar 

  33. Karelis AD, Rabasa-Lhoret R, Pompilus R, Messier V, Strychar I, Brochu M, Aubertin-Leheudre M: Relationship between the Bertin index to estimate visceral adipose tissue from dual-energy X-ray absorptiometry and cardiometabolic risk factors before and after weight loss. Obesity (Silver Spring). 2012, 20 (4): 886-890. 10.1038/oby.2011.273.

    Article  CAS  Google Scholar 

  34. Perusse L, Despres JP, Lemieux S, Rice T, Rao DC, Bouchard C: Familial aggregation of abdominal visceral fat level: Results from the Quebec family study. Metab-Clin Exp. 1996, 45 (3): 378-382. 10.1016/S0026-0495(96)90294-2.

    Article  CAS  PubMed  Google Scholar 

  35. Rice T, Despres JP, Daw EW, Gagnon J, Borecki IB, Perusse L, Leon AS, Skinner JS, Wilmore JH, Rao DC: Familial resemblance for abdominal visceral fat: the HERITAGE family study. Int J Obesity. 1997, 21 (11): 1024-1031. 10.1038/sj.ijo.0800511.

    Article  CAS  Google Scholar 

  36. Hong YL, Rice T, Gagnon J, Despres JP, Nadeau A, Perusse L, Bouchard C, Leon AS, Skinner JS, Wilmore JH: Familial clustering of insulin and abdominal visceral fat: The HERITAGE family study. J Clin Endocrinol Metab. 1998, 83 (12): 4239-4245. 10.1210/jc.83.12.4239.

    CAS  PubMed  Google Scholar 

  37. Fox CS, Massaro JM, Hoffmann U, Pou KM, Maurovich-Horvat P, Liu CY, Vasan RS, Murabito JM, Meigs JB, Cupples LA: Abdominal visceral and subcutaneous adipose tissue compartments - Association with metabolic risk factors in the Framingham Heart Study. Circulation. 2007, 116 (1): 39-48. 10.1161/CIRCULATIONAHA.106.675355.

    Article  PubMed  Google Scholar 

  38. Leslie RD: Predicting Adult-Onset Automunune Diabetes Clarity From Complexity. Diabetes. 2010, 59 (2): 330-331. 10.2337/db09-1620.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Hayashi T, Boyko EJ, Leonetti DL, McNeely MJ, Newell-Morris L, Kahn SE, Fujimoto WY: Visceral adiposity is an independent predictor of incident hypertension in Japanese Americans. Ann Int Med. 2004, 140 (12): 992-1000. 10.7326/0003-4819-140-12-200406150-00008.

    Article  PubMed  Google Scholar 

  40. Kramer CK, Von Muhlen D, Gross JL, Barrett-Connor E: A Prospective Study of Abdominal Obesity and Coronary Artery Calcium Progression in Older Adults. J Clin Endocrinol Metab. 2009, 94 (12): 5039-5044. 10.1210/jc.2009-1497.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Lear SA, Humphries KH, Kohli S, Frohlich JJ, Birmingham CL, Mancini GBJ: Visceral adipose tissue, a potential risk factor for carotid atherosclerosis - Results of the multicultural community health assessment trial (M-CHAT). Stroke. 2007, 38 (9): 2422-2429. 10.1161/STROKEAHA.107.484113.

    Article  PubMed  Google Scholar 

  42. Kim SK, Park SW, Kim SH, Cha BS, Lee HC, Cho YW: Visceral fat amount is associated with carotid atherosclerosis even in type 2 diabetic men with a normal waist circumference. Int J Obesity. 2009, 33 (1): 131-135. 10.1038/ijo.2008.222.

    Article  CAS  Google Scholar 

  43. Lin LY, Kuo HK, Hwang JJ, Lai LP, Chiang FT, Tseng CD, Lin JL: Serum bilirubin is inversely associated with insulin resistance and metabolic syndrome among children and adolescents. Atherosclerosis. 2009, 203 (2): 563-568. 10.1016/j.atherosclerosis.2008.07.021.

    Article  CAS  PubMed  Google Scholar 

  44. Giral P, Ratziu V, Couvert P, Carrie A, Kontush A, Girerd X, Chapman MJ: Plasma bilirubin and gamma-glutamyl transferase activity are inversely related in dyslipidemic patients with metabolic syndrome: Relevance to oxidative stress. Atherosclerosis. 2010, 210 (2): 607-613. 10.1016/j.atherosclerosis.2009.12.026.

    Article  CAS  PubMed  Google Scholar 

  45. Fabbrini E, Magkos F, Mohammed BS, Pietka T, Abumrad NA, Patterson BW, Okunade A, Klein S: Intrahepatic fat, not visceral fat, is linked with metabolic complications of obesity. Proc Nat Acad Sci USA. 2009, 106 (36): 15430-15435. 10.1073/pnas.0904944106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Samuel VT, Petersen KF, Shulman GI: Lipid-induced insulin resistance: unravelling the mechanism. Lancet. 2010, 375 (9733): 2267-2277. 10.1016/S0140-6736(10)60408-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Samuel VT, Liu ZX, Qu XQ, Elder BD, Bilz S, Befroy D, Romanelli AJ, Shulman GI: Mechanism of hepatic insulin resistance in non-alcoholic fatty liver disease. J Biol Chem. 2004, 279 (31): 32345-32353. 10.1074/jbc.M313478200.

    Article  CAS  PubMed  Google Scholar 

  48. Brown MS, Goldstein JL: Selective versus total insulin resistance: A pathogenic paradox. Cell Metab. 2008, 7 (2): 95-96. 10.1016/j.cmet.2007.12.009.

    Article  CAS  PubMed  Google Scholar 

  49. Targher G, Day CP, Bonora E: Risk of Cardiovascular Disease in Patients with Nonalcoholic Fatty Liver Disease. New Eng J Med. 2010, 363 (14): 1341-1350. 10.1056/NEJMra0912063.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

We would like to thank all the twin volunteers who participated who made this study possible and the staff at the Department of Twin Research who collected and processed the data for this study.

Funding

This work was supported by the Medical Research Council UK [Investigator Award 91993 to TA], the British Heart Foundation Project Grant [PG/06/032 to PC]. The TwinsUK resource is supported by the Wellcome Trust; European Community’s Seventh Framework Programme (FP7/2007-2013); and the Department of Health via the National Institute for Health Research (NIHR) comprehensive Biomedical Research Centre award to Guy’s & St Thomas’ NHS Foundation Trust in partnership with King’s College London.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Toby Andrew.

Additional information

Competing interests

The author declare that they have no competing interest.

Authors’ contributions

KD - design of the experiment, collection of data, analysis of data, writing of the manuscript; MC – collection of data; TDS, PC – collection of data, comments on the manuscript; WA, MF – analysis of data, comments on the manuscript; TA: design of the experiment, collection of data, analysis of data, writing of the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12872_2012_569_MOESM1_ESM.xls

Additional file 1: Table S1: Study sample (n = 3457) correlation and colinearity between measures of adiposity. Table S2. Type 2 diabetes (T2D) multiple regression analyses. Table S3. Residuals analysis for type 2 diabetes (T2D). Table S4. Hypertension (HT) multiple regression analyses. Table S5. Residuals analysis for hypertension (HT). Table S6. Carotid intima-media thickness (cIMT) Cox multiple regression analyses. Table S7. Residuals analysis for carotid intima-media thickness (cIMT). Table S8. Multiple regression analyses for alanine transaminase (ALT). Table S9. Residuals analysis for the quantitative trait (QT)/assay for alanine transaminase (ALT). Table S10. Multiple regression analyses for alkaline phosphatase (ALK). Table S11. Residuals analysis for the quantitative trait (QT)/assay alkaline phosphate (ALK). Table S12. Multiple regression analyses for bilirubin (BIL). Table S13. Residuals analysis for the quantitative trait (QT)/assay for bilirubin (BIL). Table S14. Multiple regression analyses for gamma-glutamyl transpeptidase (GGT). Table S15. Residuals analysis for the quantitative trait (QT)/assay gamma-glutamyl transpeptidase (GGT). (XLS 490 KB)

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Direk, K., Cecelja, M., Astle, W. et al. The relationship between DXA-based and anthropometric measures of visceral fat and morbidity in women. BMC Cardiovasc Disord 13, 25 (2013). https://doi.org/10.1186/1471-2261-13-25

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2261-13-25

Keywords