Next Article in Journal
Magnetic-Based Human Tissue 3D Cell Culture: A Systematic Review
Next Article in Special Issue
Circulating Adipokines and Hepatokines Serve as Diagnostic Markers during Obesity Therapy
Previous Article in Journal
Isomelezitose Overproduction by Alginate-Entrapped Recombinant E. coli Cells and In Vitro Evaluation of Its Potential Prebiotic Effect
Previous Article in Special Issue
Molecular and Physiological Effects of Browning Agents on White Adipocytes from Bone Marrow Mesenchymal Stromal Cells
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 20
10.3390/ijms232012680
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Macrophage Phenotypes and Gene Expression Patterns Are Unique in Naturally Occurring Metabolically Healthy Obesity

by
Alistaire D. Ruggiero
1,
Ravichandra Vemuri
1,
Masha Block
1,
Darla DeStephanis
1,
Matthew Davis
2,
Jeff Chou
3,
Abigail Williams
1,
Ashlynn Brock
1,
Swapan Kumar Das
4 and
Kylie Kavanagh
1,5,*
1
Department of Pathology, Wake Forest University School of Medicine, Winston-Salem, NC 27157, USA
2
Department of Internal Medicine, Wake Forest University School of Medicine, Winston-Salem, NC 27157, USA
3
Department of Biostatistics and Data Science, Wake Forest University School of Medicine, Winston-Salem, NC 27157, USA
4
Department of Endocrinology and Metabolism, Wake Forest University School of Medicine, Winston-Salem, NC 27157, USA
5
College of Health and Medicine, University of Tasmania, Hobart 7000, Australia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(20), 12680; https://doi.org/10.3390/ijms232012680
Submission received: 26 August 2022 / Revised: 14 October 2022 / Accepted: 18 October 2022 / Published: 21 October 2022
(This article belongs to the Special Issue New Insights into Adipose Tissue Metabolic Function and Dysfunction)
Browse Figures
Review Reports Versions Notes

Abstract

:
Obesity impacts 650 million individuals globally, often co-occurring with metabolic syndrome. Though many obese individuals experience metabolic abnormalities (metabolically unhealthy obese [MUO]), ~30% do not (metabolically healthy obese [MHO]). Conversely, >10% of lean individuals are metabolically unhealthy (MUL). To evaluate the physiologic drivers of these phenotypes, a 44-animal African green monkey cohort was selected using metabolic syndrome risk criteria to represent these four clinically defined health groups. Body composition imaging and subcutaneous adipose tissue (SQ AT) biopsies were collected. Differences in adipocyte size, macrophage subtype distribution, gene expression, vascularity and fibrosis were analyzed using digital immunohistopathology, unbiased RNA-seq, endothelial CD31, and Masson’s trichrome staining, respectively. MHO AT demonstrated significant increases in M2 macrophages (p = 0.02) and upregulation of fatty acid oxidation-related terms and transcripts, including FABP7 (p = 0.01). MUO AT demonstrated downregulation of these factors and co-occurring upregulation of immune responses. These changes occurred without differences in AT distributions, adipocyte size, AT endothelial cells, collagen I deposition, or circulating cytokine levels. Without unhealthy diet consumption, healthy obesity is defined by an increased SQ AT M2/M1 macrophage ratio and lipid handling gene expression. We highlight M2 macrophages and fatty acid oxidation as targets for improving metabolic health with obesity.

1. Introduction

More than 1.9 billion adults worldwide are overweight (BMI: 25–29.9 kg/m2), and over 650 million are obese (BMI ≥ 30 kg/m2) [1]. Obesity often coincides with metabolic syndrome (MetS), which is defined as the presence of two or more of the following: increased waist circumference, hypertriglyceridemia, low high-density lipoprotein (HDL) cholesterol, elevated fasting glucose, and/or hypertension [2]. However, roughly 30% of obese people remain free of MetS [3], while over 10% of lean individuals demonstrate metabolic abnormalities [2]. The mechanisms that sustain the metabolically healthy obese (MHO) or propagate the metabolically unhealthy lean (MUL) phenotypes are undetermined.
The MHO and MUL phenotypes are evidence that complex interactions exist between adipose tissue (AT) and metabolic regulation. Macrophages constitute 5% of white AT cell content in healthy lean states and up to 50% of cell content in obese states [4]. As the predominant leukocyte in AT, macrophages respond to stimuli in their microenvironments and range from highly pro-inflammatory to highly anti-inflammatory [5]. While heterogeneous, AT macrophages are often defined by the extremes of their inflammatory states [5]. Recruited M1, or classically activated, macrophages are prompted by Toll-like receptor (TLR) ligands [6,7], and demonstrate high antigen presentation and expression of pro-inflammatory cytokines [7]. Embryo-derived M2 macrophages, or alternatively activated macrophages, are induced by anti-inflammatory stimuli [7]. Intermediate macrophages express M1 and M2 markers, and both associate with insulin resistance and demonstrate increases in mitochondria and upregulation of fatty acid metabolism enzymes [8]. Excess M1 macrophage infiltration results in increased tissue inflammation, while an overabundance of M2 macrophages can cause aberrant fibrosis [9]. AT macrophages display variable metabolic characteristics in situ based on their phenotypes. Adipose M2 macrophages may interact beneficially with adipocytes due their ability to release catecholamines and induce local lipolysis, and accept intercellular transfer of mitochondria from neighboring adipocytes [10,11,12]. They are also better suited for fueling oxidative phosphorylation using locally produced fatty acids [10,11,12]. To date, the spectrum of macrophages present in AT has not been determined in all four metabolic health groupings. Further, macrophage types are most commonly studied as part of the high-fat and high-sugar diet-driven MUO phenotype resulting in M1 macrophages and inflammation dominating the literature [13].
The purpose of this work was to identify underlying AT differences in distributions, macrophage profiles and gene expression in these four metabolic health groups. Here, we evaluated SQ AT differences between metabolic phenotypes (metabolically healthy lean [MHL] vs. MHO, MUL, and metabolically unhealthy obese [MUO]) that occur spontaneously without unhealthy diet exposure or changes in the environment. To do this, we used a unique nonhuman primate (NHP) model, the African green monkey (Chlorocebus aethiops sabaeus). We employed standard immunohistochemical approaches to evaluate macrophage subtypes and unbiased RNA-seq to evaluate subcutaneous AT transcriptional differences between metabolic health phenotypes. In doing so, we determined that MHO SQ AT is uniquely different in its macrophage cell polarization and metabolic gene expression signatures.

2. Results

2.1. Abdominal Adipose Tissue Distribution by Depot, Adipocyte Size, and Circulating Insulin Do Not Explain Poor Metabolic Health

By design, obese animals (MHO and MUO) were selected to have increased waist circumferences (Table 1). On average, obese animals weighed 1.6 kilograms more, and their average body fat was 19 percent higher than lean animals (Table 1). One point was allotted for each MetS risk criteria (cut-off values defined in Supplementary Table S1). In line with human clinical categorizations, MHL animals had a MetS score of zero, MHO animals had a MetS score of one for their obesity status, and MUL and MUO animals had MetS scores of two or higher (Figure 1a) [14]. By design, the MHO and MUO animals demonstrated increased waist circumferences compared to the MHL and MUL groups. MUL and MUO animals demonstrated increased fasting glucose values, percent glycosylated hemoglobin), systolic and diastolic blood pressure, and triglycerides. Computed tomography (CT) data concluded that obese animals exhibited increased percent fat mass and decreased percent lean mass. However, CT image analysis revealed no differences in fat-to-lean mass ratio, or abdominal SQ or visceral (VIS) adipose volume by health status (Figure 1b–d). Obese animals demonstrated lower VIS:SQ ratios than lean animals, though values between MHO and MUO animals were not different (Table 1). Lean animals demonstrated significantly less SQ and VIS adipose. The substantially lower amounts of SQ adipose drove the high VIS:SQ ratios in lean animals. Obese animals also had decreased percentage of lean mass body weight (Table 1). While food intake and caloric expenditure were not measured, percent lean mass, a measure of primarily muscle, was comparable between the MHO and MUO groups. Additionally, no differences in liver adipose, estimated as liver attenuation values, were identified (Table 1). MHL and MHO nonhuman primate (NHP) circulating insulin concentrations were comparable (Table 1) in contrast to MUO animals, which demonstrated significantly increased circulating insulin concentrations. MUL animals demonstrated significantly increased systolic and diastolic blood pressures. Histological evaluation of SQ adipose illustrated that differences in adipocyte cell size was driven by obesity status alone, where obese animals displayed larger SQ adipocytes regardless of their health status (Figure 1e; Supplementary Figures S1 and S2). Adipocyte size did not associate with metabolic syndrome-related clinical parameters.

2.2. Metabolically Healthy Obese Subcutaneous Adipose Demonstrates Increased Anti-Inflammatory M2 Macrophage Accumulation

