Skip to main content

Principal component analysis of adipose tissue gene expression of lipogenic and adipogenic factors in obesity

Abstract

Objective

A better understanding of mechanisms regulating lipogenesis and adipogenesis is needed to overcome the obesity pandemic. We aimed to study the relationship of the transcript levels of peroxisome proliferator activator receptor γ (PPARγ), CCAAT/enhancer-binding protein alpha (C/EBP-α), liver X receptor (LXR), sterol regulatory element-binding protein-1c (SREBP-1c), fatty acid synthase (FAS), and acetyl-CoA carboxylase (ACC) in subcutaneous adipose tissue (SAT) and visceral adipose tissue (VAT) from obese and normal-weight women with a variety of anthropometric indices, metabolic and biochemical parameters, and insulin resistance.

Methods

Real‐time PCR was done to evaluate the transcript levels of the above‐mentioned genes in VAT and SAT from all participants.

Results

Using principal component analysis (PCA) results, two significant principal components were identified for adipogenic and lipogenic genes in SAT (SPC1 and SPC2) and VAT (VPC1 and VPC2). SPC1 was characterized by relatively high transcript levels of SREBP1c, PPARγ, FAS, and ACC. However, the second pattern (SPC2) was associated with C/EBPα and LXR α mRNA expression. VPC1 was characterized by transcript levels of SREBP1c, FAS, and ACC. However, the VPC2 was characterized by transcript levels of C/EBPα, LXR α, and PPARγ. Pearson’s correlation analysis showed that unlike SPC2, which disclosed an inverse correlation with body mass index, waist and hip circumference, waist to height ratio, visceral adiposity index, HOMA-IR, conicity index, lipid accumulation product, and weight-adjusted waist index, the VPC1 was positively correlated with above-mentioned obesity indices.

Conclusion

This study provided valuable data on multiple patterns for adipogenic and lipogenic genes in adipose tissues in association with a variety of anthropometric indices in obese subjects predicting adipose tissue dysfunction and lipid accumulation.

Peer Review reports

Introduction

Obesity is believed to be a principal factor for the development of insulin resistance, type 2 diabetes, non-alcoholic fatty liver, hypertension, and cardiovascular disorders. Over recent decades, the prevalence of obesity has substantially risen worldwide. In Iran, the prevalence of obesity has increased to about 26% in 2020 and based on a very recent study, the rate of obesity and overweight was reported to be higher among Iranian women than men [1,2,3]. Although white adipose tissue has been considered an active organ secreting several hormones and adipocytokines controlling lipid and glucose metabolism, it is still the principal contributor to the energy reservoir [4].

Adipose tissue dysfunction is a principal contributor to obesity and related metabolic abnormalities which primarily result from impaired lipogenesis and adipogenesis [5,6,7]. At the molecular level, adipogenesis and lipogenesis are regulated by the coordinated and highly complex transcriptional network comprising of peroxisome proliferator activator receptor γ (PPARγ), CCAAT/enhancer-binding protein alpha (C/EBP-α), liver X receptor (LXR), and sterol regulatory element-binding protein 1c (SREBP1c) [8,9,10]. PPARγ, a member of the nuclear-receptor superfamily, has been identified as the master regulator of adipogenesis as well as the main regulator of energy homeostasis in adipose tissue. C/EBPα is a family member of transcription factors. C/EBPs contribute to the regulation of adipocyte differentiation in cooperation with PPARγ and LXRα, a transcription factor belonging to the nuclear receptor superfamily which is involved in a wide range of metabolic pathways, including energy metabolism, inflammation, and particularly adipogenesis and lipogenesis [11,12,13].

Although the results are highly controversial, data from clinical studies indicated that the dysregulation of the aforementioned molecules is linked to obesity and associated metabolic abnormalities, including insulin resistance [14,15,16]. For instance, some studies demonstrated a higher expression of PPARγ and SREBP1c in obese subjects, while others have observed the opposite results in this regard. Moreover, a limited number of data showed a dysregulated expression of LXRα in the context of obesity in humans [14].

Today, a piece of considerable evidence indicated an alteration in the gene expression of the main lipogenic enzymes; fatty acid synthase (FAS) and acetyl-CoA carboxylase (ACC) which are associated with adipocyte lipid accumulation in obesity, although with contradictory findings [10, 17, 18]. For instance, Ranganathan G et al. reported the down-regulation of FAS in adipose tissues of subjects with impaired glucose tolerance [19]. In contrast, another study reported the up-regulation of FAS in adipose tissue of obese subjects which was linked to impaired insulin sensitivity and visceral fat accumulation [20].

Despite the considerable, but conflicting data on the role of dysregulated expression of lipogenic and adipogenic factors in the context of obesity in humans, there is still boundless enthusiasm in the search for this area. It is crucial to mention that the white adipose tissue is mainly distributed in two discrete anatomical depots; subcutaneous adipose tissue (SAT) and visceral adipose tissue (VAT). Remarkably, SAT and VAT are markedly distinct from each other in cellular, molecular, physiologic, clinical, and prognostic characteristics [21, 22].ence, the study of adipose tissue is principal to the understanding of metabolic abnormalities pertinent to the initiation and development of obesity. Moreover, there is still a great deal of work that must be done to unravel the precise role of adipogenic and lipogenic factors in the concept of obesity in humans.

