Skip to main content

Increased protein intake affects pro-opiomelanocortin (POMC) processing, immune function and IGF signaling in peripheral blood mononuclear cells of home-dwelling old subjects using a genome-wide gene expression approach

Abstract

Background

Adequate protein intake among older adults is associated with better health outcomes such as immune function and metabolic regulation of skeletal muscle, but conflicting results make it difficult to define the optimal intake. To further understand the impact of protein intake on metabolic processes, the aim of the study was to explore genome-wide gene expression changes in peripheral blood mononuclear cells (PBMCs) in home-dwelling old subjects after increased protein intake for 12 weeks.

Method

In a parallel double-blind randomized controlled intervention study, subjects (≥ 70 years) received a protein-enriched milk (2 × 20 g protein/day, n = 14, mean (±SD) age 76.9 ± 4.9 years) or an isocaloric carbohydrate drink (n = 17, mean (±SD) age 77.7 ± 4.8 years) for breakfast and evening meal for 12 weeks. PBMCs were isolated before and after the intervention. Microarray analysis was performed using Illumina technology. Serum levels of gut peptides and insulin growth factor (IGF)-1 were also measured.

Results

In total 758 gene transcripts were regulated after increased protein intake, and 649 gene transcripts were regulated after intake of carbohydrates (p < 0.05). Forty-two of these genes were overlapping. After adjusting for multiple testing, 27 of the 758 gene transcripts were regulated (FDR, q-value < 0.25) after protein intake. Of these 25 were upregulated and two downregulated. In particular, genes and signaling pathways involved in pro-opiomelanocortin (POMC) processing, immune function, and IGF signaling were significantly altered.

Conclusions

PBMCs can be used to study gene expression changes after long-term protein intake, as many signaling pathways were regulated after increased protein intake. The functional significance of these findings needs to be further investigated.

Trial registration

ClinicalTrials.gov, ID no. NCT02218333. The study was registered on August 18, 2014.

Background

Adequate intake of dietary protein is important to maintain muscle mass and muscle strength in older adults. Age-related loss of skeletal muscle mass, skeletal muscle strength and functional performance is linked to a number of poor health outcomes, such as impaired functionality, reduced ability to manage activities of daily life [1, 2], reduced quality of life [3], increased morbidity, but also mortality [4,5,6]. Epidemiological studies have indicated that dietary protein intake may be a modifiable risk factor for loss of skeletal muscle mass and skeletal muscle strength (sarcopenia) in older adults [7]. However, no consensus has been reached related to recommendations of optimal intake level of protein in older adults, and results from studies investigating changes in muscle mass or muscle strength after increased protein intake have been conflicting [6, 8,9,10,11].

Dietary interventions usually cause only small changes within the body and may be difficult to detect by changes in the phenotype, such as muscle mass, muscle strength or circulating inflammatory markers, especially in short-term trials [12]. Other, more sensitive methods, such as changes in gene transcripts, may detect important changes at a molecular level. These signatures may serve as biomarkers in a life-long perspective, affecting homeostatic control and the risk of developing life style related diseases [12, 13]. Protein intake and certain amino acids are known to regulate gene expression [14, 15]. The most studied pathway regulated by dietary proteins is probably the mammalian target of rapamycin (mTOR) pathway, which is involved in the regulation of muscle protein synthesis [16]. However, proteins and amino acids are also shown to play an important role in the regulation of general control nonderepressible 2 (GCN2)/activating transcription factor 4 (ATF4), which regulate protein metabolism, but may also affect lipid and glucose metabolism [14].

Peripheral blood mononuclear cells (PBMCs) can be used as a model system for studying changes in gene expression levels in dietary intervention studies [17, 18]. PBMCs include monocytes, lymphocytes and natural killer (NK) cells mainly; cells that play an important role in inflammation and in the development of chronically related diseases, such as cardiovascular diseases [19, 20]. PBMCs are circulating cells exposed to nutrients, metabolites and peripheral tissues, and PBMCs may therefore reflect whole-body health [18]. However, their role in aging and age-related conditions, such as loss of muscle mass and the decline of immune function with aging is not well described. Furthermore, how gene expression in PBMCs is affected by dietary protein intake has hardly been investigated in older subjects. We recently showed that consumption of protein-enriched milk had minor effects on a limited number of selected inflammatory genes and genes involved in muscle mass in PBMCs from older adults, when using a targeted approach [21]. In the present study we further elucidated the impact of protein intake on metabolic processes. The aim of our study was to explore genome-wide gene expression changes in PBMCs in home-dwelling old subjects with reduced physical strength and/or performance after 12 weeks with increased protein intake.

Results

Subject characteristics

There were no statistically significant differences between the study groups in either of the parameters measured at baseline, such as BMI, lean body mass, hs-CRP or the amounts of monocytes and lymphocytes (Table 1).

Table 1 Baseline characteristics

The dietary intake was calculated based on two 24-h recall interviews performed prior to, and during the study period (Table 2). As expected, the intake of protein increased in the protein group, while the intake of carbohydrates and fiber were reduced. In the carbohydrate group the intake of protein (E %), saturated fat, and monounsaturated fat were reduced, and the intake of carbohydrates increased.

Table 2 Dietary changes during the intervention in the protein group and the carbohydrate group

Gene expression profiling in PBMCs

