Hypothalamus proteomics from mouse models with obesity and anorexia reveals therapeutic targets of appetite regulation.

OBJECTIVE: This study examined the proteomic profile of the hypothalamus in mice exposed to a high-fat diet (HFD) or with the anorexia of acute illness. This comparison could provide insight on the effects of these two opposite states of energy balance on appetite regulation.
METHODS: Four to six-week-old male C56BL/6J mice were fed a normal (control 1 group; n=7) or a HFD (HFD group; n=10) for 8 weeks. The control 2 (n=7) and lipopolysaccharide (LPS) groups (n=10) were fed a normal diet for 8 weeks before receiving an injection of saline and LPS, respectively. Hypothalamic regions were analysed using a quantitative proteomics method based on a combination of techniques including iTRAQ stable isotope labeling, orthogonal two-dimensional liquid chromatography hyphenated with nanospray ionization and high-resolution mass spectrometry. Key proteins were validated with quantitative PCR.
RESULTS: Quantitative proteomics of the hypothalamous regions profiled a total of 9249 protein groups (q<0.05). Of these, 7718 protein groups were profiled with a minimum of two unique peptides for each. Hierachical clustering of the differentiated proteome revealed distinct proteomic signatures for the hypothalamus under the HFD and LPS nutritional conditions. Literature research with in silico bioinformatics interpretation of the differentiated proteome identified key biological relevant proteins and implicated pathways. Furthermore, the study identified potential pharmacologic targets. In the LPS groups, the anorexigen pro-opiomelanocortin was downregulated. In mice with obesity, nuclear factor-κB, glycine receptor subunit alpha-4 (GlyR) and neuropeptide Y levels were elevated, whereas serotonin receptor 1B levels decreased.
CONCLUSIONS: High-precision quantitative proteomics revealed that under acute systemic inflammation in the hypothalamus as a response to LPS, homeostatic mechanisms mediating loss of appetite take effect. Conversely, under chronic inflammation in the hypothalamus as a response to HFD, mechanisms mediating a sustained 'perpetual cycle' of appetite enhancement were observed. The GlyR protein may constitute a novel treatment target for the reduction of central orexigenic signals in obesity.

Introduction

Overweight and obesity as a result of positive energy balance constitute major public health burdens with significant economic and social implications.[1] The hypothalamus is the key brain site for the regulation of food intake and energy balance in mammals. Under physiological conditions, a variety of peripheral signals regulate appetite and adjust energy intake to match energy consumption requirements.[2] Systemic acute inflammatory signals can cause profound anorexia by disrupting the physiological regulation of appetite in the hypothalamus.[3] Conversely, emerging evidence suggests that hypothalamic ‘inflammatory' activation as a result of a high-fat diet (HFD) and obesity can disturb anorexigenic and thermogenic signals and promote abnormal body weight control.[4] Deciphering the molecular events of these two contradictory observations in a global and directly comparative way allows for a more causal understanding on the relations between diet, inflammatory signals and appetite regulation. Provided that nutritional intervention protocols are likely to interfere with multiple pathways inside the cell, a system-wide interrogation of the host response using non-targeted quantitative proteomics was warranted.[5]

The aim of this study was to compare the global proteomic profile of the hypothalamus in mice exposed to a HFD or with anorexia of acute illness. We hypothesise that such a comparison can shed new insight on the effects of these two opposite states of energy balance on appetite regulation.

Materials and methods

Animal model

Studies were conducted in male mice on C56BL/6J background. Mice were bred and maintained at regular housing temperatures (23 °C) and 12-h light/dark cycle starting at 0700 hours. Animals had ad libitum access to water and food and were weighed at weekly intervals. Four to six-week-old mice (n=34) were randomly divided into four groups: the control 1 (C1) group (n=7), the HFD group (n=10), the control 2 (C2) group (n=7) and the lipopolysaccharide (LPS) group (n=10). Mice in the C1 group were fed a control diet (4.5% fat, 34% starch, 5.0% sugar and 22.0% protein), whereas the HFD group was fed a HFD (24% fat, 41% carbohydrate, and 24% protein; Research diets D12451 formula) for 8 weeks after which they were anaesthetised and perfused intracardially with phosphate-buffered saline. A diet containing 24% fat (45% kcal fat) was selected because it efficiently induces obesity and better simulates a human HFD as opposed to a more extreme 60% kcal fat diet.[6, 7] Mice in the C2 and LPS groups were fed a control diet for 8 weeks. After this period, mice in the C2 group received an intraperitoneal injection of 0.3 ml saline, whereas in the LPS group received an intraperitoneal injection of LPS (Escherichia coli O111:B4 (Sigma No. L-2630)) at a dose of 120 μg per mouse dissolved in 0.3 ml saline. Mice in the C2 and LPS groups were killed 18 h after the injection as described above. As protein expression in the hypothalamus is known to be under circadian regulation,[8] mice were killed at 0700 hours. Experimental procedures were approved by the Institutional Animal Care and Use Committee at the Centre of Basic Research, Biomedical Research Foundation of the Academy of Athens.