Nowadays, a variety of anthropometric indices have been proposed to monitor and predict metabolic complications in obesity. These include body adiposity index (BAI), The abdominal volume index (AVI), visceral adiposity index (VAI), weight-adjusted waist index (WWI), conicity index (CI), waist-to-hip ratio (WHR), waist circumference (WC), and waist-to-height ratio (WHtR) [23, 24]. Thus far no study has assessed the association of lipogenic and adipogenic factors with a variety of anthropometric indices in obese subjects.

Besides, the regulation of adipose tissue expansion by targeting adipogenesis and adipogenesis has emerged as a possible therapeutic approach to obesity-associated metabolic complications. Therefore, a better understanding of mechanisms regulating lipogenesis and adipogenesis and complicated cross-talk between them is needed to develop appropriate therapeutic strategies to obesity. Noteworthy, there is a sexually dimorphic pattern of adipogenic and lipogenic genes [25, 26]. Moreover, despite increasing concern to women’s health, further studies are still needed to improve their health since obesity is more prevalent in women than men in most countries. This issue highlights an urgent need to develop relevant sex-based therapeutic avenues for obesity-related disorders and obesity. Hence, the current study was conducted to simultaneously measure the mRNA expression of PPARγ, C/EBPα, LXR, SREBP1c, FAS, and ACC in SAT and VAT from obese women with obesity and those with normal-weight. Moreover, this study aimed to evaluate the relationship of the aforementioned genes with a variety of anthropometric indices, metabolic and biochemical parameters, and insulin resistance as well as with each other.

Methods

Obese and normal-weight groups

This case–control study protocol was approved by the ethics committee of Tehran University of Medical Sciences (IR.TUMS.Medicine.REC.1397.702) in compliance with the principles of the Declaration of Helsinki. All methods were performed in accordance with the relevant guidelines and regulations. Written informed consent was obtained from each individual before participation.

A total of 46 subjects (age range 20–53 years) were studied, including 24 obese women (body mass index; BMI ≥ 30 kg/m2) who underwent bariatric surgery and 22 normal-weight women (BMI ≤ 25 kg/m2) who underwent elective surgeries, including inguinal hernia or cholecystectomy. The age range was between 20 and 53 years for normal weight women and between 20 and 48 years for women with obesity. The obese group and normal-weight group were recruited from the Bariatric Surgery Center at Erfan Hospital and from the center of advanced laparoscopic surgeries at Sina and Loqman Hakim hospitals, respectively. All participants were of Iranian ethnicity.

The exclusion criteria for obese and non-obese groups were as follows: 1) had cardiovascular diseases, diabetes, malignancy, acute or chronic inflammatory or infectious diseases, known renal and liver dysfunction; 2) treated with medications for weight loss during the previous 6 months; 3) was currently pregnant or lactating; 4) was menopause; 5) had a surgery history during the previous 6 months, and 6) was currently smoking.

Biochemical and laboratory measurements

The peripheral blood sample was drawn from the antecubital vein by a sterile venipuncture in the morning after an overnight fast before the start of the surgery. Fasting blood samples were centrifuged at 1200 × g for 10 min at 4 °C, and serum was collected and used for measuring fasting blood glucose (FBG), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), triglyceride (TG), total cholesterol (TC), alkaline phosphatase (ALP), aspartate aminotransferase (AST), alanine aminotransferase (ALT), uric acid, urea high sensitivity-reactive protein (hs-CRP), and insulin. The insulin resistance was calculated using the homeostasis model assessment of insulin resistance (HOMA-IR) formula: FBG (mg/dL) × fasting blood insulin (µU/mL) / 405. Resting blood pressure was measured three times on the right arm using a manual sphygmomanometer in seated participants.

Assessment of adiposity indices

BMI was calculated as weight in kg divided by the square of height in meters. In detail, WC was measured midway between the inferior angle of the ribs and the suprailiac crest at the end of a normal expiration. Hip circumference was measured at the maximum circumference over the buttocks using non-stretchable tape. WHR was calculated based on the ratio of WC in cm divided by hip circumference in cm. WHtR was assessed as WC in cm divided by height in cm. WHR and WHtR were also calculated as WC in cm divided by hip circumference in cm and WC in cm divided by height in cm, respectively.

BAI, another index of obesity, was calculated using the following formula:

$$=\frac{hip circumference \left(cm\right)}{{height \left(m\right)}^{1.5}}-18$$

AVI an anthropometric tool for measuring general volume, was computed from the WC and hip circumference using the following equation:

$$=\frac{[2 cm \times {(WC\left(cm\right))}^{2}+0.7cm\times {(WC(cm)-hip(cm))}^{2}]}{1000}$$

Considering the sexual difference in VAT estimations, VAI seems to be the best approach based on a mathematical model by sex, which indirectly expresses visceral adipose tissue function. This index was calculated based on the values of WC, BMI, TG, and HDL-C by using the following formula:

$$Females: VAI=(\frac{wc (cm)}{36.85+\left(1.89\times BMI(kg/{m}^{2}\right))})\times (\frac{TG (\frac{mmol}{l})}{0.81})\times (\frac{1.52}{HDL (\frac{mmol}{l})})$$

WWI is used to assess adiposity by standardizing WC for weight and was calculated as WC in cm divided by the square root of weight in kg (cm/√kg).

CI was calculated based on the values of WC, weight, and height, by using the following formula:

$$\mathrm{CI }=\frac{WC \left(m\right)}{0.109\sqrt{weight (kg)/height (m)}}-18$$

Lipid accumulation product (LAP) is another adiposity index reflecting lipid accumulation was calculated as the sum of WC and TG.