Microarray hybridization was performed on mRNA from PBMCs collected at baseline and after 12 weeks of intervention from all subjects included. From the 48,000 probe sets present at the HumanHT-12 v4 microarray chip, 12,135 unique gene transcripts were expressed in PBMCs in the present study. In total 758 gene transcripts were regulated after intake of protein intake, and 649 gene transcripts were regulated after intake of carbohydrates (p < 0.05) (Additional file 1 : Table S1 and Additional file 2: Table S2, respectively), with approximately equal number up- and downregulated gene transcripts (Fig. 1). Only 42 gene transcripts were overlapping between the two groups (Additional file 3: Table S3) (Fig. 1). Among the overlapping genes 19 and 22 gene transcripts were upregulated after intake of protein and carbohydrates, respectively, while 23 and 20 gene transcripts were downregulated in the respective groups (Additional file 3: Table S3). Nine gene transcripts were regulated in opposite direction in the two groups (Fig. 1) (Additional file 3: Table S3).

Fig. 1
figure 1

Overview of number of genes regulated after intervention. In total 48,000 probes were on the HumanHT-12 v4 microarray chip. 12,135 were defined as expressed in PBMCs. 758 gene transcripts were changed after protein intake, while 649 gene transcripts were changed after carbohydrate intake (p < 0.05). 42 gene transcripts were overlapping in both groups. After adjusting for multiple testing, 27 gene transcripts were significantly changed after protein intake, while seven were significantly changed after carbohydrate intake (FDR, q-value < 0.25). Four gene transcripts were overlapping after adjusting for multiple testing

After adjusting for multiple testing (FDR, q-value < 0.25), the expression of 27 gene transcripts were changed from baseline to end of the intervention in the protein group (Table 3), and seven gene transcripts were changed in the carbohydrate group (Table 4).

Table 3 Genes significantly regulated after intake of protein (FDR, q < 0.25)
Table 4 Genes significantly regulated after intake of carbohydrates (FDR, q < 0.25)

Among the regulated gene transcripts (FDR, q-value < 0.25), four genes were overlapping in the two groups. In both groups, the expression level of KDM5A was significantly downregulated, whereas the expression levels of RPL23A, ZNF543 and MICA/B were upregulated (Table 3 and 4). Additional gene transcripts, altered in the protein group, were primarily involved in the immune response, such as KIR2DLs, KLRC3 and CCL4L2 (Table 3).

Four of the most upregulated genes (KIR2DL1, KIR2DL4, KLRC3, and CCL4L2) in the protein group (FDR q-value < 0.25) (Table 3) were selected for confirmation by RT-qPCR analyses. The mRNA level of KIR2DL1 (p = 0.04) was significantly upregulated after protein intake. The mRNA levels of KLRC3 and KIR2DL4 were also upregulated after protein intake, but not statistically significantly from baseline (p = 0.08, and p = 0.67, respectively) (Fig. 2). Although the differences were not significant, the median effect was in the same direction as in the microarray analyses. In contrast to the microarray results, we were not able to detect any expression level of CCL4L2 with the RT-qPCR method.

Fig. 2
figure 2

qPCR validation of microarray findings after the 12-week intervention of protein-enriched milk. Change of mRNA expression (log ratio) from baseline to end of study within the protein group. The box represents the 25-75 percentiles (IQR), and the horizontal line shows the median log ratio. The upper and lower whiskers show the largest and smallest values no further than 1.5 x IQR and data beyond the end of the whiskers are outliers and are plotted individually. Differences between the time points were tested with a paired Wilcoxon-Mann-Whitney test. Significantly regulated genes are marked with * (p < 0.05)

No difference in the distribution of monocytes and lymphocytes were observed during the intervention, as previously reported (data not shown) [21].

Pathway analysis in MetaCore

To understand further the impact of increased protein intake on metabolic processes, pathway analysis was performed using MetaCore. All differently regulated genes with a nominal p value < 0.05 were imported into the software. We identified 96 differently regulated pathways in the protein group (FDR, q-value < 0.05) (Additional file 4: Table S4). Pathways related to protein folding and maturation of pro-opiomelanocortin (POMC) processing, immune response in NK cells and development of IGF-1 receptor signaling were among the most ten regulated pathways (Fig. 3). Genes involved in the folding and maturation of POMC processing are among others CAP-Gly Domain Containing Linker Protein (CLIP), joining peptide (JP), the adrenocorticotrophic hormone (ACTH), in addition to several forms of Melanocyte-Stimulating Hormones (MSH) and lipotrophin (LPH). In the present study CLIP1 and POMC were significantly downregulated, while CLIP4 were upregulated, in the protein group (nominal p < 0.05) (Additional file 1: Table S1). Pathways were then grouped according to function, e.g. all pathways having “IGF”, “mTORC” and “growth” factor in their names were grouped into one group, and all pathways having “immune” in their names were grouped into the immune response group. Pathways having both “immune” and “apoptosis/survival/signaling transduction” in their names were grouped into the immune response group, others into the apoptosis/survival group and the signal transduction group, respectively. Genes involved in lipid metabolism, and POMC-signaling were grouped into the lipid metabolism, and POMC-signaling groups, respectively (Table 5).

Fig. 3
figure 3

Ten most regulated pathways after the 12-week intervention of protein-enriched milk. -log(pValue) indicates the level of significance, the ratio the number of genes regulated in our sample material compared to the number of genes identified in this specific pathway

Table 5 Regulated pathways (FDR, q-value < 0.05) after intake of protein

