Next Article in Journal
Production of d-Tagatose by Whole-Cell Conversion of Recombinant Bacillus subtilis in the Absence of Antibiotics
Next Article in Special Issue
Angiogenin Levels and Their Association with Cardiometabolic Indices Following Vitamin D Status Correction in Saudi Adults
Previous Article in Journal
Medicinal-Cosmetic Potential of the Local Endemic Plants of Crete (Greece), Northern Morocco and Tunisia: Priorities for Conservation and Sustainable Exploitation of Neglected and Underutilized Phytogenetic Resources
Previous Article in Special Issue
Free Fatty Acid Species Differentially Modulate the Inflammatory Gene Response in Primary Human Skeletal Myoblasts
 
 
biology-logo
Article Menu

Article Menu

Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Replication of Integrative Data Analysis for Adipose Tissue Dysfunction, Low-Grade Inflammation, Postprandial Responses and OMICs Signatures in Symptom-Free Adults

by
Esther C. Gallegos-Cabriales
1,†,
Ernesto Rodriguez-Ayala
2,†,
Hugo A. Laviada-Molina
3,
Edna J. Nava-Gonzalez
4,
Rocío A. Salinas-Osornio
5,
Lorena Orozco
6,
Irene Leal-Berumen
7,
Juan Carlos Castillo-Pineda
8,
Laura Gonzalez-Lopez
5,
Claudia Escudero-Lourdes
9,
Judith Cornejo-Barrera
10,
Fabiola Escalante-Araiza
2,6,
Eira E. Huerta-Avila
6,
Fatima A. Buenfil-Rello
11,
Vanessa-Giselle Peschard
2,
Eliud Silva
2,
Rosa A. Veloz-Garza
1,
Angelica Martinez-Hernandez
6,
Francisco M. Barajas-Olmos
6,
Fernanda Molina-Segui
3,
Lucia Gonzalez-Ramirez
3,
Ruy D. Arjona-Villicaña
3,
Victor M. Hernandez-Escalante
3,
Janeth F. Gaytan-Saucedo
11,
Zoila Vaquera
11,
Monica Acebo-Martinez
9,
Areli Murillo-Ramirez
8,
Sara P. Diaz-Tena
8,
Benigno Figueroa-Nuñez
12,
Melesio E. Valencia-Rendon
5,
Rafael Garzon-Zamora
5,
Juan Manuel Viveros-Paredes
5,
Salvador B. Valdovinos-Chavez
1,
Anthony G Comuzzie
13,
Karin Haack
11,
Ashley A. Thorsell
14,
Xianlin Han
15,
Shelley A. Cole
11 and
Raul A. Bastarrachea
11,14,*
add Show full author list remove Hide full author list
1
Facultad de Enfermería, Universidad Autónoma de Nuevo León (UANL), Monterrey 64460, Mexico
2
Centro de Investigación en Ciencias de la Salud (CICSA), Facultad de Ciencias de la Salud, Universidad Anáhuac Norte, Lomas Anahuac 52786, Mexico
3
Escuela de Ciencias de la Salud, Universidad Marista de Mérida, Mérida 97300, Mexico
4
Facultad de Salud Pública y Nutrición (FASPyN), UANL, Monterrey 64460, Mexico
5
Departamento de Nutrición, Universidad del Valle de Atemajac (UNIVA), Zapopan 45050, Mexico
6
Laboratorio de Inmunogenómica y Enfermedades Metabólicas, Instituto Nacional de Medicina Genómica, SS, Ciudad de México 14610, Mexico
7
Facultad de Medicina y Ciencias Biomédicas, Universidad Autónoma de Chihuahua, Chihuahua 31125, Mexico
8
Departamento de Nutrición Humana, Universidad Latina de América, Morelia 58170, Mexico
9
Facultad de Ciencias Químicas, Universidad Autónoma de San Luis Potosí, San Luis Potosi 78240, Mexico
10
Departamento de Enseñanza, Postgrado e Investigación, Hospital Infantil de Tamaulipas, Ciudad Victoria 87150, Mexico
11
Population Health Program, Southwest National Primate Research Center (SNPRC), Texas Biomedical Research Institute, San Antonio, TX 78227-0549, USA
12
Clínica de Enfermedades Crónicas y Procedimientos Especiales (CECYPE), Morelia 58249, Mexico
13
The Obesity Society, Silver Spring, MD 20910, USA
14
Sansum Diabetes Research Institute, Santa Barbara, CA 93105, USA
15
Department of Medicine, Sam and Ann Barshop Institute for Longevity and Aging Studies, University of Texas Health San Antonio, San Antonio, TX 78229, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Biology 2021, 10(12), 1342; https://doi.org/10.3390/biology10121342
Submission received: 21 October 2021 / Revised: 11 December 2021 / Accepted: 12 December 2021 / Published: 16 December 2021

Abstract

:

Simple Summary

There is an alarming increase of cardiovascular disease (CVD) and type 2 diabetes (T2D) in Mexican nationals and Mexican Americans. Studying adipose tissue (AT) dysfunction and early biomarkers of cardiovascular and immunometabolic risk in Mexican nationals may have a strong impact on future public health policies for US-born Mexican Americans and other populations of Mexican origin in the US. The goal of this study is to evaluate the early transition towards healthy/unhealthy adipose tissue expansion to identify AT dysfunction through systemic, molecular and OMICS measures in the fasting and fed states in symptom-free volunteers with no history of age-related chronic diseases in support of precision medicine and discovery.

Abstract

We previously reported preliminary characterization of adipose tissue (AT) dysfunction through the adiponectin/leptin ratio (ALR) and fasting/postprandial (F/P) gene expression in subcutaneous (SQ) adipose tissue (AT) biopsies obtained from participants in the GEMM study, a precision medicine research project. Here we present integrative data replication of previous findings from an increased number of GEMM symptom-free (SF) adults (N = 124) to improve characterization of early biomarkers for cardiovascular (CV)/immunometabolic risk in SF adults with AT dysfunction. We achieved this goal by taking advantage of the rich set of GEMM F/P 5 h time course data and three tissue samples collected at the same time and frequency on each adult participant (F/P blood, biopsies of SQAT and skeletal muscle (SKM)). We classified them with the presence/absence of AT dysfunction: low (<1) or high (>1) ALR. We also examined the presence of metabolically healthy (MH)/unhealthy (MUH) individuals through low-grade chronic subclinical inflammation (high sensitivity C-reactive protein (hsCRP)), whole body insulin sensitivity (Matsuda Index) and Metabolic Syndrome criteria in people with/without AT dysfunction. Molecular data directly measured from three tissues in a subset of participants allowed fine-scale multi-OMIC profiling of individual postprandial responses (RNA-seq in SKM and SQAT, miRNA from plasma exosomes and shotgun lipidomics in blood). Dynamic postprandial immunometabolic molecular endophenotypes were obtained to move towards a personalized, patient-defined medicine. This study offers an example of integrative translational research, which applies bench-to-bedside research to clinical medicine. Our F/P study design has the potential to characterize CV/immunometabolic early risk detection in support of precision medicine and discovery in SF individuals.

Graphical Abstract

1. Introduction

Certain individuals can be considered metabolically healthy despite their long-standing and high degree of body fat accumulation. This effect can be thought of as healthy AT expansion. On the other hand, unhealthy AT expansion is a major contributor to the systemic metabolic disturbances that are characteristic of obesity and type 2 diabetes mellitus (T2DM) [1]. The loss of expansion capacity can occur in patients with normal weight, explaining the existence of MUH lean subjects [2,3].
Maintaining a body mass index (BMI) in the ideal range of 20.0–25.0 kg/m2 may effectively reduce the risk of early death from cardiovascular and immunometabolic disease [4]. In addition, subjects exhibiting MH obesity (MHO) are not at increased risk of cardiovascular events compared to MH people in the normal weight range. It appears that, independent of BMI, AT dysfunction [5] triggers early molecular events that ultimately lead both to insulin resistance (IR) and also to low-grade chronic subclinical inflammation (LGCSI) and subsequent morbidity and mortality. The recently developed ALR correlates with IR better than adiponectin or leptin alone [6]. This emerging biomarker decreases with increasing number of cardiometabolic risk factors, reflecting the functionality of AT, and negatively correlates with markers of LGCSI [7]. The ALR has been suggested as a significant marker of AT dysfunction and inflammation [8,9].
Cardiovascular Disease (CVD) and T2DM are serious public health problems among Mexican nationals [10] and Mexican Americans [11]. Identifying MH/MUH parameters occurring simultaneously with AT dysfunction markers in Mexican nationals will enhance the understanding of early underlying biology of T2DM and CVD [12,13]. The knowledge gained will have a strong impact on future public health policies to decrease the burden of T2DM and CVD in Hispanics, Mexican Americans and populations of Mexican origin in the US.
The GEMM family study (Genética de las Enfermedades Metabólicas en México, Genetics of Metabolic Diseases in Mexico) is a bi-national, multi-center collaborative study of cardiovascular risk phenotypes of immunometabolic origin (CVRIMO) related to the risk of T2DM and CVD. The GEMM study is a longitudinal, family-based, precision medicine research project. In addition, known relatedness provides one level of informed adjustment for admixture, an important consideration given the diverse ancestry of modern Mexicans and many other cohorts. GEMM has acquired data directly from three tissues that are highly relevant to immunometabolism: AT, SKM and blood (immune system). Our measurements in blood and tissues were taken over a time course to allow fine scale multi-OMIC profiling of individual postprandial responses. Unlike studies utilizing tissue repositories of poorly-characterized origin, our study obtained tissues from a well-characterized cohort of symptom-free adults of common ancestry [14,15].
Precision and personalized medicine, linked to the identification of early risk and prevention instead of curative pathological symptoms, is rapidly taking place in the immunometabolic and CVD field. This paper aims to identify early biomarkers of immunometabolic dysregulation in persons without clinical symptoms, which may define a future MUH phenotype, and at characterizing the AT dysfunction phenotype at a postprandial systemic and molecular level in SF individuals from the GEMM study. Our working hypothesis is that by comparing the metabolically healthy vs. unhealthy phenotype in our SF participants and their ALR levels, we will be able to screen for presence or absence of systemic immunometabolic dysfunction. We will also be able to identify in them, within the normal variation, significant differences in circulating levels of proinflammatory cytokines and the presence of impaired postprandial peaks from hormones related to the insulin–glucose axis. This paper increases and replicates previous published data from 80 individuals, with a focus on subcutaneous AT profiling [15].