We assessed tissue inflammatory profiles by characterizing SQ AT macrophage subtypes. Supplementary Figure S2 illustrates the staining combinations used to identify undefined, M1, intermediate, and M2 macrophages. SQ total AT macrophage density was comparable between all groups (Supplementary Table S2); however, MHO animals demonstrated significantly higher M2/M1 macrophage ratios (p = 0.02; Figure 2a–c). Differences in specific macrophage subtypes did not reach statistical significance (Supplementary Table S2), highlighting the importance of the M2/M1 macrophage ratio in local and systemic inflammation. Anti-inflammatory M2 macrophages were negatively associated with circulating levels of monocyte chemoattractant protein (MCP)-1 (r = −0.34, p = 0.03; Supplementary Figure S3). Though this association is moderate, it indicates that the healthy SQ AT milieu likely initiates less macrophage recruitment for harmful inflammation resolution. MUO animals demonstrated the highest percentage of both pro-inflammatory M1 and intermediate macrophages (Figure 2b; Supplementary Table S1). Notably, a large percentage of SQ adipose macrophages in each group were undefined or not yet polarized (Figure 2b).

2.3. Metabolically Healthy Obese Subcutaneous Adipose Demonstrates Upregulation of Lipid Handling-Related Transcripts without Immune Activation

We performed bulk RNA-seq on SQ AT biopsies to identify transcript expression differences between groups. In total, 15,423 transcripts were identified after the removal of non-protein coding genes and genes ambiguously mapped to chromosomes. The largest number of genes that were differentially expressed was between MHO and MHL AT (n = 1092), whereas only 56 genes were differentially expressed between the MHO and MUO groups, 154 genes were differentially expressed between the MHL and MUL groups, and 347 were differentially expressed between the MUL and MUO groups. MHO animals demonstrated upregulation of transcripts related to lipid synthesis, lipid dynamics and lipid transport and downregulation of transcripts related to inflammation initiation (Figure 3a,b), while MUO animals demonstrated the inverse pattern (Table 2 and Table 3; Supplementary Figure S5). Gene ontology term analysis revealed gene pathways specific to healthy obesity. Both the “acyl-CoA dehydrogenase activity” and “flavin adenine dinucleotide binding” pathways were upregulated in the MHO compared to the MHL, indicating that MHO animals upregulate transcripts related to fatty acid oxidation (Figure 3a). MHO animals also illustrated enrichment of the “metabolic process” pathway. Importantly, no immune response-related pathways were found to be differently enriched between the MHO and MHL groups. On the contrary, MUO animals illustrated downregulation of the “fatty acid metabolic process”, “acyl-CoA biosynthetic process”, and “lipid metabolic process” pathways compared to MHL (Figure 3b). MUO animals demonstrated increases in “immune response” and “cell chemotaxis” pathways compared to MHL, a signature not seen in MHO SQ adipose (Figure 3b). The “immune response” gene ontology term is comprised of 1973 genes, including those related to the innate immune response, B and T cell receptor signaling, and macrophage cytokine production. The “cell chemotaxis” gene ontology term is comprised on 316 genes, including those related to lymphocyte and T cell chemotaxis. Specific genes in these terms that were upregulated in the MUO compared to the MHL SQ AT included C-X-C Motif Chemokine Ligand 16 (CXCL16), C-C Motif Chemokine Ligand 11 (CCL11), CCL22, and CCL26. Pathway analysis between MUL and MUO groups identified enrichment of immune regulating pathways with combined obesity and metabolic disease (Supplementary Figure S4). MUO AT demonstrated significant downregulation of the “cellular respiration” and “tricarboxylic acid cycle” gene ontology pathways compared to MHL. This evidence of decreased mitochondrial function in the MUO but not in MHO aligns with our previous study where we identified poor mitochondrial quality control markers in MUO compared to MHO NHPs [15]. Gene term enrichment identified from comparisons of MHO and MUO, and MHL and MUL animals did not survive corrections. While they were ultimately not statistically significant, the terms “extracellular region”, “serine-type endopeptidase activity” and “transmembrane activity” were enriched in the MHO compared to the MUO (Supplementary Figure S5).
Transcripts upregulated in MHO included fatty acid binding protein 7 (FABP7), monoacylglycerol O-acyltransferase 1 (MOGAT1), ATPase phospholipid transporting 8B3 (ATP8B3), glycerol-3-phosphate acyltransferase 3 (GPAT3), and adrenoceptor beta 3 (ADRB3). Upregulation of FABP7 indicates that MHO animals have improved signaling for fatty acid movement across membranes (Table 2). Downregulation of C-C motif chemokine receptor 1 (CCR1), pancreatic lipase related protein 1 (PNLIP), macrophage scavenger receptor 1 (MSR1), and oxidized low-density lipoprotein receptor 1 (OLR1) was observed. Activation of MSR1 has been related to a phenotypic macrophage switch from anti-inflammatory to pro-inflammatory [16]. Therefore, downregulation of MSR1 in MHO adipose indicates maintenance of an anti-inflammatory M2-based macrophage milieu, which coincides with the macrophage data presented above. The MHO animals’ lower expression of CCR1, PNLIP, and OLR1 also suggests decreased tissue inflammation, atherogenic capacity and lectin activation [17,18]. The sustained anti-inflammatory environment in the MHO adipose occurred in the absence of differences in adipocyte size when compared to the MUO (Figure 1e) and radiographic adipose or SQ AT cell density (Table 1) differences.
In contrast, MUO AT demonstrated upregulated expression of transcripts related to inflammatory processes and glycogen storage (Table 3), including plasminogen (PLG) and glycogen synthase 2 (GYS2). Heats maps of differentially expressed genes between the two obese groups and the MHL group are shown in Supplementary Figure S6. Though inverse transcript expression was observed between MHO and MUO groups (Figure 3c), no differentially expressed transcripts between these groups met corrections for multiple comparisons. While they did not meet statistical inquiry, transcripts associated with inflammation, including CCL19 and CCL22, were downregulated in MHO AT compared to MUO, corroborating previous findings [19]. Transcripts associated with nutrient transport, including solute carrier family (SLC) 22A3 and SLC10A6, were upregulated in MHO AT compared to MUO (Supplementary Table S3) supporting an expression pattern consistent with low inflammation and high lipogenic capacity.
In line with these transcript expression differences between groups, the SQ AT M2/M1 ratio positively correlated with expression levels of ATP binding cassette subfamily C member 2 (ABCC2) (r = 0.315, p = 0.0448), and negatively correlated with expression levels of cyclin dependent kinase inhibitor 1A (CDKN1A) (r = −0.419, p = 0.00646), immunoglobulin heavy constant delta (IGHD) (r = −0.398, p = 0.00988), and heat shock transcription factor 4 (HSF4) (r = −0.436, p = 0.00436). These associations indicate that increases in the M2/M1 ratio correspond with decreased antigen binding, cellular senescence, and stress response, and improved movement of molecules across membranes (Supplementary Table S4), and support the MHO-specific downregulation of inflammatory transcripts and upregulation of lipid handling-related transcripts noted above.

2.4. Tissue Cytokine Burdens Match Tissue Macrophage Profiles

We found no group differences in circulating adiponectin, MCP-1, IL-6, IL-1β or PAI-1 (Supplementary Table S5), but a number of biologically relevant relationships were present. Circulating levels of PAI-1 demonstrated a trend toward a positive association with the increasing percentages of M1 macrophages (r = 0.27, p = 0.09). Circulating levels of PAI-1 also were negatively correlated with percentages of M2 macrophages (r = −0.29, p = 0.06), indicating that pro-inflammatory PAI-1 levels decrease as anti-inflammatory macrophage densities increase. Tissue levels of MCP-1 differed by obesity status (p = 0.02), where obese animals displayed higher levels. Tissue MCP-1 also demonstrated a trend toward significance by health status (p = 0.06), where unhealthy animals had higher levels (Supplementary Table S5). SQ AT levels of anti-inflammatory IL-10 were significantly different between groups (p = 0.03), where the MHL group had the highest detectable tissue levels (Supplementary Table S4). Tissue levels of IL-10 also differed significantly by obesity status (p = 0.004), where obese animals demonstrated significantly lower levels (Supplementary Table S5). Tissue levels of TGF-β exhibited (1) a trend toward significance (p = 0.09), where the MUL group displayed increased tissue levels, and (2) a trend toward a significant difference by health, where unhealthy animals illustrated higher tissue levels (p = 0.08). Additionally, abdominal SQ AT concentrations of TGF-β negatively correlated with the SQ AT M2/M1 macrophage ratio (r = −0.351, p = 0.045). No other cytokines measured correlated with the M2/M1 macrophage ratio. Accordingly, our tissue cytokines tracked with our observed macrophage phenotypic shifts.

2.5. Changes in Tissue Vasculature or Fibrosis Do Not Explain Poor Health but Relate to Adipocyte Hypertrophy

Hematoxylin and eosin staining was used to evaluate adipocyte size differences, while CD31 staining and Masson’s Trichrome staining (MTC) were used to quantitate tissue vasculature and fibrosis, respectively (Supplementary Figure S7). No differences in either CD31-positive area or percent MTC-positive staining were present between groups (Supplementary Table S6). CD31-positive area correlated positively with adipocyte size (r = 0.33, p = 0.03; Figure 4a), indicating that as adipocyte size increased, the vasculature increased, possibly to meet adipocyte nutrient requirements. Percent MTC-positive staining correlated negatively with increasing SQ adipocyte size (r = −0.34, p = 0.03; Figure 4b), highlighting that as adipocyte size increased, interstitial collagen I formation decreased. Additionally, no expression differences in other ECM transcripts, including collagen types V, VI, VII and IX, were observed in MHO AT. Collagens can be upregulated in aberrant obesity [20]. Our MUO SQ AT demonstrated upregulation of COL4A3 and downregulation of MMP7 compared to the MHL indicating deposition of fibers with decreased breakdown (Table 3). No other ECM transcripts, including other matrix metalloproteinases, laminins, or fibronectins, were differentially expressed between groups. Our results indicate that neither tissue fibrosis nor vasculature changes explain health status in lean or obese states. Rather, we posit that metabolic dysfunction in the presence and absence of obesity may be caused by other factors, such as the differences in immune cell populations and adipocyte function shifts observed here.