Furthermore, 62 differently regulated pathways were identified in the carbohydrate group (FDR, q-value < 0.05). These were related to apoptosis and survival of TNFA-induced caspase-8-signaling, apoptosis and survival of ceramides signaling pathways and IL16 signaling pathways, and are shown in Additional file 5: Table S5 and (flow chart given in Additional file 6: Figure S1).

IGF-1 in serum and gut peptides in plasma

Based on the results from the pathway analysis in the protein group, we analyzed blood concentrations of IGF-1 and gut peptides. Serum IGF-1 level increased from baseline to the end of the intervention in the protein group and was close to being significant (p = 0.08). No change was observed in the level of IGF-1 in the group receiving the isocaloric carbohydrate drink. The change between the groups was also close to being statistically significant (p = 0.07) (Table 6). Circulating levels of GIP, GLP-1, PYY, amylin and PP did not differ within or between the groups after 12 weeks (Table 6).

Table 6 Effects of protein-enriched milk and isocaloric carbohydrate drink on of IGF-1 and gut peptides

Discussion

We investigated genome-wide gene expression changes in PBMCs in older men and women (≥ 70 years) with reduced physical strength and/or performance, before and after a 12-week intervention of increased protein intake. We found that pathways related to protein folding and maturation of POMC processing, immune response in NK-cells and development of IGF-1 receptor signaling were the most regulated pathways after increased protein intake.

POMC is found in many tissues, among them leucocytes [22] and PBMCs [23], being the precursor of several molecules, among them ACTH, b-LPH [22], and α-MSH [23]. POMC-derived peptides are generally known to play an important role in regulating energy homeostasis [24] hunger and satiety [23, 25]. POMC-derived peptides may also exert different effects in different tissues, as post-translational processing events of POMC are common [22, 23]. The functional significance of POMC expression in lymphocytes is unclear, but it is assumed that it forms part of a biochemical loop linking the immune, nervous and endocrine systems [22]. To our knowledge, there has been no studies identifying how diets high in protein may affect the mRNA expression of POMC and POMC related genes in PBMCs. In the present study, CLIP1 and 4, in addition to POMC were significantly regulated (before adjusting for multiple testing), supporting the notion that intake of protein modulated this pathway. High-protein diets have been shown to promote satiety and reduce calorie intake [26] through anorexigenic gut peptides (cholecystokinin, GLP-1, and PYY), but also through hypothalamic pathways involving POMC [27]. Kinzig and colleagues showed that a diet high in protein in rats resulted in significantly increased POMC gene expression, in the hypothalamus [28]. Increased intake of leucine, present in high amounts in dairy products, is an important activator of mTORC1 [29] and has also been shown to increase the mRNA expression of POMC [27].

Among the statistically significantly upregulated genes observed after increased protein intake were Killer Cell Immunoglobulin-like receptors (KIRs) and Killer Cell Lectin Like Receptor C3 (KLRC3) which are closely related to the immune response of TYRO protein tyrosine kinase binding protein (DAP12) receptors in NK cells, as was one of the most regulated pathways within the protein group. These results indicate changes in the activity of NK-cells in the protein group [30].

Seven altered pathways related to IGF and growth hormones signaling after intake of protein were identified. Both IGF and mTORC1 pathways were regulated in PBMCs after the increased protein intake. The regulation of the mTORC1 pathway in PBMC indicated that PBMCs could be used as a model system to study gene expression changes in protein intervention studies, as the most studied pathway regulated by dietary proteins is probably the mTORC pathway [16]. The biological function of the mTORC regulation in PBMCs needs further investigation as this pathway is mainly studied in the metabolic regulation of skeletal muscle. IGF-1 serum levels were increased, but not statistically significant (p < 0.08), indicating that increased protein intake might affect metabolism and growth hormone signaling. Previous studies have shown that increased energy intake, intake of essential amino acids [31, 32], and increased intake of milk can increase circulating levels of IGF-1 [33,34,35,36]. IGF-1 is an anabolic hormone that has an important function in maintaining skeletal muscle mass across all ages [37], and concentrations of circulating IGF-1 have been shown to decline during aging [31] [38]. The level of IGF-1 and chronic low-grade inflammation may be closely linked [39]. Chronic low-grade inflammation seems to a robust predictor of disability and mortality, even in the absence of clinical disease [40, 41]. Dysregulations of the immune system, including failure in resolving inflammation, may play a role in the etiology and perpetuation of underlying inflammation [42] and sarcopenia [43]. In the present study, we observed that many genes related to the function of NK-cells were upregulated after increased protein intake and several pathways related to immune functions were altered by protein intake.

The major strengths of the present study were the double-blind randomized controlled design with the use of an isocaloric test drink and several data analysis strategies to explore their possible effects. The major limitation is that we cannot rule out that the effect of increased protein intake in the protein group, and the increased carbohydrate intake in the carbohydrate group is solely due to these changes since the subjects had other macronutrient changes during the intervention. Furthermore, we used commercially qPCR with primers that were not identical to the probes used in the microarrays. This may be one reason why we were not able to verify all gene transcripts regulated after protein intake.

Conclusions

We identified significant changes in gene transcripts and signaling pathways in PBMCs after increased protein intake. Most changes were related to protein folding and maturation of POMC processing, immune response in NK-cells and IGF-1 receptor signaling. Whether these changes in whole genome transcriptome profiles and PBMCs can affect long-term health outcomes by increased protein intake in older adults needs to be further validated.

Methods

Subjects and study design

