Transcriptomic and Proteomic Analysis Reveals the Potential Role of RBMS1 in Adipogenesis and Adipocyte Metabolism

Adipocytes play a critical role in maintaining a healthy systemic metabolism by storing and releasing energy in the form of fat and helping to regulate glucose and lipid levels in the body. Adipogenesis is the process through which pre-adipocytes are differentiated into mature adipocytes. It is a complex process involving various transcription factors and signaling pathways. The dysregulation of adipogenesis has been implicated in the development of obesity and metabolic disorders. Therefore, understanding the mechanisms that regulate adipogenesis and the factors that contribute to its dysregulation may provide insights into the prevention and treatment of these conditions. RNA-binding motif single-stranded interacting protein 1 (RBMS1) is a protein that binds to RNA and plays a critical role in various cellular processes such as alternative splicing, mRNA stability, and translation. RBMS1 polymorphism has been shown to be associated with obesity and type 2 diabetes, but the role of RBMS1 in adipose metabolism and adipogenesis is not known. We show that RBMS1 is highly expressed during the early phase of the differentiation of the murine adipocyte cell line 3T3-L1 and is significantly upregulated in the adipose tissue depots and adipocytes of high-fat-fed mice, implying a possible role in adipogenesis and adipose metabolism. Knockdown of RBMS1 in pre-adipocytes impacted the differentiation process and reduced the expression of some of the key adipogenic markers. Transcriptomic and proteomic analysis indicated that RBMS1 depletion affected the expression of several genes involved in major metabolic processes, including carbohydrate and lipid metabolism. Our findings imply that RBMS1 plays an important role in adipocyte metabolism and may offer novel therapeutic opportunity for metabolic disorders such as obesity and type 2 diabetes.


Introduction
Obesity is a complex disease characterized by an excessive amount of body fat [1,2]. It is a major public health concern worldwide and is associated with numerous health problems such as heart disease, diabetes, and certain types of cancer [3][4][5][6][7][8][9][10]. Furthermore, obese people are more likely to develop severe disease, require hospitalization, and die as a result of viral infections, including the recent SARS-CoV-2 pandemic [11][12][13]. The key factor in the pathophysiology of obesity-related metabolic abnormalities is obesity-induced adipose tissue dysfunction [14,15]. Adipose tissue dysfunction leads to the release of pro-inflammatory cytokines and adipokines, which contribute to insulin resistance and chronic low-grade inflammation [13][14][15]. Adipose tissue is a kind of connective tissue that stores energy in the form of fat and plays a crucial role in regulating metabolism, hormone production, and immune function [16]. In obesity, adipose tissue becomes inflamed and dysfunctional, leading to insulin resistance and other metabolic disturbances that contribute to the development of chronic diseases [17]. Adipocytes, the primary component of adipose tissue, are derived from multipotent mesenchymal stem cells (MSCs) through adipogenesis, a complicated and highly regulated process. Understanding the mechanisms that regulate adipogenesis and the factors that contribute to its dysregulation may provide insights into the prevention and treatment of these conditions [18,19].
The transcriptional regulation of adipogenesis has been comprehensively understood, and several transcription factors have been identified as key regulators of this process, including peroxisome proliferator-activated receptor gamma (PPARγ), CCAAT/enhancerbinding protein alpha (C/EBPα), and sterol regulatory element-binding protein 1c (SREBP-1c) [19,20]. Recent studies have highlighted the post-transcriptional regulation of adipogenesis, including the role of microRNAs and RNA-binding proteins (RBPs) in modulating adipocyte differentiation and function [21][22][23][24]. RBPs such as HuR, QK, and YBX1 have been shown to regulate adipogenesis by controlling the stability and translation of key mRNAs involved in this process, suggesting a new avenue for therapeutic intervention [25][26][27].
RNA-binding motif single-stranded interacting protein 1 (RBMS1) is a member of the RNA-binding protein family and plays key roles in the post-transcriptional regulation of gene expression. RBMS1 has been shown to be involved in various cellular processes, including mRNA splicing and stability, as well as translation initiation. Multiple studies have shown that RBMS1 expression is altered in several types of cancer, and its dysregulation has been associated with tumor progression and metastasis [28][29][30][31]. RBMS1 may play a role in the development of coronary heart disease, possibly through its involvement in lipid metabolism and inflammation [32]. In addition, RBMS1 polymorphism has been found to be significantly associated with an increased risk of developing type 2 diabetes in several populations, suggesting a possible role in metabolic regulation as well [33][34][35]. No previous study has explored the role of RBMS1 in adipogenesis and adipose metabolism. This study was undertaken to investigate the role of RBMS1 in adipocyte differentiation and lipid metabolism. We show that RBMS1 may play an important role in the regulation of adipocyte differentiation and metabolism and may offer a novel therapeutic opportunity for metabolic disorders such as obesity and type 2 diabetes.