3. Discussion

In this study, we identify features of naturally occurring obesity in primates, with the unique determination that upregulation of lipid handling-related transcripts and an increased M2/M1 macrophage ratio are associated with maintained health despite obesity. Unhealthy obesity was characterized by an opposing profile with a decrease in the SQ AT M2/M1 ratio, tissue immune activation and metabolic process downregulation. Previous characterization of metabolically healthy obese individuals determined that maintained health with obesity corresponded with upregulated mitochondrial oxidative phosphorylation and fatty acid oxidation transcript expression, decreased levels of markers of chronic inflammation, and reduced liver adipose accumulation [21]. Additionally, expression of SLC39A8 and SLC16A7, solute carrier family members, in circulation or saliva were associated with maintained metabolic health in obese women [22,23]. Our work verifies these findings and adds vital information regarding differences in local immune cell populations in healthy and unhealthy obesity. Supporting our SQ AT immune cell findings, reductions in the M2/M1 AT macrophage ratio previously correlated with increasing insulin resistance [24,25]. In our study, the M2/M1 ratio negatively correlated with expression levels of transcripts related to antigen binding, stress response, and cellular senescence, providing further evidence that an increase in the M2/M1 macrophage ratio is associated with improved local health. Additionally, human adipose spatial mapping indicates that the M2 macrophage subtype is the predominant macrophage [26] and, thus, a more important macrophage phenotype than has been appreciated to date. Even more surprising and perhaps importantly, healthy obesity correlated with higher M2 macrophage levels than seen even in their healthy lean counterparts.
Our observed macrophage and transcriptional shifts occurred without changes in abdominal adipose distribution or ectopic adipose accumulation. Western dietary components redistribute adipose to ectopic sites [27,28] and drive increases in pro-inflammatory MCP-1 and AT macrophages [29,30,31]. Our data were derived from nonhuman primates unexposed to an unhealthy diet. Though the animals consumed a healthy low-fat, low-sugar, high-fiber, and high-protein chow, we recognize that our cohort animals illustrated more abdominal VIS AT accumulation than SQ AT accumulation, which resulted in VIS:SQ ratios over 1. Previous work has shown that human females without cardiovascular disease have an abdominal VIS:SQ AT ratio under 1 [32]. While our animal model demonstrates a higher abdominal VIS:SQ AT ratio, they present the full range of metabolic abnormalities that are present in people, and they demonstrate the same MHO to MUO transition rate [33]. Additionally, our obese animals demonstrated lower VIS:SQ ratios as a result of increased SQ AT expansion with obesity onset. This corresponds with the animals’ consumption of a healthy diet and their healthy environment, and contributed to the premise of this work, which was to evaluate the differing properties of SQ AT as it expands.
Our data suggest that driving AT macrophage polarization toward an anti-inflammatory M2 profile and enhancing AT fatty acid oxidation may improve or maintain tissue health with obesity onset. MHL and MUL SQ AT demonstrated increased percentages of the undefined macrophages compared to obese animals, indicating that, in the absence of obesity and poor diet consumption, fewer inflammatory signals result in reduced macrophage polarization. In line with this, while the MHL SQ AT demonstrated a non-significant increase in intermediate macrophages, no increases in MHL SQ AT pro-inflammatory cytokine protein expression or increases in inflammation-related transcript expression were observed. Current treatment options for MUO persons rest primarily on weight loss. Our findings support the implementation of lifestyle strategies, such as healthy diets and exercise training, which have been shown to induce M2 macrophage polarization and improve AT function [34], recapitulating the MHO AT state without the need for weight reduction. AT mesenchymal stem cells regulate the tissue immune response through exosomes, as stem cells cluster near M2 macrophages and drive additional M2 polarization [26,35]. Stem cell-derived exosomes from healthy adipose were capable of inducing M2 AT macrophage polarization, highlighting a promising future direction for the treatment of MUO persons [36]. Additionally, treatment with bezafibrate induced acyl-CoA dehydrogenase in cultured fibroblasts while decreasing TNF-α [37], replicating an MHO-like physiologic state, where fatty acid degradation is increased. Along these lines, recent work identified that the use of a near-infrared fluorophore (IR-61) increases AT macrophage oxidative phosphorylation and reduces inflammation, and could be a useful agent for the treatment of obesity-related metabolic dysfunction [38].
Long-term follow-up on MHO patients determined that half of MHO individuals do not develop metabolic abnormalities within ten years [33]. Our findings validate the MHO classification and highlight differences in a key endocrine tissue that contributes to maintaining health with obesity. We used a nonhuman primate model that naturally and reliably demonstrates the four health groups and a strict definition of metabolically healthy versus unhealthy obesity, in contrast to studies which have classified MHO people and animals as obese with an additional metabolic derangement [39]. Nonhuman primates are the leading research model in this arena as, in addition to demonstrating health group distributions that are similar to people, their adipose tissue accumulation patterns [40], immune cell profiles [41], and mechanisms of lipolytic control [42] closely match those observed in humans. While the animals in this study were all female and hormones, including 17-β-estradiol, can protect against metabolic derangements [43], the selected NHP model demonstrated distinct health profiles with significant differences in MetS risk factors, which categorized them into the four clinically defined metabolic health groups.
We immunohistologically evaluated the most predominant immune cell type in SQ white adipose but acknowledge that there are significant numbers of B and T cells that reside in and play key immune modulatory and/or functional roles in maintaining metabolic homeostasis in AT [44,45]. An immunohistochemical approach was chosen, as it has been shown to identify more in situ macrophages without disrupting the tissue inflammatory state compared to other methodologies [46]. Additionally, our staining strategy allowed us to capture a spectrum of macrophage subtypes. We used positive co-staining of both CD68 and CD163 to identify all macrophages, as CD163 has been identified as being present in both resident and recruited macrophages, and CD163 cells have been noted to be increased in both healthy and unhealthy adipose [5,47]. We elected to use the transcription factors pSTAT1 and cMAF as co-stains for M1 and M2 macrophages, respectively. M1 macrophage polarization is promoted after IFN-γ/signal transducer and activator of transcription 1 (STAT1) pathway activation [48]. Similarly, cMAF acts as a metabolic checkpoint regulating the TCA cycle and UDP-GlcNAc biosynthesis, promoting M2 macrophage polarization [49,50], compared to CD206, which is a surface marker of both dendritic cells and macrophages and, thus, less specific [51].
We coupled our immunohistochemical approach with tissue cytokine measurements to secondarily assess inflammatory markers associated with macrophage types, and unbiased RNA-seq supported our inflammation-related findings. While this work was underpowered to detect significant differences in tissue cytokines, the trends observed corroborate the significant tissue immune cell profiles, and future work will include increased sample sizes. In line with this, obese animals demonstrated decreases in SQ AT IL-10 compared to lean animals as the obese animals displayed increased macrophage polarization, or response to inflammatory signaling. Specifically, MHL and MUL SQ AT showed increases in undefined macrophages while MHO SQ AT illustrated increases in M2 macrophage polarization. While RNA-seq does not allow for cell type specificity, it provides accurate gene expression quantification without requiring qPCR validation, as strong correlations have been observed between the two methods [52]. Our histological approaches concluded that there was comparable tissue vascularity but could not assess nutritive perfusion, though HIF1a and related genes were not differentially expressed, suggesting that gross hypoxia was not driving AT differences. Similarly, fibrosis by MTC staining identifies predominantly collagen I. Fibrosis and related gene transcripts were not related to health differences nor MHO M2 macrophage enrichment. However, changes in MUO AT genes, including upregulation of COL4A3 and downregulation of MMP7, may indicate dysregulation as previously described [53].
The strengths of our study include the peerless animal model, which allows MHO AT character to be explored when body adipose distributions, adipocyte size and macrophage density are all comparable. Thus, we report that a unique MHO SQ AT transcript expression profile suggests that a preferred adipocyte function co-exists with a specific immune cell profile. In conclusion, we identify M2 macrophages as an SQ AT component that corresponds with sustained metabolic health in the obese state and a novel target for improving health under such conditions.

4. Methods and Materials

4.1. Study Design and Cohort Selection