The current study was part of a previously published parallel double-blind, randomized controlled intervention trial, conducted from 2014 to 2015 at Oslo and Akershus University College of Applied Sciences, Norway [9]. Fifty home-dwelling men and women (≥ 70 years) with reduced physical strength and/or performance were found eligible for the study. Invitation letters were sent to 2820 subjects, 438 subjects met to the screening visit of which 388 did not meet the inclusion criteria. In total, 50 subjects were randomized, and 36 subjects concluded the study. Inclusion criteria were either reduced grip strength (< 20 kg in women and < 30 kg in men), gait speed < 1 m/s, time step stair test ≥ 8.4 s or timed five times to sit to stand test > 12.5 s, and willingness to keep the physical activity level stable throughout the study period. Subjects with type I and II diabetes or HbA1 ≥ 6.5%, severe inflammation, chronic obstructive pulmonary disease, high blood pressure (> 180/105 mmHg), acute cardiovascular disease within the last 6 months or a history of cancer the last 3 years were excluded. Subjects with thyroid-stimulating hormone outside reference range (0.2–10 mU/l) were only included if the thyroxine concentration was within the reference value. If a stable dosage of thyroxine treatment, hormone therapy, and antihypertensive drugs had been used prior to inclusion, these drugs were allowed during the study. Subjects were stratified by gender and smoking and allocated within each of the two groups to consume either protein-enriched milk (2 × 0.4 l/day, 2 × 20 g protein/day, 5.0% protein, 4.6% carbohydrates, < 0.1% fat, 167 kJ (39 kcal)/100 g) or an isocaloric carbohydrate drink (2 × 0.4 l/day) for breakfast and evening meal for 12 weeks. The test drinks were identical in color, labeling, and appearance, and were blinded to both the study participants and the study staff. The drinks (0.4 l) were labeled with each participant’s ID-number, and information about day and time to be consumed (morning or evening). The test drinks were picked up by the participants at the study center or delivered to their homes by the study staff. Unopened containers and leftovers were returned to the staff. None of the participants was excluded from the study due to low compliance (≥ 70%). Participants were encouraged to maintain their normal diet and physical activity levels during the study period. The participants registered all physical activity ≥ 30 min during the study period. All subjects completed the study within 12 ± 1 weeks.

Study products

The protein-enriched drink and the isocaloric carbohydrate drink were produced and provided by TINE SA, Oslo, Norway. The protein-enriched milk is commercial available for sale in Norway, but was not enriched with vitamin D when used in the study. The protein-enriched drink provided on average 167 kJ (39 kcal), 5.0% protein, 4.6% carbohydrates, < 0.1 g fat/100 g. About 80% of the milk protein was casein and the remaining protein was whey protein. The isocaloric, non-nitrogenous control drink was prepared from carbohydrates (sugar, xantan gum and maltosweet®). Calcium was added to the control drink to match the content in the protein-enriched milk, and titan dioxide was added to give the control drink a milky appearance.

Dietary assessments were performed by two 24-h-dietary recalls prior to baseline and two at the end of the intervention [9], reflecting the diet prior to, and during the intervention period. The interviews were performed using an in-house data program (KBS version 7.0) and was linked to the Norwegian food composition table. Dietary supplements were included in the analysis.

Body composition was measured using dual energy X-ray absorptiometry (DXA). The analysis was performed after an overnight fast (≥ 12 h) at the Norwegian School of Sport Sciences, Oslo, Norway, and we used enCORE Software (version 14.10.022, GE Lunar) to estimate lean muscle mass in arms and legs. Muscle mass was defines as the sum of lean body mass of the four limbs [9].

A detailed description of the protocol, participant requirement, and complete list of exclusion criteria, composition of the test drinks, how the 24-h-dietary recalls were performed and compliance are described previously [9].

All subjects provided written informed consent, and the study was conducted according to the Declaration of Helsinki. We received approval for all procedures involving human subjects by the Regional Committees for Medical and Health Research Ethics, Health Region South East, Norway. The study was registered at Clinicaltrials.gov (ID no. NCT02218333).

Blood sampling and preparation

Blood samples were drawn in the morning after an overnight fast (≥ 12 h) in BD Vacutainer® CPTTM cell preparation tubes with sodium heparin (Becton Deckenson, NJ, USA) at baseline and after 12 weeks. Using CPTTM cell preparation tubes is a well-documented and standardized method to collect mononuclear cells with high purity (above 90 %), and according to the manufacturer, approximately 80% of the cells are lymphocytes and 12% are monocytes. PBMCs were isolated, pellets were stored at − 80 °C before mRNA was extracted using RNeasy Mini Kit (Qiagen) as described elsewhere [21]. RNA quantity was measured using NanoDrop-1000 (NanoDrop Technologies, Inc., Delaware, USA), while RNA quality was checked with Aglient 2100 Bioanalyzer (Agilent Technologies, Inc., California, USA). All sample had a RIN-value above 9, except from one sample that was excluded from further analysis. We also excluded samples from four subjects from further analysis due to high-sensitive C-reactive protein (hs-CRP) > 10 Additional file 7: Figure S2. All RNA samples in the present study are the same as used in our previous study by Gjevestad et al. [21]. In addition, serum samples for the determination of hs-CRP, and EDTA-blood for the differential blood count were collected. The analyses were performed at an accredited laboratory (Fürst Laboratories, Oslo, Norway).

Microarray hybridization and processing