Determination of insulin resistance by HOMA-IR

HOMA-IR analysis was used to assess insulin resistance in HFD-fed mice. After overnight fasting, values for homeostasis model assessment of insulin resistance (HOMA-IR) were calculated from the values of fasting serum glucose (mg dl−1) and fasting serum insulin (μU ml−1) by using the following formula: HOMA-IR=fasting glucose value (mg dl−1) × fasting insulin value (μU ml−1)/405. Low HOMA-IR values indicate high insulin sensitivity, whereas high HOMA-IR values indicate low insulin sensitivity (insulin resistance). Fasting glucose concentrations were measured using a hand-held glucometer (Biorad, Hercules, CA, USA), whereas serum insulin levels were quantified by ELISA (Mercodia, Uppsala, Sweden).

Quantitative proteomics sample processing

Hypothalamic regions were removed and snap frozen at −80 °C. Specimens were dissolved in 0.5 m triethylammonium bicarbonate, 0.05% sodium dodecyl sulphate, homogenised using the Fast Prep system (Savant Bio, Cedex, France) followed by pulsed probe sonication (Misonix, Farmingdale, NY, USA). Lysates were centrifuged (16 000 g, 10 min, 4 °C) and supernatants were measured for protein content using infra-red spectroscopy (Merck Millipore, Darmstadt, Germany). For the quantitative proteomic analysis, the hypothalamic regions from three mice were used for the C1 and C2 groups and from six mice for the HFD and LPS groups (n=18 in total). Three individual protein extracts were pooled (33.3 μg from each lysate giving 100 μg final protein content) to form one sample for the C1 and C2 groups and two biological replicates for the HFD and LPS groups. Lysates were then reduced, alkylated and subjected to trypsin proteolysis. Peptides were labelled using six of the eight-plex iTRAQ reagent kit (113=C1, 114=C2, 115=HFD1, 116=HFD2, 117=LPS1, 118=LPS2) and analysed using high-precision two-dimensional liquid chromatography with nanospray ionization tandem mass spectrometry as reported previously by the authors[9, 10, 11, 12] (Figure 1a) (Supplementary Methods 1).

Figure 1

(a) Overview of the proteomics pipeline (b) Venn diagram description of the differentiated proteins between the HFD and LPS groups relative to their respective controls. (c) Hierarchical clustering analysis of modulated proteins: the HFD and LPS groups have a distinct endophenotypic portrait in the hypothalamus. (d) Gene ontology analysis using BiNGO of the commonly modulated proteins in the HFD and LPS conditions compared with their respective control showed a significant enrichment for acute inflammatory response. Coloured nodes representing GO terms correspond to those that were significant according to the P-value colour scale, whereas white nodes were not significant. The size of the node reflects the number of proteins that mapped to the corresponding GO term.

Database searching