Adipose Tissue Samples, RNA Extraction, and Real-Time Quantitative Polymerase Chain Reaction (PCR)

Paired samples of VAT and SAT were obtained during bariatric.surgery in the obese group or during the inguinal hernia or cholecystectomy in the normal-weight subjects as described previously. Following washing with sterile and cold phosphate-buffered saline, the biopsy samples were immediately frozen in liquid nitrogen and kept at 80 °C for further experiments.

Total RNA was isolated from frozen adipose tissue (100 mg) using the RNeasy Lipid Tissue Mini Kit (Qiagen, Germany).

Frozen adipose tissues (100 mg) were transferred into 1 ml QIAzol Lysis Reagent (Qiagen, Hilden, Germany), and homogenized using a pestle and porcelain mortar. Then, total RNA was isolated from the homogenate using the RNeasy Lipid Tissue Mini Kit (Qiagen, Germany). Before reverse transcription, isolated RNAs were checked for quality and quantity with a NanoDrop 2000 Spectrophotometer (Thermo Fisher Scientific) and RNA integrity was checked by running RNA samples on a 1% agarose gel. Complementary DNA (cDNA) was synthesized.from RNA (1 µg) using the PrimeScript 1st Strand cDNA Synthesis kit (Takara, Japan).

Before further gene expression studies, the standard curve was generated for each primer set using a serial dilution of cDNA synthesized from a pool of RNA extracted from adipose tissues (VAT and SAT). Accordingly, the amplification efficiency for each primer was between 90–100%.

The relative mRNA expression was assessed by real-time PCR using BioFACT™ 2X Real-Time PCR Master Mix (For SYBR Green I) in a Step-One-Plus TM real-time (ABI Applied Biosystems). The sequences of gene-specific primers for the target genes; PPARγ, C/EBPα, LXR, SREBP1c, FAS, and ACC, and the reference gene; β-actin are shown in supplementary Table 1. PCR was performed with 15 min of initial denaturation at 95 °C and then followed by 40 cycles with denaturation at 95 °C for 20 s, and annealing at 60 °C for 15 s. It should be noted that each expression was quantified in duplicate. A melting curve analysis was performed to confirm the specificity of all amplified products.

For each sample, the difference in Ct values (ΔCt) between the target gene and the reference gene was calculated. Since the efficiency (E) of amplification was from 90 to 100% in all assays, we used a 2−∆Ct method to perform relative quantification.

Statistical analysis

Data normality was checked by the Shapiro–Wilk test. Laboratory and anthropometric parameters with normal distribution were presented as mean ± standard deviation (SD), and variables without normal distribution were presented as median (interquartile ranges). All values in the figures were shown as the mean ± standard error of the mean (SEM). Log-transformation was employed for variables with non-normal distribution. The comparison of gene expression levels, as well as anthropometric and biochemical data between obese patients and normal-weight subjects was carried out by independent students t-test on log-transformed variables.

Subsequently, ANCOVA analysis was performed to remove the effects of potential confounders. Correlation coefficients were calculated using the two-tailed Pearson’s correlation analysis. It should be noted that non-normally distributed variables were log-transformed to generate a normal distribution before further analyses. Principal component analysis (PCA) with a varimax rotation was performed to reduce the adipogenic and lipogenic genes into a smaller set of principal components that account for most of the observed variations. It should be noted that the patterns were normalized using the Kaiser method to have greater interpretability. Using PCA results, two significant principal components for adipogenic and lipogenic genes in SAT (SPC1 and SPC2) and VAT (VPC1 and VPC2) were identified as shown in supplementary Table 2 and supplementary Table 3, respectively.

All statistical assessments were two-tailed and P-value < 0.05 was considered statistically significant. All data analysis was performed using SPSS 20 (SPSS, Chicago, IL, USA).

Results

Study population characteristics

Some parts of the results regarding the anthropometric, clinical, and metabolic characterizations of the study population which are presented here are related to those in earlier work (unpublished data). The anthropometric and metabolic variables of obese and normal-weight groups are displayed in Table 1. There is no statistically significant difference in terms of age, SBP, DBP, and circulating levels of FBG, urea, HDL-C, TG, AST, and ALT between the two studied groups. However, the circulating levels of LDL-C, total cholesterol, uric acid, creatinine, albumin, total protein, hs-CRP, insulin, HbA1c, and HOMA-IR value were significantly higher in the obese group compared to normal-weight subjects.

Table 1 The anthropometric, clinical, and metabolic characterizations of all participants

As for adiposity indices, the anthropometric measurements, including BMI, WC, hip circumference, BAI, AVI, CI, and WHtR were significantly higher in the obese group in comparison with the non-obese group.

Notably, remaining obesity indices including WHR (p = 0.088) and WWI (p = 0.073) showed an increasing trend albeit with a borderline statistical significance in women with obesity compared to ones with normal weight.

The mRNA Expression of PPARγ, C/EBPα, LXR, SREBP1c, FAS, and ACC in SAT and VAT from Obese and Normal-Weight Individuals

The results of the mRNA expression of PPARγ, C/EBPα, LXR, SREBP1c, FAS, and ACC in SAT and VAT from the two studied groups are shown in Fig. 1.

Fig. 1
figure 1

Transcript levels of PPARγ (a), LXR (b),C/EBPα (c), SREBP1c (d), FAS (e), and ACC (b) genes in the visceral (VAT) and subcutaneous (SAT) adipose tissues of obese (O) subjects (n = 24) and normal-weight (N) individuals (n = 22). Results were shown as mean ± standard error of mean (SEM)

