Identification of <i>SCARA3</i> with potential roles in metabolic disorders

Hui Peng; Qi Guo; Tian Su; Ye Xiao; Chang-Jun Li; Yan Huang; Xiang-Hang Luo

doi:10.18632/aging.202228

Research Paper Volume 13, Issue 2 pp 2149—2167

Identification of SCARA3 with potential roles in metabolic disorders

Hui Peng^1, , Qi Guo^1, , Tian Su^1, , Ye Xiao^1, , Chang-Jun Li^1, , Yan Huang^1, , Xiang-Hang Luo^1, ,

¹ Department of Endocrinology, Endocrinology Research Center, Xiangya Hospital of Central South University, Changsha, China

Received: July 22, 2020 Accepted: October 22, 2020 Published: December 9, 2020

https://doi.org/10.18632/aging.202228
How to Cite

Copyright: © 2020 Peng et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Obesity is characterized by the expansion of adipose tissue which is partially modulated by adipogenesis. In the present study, we identified five differentially expressed genes by incorporating two adipogenesis-related datasets from the GEO database and their correlation with adipogenic markers. However, the role of scavenger receptor class A member 3 (SCARA3) in obesity-related disorders has been rarely reported. We found that Scara3 expression in old adipose tissue-derived mesenchymal stem cells (Ad-MSCs) was lower than it in young Ad-MSCs. Obese mice caused by deletion of the leptin receptor gene (db/db) or by a high-fat diet both showed reduced Scara3 expression in inguinal white adipose tissue. Moreover, hypermethylation of SCARA3 was observed in patients with type 2 diabetes and atherosclerosis. Data from the CTD database indicated that SCARA3 is a potential target for metabolic diseases. Mechanistically, JUN was predicted as a transcriptional factor of SCARA3 in different databases which is consistent with our further bioinformatics analysis. Collectively, our study suggested that SCARA3 is potentially associated with age-related metabolic dysfunction, which provided new insights into the pathogenesis and treatment of obesity as well as other obesity-associated metabolic complications.

Introduction

Obesity, especially metabolically unhealthy obesity, has been reported to increase the level of glucose and lipids in circulation, resulting in various metabolic disorders, such as cardiovascular diseases, fatty liver and diabetes mellitus [1–4]. With aging populations, the prevalence of obesity has increased significantly worldwide over the past five decades [5–9]. Thus, obesity imposes enormous economic and heath burdens and has become a marked challenge for individuals and public health.

Obesity is characterized by the excessive accumulation of adipose tissue and is associated with adipose oxidative stress [10, 11]. The expansion of adipose tissue is modulated via two distinct mechanisms [1]. On the one hand, adipocyte hyperplasia is characterized by the increase of number of adipocytes mediated by adipogenic differentiation. On the other hand, adipocyte hypertrophy is characterized by adipocyte enlargement [8]. To date, many studies have shown that adipogenic differentiation and subsequent adipocytes maturation play crucial roles in the etiology of obesity, as well as its metabolic complications [8, 12–15]. The expansion of white adipose tissues caused by excessive adipogenesis is found in mice and humans [16–18]. In addition, adipose tissue expansion either in mice induced by a high-fat diet or in leptin receptor-deficient db/db mice, was accompanied by increased adipogenic differentiation [19–21]. Mesenchymal stem cells (MSCs) are multipotent cells with the ability to differentiate into mature cells of several mesenchymal tissues, such as fat and bone [22, 23]. Previously, we identified several factors associated with the lineage commitment of MSCs that showed important impact on metabolic disorders during aging [6, 24–26]. During aging, age-related factors generally contribute to the accumulation of adipose tissue or metabolic dysfunction [6, 26–28]. However, the pathogenesis of obesity during aging mediated by MSCs-associated adipogenesis remains unclear.

In this study, we validated several potential genes (GHR (growth hormone receptor), GPX3 (glutathione peroxidase 3), SAA1 (serum amyloid A1), SCARA3 (scavenger receptor class A member 3), and WFDC1 (WAP four-disulfide core domain 1)) during adipogenesis through the selection of differentially expressed genes (DEGs) and weighted gene co-expression network analysis (WGCNA). Mechanistically, the downregulation of SCARA3 mediated by methylation is highly likely to promote adipogenesis by regulating JUN (Jun proto-oncogene, AP-1 transcription factor subunit) involved transcription factors and peroxisome proliferator activated receptor (PPAR) signaling. Furthermore, SCARA3 was associated with metabolic dysfunction in both humans and animals in an age-dependent manner. Thus, the results of the present study suggested that SCARA3 is a potential new target for diagnosis and treatment of obesity, as well as other obesity-associated metabolic complications.

Results

Identification of the key modules through WGCNA