RBMS1 Expression Is Upregulated in Adipose Tissue Depots and Adipocytes of Mice Fed a High-Fat Diet
To understand the functional significance of RBMS1 in adipose tissue metabolism, we first evaluated the effect of high-fat diet (HFD) on the expression of Rbms1 in mouse adipose tissue depots. C57BL/6J male mice (four in each group) were fed either a regular chow diet (CD) or HFD for 18 weeks, as described in Section 4. After 18 weeks, HFD-fed mice gained more weight and fat mass compared to those fed a regular diet ( Figure 1A,B). Moreover, HFD mice developed hyperglycemia and hypertriglyceridemia, showing dietinduced obesity and insulin resistance ( Figure 1C,D). The expression of Rbms1 in adipose tissue depots in mice from the CD and HFD groups was subsequently assessed using quantitative real-time PCR using specific primers. Rbms1 expression increased significantly in the subcutaneous (scWAT) and visceral (vsWAT) white adipose tissue depots of HFD mice compared to CD animals ( Figure 1E,F). Since diet-induced obesity (DIO) is known to alter the cellular composition of adipose tissue, which may have a direct impact on the expression of genes, including Rbms1, we re-analyzed an RNA sequencing dataset (GSE #142187) to evaluate the expression of Rbms1 in the adipocyte fraction [36]. The dataset comprised RNA expression data of adipocytes extracted from white adipose tissue of mice fed a CD or HFD for 8 and 20 weeks. The mRNA expression of Rbms1 was significantly increased in the adipocytes of HFD mice as compared to the CD mice, which is consistent with our data obtained from whole adipose tissue ( Figure 1G,H).

High-Fat Diet
To understand the functional significance of RBMS1 in adipose tissue metabolism, we first evaluated the effect of high-fat diet (HFD) on the expression of Rbms1 in mouse adipose tissue depots. C57BL/6J male mice (four in each group) were fed either a regular chow diet (CD) or HFD for 18 weeks, as described in Section 4. After 18 weeks, HFD-fed mice gained more weight and fat mass compared to those fed a regular diet ( Figure 1A,B). Moreover, HFD mice developed hyperglycemia and hypertriglyceridemia, showing dietinduced obesity and insulin resistance ( Figure 1C,D). The expression of Rbms1 in adipose tissue depots in mice from the CD and HFD groups was subsequently assessed using quantitative real-time PCR using specific primers. Rbms1 expression increased significantly in the subcutaneous (scWAT) and visceral (vsWAT) white adipose tissue depots of HFD mice compared to CD animals ( Figure 1E,F). Since diet-induced obesity (DIO) is known to alter the cellular composition of adipose tissue, which may have a direct impact on the expression of genes, including Rbms1, we re-analyzed an RNA sequencing dataset (GSE #142187) to evaluate the expression of Rbms1 in the adipocyte fraction [36]. The dataset comprised RNA expression data of adipocytes extracted from white adipose tissue of mice fed a CD or HFD for 8 and 20 weeks. The mRNA expression of Rbms1 was significantly increased in the adipocytes of HFD mice as compared to the CD mice, which is consistent with our data obtained from whole adipose tissue ( Figure 1G,H). Rbms1 expression is significantly altered in adipose tissue depots and adipocytes of mice fed a high-fat diet (HFD). (A-D) Body weight, fat mass, plasma triglyceride levels, and plasma glucose levels in the two groups of mice. (E,F) Rbms1 expression in subcutaneous (E) and mesenteric adipose tissue (F) of mice fed either a CD or HFD. Each bar represents data from four mice from each group. (G,H) Rbms1 mRNA expression in adipocyte fraction of visceral adipose of mice fed a CD or HFD for 8 weeks and 20 weeks. The data were obtained through re-analysis of a dataset (GSE 142187), as described in Section 4. Dots represent each value and bars represent mean ± standard Figure 1. Rbms1 expression is significantly altered in adipose tissue depots and adipocytes of mice fed a high-fat diet (HFD). (A-D) Body weight, fat mass, plasma triglyceride levels, and plasma glucose levels in the two groups of mice. (E,F) Rbms1 expression in subcutaneous (E) and mesenteric adipose tissue (F) of mice fed either a CD or HFD. Each bar represents data from four mice from each group. (G,H) Rbms1 mRNA expression in adipocyte fraction of visceral adipose of mice fed a CD or HFD for 8 weeks and 20 weeks. The data were obtained through re-analysis of a dataset (GSE 142187), as described in Section 4. Dots represent each value and bars represent mean ± standard deviation. RT-PCR (n = 4 in each group) and RNA sequencing (n = 6 for each group). ns p > 0.05, ** p < 0.01, *** p < 0.001 and **** p < 0.0001.

RBMS1 Expression Is Downregulated during the Adipogenic Differentiation of 3T3-L1 Cells
3T3-L1 cells are derived from mouse embryos and are commonly used as a model for studying adipocyte differentiation and lipid metabolism. The 3T3-L1 cell line provides a useful tool for studying the molecular mechanisms underlying adipogenesis and for testing potential anti-obesity drugs. We subjected 3T3L1 cells (day 0) to RNA sequencing to assess the expression of the whole repertoire of RBPs in 3T3-L1 cells prior to differentiation (preadipocytes), and Rbms1 was one of the highly expressed RBPs (Supplementary Table S1). We induced the adipogenic differentiation of these cells by incubating them with differentiation cocktail of IBMX, dexamethasone, insulin, and rosiglitazone, followed by post-differentiation with insulin alone. The cells were imaged on day 7 after differentiation to confirm the differentiation process ( Figure 2A). Next, we investigated how adipogenic differentiation affects RBMS1 expression by measuring RNA and protein expression at different points during the differentiation process. Interestingly, RBMS1 was highly expressed before and during the early stages of differentiation, but declined toward the end of the differentiation process ( Figure 2B,C). For reliability and uniformity, we used multiple controls (β-actin, GAPDH, and α-tubulin) to assess RBMS1 protein expression at different stages of differentiation ( Figure 2C). We examined the subcellular localization of RBMS1 in these cells and discovered that the protein was excluded from the nucleus and predominantly expressed in the cytoplasm, consistent with its role in RNA metabolism ( Figure 2D). These findings suggest that RBMS1 may play a role in regulating adipogenesis and adipose metabolism by controlling RNA processing and transport within the cytoplasm of these cells. deviation. RT-PCR (n = 4 in each group) and RNA sequencing (n = 6 for each group). ns p > 0.05, ** p < 0.01, *** p < 0.001 and **** p < 0.0001.

RBMS1 Expression Is Downregulated during the Adipogenic Differentiation of 3T3-L1 Cells
3T3-L1 cells are derived from mouse embryos and are commonly used as a model for studying adipocyte differentiation and lipid metabolism. The 3T3-L1 cell line provides a useful tool for studying the molecular mechanisms underlying adipogenesis and for testing potential anti-obesity drugs. We subjected 3T3L1 cells (day 0) to RNA sequencing to assess the expression of the whole repertoire of RBPs in 3T3-L1 cells prior to differentiation (pre-adipocytes), and Rbms1 was one of the highly expressed RBPs (Supplementary  Table S1). We induced the adipogenic differentiation of these cells by incubating them with differentiation cocktail of IBMX, dexamethasone, insulin, and rosiglitazone, followed by post-differentiation with insulin alone. The cells were imaged on day 7 after differentiation to confirm the differentiation process ( Figure 2A). Next, we investigated how adipogenic differentiation affects RBMS1 expression by measuring RNA and protein expression at different points during the differentiation process. Interestingly, RBMS1 was highly expressed before and during the early stages of differentiation, but declined toward the end of the differentiation process ( Figure 2B,C). For reliability and uniformity, we used multiple controls (β-actin, GAPDH, and α-tubulin) to assess RBMS1 protein expression at different stages of differentiation ( Figure 2C). We examined the subcellular localization of RBMS1 in these cells and discovered that the protein was excluded from the nucleus and predominantly expressed in the cytoplasm, consistent with its role in RNA metabolism ( Figure 2D). These findings suggest that RBMS1 may play a role in regulating adipogenesis and adipose metabolism by controlling RNA processing and transport within the cytoplasm of these cells.

RBMS1 May Play an Important Role in Adipocyte Metabolism
Given the significant expression of RBMS1 before and during the early stages of the adipogenic differentiation of 3T3 L1 cells, we hypothesized that RBMS1 might play a role in the differentiation process. To test this hypothesis, we used shRNA to silence Rbms1 expression in 3T3-L1 cells using a lentiviral transduction system followed by puromycin selection. RBMS1 knockdown (RBMS1KD) was confirmed by immunoblotting with specific antibodies ( Figure 3A). Interestingly, RBMS1 knockdown had no effect on the proliferation and viability of 3T3-L1 cells, as measured using the MTT assay ( Figure 3B). Moreover, cells appeared to differentiate into adipocytes with the formation of lipid droplets similar to scrambled control cells (scr ctrl) ( Figure 3C). The quantification of oil droplets revealed no difference in the lipid content of RBMS1 knockdown cells compared to the scrambled control ( Figure 3D). Furthermore, the mRNA and protein expression of two important markers of adipogenesis, proliferator-activated receptor gamma (PPARγ) and adiponectin, were unaffected by RBMS1 knockdown ( Figure 4A-C). However, several other adipogenic markers such as adipose triglyceride lipase (ATGL), CCAAT/enhancer-binding protein alpha (C/EBPα) and hormone-sensitive lipase (HSL) were significantly downregulated in RBMS1 knockdown 3T3-L1 cells compared to the scrambled control ( Figure 4A-C). C/EBP α is a master transcription factor that is known to regulate gene expression and plays a crucial role in adipogenic differentiation and development. ATGL and HSL are two key enzymes involved in the breakdown of triglycerides into free fatty acids. They play a crucial role in regulating lipid metabolism and energy homeostasis in the body. ATGL is responsible for the initial step of triglyceride hydrolysis, while HSL acts in the second step and is also involved in the breakdown of other lipid molecules. Interestingly, ATGL expression in RBMS1KD cells is only affected at the protein level, implying that RBMS1 may be involved in the post-transcriptional regulation of ATGL and, thus, regulates lipid mobilization and metabolism.

RNA Sequencing and Downstream Analysis Reveals the Role of RBMS1 in Adipocyte Metabolism
To investigate the effect of RBMS1 deletion on adipocyte gene expression, m 3T3-L1 cells expressing scrambled control shRNA or Rbms1 shRNA were subjec RNA sequencing on the Illumina NovaSeq 6000 platform. DESeq2 and the R packag used to obtain the mRNA expression profile for these two groups. A comparison u cut-off value of FDR-adjusted p-value ≤ 0.1 revealed that 430 genes were differentia pressed in RBMS1 knockdown adipocytes as compared to the scrambled control. Th jority of the differentially expressed genes were downregulated (66%), and only 34% upregulated in RBMS1KD cells compared to the control. Below is a heat map disp the relative expression levels of the differentially expressed genes across different sa in the two groups ( Figure 5A), and the relationship between statistical significan fold change for each gene is shown in a volcano plot ( Figure 5B). The top 25 upreg and top 25 downregulated genes are shown in Supplementary Tables S2 and S3.
The downstream analysis was performed using the Ingenuity Pathway An (IPA) program to identify the biological pathways and networks that are enriched L1 adipocytes after RBMS1 knockdown. The analysis revealed that the differentia

RNA Sequencing and Downstream Analysis Reveals the Role of RBMS1 in Adipocyte Metabolism
To investigate the effect of RBMS1 deletion on adipocyte gene expression, mature 3T3-L1 cells expressing scrambled control shRNA or Rbms1 shRNA were subjected to RNA sequencing on the Illumina NovaSeq 6000 platform. DESeq2 and the R package were used to obtain the mRNA expression profile for these two groups. A comparison using a cut-off value of FDR-adjusted p-value ≤ 0.1 revealed that 430 genes were differentially expressed in RBMS1 knockdown adipocytes as compared to the scrambled control. The majority of the differentially expressed genes were downregulated (66%), and only 34% were upregulated in RBMS1KD cells compared to the control. Below is a heat map displaying the relative expression levels of the differentially expressed genes across different samples in the two groups ( Figure 5A), and the relationship between statistical significance and fold change for each gene is shown in a volcano plot ( Figure 5B). The top 25 upregulated and top 25 downregulated genes are shown in Supplementary Tables S2 and S3. 1alpha (HIF1) signaling, AMP-activated protein kinase (AMPK) signaling, insulin receptor signaling, fatty acid oxidation, and many others ( Figure 6A). Additional functional analysis revealed that RBMS1 knockdown altered cellular activities such as lipid metabolism, glucose metabolism, post-translational modification, cell death and survival, cell growth and proliferation, the inflammatory response, protein synthesis, and energy production ( Figure 6B). Our findings imply that RBMS1 regulates a variety of cellular processes in 3T3 L1 adipocytes, including energy metabolism and inflammation.  The downstream analysis was performed using the Ingenuity Pathway Analysis (IPA) program to identify the biological pathways and networks that are enriched in 3T3 L1 adipocytes after RBMS1 knockdown. The analysis revealed that the differentially expressed genes were primarily involved in pathways such as oxidative phosphorylation, mitochondrial dysfunction, glucocorticoid receptor signaling, hypoxia-inducible factor 1alpha (HIF1) signaling, AMP-activated protein kinase (AMPK) signaling, insulin receptor signaling, fatty acid oxidation, and many others ( Figure 6A). Additional functional analysis revealed that RBMS1 knockdown altered cellular activities such as lipid metabolism, glucose metabolism, post-translational modification, cell death and survival, cell growth and proliferation, the inflammatory response, protein synthesis, and energy production ( Figure 6B). Our findings imply that RBMS1 regulates a variety of cellular processes in 3T3 L1 adipocytes, including energy metabolism and inflammation.

Differentially Expressed Proteins Associated with RBMS1 Knockdown in Mature Adipocytes
To identify differentially expressed proteins as a result of RBMS1 depletion, samples were subjected to 2D-DIGE proteomic analysis..Three biological replicates from each group (RBMS1KD and scrambled control) were separated, analyzed, and used for the comparative analysis. A representative protein separation pattern is shown in Supplementary Figure S1. A total of 1206 protein spots were analyzed on the gel, and after applying a cut-off line using a p-value ≤ 0.05 and a fold change ≥ 1.5, 73 significantly different proteins were identified (Supplementary Figure S2). Subsequently, only 37 proteins could be sequenced and identified using the MALDI-TOF sequencing system combined with the SWISS-PROT database and Mascot software 2.6.00 (Supplementary Table S4). Many proteins had their expression altered only at the protein level, not at the mRNA level, indicating potential post-transcriptional regulation by RBMS1. A comparative analysis of these proteins and their mRNA expression is shown in Figure S3.

Differentially Expressed Proteins Associated with RBMS1 Knockdown in Mature Adipocytes
To identify differentially expressed proteins as a result of RBMS1 depletion, samples were subjected to 2D-DIGE proteomic analysis..Three biological replicates from each group (RBMS1KD and scrambled control) were separated, analyzed, and used for the comparative analysis. A representative protein separation pattern is shown in Supplementary Figure S1. A total of 1206 protein spots were analyzed on the gel, and after applying a cut-off line using a p-value ≤ 0.05 and a fold change ≥ 1.5, 73 significantly different proteins were identified (Supplementary Figure S2). Subsequently, only 37 proteins could be sequenced and identified using the MALDI-TOF sequencing system combined with the SWISS-PROT database and Mascot software 2.6.00 (Supplementary Table S4). Many proteins had their expression altered only at the protein level, not at the mRNA level, indicating potential post-transcriptional regulation by RBMS1. A comparative analysis of these proteins and their mRNA expression is shown in Figure S3.
The list of differentially expressed proteins was uploaded into the IPA database for The list of differentially expressed proteins was uploaded into the IPA database for downstream analysis. The top pathways affected by RBMS1 knockdown included unfolded protein response, mitochondrial dysfunction, and oxidative phosphorylation, consistent with the analysis of the RNA sequencing data ( Figure 7A). Additional functional analysis revealed that differentially expressed proteins were significantly enriched in nucleic acid metabolism, energy production, endocrine system disorders, cell cycle, cell death, and lipid metabolism ( Figure 7B). folded protein response, mitochondrial dysfunction, and oxidative phosphorylation, consistent with the analysis of the RNA sequencing data ( Figure 7A). Additional functional analysis revealed that differentially expressed proteins were significantly enriched in nucleic acid metabolism, energy production, endocrine system disorders, cell cycle, cell death, and lipid metabolism ( Figure 7B

Discussion
RBPs have emerged as important regulators of adipose metabolism and possible therapeutic targets for metabolic diseases [37][38][39]. These proteins have a distinct structure that allows them to directly bind to RNA and regulate its degradation, stability, adenylation, transport, and, ultimately, protein expression [40]. In the last few years, several members of the RBP family, including HuR, QK, and YBX1, have been shown to regulate adipogenesis by controlling the stability and translation of key mRNAs involved in this process [21][22][23]. RBMS1 is a member of the RNA-binding protein family and contains two conserved RNA binding motifs, RNP1 and RNP2. RBMS1 has been extensively studied for its role in cancer and found to regulate various cellular processes such as alternative splicing, mRNA stability, and translation, making it a potential therapeutic target for cancer treatment [28,29]. Additionally, RBMS1 expression levels have been shown to correlate with cancer progression and patient survival in certain types of cancer. Multiple studies have found that single nucleotide polymorphisms (SNPs) in RBMS1 are linked to the risk of type 2 diabetes mellitus and other obesity-related traits [37,38,41]. Moreover, Rbms1 mRNA expression shows a significant increase in the placenta from obese mothers as compared to lean mothers [42]. These observations strongly suggest a possible role of RBMS1 in metabolic regulation.
To the best of our knowledge, no prior research has looked into the role of RBMS1 in adipogenesis or adipose tissue metabolism. Therefore, the current study was initiated to investigate the potential role of RBMS1 in adipocyte differentiation and metabolism. We found that RBMS1 is highly expressed in the early stages of differentiation of the murine adipocyte cell line 3T3-L1 and decreases in the later stages. Furthermore, Rbms1 expression is significantly higher in the subcutaneous and visceral adipose tissue depots of HFD-fed mice as compared to CD-fed mice. Diet-induced obesity (DIO) is known to affect the cellular composition of adipose tissue, which may have a direct influence on the expression of genes such as Rbms1. To demonstrate that the increased Rbms1 expression in response to DIO is not attributable to changes in adipose tissue cellular composition, we re-analyzed an RNA sequencing dataset containing data from adipocytes extracted from visceral adipose tissue of HFD-and CD-fed mice. Consistent with our results from whole adipose tissue, Rbms1 expression was significantly upregulated in the adipocyte fraction of adipose tissue in response to a high-fat diet. This suggests that the upregulation of Rbms1 expression in response to a high-fat diet is not due to changes in adipose tissue cellular composition, but rather a direct response to DIO. We further showed that RBMS1 may play a crucial role in regulating the differentiation of 3T3-L1 adipocytes. The depletion of RBMS1 using specific shRNA resulted in the downregulation of some of the key adipogenic markers, including ATGL and HSL, two important lipase enzymes that regulate lipid mobilization and metabolism. Intriguingly, ATGL expression in RBMS1KD cells was only altered at the protein level, not at the mRNA level, suggesting that RBMS1 may control ATGL post-transcriptionally. ATGL is a key enzyme in the conversion of stored fat into energy, and any disruption in its regulation can result in metabolic disorders such as obesity and diabetes [43,44]. More research is needed to determine whether RBMS1 regulates ATGL translation and, if so, how. The transcriptomic and proteomic analyses further revealed that the depletion of RBMS1 led to altered expression of genes involved in key metabolic pathways, including those involved in carbohydrate and lipid metabolism. Interestingly, the mRNA expression of several differentially expressed proteins was unaffected by RBMS1 knockdown, suggesting a possible post-transcriptional regulation by RBMS1. More studies are needed to confirm these findings and to understand the mechanism by which RBMS1 controls the translation of these mRNAs.
Taken together, these findings suggest that RBMS1 may play an important role in adipose metabolism by regulating multiple metabolic pathways. More research is needed to fully understand how RBMS1 regulates adipose metabolism and the adipogenesis process. First, these findings must be replicated in animal models, preferably by creating adiposespecific RBMS1 knockout mice. Second, identifying potential RBMS1 downstream target mRNAs could provide valuable insight into the molecular mechanisms underlying its metabolic effects. Finally, identifying specific compounds capable of modulating RBMS1 expression or activity could be a promising avenue for drug development.

Cell Culture and Differentiation Induction
The 3T3-L1 mouse embryo fibroblast pre-adipocytes were purchased from the American Type Culture Collection (ATCC), (Manassas, VA, USA). Cells were maintained and differentiated in Dulbecco's modified Eagle's medium (Sigma, St. Louis, MO, USA) supplemented with 10% heat-inactivated fetal bovine serum (Thermo Fisher Scientific, Waltham, MA, USA) and 10,000 U/mL penicillin-streptomycin (Thermo Fisher Scientific, Waltham, MA, USA), and incubated at 37 • C and 5% CO 2 as described previously. Briefly, cells were grown until confluency then starved for 2-3 days before initiating the differentiation process. Differentiation was induced by incubating cells with the complete medium (DMEM + 10% FBS) containing the adipogenic cocktail (1 µM insulin, 0.25 µM dexamethasone, 0.5 mM IBMX, 3 µM Rosiglitazone). The differentiation medium was replenished after 3 days with complete medium +1 µM insulin for 24 h. Differentiated cells were maintained in complete medium and used on day 7, as described before [45].

Animals and Diet
C57Bl/6J male mice were a kind gift from Dr. Jahangir Iqbal of the King Abdullah International Medical Research Center in AL-hasaa, SA. Mice were maintained and housed at a 22 ± 0.5 • C temperature, as mentioned previously [46]. Mice were fed ad libitum either a chow diet CD (10% energy by fat, catalog #D12450B, Research Diets, Inc., New Brunswick, NJ, USA) or a high-fat diet HFD (60% energy by fat, catalog #D12492, Research Diets, Inc., New Brunswick, NJ, USA). Four mice in each group were fed the indicated diets for 18 weeks. Weekly weight measurements were recorded, and blood triglycerides, glucose, and free fatty acids were measured. At the end of the treatment, mice were euthanized, and all fat pads were removed and weighed. Experiments were managed with the approval of the KAIMRC Institutional Animal Care and Use Committee (Protocol #s RA17-013-A and RA20-005-A).

Tissue Digestion, RNA Isolation, and RT-PCR
An RNeasy mini kit was purchased from Qiagen, (Germantown, MD, USA) and utilized to isolate RNA from cell lysates. Pre-adipocytes and mature adipocytes were collected, and total RNA was extracted following the manufacturer's protocol. Meanwhile, RNA from fat tissue was extracted by utilizing Trizole-chloroform following the manufacturer's protocol, as described previously. In short, an equal amount of total RNA was used to synthesis cDNA using a High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Waltham, MA, USA) and was mixed with Power SYBR Green PCR Master Mix (Applied Biosystems, Foster City, CA, USA). The amplification reaction was performed in an Applied Biosystems QuantStudio (TM) 6 Flex real-time PCR System. Primers are listed in the Supplementary Materials (Table S6).

Plasmid, Lentivirus Preparation, and Generation of Stable RBMS1 Knockdown
3T3-L1 knockdown cells were generated using lentiviral transduction and Rbms1specific short hairpin RNA (shRNA). Lentiviral particles were packaged in human embryonic kidney 293T (HEK293T) cells following transfection with the relevant shRNA and lentiviral packaging mix. 3T3-L1 cells were plated in 6-well plates in full medium at 40-50% confluency, and lentiviral particles with transduction enhancer (Polybrene 8 µg/mL) were added to each well. The medium containing viral particles was removed after 48 h of incubation and replaced with fresh medium containing 2 µg/mL puromycin to generate stable RBMS1KD in 3T3-L1 cells.

Library Preparation and RNA Sequencing
RNA was extracted from freshly collected cell pellets as described above. Only samples with a high-quality RNA integrity (RIN ≥ 9) on an Agilent 2100 bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) were judged to be suitable for RNA library preparation prior to whole-transcriptome sequencing. Three biological replicates from each group were utilized for RNA sequencing at the King Abdullah International Medical Research Center's medical genomics laboratory (Riyadh, Saudi Arabia). In short, the library was created using the Illumina Stranded Total RNA Prep protocol. After obtaining purified total RNA at a concentration of 1000 ng per sample, ribosomal RNA (rRNA) was enzymatically depleted using an Illumina Ribo-zero plus rRNA depletion kit, followed by fragmentation and random hexanucleotide treatment to prime the sample for cDNA synthesis. RNA fragments primed with hexamer were then used to create double-stranded cDNA. To prevent them from ligating each other, adenylate 3' ends were added, and a ligate anchor was added to prepare a cDNA fragment for dual indexing primers (i7, i5) using RNA prep ligation, IDT DNA/RNA UD indexes, and RNA index anchor kits (New England BioLabs, Ipswich, MA, USA). The library was amplified using PCR and cleaned with Agencourt AMPureXP (Beckman, CA, USA). The final library was then normalized, pooled, and loaded into an SP flow cell for up to 400 million reads of coverage. On an Illumina NovaSeq 6000 platform, the library was sequenced (200 cycles).

RNA Sequencing Data Processing
The generated raw sequence reads were trimmed using Trim Galore (version 0.6.7) [47] to remove adapter sequences and low-quality bases. Quality assessment was performed using FastQC (version 0.11.9) to evaluate the overall quality of the reads, including measures such as per-base sequence quality scores, sequence length distribution, and GC content. Processed reads were aligned to the Ensembl reference genome GRCm38 using STAR (version 2.7.10a) [48], a spliced transcripts aligner, with default parameters. Gene-level counts were generated from the aligned reads using Salmon (version 1.9.0) [49], a fast and accurate transcript quantification tool, with default parameters. Differentially expressed genes (DEGs) between the study groups were identified using DESeq2 (version 1.38.0) [50], an R package for RNA-seq data analysis. DESeq2 was used to calculate the fold change, p-value, and false discovery rate (FDR) for each gene, with the threshold for significance set at an FDR-adjusted p-value < 0.05. Batch effects were investigated using principal component analysis (PCA). Logarithmic scale (base 2) fold change was used to determine the magnitude of gene expression changes. DEG lists for each group were uploaded to Ingenuity Pathways Analysis (IPA), a bioinformatics tool for functional mapping and pathway analysis. IPA was used to identify enriched biological functions, pathways, and networks associated with the DEGs, using default settings and a significance threshold of p-value < 0.05. All statistical analyses were performed using R (version 4.2.1) unless otherwise stated. Heat maps were generated using the integrated Differential Expression and Pathway analysis (iDEP.96) platform [51]. TBtools 1.098 software was used to extract data of the entire repertoire of RNA-binding proteins and to generate volcano plots from RNA sequencing data [52]. All raw sequencing data has been deposited to Gene Expression Omnibus record number GSE233256.

RNA Sequencing Re-Analysis
The GSE142187 dataset was analyzed using the GEO RNA-seq Experiments Interactive Navigator (GREIN) platform [53]. Dataset GSE142187 contained RNA expression data from adipocytes isolated from the white adipose tissue of mice fed either a regular chow or high-fat diet for 8 and 20 weeks [36]. CPM-normalized mRNA expression was obtained from the GREIN platform, and differential gene expression was estimated. Raw p values were adjusted for multiple testing using the Benjamini-Hochberg procedure and only adjusted p (Padj) values lower than 0.05 were considered significant.

Cellular Localization through Confocal Microscopy
Pre-and post-differentiated 3T3-L1 cells were plated in glass-bottom dishes, washed with PBS, and fixed with 4% paraformaldehyde. Cells were permeabilized by incubating with PBS containing 0.1% Triton X-100 and incubated in blocking solution (1% BSA, 22.52 mg/mL glycine) in PBST (PBS + 0.1% Tween 20). Cells were then incubated with RBMS1 antibody for 1 h at room temperature followed by washing and incubation with goat anti-rabbit IgG (H+L) highly cross-adsorbed secondary antibody for 1 h in the dark. 40,6-diamidino-2-phenylindole (DAPI, Cell Signaling Technology, Danvers, MA, USA) was used to stain nuclei in cell cultures. The cells were imaged using an Eclipse TE2000 confocal laser scanning microscope (CLSM) (Nikon Corporation, Tokyo, Japan).

Measurement of Cell Growth and Viability
To investigate the impact of the depletion of RBMS1 on the proliferation and viability of 3T3 L1 cells, a 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay was performed. Briefly, cells were seeded at a density of 2 × 10 3 cells/well in 96well plates and incubated at 37 • C in 5% CO 2 . Then, 5 uL of 0.5 mg/mL MTT (Thermo Fisher Scientific, Waltham, MA, USA) was added to each well after 24, 48, and 72 h of incubation and they were incubated for 4 h in 5% CO 2 at 37 • C. To dissolve the formazan crystals, 100 uL DMSO (Sigma, St. Louis, MO, USA) was used, and optical density (OD) at 562 nm was measured using a SpectraMax ® Microplate Reader (Molecular Device, San Jose, CA, USA).

Oil Red O Stain and Lipid Quantification
The impact of RBMS1KD on lipid accumulation was measured through oil red O staining and lipid quantification. Briefly, cells were seeded in 12-well plates and were differentiated after they reached confluency, as described before [45]. On day 7, media were removed, and cells were washed and fixed in 4% paraformaldehyde methanol free for 1 h. A freshly made oil red O stain was filtered through a 0.2 µm filter and used to stain cells for 20 min at room temperature. Cells were rinsed with 60% isopropanol for 5 min to remove any unwanted background. The cells were imaged on a phase contrast microscope. To quantify the lipid droplets, 100% isopropanol was used to extract the stain from the droplets, and absorbance was measured at 562 nm using a SpectraMax ® Microplate Reader (Molecular Device, San Jose, CA, USA).

Western Blotting
Proteins from cell pellets were extracted using RIPA buffer (Thermo Fisher Scientific, Waltham, MA, USA) supplemented with a protease inhibitor cocktail (Thermo Fisher Scientific, Waltham, MA, USA). Then, 15 µg protein was separated through 10% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to a nitrocellulose membrane (0.2 µm, Bio-Rad Laboratories, Philadelphia, PA, USA) using a wet transfer system at 60 V for 1 h at 4 • C. Membranes were blocked using 5% fat-free milk and probed with the indicated primary antibodies, followed by the suitable anti-rabbit/mouse horseradish peroxidase (HRP)-conjugated secondary antibodies. Protein bands were developed and visualized using a ChemiDoc XRS imager (Bio-Rad) with Image Lab software (version 4.1, Bio-Rad).

Protein Extraction and Labelling for Two-Dimensional Electrophoresis and Image Scanning
Three separated replicates from each experimental group were used for proteomics analysis at the Obesity Research Center, College of Medicine at King Saud University (Riyadh, Saudi Arabia). Briefly, mature adipocytes were lysed and sonicated in buffer (1.5 mL, pH 8.8, 30 mM Tris buffer containing 7 M urea, 2 M thiourea, 2% Chaps, 1× protease inhibitor mix), as described previously [54]. In brief, after the sonication, 50 mM dithiothreitol (DTT) was added to stabilize the proteins, followed by 40 min centrifugation at 20,000× g at 4 • C. Supernatant was used and then precipitated using a 2D clean-up kit (GE Healthcare, USA) according to the manufacturer's protocol. A buffer (7 M urea, 2 M thiourea, 4% CHAPS, and 30 mM Tris) was prepared and utilized to re-suspend the protein pellets, followed by protein pH adjustment to 8.5. Protein concentrations were measured in triplicate using a 2D-Quantkit (GE Healthcare, Danderyd, Sweden). After that, proteins were labeled with CyDye™ DIGE Fluor minimal dye (400 pmol/50 µg) (GE Healthcare, Danderyd, Sweden) according to the manufacturer's recommendation. Each biological replicate was labeled with either cy3 or cy5. Moreover, a pooled standard of 50 µg protein from each of the 6 samples was labeled with cy2. Accordingly, samples were combined and separated by first dimension using Immobiline Dry Strips (24 cm, pH 4-7; GE Healthcare, Danderyd, Sweden). These separated proteins were subsequently separated using second-dimension sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) (Ettan DALT six vertical units, GE Healthcare, Danderyd, Sweden), scanned with a Sapphire Biomolecular Imager (Azure Bio systems, Dublin, OH, USA), and digitalized via the image analysis software Sapphire Capture system Version 1.12.0921.0 (Azure Biosystems, Dublin, OH, USA). In order to identify the statistically significant protein spots, cy2 labeling (internal standard) was included to enable normalization across gels and quantitative differential analysis of the protein levels. An ANOVA parameter was used to determine the significant changes in protein abundance by applying a cut-off line p-value ≤ 0.05 and fold change ≥ 1.5. Subsequently, a total protein of 1 mg was created by pooling equal protein amounts from 6 samples to prepare a preparative gel (3 knockdown, 3 controls). The gel was stained for five days with Coomassie stains (Thermo Fisher Scientific, Waltham, MA, USA), followed by a gentle rinse with distilled water. The stained gel was kept at 4 • C until the protein spots could be selected, picked, and identified using MS.

MALDI-TOF Mass Spectrometry and Data Processing
Statistically significant protein spots were collected manually, washed, and digested according to previously described methods [55,56]. Trypsin was used to cleave protein peptides at 37 • C overnight. Next, 0.8 µL of tryptic peptides was loaded into a MALDI target plate (384 MTP Anchorchip, 800 µm Anchorchip, Bruker Daltonics, Bremen, Germany). MALDI-MS/MS spectra were measured on an UltraflexTerm TOF mass spectrometer connected to a LIFT-MS/MS device (Bruker Daltonics). MS data were translated using BioTools v3.2 (Bruker Daltonics). Each peptide sequence was aligned to the database to identify proteins by utilizing the Mascot engine (v2.0.04, updated on 9 May 2021; Matrix Science Ltd., London, UK). Only those proteins that showed a Mascot score greater than 56 with p ≤ 0.05 were considered for Ingenuity Pathway Analysis (IPA).

Pathways Analysis and Functional Ontology
The functional analysis and pathway analysis were performed through the use of IPA (QIAGEN Inc. https://www.qiagenbioinformatics.com/products/ingenuity-pathwayanalysis (accessed on 8 March 2023)). Only differentially expressed genes and proteins were uploaded for core analysis. In short, a list of DE genes containing the ratio-transformed log 2 fold change and adjusted p-value for each gene was imported, and core analysis was run with consideration of fold change less than or equal to 1.5 and p-value 0.1 and 0.05 for proteomics. Bio-functional analysis was carried out to characterize their functions.

Statistical Analysis
Each experiment was carried out three independent times for two experimental groups. A total of at least 3 technical/biological replicates were statistically analyzed using the twotailed unpaired t-test or analysis of variance (ANOVA) using Prism GraphPad 9 software (GraphPad Software, Inc., La Jolla, CA, USA). Only values that showed a p-value less than or equal to 0.05 were considered significant.

Conclusions
In conclusion, the current study identifies RBMS1 as an important regulator of adipogenesis and adipose metabolism. These preliminary findings must be confirmed in animal models, and the use of adipose-specific RBMS1 knockout mice could provide valuable insight into the functional significance of this protein.

Informed Consent Statement: Not applicable.
Data Availability Statement: The datasets generated during the current study has been deposited to the Gene Expression Omnibus record number GSE138242.