2. Methods

Recruitment of study participants: As of February 2020, GEMM had performed the meal challenge, tissue biopsies and sample collection for 225 participants [14,15]. Unfortunately, the COVID-19 pandemic forced us to temporarily stop recruitment. However, we were able to bring samples from 124 study participants to the US in December 2019 and have obtained data from these recruited participants (81 females and 43 males). Our cohort is comprised of SF volunteers that were carefully screened to rule out presence of any chronic or acute inflammatory disease, as described in detail elsewhere [14,15,16].
Mixed Meal Challenge: GEMM participants were given an innovative balanced mixed-meal challenge [17]. We provided a mixed meal in the form of Ensure® Plus liquid shake from Abbott Laboratories (macronutrient composition: 65% carbohydrate, 15% protein and 20% fat) dosed at 30% of each participant’s resting energy expenditure in Kcal/day. This dosage mimics the first morning meal after an overnight fast and allows each participant to receive his/her own ideal calorie amount as dictated by his/her own resting metabolic rate and expenditure. Our goal was to compare the metabolic response to a fixed macronutrient composition among each participant [18]. We chose a liquid shake instead of a solid meal for easier means to administer the mixed meal challenge.
Phenotyping Methodology: We recruited study participants according to established guidelines and strategies [19]. We chose to initially classify our SF participants with the presence/absence of adipose tissue dysfunction into two groups according to their ALR cut-off point > 1 or <1 [7,9]. The different experimental groups were established through adiponectin and leptin circulating measures to determine the adiponectin/leptin ratio (ALR) trait and to characterize adipose tissue dysfunction. The rationale for this approach is that a low ALR as a marker of AT dysfunction is characterized by a lower secretion of adiponectin in relation to leptin levels, triggering unhealthy adipose tissue hypoxia, proinflammatory macrophage polarization traits, altered adipokine profile and insulin resistance (IR). Bioimpedance was measured by Tanita BC—418 Body Composition Analyzer, and body composition by dual energy X-ray absorptiometry. A wide range of clinical biochemistries, hormones, cytokines and endophenotypes were measured for F/P metabolic assessments through analyte categories, including β-cell and the insulin–glucose axis, inflammation endophenotypes and lipid–lipoprotein metabolism. These biochemical phenotypes were analyzed on a Luminex 100IS platform and an Immulite 1000 which runs enzyme-linked immunosorbent assay (ELISA) and radioimmunoassay (RIA) analyses.
Shotgun Lipidomics: Lipidomics analysis serves as a powerful tool for understanding the biochemical and cellular mechanisms underlying lipid-related immunometabolic disease processes by quantifying the changes of individual lipid classes, subclasses and molecular species [20]. Plasma samples were transported to Dr. Xianlin Han’s laboratory at the University of Texas Health Sciences Center San Antonio [21]. Lipids were extracted with chloroform/methanol in the presence of a cocktail of internal standards [22] and analyzed by multi-dimensional mass spectrometry-based shotgun lipidomics (MDMS-SL), as previously described [14,15,23].
miRNA Sequencing from plasma exosomes: Exosomal miRNA was isolated from 1 mL of plasma using the Plasma/Serum Exosome Purification and RNA Isolation Mini Kit (Norgen, Thorold, ON, Canada). The quality and quantity of miRNA were assessed using the Low Abundance RNA Quantification Kit (Norgen). Conversion of miRNA into a cDNA library was performed using the NEXTflex Illumina Small RNA-Seq Kit v3 (Bio Scientific, Phoenix, AZ, USA). The purification and enrichment of cDNA by PCR was performed to create cDNA libraries. Forty-eight sample libraries were pooled and quantified using the Kapa Library ABI Prism Quantification kit. Library pools were sequenced using a commercial service.
Methodology RNA isolation: Total RNA was isolated from pre- and post-prandial samples for each tissue type. Isolation from adipose and muscle shock frozen tissue used a Direct-zol RNA Miniprep Plus kit w/TRI reagent (Zymo Research, Tustin, CA, USA). Tissue was homogenized in TRI reagent using a Beadbeater (Biospec Products, Bartlesville, OK, USA). Total RNA was isolated from blood samples using the TempusTM spin RNA isolation kit (ThermoFisher Scientific, Waltham, MA, USA). RNA quantity was determined using a QubitTM RNA BR assay kit (ThermoFisher Scientific), and RNA quality/integrity was assessed using the RNA screen tape assay (Agilent Technologies, Santa Clara, CA, USA). The GLOBINclearTM kit (ThermoFisher Scientific) was used to isolate globin mRNA from total RNA in blood.
RNASeq of RNA from fat and muscle: mRNA (from total RNA) was converted into a cDNA library using the Kapa mRNA HyperPrep kit for Illumina Platforms and unique dual-indexed adapter kit (Kapa Biosystem, Wilmington, MA, USA). Library quality was checked using D1000 Screen Tapes (Agilent Technologies). Samples were quantified by qPCR using the Kapa Library ABI Prism Quantification kit (Kapa Biosystems) before pooling, and pools were quantified using the same kit after pooling 24 samples to load on one flow cell lane. Library pools were sequenced using a commercial service.
Multidimensional analysis platform development: For this paper, we used FunRich, an open-access, standalone functional enrichment and network analysis tool to perform analysis on background databases that are integrated from heterogeneous genomic and proteomic resources (>1.5 million annotations) [24,25] and the statistical software R version 4.1.1 (R Core Team, 2021, https://www.R-project.org/, accessed on 15 October 2021).
Statistical Analysis: Data are presented as mean ± standard deviation (SD), unless otherwise indicated. Differences between groups are analyzed by two-tailed unpaired Student’s t-tests, as appropriate. The calculations were performed using JMP® 14.2.0 (SAS Institute Inc., Cary, NC, USA) and GraphPad Prism 6 (GraphPad Software, Inc., La Jolla, CA, USA). Multivariate linear regression analysis was used to measure the strength of the relationship between our immunometabolic F/P phenotypes with the ALR. For statistical analysis of independent multiple testing, we used a one-way analysis of variance (ANOVA) and Tukey’s multiple comparison test. Our models incorporated the random effect of kinship and included age and sex as covariates. The analysis of the relative gene expression data was conducted by using the delta-delta Ct method to analyze the relative changes in gene expression from our real-time quantitative PCR experiments. This approach relates to the PCR signal of the target transcript in a treatment group to that of another sample such as an untreated control as has been described by Livak and Schmittgen [26]. One-way ANOVA was used to study the association between serum stimulation time and the level of gene expression. Fold increase in mRNA was the dependent variable.

3. Results

We chose to initially classify our SF participants with the presence/absence of adipose tissue dysfunction into two groups according to their ALR cut-off point > 1 or <1 [7,9]. We found interesting differences between the high (H)ALR and low (L)ALR groups. Table 1 and Table 2 show the demographic characteristics and immunometabolic parameters for the participants with a (H)ALR or a (L)ALR in our symptom-free (SF) GEMM cohort. Eighty-four individuals had an ALR > 1, and 40 individuals had an ALR < 1. There were significant and striking differences of non-traditional immunometabolic biomarkers [27,28] between individuals with a (H)ALR vs. a (L)ALR. Systemic levels of leptin, weight, waist circumference, % body fat, circulating proinflammatory levels of IL-6, plasminogen activator inhibitor (PAI)-1, hsCRP and fibrinogen and the 5 h postprandial triglyceride curves were significantly increased in the group with AT dysfunction. Our data also showed significantly lower levels of adiponectin. All phenotypes pertaining to the insulin–glucose axis (fasting and 2 h insulin levels, 2 h postprandial glucose levels, fasting glucagon, the Matsuda and the HOMA-IR indexes and curves of 5 h glucose and GLP-1) in the (L)ALR category were impaired when compared with the (H)ALR group (Table 1 and Table 2).
We estimated the presence/absence of metabolic health (MH and MUH) among the SF GEMM subjects (according to their ALR cut-off point > 1 or <1) using the proposed criteria for defining MH/MUH. This was based on the number of metabolic syndrome parameters present and an assessment of insulin sensitivity (using the Matsuda index score [29]) and LGCSI [30] as determined by hsCRP levels [4,31,32] (Figure 1). We identified three subgroups of symptom-free individuals (Figure 1): (i) Group 1: subjects with a (H)ALR having <2 of the metabolic syndrome risk factors and the absence of IR and LGCSI; (ii) Group 2: subjects with a (H)ALR having ≥2 MUH phenotypes for metabolic risk and the presence of IR and LGCSI; and (iii) Group 3: the same characteristics as Group 2 but with a (L)ALR (Table 3). The proposed criteria to identify MH individuals was based on cutoff values for a healthy cardiometabolic profile [4,29,32,33,34,35]. Reference values in Mexican Americans have been determined with a HOMA-IR > 3.80 as having clear correlates of insulin resistance [36].
We performed shotgun lipidomics and a comprehensive sequencing of exosomal miRNA study using fasting and postprandial plasma samples from a subgroup of 14 females from our SF cohort. Table 4 shows their clinical characteristics. We carefully matched this subgroup by age and fat% from our total SF female cohort (N = 81). We performed whole shotgun lipidomics in fasting and 3 h postprandial plasma. Using these lipidomics results, we classified the 14 females having a high (H; mean 2.2) or low (L; mean 0.5) ALR, and compared the differential signatures of lipid classes, subclasses and molecular species. This strategy allowed identification and quantification of an immense quantity of individual lipid species and classes, including the following: triglycerides, diglycerides, monoglycerides, cholesterol and cholesteryl esters, oxysterols, acylcarnitines, acyl-CoAs, phospholipids (including cardiolipin), lysophospholipids, eicosanoids, 4-hydroxy alkenals and retinoic acid [37]. Shotgun lipidomics and miRNA differential expression profiles are shown in Figure 2, Figure 3, Figure 4 and Figure 5. The implications of these data are explained in the Section 4.
We compared F/P gene expression profiles of two key tissues (SKM and SQAT) between MH and MUH individuals with or without AT dysfunction in 6 SF volunteers (Table 5), chosen from our previous 14 SF adult females, matched by gender and fat% (Table 4).
We assigned participants with a (L)ALR (mean 0.52) < 1 as the at-risk group (RG) and participants with a (H)ALR (mean 1.98) > 1 as the control (low risk) group (CG) (Table 5). The contrast linear combination of variables between the CG and RG are shown in Table 6. Figure 6 and Figure 9 report differentially expressed genes in SKM (twenty genes) and AT (twenty genes) between the CG vs. RG. By comparing the magnitude of change vs. the statistical significance (p-value) we were able to allegedly identify the most biologically significant genes with large fold changes that are statistically significant. Figure 7 and Figure 10 show heat maps with differences in gene expression from SKM and AT signatures from the 20 genes identified for each tissue as shown in Figure 6 and Figure 9. Figure 8 and Figure 11 show the enrichment analysis and interaction networks obtained through the FunRich database analysis tool. The explanation of the associations and gene interactions is written in the Section 4.

4. Discussion

In this study, we used a precision medicine screening approach anchored by measurements of AT function and MH/MUH phenotypes among active, symptom-free (SF) adults to identify integrative strategies for detection of early immunometabolic and cardiovascular risk. Our cohort was chosen with a careful exclusion criterion applied for participants with acute illness, activity-limiting unexplained illness, absence of metabolic syndrome criteria, hypertension, dyslipidemia, prevalent diabetes, known cardiovascular or chronic lung disease, cancer or renal failure [16]. We decided to characterize our SF cohort with the presence/absence of AT dysfunction through a novel biomarker, the adiponectin/leptin ratio (ALR) [6,7,9].
As shown in Table 1 and Table 2, we found that one group (N = 84) had a high (H)ALR (6.7 ± 10.0) that translated to a healthy AT expansion and the second group (N = 40) had a low (L)ALR (0.6 ± 0.3) that translated to an unhealthy AT expansion. We found interesting differences between (H) vs. (L)ALR when we characterized AT dysfunction. The anthropometric variables, the β-cell and insulin–glucose axis phenotypes, the adipose tissue hormones, the lipid–lipoprotein and the proinflammatory phenotypes were significantly altered in the (L)ALR group. These findings represent a solid indication that there is an AT metabolic imbalance measured through the plasma ALR in our GEMM SF subjects. However, we should keep in mind that males and females may exhibit differences in adiponectin/leptin ratios and metabolic dysfunction, given that 81 participants were females out of our 124 recruited SF subjects. Our findings also highlight the importance of postprandial immunometabolism through measurements of dynamic phenotypes. After careful screening of AT dysfunction through the ALR, we estimated the presence/absence of MH and MUH risk among our SF GEMM subjects as shown in Figure 1 and Table 3. We were able to combine the MH/MUH phenotype with the AT dysfunction phenptype in our SF cohort. In order to clasify our 124 SF subjects as metabolically healthy, we used well-defined criteria: absence of abdominal obesity based on waist circumference, absence of metabolic syndrome components, normal lipid values, normal fasting glucose concentrations, normal C-reactive protein concentrations and absence of IR as evaluated by the Matsuda index (Table 3) [4,32,38,39].
Most studies regarding metabolic health have focused their reports on individuals with obesity. They have designated that metabolically healthy obesity (MHO) is applied to those who meet the BMI cutoff point of ≥30 kg/m2, but do not have other major cardiovascular risk factors and who are not at higher cardiovascular risk than nonobese individuals [40]. One major weakness with respect to this phenotype is that perhaps it solely represents a subgroup of obese subjects with just a lower quantity of metabolic risk than expected for their degree of BMI [41]. Moreover, there are subsets of individuals who are considered “normal/healthy” weight, categorized by a low BMI, but with an increased metabolic/cardiovascular risk. This subset, characterized by elevated insulin levels, elevated triglyceride levels and insulin resistance are much harder to characterize than the MHO group and are predisposed to subsequent development of T2DM and coronary artery disease [42,43,44]. Despite these controversies, we decided to estimate the presence/absence of metabolic health (MH/MUH) among the SF GEMM subjects according to their ALR cut-off points characterizing AT dysfunction. The analysis of the phenotypes associated with MH/MUH are combined for males and females due to their similarity.
Three subgroups of SF individuals were identified: one group with only one metabolic risk phenotype and the absence of IR and LGCSI with a high (H) mean ALR of 9.5 ± 13.0, and a second and third group with altered fasting glucose levels (>100 mg/dL), low HDL-C values, high triglyceride levels, visceral obesity and isulin resistance with a mean ALR of 4.0 ± 4.5 and 0.6 ± 0.3, respectively (Table 3, Figure 1). Twenty-eight individuals fulfilled the criteria of MH with the presence of only one metabolic risk phenotype. We identified a mid (M) “gray zone” where a second group of GEMM SF individuals (N = 43) had a mean ALR of 4.0 ± 4.5 and presented between 2 to 7 metabolic risk phenotypes. For the third group (N = 40) with a mean (L)ALR of 0.6 ± 0.3 the participants presented between 2–8 metabolic risk factors. Table 3 and Figure 1 demonstrate these data by presenting the risk criteria cut-offs, the number of cardiometabolic risk phenotypes present in everyone and the percentage of these risk phenotypes within each of the three groups with AT dysfunction. Compared to the group with a (H)ALR, subjects from the (M) and (L)ALR groups had a much larger and more significant prevalence of IR and prediabetes.
When comparing the three different groups for the presence of AT dysfunction through a (H), (M) and (L)ALR, a clear picture emerges. There is a gradual and ascendant increase of the prevalence of their dysglycemic state, their IR state, increased triglycerides, their proinflammatory state and prediabetic state that can be easily correlated with their AT function. Through the biology and molecular actions of adiponectin and leptin, the ALR seems to offer a means to translate clinical biochemistry features into cellular steps for the transition between a healthy and unhealthy AT expansion in SF individuals. We could speculate that the (H)ALR group in our GEMM SF cohort represents a mild physiologic hypoxia during early healthy AT expansion, appropriate signaling for sufficient angiogenesis and small initial inflammatory cell infiltration. This induces an organized extracellular matrix (ECM) remodeling process for adequate expansion, involving functional adipocytes and the presence of predominantly M2 macrophages without fibrosis (Graphical Abstract).
In the presence of an excessive calorie-dense, palatable macronutrient intake there is adipocyte growth, which ultimately changes the healthy AT expansion through a progressive increase in ECM component synthesis, leading to a discrete and polarized macrophage phenotypic conversion, which characterizes the (M)ALR group. Chronic high caloric intake over time results in adipocyte hypertrophy, hypoxia and impaired ECM synthesis. The deleterious inflammatory response is fully triggered, predominantly by the presence of M1-polarized macrophages secreting a vast array of proinflammatory cytokines (mainly IL-6 and TNF-α). Due to severe hypoxia, a vicious cycle of persistent injury signaling and resultant impaired angiogenesis limits oxygen supply, subsequently triggering CLGSCI and massive proinflammatory macrophage infiltration, which characterizes the (L)ALR group (Graphical Abstract).
This scenario sets the landscape for excessive fibrosis, promoting significant AT hypertrophy, triggering the associated systemic immunometabolic dysfunction and leading to deleterious metabolic consequences such as T2DM, impaired angiogenesis, systemic inflammation, IR, endothelial dysfunction, altered lipid–lipoprotein metabolism and atherosclerosis [45,46] This same scenario could also be applicable to SF individuals with a (L)ALR independent of their AT accumulation and BMI measurements. It seems clear that AT dysfunction as a proxy for AT unhealthy expansion, measured through the ALR, is the common denominator that dictates early transition to immunometabolic impairments such as LGCSI, IR, postprandial proinflammatory response leading to dysglycemic states, prediabetes, lipid–lipoprotein abnormalities, and endothelial dysfunction phenotypes (Graphical Abstract). It appears that these phenotypes pertain to an approach that has not been exploited in full within the routine clinical setting [47,48]. These phenotypes are also considered CVD risk factors of immunometabolic origin [15], termed as nontraditional CVD risk factors (Table 1 and Table 2) [48,49]. It is also worth mentioning that this scenario may have individualized permissive factors and a strong genetic background to affect different individuals with different thresholds of fat accumulation. This has been implied in animal models where scientists have shown that there are potential differences in lipolysis by demonstrating that leptin-induced lipolysis may be regulated by nitric oxide (NO) [50].
As described in our results, we used heat maps to provide a descriptive representation of the postprandial comprehensive sequencing of miRNA in F/P plasma exosomal samples and a shotgun lipidomics dynamic process at 180′ to compare differences between female participants (N = 14) with (H) and (L)ALR (Table 4). We are only presenting the postprandial results 3 h after food ingestion, representing the peak of digestion. The lipidomic and miRNA heat maps show a fold increase of molecular lipid species and miRNA differential expression at time point 180′ for (L)ALR individuals as compared to (H)ALR ones (Figure 2 and Figure 4). Our preliminary results interestingly showed that the “hot brown” highest values in shotgun lipidomics and miRNA differential expression were found in participant 16 with a (L)ALR, and participant 21 with a (H)ALR. We performed volcano plots that enabled identification of lipid signatures with large fold changes, with the most biologically significant lipid species differentially expressed, also identifying miRNAs signatures with a trend for downregulation (Figure 3 and Figure 5). Figure 3 shows the specific classes of lipids differentially expressed between the (H) and (L)ALR at 180 min. These classes were PI [51], AC [52], Cer [53] and TAG [54]. The miRNA identified were hsa-miR-335-3p [51], hsa-miR-483-5p [52], hsa-miR-505-3p [53] and hsa-miR-584-5p [54] (Figure 5).
Female participants 16 and 21 seemed to have had the highest lipidomics and miRNA activities and expression profiles (Figure 2 and Figure 4). Participant 16 is a symptom-free metabolically unhealthy 35 year old female with AT dysfunction (ALR = 0.20), insulin resistance (Matsuda index of 2.7, HOMA-IR of 4.06), chronic subclinical inflammation (hsCRP of 51.8 mg/L) and HDL-C of 45 mg/dL, with a high risk of developing T2DM. Participant 21 is a symptom-free metabolically unhealthy 30 year old female with an ALR = 1.41, a Matsuda index of 2.6, fasting glucose levels of 101 mg/dL and 2 h-postprandial glucose levels of 155 mg/dL and fasting triglycerides of 166 mg/dL, which translates to a strong predisposition of developing a dysglycemic state [33]. Although our sample is too small for rigorous conclusions, these results illustrate examples of how we can characterize our future full dataset to screen plasma circulating molecular biomarkers and identify early signatures of risk for immunometabolic and cardiovascular disease and achieve a translational (from bench to bedside) and personalized medicine approach.
Finally, we obtained SKM and SQAT biopsies from six SF volunteers chosen from our previous fourteen SF adult females to perform whole-transcriptome analysis with total RNA sequencing (RNA-Seq). We obtained the F/P SKM (vastus lateralis) and SQAT biopsies from the thigh of each SF participant. We obtained our subcutaneous (SQ) tissue biopsy from the thigh instead of the abdominal SQ depot, and it has been documented that adipocyte formation is higher in the gluteofemoral compared to the abdominal SQ adipose depot when in vivo adipogenesis is directly measured [55,56]. In addition, thigh SQAT mass (but not abdominal SQAT mass) is weakly associated with a lower prevalence of metabolic syndrome and cardiometabolic risk [57]. Support for these findings is provided by studies showing that the relative release of palmitoleate, an insulin-sensitizing lipokine, is markedly higher from SQ gluteofemoral AT compared with SQ abdominal AT. Furthermore, by directly measuring adipogenesis, important information has been documented regarding preadipocyte and adipocyte formation being higher in the femoral compared to the abdominal SQ adipose depot [58,59,60]. Although this study comprises a small sample size (N = 6), the purpose was to obtain proof of concept data, which we performed by repeated measurements of key immunometabolic phenotypes that are related to body composition: adipose tissue biology, the insulin–glucose axis proinflammatory biomarkers and postprandial metabolism. This is a subset of only six female participants matched by gender and percentage of body fat, independent of BMI, weight or age but divided by (H) vs. (L)ALR; therefore, a weakness of the analysis is that perhaps a power analysis should have been done instead given the small sample size and the number of phenotypes compared between the (H) and (L)ALR groups. Although our results showed several not significant p values and some phenotypes, specifically BMI, showed no differences between the compared groups, we speculate that with a larger sample size, differences could have been detected.
The contrast of the linear combination of phenotypes between the mean (H) and (L)ALR groups showed that the ALR, IL-6, hsCRP and the AUC for insulin were significantly higher in the L(ALR) group, perhaps indicating the presence of a higher proinflammatory activity in the (L)ALR group (Table 6). We assigned participants with a (L)ALR (N = 2, mean 0.52, Table 5) < 1 as the at-risk group (RG) and participants with a (H)ALR (N = 4, mean 1.98, Table 5) > 1 as the control group (CG). The logarithmic fold change (FC) for RNA was calculated for 0′ and 180′ (F/P). We also contrasted the two groups (RG vs. CG) using the natural logarithmic FC to identify genes ranked by unadjusted p-values to obtain a probability measurement. Figure 6 shows that the three most relevant genes differentially expressed in SKM were CCT6A [61,62], PMP22 [63] and TSC22D1 [64]. Interestingly, Figure 7 shows that the hottest activity detected in the heatmap for SKM gene expression was CCT6A for adult pathways involving the immune system (B cells, CDK4, CDK8, monocytes and NK cells). Figure 8 shows that TSC22D1 was the key centered gene associating a vast array of interactive pathways obtained through FunRich annotations for SKM transcript expression.
Figure 9 shows that the two most relevant genes differentially expressed in AT were WNT5A [65] and FNDC1 [66,67]. Figure 10 demonstrates that the MCM4 gene, involved with adult pathways involving the immune system (B cells, CDK4, CDK8, monocytes and NK cells), had the hottest activity detected in the heatmap for SQAT gene expression. Figure 11 shows that MCM4 was the key centered gene associating a vast array of interactive pathways obtained through FunRich annotations for SQAT transcript expression [68]. Overall, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 demonstrate a molecular workflow for bioinformatics analysis and provide a stepwise strategy to translate key transcriptomic/lipidomic expression into biological pathways and gene–gene expression interactions and associations of SKM and AT in the fasting and postprandial state.

5. Conclusions

Our study characterized F/P nonconventional biomarkers in blood and directly measured genomic and molecular pathways of F/P AT inflammation and other features of dysfunctional AT. We combined integrated, clinical and multiOMIC profiles from blood, AT and SKM, three tissues key to immunometabolic biology, and outlined novel methods and procedures in SF adults to identify biochemical pathways and potential regulatory networks. However, our data are cross-sectional, not longitudinal, therefore, it remains to be seen that the identification of early biomarkers are truly associated with risk factors with future abnormal metabolic symptoms. Nevertheless, by characterizing some of the earliest biomarkers of immunometabolic dysregulation, we can identify which subjects may have an unhealthy metabolism, therefore identifying risk factors that lead to the development of T2DM and CVD.

Author Contributions

Conceptualization, E.C.G.-C., E.R.-A., H.A.L.-M., S.A.C. and R.A.B.; Data curation, E.S., A.M.-H., V.M.H.-E., J.M.V.-P. and K.H.; Formal analysis, E.R.-A., E.S., A.M.-H., F.M.B.-O., V.M.H.-E., J.M.V.-P. and K.H.; Funding acquisition, H.A.L.-M., F.M.-S. and R.A.B.; Investigation, E.C.G.-C., E.R.-A., H.A.L.-M., E.J.N.-G., L.O., I.L.-B., J.C.C.-P., L.G.-L., R.A.S.-O., C.E.-L., J.C.-B., F.E.-A., E.E.H.-A., F.A.B.-R., V.-G.P., A.M.-H., F.M.B.-O., F.M.-S., L.G.-R., R.D.A.-V., V.M.H.-E., Z.V., M.A.-M., A.M.-R., S.P.D.-T., B.F.-N., M.E.V.-R., R.G.-Z., S.B.V.-C., A.G.C., X.H., S.A.C. and R.A.B.; Methodology, E.C.G.-C., E.R.-A., E.J.N.-G., L.O., J.F.G.-S. and R.A.B.; Project administration, R.A.V.-G., M.E.V.-R. and S.B.V.-C.; Resources, J.C.C.-P., F.A.B.-R., L.G.-R., S.P.D.-T., B.F.-N., M.E.V.-R., R.G.-Z., J.M.V.-P. and S.B.V.-C.; Software, K.H.; Supervision, E.C.G.-C., L.O., A.G.C. and R.A.B.; Validation, E.E.H.-A., V.-G.P., L.G.-R., R.D.A.-V. and R.G.-Z.; Writing—original draft, R.A.B.; Writing—review & editing, A.A.T. All authors have read and agreed to the published version of the manuscript.

Funding

Research reported in this publication was supported by the National Institute of Diabetes and Digestive and Kidney Diseases of the National Institutes of Health under Award Number R56DK114703-01. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Institutional Review Board Statement

All procedures for this family-based study were approved by the Ethical Committees (Institutional Review Boards) of the respective University-based Centers and is currently approved by the Institutional Review Board at the University of Texas Health Sciences Center at San Antonio, TX (IRB Number HSC20170448H). The GEMM study was conducted according to the guidelines and principles expressed in the Declaration of Helsinki.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

No potential conflict of interest was reported by the authors.

References

  1. Stefan, N.; Haring, H.U.; Hu, F.B.; Schulze, M.B. Metabolically healthy obesity: Epidemiology, mechanisms, and clinical implications. Lancet Diabetes Endocrinol. 2013, 1, 152–162. [Google Scholar] [CrossRef]
  2. Virtue, S.; Vidal-Puig, A. It’s not how fat you are, it’s what you do with it that counts. PLoS Biol. 2008, 6, e237. [Google Scholar] [CrossRef] [PubMed]
  3. Lionetti, L.; Mollica, M.P.; Lombardi, A.; Cavaliere, G.; Gifuni, G.; Barletta, A. From chronic overnutrition to insulin resistance: The role of fat-storing capacity and inflammation. Nutr. Metab. Cardiovasc. Dis. 2009, 19, 146–152. [Google Scholar] [CrossRef]
  4. Stefan, N.; Schick, F.; Haring, H.U. Causes, Characteristics, and Consequences of Metabolically Unhealthy Normal Weight in Humans. Cell Metab. 2017, 26, 292–300. [Google Scholar] [CrossRef] [PubMed]
  5. Longo, M.; Zatterale, F.; Naderi, J.; Parrillo, L.; Formisano, P.; Raciti, G.A.; Beguinot, F.; Miele, C. Adipose Tissue Dysfunction as Determinant of Obesity-Associated Metabolic Complications. Int. J. Mol. Sci. 2019, 20, 2358. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Bravo, C.; Cataldo, L.R.; Galgani, J.; Parada, J.; Santos, J.L. Leptin/Adiponectin Ratios Using Either Total Or High-Molecular-Weight Adiponectin as Biomarkers of Systemic Insulin Sensitivity in Normoglycemic Women. J. Diabetes Res. 2017, 2017, 9031079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Fruhbeck, G.; Catalan, V.; Rodriguez, A.; Ramirez, B.; Becerril, S.; Salvador, J.; Colina, I.; Gomez-Ambrosi, J. Adiponectin-Leptin Ratio Is a Functional Biomarker of Adipose Tissue Inflammation. Nutrients 2019, 11, 454. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Fruhbeck, G.; Catalan, V.; Rodriguez, A.; Ramirez, B.; Becerril, S.; Salvador, J.; Portincasa, P.; Colina, I.; Gomez-Ambrosi, J. Involvement of the leptin-adiponectin axis in inflammation and oxidative stress in the metabolic syndrome. Sci. Rep. 2017, 7, 6619. [Google Scholar] [CrossRef] [PubMed]
  9. Fruhbeck, G.; Catalan, V.; Rodriguez, A.; Gomez-Ambrosi, J. Adiponectin-leptin ratio: A promising index to estimate adipose tissue dysfunction. Relation with obesity-associated cardiometabolic risk. Adipocyte 2018, 7, 57–62. [Google Scholar] [CrossRef]
  10. Aguilar-Salinas, C.A.; Velazquez Monroy, O.; Gomez-Perez, F.J.; Gonzalez Chavez, A.; Esqueda, A.L.; Molina Cuevas, V.; Rull-Rodrigo, J.A.; Tapia Conyer, R. Characteristics of patients with type 2 diabetes in Mexico: Results from a large population-based nationwide survey. Diabetes Care 2003, 26, 2021–2026. [Google Scholar] [CrossRef] [Green Version]
  11. Harris, M.I.; Flegal, K.M.; Cowie, C.C.; Eberhardt, M.S.; Goldstein, D.E.; Little, R.R.; Wiedmeyer, H.M.; Byrd-Holt, D.D. Prevalence of diabetes, impaired fasting glucose, and impaired glucose tolerance in U.S. adults. The Third National Health and Nutrition Examination Survey, 1988–1994. Diabetes Care 1998, 21, 518–524. [Google Scholar] [CrossRef] [PubMed]
  12. Morales, L.S.; Flores, Y.N.; Leng, M.; Sportiche, N.; Gallegos-Carrillo, K.; Salmeron, J. Risk factors for cardiovascular disease among Mexican-American adults in the United States and Mexico: A comparative study. Salud Publica Mex. 2014, 56, 197–205. [Google Scholar] [CrossRef] [Green Version]
  13. Barquera, S.; Campos, I.; Rivera, J.A. Mexico attempts to tackle obesity: The process, results, push backs and future challenges. Obes. Rev. 2013, 14, 69–78. [Google Scholar] [CrossRef] [Green Version]
  14. Bastarrachea, R.A.; Laviada-Molina, H.A.; Nava-Gonzalez, E.J.; Leal-Berumen, I.; Escudero-Lourdes, C.; Escalante-Araiza, F.; Peschard, V.G.; Veloz-Garza, R.A.; Haack, K.; Martinez-Hernandez, A.; et al. Deep Multi-OMICs and Multi-Tissue Characterization in a Pre- and Postprandial State in Human Volunteers: The GEMM Family Study Research Design. Genes 2018, 9, 532. [Google Scholar] [CrossRef] [Green Version]
  15. Rodriguez-Ayala, E.; Gallegos-Cabrales, E.C.; Gonzalez-Lopez, L.; Laviada-Molina, H.A.; Salinas-Osornio, R.A.; Nava-Gonzalez, E.J.; Leal-Berumen, I.; Escudero-Lourdes, C.; Escalante-Araiza, F.; Buenfil-Rello, F.A.; et al. Towards precision medicine: Defining and characterizing adipose tissue dysfunction to identify early immunometabolic risk in symptom-free adults from the GEMM family study. Adipocyte 2020, 9, 153–169. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Perkins, B.A.; Caskey, C.T.; Brar, P.; Dec, E.; Karow, D.S.; Kahn, A.M.; Hou, Y.C.; Shah, N.; Boeldt, D.; Coughlin, E.; et al. Precision medicine screening using whole-genome sequencing and advanced imaging to identify disease risk in adults. Proc. Natl. Acad. Sci. USA 2018, 115, 3686–3691. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. San-Cristobal, R.; Navas-Carretero, S.; Martínez-González, M.; Ordovas, J.M.; Martínez, J.A. Contribution of macronutrients to obesity: Implications for precision nutrition. Nat. Rev. Endocrinol. 2020, 16, 305–320. [Google Scholar] [CrossRef] [PubMed]
  18. Ma, Y.; Olendzki, B.C.; Pagoto, S.L.; Hurley, T.G.; Magner, R.P.; Ockene, I.S.; Schneider, K.L.; Merriam, P.A.; Hébert, J.R. Number of 24-hour diet recalls needed to estimate energy intake. Ann. Epidemiol. 2009, 19, 553–559. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Bergmann, S.; Keitel-Korndorfer, A.; Herfurth-Majstorovic, K.; Wendt, V.; Klein, A.M.; von Klitzing, K.; Grube, M. Recruitment strategies in a prospective longitudinal family study on parents with obesity and their toddlers. BMC Public Health 2017, 17, 145. [Google Scholar] [CrossRef] [Green Version]
  20. Han, X. Lipidomics for studying metabolism. Nat. Rev. Endocrinol. 2016, 12, 668–679. [Google Scholar] [CrossRef] [PubMed]
  21. Frost, P.A.; Chen, S.; Rodriguez-Ayala, E.; Laviada-Molina, H.A.; Vaquera, Z.; Gaytan-Saucedo, J.F.; Li, W.-H.; Haack, K.; Grayburn, P.A.; Sayers, K.; et al. Research methodology for in vivo measurements of resting energy expenditure, daily body temperature, metabolic heat and non-viral tissue-specific gene therapy in baboons. Res. Vet. Sci. 2020, 133, 136–145. [Google Scholar] [CrossRef] [PubMed]
  22. Bligh, E.G.; Dyer, W.J. A rapid method of total lipid extraction and purification. Can. J. Biochem. Physiol. 1959, 37, 911–917. [Google Scholar] [CrossRef]
  23. Wang, M.; Wang, C.; Han, R.H.; Han, X. Novel advances in shotgun lipidomics for biology and medicine. Prog. Lipid Res. 2016, 61, 83–108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Pathan, M.; Keerthikumar, S.; Ang, C.S.; Gangoda, L.; Quek, C.Y.; Williamson, N.A.; Mouradov, D.; Sieber, O.M.; Simpson, R.J.; Salim, A.; et al. FunRich: An open access standalone functional enrichment and interaction network analysis tool. Proteomics 2015, 15, 2597–2601. [Google Scholar] [CrossRef] [PubMed]
  25. Dong, Z.; Lei, X.; Kujawa, S.A.; Bolu, N.; Zhao, H.; Wang, C. Identification of core gene in obese type 2 diabetes patients using bioinformatics analysis. Adipocyte 2021, 10, 310–321. [Google Scholar] [CrossRef] [PubMed]
  26. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  27. Du, T.; Zhang, J.; Yuan, G.; Zhang, M.; Zhou, X.; Liu, Z.; Sun, X.; Yu, X. Nontraditional risk factors for cardiovascular disease and visceral adiposity index among different body size phenotypes. Nutr. Metab. Cardiovasc. Dis. 2015, 25, 100–107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Mahalle, N.P.; Garg, M.K.; Kulkarni, M.V.; Naik, S.S. Differences in traditional and non-traditional risk factors with special reference to nutritional factors in patients with coronary artery disease with or without diabetes mellitus. Indian J. Endocrinol. Metab. 2013, 17, 844–850. [Google Scholar] [CrossRef] [PubMed]
  29. Furugen, M.; Saitoh, S.; Ohnishi, H.; Akasaka, H.; Mitsumata, K.; Chiba, M.; Furukawa, T.; Miyazaki, Y.; Shimamoto, K.; Miura, T. Matsuda-DeFronzo insulin sensitivity index is a better predictor than HOMA-IR of hypertension in Japanese: The Tanno-Sobetsu study. J. Hum. Hypertens. 2012, 26, 325–333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Pietzner, M.; Kaul, A.; Henning, A.K.; Kastenmuller, G.; Artati, A.; Lerch, M.M.; Adamski, J.; Nauck, M.; Friedrich, N. Comprehensive metabolic profiling of chronic low-grade inflammation among generally healthy individuals. BMC Med. 2017, 15, 210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Iacobini, C.; Pugliese, G.; Blasetti Fantauzzi, C.; Federici, M.; Menini, S. Metabolically healthy versus metabolically unhealthy obesity. Metabolism 2019, 92, 51–60. [Google Scholar] [CrossRef]
  32. Smith, G.I.; Mittendorfer, B.; Klein, S. Metabolically healthy obesity: Facts and fantasies. J. Clin. Investig. 2019, 129, 3978–3989. [Google Scholar] [CrossRef] [Green Version]
  33. Mechanick, J.I.; Garber, A.J.; Grunberger, G.; Handelsman, Y.; Garvey, W.T. Dysglycemia-based chronic disease: An american association of clinical endocrinologists position statement. Endocr. Pract. 2018, 24, 995–1011. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Wilson, M.L. Prediabetes: Beyond the Borderline. Nurs. Clin. N. Am. 2017, 52, 665–677. [Google Scholar] [CrossRef]
  35. Misra, A.; Wasir, J.S.; Vikram, N.K. Waist circumference criteria for the diagnosis of abdominal obesity are not applicable uniformly to all populations and ethnic groups. Nutrition 2005, 21, 969–976. [Google Scholar] [CrossRef] [PubMed]
  36. Qu, H.Q.; Li, Q.; Rentfro, A.R.; Fisher-Hoch, S.P.; McCormick, J.B. The definition of insulin resistance using HOMA-IR for Americans of Mexican descent using machine learning. PLoS ONE 2011, 6, e21041. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Han, X.; Yang, K.; Gross, R.W. Multi-dimensional mass spectrometry-based shotgun lipidomics and novel strategies for lipidomic analyses. Mass Spectrom. Rev. 2012, 31, 134–178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Bagheri, M.; Farzadfar, F.; Qi, L.; Yekaninejad, M.S.; Chamari, M.; Zeleznik, O.A.; Kalantar, Z.; Ebrahimi, Z.; Sheidaie, A.; Koletzko, B.; et al. Obesity-Related Metabolomic Profiles and Discrimination of Metabolically Unhealthy Obesity. J. Proteome Res. 2018, 17, 1452–1462. [Google Scholar] [CrossRef]
  39. Caleyachetty, R.; Thomas, G.N.; Toulis, K.A.; Mohammed, N.; Gokhale, K.M.; Balachandran, K.; Nirantharakumar, K. Metabolically Healthy Obese and Incident Cardiovascular Disease Events among 3.5 Million Men and Women. J. Am. Coll. Cardiol. 2017, 70, 1429–1437. [Google Scholar] [CrossRef]
  40. Mathew, H.; Farr, O.M.; Mantzoros, C.S. Metabolic health and weight: Understanding metabolically unhealthy normal weight or metabolically healthy obese patients. Metabolism 2016, 65, 73–80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Ruderman, N.; Chisholm, D.; Pi-Sunyer, X.; Schneider, S. The metabolically obese, normal-weight individual revisited. Diabetes 1998, 47, 699–713. [Google Scholar] [CrossRef]
  42. Shea, J.L.; King, M.T.; Yi, Y.; Gulliver, W.; Sun, G. Body fat percentage is associated with cardiometabolic dysregulation in BMI-defined normal weight subjects. Nutr. Metab. Cardiovasc. Dis. 2012, 22, 741–747. [Google Scholar] [CrossRef] [PubMed]
  43. Roberson, L.L.; Aneni, E.C.; Maziak, W.; Agatston, A.; Feldman, T.; Rouseff, M.; Tran, T.; Blaha, M.J.; Santos, R.D.; Sposito, A.; et al. Beyond BMI: The “Metabolically healthy obese” phenotype & its association with clinical/subclinical cardiovascular disease and all-cause mortality—A systematic review. BMC Public Health 2014, 14, 14. [Google Scholar]
  44. Franzosi, M.G. Should we continue to use BMI as a cardiovascular risk factor? Lancet 2006, 368, 624–625. [Google Scholar] [CrossRef]
  45. Crewe, C.; An, Y.A.; Scherer, P.E. The ominous triad of adipose tissue dysfunction: Inflammation, fibrosis, and impaired angiogenesis. J. Clin. Investig. 2017, 127, 74–82. [Google Scholar] [CrossRef] [Green Version]
  46. Vishvanath, L.; Gupta, R.K. Contribution of adipogenesis to healthy adipose tissue expansion in obesity. J. Clin. Investig. 2019, 129, 4022–4031. [Google Scholar] [CrossRef] [PubMed]
  47. Despres, J.P. Body fat distribution and risk of cardiovascular disease: An update. Circulation 2012, 126, 1301–1313. [Google Scholar] [CrossRef] [Green Version]
  48. Lin, J.S.; Evans, C.V.; Johnson, E.; Redmond, N.; Coppola, E.L.; Smith, N. Nontraditional Risk Factors in Cardiovascular Disease Risk Assessment: Updated Evidence Report and Systematic Review for the US Preventive Services Task Force. JAMA 2018, 320, 281–297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Balagopal, P.B.; de Ferranti, S.D.; Cook, S.; Daniels, S.R.; Gidding, S.S.; Hayman, L.L.; McCrindle, B.W.; Mietus-Snyder, M.L.; Steinberger, J. Nontraditional risk factors and biomarkers for cardiovascular disease: Mechanistic, research, and clinical considerations for youth: A scientific statement from the American Heart Association. Circulation 2011, 123, 2749–2769. [Google Scholar] [CrossRef] [Green Version]
  50. Frühbeck, G.; Gómez-Ambrosi, J. Modulation of the leptin-induced white adipose tissue lipolysis by nitric oxide. Cell Signal 2001, 13, 827–833. [Google Scholar] [CrossRef]
  51. Asano, T.; Fujishiro, M.; Kushiyama, A.; Nakatsu, Y.; Yoneda, M.; Kamata, H.; Sakoda, H. Role of phosphatidylinositol 3-kinase activation on insulin action and its alteration in diabetic conditions. Biol. Pharm. Bull. 2007, 30, 1610–1616. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Makarova, E.; Makrecka-Kuka, M.; Vilks, K.; Volska, K.; Sevostjanovs, E.; Grinberga, S.; Zarkova-Malkova, O.; Dambrova, M.; Liepinsh, E. Decreases in Circulating Concentrations of Long-Chain Acylcarnitines and Free Fatty Acids During the Glucose Tolerance Test Represent Tissue-Specific Insulin Sensitivity. Front. Endocrinol. 2019, 10, 870. [Google Scholar] [CrossRef] [Green Version]
  53. Chung, J.O.; Koutsari, C.; Blachnio-Zabielska, A.U.; Hames, K.C.; Jensen, M.D. Effects of meal ingestion on intramyocellular ceramide concentrations and fractional de novo synthesis in humans. Am. J. Physiol. Endocrinol. Metab. 2018, 314, E105–E114. [Google Scholar] [CrossRef] [PubMed]
  54. Lambert, J.E.; Parks, E.J. Postprandial metabolism of meal triglyceride in humans. Biochim. Biophys. Acta 2012, 1821, 721–726. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. White, U.A.; Fitch, M.D.; Beyl, R.A.; Hellerstein, M.K.; Ravussin, E. Differences in In Vivo Cellular Kinetics in Abdominal and Femoral Subcutaneous Adipose Tissue in Women. Diabetes 2016, 65, 1642–1647. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Bahmad, H.F.; Daouk, R.; Azar, J.; Sapudom, J.; Teo, J.C.M.; Abou-Kheir, W.; Al-Sayegh, M. Modeling Adipogenesis: Current and Future Perspective. Cells 2020, 9, 2326. [Google Scholar] [CrossRef]
  57. Goodpaster, B.H.; Krishnaswami, S.; Harris, T.B.; Katsiaras, A.; Kritchevsky, S.B.; Simonsick, E.M.; Nevitt, M.; Holvoet, P.; Newman, A.B. Obesity, regional body fat distribution, and the metabolic syndrome in older men and women. Arch. Intern. Med. 2005, 165, 777–783. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Cao, H.; Gerhold, K.; Mayers, J.R.; Wiest, M.M.; Watkins, S.M.; Hotamisligil, G.S. Identification of a lipokine, a lipid hormone linking adipose tissue to systemic metabolism. Cell 2008, 134, 933–944. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Stefan, N.; Kantartzis, K.; Celebi, N.; Staiger, H.; Machann, J.; Schick, F.; Cegan, A.; Elcnerova, M.; Schleicher, E.; Fritsche, A.; et al. Circulating palmitoleate strongly and independently predicts insulin sensitivity in humans. Diabetes Care 2010, 33, 405–407. [Google Scholar] [CrossRef] [Green Version]
  60. Pinnick, K.E.; Neville, M.J.; Fielding, B.A.; Frayn, K.N.; Karpe, F.; Hodson, L. Gluteofemoral adipose tissue plays a major role in production of the lipokine palmitoleate in humans. Diabetes 2012, 61, 1399–1403. [Google Scholar] [CrossRef] [Green Version]
  61. Huang, K.; Zeng, Y.; Xie, Y.; Huang, L.; Wu, Y. Bioinformatics analysis of the prognostic value of CCT6A and associated signalling pathways in breast cancer. Mol. Med. Rep. 2019, 19, 4344–4352. [Google Scholar] [CrossRef] [Green Version]
  62. Rome, S.; Clément, K.; Rabasa-Lhoret, R.; Loizon, E.; Poitou, C.; Barsh, G.S.; Riou, J.P.; Laville, M.; Vidal, H. Microarray profiling of human skeletal muscle reveals that insulin regulates approximately 800 genes during a hyperinsulinemic clamp. J. Biol. Chem. 2003, 278, 18063–18068. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Shan, B.; Wang, X.; Wu, Y.; Xu, C.; Xia, Z.; Dai, J.; Shao, M.; Zhao, F.; He, S.; Yang, L.; et al. The metabolic ER stress sensor IRE1α suppresses alternative activation of macrophages and impairs energy expenditure in obesity. Nat. Immunol. 2017, 18, 519–529. [Google Scholar] [CrossRef] [PubMed]
  64. Jäger, J.; Greiner, V.; Strzoda, D.; Seibert, O.; Niopek, K.; Sijmonsma, T.P.; Schäfer, M.; Jones, A.; De Guia, R.; Martignoni, M.; et al. Hepatic transforming growth factor-β 1 stimulated clone-22 D1 controls systemic cholesterol metabolism. Mol. Metab. 2014, 3, 155–166. [Google Scholar] [CrossRef]
  65. Akoumianakis, I.; Sanna, F.; Margaritis, M.; Badi, I.; Akawi, N.; Herdman, L.; Coutinho, P.; Fagan, H.; Antonopoulos, A.S.; Oikonomou, E.K.; et al. Adipose tissue-derived WNT5A regulates vascular redox signaling in obesity via USP17/RAC1-mediated activation of NADPH oxidases. Sci. Transl. Med. 2019, 11, eaav5055. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Boström, P.; Wu, J.; Jedrychowski, M.P.; Korde, A.; Ye, L.; Lo, J.C.; Rasbach, K.A.; Boström, E.A.; Choi, J.H.; Long, J.Z.; et al. A PGC1-α-dependent myokine that drives brown-fat-like development of white fat and thermogenesis. Nature 2012, 481, 463–468. [Google Scholar] [CrossRef] [PubMed]
  67. Frühbeck, G.; Fernández-Quintana, B.; Paniagua, M.; Hernández-Pardos, A.W.; Valentí, V.; Moncada, R.; Catalán, V.; Becerril, S.; Gómez-Ambrosi, J.; Portincasa, P.; et al. FNDC4, a novel adipokine that reduces lipogenesis and promotes fat browning in human visceral adipocytes. Metabolism 2020, 108, 154261. [Google Scholar] [CrossRef] [PubMed]
  68. Leontovich, A.A.; Intine, R.V.; Sarras, M.P., Jr. Epigenetic Studies Point to DNA Replication/Repair Genes as a Basis for the Heritable Nature of Long Term Complications in Diabetes. J. Diabetes Res. 2016, 2016, 2860780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Number of cardiometabolic risk phenotypes among our symptom-free participants (N = 124) stratified by a very high ALR (9.5, N = 41), high ALR (4.0, N = 43) and low ALR (0.6, N = 40).
Figure 1. Number of cardiometabolic risk phenotypes among our symptom-free participants (N = 124) stratified by a very high ALR (9.5, N = 41), high ALR (4.0, N = 43) and low ALR (0.6, N = 40).
Biology 10 01342 g001
Figure 2. Heat maps with differences in lipidomics signatures for Time 180′ (peak of postprandial response) after a mixed meal among participants with (H) vs. (L)ALR.
Figure 2. Heat maps with differences in lipidomics signatures for Time 180′ (peak of postprandial response) after a mixed meal among participants with (H) vs. (L)ALR.
Biology 10 01342 g002
Figure 3. Volcano plot with differences in lipidomic signatures for Time 180′ (peak of postprandial response) after a mixed meal among participants with (H) vs. (L)ALR. Individual lipid classes and species are shown in boxes (right).
Figure 3. Volcano plot with differences in lipidomic signatures for Time 180′ (peak of postprandial response) after a mixed meal among participants with (H) vs. (L)ALR. Individual lipid classes and species are shown in boxes (right).
Biology 10 01342 g003
Figure 4. Heat maps with differences in miRNA signatures for Time 180′ (peak of postprandial response) after a mixed meal among participants with (H) vs. (L)ALR.
Figure 4. Heat maps with differences in miRNA signatures for Time 180′ (peak of postprandial response) after a mixed meal among participants with (H) vs. (L)ALR.
Biology 10 01342 g004
Figure 5. Volcano plot with differences in miRNA signatures for Time 180′ (peak of postprandial response) after a mixed meal among participants with (H) vs. (L)ALR.
Figure 5. Volcano plot with differences in miRNA signatures for Time 180′ (peak of postprandial response) after a mixed meal among participants with (H) vs. (L)ALR.
Biology 10 01342 g005
Figure 6. Twenty logarithmic fold change (FC) > 2 and significant p-values for skeletal muscle (SKM) differentially expressed genes (DEGs). Volcano plot of FC differences of SKM RNA expression between RG and CG. The vertical lines correspond to FC > 2.
Figure 6. Twenty logarithmic fold change (FC) > 2 and significant p-values for skeletal muscle (SKM) differentially expressed genes (DEGs). Volcano plot of FC differences of SKM RNA expression between RG and CG. The vertical lines correspond to FC > 2.
Biology 10 01342 g006
Figure 7. Heat map with differences in gene expression from SKM signatures from the 20 genes identified as shown in Figure 6. Gene CCT6A shows “hot red” values perhaps indicating higher activity in several fetal and adult pathways and, interestingly, from the immune systems through B cells, CDK4, CDK8, monocytes and NK cells.
Figure 7. Heat map with differences in gene expression from SKM signatures from the 20 genes identified as shown in Figure 6. Gene CCT6A shows “hot red” values perhaps indicating higher activity in several fetal and adult pathways and, interestingly, from the immune systems through B cells, CDK4, CDK8, monocytes and NK cells.
Biology 10 01342 g007
Figure 8. Graphical representation of enrichment analysis and interaction networks obtained through the FunRich database analysis tool for gene annotation pertaining to transcript expression associations in SKM tissue.
Figure 8. Graphical representation of enrichment analysis and interaction networks obtained through the FunRich database analysis tool for gene annotation pertaining to transcript expression associations in SKM tissue.
Biology 10 01342 g008
Figure 9. Twenty logarithmic fold change (FC) > 2 and significant p-values for adipose tissue (AT) differentially expressed genes (DEGs). Volcano plot of FC differences of AT RNA expression between RG and CG. The vertical lines correspond to FC > 2.
Figure 9. Twenty logarithmic fold change (FC) > 2 and significant p-values for adipose tissue (AT) differentially expressed genes (DEGs). Volcano plot of FC differences of AT RNA expression between RG and CG. The vertical lines correspond to FC > 2.
Biology 10 01342 g009
Figure 10. Heat map with differences in gene expression from adipose tissue signatures from the 20 genes identified as shown in Figure 9. Gene MCM4 shows “hot red” values perhaps indicating higher activity in several fetal pathways and, interestingly, from the immune systems through B cells, CDK4, CDK8, monocytes and NK cells.
Figure 10. Heat map with differences in gene expression from adipose tissue signatures from the 20 genes identified as shown in Figure 9. Gene MCM4 shows “hot red” values perhaps indicating higher activity in several fetal pathways and, interestingly, from the immune systems through B cells, CDK4, CDK8, monocytes and NK cells.
Biology 10 01342 g010
Figure 11. Graphical representation of enrichment analysis and interaction networks obtained through the FunRich database analysis tool for gene annotation pertaining to transcript expression associations in adipose tissue.
Figure 11. Graphical representation of enrichment analysis and interaction networks obtained through the FunRich database analysis tool for gene annotation pertaining to transcript expression associations in adipose tissue.
Biology 10 01342 g011
Table 1. Demographic characteristics and clinical parameters for the adipose tissue dysfunction phenotype in our symptom-free GEMM volunteers.
Table 1. Demographic characteristics and clinical parameters for the adipose tissue dysfunction phenotype in our symptom-free GEMM volunteers.
Demographic CharacteristicsPhenotypes and Cut-Offs (Mean ± SD)
N = 124Adipose Tissue (Dys)function
Sympton-Free Participants
(Metabolic Risk Criteria)
ALR > 1 (N = 84)ALR < 1 (N = 40)p
ALR6.7 ± 10.00.6 ± 0.3<0.0001
Age37 ± 1439 ± 140.3887
Weight (kg)70 ± 1584 ± 16<0.0001
Waist Circunference (cm)87 ± 1299 ± 15<0.0001
BMI (kg/m2)27 ± 532 ± 6<0.0001
% Fat Total34 ± 940 ± 8<0.0001
Systolic Pressure (mmHg) ≥ 130111 ± 11111 ± 140.7483
Diastolic Pressure (mmHg) ≥ 8570 ± 1072 ± 100.3103
Fasting Glucose (mg/dL) ≥ 10097 ± 22110 ± 290.0419
2-h Glucose (mg/dL) ≥ 140127 ± 42145 ± 540.0935
Triglycerides (mg/dL) ≥ 150118 ± 51150 ± 570.0007
HDL-cholesterol (mg/dL) < 40, Men; <50 Women48 ± 1442 ± 110.0313
High sensitive C-reactive protein (mg/L) (>90th percentile) *10.0 ± 13.819.6 ± 22.20.0027
Whole body insulin resistance (Matsuda index) **4.6 ± 4.02.7 ± 2.9<0.0001
Absence of diagnosis or therapy of cardiometabolic diseases: T2D, hypertension, dyslipidemia, NAFLD, CKD, or CVD; or treatment with blood pressure, lipid, or diabetes medications.* hsCRP 90% percentile: >35.72 mg/L; ** Whole body insulin resistant: equal or lower than 2.5.
Table 2. Immunometabolic and postprandial parameters for the adipose tissue dysfunction phenotype in our symptom-free GEMM volunteers.
Table 2. Immunometabolic and postprandial parameters for the adipose tissue dysfunction phenotype in our symptom-free GEMM volunteers.
Immunometabolic and Postprandial PhenotypesPhenotypes (Mean ± SD)
N = 124Adipose Tissue (Dys)function
Sympton-Free ParticipantsALR > 1 (N = 84)ALR < 1 (N = 40)p
Adiponectin (μg/mL)32.7 ± 33.010.1 ± 5.2<0.0001
Leptin (ng/mL)8.7 ± 6.521.0 ± 18.4<0.0001
TNFa (pg/mL)5.0 ± 2.95.2 ± 2.80.6290
IL-6 (pg/mL)9.1 ± 22.68.6 ± 9.50.0014
MCP-1 pg/mL108.2 ± 143.9107.1 ± 41.80.0544
PAI-1 (pg/mL) 34,169.0 ± 29,375.539,182.1 ± 30,956.30.2658
Fibrinogen (mg/dL)107.7 ± 40.8131.4 ± 73.40.0355
Fasting Insulin (microU/mL)14.6 ± 18.625.2 ± 20.3<0.0001
C-peptide (pg/mL)1.2 ± 0.51.9 ± 0.9<0.0001
Insulin 120′ [microU/mL]64.8 ± 44.4120.9 ± 74.4<0.0001
HOMA-IR3.6 ± 4.77.2 ± 6.5<0.0001
Glucagon (pg)mL)47.4 ± 64.246.5 ± 45.30.5168
GLP-1 (pg/mL)13.0 ± 36.317.3 ± 44.30.8721
Leptin AUC (5h)259861940.0002
Insulin AUC (5h)16,84728,1580.0821
Glucose AUC (5h)36,48640,8830.1124
Triglycerides AUC 5h48,22059,7980.0025
GLP-1 AUC 5h474353460.1306
Glucagon AUC 5h18,48717,4350.1988
C-peptide AUC 5h916.412460.0821
Table 3. Prevalence and percentages of parameters and cut-offs widely used to define metabolic risk and a healthy/unhealthy cardiometabolic profile (MH/MUH) in subjects with adipose tissue dysfunction determined by a mean (H) or (L)ALR.
Table 3. Prevalence and percentages of parameters and cut-offs widely used to define metabolic risk and a healthy/unhealthy cardiometabolic profile (MH/MUH) in subjects with adipose tissue dysfunction determined by a mean (H) or (L)ALR.
Prevalence and Percentage of Individuals with Risk Phenotypes
Group 1 (N = 41)Group 2 (N = 43)Group 3 (N = 40)
MH/MUH risk criteria and cut-offsMean ALR 9.5 ± 13.0Mean ALR 4.0 ± 4.5Mean ALR 0.6 ± 0.3
Diabetic A1c > 6.50 (0.0%)3 (7.0%)4 (10.0%)
Prediabetic A1c 5.7–6.4 1 (2.4%)6 (14.0%)7 (17.5%)
Matsuda Index < 2.54 (9.8%)17 (39.5%)25 (62.5%)
HOMA-IR > 2.612 (29.3)16 (37.2%)28 (70.0%)
hsCRP > 35.70 (0.0%)4 (9.3%)8 (20.0%)
Glucose > 1000 (0.0%)28 (65.1%)16 (40.0%)
Triglycerides > 1503 (7.3%)17 (39.5%)16 (40.0%)
HDL < 40 Men < 50 Women16 (39.0%)36 (83.7%)33 (82.5%)
Dias BP > 851 (2.4%)5 (11.6%)4 (10.0%)
Sys BP > 1300 (0.0%)5 (11.6%)4 (10.0%)
Waist > 88 Women - >102 Men3 (7.3%)26 (60.5%)30 (75.0%)
Table 4. Demographics and clinical biochemistries in symptom-free adult volunteers (N = 14).
Table 4. Demographics and clinical biochemistries in symptom-free adult volunteers (N = 14).
Demographic CharacteristicsMean Values for 14 Females
(H) ALR (N = 9)SD (±)(L) ALR (N = 5)SD (±)
Adipo/Lep Ratio2.21.10.50.4
Age (Yr)38.511.832.613.0
% Fat Total42.76.246.03.2
Weight (kg)62.15.777.617.8
Waist Circumference (cm)83.99.993.916.7
BMI (kg/m2)26.83.132.87.8
Triglycerides (mg/dL)119.842.2145.060.1
HDL-Cholesterol (mg/dL)41.49.647.89.4
Adiponectin (μg/mL)24.915.57.02.8
Leptin (ng/mL)11.13.921.011.3
Fasting Glucose (mg/dL)87.06.986.28.4
Glucose 120′ (mg/dL)121.221.8107.07.0
Fasting Insulin (microU/mL)7.14.217.39.5
Insulin 120′ (microU/mL)59.542.4101.742.2
Matsuda Index6.64.12.91.3
HOMA-IR1.50.93.72.0
PAI-1 (pg/mL)42,867.531,186.856,408.457,940.0
IL-6 (pg/mL)1.61.23.12.5
TNFa (pg/mL)2.51.53.42.2
MCP-1 pg/mL109.838.3133.861.2
Table 5. Demographics and ALR measurements of each participant (N = 6).
Table 5. Demographics and ALR measurements of each participant (N = 6).
Female ID (N = 6)GenderAge % FatAdiponectin
(μg/mL)
Leptin
(ng/mL)
Adipo/Lep Ratio
MTY0003F4934.518.377.412.48
MTY0007F3836.212.135.422.24
MTY0006F3350.323.7311.502.06
MTY0014F4549.916.9615.031.13
Mean 4142.717.809.841.98
SD (±) 78.54.774.290.59
MTY0009F2043.28.468.930.95
MTY0010F35473.1235.270.09
Mean 2845.15.7922.100.52
SD (±) 102.73.7818.620.61
Table 6. Contrast linear combination of variables between the CG and RG (N = 6).
Table 6. Contrast linear combination of variables between the CG and RG (N = 6).
(H)ALR > 1 (L)ALR < 1
VariableMeanSD (±)SEMeanSD (±)SEDiff.SE Diff.p-Value
Waist Circumference (cm)82.57.54.339018.0810.44−7.511.30.543
BMI (kg/m2)25.961.8921.09232.48.4124.856−6.434.9780.265
% total fat40.338.6735.00746.73.361.939−6.365.370.301
Fat Mass kg23.87.94.56132.8810.826.248−9.087.7360.305
Muscle Mass kg34.281.7180.99137.1811.056.384−2.96.460.676
Triglycerides (mg/dL)125.39.0185.206155.686.4349.9−30.350.170.578
Creatinin (mg/dL)0.60.2640.1520.6660.1520.088−0.060.1760.724
Uric acid (mg/dL)3.7331.7951.0365.9330.650.375−2.21.1020.116
BUN (mg/dL)842.3099.3330.5770.333−1.332.3330.598
Total Cholesterol (mg/dL)162.380.2246.3115611.136.4296.33346.760.898
HDL (mg/dL)3514.798.54447.6613.277.666−12.611.470.331
LDL (mg/dL)11059.0234.07799.5395.5073134.520.419
VLDL (mg/dL)17.336.4293.71129.33179.82−1210.490.316
Alt (U/L)16.665.8593.38228.3312.747.356−11.68.0960.223
Ast (U/L)33.3312.897.44655.3318.8210.86−2213.170.17
Alk phos (U/L)73.3324.5414.1649.6632.1218.5523.6623.340.367
Adipo/lep Ratio2.8041.4480.8360.6270.2130.1232.1770.8450.061 *
PAI-1 (pg/mL)21,09528,05716,19868,89646,81027,026−478031,5090.203
MCP-1 (pg/mL)98.320.5511.86124.733.4119.28−26.422.640.307
IL-6 (pg/mL)0.6560.7050.4072.8330.6670.385−2.170.560.017 **
TNF-a (pg/mL)2.1691.470.8493.0351.1150.644−0.861.0650.461
hsCRP (mg/L)0.0580.0290.0170.3040.1650.095−0.240.0960.063 *
Matsuda Index8.9915.4053.123.3131.6670.9625.6783.2650.157
HOMA-IR0.9050.3930.2263.392.881.663−2.481.6780.212
AUC Glucose 35061.5535.53310.92.1551.24439.0835.560.333
AUC Insulin134.861.8635.71282.887.2250.36−14861.740.074 *
AUC GLP-1304.852.3230.21470.3236.1136.3−165139.60.301
Statistically Significant p-Value: **; Marginally Significant p-value: *.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gallegos-Cabriales, E.C.; Rodriguez-Ayala, E.; Laviada-Molina, H.A.; Nava-Gonzalez, E.J.; Salinas-Osornio, R.A.; Orozco, L.; Leal-Berumen, I.; Castillo-Pineda, J.C.; Gonzalez-Lopez, L.; Escudero-Lourdes, C.; et al. Replication of Integrative Data Analysis for Adipose Tissue Dysfunction, Low-Grade Inflammation, Postprandial Responses and OMICs Signatures in Symptom-Free Adults. Biology 2021, 10, 1342. https://doi.org/10.3390/biology10121342