A 44-animal age-matched cohort was selected from the Vervet Research Colony (VRC) at Wake Forest University School of Medicine (WFSM). We used a modified version of the National Cholesterol Education Program (NCEP) ATP III definition [54] to categorize our animals, as these colony-raised monkeys exclusively consume a healthy low-fat, low-cholesterol and low-sugar diet and live in a controlled indoor/outdoor setting for the entirety of their lives. We used metabolic syndrome (MetS) criteria, with criteria including the waist measure adjusted for NHPs (waist circumference ≥ 40 cm), HbA1c > 6, fasting glucose ≥ 100 mg/dL, systolic blood pressure > 135 mmHg, diastolic blood pressure > 85 mmHg, triglycerides ≥ 125 mg/dL, HDL cholesterol ≤ 50 mg/dL to identify animals in each health group [15]. MetS scores were calculated as the sum of each MetS criterion met [15]. MHL animals had a MetS score of zero, MHO animals had a MetS score of one for their obesity status only, and MUL and MUO animals had MetS scores of two or higher. We have previously described that African green monkeys demonstrate metabolic health group distributions and MetS risk factors similar to people [15,55,56]. MetS in people is typically defined as the presence of three or more MetS risk factors. However, we found the percentage of African green monkeys with two or more metabolic syndrome criteria to be 29%, which aligns with the percentage of Americans identified as having MetS [57] and formed the basis of our groupings. Given that our animals are born and housed in such healthy environments, we found that 36% of screened African green monkeys demonstrate zero MetS risk factors [15], which is three-fold higher than the prevalence of American adults without MetS risk [56]. Accordingly, a MetS score of two or higher indicates overt unhealthiness in this animal model.
Waist circumference determined obesity. A waist circumference of ≥40 cm corresponds with the upper 20th percentile of the VRC animals [55]. This percentile cutoff was modeled according to where the Adult Treatment Panel III risk waist definition falls for males and females in normal human populations [55,58,59,60]. A waist circumference of 40 cm represents mild obesity, while a waist circumference of ≥44.5 cm represents severe obesity and the upper 10th percentile of the VRC animals. Three animals in each obese group were classified as severely obese. Groups were selected to balance their relatedness coefficients (Supplementary Table S7) to avoid family overrepresentation in any phenotypic group. All animals were socially housed and consumed a low-fat, low-cholesterol, high protein and fiber chow with no added sugar (Lab Diet 5038; LabDiet).
We used the animals’ group placements at the time of biopsy when conducting our analyses, allowing for 9–13 observations per group for each endpoint. As with many NHP studies, selection from the breeding colony led to our study including only females.
All procedures were performed in accordance with the Guide for Care and Use of Laboratory Animals. Protocols for the avoidance of pain and discomfort were adhered to and conducted in compliance with the WFSM Institutional Animal Care and Use Committee.

4.2. Blood and Metabolic Measures

Animals were fasted overnight and anesthetized with intramuscular ketamine injections (10–15 mg/kg). Animals were weighed and had waist circumference measured with a flexible tape measure at the level of the umbilicus. Blood samples were obtained via percutaneous femoral venipuncture into ethylenediaminetetraacetic acid (EDTA) blood tubes and placed on ice until processed for plasma, which was stored at −80 °C until analyzed. Fasting glucose was determined by the glucose oxidase method, and fasting insulin concentrations were determined by enzyme-linked immunosorbent assay performed using isolated plasma (ELISA; Mercodia, Winston-Salem, NC, USA; cat. no. 10-1113-01). Whole blood was used to measure glycated hemoglobin using high-performance liquid chromatography methodology (Afinion Alere, Abbott, Chicago, IL, USA). Triglyceride, high-density lipoprotein cholesterol (HDLc), circulating free fatty acid (Abcam, Cambridge, UK; cat. no. ab65341), and total cholesterol concentrations were measured enzymatically using isolated plasma [61]. Blood pressure was measured indirectly by sphygmomanometer as previously described [61]. Supplementary Figure S8 illustrations the methods performed on collected blood samples and isolated plasma.

4.3. Body Composition

Computed tomography (CT) scans were performed on a GE 0.625-mm slice Discovery MI DR scanner (GE Healthcare, Chicago, IL, USA). Images were reconstructed with TeraRecon Aquarius Intuition software (Durham, NC, USA) and converted into dicom format for analysis. Fat- and lean-containing voxels were used to calculate volume results and expressed as a percentage of the animal’s total body weight [62]. Abdominal fat distribution was assessed by regional analysis between the thoracolumbar junction and the fourth lumbar vertebrae using automated differentiation of deposits within and outside the body wall [62]. Liver attenuation values were collected in triplicate [62], where liver attenuation, or density, was used to determine liver adipose accumulation.

4.4. Tissue Collections

Animals under isoflurane anesthesia had SQ AT collected from the ventral aspect of the caudal abdomen at the umbilicus. Upon collection, tissues were washed twice in Dulbecco’s phosphate-buffered saline (DPBS) and either flash frozen in liquid nitrogen or placed in 4% paraformaldehyde (PFA). Supplementary Figure S8 depicts the methods performed on collected adipose tissues.

4.5. Cytokine Evaluations

Protein was extracted from ~100 mg of frozen tissue using the Invent Minute Total Protein Extraction Kit for Adipose Tissues (Plymouth, MN, USA; cat. no. AT-022) according to the manufacturer’s instructions. Extracted protein was quantitated by bicinchoninic acid assay. A total of 75 μg of protein was used to analyze monocyte chemoattractant protein-1 (MCP-1), transforming growth factor beta (TGF-β), interleukin (IL)-6, (R & D Systems, Minneapolis, MN, USA; cat. nos. DCP00, DB100B, and D6050) and IL-10 by ELISA (U-Cytech Biosciences, Utrecht, The Netherlands; cat. no. CT146). Frozen plasma samples were used to evaluate circulating MCP-1, IL-6, interleukin (IL)-1β (cat. no. DLB50), plasminogen activator inhibitor (PAI)-1 (cat. no. DSE100), and adiponectin (cat. no. DRP300) via ELISA (R & D Systems, Minneapolis, MN, USA). Assays were run in triplicate using predetermined dilutions.

4.6. Histology and Immunohistochemistry

PFA-fixed, paraffin-embedded tissues were used to create histologic sections. Sections were stained with fluorescent CD68 and CD163 to evaluate total macrophage infiltration (AbD Serotec Bio-Rad, Hercules, CA, USA; 1 μg/mL), and co-stained with pSTAT1 (Cell Signaling Technology, Danvers, MA, USA, cat. no. 9167S; 1:250 dilution) and CMAF (Abcam, Cambridge, UK; cat. no. ab199424; 1:250 dilution) identify undefined, M1, intermediate, and M2 macrophages [47]. Undefined macrophages were identified as CD68+CD163+, M1s were identified as CD68+CD163+pSTAT1+, intermediate macrophages were identified as CD68+CD163+pSTAT1+CMAF+, and M2s were identified as CD68+CD163+CMAF+. A Visiopharm software macro (v.2019.02.1.6005) distinguished between macrophage types. Macrophage subtypes were identified across an entire tissue section. Macrophage subtype percentages were determined by taking the number of the specific macrophage subtype and dividing by the total number of macrophages identified in a tissue section. Macrophage density was determined by taking the total number of macrophages for a tissue section and dividing by the tissue area, which was calculated using a specialized Visiopharm macro. Unstained sections and lung tissue were used as negative and positive controls, respectively. Adipocyte analyses were performed using a Visiopharm (v.2019.02.1.6005) macro on hemoxylin and eosin-stained slides that determined total tissue area, total adipocyte number, total adipocyte area, and individual adipocyte area. Adipose sections were stained with Masson’s Trichrome (Abcam, Cambridge, UK) and CD31 (Biocare Medical, Pacheco, CA, USA; cat. no. CM131C) and evaluated using Visiopharm.

4.7. Transcriptomic Analyses

