Study population
The present study had two phases. The first phase was aimed to adapt VAI to DAI and evaluate its association with cardiometabolic diseases. For that purpose, we selected healthy subjects without cardiovascular risk factors from the Genetics of Atherosclerotic Disease study (GEA Spanish acronym). The GEA study was designed at the National Institute of Cardiology Ignacio Chávez, to examine the genomic basis of coronary heart disease (CHD) and to evaluate its relation with traditional and emerging cardiovascular risk factors in an adult Mexican population. The sample included 1200 patients with established premature CHD, and a control group of 1600 subjects aged 35 to 70 years, with no clinical or family history of coronary arterial disease, selected from the blood donors attending blood bank of the National Institute of Cardiology, or recruited by advertisements posted in social service centers from June 2008 through November 2012.[27] Participants of the GEA study were extensively characterized, including imaging studies, biochemical measurements, anthropometry, medical history, and sociodemographic and nutritional information collected through validated questionnaires.[28] With the purpose to estimates the constant values needed for the DAI equation, a subgroup of 340 healthy subjects without diabetes and with glucose < 5.6 mmol/L, systolic/diastolic blood pressure < 140/90 mmHg, TG < 1.69 mmol/L, HDL-C ≥ 1.03/1.29 mmol/L for men/women, and BMI < 30 kg/m2, was selected from the control GEA participants. Additionally, 1418 subjects from the control group of the GEA study were used to estimate the prevalence of cardiometabolic abnormalities and its association with DAI, excluding those participants with TG > 3.15 mmol/L or BMI > 40 kg/m2 that are potential confounders of the original visceral adiposity index (Fig. 1a). [24] In the second phase of the present study we analyzed subcutaneous adipose tissue biopsies and evaluated body composition and circulating adipokines and cytokines in a sample of healthy subjects, to validate the usefulness of DAI as a marker of the morpho-functional state of the adipose tissue. Briefly, 350 healthy individuals without cardiovascular disease, dyslipidemia (TG > 3.15 mmol/L), infection disease, diabetes, cancer, or any autoimmune disease, were invited through phone calls and personal interviews (from February 2018 to May 2019), to participate in the study to investigate the morpho-functionality of their adipose tissue. After invitations, 52 healthy subjects voluntarily accepted to participate in the study and were enrolled, but only 36 participants were included in the analyses due to incomplete subcutaneous adipose tissue biopsy, biochemical measurements, anthropometry, and clinical information from some participants (Fig. 1b). Both studies were approved by the research and bioethics committee of the National Institute of Cardiology and were performed following the guidelines of the Helsinki Declaration. All participants voluntarily signed the institutional informed consent form.
Clinical and anthropometrical evaluation
Participants in both studies were interviewed by a trained research staff and completed questionnaires detailing medical history, demographic characteristics, CHD history, medication, alcohol, and tobacco use. Height, weight, and waist circumference were measured; BMI was estimated based on weight (kg) divided by height (m2). Waist circumference was measured at the midpoint between the top of the iliac crest and the lower margin of the last palpable rib in the midaxillary line, with the patient in the standing position. After a 10 minutes rest period, blood pressure was measured three times; the average of the second and third blood pressure measurements was used for the analysis.
Biochemical measurements
In both studies, blood samples were obtained from an antecubital vein of each patient after a 12 hours overnight fast and 20 minutes in a sitting position. Plasma glucose, total cholesterol, TG, and HDL-C were measured using standardized procedures (Roche Diagnostics GmbH, Mannheim, Germany). Low density lipoprotein cholesterol was estimated by using the De Long et al. formula.[29] Accuracy and precision of lipid measurements in our laboratory are under periodic surveillance by the Center for Disease Control and Prevention service (Atlanta, GA, USA). Inter-assay coefficients of variation were < 6% for all these assays. High-sensitivity C-reactive protein (hs-CRP) was determined by immunonephelometry on a BN Pro Spec nephelometer (Dade Behring, Marburg, Hesse, Germany), according to the manufacturer’s method, with intra- and inter-assay variation coefficients below 3%. Adiponectin, leptin, interleukin 6 (IL-6), interleukin 1β (IL-1β), monocyte chemoattractant protein-1 (MCP-1), plasminogen activator inhibitor-1 (PAI-1), and insulin were quantified using a Bio-Plex system (Bio-Rad Inc, Hercules, CA, USA), with intra- and inter-assay variation coefficients below 4% and 5% (respectively), for all these assays. Adiponectin/leptin ratio was calculated as indicator of adipose tissue functionality. [22, 30] The homeostatic model assessment of insulin resistance (HOMA-IR) index was calculated using the formula: HOMA-IR = (Glucose [mmol/l] × Insulin [µIU/l]) / (22.5).
Computed tomography (CT)
Participants in the control group from the GEA study underwent CT to determine the presence of fatty liver or arterial calcium score > 100 Agatston units as a subrogated of subclinical atherosclerosis. CT is a validated method for measuring coronary artery calcium [31] and non-alcoholic fatty liver disease.[32] CT of the chest was performed using a 64- channel multi-detector helical computed tomography system (Somatom Sensation, Siemens), and images were interpreted by experienced radiologists. Scans were read to determine coronary artery calcification scores using the Agatston method. [31] To determine the liver and spleen attenuation, a single slice CT scan was obtained at the level of T11–T12 or T12–L1. [32]
Electrical bioimpedance
Body composition was determined in participants of the second phase of the study, using a 6-frequency electrical bioimpedance analyzer (InBodyS10, Korea), in the supine position. Total body fat percentage and visceral fat area (cm2) were obtained. All study participants were asked to adhere to a fast of at least 4 hours before measurement and to have no metal or electronic equipment with them during the test.
Cardiometabolic abnormalities definition
In all subjects, diabetes was defined as glucose > 7.0 mmol/L, hypoglycemic treatment or previous diagnosis, [33] non-alcoholic fatty liver disease as spleen-liver attenuation ratio < 1.0, [32] hypertension as self-reported treatment with antihypertensive medications or systolic/diastolic blood pressure ≥ 140/90 mmHg [34] and subclinical atherosclerosis as coronary artery calcium > 100 Agatston units. [31]
Adipose tissue biopsies
In the 68 healthy volunteers enrolled in the study of the morpho-functionality of adipose tissue, a subcutaneous white adipose tissue sample was obtained from periumbilical fat, with surgical technique under local anesthesia (2% lidocaine), and after an overnight fast.[35] Biopsies were immediately rinsed with sterile saline and visible blood vessels were removed with sterile tweezers. Adipose tissue biopsies were fractionated in three pieces. Two pieces were stored in cryotubes and immediately frozen in liquid nitrogen. The last piece was immediately fixed in PBS-buffered 4% paraformaldehyde for histological analyses.
Analysis of the number of adipocytes per field and adipocytes mean area
After 24 hours, fixed tissues were dehydrated in ethanol, cleared in xylene, embedded in paraffin, sectioned at 4 µm, and stained with hematoxylin and eosin. Digital images were obtained using a digital camera (Leica ICC50 HD) coupled to Leica DM750 microscope using a 20X lens at a resolution of 2048 × 1536 pixels using LAS EX V 3.0 software (Leica Microsystems, Heerbrugg Switzerland). Five fields of view image captures were taken per slide in varying parts of the fat biopsies of each patient. Adipocytes size was measured by converting pixels into microns utilizing Adiposoft (ImageJ) software with the following parameters: minimum diameter 10 µm, maximum diameter 1000 µm and microns per pixel 0.439 as previously described. [36] The mean adipocytes area (µm2) and number of adipocytes (per field) were determined in a blinded manner by taking into account all the adipocytes of the 5 fields of each patient. After completing the automated analyzes of the counting and area adipocytes measurements, each value in the data table was checked manually to ensure that it represented a single adipocyte in order to prevent errors that may occur during automated analyzes.
Dysfunctional adiposity index (DAI)
The original formula of VAI is separately calculated in men or women and includes two sections.[23] The first section represents the proportion of central adipose tissue through waist circumference with respect to BMI as a measure of generalized fat. With the analysis of the linear regression between both adiposity markers, the constants of the intercepts and the slopes can be determined, which are included in the formula: WC/ [constant intercept + (slope constant ∗ BMI)]. The second section represents AT function through TG and HDL-C plasma concentration assessment. In this case, median values of TG (mmol/L) and HDL-C (mmol/L) from healthy subjects without cardiometabolic risk factors are used as a reference. Thus, the final equation for DAI estimation by gender is represented as:
$$\text{DAI}\text{ =}\left(\frac{\text{WC}}{\text{constant intercept + }\left(\text{constant slope * BMI}\right)}\right)\text{ }\left(\frac{\text{TG}}{\begin{array}{c}\text{ }\text{ }\text{T}\text{G}\text{ }\text{m}\text{e}\text{d}\text{i}\text{a}\text{n}\ \text{r}\text{e}\text{f}\text{e}\text{r}\text{e}\text{n}\text{c}\text{e}\text{ }\text{ }\text{ }\end{array}}\right)\left(\frac{\begin{array}{c}\text{H}\text{D}\text{L}\text{-}\text{C}\text{ }\text{m}\text{e}\text{d}\text{i}\text{a}\text{n}\ \text{r}\text{e}\text{f}\text{e}\text{r}\text{e}\text{n}\text{c}\text{e}\end{array}}{\text{ HDL-C }}\right)$$
Data are presented as means ± standard deviation, median (interquartile range), and the number of subjects (percentages). In participants of the morpho-functionality adipose tissue study, HOMA-IR was calculated and subjects were stratified according to those below or above the median HOMA-IR value. A comparison of these groups was accomplished with Student’s t, Mann–Whitney U, or Chi-square tests; respectively. Spearman coefficient correlation was calculated to evaluate the relationship of DAI with inflammatory and morpho-functional characteristics of adipose tissue. The discriminative power of the DAI and other adiposity measurements, to identify subjects with HOMA-IR above the median, was tested by calculation of the area under the receiver operating characteristic (ROC) curve (area under the curve [AUC]); the optimal cut-off point was determined by the maximal Youden index. In order to investigate the association of high DAI values with cardiometabolic abnormalities in the control group of the GEA study, a multiple logistic regression stepwise forward analysis was performed adjusting by potential confounders. The values of regression analysis are shown as odds ratio (95% interval of confidence). All p values < 0.05 were considered statistically significant. The statistical analyses were performed using 15.0 software SPSS (Chicago, IL, USA) and SAS JMP ® Trial version 15.1.0 (Cary, NC, USA).