AMA Style

Gallegos-Cabriales EC, Rodriguez-Ayala E, Laviada-Molina HA, Nava-Gonzalez EJ, Salinas-Osornio RA, Orozco L, Leal-Berumen I, Castillo-Pineda JC, Gonzalez-Lopez L, Escudero-Lourdes C, et al. Replication of Integrative Data Analysis for Adipose Tissue Dysfunction, Low-Grade Inflammation, Postprandial Responses and OMICs Signatures in Symptom-Free Adults. Biology. 2021; 10(12):1342. https://doi.org/10.3390/biology10121342

Chicago/Turabian Style

Gallegos-Cabriales, Esther C., Ernesto Rodriguez-Ayala, Hugo A. Laviada-Molina, Edna J. Nava-Gonzalez, Rocío A. Salinas-Osornio, Lorena Orozco, Irene Leal-Berumen, Juan Carlos Castillo-Pineda, Laura Gonzalez-Lopez, Claudia Escudero-Lourdes, and et al. 2021. "Replication of Integrative Data Analysis for Adipose Tissue Dysfunction, Low-Grade Inflammation, Postprandial Responses and OMICs Signatures in Symptom-Free Adults" Biology 10, no. 12: 1342. https://doi.org/10.3390/biology10121342

APA Style

Gallegos-Cabriales, E. C., Rodriguez-Ayala, E., Laviada-Molina, H. A., Nava-Gonzalez, E. J., Salinas-Osornio, R. A., Orozco, L., Leal-Berumen, I., Castillo-Pineda, J. C., Gonzalez-Lopez, L., Escudero-Lourdes, C., Cornejo-Barrera, J., Escalante-Araiza, F., Huerta-Avila, E. E., Buenfil-Rello, F. A., Peschard, V. -G., Silva, E., Veloz-Garza, R. A., Martinez-Hernandez, A., Barajas-Olmos, F. M., ... Bastarrachea, R. A. (2021). Replication of Integrative Data Analysis for Adipose Tissue Dysfunction, Low-Grade Inflammation, Postprandial Responses and OMICs Signatures in Symptom-Free Adults. Biology, 10(12), 1342. https://doi.org/10.3390/biology10121342

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop