Abstract

The hypothalamus plays a central role in the integrated regulation of feeding and energy homeostasis. The hypothalamic arcuate nucleus (ARC) contains a population of neurons that express orexigenic and anorexigenic factors and is thought to control feeding behavior via several neuronal circuits. In this study, a comparative proteomic analysis of low-fat control diet- (LFD-) and high-fat diet- (HFD-) induced hypothalamic ARC was performed to identify differentially expressed proteins (DEPs) related to changes in body weight. In the ARC in the hypothalamus, 6621 proteins () were detected, and 178 proteins were categorized as DEPs (89 upregulated and 89 downregulated in the HFD group). Among the Gene Ontology molecular function terms associated with the DEPs, protein binding was the most significant. Fibroblast growth factor receptor substrate 2 (Frs2) and SHC adaptor protein 3 (Shc3) were related to protein binding and involved in the neurotrophin signaling pathway according to Kyoto Encyclopedia of Genes and Genomes analysis. Furthermore, high-precision quantitative proteomic analysis revealed that the protein profile of the ARC in mice with HFD-induced obesity differed from that in LFD mice, thereby offering insight into the molecular basis of feeding regulation and suggesting Frs2 and Shc3 as novel treatment targets for central anorexigenic signal induction.

1. Introduction

Obesity is a common issue worldwide and a cause of various disorders, including hypertension, diabetes, and cardiovascular disorders [1]. Several neural circuits in the brain, including the hypothalamus, solitary nucleus, midbrain, and limbic system, are involved in the maintenance of body weight and appetite [2]. Leptin and insulin regulate feeding behavior and energy homeostasis in the hypothalamus [3, 4], and the arcuate nucleus (ARC) of the hypothalamus is involved in feeding and energy expenditure [5].

The ARC is anatomically located at the lateral wall of the third ventricle and lies directly above the pituitary gland upon the median eminence in the hypothalamus. Several types of neurons in the feeding circuit of the ARC have been discovered, including those expressing agouti-related protein/neuropeptide Y (agouti-related protein [AgRP]/NPY) or proopiomelanocortin (POMC) derivatives [6, 7]. Orexigenic AgRP neurons are negatively regulated by adipocyte-derived leptin, whereas anorexigenic POMC neurons are triggered by leptin to release mainly α-MSH derivatives that affect nuclei, including the paraventricular nucleus (PVN) of the brain [8, 9]. Aside from several peptides and neural pathways, the molecular mechanisms underlying the regulation of the feeding circuit in the ARC have not been fully elucidated. Considering its pivotal role in food intake and satiety, we hypothesized that proteome profiling of the ARC would uncover critical circuits related to obesity. Although proteomic analyses of the hypothalamus in a high-fat diet- (HFD-) induced obesity (DIO) model have been performed [9], few studies have focused specifically on the ARC.

In this study, we characterized the ARC proteome in order to identify protein alterations in mice exposed to an HFD and provide new insights into the molecular basis of appetite regulation and energy metabolism.

2. Materials and Methods

2.1. Animals

Male C57BL/6 mice (7 weeks old; 20–25 g; ) were purchased from Orient Bio (Seongnam, Korea) and maintained at under a 12 h light/dark cycle starting at 8:00 A.M., with ad libitum access to the designated type of food and water in a breeding room. Body weights were measured weekly.

2.2. Ethics Statement

This study was carried out in accordance with the Guide for the Care and Use of Laboratory Animals of Seoul National University Hospital. The protocol was approved by the Institutional Animal Care and Use Committee of the Seoul National University Hospital (protocol 18-0154-S1A1).

2.3. DIO Model

All mice underwent a week of adjustment to the housing environment and were fed a chow diet. The mice were then randomly assigned to two groups (/group) for the 8-week experimental period. One group was fed a low-fat control diet (LFD; rodent diet with 10% of energy from fat and 3.8 kcal/g), and the other was fed an HFD (60% of energy from fat and 5.2 kcal/g) ad libitum for 8 weeks (D12450B and D12492, respectively; Research Diets, Brunswick, NJ, USA). Mean body weights were measured weekly in both groups to evaluate the reliability of the DIO model.