The mRNA expression of PPARγ was significantly lower in VAT (p = 0.033) of obese women compared to normal-weight subjects (Fig. 1a). Moreover, we found that the mRNA expression of LXRα was significantly lower in both VAT (p = 0.009) and SAT (p = 0.012) from the obese group compared to the non-obese group (Fig. 1b). Although the expression of C/EBPα in the VAT (p = 0.75) was not significantly different between the two studied groups, it was lower in SAT (p = 0.003) in obese subjects compared to the non-obese group (Fig. 1c).

Conversely, our results revealed higher transcript levels of SREBP1c and its target genes; FAS and ACC in the VAT from obese subjects compared to normal-weight ones (p = 0.050, p = 0.017, and p = 0.02, respectively), whereas the mRNA expression of these genes in the SAT was not significantly different between the two groups (p = 0.725 for SREBP1c, p = 0.598 for FAS, and p = 0.742 for ACC) (Fig. 1d-e). Therefore, it is likely that the expression of PPARγ, C/EBPα, and LXR was lower while the expression of SREBP1c, FAS, and ACC was higher in the obese group compared to the non-obese counterparts.

ANCOVA was performed to remove the possible effect of age and HOMA-IR on the gene expression of all studied genes in VAT and SAT. The results indicated that the decrease in the expression of LXRα both in VAT (p = 0.004) and SAT (p = 0.002) of the obese group was independent of age and HOMA-IR. Moreover, the C/EBPα transcript level (P = 0.02) in the SAT of the obese group remained significantly different from that of the normal-weight group after adjustment for age and HOMA-IR. However, the difference in PPARγ (P = 0.131), SREBP1c (P = 0.46), FAS (P = 0.24), and ACC (P = 0.14) in VAT did not remain statistically significant after correction for age and HOMA-IR.

Univariate Correlations of mRNA Expression of All Studied Genes with Adiposity Indices and Metabolic Profile

We also assessed the Pearson correlation coefficient of PPARγ, C/EBPα, LXR, SREBP1c, FAS, and ACC mRNA levels with adiposity indices and metabolic profiles in SAT (Table 2) and VAT (Table 3) of the whole population.

Table 2 Pearson correlation of PPARγ, C/EBPα, LXR, SREBP1c, FAS, and ACC mRNA levels in subcutaneous adipose tissue of whole participants with adiposity indices and metabolic profile
Table 3 Pearson correlation of PPARγ, C/EBPα, LXR, SREBP1c, FAS, and ACC mRNA levels in visceral adipose tissue of whole participants with adiposity indices and metabolic profile

Our results indicated that the mRNA expression of PPARγ in the SAT is positively correlated with insulin levels (r = 0.365, p = 0.013).

Moreover, SAT mRNA expression of LXR showed an inverse correlation with adiposity indices including BMI (r = -0.330, p = 0.025), WC (r = -0.413, p = 0.004), HC (r = -0.376, p = 0.010), BAI (r = -0.380, p = 0.009), CI (r = -0.457, p = 0.001), AVI (r = -0.413, p = 0.004), WHtR (r = -0.467, p = 0.001), WWI (r = -0.317, p = 0.032), and LAP (r = -0.354, p = 0.016).

Similarly, SAT mRNA expression of C/EBPα was inversely correlated with BMI (r = -0.410, p = 0.005), WC (r = -0.410,p = 0.005), HC (r = -0.462, p = 0.001), BAI (r = -0.544, p = 0.0001), CI (r = -0.421, p = 0.004), AVI (r = -0.412, p = 0.004), WHtR (r = -0.425, p = 0.003), WWI (r = -0.374, p = 0.010), and LAP (r = -0.345, p = 0.019), and HOMA-IR (r = -0.373, p = 0.011).

In this study, according to correlation analysis, mRNA expression of PPARγ in VAT showed a significant inverse relationship with the BMI (r = -0.299, p = 0.044), VAI (r = -0.355, p = 0.015), LAP (r = -0.367, p = 0.012), HC (r = -0.368, p = 0.012), and BAI (r = -0.361, p = 0.014) in the whole study population. Our results also revealed a significant inverse correlation of VAT LXRα transcript levels with LAP (r = -0.288, p = 0.043), BAI (r = -0.329, p = 0.026), and WHtR (r = -0.297, p = 0.045). Moreover, there was an inverse correlation between VAT C/EBPα transcript levels and VAI (r = -0.336, p = 0.023).

On the other hand, VAT SREBP1c mRNA expression showed a significant positive correlation with BMI (r = 0.325, p = 0.027), WC (r = 0.346, p = 0.019), HC (r = 0.384, p = 0.008), WHtR (r = 0.328, p = 0.026), BAI (r = 0.333, p = 0.024), AVI (r = 0.349, p = 0.017), CI (r = 0.293, p = 0.048), and HOMA-IR (r = 0.349, p = 0.017).

FAS gene expression in the VAT also correlated positively with LAP (r = 0.424, p = 0.003), BMI (r = 0.387, p = 0.008), WC (r = 0.496, p = 0.0001), HC (r = 0.435, p = 0.002), WHR (r = 0.322, p = 0.029), BAI (r = 0.350, p = 0.017),CI(r = 0.473, p = 0.001), AVI (r = 0.494, p = 0.0001), WHtR (r = 0.438, p = 0.002), and HOMA-IR (r = 0.305, p = 0.039). We also found a positive correlation between ACC gene expression in VAT and BMI (r = 0.344, p = 0.019), and insulin levels (r = 0.296, p = 0.046).