RNA was extracted from 100mg adipose biopsies (n = 43) using the Qiagen RNeasy Lipid Mini tissue kit (Germantown, MD, USA; cat. no. 74804) according to the manufacturer’s protocol. Genomic DNA was removed using the Qiagen RNase-Free DNase set according to the manufacturer’s protocol. Total RNA was used to prepare cDNA libraries using the Illumina® TruSeq Stranded mRNA Library Prep and IDT for Illumina-TruSeq RNA UD Indexes. RNA quality was assessed using the Agilent TapeStation 4150. RIN values ranged from 7.2 to 9.2. Poly-A selection using oligo-dT magnetic beads occurred and was followed by enzymatic fragmentation, reverse-transcription and double-stranded cDNA purification using AMPure XP magnetic beads. The cDNA was end repaired, 3′ adenylated with Illumina sequencing adaptors ligated onto the fragment ends, and the stranded libraries were pre-amplified with PCR. The library size distribution was validated and quality inspected using an Agilent 2100 Bioanalyzer. The quantity of each cDNA library was measured using the Qubit 3.0 (Thermo Fisher, Waltham, MA, USA). The libraries were pooled and sequenced to a target read depth of 25M reads per library using single-end, NovaSeq 6000 SP Reagent Kit (100 cycles) (Illumina Inc., San Diego, CA, USA) on the NovaSeq 6000. RNA-seq data were screened using FastQC for quality analysis. The STAR sequence aligner and the vervet monkey genome (Chlorocebus aethiops) were used to align clean reads [63]. 32.52 ± 5.16 million reads were uniquely mapped. “featureCounts” quantified RNA-seq data as counts [64]. Differentially expressed genes (DEGs) were analyzed using the Limma-voom tool in R, where gene counts were used as inputs [65]. Database for Annotation, Visualization and Integrated Discovery (DAVID) was to identify biologically enriched gene ontology (GO) terms (https://david.ncifcrf.gov/) (accessed on 1 November 2021).

4.8. Statistical Methods

The Shapiro–Wilk and Levene’s tests were used to assess normality and homogeneity of variance for histological and protein endpoints. Adipocytes under 400 um2 in size were removed prior to downstream analyses. Endpoints that did not meet assumptions underwent logarithmic transformation. Analysis of variance was performed on all histological and protein endpoints to determine group differences. Factorial analysis of variance determined the main effects of obesity and health status. Pearson’s correlation was used to determine associations. A p-value of <0.05 was used to determine significance and a p-value of <0.1 was used to determine a trend. Significant DEGs were defined as p < 0.01 and log2 fold change > |0.585| (1.5-fold change).

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232012680/s1.

Author Contributions

Authors involved in data generation and data analysis: A.D.R., R.V., J.C., A.B., A.W., M.D., D.D., S.K.D. and K.K. Authors involved in sample collection and trial conductance: A.D.R., R.V., M.B., M.D. and K.K. Authors involved in manuscript writing and editing: A.D.R., R.V. and K.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was support by the National Institutes of Health [R01HL142930, P40OD010965, ULTR001420, T32AI007401].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All transcriptomic data will be made available upon request.

Acknowledgments

The authors wish to acknowledge our funding sources (R01HL142930, P40OD010965, ULTR001420, T32AI007401), all technicians involved in the care and evaluation of the study animals, and the nonhuman primates that made this study possible.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization. Obesity and Overweight. Available online: https://www.who.int/news-room/fact-sheets/detail/obesity-and-overweight (accessed on 18 October 2022).
  2. Aung, K.; Lorenzo, C.; Hinojosa, M.A.; Haffner, S.M. Risk of developing diabetes and cardiovascular disease in metabolically unhealthy normal-weight and metabolically healthy obese individuals. J. Clin. Endocrinol. Metab. 2014, 99, 462–468. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Wahba, I.M.; Mak, R.H. Obesity and obesity-initiated metabolic syndrome: Mechanistic links to chronic kidney disease. Clin. J. Am. Soc. Nephrol. 2007, 2, 550–562. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Weisberg, S.P.; McCann, D.; Desai, M.; Rosenbaum, M.; Leibel, R.L.; Ferrante, A.W. Obesity is associated with macrophage accumulation in adipose tissue. J. Clin. Investig. 2003, 112, 1796–1808. [Google Scholar] [CrossRef] [PubMed]
  5. Russo, L.; Lumeng, C.N. Properties and functions of adipose tissue macrophages in obesity. Immunology 2018, 155, 407–417. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Chylikova, J.; Dvorackova, J.; Tauber, Z.; Kamarad, V. M1/M2 macrophage polarization in human obese adipose tissue. Biomed. Pap. Med. Fac. Univ. Palacky Olomouc Czech. Repub. 2018, 162, 79–82. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Wang, Y.; Smith, W.; Hao, D.; He, B.; Kong, L. M1 and M2 macrophage polarization and potentially therapeutic naturally occurring compounds. Int. Immunopharmacol. 2019, 70, 459–466. [Google Scholar] [CrossRef]
  8. Wentworth, J.M.; Naselli, G.; Brown, W.A.; Doyle, L.; Phipson, B.; Smyth, G.K.; Wabitsch, M.; O’Brien, P.E.; Harrison, L.C. Pro-inflammatory CD11c+ CD206+ adipose tissue macrophages are associated with insulin resistance in human obesity. Diabetes 2010, 59, 1648–1656. [Google Scholar] [CrossRef] [Green Version]
  9. Wynn, T.A.; Vannella, K.M. Macrophages in tissue repair, regeneration, and fibrosis. Immunity 2016, 44, 450–462. [Google Scholar] [CrossRef] [Green Version]
  10. Brestoff, J.R.; Wilen, C.B.; Moley, J.R.; Li, Y.; Zou, W.; Malvin, N.P.; Rowen, M.N.; Saunders, B.T.; Ma, H.; Mack, M.R. Intercellular mitochondria transfer to macrophages regulates white adipose tissue homeostasis and is impaired in obesity. Cell Metab. 2021, 33, 270–282.e278. [Google Scholar] [CrossRef]
  11. Nguyen, K.D.; Qiu, Y.; Cui, X.; Goh, Y.; Mwangi, J.; David, T.; Mukundan, L.; Brombacher, F.; Locksley, R.M.; Chawla, A. Alternatively activated macrophages produce catecholamines to sustain adaptive thermogenesis. Nature 2011, 480, 104–108. [Google Scholar] [CrossRef]
  12. Gao, J.; Qiao, Z.; Liu, S.; Xu, J.; Wang, S.; Yang, X.; Wang, X.; Tang, R. A small molecule nanodrug consisting of pH-sensitive ortho ester–dasatinib conjugate for cancer therapy. Eur. J. Pharm. Biopharm. 2021, 163, 188–197. [Google Scholar] [CrossRef] [PubMed]
  13. Wellen, K.E.; Hotamisligil, G.S. Obesity-induced inflammatory changes in adipose tissue. J. Clin. Investig. 2003, 112, 1785–1788. [Google Scholar] [CrossRef] [PubMed]
  14. Paczkowska-Abdulsalam, M.; Kretowski, A. Obesity, metabolic health and omics: Current status and future directions. World J. Diabetes 2021, 12, 420. [Google Scholar] [CrossRef] [PubMed]
  15. Kavanagh, K.; Davis, A.T.; Peters, D.E.; LeGrand, A.C.; Bharadwaj, M.S.; Molina, A.J. Regulators of mitochondrial quality control differ in subcutaneous fat of metabolically healthy and unhealthy obese monkeys. Obesity 2017, 25, 689–696. [Google Scholar] [CrossRef] [Green Version]
  16. Guo, M.; Härtlova, A.; Gierliński, M.; Prescott, A.; Castellvi, J.; Losa, J.H.; Petersen, S.K.; Wenzel, U.A.; Dill, B.D.; Emmerich, C.H. Triggering MSR1 promotes JNK-mediated inflammation in IL-4-activated macrophages. EMBO J. 2019, 38, e100299. [Google Scholar] [CrossRef]
  17. De Oliveira, C.; Khatua, B.; Noel, P.; Kostenko, S.; Bag, A.; Balakrishnan, B.; Patel, K.S.; Guerra, A.A.; Martinez, M.N.; Trivedi, S. Pancreatic triglyceride lipase mediates lipotoxic systemic inflammation. J. Clin. Investig. 2020, 130, 1931–1947. [Google Scholar] [CrossRef] [Green Version]
  18. Mehta, J.L.; Chen, J.; Hermonat, P.L.; Romeo, F.; Novelli, G. Lectin-like, oxidized low-density lipoprotein receptor-1 (LOX-1): A critical player in the development of atherosclerosis and related disorders. Cardiovasc. Res. 2006, 69, 36–45. [Google Scholar] [CrossRef]
  19. Das, S.K.; Ma, L.; Sharma, N.K. Adipose tissue gene expression and metabolic health of obese adults. Int. J. Obes. 2015, 39, 869–873. [Google Scholar] [CrossRef] [Green Version]
  20. 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]
  21. Naukkarinen, J.; Heinonen, S.; Hakkarainen, A.; Lundbom, J.; Vuolteenaho, K.; Saarinen, L.; Hautaniemi, S.; Rodriguez, A.; Frühbeck, G.; Pajunen, P. Characterising metabolically healthy obesity in weight-discordant monozygotic twins. Diabetologia 2014, 57, 167–176. [Google Scholar] [CrossRef]
  22. Schlauch, K.A.; Kulick, D.; Subramanian, K.; De Meirleir, K.L.; Palotás, A.; Lombardi, V.C. Single-nucleotide polymorphisms in a cohort of significantly obese women without cardiometabolic diseases. Int. J. Obes. 2019, 43, 253–262. [Google Scholar] [CrossRef] [PubMed]
  23. Gharooi Ahangar, O.; Javanrouh, N.; Daneshpour, M.S.; Barzin, M.; Valizadeh, M.; Azizi, F.; Hosseinpanah, F. Genetic markers and continuity of healthy metabolic status: Tehran cardio-metabolic genetic study (TCGS). Sci. Rep. 2020, 10, 1–8. [Google Scholar]
  24. Fujisaka, S. The role of adipose tissue M1/M2 macrophages in type 2 diabetes mellitus. Diabetol. Int. 2021, 12, 74–79. [Google Scholar] [CrossRef] [PubMed]
  25. Miyachi, Y.; Tsuchiya, K.; Shiba, K.; Mori, K.; Komiya, C.; Ogasawara, N.; Ogawa, Y. A reduced M1-like/M2-like ratio of macrophages in healthy adipose tissue expansion during SGLT2 inhibition. Sci. Rep. 2018, 8, 1–13. [Google Scholar] [CrossRef] [Green Version]
  26. Bäckdahl, J.; Franzén, L.; Massier, L.; Li, Q.; Jalkanen, J.; Gao, H.; Andersson, A.; Bhalla, N.; Thorell, A.; Rydén, M. Spatial mapping reveals human adipocyte subpopulations with distinct sensitivities to insulin. Cell Metab. 2021, 33, 1869–1882.e1866. [Google Scholar] [CrossRef]
  27. Kavanagh, K.; Jones, K.L.; Sawyer, J.; Kelley, K.; Carr, J.J.; Wagner, J.D.; Rudel, L.L. Trans fat diet induces abdominal obesity and changes in insulin sensitivity in monkeys. Obesity 2007, 15, 1675–1684. [Google Scholar] [CrossRef]
  28. Kavanagh, K.; Wylie, A.T.; Tucker, K.L.; Hamp, T.J.; Gharaibeh, R.Z.; Fodor, A.A.; Cullen, J.M. Dietary fructose induces endotoxemia and hepatic injury in calorically controlled primates. Am. J. Clin. Nutr. 2013, 98, 349–357. [Google Scholar] [CrossRef] [Green Version]
  29. Nappo, F.; Esposito, K.; Cioffi, M.; Giugliano, G.; Molinari, A.M.; Paolisso, G.; Marfella, R.; Giugliano, D. Postprandial endothelial activation in healthy subjects and in type 2 diabetic patients: Role of fat and carbohydrate meals. J. Am. Coll. Cardiol. 2002, 39, 1145–1150. [Google Scholar] [CrossRef] [Green Version]
  30. Spencer, M.; Finlin, B.S.; Unal, R.; Zhu, B.; Morris, A.J.; Shipp, L.R.; Lee, J.; Walton, R.G.; Adu, A.; Erfani, R. Omega-3 fatty acids reduce adipose tissue macrophages in human subjects with insulin resistance. Diabetes 2013, 62, 1709–1717. [Google Scholar] [CrossRef] [Green Version]
  31. Brykczynska, U.; Geigges, M.; Wiedemann, S.J.; Dror, E.; Böni-Schnetzler, M.; Hess, C.; Donath, M.Y.; Paro, R. Distinct transcriptional responses across tissue-resident macrophages to short-term and long-term metabolic challenge. Cell Rep. 2020, 30, 1627–1643.e1627. [Google Scholar] [CrossRef]
  32. Fox, C.S.; Massaro, J.M.; Hoffmann, U.; Pou, K.M.; Maurovich-Horvat, P.; Liu, C.-Y.; Vasan, R.S.; Murabito, J.M.; Meigs, J.B.; Cupples, L.A. Abdominal visceral and subcutaneous adipose tissue compartments: Association with metabolic risk factors in the Framingham Heart Study. Circulation 2007, 116, 39–48. [Google Scholar] [CrossRef] [PubMed]
  33. Bell, J.A.; Hamer, M.; Batty, G.D.; Singh-Manoux, A.; Sabia, S.; Kivimäki, M. Incidence of metabolic risk factors among healthy obese adults: 20-year follow-up. J. Am. Coll. Cardiol. 2015, 66, 871–873. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Kolahdouzi, S.; Talebi-Garakani, E.; Hamidian, G.; Safarzade, A. Exercise training prevents high-fat diet-induced adipose tissue remodeling by promoting capillary density and macrophage polarization. Life Sci. 2019, 220, 32–43. [Google Scholar] [CrossRef]
  35. Heo, J.S.; Choi, Y.; Kim, H.O. Adipose-derived mesenchymal stem cells promote M2 macrophage phenotype through exosomes. Stem Cells Int. 2019, 2019, 7921760. [Google Scholar] [CrossRef] [PubMed]
  36. Zhao, H.; Shang, Q.; Pan, Z.; Bai, Y.; Li, Z.; Zhang, H.; Zhang, Q.; Guo, C.; Zhang, L.; Wang, Q. Exosomes from adipose-derived stem cells attenuate adipose inflammation and obesity through polarizing M2 macrophages and beiging in white adipose tissue. Diabetes 2018, 67, 235–247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Cabrero, A.; Alegret, M.; Sánchez, R.M.; Adzet, T.; Laguna, J.C.; Vázquez, M. Bezafibrate reduces mRNA levels of adipocyte markers and increases fatty acid oxidation in primary culture of adipocytes. Diabetes 2001, 50, 1883–1890. [Google Scholar] [CrossRef]
  38. Wang, Y.; Tang, B.; Long, L.; Luo, P.; Xiang, W.; Li, X.; Wang, H.; Jiang, Q.; Tan, X.; Luo, S. Improvement of obesity-associated disorders by a small-molecule drug targeting mitochondria of adipose tissue macrophages. Nat. Commun. 2021, 12, 1–16. [Google Scholar] [CrossRef]
  39. Camhi, S.M.; Must, A.; Gona, P.N.; Hankinson, A.; Odegaard, A.; Reis, J.; Gunderson, E.P.; Jacobs, D.R.; Carnethon, M.R. Duration and stability of metabolically healthy obesity over 30 years. Int. J. Obes. 2019, 43, 1803–1810. [Google Scholar] [CrossRef] [PubMed]
  40. Comuzzie, A.G.; Cole, S.A.; Martin, L.; Carey, K.D.; Mahaney, M.C.; Blangero, J.; VandeBerg, J.L. The baboon as a nonhuman primate model for the study of the genetics of obesity. Obes. Res. 2003, 11, 75–80. [Google Scholar] [CrossRef]
  41. Laparra, A.; Tricot, S.; Le Van, M.; Damouche, A.; Gorwood, J.; Vaslin, B.; Favier, B.; Benoist, S.; Ho Tsong Fang, R.; Bosquet, N. The frequencies of immunosuppressive cells in adipose tissue differ in human, non-human primate, and mouse models. Front. Immunol. 2019, 10, 117. [Google Scholar] [CrossRef]
  42. Bousquet-Mélou, A.; Galitzky, J.; Lafontan, M.; Berlan, M. Control of lipolysis in intra-abdominal fat cells of nonhuman primates: Comparison with humans. J. Lipid Res. 1995, 36, 451–461. [Google Scholar] [CrossRef]
  43. Ziaei, S.; Mohseni, H. Correlation between hormonal statuses and metabolic syndrome in postmenopausal women. J. Fam. Reprod. Health 2013, 7, 63. [Google Scholar]
  44. Misumi, I.; Starmer, J.; Uchimura, T.; Beck, M.A.; Magnuson, T.; Whitmire, J.K. Obesity expands a distinct population of T cells in adipose tissue and increases vulnerability to infection. Cell Rep. 2019, 27, 514–524.e515. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Wu, L.; Parekh, V.V.; Hsiao, J.; Kitamura, D.; Van Kaer, L. Spleen supports a pool of innate-like B cells in white adipose tissue that protects against obesity-associated insulin resistance. Proc. Natl. Acad. Sci. USA 2014, 111, E4638–E4647. [Google Scholar] [CrossRef] [Green Version]
  46. Morgan-Bathke, M.; Harteneck, D.; Jaeger, P.; Sondergaard, E.; Karwoski, R.; Espinosa De Ycaza, A.; Carranza-Leon, B.G.; Faubion, W.A., Jr.; Oliveira, A.M.; Jensen, M.D. Comparison of methods for analyzing human adipose tissue macrophage content. Obesity 2017, 25, 2100–2107. [Google Scholar] [CrossRef] [PubMed]
  47. Barros, M.H.M.; Hauck, F.; Dreyer, J.H.; Kempkes, B.; Niedobitek, G. Macrophage polarisation: An immunohistochemical approach for identifying M1 and M2 macrophages. PLoS ONE 2013, 8, e80908. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Leopold Wager, C.M.; Hole, C.R.; Wozniak, K.L.; Olszewski, M.A.; Mueller, M.; Wormley, F.L., Jr. STAT1 signaling within macrophages is required for antifungal activity against Cryptococcus neoformans. Infect. Immun. 2015, 83, 4513–4527. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Liu, M.; Tong, Z.; Ding, C.; Luo, F.; Wu, S.; Wu, C.; Albeituni, S.; He, L.; Hu, X.; Tieri, D. Transcription factor c-Maf is a checkpoint that programs macrophages in lung cancer. J. Clin. Investig. 2020, 130, 2081–2096. [Google Scholar] [CrossRef] [Green Version]
  50. Chen, L.; Wang, J.; Li, S.; Yu, Z.; Liu, B.; Song, B.; Su, Y. The clinical dynamic changes of macrophage phenotype and function in different stages of human wound healing and hypertrophic scar formation. Int. Wound J. 2019, 16, 360–369. [Google Scholar] [CrossRef]
  51. Xu, Z.-J.; Gu, Y.; Wang, C.-Z.; Jin, Y.; Wen, X.-M.; Ma, J.-C.; Tang, L.-J.; Mao, Z.-W.; Qian, J.; Lin, J. The M2 macrophage marker CD206: A novel prognostic indicator for acute myeloid leukemia. Oncoimmunology 2020, 9, 1683347. [Google Scholar] [CrossRef] [Green Version]
  52. Everaert, C.; Luypaert, M.; Maag, J.L.; Cheng, Q.X.; Dinger, M.E.; Hellemans, J.; Mestdagh, P. Benchmarking of RNA-sequencing analysis workflows using whole-transcriptome RT-qPCR expression data. Sci. Rep. 2017, 7, 1–11. [Google Scholar]
  53. Ruiz-Ojeda, F.J.; Méndez-Gutiérrez, A.; Aguilera, C.M.; Plaza-Díaz, J. Extracellular matrix remodeling of adipose tissue in obesity and metabolic diseases. Int. J. Mol. Sci. 2019, 20, 4888. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Huang, P.L. A comprehensive definition for metabolic syndrome. Dis. Model. Mech. 2009, 2, 231–237. [Google Scholar] [CrossRef]
  55. Kavanagh, K.; Fairbanks, L.A.; Bailey, J.N.; Jorgensen, M.J.; Wilson, M.; Zhang, L.; Rudel, L.L.; Wagner, J.D. Characterization and heritability of obesity and associated risk factors in vervet monkeys. Obesity 2007, 15, 1666–1674. [Google Scholar] [CrossRef]
  56. Araújo, J.; Cai, J.; Stevens, J. Prevalence of optimal metabolic health in American adults: National Health and Nutrition Examination Survey 2009–2016. Metab. Syndr. Relat. Disord. 2019, 17, 46–52. [Google Scholar] [CrossRef]
  57. Moore, J.X.; Chaudhary, N.; Akinyemiju, T. Peer reviewed: Metabolic syndrome prevalence by race/ethnicity and sex in the United States, National Health and Nutrition Examination Survey, 1988–2012. Prev. Chronic Dis. 2017, 14, E24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. National Cholesterol Education Program (US); Expert Panel on Detection, Treatment of High Blood Cholesterol in Adults. Third Report of the National Cholesterol Education Program (NCEP) Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol in Adults (Adult Treatment Panel III); The Program: Bethesda, MD, USA, 2002. [Google Scholar]
  59. Lean, M.; Han, T.; Morrison, C. Waist circumference as a measure for indicating need for weight management. BMJ 1995, 311, 158–161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Palaniappan, L.; Carnethon, M.R.; Wang, Y.; Hanley, A.J.; Fortmann, S.P.; Haffner, S.M.; Wagenknecht, L. Predictors of the incident metabolic syndrome in adults: The Insulin Resistance Atherosclerosis Study. Diabetes Care 2004, 27, 788–793. [Google Scholar] [CrossRef] [Green Version]
  61. Vemuri, R.; Ruggiero, A.; Whitfield, J.B.; Dugan, G.O.; Cline, J.M.; Block, M.R.; Guo, H.; Kavanagh, K. Hypertension drives microbial translocation and shifts in the fecal microbiome of non-human primates. bioRxiv 2021. [Google Scholar] [CrossRef]
  62. Bacarella, N.; Ruggiero, A.; Davis, A.T.; Uberseder, B.; Davis, M.A.; Bracy, D.P.; Wasserman, D.H.; Cline, J.M.; Sherrill, C.; Kavanagh, K. Whole body irradiation induces diabetes and adipose insulin resistance in nonhuman primates. Int. J. Radiat. Oncol. Biol. Phys. 2020, 106, 878–886. [Google Scholar] [CrossRef]
  63. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  64. Liao, Y.; Smyth, G.K.; Shi, W. featureCounts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014, 30, 923–930. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef] [PubMed]