2.4. Sample Preparation

After 8 weeks, each mouse was anesthetized with a ketamine-xylazine mixture (80–120 mg/kg) and intracardially perfused with sterilized cold saline and 4% paraformaldehyde. Brains were extracted, snap-frozen over liquid nitrogen, and stored at −80°C. A punch biopsy was performed, as described previously, to excise ARC tissue samples from the brain [10, 11], and the samples were deep-frozen immediately after sampling.

2.5. Protein Quantification

Each ARC sample was collected, as described previously [10, 11], and left and right samples were pooled and subjected to BCA/Bradford peptide assay for peptide extraction before quantification. Peptide sequences in the extracts were profiled in depth using a quadrupole-orbitrap mass spectrometer (Q Exactive Plus; Thermo Fisher Scientific, Waltham, MA, USA). Peptides were identified using the MaxQuant database (Max Planck Institute of Biochemistry, Munich, Germany) and assigned with a [12]. The level of peptide expression was calculated using intensity-based absolute quantification (iBAQ) [13]. Peptides identified in at least one sample were included in the final list of valid proteins, which was imported to the Perseus platform (Max Planck Institute of Biochemistry) for statistical analyses and detection of differentially expressed proteins (DEPs).

2.6. Functional Pathway Analysis

The list of DEPs was uploaded via the STRING 11.0 tool to the STRING database (http://string-db.org) to generate a protein-protein interaction network (ELIXIR Core Data Resource) [14], and DEPs were searched against the Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.kegg.jp/kegg/) and Gene Ontology (GO) (http://geneontology.org/) databases to evaluate functional enrichment. The overall characteristics of DEPs were examined by GO enrichment analysis [15] using the Biological Networks Gene Ontology (BiNGO) plug-in of Cytoscape [16]. Biochemical pathways involving two or more DEPs were determined [17].

2.7. Statistical Analysis

Bodyweight data are shown as and analyzed with two-way repeated measures analysis of variance implemented in SPSS (v.24.0; IBM Corp., Armonk, NY, USA).

For proteomic analysis for each protein, iBAQ intensities were transformed to log2 values, and missing values were imputed based on the Gaussian distribution. We then performed width adjustment (a process to equalize the quartiles of intensity values between the samples) for normalization. Student’s -tests were used for comparisons between the normalized iBAQ intensities of the HFD and LFD groups, and fold change values (defined as the ratio of the normalized mean iBAQ intensities of group HFD to group LFD) were then calculated. A DEP was defined as a protein with a and either a (upregulated protein) or <1.0/1.5 (downregulated protein). A volcano plot was generated for the DEPs based on values and fold changes.

3. Results

3.1. Mouse Model of HFD-Induced Obesity

After unrestricted administration of an HFD for 8 weeks, a DIO model was successfully established. From 3 weeks (LFD vs. HFD: vs. , ) to 8 weeks (LFD vs. HFD: vs. , ) after the adjustment period, the mean weight in the HFD group was significantly higher than that in the LFD group (Figure 1(a)). After 15 weeks, the HFD group showed a significantly higher body weight than the LFD group () (Figure 1(b)).

3.2. Protein Quantification

Left and right ARC tissue samples were pooled and used for peptide extraction and quantification. In total, 6,621 proteins were identified in the LFD and HFD groups. A volcano plot was generated to visualize the proteins and display the distribution and significant proteins (determined by −log10 values for the mean normalized iBAQ intensity). In total, 336 proteins were identified at a threshold value of 0.05. After filtering with the fold change threshold, 178 proteins were identified as DEPs, including 89 upregulated DEPs () and 89 downregulated DEPs () (Figure 2(a), Table S1). Hierarchical clustering revealed a distinct proteomic signature in the hypothalamus under HFD and LFD. The normalized iBAQ intensities of the 178 DEPs and 6,621 valid proteins for each sample (five samples each for the LFD and HFD groups) were visualized using a heatmap (Figures 2(b) and 2(c)).

3.3. iBAQ Analysis of the Neurotrophin Signaling Pathway

We evaluated a candidate pathway (the neurotrophin signaling pathway) based on previous results demonstrating its important roles in neurogenesis, axogenesis, and synaptic plasticity in the hypothalamus and feeding and energy imbalances due to obesity [18, 19]. Among the DEPs in our analysis, we detected two proteins (fibroblast growth factor (FGF) receptor substrate 2 (Frs2) and SHC adaptor protein 3 (Shc3)) involved in the neurotrophin signaling pathway (Figure 3(a)). The intensity-based quantification confirmed that Frs2 and Shc3 levels were lower in the HFD group than in the control group (Figures 3(b) and 3(c)).

3.4. KEGG Pathway Analysis

The DEPs were further evaluated by KEGG pathway analysis, with a KEGG pathway diagram generated after deleting unrelated molecules and redundant data. Among 45 DEPs related to “protein binding,” five KEGG pathways that included ≥3 DEPs related to this term were further explored. The threshold was set to three DEPs, because the number of KEGG pathways was high () when the threshold was set to two DEPs, which increased the risk of false-positive results. The five selected KEGG pathways and their corresponding genes were as follows: mitogen-activated protein kinase (MAPK) signaling pathway (Mapk14, Map3k12, and Rps6ka4), thermogenesis (Mapk14, Frs2, Prkab2, and Uqcrh), proteoglycans in cancer (Mapk14, Frs2, and Lum), Salmonella infection (Mapk14, Arpc5, and Wasf2), and neurotrophin signaling pathway (Mapk14, Frs2, and Shc3). The roles of the MAPK signaling pathway and the thermogenesis pathway are already clearly established; therefore, these pathways were excluded from further analyses. In addition, proteoglycans in cancer and Salmonella infection pathways were excluded, because they were highly unlikely to have therapeutic value in obesity. Accordingly, the neurotrophin signaling pathway (including Mapk14, Frs2, and Shc3) was selected as a key pathway for further analyses (Figure 4). The modified pathway diagram suggests that these proteins might suppress axonal outgrowth (Figure 5).

3.5. GO Functional Enrichment Analysis

According to GO enrichment analysis and the BiNGO network (Tables 13), DEPs were enriched for various GO terms and particularly in two categories: biological processes and cellular components (Figures S1 and S2). BiNGO enabled graphical visualization of several DEPs and values for each GO term (Figure 4). The results showed that enriched terms in the molecular function category included “binding” (, 90 DEPs) and catalytic activity (, 47 DEPs) (Figure 4). Compared with catalytic activity, “binding” included more DEPs and was more highly significant. Among subgroups of the “binding” term, “protein binding” (, 45 DEPs) was the most highly significant and included the most DEPs. The DEPs related to “protein binding” were further evaluated for involvement in KEGG pathways, revealing one pathway related to the hypothalamic ARC (i.e., the neurotrophin signaling pathway) (Figure 4).

4. Discussion

This study identified DEPs in the ARC between mice with HFD-induced obesity and mice fed an LFD. Glutamate transporters (e.g., GLT and GLAST) are among the exceptionally abundant and important proteins in the central nervous system (~>2.5% of the brain protein mass exists in the form of glutamate transporters) and critical for the function of glutamatergic signaling (the balance between glutamate excitation and GABA inhibition is particularly important and relevant for body mass control mechanisms in the hypothalamus). Although proteomic changes in the hypothalamus of an obesity model have been studied [9], our comparison of the ARC in LFD and HFD groups was novel. Moreover, in-depth profiling of the global proteome targeting the ARC using the punch biopsy method revealed various proteins that have not previously been reported in studies of obesity or energy metabolism.

Recent studies suggest that tumor necrosis factor α and interleukin 1β, especially in AgRP-producing neurons, might contribute to inflammation-induced anorexia during acute inflammatory conditions [20]. Additionally, single-cell RNA sequencing of the hypothalamic ARC after HFD-induced obesity demonstrated selective changes in AgRP neurons via neuron-astrocyte interactions, which contributed to the exaggerated sympathoexcitation observed in obese rats [21, 22]. These effects might be mediated by the regulatory effects of leptin on astrocytic glutamate transporters within the ARC of the hypothalamus [23].

GO enrichment analysis identified binding and catalytic activity as key molecular functions. Notably, the DEPs related to protein binding were involved in the neurotrophin signaling pathway. These results were consistent with those of previous reports indicating that the neurotrophin signaling pathway contributes to neurogenesis, axogenesis, and synaptic plasticity in the hypothalamus and dysfunctions in feeding and energy balance due to obesity [18, 19]. “Binding” was a highly significant term and included a number of DEPs, and the subgroup “protein binding” was further identified as a significant molecular function. DEPs associated with “protein binding,” particularly Frs2 and Shc3, were involved in the neurotrophin signaling pathway. This pathway is generally involved in retrograde transport, cell survival, cellular differentiation, and axonal outgrowth of neurons and triggered by the binding of extracellular signaling molecules, such as nerve growth factor, brain-derived neurotrophic factor (BDNF), and neurotrophins 4 and 3. Furthermore, a previous study reported that the distribution of tropomyosin receptor kinase B- (TrkB-) expressing neurons is altered in response to an HFD in the mediobasal hypothalamus [24].

These results showed that Frs2 and Shc3 levels and the related neurotrophin pathway were altered in the ARC of a mouse model of HFD-induced obesity. Neurotrophins, which mediate receptors and signaling molecules, play multiple roles in survival and development [25]. FRS2 functions by transmitting FGF and/or neurotrophin signaling during dorsal forebrain development, with a critical role in the formation of the hippocampal dentate gyrus [26]. Additionally, ARCTrkB neurons regulate food intake by projecting to BDNF-expressing neurons in the PVN [27].

A previous study reported that BDNF deficiency signaled a reduced number of axonal projections of ARCTrkB to PVN, subsequently inducing hyperphagia [28]. SHC molecules and neurotrophins bind to Trks and the neurotrophin receptor [25, 29]. In the present study, the DEPs included Frs2 and Shc3, which represent novel candidate pharmaceutical targets for obesity based on their close links to pathways involved in obesity. Furthermore, these two proteins identified by in-depth profiling of the ARC using the punch biopsy method have not previously been reported in studies of obesity or energy metabolism. However, the central effects of Frs2 or Shc3 remain unknown. To fully understand their roles in the regulation of energy expenditure and food intake, further studies are needed, such as an intracerebroventricular injection study along with analyses of hourly and daily food intake and metabolic parameters.

5. Conclusions

We identified 178 DEPs in the ARC of a reliable mouse model of HFD-induced obesity. These DEPs play key roles in protein binding and the neurotrophin signaling pathway and include Frs2 and Shc3 as novel candidates involved in appetite and energy metabolism.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Chang Yeon Kim and Jang Ho Ahn contributed equally to this study.

Acknowledgments

We thank Ki Soon Dan and Hye-Yeon Kim from the Proteomics & Biomarker Lab (Seoul National University College of Medicine, Seoul, Korea) for their contributions to the proteomic analysis. We also thank Min Sun Kim for providing technical assistance. This work was supported by the National Research Foundation of Korea (NRF) Grant funded by the Korean Government (MSIT) (grant no. NRF-2018R1A5A2025964) and Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (grant no. 2017R1A6A3A1131814) and the Seoul National University Hospital Research Fund (grant no. 26-2016-0030).

Supplementary Materials

Supplementary 1. Figure S1: GO pathway—biological process. Figure S2: GO pathway—cellular component.

Supplementary 2. Table S1: DEPs and iBAQ intensity values in the HFD group relative to the LFD group.