Altogether, it seems that the expression level of SREBP1c, FAS, and ACC had a positive correlation with obesity indices as well as HOMA-IR while the expression level of PPARγ, C/EBPα, and LXR had an inverse correlation with adiposity indices in all subjects.

Correlation Analysis of the Transcript Level Patterns of Adipogenic and Lipogenic Genesin SAT and VAT of Whole Participants with Adiposity Indices and Metabolic Profile

Using PCA results, we identified two significant principal components (SPC1 and SPC2) for adipogenic and lipogenic genes in SAT. SPC1 was characterized by relatively high transcript levels of SREBP1c, PPARγ, FAS, and ACC. However, the second pattern (SPC2) was associated with C/EBPα and LXR α mRNA expression.

Similarly, two patterns for adipogenic and lipogenic genes were identified (VPC1 and VPC2) in VAT. Specifically, the first pattern was characterized by transcript levels of SREBP1c, FAS, and ACC. The second one was characterized by C/EBPα, LXR α, and PPARγ mRNA expression.

Next, Pearson’s correlation analysis was used to assess the possible association of each lipogenic and adipogenic pattern score with obesity indices and insulin resistance indices (Table 4).

Table 4 Pearson correlation of pattern of the transcript level of adipogenic and lipogenic genes in subcutaneous and visceral adipose tissue of whole participants with adiposity indices and metabolic profile

As shown in Table 4, SPC1 was positively correlated with insulin level (r = 0.301, p = 0.042) in the whole study population. However, SPC2 was inversely correlated with BMI (r = -0.492, p = 0.001), WC (r = -0.530,p = 0.0001), HC(r = -0.549, p = 0.0001), WHtR (r = -0.579, p = 0.0001), VAI (r = -0.445, p = 0.002), LAP (r = -0.615, p = 0.0001), AVI (r =—0.429, p = 0.003), WWI (r =—0.554, p = 0.0001), CI (r = -0.531, p = 0.0001), HOMA-.IR (r = -0.374, p = 0.011) in the whole study population.

The VPC1 was significantly correlated with BMI (r = 0.398, p = 0.006), WC (r = 0.433,p = 0.003), HC(r = 0.440, p = 0.002), WHtR (r = 0.396, p = 0.012), VAI (r = 0.355, p = 0.015), LAP (r = 0.3035, p = 0.023), WWI (r = 0.329, p = 0.026), CI (r = 0.434, p = 0.003), insulin level (r = 0.387, p = 0.008), and HOMA-IR (r = 0.306, p = 0.038) in the whole study population. However, BAI (r = -0.380, p = 0.009), VAI (r = -0.311, p = 0.035), and LAP (r = -0.305, p = 0.039) showed an inverse correlation with VPC2.

Remarkably, the observed results based on PCA data was almost in parallel to individual expression levels of adipogenic and lipogenic factors and their correlation with inulin value and adiposity indices. Indeed, the association of the identified patterns, especially SPC2 and VPC1 with adiposity parameters as well as insulin resistance demonstrated significant p-values.

Discussion

The alteration in the gene expression of lipogenic and adipogenic factors in relation with obesity and associated metabolic abnormalities has been reported in numerous pieces of evidence although with inconsistent results. Moreover, this association has been mostly investigated individually [27, 28]. Therefore, the present study aimed to investigate the association of the gene expression of lipogenic and adipogenic factors pattern with insulin resistance and a variety of anthropometric indices in the context of obesity in humans. In particular, here, we explain it.from two perspectives; (i) individually and [2] PCA-derived findings.

According to the first perspective, we found a reduced transcript level of PPARγ, C/EBPα, and LXRwhile the expression level of SREBP1c, FAS, and ACC was higher in adipose depots of subjects with obesity compared to the non-obese counterparts. This association was further substantiated by correlation analysis, as revealed that transcript levels of SREBP1c, FAS, and ACC had a positive correlation with a variety of adiposity indices and HOMA-IR values. However, the transcript level of PPARγ, C/EBPα, and LXR displayed an inverse correlation with the indicated parameters. Of great interest, the involvement of adipogenic and lipogenic has been reported in the context of obesity, albeit with contrary results. For instance, some studies reported a higher expression of PPARγ and SREBP1c in obese subjects [29, 30], while others observed the opposite results in this regard [15, 31]. Also, a limited number of data showed a dysregulated expression of LXRα in the context of obesity in humans [14]. In contrast to our findings, Ranganathan G et al. reported the down-regulation of FAS in adipose tissues of subjects with impaired glucose tolerance [19]. Nevertheless, another study reported the up-regulation of FAS in adipose tissue of obese subjects which was linked to impaired insulin sensitivity and visceral fat accumulation [20].

According to the second perspective based PCA-derived data, two separate patterns for adipogenic and lipogenic genes in SAT and VAT were obtained. In SAT, the first pattern was characterized by relatively high transcript levels of SREBP1c, PPARγ, FAS, and ACC. However, the second pattern was associated with C/EBPα and LXR α mRNA expression. In VAT, one pattern included transcript levels of SREBP1c, FAS, and ACC while, the second one was linked to C/EBPα, LXR α, and PPARγ mRNA expression. Interestingly, the link between the identified patterns, especially the second pattern in SAT and the first pattern in VAT and adiposity parameters, and insulin resistance yielded significant p-values.