Ijms 23 12680 g001 550
Figure 1. Nonhuman primate cohort body composition and adipocyte sizes. All values are presented as means ± SEM. Factorial ANOVA and ANCOVA with Tukey post hoc analyses were used to determine differences by obesity status and health status. Unlike letters denote statistically significant differences (p < 0.05) by health status or obesity status, specifically between groups labeled “a” versus “b”. Groups labeled with the same letter were not statistically different. (a) Metabolic syndrome scores for each health group (p = 0.00016; MHL, n = 12; MHO, n = 10; MUL, n = 9; MUO, n = 13). Differences in (b) fat to lean mass ratio (p ≤ 0.0001), (c) abdominal subcutaneous (SQ) adipose volume (p = 0.00023), and (d) abdominal visceral (VIS) adipose volume (p = 0.0095) between health groups. (e) Histological quantification of adipocyte sizes (p = 0.00056). The violin plot provided illustrates the distribution of adipocyte sizes in each health group. The cross bars indicate the mean adipocyte size value for each group.
Figure 1. Nonhuman primate cohort body composition and adipocyte sizes. All values are presented as means ± SEM. Factorial ANOVA and ANCOVA with Tukey post hoc analyses were used to determine differences by obesity status and health status. Unlike letters denote statistically significant differences (p < 0.05) by health status or obesity status, specifically between groups labeled “a” versus “b”. Groups labeled with the same letter were not statistically different. (a) Metabolic syndrome scores for each health group (p = 0.00016; MHL, n = 12; MHO, n = 10; MUL, n = 9; MUO, n = 13). Differences in (b) fat to lean mass ratio (p ≤ 0.0001), (c) abdominal subcutaneous (SQ) adipose volume (p = 0.00023), and (d) abdominal visceral (VIS) adipose volume (p = 0.0095) between health groups. (e) Histological quantification of adipocyte sizes (p = 0.00056). The violin plot provided illustrates the distribution of adipocyte sizes in each health group. The cross bars indicate the mean adipocyte size value for each group.
Ijms 23 12680 g001
Ijms 23 12680 g002 550
Figure 2. Metabolically healthy obese subcutaneous adipose demonstrates increases in anti-inflammatory M2 macrophages. (a) The ratio of M2/M1 macrophages between metabolic health groups (overall analysis of variance (ANOVA) p = 0.02; health effect p = 0.007). M2/M1 ratios are presented as means ± SEM (metabolically healthy lean (MHL), n = 12; metabolically healthy obese (MHO), n = 10; metabolically unhealthy lean (MUL), n = 9; metabolically unhealthy obese (MUO), n = 13). * Indicates Tukey post hoc significance of p < 0.05 compared to all other groups. (b) Stacked bars of the macrophage subtype percentages found in the SQ AT of each group. (c) Representative M2 macrophage AT heat maps, where increasing green to red coloration indicates higher densities of M2 macrophages.
Figure 2. Metabolically healthy obese subcutaneous adipose demonstrates increases in anti-inflammatory M2 macrophages. (a) The ratio of M2/M1 macrophages between metabolic health groups (overall analysis of variance (ANOVA) p = 0.02; health effect p = 0.007). M2/M1 ratios are presented as means ± SEM (metabolically healthy lean (MHL), n = 12; metabolically healthy obese (MHO), n = 10; metabolically unhealthy lean (MUL), n = 9; metabolically unhealthy obese (MUO), n = 13). * Indicates Tukey post hoc significance of p < 0.05 compared to all other groups. (b) Stacked bars of the macrophage subtype percentages found in the SQ AT of each group. (c) Representative M2 macrophage AT heat maps, where increasing green to red coloration indicates higher densities of M2 macrophages.
Ijms 23 12680 g002
Ijms 23 12680 g003a 550Ijms 23 12680 g003b 550
Figure 3. Metabolically healthy obese subcutaneous adipose illustrates increases in fatty acid oxidation-related gene ontology terms and downregulation of pro-inflammatory transcripts. Red coloration indicates pathway upregulation while green coloration indicates pathway downregulation. Pathways shown all survived false discovery rate and multiple comparisons corrections, and reached statistical significance (p < 0.05). (a) Metabolically healthy obese (MHO; n = 10) subcutaneous (SQ) adipose tissue (AT) pathway increases compared to metabolically healthy lean (MHL; n = 12). (b) Metabolically unhealthy obese (MUO, n = 13) SQ AT pathway increases and decreases compared to MHL (n = 12). (c) Heat map of differentially expressed genes in MHO SQ AT compared to MUO (MUO n = 13; MHO n = 10). Specific transcripts are not denoted in this heat map.
Figure 3. Metabolically healthy obese subcutaneous adipose illustrates increases in fatty acid oxidation-related gene ontology terms and downregulation of pro-inflammatory transcripts. Red coloration indicates pathway upregulation while green coloration indicates pathway downregulation. Pathways shown all survived false discovery rate and multiple comparisons corrections, and reached statistical significance (p < 0.05). (a) Metabolically healthy obese (MHO; n = 10) subcutaneous (SQ) adipose tissue (AT) pathway increases compared to metabolically healthy lean (MHL; n = 12). (b) Metabolically unhealthy obese (MUO, n = 13) SQ AT pathway increases and decreases compared to MHL (n = 12). (c) Heat map of differentially expressed genes in MHO SQ AT compared to MUO (MUO n = 13; MHO n = 10). Specific transcripts are not denoted in this heat map.
Ijms 23 12680 g003aIjms 23 12680 g003b
Ijms 23 12680 g004 550
Figure 4. Endothelial area increases with increasing subcutaneous adipose adipocyte size, while fibrosis decreases with increasing subcutaneous adipose adipocyte size. (a) Correlation between CD31-positive staining area and SQ adipocyte size (Pearson’s r = 0.33, p = 0.03, n = 43). (b) Correlation between transformed percent MTC-positive area and SQ adipocyte area (Pearson’s r = −0.34, p = 0.03, n = 43).
Figure 4. Endothelial area increases with increasing subcutaneous adipose adipocyte size, while fibrosis decreases with increasing subcutaneous adipose adipocyte size. (a) Correlation between CD31-positive staining area and SQ adipocyte size (Pearson’s r = 0.33, p = 0.03, n = 43). (b) Correlation between transformed percent MTC-positive area and SQ adipocyte area (Pearson’s r = −0.34, p = 0.03, n = 43).
Ijms 23 12680 g004
Table
Table 1. Nonhuman primate cohort demographic information. Values are presented as means with SEM in parentheses. ANOVA and factorial ANOVA were used to identify differences between groups and by health and obesity status. A total of 44 female African green monkeys were selected for investigation. Metabolically healthy lean (MHL; n = 12); metabolically healthy obese (MHO; n = 10); metabolically unhealthy lean (MUL; n = 9); metabolically unhealthy obese (MUO; n = 13).
Table 1. Nonhuman primate cohort demographic information. Values are presented as means with SEM in parentheses. ANOVA and factorial ANOVA were used to identify differences between groups and by health and obesity status. A total of 44 female African green monkeys were selected for investigation. Metabolically healthy lean (MHL; n = 12); metabolically healthy obese (MHO; n = 10); metabolically unhealthy lean (MUL; n = 9); metabolically unhealthy obese (MUO; n = 13).
MeasurementMHLMHOMULMUOANOVA p-ValueHealth p-ValueObesity p-Value
Age (yrs)15.15 (1.30)15.14 (1.55)14.32 (1.99)15.87 (0.92)0.740.690.26
Body Weight (kg)5.25 (0.26)6.47 (0.24)4.95 (0.24)6.91 (0.26)<0.001 0.190.01
Waist Circumference (cm)34.29 (0.97)42.86 (1.07)34.80 (0.78)43.23 (0.85)<0.001 0.02<0.001
Fasting Glucose (mg/dL)73.5 (3.20)67.2 (3.59) 108.67 (19.15)104.32 (12.37)0.01 0.360.86
Glycosylated Hemoglobin A1c (%)4.07 (0.05)4.07 (0.07) 4.83 (0.58)5.38 (0.57)0.009 0.012 0.31
Systolic Blood Pressure (mmHg)114.44 (3.79)110.10 (5.43)144.78 (4.95)111.62 (5.06)<0.001 0.480.31
Diastolic Blood Pressure (mmHg)67.53 (3.38)61.57 (3.80)89.22 (4.70) 66.79 (2.46)<0.001 0.250.42
High Density Lipoprotein Cholesterol (mg/dL)68.08 (3.37)56.10 (5.76)73.44 (10.38)62.50 (4.39)0.240.330.06
Triglycerides (mg/dL)56.34 (3.46)55.90 (5.89)79.10 (15.67)90.52 (10.80)0.030.240.76
Liver Density (HU)60.37 (2.26)55.71 (4.03)59.54 (2.84)57.59 (1.94)0.860.690.83
% Fat Mass18.35 (2.69)36.59 (2.81)16.70 (2.07)35.98 (3.05)0.01 0.940.06
% Lean Mass54.13 (1.56) 47.29 (2.89)55.03 (1.52)45.63 (1.60)0.02 0.910.07
Visceral:Subcutaneous Adipose Area (cm2)3.44 (0.46)1.78 (0.19)3.86 (0.83)1.73 (0.30)0.310.880.19
Plasma Insulin (mU/L)16.80 (4.15)19.98 (3.87)8.84 (2.08)51.12 (11.18)0.0055 0.710.0019
Table
Table 2. Top false discovery rate correction-surviving up- and downregulated transcripts illustrated in the metabolically healthy obese (n = 10) subcutaneous adipose tissue compared to the metabolically healthy lean (n = 12).
Table 2. Top false discovery rate correction-surviving up- and downregulated transcripts illustrated in the metabolically healthy obese (n = 10) subcutaneous adipose tissue compared to the metabolically healthy lean (n = 12).
TranscriptLog Fold ChangeAdjusted p-Value
Fatty Acid Binding Protein 7 (FABP7)1.510.010
Prostaglandin D2 Synthase (PTGDS)1.530.015
Aldehyde Dehydrogenase 1 Family Member B1 (ALDH1B1)1.560.044
Transferrin (TF)1.660.024
Thyroid Stimulating hormone receptor (TSHR)1.700.0036
Lactate Dehydrogenase D (LDHD)1.740.00016
Phosphoenolpyruvate Carboxykinase 2 (PCK2)1.770.0063
ATPase Phospholipid Transporting 8B3 (ATP8B3)1.790.0063
Adrenoceptor Beta 3 (ADRB3)1.790.010
Trefoil Factor 3 (TFF3)1.880.027
ATP Binding Cassette Subfamily C Member 6 (ABCC6)1.890.0051
Phosphoenolpyruvate Carboxykinase 1 (PCK1)2.050.0075
Monoacylglycerol O-Acyltransferase 1 (MOGAT1)2.120.00071
Glycerol-3-Phosphate Acyltransferase 3 (GPAT3)2.160.015
Protein Phosphatase 1 Regulatory Inhibitor Subunit 1B (PPP1R1B)2.380.034
ADAM-Like Decysin 1 (ADAMDEC1)−4.110.00011
Cartilage Oligomeric Matrix Protein (COMP)−3.130.00022
Sterile Alpha Motif Domain Containing 5 (SAMD5)−2.980.000056
Dendrocyte Expressed Seven Transmembrane Protein (DCSTAMP)−2.930.014
R-Spondin 1 (RSPO1)−2.760.0034
Thymic Stromal Lymphopoietin (TSLP)−2.520.00067
Osteoclast Stimulatory Transmembrane Protein (OCSTAMP)−2.460.045
CUB and Sushi Multiple Domains 2 (CSMD2)−2.310.011
ETS Homologous Factor (EHF)−2.280.027
Oxidized Low Density Lipoprotein Receptor 1 (OLR1)−2.230.0063
Carboxypeptidase Z (CPZ)−2.200.029
C-C Motif Chemokine Receptor 1 (CCR1)−2.170.00070
Serpin Family A Member 5 (SERPINA5)−2.110.0020
Pancreatic Lipase Related Protein 1 (PNLIPRP1)−2.100.020
Macrophage Scavenger Receptor 1 (MSR1)−2.080.014
Table
Table 3. Top false discovery rate correction-surviving up- and downregulated transcripts in metabolically unhealthy obese (n = 13) subcutaneous adipose tissue compared to metabolically healthy lean (n = 12).
Table 3. Top false discovery rate correction-surviving up- and downregulated transcripts in metabolically unhealthy obese (n = 13) subcutaneous adipose tissue compared to metabolically healthy lean (n = 12).
TranscriptLog Fold ChangeAdjusted p-Value
Glycogen Synthase 2 (GYS2)3.720.0059
Collagen Type IV Alpha 3 Chain (COL4A3)2.690.0084
Family with Sequence Similarity 131 Member C (FAM131C)2.620.0067
Serum Amyloid A4, Constitutive (SAA4)2.450.028
SH3 Domain Binding Kinase Family Member 2 (SBK2)2.410.0011
Glycine-N-Acyltransferase (GLYAT)2.400.00033
Plasminogen (PLG)2.330.000012
Adrenoceptor Beta 3 (ADRB3)2.280.0019
Stathmin 4 (STMN4)2.230.0456
C-type Lectin Domain Family 4 Member G (CLEC4G)2.220.0060
Thyroid Stimulating Hormone Receptor (TSHR)2.210.000049
Monoacylglycerol O-Acyltransferase 1 (MOGAT1)2.190.00013
Polypeptide N-Acetylgalactosaminyltransferase 9 (GALNT9)2.170.0016
V-Set Furthermore, Transmembrane Domain Containing 2 Like (VSTM2L)2.170.0043
Cat Eye Syndrome Chromosome Region, Candidate 2(CECR2)2.110.0017
Glycoprotein Hormones, Alpha Polypeptide (CGA)−2.620.0054
Osteoclast Stimulatory Transmembrane Protein (OCSTAMP)−2.630.013
WAP Four-Disulfide Core Domain 2 (WFDC2)−2.670.00025
Solute Carrier Family 5 Member 9 (SLC5A9)−2.680.00038
Glutamate Ionotropic Receptor AMPA Type Subunit 1 (GRIA1)−2.910.000011
Cadherin 1 (CDH1)−2.980.0047
Sterile Alpha Motif Domain Containing 5 (SAMD5)−2.990.000010
ETS Homologous Factor (EHF)−3.000.00013
Dendrocyte Expressed Seven Transmembrane Protein (DCSTAMP)−3.020.0040
WD Repeat Domain 64 (WDR64)−3.120.0000085
Matrix Metallopeptidase 7 (MMP7)−3.290.035
Solute Carrier Family 24 Member 2 (SLC24A2)−3.320.00021
S100 Calcium Binding Protein A7A (S100A7A)−3.320.00040
ADAM Like Decysin 1 (ADAMDEC1)−4.060.000028
Death Associated Protein Like 1 (DAPL1)−4.270.000056
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ruggiero, A.D.; Vemuri, R.; Block, M.; DeStephanis, D.; Davis, M.; Chou, J.; Williams, A.; Brock, A.; Das, S.K.; Kavanagh, K. Macrophage Phenotypes and Gene Expression Patterns Are Unique in Naturally Occurring Metabolically Healthy Obesity. Int. J. Mol. Sci. 2022, 23, 12680. https://doi.org/10.3390/ijms232012680