Following RNA preparation and amplification, using the Illumina Total Prep RNA Amplification Kit (Illumina Inc., California, USA), gene expression measurements were performed by hybridizing the amplified RNA to Illumina HumanHT-12 v4 Expression BeadChip (Illumina Inc., California, USA) according to the manufacturer’s instructions. Samples were scanned using Illumina HiScan System (Illumina Inc., California, USA). The Illumina HumanHT-12 v4 Expression BeadChip provides genome-wide measurements of the expression of over 48,000 probe sets. IlluminaGenome Studio was used to calculate and report a detection p value, which represents the confidence that a given transcript is expressed above the background. A gene was defined as expressed when relevant probes with a p value below 0.01 were found in more than five samples. After hybridization and scanning, a manual quality control was performed investigating density plots and hierarchical clustering of raw probe densities. One probe per gene (max IQR) was selected for further analysis. The microarray experiments were conducted according to the MIAME (Minimum Information about a Microarray Experiment) guidelines.

Microarray data analysis

After correcting for background noise, using normexp background correction (neqc filtration, Limma), quantile normalization of the data was performed using the Illumina GenomeStudio software, version 1.7.0. Data were log2-transformed and exported raw (non-normalized) to R (http://www.r-project.org/) for biostatistical analysis using the Linear Models for Microarray Data (Limma) Bioconductor package version 1.1.0. Differential gene expression was estimated by a moderated paired t test (Limma) by comparing the relative change from baseline to after the intervention using R software. Gene transcripts that were significantly regulated during the intervention (nominal p value < 0.05) were subjected to further gene pathway analysis using MetaCoreTM (GeneGo, Thomson Reuters, Michigan, USA). Pathways identified in MetaCoreTM with a FDR q-value < 0.05 were considered significantly modulated.

RT-qPCR

RNA was reversely transcribed by a high-capacity cDNA reverse transcription kit (Applied Biosystems). RT-qPCR was performed on an BioRad CFX96 (Bio–Rad Laboratories) with inventoried TaqMan gene expression assays for Killer cell immunoglobulin like receptor, two Ig domains and long cytoplasmic tail 1 (KIR2DL1, Hs 04961778_gH), KIR2DL4 (Hs00427106_m1), Killer cell lectin like receptor C3 (KLRC3, Hs01652462_m1), and C-C motif chemokine ligand 4 like 2 (CCL4L2, Hs04400556_m1) (Thermo Fisher Scientific). TATA-binding protein (TBP, Hs00427620_m1) and glucuronidase beta (GUSB, Hs00939627_m1) were chosen as reference genes due to previous experience with these genes in PBMCs [44]. The assays used for the selected genes were chosen due to best coverage according to Thermo Fischer. The relative mRNA level for each transcript was calculated by the ∆∆cycle threshold (Ct) method [45]. Ct values for each target gene were normalized to the average Ct value of the reference genes (Ctreference − Cttarget = ∆Ct) and the relative change from baseline to the end of study visits was calculated and expressed as log ratio (∆Ctend of study − ∆Ctbaseline = ∆∆Ct).

Measurements of insulin-like growth factor 1

Serum levels of IGF-1 were measured using an ELISA-kit (R&D Systems Inc., Minneapolis, USA) in accordance with the protocol provided. All samples were measured in duplicates.

Measurements of gut peptides

Plasma concentrations of glucose-dependent insulinotropic polypeptide (GIP), glucagon-like peptide-1 (GLP-1), pancreatic polypeptide (PP), amylin and peptide YY (PYY) were analyzed using Milliplex Map Kit for human metabolic hormone magnetic bead panel (Cat. no. HMHEMAG-34 K, EMD Millipore Corporation MA, USA). All samples were measured in duplicate along with controls using Bio-Plex 200 system, based on Luminex xMAP technology (Bio-Rad Laboratories Inc., CA, USA) [46].

Statistical analysis

For normally distributed data, differences between study-groups were performed using the independent samples t test, and Wilcoxon-Mann-Whitney test was performed on not normally distributed data. Changes within each study group were performed using the paired sample t test or Wilcoxon-Mann-Whitney test. We considered a p value of < 0.05 statistically significant when comparing within and between groups. SPSS statistical software, version 25 from Microsoft (SPSS, Inc., USA) and R were used for statistical analyses.

Availability of data and materials

The datasets generated and analyzed during the current study are available from the corresponding author on reasonable request, pending the authorization to deposit them in a public repository.

Abbreviations

ACTH:

Adrenocorticotrophic hormone

ATF4:

Activating transcription factor 4

BMI:

Body mass index

CLIP:

CAP-Gly Domain Containing Linker Protein

DAP12:

TYRO protein tyrosine kinase binding protein

DXA:

Dual energy X-ray absorptiometry

FDR:

False discover rate

GCN2:

General control nonderepressible 2

GIP:

Glucose-dependent insulinotropic polypeptide

GLP-1:

Glucagon-like peptide-1

hs-CRP:

High-sensitive C-reactive protein

IGF-1:

Insulin-like growth factor 1

IL:

Interleukin

JP:

Joining peptide

KDM5A:

lysine-specific demethylase 5A

KIRs:

killer cell immunoglobulin-like receptors

KLRC3:

Killer Cell Lectin Like Receptor C3

LPH:

Lipotrophin

MICA/B:

MHC Class I Polypeptide-Related Sequence A/B

MSH:

Melanocyte-stimulating hormone

mTORC11:

Mechanistic target of rapamycin complex 1

NK-cells:

Natural killer cells

NKG2D:

Natural killer group 2, member D

PBMC:

Peripheral blood mononuclear cells

POMC:

Pro-opiomelanocortin

PP:

Pancreatic polypeptide

PRL23A:

Ribosomal protein L23a

PYY:

Amylin and peptide YY

TNFA:

Tumor necrosis factor alpha

TNFRSF1A:

TNF Receptor Superfamily Member 1A

ZNF543:

Zinc finger protein 543

References

  1. Rennie MJ, Wackerhage H, Spangenburg EE, Booth FW. Control of the size of the human muscle mass. Annu Rev Physiol. 2004;66:799–828.

    Article  CAS  PubMed  Google Scholar 

  2. Roubenoff R. Physical activity, inflammation, and muscle loss. Nutrition reviews. 2007;65(12 Pt 2):S208–12.

    Article  PubMed  Google Scholar 

  3. da Silva AT, de Oliveira Duarte YA, Ferreira Santos JL, Wong R, Lebrao ML. Sarcopenia according to the european working group on sarcopenia in older people (EWGSOP) versus Dynapenia as a risk factor for disability in the elderly. The journal of nutrition, health & aging. 2014;18(5):547–53.

    Article  Google Scholar 

  4. Ambrose AF, Paul G, Hausdorff JM. Risk factors for falls among older adults: a review of the literature. Maturitas. 2013;75(1):51–61.

    Article  PubMed  Google Scholar 

  5. Bauer J, Biolo G, Cederholm T, Cesari M, Cruz-Jentoft AJ, Morley JE, et al. Evidence-based recommendations for optimal dietary protein intake in older people: a position paper from the PROT-AGE Study Group. Journal of the American Medical Directors Association. 2013;14(8):542–59.

    Article  PubMed  Google Scholar 

  6. Baum JI, Kim IY, Wolfe RR. Protein Consumption and the Elderly: What Is the Optimal Level of Intake? Nutrients. 2016;8(6).

    Article  PubMed Central  CAS  Google Scholar 

  7. Houston DK, Nicklas BJ, Ding J, Harris TB, Tylavsky FA, Newman AB, et al. Dietary protein intake is associated with lean mass change in older, community-dwelling adults: the Health, Aging, and Body Composition (Health ABC) Study. The American journal of clinical nutrition. 2008;87(1):150–5.

    Article  CAS  PubMed  Google Scholar 

  8. Kim H, Suzuki T, Kim M, Kojima N, Yoshida Y, Hirano H, et al. Incidence and predictors of sarcopenia onset in community-dwelling elderly Japanese women: 4-year follow-up study. Journal of the American Medical Directors Association. 2015;16(1):85 e1-8.

    Article  Google Scholar 

  9. Ottestad I, Lovstad AT, Gjevestad GO, Hamarsland H, Saltyte Benth J, Andersen LF, et al. Intake of a Protein-Enriched Milk and Effects on Muscle Mass and Strength. A 12-Week Randomized Placebo Controlled Trial among Community-Dwelling Older Adults. The journal of nutrition, health & aging. 2017;21(10):1160–9.

    Article  CAS  Google Scholar 

  10. Mamerow MM, Mettler JA, English KL, Casperson SL, Arentson-Lantz E, Sheffield-Moore M, et al. Dietary protein distribution positively influences 24-h muscle protein synthesis in healthy adults. The Journal of nutrition. 2014;144(6):876–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Tieland M, Beelen J, Laan ACM, Poon S, de Groot L, Seeman E, et al. An Even Distribution of Protein Intake Daily Promotes Protein Adequacy but Does Not Influence Nutritional Status in Institutionalized Elderly. Journal of the American Medical Directors Association. 2018;19(1):33–9.

    Article  PubMed  Google Scholar 

  12. Albers R, Bourdet-Sicard R, Braun D, Calder PC, Herz U, Lambert C, et al. Monitoring immune modulation by nutrition in the general population: identifying and substantiating effects on human health. The British journal of nutrition. 2013;110(Suppl 2):S1–30.

    Article  CAS  PubMed  Google Scholar 

  13. Afman L, Muller M. Nutrigenomics: from molecular nutrition to prevention of disease. Journal of the American Dietetic Association. 2006;106(4):569–76.

    Article  CAS  PubMed  Google Scholar 

  14. Haro D, Marrero PF, Relat J. Nutritional Regulation of Gene Expression: Carbohydrate-, Fat- and Amino Acid-Dependent Modulation of Transcriptional Activity. Int J Mol Sci. 2019;20(6).

    Article  PubMed Central  Google Scholar 

  15. Bruhat A, Fafournoux P. Recent advances on molecular mechanisms involved in amino acid control of gene expression. Current opinion in clinical nutrition and metabolic care. 2001;4(5):439–43.

    Article  CAS  PubMed  Google Scholar 

  16. Vary TC, Lynch CJ. Nutrient signaling components controlling protein synthesis in striated muscle. The Journal of nutrition. 2007;137(8):1835–43.

    Article  CAS  PubMed  Google Scholar 

  17. Afman L, Milenkovic D, Roche HM. Nutritional aspects of metabolic inflammation in relation to health--insights from transcriptomic biomarkers in PBMC of fatty acids and polyphenols. Molecular nutrition & food research. 2014;58(8):1708–20.

    Article  CAS  Google Scholar 

  18. de Mello VD, Kolehmanien M, Schwab U, Pulkkinen L, Uusitupa M. Gene expression of peripheral blood mononuclear cells as a tool in dietary intervention studies: What do we know so far? Molecular nutrition & food research. 2012;56(7):1160–72.

    Article  CAS  Google Scholar 

  19. Visvikis-Siest S, Marteau JB, Samara A, Berrahmoune H, Marie B, Pfister M. Peripheral blood mononuclear cells (PBMCs): a possible model for studying cardiovascular biology systems. Clin Chem Lab Med. 2007;45(9):1154–68.

    Article  CAS  PubMed  Google Scholar 

  20. Pasterkamp G, Daemen M. Circulating cells: the biofactory for markers of atherosclerotic disease. European heart journal. 2008;29(22):2701–2.

    Article  PubMed  Google Scholar 

  21. Gjevestad GO, Ottestad I, Biong AS, Iversen PO, Retterstol K, Raastad T, et al. Consumption of protein-enriched milk has minor effects on inflammation in older adults-A 12-week double-blind randomized controlled trial. Mech Ageing Dev. 2017;162:1–8.

    Article  CAS  PubMed  Google Scholar 

  22. Bicknell AB. The tissue-specific processing of pro-opiomelanocortin. J Neuroendocrinol. 2008;20(6):692–9.

    Article  CAS  PubMed  Google Scholar 

  23. Ehrlich S, Weiss D, Burghardt R, Infante-Duarte C, Brockhaus S, Muschler MA, et al. Promoter specific DNA methylation and gene expression of POMC in acutely underweight and recovered patients with anorexia nervosa. J Psychiatr Res. 2010;44(13):827–33.

    Article  PubMed  Google Scholar 

  24. Wardlaw SL. Hypothalamic proopiomelanocortin processing and the regulation of energy balance. Eur J Pharmacol. 2011;660(1):213–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Valassi E, Scacchi M, Cavagnini F. Neuroendocrine control of food intake. Nutrition, metabolism, and cardiovascular diseases : NMCD. 2008;18(2):158–68.

    Article  CAS  PubMed  Google Scholar 

  26. Davidenko O, Darcel N, Fromentin G, Tome D. Control of protein and energy intake - brain mechanisms. European journal of clinical nutrition. 2013;67(5):455–61.

    Article  CAS  PubMed  Google Scholar 

  27. Journel M, Chaumontet C, Darcel N, Fromentin G, Tome D. Brain responses to high-protein diets. Advances in nutrition. 2012;3(3):322–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Kinzig KP, Hargrave SL, Hyun J, Moran TH. Energy balance and hypothalamic effects of a high-protein/low-carbohydrate diet. Physiology & behavior. 2007;92(3):454–60.

    Article  CAS  Google Scholar 

  29. Dardevet D, Sornet C, Balage M, Grizard J. Stimulation of in vitro rat muscle protein synthesis by leucine decreases with age. The Journal of nutrition. 2000;130(11):2630–5.

    Article  CAS  PubMed  Google Scholar 

  30. Yokoyama WM, Plougastel BF. Immune functions encoded by the natural killer gene complex. Nat Rev Immunol. 2003;3(4):304–16.

    Article  CAS  PubMed  Google Scholar 

  31. Maggio M, De Vita F, Lauretani F, Butto V, Bondi G, Cattabiani C, et al. IGF-1, the cross road of the nutritional, inflammatory and hormonal pathways to frailty. Nutrients. 2013;5(10):4184–205.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Levine ME, Suarez JA, Brandhorst S, Balasubramanian P, Cheng CW, Madia F, et al. Low protein intake is associated with a major reduction in IGF-1, cancer, and overall mortality in the 65 and younger but not older population. Cell metabolism. 2014;19(3):407–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Cadogan J, Eastell R, Jones N, Barker ME. Milk intake and bone mineral acquisition in adolescent girls: randomised, controlled intervention trial. Bmj. 1997;315(7118):1255–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Ben-Shlomo Y, Holly J, McCarthy A, Savage P, Davies D, Davey SG. Prenatal and postnatal milk supplementation and adult insulin-like growth factor I: long-term follow-up of a randomized controlled trial. Cancer Epidemiol Biomarkers Prev. 2005;14(5):1336–9.

    Article  CAS  PubMed  Google Scholar 

  35. Heaney RP, McCarron DA, Dawson-Hughes B, Oparil S, Berga SL, Stern JS, et al. Dietary changes favorably affect bone remodeling in older adults. Journal of the American Dietetic Association. 1999;99(10):1228–33.

    Article  CAS  PubMed  Google Scholar 

  36. Zhu K, Du X, Cowell CT, Greenfield H, Blades B, Dobbins TA, et al. Effects of school milk intervention on cortical bone accretion and indicators relevant to bone metabolism in Chinese girls aged 10-12 y in Beijing. The American journal of clinical nutrition. 2005;81(5):1168–75.

    Article  CAS  PubMed  Google Scholar 

  37. Sharples AP, Hughes DC, Deane CS, Saini A, Selman C, Stewart CE. Longevity and skeletal muscle mass: the role of IGF signalling, the sirtuins, dietary restriction and protein intake. Aging Cell. 2015;14(4):511–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Sandri M, Barberi L, Bijlsma AY, Blaauw B, Dyar KA, Milan G, et al. Signalling pathways regulating muscle mass in ageing skeletal muscle: the role of the IGF1-Akt-mTOR-FoxO pathway. Biogerontology. 2013;14(3):303–23.

    Article  CAS  PubMed  Google Scholar 

  39. Barbieri M, Ferrucci L, Ragno E, Corsi A, Bandinelli S, Bonafe M, et al. Chronic inflammation and the effect of IGF-I on muscle strength and power in older persons. American journal of physiology Endocrinology and metabolism. 2003;284(3):E481–7.

    Article  CAS  PubMed  Google Scholar 

  40. Beavers KM, Hsu FC, Isom S, Kritchevsky SB, Church T, Goodpaster B, et al. Long-term physical activity and inflammatory biomarkers in older adults. Medicine and science in sports and exercise. 2010;42(12):2189–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Singh T, Newman AB. Inflammatory markers in population studies of aging. Ageing research reviews. 2011;10(3):319–29.

    Article  CAS  PubMed  Google Scholar 

  42. Franceschi C. Inflammaging as a major characteristic of old people: can it be prevented or cured? Nutrition reviews. 2007;65(12 Pt 2):S173–6.

    Article  PubMed  Google Scholar 

  43. Schaap LA, Pluijm SM, Deeg DJ, Harris TB, Kritchevsky SB, Newman AB, et al. Higher inflammatory marker levels in older persons: associations with 5-year change in muscle mass and muscle strength. The journals of gerontology Series A, Biological sciences and medical sciences. 2009;64(11):1183–9.

    Article  PubMed  CAS  Google Scholar 

  44. Myhrstad MC, Narverud I, Telle-Hansen VH, Karhu T, Lund DB, Herzig KH, et al. Effect of the fat composition of a single high-fat meal on inflammatory markers in healthy young women. Br J Nutr. 2011;106(12):1826–35.

    Article  CAS  PubMed  Google Scholar 

  45. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

  46. Ingerslev AK, Mutt SJ, Laerke HN, Hedemann MS, Theil PK, Nielsen KL, et al. Postprandial PYY increase by resistant starch supplementation is independent of net portal appearance of short-chain fatty acids in pigs. PloS one. 2017;12(10):e0185927.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

We want to thank all the participants in the trial and Nadia Alazraq for performing RT-qPCR.

Funding

The Research Council of Norway (225258/E.40), Throne Holst Foundation for Nutrition Research, University of Oslo, and TINE SA supported this work.

Author information

Authors and Affiliations

Authors

Contributions

GOG, KBH, IO, and SMU designed the study. GOG, KBH, AR, KK, SJM, KHH, IO, and SMU. conducted the research. AF performed the statistical analysis of the microarray data, while AR, KK and IO performed statistical analysis of other variables. GOG and MM performed the pathway analysis, GOG, KBH, IO, MM and SMU were involved in the interpretation of the data, GOG, KBH and SMU wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Stine M. Ulven.

Ethics declarations

Ethics approval and consent to participante

All subjects provided written informed consent, and the study was conducted according to the Declaration of Helsinki. We received approval for all procedures involving human subjects by the Regional Committees for Medical and Health Research Ethics, Health Region South East, Norway. The study was registered at Clinicaltrials.gov (ID no. NCT02218333).

Consent for publication

Not applicable.

Competing interests

The protein-enriched milk and the isocaloric carbohydrate drink were provided by TINE SA, Oslo, Norway, where G.O.G. is employed. She has no financial interest to declare. K.B.H. and S.M.U have received funding from TINE SA. A.R, A.F, I.O, M.M, K.K, S.J.M, K.H.H report no conflict of interest.

K.B.H. has received research grant from Mills DA, Olympic Seafood, Amgen, Sanofi, Kaneka and Pronova. S.M.U. has received research grant from Mills DA and Olympic Seafood. None of these grants or honoraria is related to the contents of this manuscript.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1: Table S1.

Genes that were differently expressed (p<0.05) after intake of protein. LogFoldChange is calculated by substracting Log2 ratios.

Additional file 2: Table S2.

Genes that were differently expressed (p<0.05) after intake of carbohydrates. LogFoldChange is calculated by substracting Log2 ratios.

Additional file 3: Table S3.

Overlapping genes regulated in the protein and carbohydrate groups-those in bold were regulated in opposite direction.

Additional file 4: Table S4.

All pathway maps regulated after intake of protein (FDR, q-value <0.05). The ratio indicates the number of regulated genes in our gene set compared to the total number of genes included in the pathway. Examples of genes included in each pathway are listed in the column to the right. All pathways are manually classified into selected biological processes (immune response, apoptosis and survival, signal transduction and others)

Additional file 5: Table S5.

All pathways regulated within the carbohydrate group (FDR, q<0.05). The ratio indicates the number of regulated genes in our gene set compared to the total number of genes included in the pathway. Examples of genes included in each pathway are listed in the column to the right. All pathways are manually classified into selected biological processes (immune response, apoptosis and survival, signal transduction and others.

Additional file 6: Figure S1.

Flow chart

Additional file 7: Figure S2.

Pathways altered within the carbohydrate group

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gjevestad, G.O., Holven, K.B., Rundblad, A. et al. Increased protein intake affects pro-opiomelanocortin (POMC) processing, immune function and IGF signaling in peripheral blood mononuclear cells of home-dwelling old subjects using a genome-wide gene expression approach. Genes Nutr 14, 32 (2019). https://0-doi-org.brum.beds.ac.uk/10.1186/s12263-019-0654-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12263-019-0654-6

Keywords