To further explain the second perspective, it is crucial to mention that adipose tissue dysfunction is a principal contributor to obesity and related metabolic abnormalities which primarily results from impaired transcriptional regulation of the key factors that control lipogenesis and adipogenesis. A complex and finely well-organized transcriptional network (cascade) controls adipogenesis and lipogenesis. Hence, it seems that the evaluation of a combination of adipogenic and lipogenic factors could, at least partly, clarify obscure issues in the pathogenesis of obesity [27]. In this way, PCA is not only used as a powerful method to recognize the clustering of adipogenic and lipogenic genes in the context of human obesity, but is also applied to the data reduction of the existing 6 orthogonal variables. Subsequently, linear regression was performed to provide detailed information about quantitative associations between the PCA-obtained pattern and adiposity indices.

Induction of adipocyte differentiation is mediated by a combination of adipogenic transcription factors; PPARγ, CEBPα, and LXRα, and the subsequent lipid storage by a lipogenic transcription factor; SREBP1c [29, 32]. Adipogenesis is the process of pre-adipocyte differentiation to lipid-laden adipocytes and has a central role in systemic energy homeostasis. However, lipogenesis is a process that happens preferentially in adipose tissue and encompasses the synthesis of fatty acids and triglycerides as energy reserves. A well-functioning adipose tissue can sequester nutritional overload in the form of triglycerides and therefore have a protective role against detrimental effects of lipid storage in other organs [8]. The expression of adipogenic genes usually shows that how the adipose tissue expands [33]. In line with this, a decreased expression of PPARγ, CEBPα, and LXRα in adipose tissue of women with obesity can reflect low adipose tissue expansion in the present study. Remarkably, previous evidence points out that C/EBPα and PPARγ are principal players in adipogenesis and participate in a single pathway of adipocyte formation with PPARγ being the master regulator of adipogenesis [34]. Remarkably, the failure of fat accumulation was reported before in CEBPα knockout mice immediately after birth. Moreover, male PPARγ knockout mice were more insulin resistant on a normal diet in comparison with the wild-type ones. Indeed, this study indicated that impaired adipogenesis could trigger metabolic abnormalities, such as increased adiposity markers and insulin resistance in obese cases [35]. It was also suggested that adipose tissue-specific ablation of PPARγ causes progressive lipodystrophy [36]. Intriguingly, VAI and LAP have been proposed as substitute indicators of dysfunction and distribution of adipose tissue that indirectly predict cardiometabolic risk [37, 38]. Here, the SPC2 pattern score characterized by C/EBPα and LXRα expression was also associated with VAI and LAP expressing lipid accumulation and adipose tissue dysfunction. It is important to note that.an impaired fat deposition in adipose tissue is closely linked to lipotoxicity and subsequently lipid accumulation in the liver, skeletal muscle, and pancreas, which in turn lead to the development of insulin resistance and other abnormalities pertinent to obesity [39].

This concept was strengthened by increased expression of lipogenic genes; SREBP1c, FAS, and ACC too. Indeed,

SREBP1c which promotes lipogenesis by inducing lipogenic enzymes; FAS and ACC and in turn, leads to increased lipid accumulation within the adipocyte [40] was reported to be overexpressed in obese subjects compared to the non-obese ones. To support this notion, VPC1 reflecting high transcript levels of SREBP1c, FAS, and ACC was significantly correlated with LAP and VAI as well.

Here, we identified several potential patterns of adipogenic and lipogenic genes in adipose tissues using PCA for future investigations. However, it should be noted that this research area needs a more comprehensive study considering other factors in adipose tissue function. Also, it should be noted that the measurement of adiponectin transcript level, as a major component of adiposity can provide valuable information on in-depth interpretation of our results. Finally, we would like to stress that further studies on both sexes and with a larger sample size are required to strengthen the results.

Conclusion

The present study used PCA to simultaneously evaluate the expression levels of the key genes involved in the complex transcriptional network regulating adipogenesis and lipogenesis; comprising of PPARγ, C/EBPα, LXR, SREBP1c, FAS, and ACC in SAT and VAT from obese and normal-weight subjects to decrease the adipogenic and lipogenic genes into a smaller set of principal components that account for most of the observed variations. Interestingly, this study provided valuable information on multiple patterns for adipogenic and lipogenic genes in adipose tissues with insulin resistance and a variety of anthropometric indices predicting adipose tissue dysfunction and lipid accumulation in the context of obesity.

Availability of data and material

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