Unprocessed raw files were submitted to Proteome Discoverer 1.4 for target decoy searching against the UniProtKB/TrEMBL mus musculus database comprised of 57 475 entries (release date 03 September 2014), allowing for up to two missed cleavages, a precursor mass tolerance of 10 ppm, a minimum peptide length of six and a maximum of two variable (one equal) modifications of; iTRAQ 8-plex (Y), oxidation (M), deamidation (N, Q) or phosphorylation (S, T, Y). Methylthio (C) and iTRAQ (K, N terminus) were set as fixed modifications. Reporter ion ratios from unique peptides only were taken into consideration for the quantitation of the respective protein. Quantification ratios were median-normalised and log2 transformed. Proteins were grouped using the protein grouping inference algorithm in the Proteome Discoverer 1.4 application (Supplementary Methods 2). A protein was considered modulated in the HFD or LPS group relative to control when its log2 ratio was above or below±1 s.d. across both biological replicates.[13] In adherence to the Paris Publication Guidelines for the analysis and documentation of peptide and protein identifications (http://www.mcponline.org/site/misc/ParisReport_Final.xhtml) and the recently published guidelines for high-confidence protein identification by Omenn et al.,[14] only proteins identified with at least two unique peptides were further subjected to bioinformatics analysis. Proteins of biological relevance, two of which were identified with one unique peptide, were validated using targeted quantitative PCR. Proteomics data were deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD003475.

Bioinformatics analysis

Heatmap construction of differentially expressed proteins in the HFD and LPS groups compared with their respective control was generated using Cluster 3.0 (http://bonsai.hgc.jp/~mdehoon/software/cluster/software.htm) and Java Treeview (http://jtreeview.sourceforge.net). MetaCore (GeneGo, St Joseph, MI, USA) and BiNGO were applied to differentially expressed proteins analysed with at least two unique peptides to identify over-represented processes in the modulated proteome of each group compared with their respective control. FDR-corrected P-values <0.05 were considered significant.

Total RNA isolation and cDNA synthesis

For the quantitative PCR analysis the hypothalamic regions from four mice were utilised per group (C1, C2, HFD, LPS) (n=16 in total). Total RNA from mouse hypothalamus was extracted using an RNA Isolation Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. RNA concentration and purity were determined by NanoDrop 2000c Spectrophotometer (Thermo Scientific, Waltham, MA, USA). Ribosomal RNA band integrity was evaluated by conventional 1% agarose gel electrophoresis and the Agilent Bioanalyzer with the RNA 6000 Nano Kit (Agilent, Santa Clara, CA, USA).

For complementary DNA (cDNA) synthesis, 300 ng total RNA from each sample was reverse transcribed into cDNA using M-MLV Reverse Transcriptase (Thermo Scientific) according to the manufacturer's instructions. All the cDNA samples were stored at −20 °C until quantitative PCR analyses.

Quantification of mRNA

Quantification of mRNA was performed for the following proteins: NFκ-B, neuropeptide Y (NPY), 5-hydroxytryptamine (serotonin) receptor 1B, glycine receptor alpha-4 subunit, pro-opiomelanocortin. The PCR mixture contained 1 μl diluted cDNA, 10 μM gene-specific primer (forward and reverse mixed together) and 10 μl of 2 × Fast SYBR Green Master Mix (Roche Diagnostics, Rotkreuz, Switzerland) in a total volume of 20 μl. Amplification was performed in 96-well optical reaction plates (Roche Diagnostics) on LightCycler 480 (Roche Diagnostics) using the following programme: 94 °C for 3 min to activate polymerase, 40 cycles at 94 °C for 20 s, 60 °C for 20 s and 72 °C for 20 s; melting curve analysis was performed after every run by heating up to 95 °C to monitor presence of unspecific products. Two negative controls were included in each assay run, with water instead of template. Three replicate measurements for each sample were performed. Primers were designed and checked with Primer Quest Tool (IDT) and NCBI primer BLAST tool and synthesised by Macrogen (Seoul, South Korea). Primer sequences are listed in Supplementary Table 1.

Quantitative PCR data analysis

The mRNA expression of genes in the hypothalamus of LPS-treated and HFD-fed mice was calculated relative to the expression in their respective control mice, according to the delta–delta Ct method (2−ΔΔCt) using the most and least stable reference genes found, as well as the most commonly used glyceraldehyde-3-phosphate dehydrogenase and beta-actin.

Results

Mice fed a HFD weighed significantly more than those fed a control diet (Supplementary Figure 1A). Furthermore, they were insulin resistant with significantly higher HOMA-IR values (Supplementary Figure 1B). The proteomic analysis of the hypothalamic regions resulted in the profiling of 9249 protein groups (q<0.05) (Supplementary Tables 2 and 3), of which 7718 were identified with at least two unique peptides. Among the proteins identified with at least two unique peptides, 201 were upregulated in the HFD and 193 in the LPS groups compared with their respective control, of which 30 were common between the two conditions. Furthermore, 169 proteins were downregulated in the HFD and 194 in the LPS groups, of which 18 were common between the two conditions (Figure 1b, Supplementary Tables 4–6). The R2 value between biological replicates was 0.69 and 0.89 for the differentially expressed proteins of the HFD and LPS groups, respectively. Hierarchical clustering of the differentiated proteome revealed a distinct proteomic signatures of the hypothalamus under the HFD and LPS (Figure 1c) nutritional conditions.

Gene ontology analysis using BiNGO of the 48 commonly modulated proteins showed a significant enrichment for acute inflammatory response (Figure 1d). In the hypothalamus of the LPS-exposed mice protein folding and maturation_posttranslational processing of neuroendocrine peptides (FDR-corrected P-value=1.2E-6) (Figure 2a) and apoptosis and survival_caspase cascade (FDR-corrected P-value=1.8E-2) (Figure 2b) were significantly enriched.

Figure 2

In silico analysis using MetaCore showed that (a) protein folding and maturation_posttranslational processing of neuroendocrine peptides (FDR-corrected P-value=1.2E-6) and (b) apoptosis and survival_caspase cascade (FDR-corrected P-value=1.8E-2) were significantly enriched in the differentially expressed proteins of the LPS groups compared with control. Analysed proteins are denoted with a circle (red=upregulation, blue=downregulation).

In the differentially expressed proteins of the HFD groups, inflammation_IL-6 signalling (FDR-corrected P-value=6.5E-6) (Figure 3a), inflammation_kallikrein-kinin system (FDR-corrected P-value=9.3E-6), and inflammation_protein C signalling (FDR-corrected P-value=9.2E-3) were significantly over-represented. A nodal protein in all three inflammatory pathways was NF-κB, analysed to increase in the HFD groups compared with control (Figure 1c). Gene ontology analysis using BiNGO confirmed that both acute and chronic inflammatory responses were significantly over-represented in the modulated proteome of the hypothalamus in HFD-fed mice compared with control (Figure 3b).

Figure 3

(a) In silico analysis using MetaCore showed that inflammation_IL-6 signalling was significantly enriched in the differentially expressed proteins of the HFD groups compared with control (FDR-corrected P-value=6.5E-6). Analysed proteins are denoted with a circle (red=upregulation, blue=downregulation). (b) Gene ontology analysis using BiNGO confirmed that acute and chronic inflammatory responses are significantly enriched in the differentially expressed hypothalamic proteins of HFD-fed mice compared with control. Coloured nodes representing GO terms correspond to those that were significant according to the P-value colour scale, whereas white nodes were not significant. The size of the node reflects the number of proteins that mapped to the corresponding GO term.

Of potential biomedical interest to be discussed, nuclear factor-κB (NF-κB), pro-NPY and glycine receptor alpha-4 subunit levels increased, whereas serotonin receptor 1B levels were reduced in the HFD groups relative to control. In the LPS groups compared with control, pro-opiomelanocortin was profiled to be downregulated. The relative quantitation of these proteins at the mRNA level was validated using quantitative PCR (Figure 4).

Figure 4

Validation of key proteomic findings using qPCR: NF-κB, pro-neuropeptide Y and glycine receptor alpha-4 subunit levels increased, whereas serotonin receptor 1B levels were reduced in the HFD groups relative to control. In the LPS groups compared with control, pro-opiomelanocortin was profiled to be downregulated.

Discussion

Molecular perturbations in the hypothalamus as a result of obesity have been studied primarily at the transcriptome but not the proteome level.[15, 16] This study compared the global proteomic profile of the mouse hypothalamus between two experimental groups in opposite states of energy balance, that is, mice with increased caloric intake through a HFD, versus mice with severely compromised food intake due to acute systemic inflammation. Acute inflammatory response was significantly enriched in the commonly modulated proteins of both groups (Figure 1d), in agreement with a number of studies indicating the altered ‘inflammatory' state in these experimental models.[3, 17]

It is well established that LPS induces anorexia in mouse models.[18] Our results showed that the acute systemic inflammatory response caused by LPS treatment specifically decreased pro-opiomelanocortin (Figure 4) and the related ACTH and alpha-MSH peptides with well-shown anorexigenic effects (Figure 2a). This paradoxical downregulation of anorexigenic signals possibly highlights the presence of a negative feedback loop, aiming to re-establish homoeostasis. Furthermore, hypothalamic sampling was done 18 h after LPS administration, as the mice were recovering from the endotoxin effects, as the administered LPS dose was at a medium range. Along this theme, reduced expression levels of caspase isoforms 3, 6 and 7 in the hypothalamus of LPS-treated mice was observed (Figure 2b). This finding suggests control of apoptotic processes and protection of neurons during re-establishment of homoeostasis, supporting the theory that inflammatory and homoeostatic control mechanisms may antagonise in a time- and immune activation phase-dependent manner.[19]

In the modulated proteome of the HFD groups, a number of inflammatory pathways, including the IL-6 pathway (Figure 3a), kallikrein-kinin system and protein C signalling, were significantly enriched, reflecting a more pronounced, chronic inflammatory response, also evident in the gene ontology analysis (Figure 3b). Upregulation of pro-inflammatory mediators in the hypothalamus is a hallmark of obesity and constitutes a major component of insulin resistance pathogenesis.[20] The exact mechanisms mediating this phenotype remain a subject of intense investigation with endoplasmic reticulum stress, oxidative stress or direct activation of inflammatory-related transcription factors among the implicated processes.

A central node in the hypothalamic inflammatory processes was NF-κB, found to be upregulated as a result of HFD (Figure 4). NF-κB, an important mediator of metabolic inflammation, is enriched in hypothalamic neurons but normally remains inactive.[21] Over-nutrition alone, even before the onset of obesity, activates hypothalamic NF-κB partly through increased endoplasmic reticulum stress.[21] Furthermore, activation of hypothalamic NF-κB interrupts central insulin and leptin signalling and actions, whereas local suppression in the mediobasal hypothalamus, or more specifically in hypothalamic Agouti-related peptide neurons, significantly protects against obesity and glucose intolerance.[21]

NPY was profiled to increase in the hypothalamus of HFD-fed mice compared with control (Figure 4). Central administration of NPY, a well-established orexigenic peptide in the hypothalamus,[22] has been found to increase food intake and body weight, with chronic administration leading to development of obesity.[23, 24] Furthermore, NPY overexpression in the dorsomedial hypothalamus induces hyperphagia and obesity,[25, 26] whereas knockdown of NPY improves these changes.[27, 28] A recent study showed that dorsomedial hypothalamus NPY knockdown rats on HFD had normal food intake, body weight and glucose tolerance/insulin sensitivity indexes as observed in lean control rats.[28] These authors suggest that dorsomedial hypothalamus NPY could be a potential therapeutic target for obesity and diabetes.

In the hypothalamus of HFD-fed mice compared with control, levels of serotonin receptor 1B (5-hydroxytryptamin receptor 1B; 5-HT1BR), a key modulator of energy homeostasis and a therapeutic target for obesity, were analysed to decrease (Figure 4). Serotonin reduces the activity of the orexigenic AgRP/NPY neurons via the 5-HT1BR in mice (or 5-HT1DBR in humans).[29] Studies have shown that mice lacking the 5-HT1BR develop hyperphagia, resulting in body weight increase.[30] In addition, administration of CP-94,253, a highly selective 5-HT1BR agonist, significantly decreases food intake[31] and potentiates the anorectic effects of pro-opiomelanocortin.

GlyR was analysed to be upregulated in the hypothalamus of mice with obesity compared with control (Figure 4). Although glycine is a well-established inhibitory neurotransmitter in the spinal cord and brain stem,[32] studies examining the function of glycine receptors in higher brain structures including the hypothalamus are limited. A study by Karnani et al.[33] showed that glycine receptor agonists exert a hyperpolarizing action on mature hypocretin/orexin (hcrt/orx) neurons in the hypothalamus, important regulators of arousal and reward seeking behaviour. To the best of our knowledge, our study is the first to identify an increase of GlyR levels in the hypothalamus of mice exposed to HFD. GlyR receptor antagonists, could be tested for their efficacy to reduce central orexigenic signals in obesity.

Conclusion

High-precision quantitative proteomics revealed that under acute systemic inflammation in the hypothalamus as a response to LPS, homeostatic mechanisms mediating loss of appetite take effect. Conversely, under chronic inflammation in the hypothalamus as a response to HFD, mechanisms mediating a sustained ‘perpetual cycle' of appetite enhancement were observed. The combinatorial protein level modulation of the NF-κB, 5-HT1BR, GlyR and pro-NPY provides strong evidence for the enhancement of orexigenic signals as a response to HFD. Our study is the first to identify an increase of GlyR levels in the hypothalamus of mice exposed to HFD. Antagonists for this receptor could be tested for their efficacy to reduce central orexigenic signals in obesity.