AMA Style

Ruggiero AD, Vemuri R, Block M, DeStephanis D, Davis M, Chou J, Williams A, Brock A, Das SK, Kavanagh K. Macrophage Phenotypes and Gene Expression Patterns Are Unique in Naturally Occurring Metabolically Healthy Obesity. International Journal of Molecular Sciences. 2022; 23(20):12680. https://doi.org/10.3390/ijms232012680

Chicago/Turabian Style

Ruggiero, Alistaire D., Ravichandra Vemuri, Masha Block, Darla DeStephanis, Matthew Davis, Jeff Chou, Abigail Williams, Ashlynn Brock, Swapan Kumar Das, and Kylie Kavanagh. 2022. "Macrophage Phenotypes and Gene Expression Patterns Are Unique in Naturally Occurring Metabolically Healthy Obesity" International Journal of Molecular Sciences 23, no. 20: 12680. https://doi.org/10.3390/ijms232012680

APA Style

Ruggiero, A. D., Vemuri, R., Block, M., DeStephanis, D., Davis, M., Chou, J., Williams, A., Brock, A., Das, S. K., & Kavanagh, K. (2022). Macrophage Phenotypes and Gene Expression Patterns Are Unique in Naturally Occurring Metabolically Healthy Obesity. International Journal of Molecular Sciences, 23(20), 12680. https://doi.org/10.3390/ijms232012680

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