References

  1. Pourfarzi F, Sadjadi A, Poustchi H, Amani F. Prevalence of overweight and obesity in Iranian population: A population-based study in northwestern of Iran. J Public Health Res. 2021;11(1).

  2. Okati-Aliabad H, Ansari-Moghaddam A, Kargar S, Jabbari N. Prevalence of Obesity and Overweight among Adults in the Middle East Countries from 2000 to 2020: A Systematic Review and Meta-Analysis. J Obes. 2022;2022:8074837.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Esteghamati A, Khalilzadeh O, Mohammad K, Meysamie A, Rashidi A, Kamgar M, et al. Secular trends of obesity in Iran between 1999 and 2007: National Surveys of Risk Factors of Non-communicable Diseases. Metab Syndr Relat Disord. 2010;8(3):209–13.

    Article  PubMed  Google Scholar 

  4. Longo M, Zatterale F, Naderi J, Parrillo L, Formisano P, Raciti GA, et al. Adipose Tissue Dysfunction as Determinant of Obesity-Associated Metabolic Complications. International journal of molecular sciences. 2019;20(9).

  5. Smith U, Kahn BB. Adipose tissue regulates insulin sensitivity: role of adipogenesis, de novo lipogenesis and novel lipids. J Intern Med. 2016;280(5):465–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Luo L, Liu M. Adipose tissue in control of metabolism. J Endocrinol. 2016;231(3):R77-r99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. White U, Ravussin E. Dynamics of adipose tissue turnover in human metabolic health and disease. Diabetologia. 2019;62(1):17–23.

    Article  PubMed  Google Scholar 

  8. Bjorndal B, Burri L, Staalesen V, Skorve J, Berge RK. Different adipose depots: their role in the development of metabolic syndrome and mitochondrial response to hypolipidemic agents. Journal of obesity. 2011;2011: 490650.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Feng B, Zhang T, Xu H. Human adipose dynamics and metabolic health. Ann N Y Acad Sci. 2013;1281:160–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Lenhard JM. Lipogenic enzymes as therapeutic targets for obesity and diabetes. Curr Pharm Des. 2011;17(4):325–31.

    Article  CAS  PubMed  Google Scholar 

  11. Kersten S. Mechanisms of nutritional and hormonal regulation of lipogenesis. EMBO Rep. 2001;2(4):282–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Laurencikiene J, Rydén M. Liver X receptors and fat cell metabolism. International journal of obesity (2005). 2012;36(12):1494–502.

  13. Moreno-Indias I, Tinahones FJ. Impaired adipose tissue expandability and lipogenic capacities as ones of the main causes of metabolic disorders. J Diabetes Res. 2015;2015: 970375.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Dahlman I, Nilsson M, Jiao H, Hoffstedt J, Lindgren CM, Humphreys K, et al. Liver X receptor gene polymorphisms and adipose tissue expression levels in obesity. Pharmacogenet Genomics. 2006;16(12):881–9.

    Article  CAS  PubMed  Google Scholar 

  15. Kolehmainen M, Vidal H, Alhava E, Uusitupa MI. Sterol regulatory element binding protein 1c (SREBP-1c) expression in human obesity. Obes Res. 2001;9(11):706–12.

    Article  CAS  PubMed  Google Scholar 

  16. Giusti V, Verdumo C, Suter M, Gaillard RC, Burckhardt P, Pralong F. Expression of peroxisome proliferator-activated receptor-gamma1 and peroxisome proliferator-activated receptor-gamma2 in visceral and subcutaneous adipose tissue of obese women. Diabetes. 2003;52(7):1673–6.

    Article  CAS  PubMed  Google Scholar 

  17. Kusunoki J, Kanatani A, Moller DE. Modulation of fatty acid metabolism as a potential approach to the treatment of obesity and the metabolic syndrome. Endocrine. 2006;29(1):91–100.

    Article  CAS  PubMed  Google Scholar 

  18. Menendez JA, Vazquez-Martin A, Ortega FJ, Fernandez-Real JM. Fatty acid synthase: association with insulin resistance, type 2 diabetes, and cancer. Clin Chem. 2009;55(3):425–38.

    Article  CAS  PubMed  Google Scholar 

  19. Ranganathan G, Unal R, Pokrovskaya I, Yao-Borengasser A, Phanavanh B, Lecka-Czernik B, et al. The lipogenic enzymes DGAT1, FAS, and LPL in adipose tissue: effects of obesity, insulin resistance, and TZD treatment. J Lipid Res. 2006;47(11):2444–50.

    Article  CAS  PubMed  Google Scholar 

  20. Berndt J, Kovacs P, Ruschke K, Kloting N, Fasshauer M, Schon MR, et al. Fatty acid synthase gene expression in human adipose tissue: association with obesity and type 2 diabetes. Diabetologia. 2007;50(7):1472–80.

    Article  CAS  PubMed  Google Scholar 

  21. Lee MJ, Wu Y, Fried SK. Adipose tissue heterogeneity: implication of depot differences in adipose tissue for obesity complications. Mol Aspects Med. 2013;34(1):1–11.

    Article  CAS  PubMed  Google Scholar 

  22. Guglielmi V, Sbraccia P. Obesity phenotypes: depot-differences in adipose tissue and their clinical implications. Eating and weight disorders : EWD. 2018;23(1):3–14.

    Article  PubMed  Google Scholar 

  23. de Oliveira CC, Roriz AK, Ramos LB, Gomes NM. Indicators of Adiposity Predictors of Metabolic Syndrome in the Elderly. Ann Nutr Metab. 2017;70(1):9–15.

    Article  PubMed  Google Scholar 

  24. Gateva A, Assyov Y, Kamenov Z. Usefulness of different adiposity indexes for identification of metabolic disturbances in patients with obesity. Archives of physiology and biochemistry. 2021:1–6.

  25. Moradi N, Fadaei R, Emamgholipour S, Kazemian E, Panahi G, Vahedi S, et al. Association of circulating CTRP9 with soluble adhesion molecules and inflammatory markers in patients with type 2 diabetes mellitus and coronary artery disease. PloS one. 2018;13(1):e0192159-e.

  26. Chang E, Varghese M, Singer K. Gender and Sex Differences in Adipose Tissue. Curr Diab Rep. 2018;18(9):69.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Longo M, Zatterale F, Naderi J, Parrillo L, Formisano P, Raciti GA, et al. Adipose tissue dysfunction as determinant of obesity-associated metabolic complications. Int J Mol Sci. 2019;20(9):2358.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Halliwell KJ, Fielding BA, Samra JS, Humphreys SM, Frayn KN. Release of individual fatty acids from human adipose tissue in vivo after an overnight fast. J Lipid Res. 1996;37(9):1842–8.

    Article  CAS  PubMed  Google Scholar 

  29. Crewe C, Zhu Y, Paschoal VA, Joffin N, Ghaben AL, Gordillo R, et al. SREBP-regulated adipocyte lipogenesis is dependent on substrate availability and redox modulation of mTORC1. JCI Insight. 2019;5(15).

  30. Ebrahimi R, Toolabi K, Jannat Ali Pour N, Mohassel Azadi S, Bahiraee A, Zamani-Garmsiri F, et al. Adipose tissue gene expression of long non-coding RNAs; MALAT1, TUG1 in obesity: is it associated with metabolic profile and lipid homeostasis-related genes expression? Diabetol Metab Syndr. 2020;12:36.

  31. Torres J-L, Usategui-Martín R, Hernández-Cosido L, Bernardo E, Manzanedo-Bueno L, Hernández-García I, et al. PPAR-γ Gene Expression in Human Adipose Tissue Is Associated with Weight Loss After Sleeve Gastrectomy. J Gastrointest Surg. 2022;26(2):286–97.

    Article  PubMed  Google Scholar 

  32. Madsen MS, Siersbæk R, Boergesen M, Nielsen R, Mandrup S. Peroxisome proliferator-activated receptor γ and C/EBPα synergistically activate key metabolic adipocyte genes by assisted loading. Mol Cell Biol. 2014;34(6):939–54.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Rodríguez-Acebes S, Palacios N, Botella-Carretero JI, Olea N, Crespo L, Peromingo R, et al. Gene expression profiling of subcutaneous adipose tissue in morbid obesity using a focused microarray: Distinct expression of cell-cycle- and differentiation-related genes. BMC Med Genomics. 2010;3(1):61.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Rosen ED, Hsu CH, Wang X, Sakai S, Freeman MW, Gonzalez FJ, et al. C/EBPalpha induces adipogenesis through PPARgamma: a unified pathway. Genes Dev. 2002;16(1):22–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Rosen ED, Walkey CJ, Puigserver P, Spiegelman BM. Transcriptional regulation of adipogenesis. Genes Dev. 2000;14(11):1293–307.

    Article  CAS  PubMed  Google Scholar 

  36. Wang F, Mullican SE, DiSpirito JR, Peed LC, Lazar MA. Lipoatrophy and severe metabolic disturbance in mice with fat-specific deletion of PPARγ. Proc Natl Acad Sci U S A. 2013;110(46):18656–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Du T, Yu X, Zhang J, Sun X. Lipid accumulation product and visceral adiposity index are effective markers for identifying the metabolically obese normal-weight phenotype. Acta Diabetol. 2015;52(5):855–63.

    Article  CAS  PubMed  Google Scholar 

  38. Ahn N, Baumeister SE, Amann U, Rathmann W, Peters A, Huth C, et al. Visceral adiposity index (VAI), lipid accumulation product (LAP), and product of triglycerides and glucose (TyG) to discriminate prediabetes and diabetes. Sci Rep. 2019;9(1):9693.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Ahmed B, Sultana R, Greene MW. Adipose tissue and insulin resistance in obese. Biomed Pharmacother. 2021;137: 111315.

    Article  CAS  PubMed  Google Scholar 

  40. Boizard M, Le Liepvre X, Lemarchand P, Foufelle F, Ferré P, Dugail I. Obesity-related Overexpression of Fatty-acid Synthase Gene in Adipose Tissue Involves Sterol Regulatory Element-binding Protein Transcription Factors*. J Biol Chem. 1998;273(44):29164–71.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This study was supported by Tehran University of Medical Sciences, Tehran, Iran (grant number: 97.01-30-37421).