The workflow of the present study, comprising various bioinformatic data analyses and further validation in vivo is shown in Figure 1. To identify the key modules closely associated with adipogenesis, we firstly utilized WGCNA to analyze all genes from the adipogenesis-related dataset (GEO100748), which comprised three different stages of adipogenesis in mesenchymal stem cells (MSCs), such as the undifferentiated stage, day 7 of the differentiated stage and day 21 of the differentiated stage [29]. The sample correlation showed there was no obvious batch effect (Figure 2A, 2B). This dataset has 37 samples, and we selected the soft-thresholding power as 14, according to the instructions of WGCNA (https://horvath.genetics.ucla.edu/html/CoexpressionNetwork/Rpackages/WGCNA/) (Figure 2C). Twenty-four modules were identified based on the average linkage hierarchical clustering and the soft-thresholding power (Figure 2D). The top four modules with most close association with adipogenic stages were grey 60, turquoise, green and red modules (Figure 2E). Thus, we selected these top four modules for further analysis. The four modules contain 276, 1459, 853 and 835 genes, respectively (Supplementary Table 1–4). The grey 60 and turquoise modules were highly positively correlated with the stages of adipogenesis, while the modules of green and red are highly negatively correlated with the stages of adipogenesis (Grey 60: r = -0.74, p = 4.2e-49; Turquoise: r = -0.9, p=1e-200; Green: r = 0.75, p = 6.3e-155; Red: r = 0.99, p = 1e-200) (Figure 3A–3D). The enrichments for the genes in each module were evaluated through GO analysis (Figure 3E–3H).

Figure 1. Workflow of the present study.

Identification of key modules that correlate with adipogenesis in the GEO dataset through WGCNA. (A) Sample dendrogram and trait heatmap. (B) Sample correlations among three different stages. (C) Analysis of the scale-free fit index (left) and the mean connectivity (right) for various soft-thresholding powers based on a scale-free R2 (R2 = 0.686, power = 14). (D) Dendrogram of all genes clustered based on a dissimilarity measure (1-TOM). Each branch in the dendrogram represents one gene and each module color represents one co-expression module. (E) Heatmap of the correlation between epigengene module and traits of adipogenesis. Each group contains the correlation coefficient and P value. The digits in the brackets on the left side represent the number of genes in the corresponding epigengene module.

Figure 2. Identification of key modules that correlate with adipogenesis in the GEO dataset through WGCNA. (A) Sample dendrogram and trait heatmap. (B) Sample correlations among three different stages. (C) Analysis of the scale-free fit index (left) and the mean connectivity (right) for various soft-thresholding powers based on a scale-free R² (R² = 0.686, power = 14). (D) Dendrogram of all genes clustered based on a dissimilarity measure (1-TOM). Each branch in the dendrogram represents one gene and each module color represents one co-expression module. (E) Heatmap of the correlation between epigengene module and traits of adipogenesis. Each group contains the correlation coefficient and P value. The digits in the brackets on the left side represent the number of genes in the corresponding epigengene module.

Identification of functional annotation of the WGCNA module that correlated highly with adipogenesis. (A–D) Scatter plot of module Eigengenes in the grey 60 module (A), turquoise module (B), green module (C), and red module (D). (E–H) Biological process GO terms for genes in the grey 60 module (E), turquoise module (F), green module (G), and red module (H). P-value cutoff = 0.01; q-value cutoff = 0.05. GO, Gene Ontology.

Figure 3. Identification of functional annotation of the WGCNA module that correlated highly with adipogenesis. (A–D) Scatter plot of module Eigengenes in the grey 60 module (A), turquoise module (B), green module (C), and red module (D). (E–H) Biological process GO terms for genes in the grey 60 module (E), turquoise module (F), green module (G), and red module (H). P-value cutoff = 0.01; q-value cutoff = 0.05. GO, Gene Ontology.

Identification of differentially expressed genes (DEGs)

We selected two available adipogenic differentiation datasets from the Gene Expression Omnibus (GEO (GEO100748 and GSE80614) to identify DEGs [29, 30]. The data from GEO100748 showed 1905 DEGs (Supplementary Table 5), and the data from GSE80614 showed 901 DEGs (Supplementary Table 6) (Adjusted p-value ≤ 0.05, |logFC| ≥ 1). We kept retained 459 genes from GEO100748 and 225 genes from GSE80614 whose absolute logFC values were greater than 2 (Figure 4A–4D). The overlapping genes among the key modules analyzed by WGCNA, during adipogenic differentiation from undifferentiating stage to day 4 or day 21, which identified 12 shared DEGs (Figure 4E). KEGG analysis of 459 DEGs from GEO100748 or GSE80614 indicated that the DEGs are highly related to cell cycle or PPAR signaling (Figure 4F, 4G). CPA4 was removed because of its controversial expression in the two different datasets. The expression of the other 11 genes was demonstrated during the different stages of adipogenic differentiation. The results identified that five shared downregulated genes during adipogenesis including TK1, SCARA3, PRC1, CENPF and CDC20, as well as six shared upregulated genes SAA1, WFDC1, GPX3, GHR, DDIT4L and COMP (Figure 5A, 5B).

Identification of shared genes during adipogenesis in two datasets. (A, B) Heatmap (A) and volcano plot (B) of the DEGs from GSE100748 (Adjusted p-value ≤ 0.05 and |logFC| ≥ 2). (C, D) Heatmap (C) and volcano plot (D) of the DEGs from GSE80614 (Adjust p-value ≤ 0.05 and |logFC| ≥ 2). (E) Venn diagram for shared genes among DEGs and key modules of GSE100748 and DEGs of GSE80614. (F) KEGG analysis of DEGs from GSE100748. (G) KEGG analysis of DEGs from GSE80614.

Figure 4. Identification of shared genes during adipogenesis in two datasets. (A, B) Heatmap (A) and volcano plot (B) of the DEGs from GSE100748 (Adjusted p-value ≤ 0.05 and |logFC| ≥ 2). (C, D) Heatmap (C) and volcano plot (D) of the DEGs from GSE80614 (Adjust p-value ≤ 0.05 and |logFC| ≥ 2). (E) Venn diagram for shared genes among DEGs and key modules of GSE100748 and DEGs of GSE80614. (F) KEGG analysis of DEGs from GSE100748. (G) KEGG analysis of DEGs from GSE80614.

Expression of 11 shared genes during adipogenesis. (A) Heatmap of the expression of 11 shared genes at different adipogenic stages including the undifferentiated stage, 1 hour, 2 hours, 3 hours, 5 hours, 6 hours, 12 hours, day 1, day 2, day 3, and day 4 (GSE80614). (B) Heatmap of the expression of 11 shared genes at different adipogenic stages, including the undifferentiated stage, day 7, and day 21 (GSE100748).

Figure 5. Expression of 11 shared genes during adipogenesis. (A) Heatmap of the expression of 11 shared genes at different adipogenic stages including the undifferentiated stage, 1 hour, 2 hours, 3 hours, 5 hours, 6 hours, 12 hours, day 1, day 2, day 3, and day 4 (GSE80614). (B) Heatmap of the expression of 11 shared genes at different adipogenic stages, including the undifferentiated stage, day 7, and day 21 (GSE100748).

The selection of shared key genes during adipogenesis

Based on the fact that PPARG, FABP4, LPL, and PPARA are crucial markers during adipogenesis [31, 32], we observed the correlation between the retained eleven genes and these four adipogenic markers. We calculated their co-expression based on the adipose tissue data from GTEX and found that only five genes (GPX3, GHR, SAA1, SCARA3, and WFDC1) correlated significantly with these markers. Moreover, GPX3, GHR and SAA1 correlated positively with adipogenic differentiation, while SCARA3 and WFDC1 correlated negatively with adipogenic differentiation (Figure 6A–6T). GPX3, GHR and SAA1 were removed from further study because they have been reported previously to be associated with obesity [33–35]. WFDC1 was filtered out because of the controversial result between its gradually increased expression during adipogenesis and its reverse correlation with adipogenic markers. The association of Scara3 with other four genes was showed in Supplementary Figure 1. During adipogenesis, although the deletion of Scara3 gene in AD-MSCs did not elevate the expression of Fabp4 compared to that in control group, it enhanced other adipogenic markers, such as Pparα, Pparγ, and Lpl (Figure 7A–7E). Together, these results verified the prediction that Scara3 inhibits adipogenesis. Thus, we focused on SCARA3 (included in turquoise model) for further analysis.

Co-expression of four shared genes with adipogenic markers in adipose tissue. (A–D) Correlation of GHR with PPARA (A), PPARG (B), FABP4 (C), and LPL (D) expression in adipose tissue, based on data from the Genotype Tissue Expression (GTEx) databases, respectively. (E–H) Correlation of GPX3 with PPARA (E), PPARG (F), FABP4 (G), and LPL (H) expression in adipose tissue, based on data from the GTEx databases, respectively. (I–L) Correlation of SAA1 with PPARA (I), PPARG (J), FABP4 (K), and LPL (L) expression in adipose tissue, based on data from the GTEx databases, respectively. (M–P) Correlation of SCARA3 with PPARA (M), PPARG (N), FABP4 (O), and LPL (P) expression in adipose tissue, based on data from GTEx databases, respectively. (Q–T) Correlation of WFDC1 with PPARA (Q), PPARG (R), FABP4 (S), and LPL (T) expression in adipose tissue, based on data from the GTEx databases, respectively.

Figure 6. Co-expression of four shared genes with adipogenic markers in adipose tissue. (A–D) Correlation of GHR with PPARA (A), PPARG (B), FABP4 (C), and LPL (D) expression in adipose tissue, based on data from the Genotype Tissue Expression (GTEx) databases, respectively. (E–H) Correlation of GPX3 with PPARA (E), PPARG (F), FABP4 (G), and LPL (H) expression in adipose tissue, based on data from the GTEx databases, respectively. (I–L) Correlation of SAA1 with PPARA (I), PPARG (J), FABP4 (K), and LPL (L) expression in adipose tissue, based on data from the GTEx databases, respectively. (M–P) Correlation of SCARA3 with PPARA (M), PPARG (N), FABP4 (O), and LPL (P) expression in adipose tissue, based on data from GTEx databases, respectively. (Q–T) Correlation of WFDC1 with PPARA (Q), PPARG (R), FABP4 (S), and LPL (T) expression in adipose tissue, based on data from the GTEx databases, respectively.

The deficiency of Scara3 inhibits adipogenic differentiation in vitro. The Ad-MSCs were transfected with NC siRNA and Scara3 SiRNA followed by adipogenic differentiation. (A) qRT-PCR analysis of depletion of Scara3. (B–E) qRT-PCR analysis of the relative levels of Pparα (B), PPARγ (C), Fabp4 (D), and Lpl (E). Error bars show standard deviation. *P

Figure 7. The deficiency of Scara3 inhibits adipogenic differentiation in vitro. The Ad-MSCs were transfected with NC siRNA and Scara3 SiRNA followed by adipogenic differentiation. (A) qRT-PCR analysis of depletion of Scara3. (B–E) qRT-PCR analysis of the relative levels of Pparα (B), PPARγ (C), Fabp4 (D), and Lpl (E). Error bars show standard deviation. *P < 0.05, **P < 0.01.

The prediction of transcription factors of SCARA3

To investigate the regulation of SCARA3, we constructed a Protein–protein interaction (PPI) network of adipogenic differentiation-associated DEGs (GSE100748) in the STRING database, and the top 10 hub genes were calculated in Cytoscape (Figure 8A). The result showed that the hub genes are mostly related to the AP-1 complex and JUN was ranked top. The PROMO database identified 62 transcription factors for SCARA3 (Maximum matrix dissimilarity rate = 10%) (Supplementary Figure 2). The top five transcription factors for SCARA3 by QIAGEN company were AP-1, ATF-2, GATA-1, PPAR-α and PPAR-γ. AP-1 is a heterodimer composed of proteins belonging to the c-FOS, c-Jun, ATF, and JDP families [36]. Thus, there were three shared transcription factors in PROMO database, including c-Jun, GATA-1, and PPAR-α (Figure 8B). Moreover, Gene Set Enrichment Analysis (GSEA) suggested that the signaling, “TNFA signaling via NFKB”, was inhibited during adipogenesis, and the key proteins identified in this signaling pathway also include various components of AP-1 (Figure 8C). We also found that SCARA3 correlated significantly and positively with JUN (Pearson r = 0.48, p value = 0.0) (Figure 8D). Negative correlations of JUN with PPARA further supported its potential role in adipogenesis (Figure 8E). To further confirm the correlation between SCARA3 with other genes associated with the signaling regulation, we interfered the Scara3 gene in AD-MSCs (Figure 8F). Quantitative real-time reverse transcription PCR (qRT-PCR) results were in line with the predicted correlation. With the deficiency of Scara3, the expression of Jun and Ppara were both decreased (Figure 8G, 8H).

Selection of the transcription factors of SCARA3. (A) PPI network of the DEGs from GSE100748. (B) Venn diagram for shared genes through the PROMO (a) and GENECARD (b) databases. (C) GSEA analysis of the DEGs from GSE100748. (D) Correlation of SCARA3 with JUN in expression in adipose tissue, based on data from GTEx databases. (E) Correlation of JUN with PPARA expression in adipose tissue, based on data from GTEx databases. (F) qRT-PCR analysis of depletion of Scara3. (G–H) qRT-PCR analysis of the relative levels of Pparα (G) and Jun (H). Error bars show standard deviation. **P

Figure 8. Selection of the transcription factors of SCARA3. (A) PPI network of the DEGs from GSE100748. (B) Venn diagram for shared genes through the PROMO (a) and GENECARD (b) databases. (C) GSEA analysis of the DEGs from GSE100748. (D) Correlation of SCARA3 with JUN in expression in adipose tissue, based on data from GTEx databases. (E) Correlation of JUN with PPARA expression in adipose tissue, based on data from GTEx databases. (F) qRT-PCR analysis of depletion of Scara3. (G–H) qRT-PCR analysis of the relative levels of Pparα (G) and Jun (H). Error bars show standard deviation. **P < 0.01, ***P < 0.001.

The indications of SCARA3 for phenotypic traits

It is well established that db/db mice and mice fed by a high-fat diet (HFD) presents obvious obesity [21, 31, 37–39]. To determine the role of SCARA3in vivo, we collected the inguinal white adipose tissue (iWAT) from db/db mice and mice fed a HFD or chow control. qRT-PCR and Western-Blot indicated that the expression of Scara3 decreased considerably in the iWAT in db/db mice compared with that in the controls (Figure 9A–9C). The similar decrease was observed in mice fed by HFD (Figure 9D–9F). To investigate the role of adipose tissue-derived mesenchymal stem cells (Ad-MSCs) in obesity, we analyzed the expression of Scara3 in Ad-MSCs from GSE115068. Compared with the expression of Scara3 in young Ad-MSCs, Scara3 expression was lower in old Ad-MSCs (Figure 9G). HFD stimulus contributed to the downregulation of Scara3 in young Ad-MSCs (Figure 9H). A previous study reported increased methylation of SCARA3 in patients with type 2 diabetes mellitus [40]. Here, we also found that SCARA3 methylation had some correlations with metabolic diseases, such as atherosclerotic lesions and type 2 diabetes via a search on the DiseaseMeth version 2.0 database (Figure 9I, 9J). Thus, we hypothesized that the decreased expression of SCARA3 during adipogenesis might be caused by its methylation. The Comparative Toxicogenomics Database (CTD; http://ctdbase.org/) provides information about interactions between environmental chemicals and gene products and their relationships to diseases (Figure 9K). We searched the top fifty relationships between SCARA3 with diseases. Intriguingly, the results implied that SCARA3 is associated with several metabolic disorders, such as weight loss, weight gain, glucose intolerance and insulin resistance (Figure 9K).

Phenotypical traits of SCARA3. (A) mRNA expression of Scara3 in iWAT of db/db mice. 8 mice were included in each group. (B, C) Protein expression of Scara3 in iWAT of db/db mice. 8 mice were included in in each group. (D) mRNA expression of Scara3 in iWAT of mice fed by HFD. 9 mice were included in each group. (E, F) Protein expression of Scara3 in iWAT of mice fed by HFD. 8 mice were included in each group. (G) Expression of Scara3 in white adipose tissue-derived mesenchymal stem cells (Ad-MSCs) of young and old mice. Each group has 2 mice, respectively. Data were obtained from GEO database (GSE115068). (H) Expression of Scara3 in Ad-MSCs of young mice fed a HFD. Each group has 2 mice, respectively. Data were obtained from GEO database (GSE115068). (I, J) Methylation of SCARA3 in metabolic disorders, including atherosclerotic lesions (I) and type 2 diabetes (J). (K) Relationship of weight-associated diseases with SCARA3 via a curated chemical interaction, based on the CTD database.

Figure 9. Phenotypical traits of SCARA3. (A) mRNA expression of Scara3 in iWAT of db/db mice. 8 mice were included in each group. (B, C) Protein expression of Scara3 in iWAT of db/db mice. 8 mice were included in in each group. (D) mRNA expression of Scara3 in iWAT of mice fed by HFD. 9 mice were included in each group. (E, F) Protein expression of Scara3 in iWAT of mice fed by HFD. 8 mice were included in each group. (G) Expression of Scara3 in white adipose tissue-derived mesenchymal stem cells (Ad-MSCs) of young and old mice. Each group has 2 mice, respectively. Data were obtained from GEO database (GSE115068). (H) Expression of Scara3 in Ad-MSCs of young mice fed a HFD. Each group has 2 mice, respectively. Data were obtained from GEO database (GSE115068). (I, J) Methylation of SCARA3 in metabolic disorders, including atherosclerotic lesions (I) and type 2 diabetes (J). (K) Relationship of weight-associated diseases with SCARA3 via a curated chemical interaction, based on the CTD database.

Discussion

Obesity generally results in predisposition towards other metabolic diseases, which can all result in a decline of life quality and life expectancy [5, 41, 42]. Thus, researchers need to find some therapeutic targets to prevent the development of obesity. Based on the important effect of adipogenesis on obesity, we aimed to find the causes of obesity from the perspective of adipogenesis. Firstly, we identified the overlapping genes in two adipogenic differentiation-associated datasets by performing WGCNA and DEGs selection. We identified five DEGs: GHR, GPX3, SAA1, SCARA3, and WFDC1. Except for SCARA3, other the four genes have already been implicated in the regulation of the pathogenesis of obesity [33–35]. Through further functional analysis, we revealed more indications that SCARA3 is an important gene during adipogenesis. Combining data form several databases, we further confirmed the possible mechanism mediated by SCARA3 and potential therapeutic effect of SCARA3 for obese individuals.

SCARA3 encodes a macrophage scavenger receptor-like protein. SCARA3 is ubiquitously expressed in human tissues, although it shows a relatively low expression in the liver and peripheral blood leukocytes [43]. Its ubiquitous expression was also revealed by the data from the GTEX database, which is a comprehensive public repository used to study tissue-specific gene expression and regulation. Cellular stress stimuli at low doses, such as UV radiation and hydrogen peroxide, both remarkably induced the expression of SCARA3 in human fibroblasts. However, excessive oxidative stress contributed to lower SCARA3 expression [43]. Thus, it has been reported to be a cellular stress response (CSR) gene with the role of protecting cells from oxidative stress by scavenging oxidative molecules or harmful products of oxidation [44, 45]. In addition, Yu et al. observed that SCARA3 might be a tumor suppressor-related gene, because down-regulation of SCARA3 was found in prostate cancer tissues and it is involved in cancers metastases and progression [46]. Obesity and its related disorders, such as diabetes and atherosclerosis, are thought to be related to impaired oxidative defense [10, 47, 48]. A previous study indicated that elevated methylation of SCARA3 probably disturbed the oxidative stress protection in type 2 diabetes mellitus (T2DM) [40]. Consistent with these observations, our analysis showed that hypermethylation of SCARA3 was associated with metabolic disorders, including T2DM and atherosclerosis lesions. Considering obesity can result from the environmental stimulus except for inheritable factors, epigenetic alterations, such as methylation, are a considerable cause of obesity [49]. Furthermore, previous studies revealed that promoter methylation contributed to the down-regulation of SCARA3 in prostate cancer [46, 50]. Based on this evidence, we assumed that the downregulation of SCARA3 during adipogenesis probably resulted from its methylation. The deletion of Scara3 in Ad-MSCs caused the downregulated expression of adipogenic markers. Notably, we verified the lower expression of Scara3 in db/db mice and in HFD-fed mice compared with that in the control groups. Consistently, downregulated expression of Scara3 was tested in Ad-MSCs isolated from young mice fed a HFD. Lower Scara3 expression was observed in young mice than in old mice. In addition, the data from the CTD database showed that SCARA3 is associated with several metabolic disorders, such as weight loss, weight gain, glucose intolerance and insulin resistance. Taken together, these results revealed the potentially important roles of SCARA3 in metabolic disorders.

Previous studies showed that scavenger receptor A gene regulatory elements targeted gene expression to macrophages and foam cells of atherosclerotic lesions [51, 52]. Furthermore, Horvai et al. demonstrated that scavenger receptor A gene regulatory elements contain binding sites for PU.1 and AP-1 in bone marrow progenitor cells [52]. In the present study, we predicted the transcription factors for SCARA3 using different databases and found the activator protein 1 (AP-1) was shared in the different databases. AP-1 is a heterodimer composed of proteins belonging to the c-FOS, c-Jun, ATF, and JDP families, which is involved in various cellular processes, such as cell differentiation, proliferation and apoptosis [36, 53]. PPAR signaling and nuclear factor κB (NF-κB) pathway were both reported to be related to obesity [42, 54–56]. The filtered DEGs expressed during adipogenesis enriched in “TNFA signaling via NFKB” by GSEA, and several genes at the leading edge (FOS, JUN, JUNB and ATF3) mostly belonged to the AP-1 heterodimer. Moreover, the hub genes selected from the PPI network, such as JUN, FOS and ATF3, also suggested that the AP-1 heterodimer plays a major role in adipogenesis. C-Jun levels were observed to be decreased during adipogenic differentiation in MC3T3 cells and inhibited adipogenic differentiation [57, 58]. Based on this evidence, we supposed that SCARA3 could bind to c-Jun and AP-1, thereby controlling adipogenic differentiation.

Collectively, we hypothesized that SCARA3 contributes to obesity and obesity-related metabolic complications in an age-dependent manner. Specifically, SCARA3 downregulation was possibly a consequence of its methylation. Mechanistically, the dysfunction of adipogenic differentiation might be associated with impaired modulation between SCARA3 and JUN. Further studies will be worthy of being performed to validate the role of SCARA3, such as the phenotypical observation of Scara3 knock-out mice and the methylation of SCARA3 promoter in mice with obesity. Other function studies, such as Chromatin Immunoprecipitation (ChIP)-PCR and luciferase reporter assay, can further verify SCARA3 gene regulates transcriptional factor JUN and PPAR signaling. Together, the present study combining bioinformatic analysis with the several biochemical functional analyses in vitro and in vivo, provided strong supports for the hypothesis that SCARA3 is a potential target for the treatment of obesity and other metabolic disorders.

Materials and Methods

Identification of key genes using WGCNA method

The R package “WGCNA” was used to find the key modules and genes related to adipogenesis [59]. Firstly, we input the raw expression of matrix from GSE100748 and transformed it into the log2 format. Then, we tested the missing values and removed the data with low quality. After the genes with zero variance were filtered out, we generated the sampleTree to check the outliers. Secondly, we normalized the data through quantile normalization and generated a heatmap of sample correlations. Thirdly, we constructed the network and detected the modules using an automatic, one-step method. The adjacency matrix was transformed into topological overlap matrix (TOM). According to the TOM-based dissimilarity measure, we set the soft-thresholding power as 14 (scale free R2 = 0.686), the minimal module size as 30, and cut height as 0.25 to identify key modules. Although the scale free R2 did not reach 0.8 when the power was 14, we considered this was because of the obvious biological difference between different stages of adipogenesis. Based on the number of samples (37 samples in total), we finally set the power as 14. Finally, we calculated the correlation between modules and traits during adipogenesis. Further data analysis based on the top four modules (grey 60, turquoise, green and red modules) during adipogenesis was performed.

Identification of DEGs

The series matrix files of GSE100748 and GSE80614 from GEO were downloaded. The data at different adipogenic stages from GSE100748 including the undifferentiated stage, day 7 and day 21 were retained. The data at different adipogenic stages from GSE80614 including the undifferentiated stage, and hour 1-6, hour 12 and day 1-4 were retained. To filter the DEGs, we only analyzed the expression at the undifferentiated stage and day 4. Firstly, we input the selected raw expression data of the matrix and transform it into the log2 format if needed. Secondly, we normalized data through quantile normalization. The R package “limma” [60] was applied to normalize the data and identify DEGs. The DEGs was retained if they had an adjust p-value ≤ 0.05 and |logFC| ≥ 1. Further data analysis we used the DEGs with an adjust p-value ≤ 0.05 and |logFC| ≥ 2 except for the GSEA analysis. Then, we exported the heatmap and volcano plot of DEGs with adjusted p-value ≤ 0.05 and |logFC| ≥ 2.

Functional enrichment analyses

Gene Ontology (GO) analyses and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were carried out using the R package “Clusterprofiler” [61]. We used DEGs (adjusted p-value <= 0.05 and |logFC| ≥ 2) from GSE80614, and transformed them into the format of ENTREZ ID. Then we use the R package “Clusterprofiler” to evaluate the enriched signaling (P value Cutoff = 0.05, Show Category = 5). R/Bioconductor package “Pi” was used to perform xPierGSEA analysis to find the pathways most associated with the hub genes. Here, the DEGs with |logFC| ≤ 1 were selected as our targeted gene set. The range size was set from 20 to 5000.

Construction of protein–protein interaction (PPI) networks of DEGs

STRING (http://string-db.org/) is a database of known and predicted protein-protein interactions. It was applied to predict the interacting genes among the DGEs using the default settings. The associated file of the network was imported into the Cytoscape software to visualize the network. The hub genes were identified using the plug-in cytohubba by setting the top node as 10 ranked by Degree.

Co-expression calculation

The raw expression data of genes and clinical phenotypes were downloaded from the GTEx program (https://www.gtexportal.org/). Then, we removed the data without associated clinical information. Further analysis then focused on correlations in adipose tissue. The Pearson correlation (r) and p-value were calculated.

Methylation and gene expression analyses

The human disease methylation database, DiseaseMeth version 2.0, is a web-based resource focusing on the aberrant methylomes of human diseases [62] (http://bioinfo.hrbmu.edu.cn/diseasemeth/). We utilized DiseaseMeth 2.0 to search the methylation levels of SCARA3 related to metabolic disorders.

CTD

The Comparative Toxicogenomics Database (CTD) illuminates how environmental chemicals affect human health [63] (http://ctdbase.org/). This database was used to validate SCARA3 associations with weight-related diseases and metabolic disorders. Here, we selected our diseases of interested only among the top 50 relationships, such as weight loss, weight gain, glucose intolerance, and insulin resistance.

Mice

Male C57BL/6J mice at 3 months of age were fed a HFD or chow control for four months before the start of the experiment to induce obesity. Db/db mice were used in our experiments at 12 weeks of age. All mice were maintained in a standard, specific pathogen-free facility of the Laboratory Animal Research Center of Central South University. All procedures involving mice were approved by the Animal Ethics Committee of the Central South University.

qRT-PCR analysis

Total RNA from the iWAT of mice or cells was extracted using Trizol reagent (Takara) as described previously [6, 64]. RNA (1000 ng) was reverse-transcribed into first-strand cDNA using a Reverse Transcription Kit (Takara). qPCR was then performed using SYBR Green PCR Master Mix (Takara) and mRNA expression was normalized to the expression of the reference gene Gapdh.

Western blot

1 mg adipose tissue was lysed in the mixture of RIPA (100 ul) and protease inhibitor (1:100). Western Blot was performed according to the previous described method [26]. The primary antibodies, Monoclonal Anti-SCARA3 antibody produced in mouse (#WH0051435M1-100UG; Sigma-Aldrich, MO, USA), and GAPDH (TA802519; ORIGENE), were incubated overnight at 4° C, then incubated with appropriate secondary antibodies for 1 hour at room temperature. The blots were visualized using ECL detection reagents.

Cell isolation and culture

Ad-MSCs were isolated from iWAT of 1-month male C57BL/6J mice. Firstly, we isolated the iWAT from mice. Then, the tissue was cut and digested using 0.1% I type Collagenase (Gibco,17108-029), a mixture of the collagenases I in PBS. Then it was incubated in 37° C for 30 min and followed by centrifuge at 1000 rpm 5 min. The cells were then cultured in DMEM containing 10% fetal bovine serum, 100 U/mL penicillin, and 100 μg/ mL streptomycin. Culture medium was changed every other day. Cells were collected for cell transfection, adipogenic differentiation or RNA extraction.

Cell transfection

The Scara3 siRNA and the negative control (NC) were purchased from Ribibio (Guangzhou, China). The siRNAs were transfected at the concentration of 100 nM using lipofectamine 2000 (Invitrogen, USA) according to manufacturer's recommendations.

Adipogenic differentiation

Adipogenic differentiation was conducted as described previously [25]. MSCs were treated with DMEM containing 10% fetal bovine serum, 100 U/mL penicillin, 100 μg/ mL streptomycin, 0.5 mM 3-isobutyl-1-methylxanthine, 5 μg/ml insulin, and 1 μM dexamethasone for 4 days. The medium was changed every other day.

Statistical analysis

All data were analyzed using R and R studio software (Version 4.0.1) except for GSEA analysis (Version 3.6.3). Key modules were selected with the power as 14 (scale free R2 = 0.686), the minimal module size as 30, and the cut height as 0.25 using package “WGCNA”. DEGs were selected with adjusted P < 0.05 and |logFC| > 1 through the package “limma”. GO terms or KEGG pathways with adjusted P < 0.05 were considered statistically significant. GSEA was generated using the package “Pi”. Pearson correlation analysis was adopted to determine the linear relationship between the two groups. For the animal experiment, each group comprised ten mice. The animals were randomly divided into two groups including mice fed a HFD and those fed a control chow diet. qRT-PCR was performed independently three times and all results are expressed as means ± standard deviation (SD). Quantitative data from qRT-PCR, Western Blot or the expression of Scara3 in white adipose tissue derived mesenchymal stem cells (Ad-MSCs) from the GEO database were imported into a spreadsheet and scaled and normalized to their appropriate controls using R. Two-way, paired or unpaired t-tests were performed using the package “ggplot”.

Availability of data and materials

The raw data of the two adipogenesis-related microarray datasets are available in the GEO database (https://www.ncbi.nlm.nih.gov/geo/. Accession number: GSE100748 and GSE80614.). The data for Scara3 expression in Ad-MSCs are available in the GEO database (Accession number: GSE115068). The transcription factors were predicted in the PROMO database (http://alggen.lsi.upc.es/cgi-bin/promo_v3/promo/promoinit.cgi?dirDB=TF_8.3). Associations of SCARA3 with metabolic disorders are available in CTD (http://ctdbase.org/). The methylation of SCARA3 in patients with metabolic disorders is available in DiseaseMeth version 2.0 (http://bio-bigdata.hrbmu.edu.cn/diseasemeth/).

Supplementary Materials

Supplementary Figures

Supplementary Table 1

Supplementary Table 2

Supplementary Table 3

Supplementary Table 4

Supplementary Table 5

Supplementary Table 6

Abbreviations

GEO: Gene Expression Omnibus; DEGs: Differentially expressed genes; WGCNA: Weighted Gene Co-Expression Network Analysis; TOM: topological overlap matrix; PPI: Protein–protein interaction; GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; AP-1: activator protein 1; CTD: Comparative Toxicogenomics Database; GSEA: Gene Set Enrichment Analysis; HFD: High-fat diet; iWAT: inguinal white adipose tissue; Ad-MSCs: adipose tissue-derived mesenchymal stem cells; MSCs: Mesenchymal stem cells; GHR: growth hormone receptor; GPX3: glutathione peroxidase 3; SAA1: serum amyloid A1; SCARA3: scavenger receptor class A member 3; WFDC1: WAP four-disulfide core domain 1; qRT-PCR: Quantitative real-time reverse transcription PCR; NC: Negative control.

Author Contributions

Hui Peng designed the study and carried out bioinformatical data analysis; Hui Peng and Yan Huang performed most of biological experiments; Yan Huang, Qi Guo, Ye Xiao, Tian Su, Chang-Jun Li fed the animals and collected the samples; Hui Peng wrote the manuscript; Yan Huang and Xiang-Hang Luo supervised the data analysis and revised the manuscript. All authors reviewed and approved the final manuscript.

Acknowledgments

The authors thank ELIXIGEN Co. for English editing of the final manuscript. This work was supported by the grant of the Key Program of the National Natural Science Foundation of China (No. 81930022), Projects of International Cooperation and Exchanges the National Natural Science Foundation of China (No. 81520108008), Training Program of the Major Research Plan of the National Natural Science Foundation of China (No. 91749105), Young Scientists Fund of the National Natural Science Foundation of China (No. 81700785) and National Natural Science Foundation of China (No. 82000811).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Funding

This work was supported by the grant of the Key Program of the National Natural Science Foundation of China (No. 81930022), Projects of International Cooperation and Exchanges the National Natural Science Foundation of China (No. 81520108008), Training Program of the Major Research Plan of the National Natural Science Foundation of China (No. 91749105), Young Scientists Fund of the National Natural Science Foundation of China (No. 81700785) and National Natural Science Foundation of China (No. 82000811).

References

1. Ghaben AL, Scherer PE. Adipogenesis and metabolic health. Nat Rev Mol Cell Biol. 2019; 20:242–58. https://doi.org/10.1038/s41580-018-0093-z [PubMed]
2. Heindel JJ, Newbold R, Schug TT. Endocrine disruptors and obesity. Nat Rev Endocrinol. 2015; 11:653–61. https://doi.org/10.1038/nrendo.2015.163 [PubMed]
3. Su T, Xiao Y, Xiao Y, Guo Q, Li C, Huang Y, Deng Q, Wen J, Zhou F, Luo XH. Bone marrow mesenchymal stem cells-derived exosomal MiR-29b-3p regulates aging-associated insulin resistance. ACS Nano. 2019; 13:2450–62. https://doi.org/10.1021/acsnano.8b09375 [PubMed]
4. Al Batran R, Gopal K, Capozzi ME, Chahade JJ, Saleme B, Tabatabaei-Dakhili SA, Greenwell AA, Niu J, Almutairi M, Byrne NJ, Masson G, Kim R, Eaton F, et al. Pimozide alleviates hyperglycemia in diet-induced obesity by inhibiting skeletal muscle ketone oxidation. Cell Metab. 2020; 31:909–19.e8. https://doi.org/10.1016/j.cmet.2020.03.017 [PubMed]
5. Blüher M. Obesity: global epidemiology and pathogenesis. Nat Rev Endocrinol. 2019; 15:288–98. https://doi.org/10.1038/s41574-019-0176-8 [PubMed]
6. Huang Y, Xiao Y, Liu Y, Guo M, Guo Q, Zhou F, Liu T, Su T, Xiao Y, Luo XH. MicroRNA-188 regulates aging-associated metabolic phenotype. Aging Cell. 2020; 19:e13077. https://doi.org/10.1111/acel.13077 [PubMed]
7. Adlung L, Amit I, Elinav E. Embrace the fat when getting old. Aging (Albany NY). 2019; 11:8730–32. https://doi.org/10.18632/aging.102341 [PubMed]
8. Lettieri Barbato D, Aquilano K. Feast and famine: adipose tissue adaptations for healthy aging. Ageing Res Rev. 2016; 28:85–93. https://doi.org/10.1016/j.arr.2016.05.007 [PubMed]
9. Gao H, Arner P, Beauchef G, Guéré C, Vie K, Dahlman I, Mejhert N, Rydén M. Age-induced reduction in human lipolysis: a potential role for adipocyte noradrenaline degradation. Cell Metab. 2020; 32:1–3. https://doi.org/10.1016/j.cmet.2020.06.007 [PubMed]
10. Okuno Y, Fukuhara A, Hashimoto E, Kobayashi H, Kobayashi S, Otsuki M, Shimomura I. Oxidative stress inhibits healthy adipose expansion through suppression of SREBF1-mediated lipogenic pathway. Diabetes. 2018; 67:1113–27. https://doi.org/10.2337/db17-1032 [PubMed]
11. Redman LM, Smith SR, Burton JH, Martin CK, Il’yasova D, Ravussin E. Metabolic slowing and reduced oxidative damage with sustained caloric restriction support the rate of living and oxidative damage theories of aging. Cell Metab. 2018; 27:805–15.e4. https://doi.org/10.1016/j.cmet.2018.02.019 [PubMed]
12. Rosen ED, Spiegelman BM. What we talk about when we talk about fat. Cell. 2014; 156:20–44. https://doi.org/10.1016/j.cell.2013.12.012 [PubMed]
13. Hammarstedt A, Gogg S, Hedjazifar S, Nerstedt A, Smith U. Impaired adipogenesis and dysfunctional adipose tissue in human hypertrophic obesity. Physiol Rev. 2018; 98:1911–41. https://doi.org/10.1152/physrev.00034.2017 [PubMed]
14. Jeffery E, Wing A, Holtrup B, Sebo Z, Kaplan JL, Saavedra-Peña R, Church CD, Colman L, Berry R, Rodeheffer MS. The Adipose Tissue Microenvironment Regulates Depot-Specific Adipogenesis in Obesity. Cell Metab. 2016; 24:142–50. https://doi.org/10.1016/j.cmet.2016.05.012 [PubMed]
15. Lettieri-Barbato D, Giovannetti E, Aquilano K. Effects of dietary restriction on adipose mass and biomarkers of healthy aging in human. Aging (Albany NY). 2016; 8:3341–55. https://doi.org/10.18632/aging.101122 [PubMed]
16. Jeffery E, Church CD, Holtrup B, Colman L, Rodeheffer MS. Rapid depot-specific activation of adipocyte precursor cells at the onset of obesity. Nat Cell Biol. 2015; 17:376–85. https://doi.org/10.1038/ncb3122 [PubMed]
17. Arner P, Andersson DP, Thörne A, Wirén M, Hoffstedt J, Näslund E, Thorell A, Rydén M. Variations in the size of the major omentum are primarily determined by fat cell number. J Clin Endocrinol Metab. 2013; 98:E897–901. https://doi.org/10.1210/jc.2012-4106 [PubMed]
18. Wang W, Ishibashi J, Trefely S, Shao M, Cowan AJ, Sakers A, Lim HW, O’Connor S, Doan MT, Cohen P, Baur JA, King MT, Veech RL, et al. A PRDM16-driven metabolic signal from adipocytes regulates precursor cell fate. Cell Metab. 2019; 30:174–89.e5. https://doi.org/10.1016/j.cmet.2019.05.005 [PubMed]
19. Wang QA, Tao C, Gupta RK, Scherer PE. Tracking adipogenesis during white adipose tissue development, expansion and regeneration. Nat Med. 2013; 19:1338–44. https://doi.org/10.1038/nm.3324 [PubMed]
20. Gross B, Pawlak M, Lefebvre P, Staels B. PPARs in obesity-induced T2DM, dyslipidaemia and NAFLD. Nat Rev Endocrinol. 2017; 13:36–49. https://doi.org/10.1038/nrendo.2016.135 [PubMed]
21. Liu L, Li T, Liao Y, Wang Y, Gao Y, Hu H, Huang H, Wu F, Chen YG, Xu S, Fu S. Triose kinase controls the lipogenic potential of fructose and dietary tolerance. Cell Metab. 2020; 32:605–18.e7. https://doi.org/10.1016/j.cmet.2020.07.018 [PubMed]
22. Matsushita K, Dzau VJ. Mesenchymal stem cells in obesity: insights for translational applications. Lab Invest. 2017; 97:1158–1166. https://doi.org/10.1038/labinvest.2017.42 [PubMed]
23. Chen Q, Shou P, Zheng C, Jiang M, Cao G, Yang Q, Cao J, Xie N, Velletri T, Zhang X, Xu C, Zhang L, Yang H, et al. Fate decision of mesenchymal stem cells: adipocytes or osteoblasts? Cell Death Differ. 2016; 23:1128–39. https://doi.org/10.1038/cdd.2015.168 [PubMed]
24. Li CJ, Cheng P, Liang MK, Chen YS, Lu Q, Wang JY, Xia ZY, Zhou HD, Cao X, Xie H, Liao EY, Luo XH. MicroRNA-188 regulates age-related switch between osteoblast and adipocyte differentiation. J Clin Invest. 2015; 125:1509–22. https://doi.org/10.1172/JCI77716 [PubMed]
25. Li CJ, Xiao Y, Yang M, Su T, Sun X, Guo Q, Huang Y, Luo XH. Long noncoding RNA Bmncr regulates mesenchymal stem cell fate during skeletal aging. J Clin Invest. 2018; 128:5251–66. https://doi.org/10.1172/JCI99044 [PubMed]
26. Xiao YZ, Yang M, Xiao Y, Guo Q, Huang Y, Li CJ, Cai D, Luo XH. Reducing hypothalamic stem cell senescence protects against aging-associated physiological decline. Cell Metab. 2020; 31:534–48.e5. https://doi.org/10.1016/j.cmet.2020.01.002 [PubMed]
27. Wu D, Wong CK, Han JM, Orban PC, Huang Q, Gillies J, Mojibian M, Gibson WT, Levings MK. T reg-specific insulin receptor deletion prevents diet-induced and age-associated metabolic syndrome. J Exp Med. 2020; 217:e20191542. https://doi.org/10.1084/jem.20191542 [PubMed]
28. Gnad T, Navarro G, Lahesmaa M, Reverte-Salisa L, Copperi F, Cordomi A, Naumann J, Hochhäuser A, Haufs-Brusberg S, Wenzel D, Suhr F, Jespersen NZ, Scheele C, et al. adenosine/A2B receptor signaling ameliorates the effects of aging and counteracts obesity. Cell Metab. 2020; 32:56–70.e7. https://doi.org/10.1016/j.cmet.2020.06.006 [PubMed]
29. Vaz C, Bee Tee BT, Yong D, Yi Lee Q, Tanavde V. Mesenchymal Stromal Cells Derived from Human Embryonic Stem Cells, Fetal Limb and Bone Marrow Share a Common Phenotype but Are Transcriptionally and Biologically Different. Stem Cell Discov. 2017; 7:1–26. https://doi.org/10.4236/scd.2017.71001
30. van de Peppel J, Strini T, Tilburg J, Westerhoff H, van Wijnen AJ, van Leeuwen JP. Identification of three early phases of cell-fate determination during osteogenic and adipogenic differentiation by transcription factor dynamics. Stem Cell Reports. 2017; 8:947–60. https://doi.org/10.1016/j.stemcr.2017.02.018 [PubMed]
31. Moseti D, Regassa A, Kim WK. Molecular regulation of adipogenesis and potential anti-adipogenic bioactive molecules. Int J Mol Sci. 2016; 17:124. https://doi.org/10.3390/ijms17010124 [PubMed]
32. Garin-Shkolnik T, Rudich A, Hotamisligil GS, Rubinstein M. FABP4 attenuates PPARγ and adipogenesis and is inversely correlated with PPARγ in adipose tissues. Diabetes. 2014; 63:900–11. https://doi.org/10.2337/db13-0436 [PubMed]
33. Chiba T, Han CY, Vaisar T, Shimokado K, Kargi A, Chen MH, Wang S, McDonald TO, O’Brien KD, Heinecke JW, Chait A. Serum amyloid A3 does not contribute to circulating SAA levels. J Lipid Res. 2009; 50:1353–62. https://doi.org/10.1194/jlr.M900089-JLR200 [PubMed]
34. Yang X, Deignan JL, Qi H, Zhu J, Qian S, Zhong J, Torosyan G, Majid S, Falkard B, Kleinhanz RR, Karlsson J, Castellani LW, Mumick S, et al. Validation of candidate causal genes for obesity that affect shared metabolic pathways and networks. Nat Genet. 2009; 41:415–23. https://doi.org/10.1038/ng.325 [PubMed]
35. List EO, Sackmann-Sala L, Berryman DE, Funk K, Kelder B, Gosney ES, Okada S, Ding J, Cruz-Topete D, Kopchick JJ. Endocrine parameters and phenotypes of the growth hormone receptor gene disrupted (GHR-/-) mouse. Endocr Rev. 2011; 32:356–86. https://doi.org/10.1210/er.2010-0009 [PubMed]
36. Curran T, Franza BR Jr. Fos and Jun: the AP-1 connection. Cell. 1988; 55:395–97. https://doi.org/10.1016/0092-8674(88)90024-4 [PubMed]
37. Xu J, Bartolome CL, Low CS, Yi X, Chien CH, Wang P, Kong D. Genetic identification of leptin neural circuits in energy and glucose homeostases. Nature. 2018; 556:505–09. https://doi.org/10.1038/s41586-018-0049-7 [PubMed]
38. Martinez-Botas J, Anderson JB, Tessier D, Lapillonne A, Chang BH, Quast MJ, Gorenstein D, Chen KH, Chan L. Absence of perilipin results in leanness and reverses obesity in Lepr(db/db) mice. Nat Genet. 2000; 26:474–79. https://doi.org/10.1038/82630 [PubMed]
39. Kim JD, Yoon NA, Jin S, Diano S. Microglial UCP2 mediates inflammation and obesity induced by high-fat feeding. Cell Metab. 2019; 30:952–62.e5. https://doi.org/10.1016/j.cmet.2019.08.010 [PubMed]
40. Karachanak-Yankova S, Dimova R, Nikolova D, Nesheva D, Koprinarova M, Maslyankov S, Tafradjiska R, Gateva P, Velizarova M, Hammoudeh Z, Stoynev N, Toncheva D, Tankova T, Dimova I. Epigenetic alterations in patients with type 2 diabetes mellitus. Balkan J Med Genet. 2016; 18:15–24. https://doi.org/10.1515/bjmg-2015-0081 [PubMed]
41. Neeland IJ, Poirier P, Després JP. Cardiovascular and Metabolic Heterogeneity of Obesity: Clinical Challenges and Implications for Management. Circulation. 2018; 137:1391–406. https://doi.org/10.1161/CIRCULATIONAHA.117.029617 [PubMed]
42. El Ouarrat D, Isaac R, Lee YS, Oh DY, Wollam J, Lackey D, Riopel M, Bandyopadhyay G, Seo JB, Sampath-Kumar R, Olefsky JM. TAZ is a negative regulator of PPARγ activity in adipocytes and TAZ deletion improves insulin sensitivity and glucose tolerance. Cell Metab. 2020; 31:162–73.e5. https://doi.org/10.1016/j.cmet.2019.10.003 [PubMed]
43. Han HJ, Tokino T, Nakamura Y. CSR, a scavenger receptor-like protein with a protective role against cellular damage causedby UV irradiation and oxidative stress. Hum Mol Genet. 1998; 7:1039–46. https://doi.org/10.1093/hmg/7.6.1039 [PubMed]
44. Guadagno NA, Moriconi C, Licursi V, D’Acunto E, Nisi PS, Carucci N, De Jaco A, Cacci E, Negri R, Lupo G, Miranda E. Neuroserpin polymers cause oxidative stress in a neuronal model of the dementia FENIB. Neurobiol Dis. 2017; 103:32–44. https://doi.org/10.1016/j.nbd.2017.03.010 [PubMed]
45. Cao M, Cai J, Yuan Y, Shi Y, Wu H, Liu Q, Yao Y, Chen L, Dang W, Zhang X, Xiao J, Yang K, He Z, et al. A four-gene signature-derived risk score for glioblastoma: prospects for prognostic and response predictive analyses. Cancer Biol Med. 2019; 16:595–605. https://doi.org/10.20892/j.issn.2095-3941.2018.0277 [PubMed]
46. Yu G, Tseng GC, Yu YP, Gavel T, Nelson J, Wells A, Michalopoulos G, Kokkinakis D, Luo JH. CSR1 suppresses tumor growth and metastasis of prostate cancer. Am J Pathol. 2006; 168:597–607. https://doi.org/10.2353/ajpath.2006.050620 [PubMed]
47. Lee JM, Govindarajah V, Goddard B, Hinge A, Muench DE, Filippi MD, Aronow B, Cancelas JA, Salomonis N, Grimes HL, Reynaud D. Obesity alters the long-term fitness of the hematopoietic stem cell compartment through modulation of Gfi1 expression. J Exp Med. 2018; 215:627–44. https://doi.org/10.1084/jem.20170690 [PubMed]
48. Cienfuegos S, Gabel K, Kalam F, Ezpeleta M, Wiseman E, Pavlou V, Lin S, Oliveira ML, Varady KA. Effects of 4- and 6-h time-restricted feeding on weight and cardiometabolic health: a randomized controlled trial in adults with obesity. Cell Metab. 2020; 32:366–78.e3. https://doi.org/10.1016/j.cmet.2020.06.018 [PubMed]
49. Ling C, Rönn T. Epigenetics in human obesity and type 2 diabetes. Cell Metab. 2019; 29:1028–44. https://doi.org/10.1016/j.cmet.2019.03.009 [PubMed]
50. Bock AJ, Nymoen DA, Brenne K, Kærn J, Davidson B. SCARA3 mRNA is overexpressed in ovarian carcinoma compared with breast carcinoma effusions. Hum Pathol. 2012; 43:669–74. https://doi.org/10.1016/j.humpath.2011.06.003 [PubMed]
51. Cozen AE, Moriwaki H, Kremen M, DeYoung MB, Dichek HL, Slezicki KI, Young SG, Véniant M, Dichek DA. Macrophage-targeted overexpression of urokinase causes accelerated atherosclerosis, coronary artery occlusions, and premature death. Circulation. 2004; 109:2129–35. https://doi.org/10.1161/01.CIR.0000127369.24127.03 [PubMed]
52. Horvai A, Palinski W, Wu H, Moulton KS, Kalla K, Glass CK. Scavenger receptor a gene regulatory elements target gene expression to macrophages and to foam cells of atherosclerotic lesions. Proc Natl Acad Sci USA. 1995; 92:5391–95. https://doi.org/10.1073/pnas.92.12.5391 [PubMed]
53. Hess J, Angel P, Schorpp-Kistner M. AP-1 subunits: quarrel and harmony among siblings. J Cell Sci. 2004; 117:5965–73. https://doi.org/10.1242/jcs.01589 [PubMed]
54. Malle EK, Zammit NW, Walters SN, Koay YC, Wu J, Tan BM, Villanueva JE, Brink R, Loudovaris T, Cantley J, McAlpine SR, Hesselson D, Grey ST. Nuclear factor κB-inducing kinase activation as a mechanism of pancreatic β cell failure in obesity. J Exp Med. 2015; 212:1239–54. https://doi.org/10.1084/jem.20150218 [PubMed]
55. Sui Y, Park SH, Xu J, Monette S, Helsley RN, Han SS, Zhou C. IKKβ links vascular inflammation to obesity and atherosclerosis. J Exp Med. 2014; 211:869–86. https://doi.org/10.1084/jem.20131281 [PubMed]
56. Nakamura M, Liu T, Husain S, Zhai P, Warren JS, Hsu CP, Matsuda T, Phiel CJ, Cox JE, Tian B, Li H, Sadoshima J. Glycogen synthase kinase-3α promotes fatty acid uptake and lipotoxic cardiomyopathy. Cell Metab. 2019; 29:1119–34.e12. https://doi.org/10.1016/j.cmet.2019.01.005 [PubMed]
57. Lee DS, Choi H, Han BS, Kim WK, Lee SC, Oh KJ, Bae KH. C-jun regulates adipocyte differentiation via the KLF15-mediated mode. Biochem Biophys Res Commun. 2016; 469:552–58. https://doi.org/10.1016/j.bbrc.2015.12.035 [PubMed]
58. Wang H, Scott RE. Adipocyte differentiation selectively represses the serum inducibility of c-jun and junB by reversible transcription-dependent mechanisms. Proc Natl Acad Sci USA. 1994; 91:4649–53. https://doi.org/10.1073/pnas.91.11.4649 [PubMed]
59. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008; 9:559. https://doi.org/10.1186/1471-2105-9-559 [PubMed]
60. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015; 43:e47. https://doi.org/10.1093/nar/gkv007 [PubMed]
61. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012; 16:284–87. https://doi.org/10.1089/omi.2011.0118 [PubMed]
62. Xiong Y, Wei Y, Gu Y, Zhang S, Lyu J, Zhang B, Chen C, Zhu J, Wang Y, Liu H, Zhang Y. DiseaseMeth version 2.0: a major expansion and update of the human disease methylation database. Nucleic Acids Res. 2017; 45:D888–95. https://doi.org/10.1093/nar/gkw1123 [PubMed]
63. Davis AP, Grondin CJ, Johnson RJ, Sciaky D, McMorran R, Wiegers J, Wiegers TC, Mattingly CJ. The comparative toxicogenomics database: update 2019. Nucleic Acids Res. 2019; 47:D948–54. https://doi.org/10.1093/nar/gky868 [PubMed]
64. Yang M, Guo Q, Peng H, Xiao YZ, Xiao Y, Huang Y, Li CJ, Su T, Zhang YL, Lei MX, Chen HL, Jiang TJ, Luo XH. Krüppel-like factor 3 inhibition by mutated lncRNA Reg1cp results in human high bone mass syndrome. J Exp Med. 2019; 216:1944–64. https://doi.org/10.1084/jem.20181554 [PubMed]