Funding

This study was supported by Tehran University of Medical Sciences, Tehran, Iran (grant number: 97.01–30-37421).

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization of the project: S.E, Laboratory procedures: N.J, R.E, Data curation: R.M, To collect tissues and clinical data: H.Z, Data analysis: S.E, M.Y, Interpretation of data: M.Y, S.E, Writing an original draft: S.E, N.A, SM.H, Reviewing & editing: S.E, R.M All authors read and approved the manuscript.

Corresponding author

Correspondence to Solaleh Emamgholipour.

Ethics declarations

Ethics approval and consent to participate

This case–control study protocol was approved by the ethics committee of Tehran University of Medical Sciences (IR.TUMS.Medicine.REC.1397.702) in compliance with the principles of the Declaration of Helsinki. All methods were performed in accordance with the relevant guidelines and regulations. Written informed consent was obtained from each individual before participation.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Forward and reverse primers used for real-time PCR. Table S2. Principal factor loading of transcript levels of adipogenic and lipogenic genes Table S3. Principal factor loading of transcript levels of adipogenic and lipogenic genes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jannat Ali Pour, N., Zabihi-Mahmoudabadi, H., Ebrahimi, R. et al. Principal component analysis of adipose tissue gene expression of lipogenic and adipogenic factors in obesity. BMC Endocr Disord 23, 94 (2023). https://doi.org/10.1186/s12902-023-01347-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12902-023-01347-w

Keywords