Abstract
Ghrelin is a gastrointestinal peptide hormone that is present in blood mostly in a non-posttranslationally modified form, with a minor proportion acylated at Ser3. Both ghrelin forms were initially assigned a role in the control of food intake but there is accumulating evidence for their involvement in fat allocation and utilization. We investigated changes in the ghrelin system in dairy cows, exhibiting differences in body fat mobilization and fatty liver, from late pregnancy to early lactation. Sixteen dairy cows underwent liver biopsy and were retrospectively grouped based on high (H) or low (L) liver fat content post-partum. Both groups had a comparable feed intake in week −6 (before parturition) and week 2 (after parturition). Only before parturition was preprandial total ghrelin concentration higher in L than in H cows and only after parturition was the basal plasma concentration of non-esterified fatty acids higher in H than in L cows. Both before and after parturition, H cows had higher preprandial plasma concentrations of acyl ghrelin, a higher acyl:total ghrelin ratio, lower plasma triacylglyceride concentrations and a lower respiratory quotient compared with L cows. These group differences could not be attributed to an allelic variant of the acyl ghrelin receptor. Rather, the ratio of acyl:total ghrelin correlated with several aspects of fat metabolism and with respiratory quotient but not with feed intake. These results show that endogenous ghrelin forms are associated with fat allocation, fatty liver, and utilization of fat during the periparturient period.
Introduction
Ghrelin is a 28 amino acid peptide in monogastric species and a 27 amino acid peptide in ruminants (Dickin et al. 2004). It is primarily produced by oxyntic cells in the stomach or abomasum respectively with smaller amounts produced by the intestine and pancreas (Hayashida et al. 2001, Hosoda et al. 2006). A small proportion of ghrelin is acylated (mostly octanoylated) by ghrelin O-acyltransferase (GOAT) at Ser3 (Yang et al. 2008), resulting in the secretion of both acyl and non-acyl ghrelin into the blood stream (Hosoda et al. 2000, ThidarMyint et al. 2006). Both forms are involved in the maintenance of energy homeostasis (ThidarMyint et al. 2006), but only the acyl form of ghrelin exhibits endocrine activity by activating GH secretagogue receptor (GHSR) signalling and GH release (Itoh et al. 2005, Iqbal et al. 2006, Bradford & Allen 2008). In ruminants, plasma concentrations of acyl ghrelin increase before meals, in response to fasting, and with decreasing rumen fill (Sugino et al. 2002a, Wertz-Lutz et al. 2006, Gregorini et al. 2009), whereas feeding suppresses the ghrelin surge (Wertz-Lutz et al. 2006). Thus, acyl ghrelin concentration is influenced by nutritional state in situations of positive energy balance (EB). Wertz-Lutz et al. (2006) found a weak increase in feed intake upon acyl ghrelin injection; however, neither s.c. (Roche et al. 2008) nor centrally infused acyl ghrelin (Iqbal et al. 2006) leads to a significant increase in feed intake, so an orexigenic role of acyl ghrelin in ruminant species appears questionable. Acyl ghrelin is elevated during early lactation (Itoh et al. 2005, Bradford & Allen 2008), a period characterized by high body fat mobilization. This suggests that, in addition to its role in preventing starvation, ghrelin exerts a prominent role in the regulation of fat metabolism. For example, feeding rumen-protected long-chain fatty acids increases the prandial acyl ghrelin surge in cows (Fukumori et al. 2012) and sheep with high levels of adiposity have higher ghrelin levels compared with lean animals after overnight fasting (Kurose et al. 2005), indicating that the ghrelin system is triggered by lipids. On the other hand, infusion of acyl ghrelin increases plasma non-esterified fatty acid (NEFA) concentrations in cows (ThidarMyint et al. 2006, Roche et al. 2008), although these circulating fatty acids do not appear to be released from adipose tissue (Grala et al. 2010). Studies on rodents and humans suggest that the lipolytic effect of acyl ghrelin may occur in muscle tissue (Barazzoni et al. 2005, Vestergaard et al. 2011). In different adipose tissues and in the liver, acyl ghrelin reduces lipid export and increases the expression of lipogenic enzymes, thereby favouring hepatic triacylglyceride (TG) deposition (Barazzoni et al. 2005, Davies et al. 2009, Grala et al. 2010). These studies indicate that acyl ghrelin contributes to the regulation of whole-body fat distribution (Barazzoni et al. 2005), while desacyl ghrelin promotes fatty acid uptake and lipid accumulation in adipose tissue (Miegueu et al. 2011).
Dairy cows adapt to early lactation by mobilizing enormous amounts of adipose tissue (up to 10% of body weight, Komaragiri et al. (1998)), which results in increased plasma NEFA concentrations and liver TG accumulation (Bobe et al. 2004, Vernon 2005, Ingvartsen 2006). Notably, cows with higher plasma NEFA concentrations develop a higher liver fat content (LFC) post-partum (pp) compared with cows exhibiting a lower NEFA level (Hammon et al. 2009). Thus, we hypothesized that the ghrelin system plays a pivotal role in fat allocation and the development of fatty liver in cows entering the early lactation period. Following this idea, we took liver biopsies throughout the periparturient period and grouped the cows with the lowest and highest LFC pp. The results of the study show that animals with a higher LFC have higher preprandial plasma concentrations of acyl ghrelin, higher acyl:total ghrelin ratios, lower plasma TG concentrations and a lower respiratory quotient, both before and after parturition. The ratio of acyl:total ghrelin correlated with several aspects of fat metabolism and with the respiratory quotient but not with feed intake.
Materials and methods
Animal experiment
All experimental procedures were approved by an ethics committee of the State Government in Mecklenburg-West Pommerania (Registration No. LALLF M-V/TSD/7221.3-2.1-021/09). Sixteen of 20 multiparous high-yielding dairy cows (at the end of their second to forth lactation, each >10.000 kg milk/305d in an earlier lactation) were bought from one commercial dairy farm and kept in tie stalls in the institutional facility. To reduce genetically based variation in fat metabolism, all cows were selected to be heterozygous at a polymorphic locus in acyl-CoA: diacylglycerol-acyltransferase 1 gene (DGAT1 K232A: amino acid lysine or alanine at position 232), which is known to affect body fat (Thaller et al. 2003) and milk fat (Kühn et al. 2007) content. Animals were not milked within the 7 weeks before the expected calving date. They received a total mixed ration twice daily (0700 and 1500 h) during weeks 6 to 4 before parturition consisting of grass silage, corn silage, barley straw, hay, concentrate and mineral feed (net energy lactation (NEL)=5.9 MJ/kg of dry matter (DM), utilizable protein=128 g/kg of DM) according to the recommendations of the German Society of Nutrition Physiology (2001). During lactation, animals received a total mixed ration consisting of grass silage, corn silage, barley straw, hay, concentrate, extracted rapeseed meal, molassed sugar beet pulp, mineral feed and propylene glycol (NEL=7.1 MJ/kg of DM, utilizable protein=163 g/kg of DM). Mineral salt blocks and water were available ad libitum. After parturition, animals were milked twice daily (0630 and 1630 h). Liver biopsies were conducted using a custom-made biopsy needle (diameter 6 mm) on day −34 and day +18 relative to parturition (week −5 and +2 respectively) and liver tissue was frozen at −80 °C before analysis of LFC. Based on LFC on day 18 pp, cows were retrospectively grouped (n=8 in each group) according to low (L) or high (H) LFC (for L: LFC<31%; for H: LFC>32% total fat/g DM).
On days −40 and +9 relative to parturition, body weight was determined and back fat thickness was measured weekly using ultrasound (Aloka SSK-500, PPG Hellige GmbH, Freiburg, Germany). On day −39 (week −6) before parturition, and again on day 10 (week 2) after parturition, cows were equipped with an indwelling jugular catheter and kept in respiratory chambers (15 °C, dark–light cycle 0600 to 1900 h; Derno et al. (2009)) to which they had been well adapted before the 3-day trial (habituation at least three times: eating, ruminating and lying down). O2 consumption and the production of CO2 and CH4 were measured at intervals of 6 min and the respiratory quotient was calculated as described previously (Derno et al. 2009). On the first day of the trial, starting at 0630 h (≙0 h), cows were fed ad libitum for 24 h. On the second day of the trial, animals were deprived of feed for 10 h starting at 0630 h. At 1630 h on the second day of the trial, animals were re-fed ad libitum until 0630 h of the third day of the trial (to determine the 1- and 14-h compensatory feed intake). Feed intake was determined every 6 min by feed disappearance measured by a scale connected to an electronic registration device. Cows were milked at 0630 and 1630 h on all 3 days of trial. Milk yield was recorded and milk samples were taken for analysis (see below).
Blood samples were taken every 60 min (via extension tubing, i.d. 1 mm, 4 m) from the outside of the chamber starting at 0920 h (≙3 h) until 1420 h (≙8 h) on day 1 (ad libitum feeding) and from 0630 h (≙0 h) until re-feeding at 1630 h (≙10 h) on day 2 of the trial. On the third day (14-h re-feeding period), one blood sample was taken at 0630 h (≙0 h). Blood samples for acyl and total ghrelin were dispensed into 9 ml EDTA-containing monovettes (Sarstedt, Nürnbrecht, Germany) supplemented with 50 μl protease inhibitor (Trasylol, Bayer) and were immediately placed on ice. Samples for NEFA and TGs were drawn into 2.7 ml fluoride/EDTA-monovettes (Sarstedt). Blood samples were centrifuged at 1565 g for 20 min at 4 °C and plasma stored at −80 °C before analysis.
Milk analyses
Milk samples were analysed for fat, protein and lactose contents (Landeskontrollverband für Leistungs- und Qualitätsprüfung Mecklenburg-Vorpommern E.V.). Energy-corrected milk yield ((0.038×g crude fat+0.024×g crude protein+0.017×g lactose)×kg milk/3.14) was calculated in kg according to Reist et al. (2002). EB was calculated (EB=NEL intake−(energy-corrected milk yield×3.14+0.293×kg body weight0.75)) (Reist et al. 2002) and expressed in MJ/day.
Analyses of liver fat, glycogen and protein content
Liver tissue (100 mg) was dried in a muffle-type furnace at 105 °C for 3 h to determine the DM content. Dried samples underwent elementary carbon and nitrogen analysis (CNS–2000; LECO Instrumente GmbH, Mönchengladbach, Germany). Fat and protein content was calculated according to the methods described previously (Kuhla et al. 2004). Liver glycogen content was determined using a commercially available photometric test based on amyloglucosidase-catalysed release of glucose (no. 10207748035, Boehringer Mannheim).
Analyses of plasma hormones and metabolites
Plasma NEFA and TG concentrations were determined by an automatic analyzer (Abx Pentra 400; Horiba, Montpellier, France, Clinic for Cattle, Stiftung Tierärztliche Hochschule Hannover, Hannover, Germany) using NEFA kit #600-215S (Wako Chemicals distributed by mti diagnostics, Neuss, Germany) and TG kit A11A01640 (ABX-Horiba). Acyl ghrelin was measured as described previously by Relling et al. (2010) using an RIA kit (#GHRA-88HK; Linco Research, St Charles, Billerica, MO, USA). To this end, 200 μl plasma was acidified with 2 μl of 1 mol/l HCl during thawing. Total ghrelin (acyl+desacyl ghrelin) was measured in 400 μl freeze-dried plasma using the RIA method described previously (ThidarMyint et al. 2006). Due to limited amounts of plasma, this was not done on the sample obtained at 0630 h on the second day of the trial.
Comparative sequencing of the GHSR gene
Several single nucleotide polymorphisms (SNPs) and two different alternative transcript variants (GHSR1a and GHSR1b) have been described for the bovine GHSR gene (Colinet et al. 2009, Komatsu et al. 2010). Primers for PCR amplification and re-sequencing of the GHSR gene were designed based on the reference mRNA sequence (NM_001143736.1) and genomic sequences (AM931584.1 and NW_003103817.1). Amplicons covered the complete protein coding region and the 3′ and 5′ UTRs of the GHSR gene (Supplementary Table 1, see section on supplementary data given at the end of this article). Comparative sequencing of the targeted gene regions (1365 bp) was performed on genomic DNA and cDNA was prepared from liver samples. Total RNA was extracted using the NucleoSpin II RNA kit (Macherey & Nagel, Düren, Germany). cDNA was synthesized by RT using Superscript III RNase H-reverse transcriptase (Invitrogen) with gene-specific primers (GHSR_R1, GHS_R4; Supplementary Table 1) according to the manufacturer's instructions. The PCR amplicons were isolated from agarose gels using the NucleoSpin Extract II kit (Macherey & Nagel) and sequenced with PCR primers on a capillary MEGABACE sequencer (GE Healthcare, Freiburg, Germany). To identify variable DNA positions, the sequences were analysed by visual inspection of sequencing profiles and sequence alignment to the respective bovine reference sequences.
Statistical analyses
One cow from each period (ante-partum (ap) and pp respectively) had to be excluded from the analysis due to illness. Descriptive statistics were calculated and tests for normality conducted (using the UNIVARIATE procedure of Base SAS software). The response variables were approximately normal. Data were evaluated by repeated measures ANOVA (using the MIXED procedure of SAS/STAT software: SAS Institute, Inc. 2009. SAS/STAT 9.2 User's Guide, Second Edition, Cary, NC, USA: SAS Institute, Inc.).
The statistical model for the zootechnical data contained the fixed effects LFC, period (ap and pp) and their interaction. Repeated measures on the same animal were taken into account (by the repeated statement in proc mixed using a compound symmetry residual covariance structure). The statistical model for the plasma data contained the fixed effects LFC, time (hours) and their interaction. Repeated measures on the same cow were taken into account in the same way as described earlier using an autoregressive block diagonal structure of the residual covariance matrix. Least square means (LSM) and s.e.m. were calculated and pairwise tested for each effect in each model using the Tukey–Kramer procedure for pairwise multiple comparisons. P<0.05 was considered significant. Pearson correlation coefficients between ghrelin forms and dry matter intake (DMI), and between ghrelin forms and parameters reflecting fat metabolism and EB, were calculated and tested (using the CORR procedure of Base SAS software).
Results
Before parturition, the groups did not differ in body weight or EB but H cows had significantly greater back fat thickness than L cows (Table 1), which further increased until parturition (data not shown). Animals in both groups were in negative EB 10 days after parturition, with a significant difference between H and L cows (P<0.05; Table 1). Although H cows had a significantly higher milk fat content (P=0.02) compared with L cows on both days, both groups had a comparable energy-corrected milk yield (Table 1) and comparable ad libitum and 14 h compensatory feed intake (Table 1) in week −6 before parturition and week 2 after parturition. Also, liver fat, protein and glycogen contents determined in the ap period did not differ between groups (Fig. 1). From pregnancy to early lactation, LFC increased threefold in H cows (P<0.01) and twofold in L cows (P<0.01; Fig. 1A), while hepatic glycogen decreased in both groups: sevenfold in H cows (P<0.01) and twofold in L animals (P=0.04; Fig. 1B). Moreover, liver protein content decreased after parturition only in cows developing a high LFC (P<0.01; Fig. 1C).
Zootechnical data before and after parturition from cows with a high (H) and a low (L) liver fat content pp (LSM±s.e.m.; P, Tukey–Kramer test)
Ante-partum (ap) | Post-partum (pp) | P (ap vs pp) | ||||||
---|---|---|---|---|---|---|---|---|
L | H | P | L | H | P | L | H | |
Energy balance (MJ/day) | 20.0±4.8 | 18.0±4.8 | 0.99 | −63.4±5.1 | −85.7±4.8 | 0.03 | <0.01 | <0.01 |
Body weight (kg) | 709±18.2 | 768±18.2 | 0.15 | 669±18.5 | 713±18.2 | 0.35 | <0.01 | <0.01 |
Ad libitum DMI (kg) | 10.7±0.8 | 10.6±0.8 | 1 | 17.0±0.9 | 15.9±0.8 | 0.79 | <0.01 | <0.01 |
1-h compensatory DMI (kg) | 3.1±0.4 | 2.9±0.3 | 0.97 | 4.4±0.4 | 5.2±0.3 | 0.39 | 0.03 | <0.01 |
14-h compensatory DMI (kg) | 9.2±0.9 | 9.9±0.9 | 0.94 | 14.0±0.9 | 13.7±0.9 | 1 | <0.01 | <0.01 |
Back fat thickness (mm)a | 12.1±1.8 | 21.1±1.8 | 0.01 | 13.0±1.8 | 23.5±1.8 | <0.01 | <0.01 | 0.28 |
Energy-corrected milk yield (kg)b | 46.0±2.1 | 49.5±2 | 0.25 | |||||
Milk fat content (%) | 4.0±0.4 | 5.5±0.35 | 0.02 |
DMI, dry matter intake.
Note: back fat thickness still increased during ap period until parturition (data not shown) above pp values presented herein.
Mean between day 1 (ad libitum feeding) and 2 (10-h food deprivation).
Because plasma metabolite and hormone concentrations may vary periprandially, we took hourly blood samples during the ad libitum period as well as after feed deprivation. During both the ap and the pp periods, plasma NEFA concentrations increased with the duration of feed deprivation and reduced after re-feeding (P<0.01). In the ap period, plasma NEFA concentrations did not differ between groups, neither during ad libitum feeding nor in response to feed deprivation (Fig. 2A). By contrast, in the pp period, plasma NEFA concentrations were significantly increased in both groups when compared with ap ad libitum feeding or the ap feed deprivation respectively (P<0.01). Furthermore, H cows had significantly higher plasma NEFA concentrations than L cows (P<0.01), but these differences were not evident during feed deprivation (Fig. 2B).
Plasma TG concentrations were not different between L and H cows during ad libitum feeding, neither in the ap nor in the pp period (Fig. 2C and D). However, in both groups, plasma TG concentrations were significantly lower after parturition compared with before parturition (P<0.05). Before parturition, feed deprivation provoked no changes in plasma TG concentrations (Fig. 2C), while after parturition, feed deprivation led to increasing plasma TG concentrations in both groups (P<0.01) and to a greater extent in L than in H cows (P<0.01; Fig. 2D).
Because H and L cows differed markedly in their fat metabolism during the pp period, we next examined whether the ghrelin system was involved in the differences in fat allocation and development of fatty liver. Plasma concentrations of acyl ghrelin were comparable between groups during both ap and pp ad libitum feeding (Fig. 3A and B). However, the pp ad libitum acyl ghrelin concentrations tended to be higher in H than in L cows (P=0.07). Acyl ghrelin decreased from the ap to the pp period in H (P=0.03) but not in L animals (P<0.13). During times of feed deprivation, acyl ghrelin concentrations continuously increased in both groups ap and pp (P<0.01). Moreover, feed-deprived H cows secreted significantly more acyl ghrelin into the circulation than L cows (P<0.02), whereas the extent in each group was higher in the pp than in the ap period (P<0.01; Fig. 3A and B).
To examine whether the higher acyl ghrelin concentrations in H cows observed after feed deprivation may be due to greater production of total (acyl and desacyl) ghrelin, we next measured the concentration of total ghrelin in plasma. During ad libitum feeding, total ghrelin concentrations did not differ between groups ap (Fig. 3C and D) but increased after parturition in H (P=0.03) but not in L cows. Also, total ghrelin continuously increased with duration of feed deprivation in both groups, both before and after parturition (P<0.02). Interestingly, L cows responded to feed deprivation in the ap period with a significantly higher total ghrelin release compared with H cows (P=0.03; Fig. 3C), while this effect disappeared during the pp period (Fig. 3D).
During ad libitum feeding, the ratio between acyl and total ghrelin did not differ between groups, nor between the ap and the pp period. However, during feed deprivation, the ratio decreased from ap to pp in H (P<0.01) but not in L cows. Also, H cows had a significantly higher acyl:total ghrelin ratio during feed deprivation than L cows, before (P=0.02) as well as after (P<0.01) parturition (Fig. 4A and B).
To determine whether differences in L and H phenotypes can be attributed to genetic variation in the ghrelin receptor gene, GHSR, we sequenced the complete coding region as well as the 5′ and 3′ UTRs of the GHSR gene. Comparative sequence analysis revealed that the known polymorphisms leading to amino acid variation in the GHSR protein (L24V, D194N and DelR242) displayed no variability in the cows included in our study. Two variable SNPs (c.−7 A>C and c. −128 T>G) located in the 5′ UTR of the GHSR gene were identified compared with the reference sequences (NM_001143736.1, AM931584.1 and NW_003103817.1, data not shown). However, a clear relationship between allelic variants of these SNPs and cow phenotypes could not be observed.
Group differences in the release of both ghrelin forms may be associated with the type of preferred nutrient utilization, which can be assessed by indirect calorimetry. The respiratory quotient was significantly lower in H compared with L cows (P<0.05) both before and after parturition, as well as during ad libitum feeding and feed deprivation (Fig. 5), although there were no differences in ad libitum or compensatory feed intake between groups (Table 1).
In order to examine whether the plasma ghrelin system plays a role in feed intake control or the regulation of fat metabolism, we performed Pearson correlation studies. To this end, we used the pp ghrelin peak during feed deprivation as a model for the preprandial ghrelin surge. We found no correlation between ad libitum nor between 1- or 14-h compensatory DMI and acyl ghrelin, total ghrelin or the ratio between both forms (Fig. 6A and B). Also, LFC did not correlate with acyl or total ghrelin (data not shown), while we observed a significant positive correlation between the ratio of acyl:total ghrelin and LFC (Fig. 6C). In addition to this, the ratio between acyl and total ghrelin correlated with milk fat content (Fig. 6D), but again no correlations were found between acyl or total ghrelin and milk fat content. Moreover, we found a negative correlation between the acyl:total ghrelin ratio and respiratory quotient (Fig. 6E), as well as a positive correlation between the acyl:total ghrelin ratio and back fat thickness respectively (Fig. 6F) but not with plasma NEFA concentrations.
Discussion
Acyl ghrelin and feed intake
The transition from pregnancy to lactation is characterized by an increase in feed intake and the allocation of mobilized body fat reserves to the mammary gland to meet the energy requirements of milk production. It has been suggested that the ghrelin system is involved in both the control of feed intake and the regulation of fat metabolism, which led us to examine whether ghrelin is involved in the adaptation to lactation in dairy cows. Pioneer studies suggested that the desacyl ghrelin form depressed feed intake (Asakawa et al. 2005) and diminished lipolysis (Muccioli et al. 2004, Thompson et al. 2004, Theander-Carrillo et al. 2006), while acyl ghrelin exerted the opposite effect on food intake (Tschöp et al. 2000). However, these results are currently the subject of controversial debate (Wortley et al. 2004, Barazzoni et al. 2007). In ruminants, studies have shown that plasma acyl ghrelin concentrations decrease after feeding and increase before meals and in response to fasting (Sugino et al. 2002a, Wertz-Lutz et al. 2006), using an intervention strategy to model the preprandial acyl ghrelin surge (Bradford & Allen 2008) and synchronize individuals for studying short-term compensatory feed intake. From these observations, it appears that ghrelin might play a role in the control of feed intake in ruminants, as suggested in rodents (Tschöp et al. 2000, Nakazato et al. 2001). In this study, we measured plasma acyl ghrelin concentrations at hourly intervals during ad libitum feeding and feed deprivation. Probably because of different individual feed intakes relative to sampling, we did not observe preprandial acyl ghrelin surges during ad libitum feeding, neither before nor after parturition. This is in agreement with findings in sheep (Sugino et al. 2002b) and our data also support the observation that acyl ghrelin concentrations increase in response to feed deprivation (Wertz-Lutz et al. 2006). The current data are in contrast to the findings of Bradford & Allen (2008), showing that late-lactating cows (positive EB) do not respond with increasing acyl ghrelin concentration after feed withdrawal. However, our cows investigated 6 weeks ap developed only a twofold increase while cows in early lactation responded with a fourfold increase in plasma acyl ghrelin after 9-h feed deprivation. The greater increase in acyl ghrelin in H cows pp compared with ap is in line with earlier findings (Itoh et al. 2005, Bradford & Allen 2008), stating that negative EB increases periprandial acyl ghrelin concentrations (Bradford & Allen 2008). During ad libitum feeding, only H cows (not L cows) had a significantly higher acyl ghrelin concentration pp when compared with the ap period. Only after feed deprivation did both groups exhibit significantly different plasma acyl ghrelin concentrations, before and after parturition, suggesting that it is not negative EB per se (ad libitum-fed early-lactating cows are already in negative EB) but rather reduced feed intake in early lactation that increases preprandial acyl ghrelin concentrations. In addition, if the magnitude of the negative EB per se accounted for the preprandial acyl ghrelin pulse, large differences in acyl ghrelin between groups would be expected during ad libitum feeding because both groups differ significantly in their EB; we did not find such differences, and there was no correlation between preprandial acyl ghrelin peak and negative EB.
There was also no correlation between acyl ghrelin and the amount of DMI during ad libitum feeding or after re-feeding. The latter finding leads to doubt regarding the orexic characteristics of acyl ghrelin in dairy cows, supporting earlier findings that peripheral (Roche et al. 2008) or central (Iqbal et al. 2006) infusion of acyl ghrelin failed to significantly increase feed intake in ruminants.
Desacyl ghrelin, the ratio between both ghrelin forms and feed intake
One reason for the proposed non-orexic character of acyl ghrelin in ruminants might be due to a masking effect of desacyl ghrelin, which is not only the dominant ghrelin form but may also exert inverse effects, at least in rodents (Asakawa et al. 2005). As already observed for the acyl form, total ghrelin concentration did not change during ad libitum feeding but increased from ap to pp only in H cows; during the pp period, ad libitum-fed H cows tended to have higher total ghrelin concentrations than L cows. However, the ratio calculated between ad libitum acyl and total ghrelin concentrations revealed no differences between groups or between ap and pp, which is in line with the comparable DMI of both groups. During deprivation, total ghrelin concentrations increased with time in both groups. Interestingly, rats also respond to fasting with an increase in desacyl ghrelin and no change in acyl ghrelin concentration (Kirchner et al. 2009). In fasted compared with fed humans, increased desacyl but decreased acyl ghrelin concentrations have been described while total ghrelin concentrations remained unchanged (Liu et al. 2008). Thus, as in monogastric species, ruminants respond to fasting with increased release of desacyl ghrelin. The fasting-induced increase in circulating desacyl ghrelin may be in line with its orexic effect described by Toshinai et al. (2006), but the majority of studies performed in rodents report reduced feed intake after central desacyl ghrelin injection (Asakawa et al. 2005, Chen et al. 2005). To the best of our knowledge, there are no data on the effect of desacyl ghrelin on feed intake in ruminants. However, in rats, it has been suggested to be involved in the regulation of energy homeostasis because of its ability to cross the blood–brain barrier and to induce increased neuronal activity in the hypothalamic arcuate nucleus (Inhoff et al. 2009).
In the ap period, total ghrelin was higher in feed-deprived L cows, while at the same time acyl ghrelin was higher in H cows. Accordingly, the acyl:total ghrelin ratio was higher in feed-deprived H than in L cows. The latter finding can be seen in both the ap and the pp periods, suggesting that H cows possess either higher GOAT or lower esterase activity compared with L cows. Moreover, because the ratio did not change in either of the groups with progressive duration of feed deprivation, it may be possible that the granula of the X/A cells of the abomasum (Hayashida et al. 2001) store high enough amounts of both ghrelin forms to release them in equal amounts slowly over time.
Ghrelin and fat metabolism; ghrelin and respiratory quotient
There is increasing evidence from recent literature that acyl ghrelin is more involved in lipid mobilization, fat distribution between tissues and substrate utilization than in the control of feed intake (Wortley et al. 2004, Barazzoni et al. 2005, Kurose et al. 2005, Rigault et al. 2007, Huda et al. 2011, Vestergaard et al. 2011). For example, acyl ghrelin administration into the human femoral artery elevates venous NEFA concentrations, suggesting that ghrelin stimulates lipolysis in skeletal muscle (Vestergaard et al. 2011). Similarly, acyl ghrelin infusion into cows leads to an elevation of plasma NEFA concentration (ThidarMyint et al. 2006, Roche et al. 2008). These findings are consistent with our results showing that early-lactating H cows with the higher preprandial acyl ghrelin surge have a higher lipolysis rate (as indicated by higher basal plasma NEFA concentrations). On the other hand, long-chain fatty acids fed in a rumen-protected form may trigger the acylation of ghrelin (Fukumori et al. 2012). Whether endogenously released fatty acids from adipose tissue contribute to the higher preprandial acyl ghrelin concentration during early lactation remains to be determined.
The different magnitude of increase of both ghrelin forms in response to feed deprivation might also be associated with a different sensitivity of ghrelin receptors. Whereas the receptor for desacyl ghrelin is not yet known, acyl ghrelin binds to GHSR1a (Komatsu et al. 2010, 2012). SNPs in the GHSR gene were found to be associated with growth and carcass traits in cattle (Zhang et al. 2009, Komatsu et al. 2011) and with growth, eating behaviour and obesity in humans (Baessler et al. 2005, Gueorguiev et al. 2009). In the current study, the allele distribution of the identified variable GHSR SNPs did not differ between groups, suggesting that the differences in fat distribution and respiratory quotient between H and L cows cannot be attributed to variations in the GHSR gene. However, differences in expression of the GHSR1a and GHSR1b genes, as is found in several tissues and is altered in an age-dependent manner (Komatsu et al. 2012), may be associated with the H and L phenotype. This question will be addressed in future studies.
Administration of acyl ghrelin in rats also diminishes the TG content in muscle (Barazzoni et al. 2005) and reduces hepatic mitochondrial β-oxidation (Rigault et al. 2007), leading to an increase in liver TGs (Barazzoni et al. 2005). This tissue partitioning effect may be present in periparturient dairy cows because H cows (with higher preprandial acyl ghrelin concentrations) developed higher liver and milk fat contents in early lactation compared with L cows. This assumption is further strengthened by the fact that H cows exhibited a higher ratio between acyl and total ghrelin. Desacyl ghrelin stimulates fatty acid uptake in cultured adipocytes (Miegueu et al. 2011), reduces gastric emptying and thus leads to negative EB in rodents (Asakawa et al. 2005). Interestingly, ad libitum-fed H cows, which tended to have higher total ghrelin concentrations, were indeed in more negative EB. The association between higher desacyl ghrelin concentration and lower feed passage rate in ad libitum-fed H cows was not investigated. However, in obese fasted humans, desacyl ghrelin was lower and acyl ghrelin was elevated compared with non-obese counterparts (Rodríguez et al. 2009); this supports our findings in feed-deprived H and L cows before parturition. The ratio of acyl:total ghrelin is also higher in obese than in non-obese humans (Barazzoni et al. 2007, Rodríguez et al. 2009), suggesting this ratio as a suitable indicator for body fat. Following this idea, we found strong direct correlations between the acyl:total ghrelin ratio and parameters reflecting the extent of body fat mobilization in early lactation, such as liver and milk fat contents and back fat thickness. Additionally, we found a strong indirect correlation between the acyl:total ghrelin ratio and the respiratory quotient, the latter a reciprocal measure of whole-body fat oxidation. Surprisingly, H cows had a lower respiratory quotient not only in early lactation but also in late gestation (when fat mobilization does not occur), indicating a generally greater utilization of fat in H than in L cows, independent of the fatty acid load in early lactation. Consistently, H cows with a higher acyl:total ghrelin ratio and higher back fat thickness had a greater propensity to mobilize and oxidize body fat reserves, which confirms the suggested role of ghrelin in influencing metabolic fuel preference (Wortley et al. 2004).
Taking the results of this study together, we conclude that neither acyl nor desacyl ghrelin concentrations refer to the amount of DMI of high-yielding dairy cows in early lactation. By contrast, due to the high correlation of the acyl:total ghrelin ratio with back fat thickness, liver and milk fat contents, but also whole-body fat oxidation, the ghrelin system links mobilized tissue fat distribution, fat retention in the liver and metabolic fuel preference to ensure proper metabolic adaptation to lactation.
Supplementary data
This is linked to the online version of the paper at http://dx.doi.org/10.1530/JOE-12-0384.
Declaration of interest
The authors declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported.
Funding
This study was supported by Deutsche Forschungsgemeinschaft (DFG; KU 1956/4-1). This research did not receive any grant from the commercial sector.
Acknowledgements
The authors thank U Wiedemuth and S Wöhl for biochemical analyses and the staff at the FBN ‘Tiertechnikum’ for assistance with animal care. They also further acknowledge the help of the Cattle Breeding Organisation Mecklenburg-West Pommerania (RMV) and the Griepentrog farm for assortment of cows. Linguistic refinements of the text by Annelise Myers and Kim Russell are gratefully acknowledged.
References
Asakawa A, Inui A, Fujimiya M, Sakamaki R, Shinfuku N, Ueta Y, Meguid MM & Kasuga M 2005 Stomach regulates energy balance via acylated ghrelin and desacyl ghrelin. Gut 54 18–24. (doi:10.1136/gut.2004.038737)
Baessler A, Hasinoff MJ, Fischer M, Reinhard W, Sonnenberg GE, Olivier M, Erdmann J, Schunkert H, Doering A & Jacob HJ et al. 2005 Genetic linkage and association of the growth hormone secretagogue receptor (ghrelin receptor) gene in human obesity. Diabetes 54 259–267. (doi:10.2337/diabetes.54.1.259)
Barazzoni R, Bosutti A, Stebel M, Cattin MR, Roder E, Visintin L, Cattin L, Biolo G, Zanetti M & Guarnieri G 2005 Ghrelin regulates mitochondrial-lipid metabolism gene expression and tissue fat distribution in liver and skeletal muscle. American Journal of Physiology. Endocrinology and Metabolism 288 E228–E235. (doi:10.1152/ajpendo.00115.2004)
Barazzoni R, Zanetti M, Ferreira C, Vinci P, Pirulli A, Mucci M, Dore F, Fonda M, Ciocchi B & Cattin L et al. 2007 Relationships between desacylated and acylated ghrelin and insulin sensitivity in the metabolic syndrome. Journal of Clinical Endocrinology and Metabolism 92 3935–3940. (doi:10.1210/jc.2006-2527)
Bobe G, Young JW & Beitz DC 2004 Pathology, ethiology, prevention and treatment of fatty liver in dairy cows. Journal of Dairy Science 87 3105–3124. (doi:10.3168/jds.S0022-0302(04)73446-3)
Bradford B & Allen MS 2008 Negative energy balance increases periprandial ghrelin and growth hormone concentrations in lactating dairy cows. Domestic Animal Endocrinology 34 196–203. (doi:10.1016/j.domaniend.2007.02.002)
Chen CY, Chao Y, Chang FY, Chien EJ, Lee SD & Doong ML 2005 Intracisternal des-acyl ghrelin inhibits food intake and non-nutrient gastric emptying in conscious rats. International Journal of Molecular Medicine 16 695–699.
Colinet FG, Vanderick S, Charloteaux B, Eggen A, Gengler N, Renaville B, Brasseur R, Portetelle D & Renaville R 2009 Genomic location of the bovine growth hormone secretagogue receptor (GHSR) gene and investigation of genetic polymorphism. Animal Biotechnology 20 28–333. (doi:10.1080/10495390802602926)
Davies JS, Kotokorpi P, Eccles SR, Barnes SK, Tokarczuk PF, Allen SK, Whitworth HS, Guschina IA, Evans BA & Mode A et al. 2009 Ghrelin induces abdominal obesity via GHS-R-dependent lipid retention. Molecular Endocrinology 23 914–924. (doi:10.1210/me.2008-0432)
Derno M, Elsner EA, Scholze EA & Schweigel M 2009 Technical note: a new facility for continuous respiration measurements in lactating cows. Journal of Dairy Science 92 2804–2808. (doi:10.3168/jds.2008-1839)
Dickin JC, Thue TD & Buchanan FC 2004 An alternative splice site in ghrelin is missing in ruminants. Animal Genetics 35 411–412. (doi:10.1111/j.1365-2052.2004.01175.x)
Fukumori R, Sugino T, Shingu H, Moriya N, Hasegawa Y, Kofima M, Kangawa K, Obitsu T, Kushibiki S & Taniguchi K 2012 Effects of calcium salts of long-chain fatty acids and rumen-protected methionine on plasma concentrations of ghrelin, glucagon-like peptide-1 (7 to 36) amide and pancreatic hormones in lactating cows. Domestic Animal Endocrinology 42 74–82. (doi:10.1016/j.domaniend.2011.09.005)
German Society of Nutrition Physiology (Gesellschaft für Ernährungsphysiologie) 2001 Empfehlungen zur Energie- und Nährstoffversorgung der Milchkühe und Aufzuchtrinder. Frankfurt am Main, Germany: DLG-Verlag.
Grala TM, Kay JK, Walker CG, Sheahan AJ, Littlejohn MD, Lucy MC & Roche JR 2010 Expression analysis of key somatotropic axis and liporegulatory genes in ghrelin- and obestatin-infused dairy cows. Domestic Animal Endocrinology 39 76–83. (doi:10.1016/j.domaniend.2010.02.004)
Gregorini P, Soder KJ & Kensinger RS 2009 Effects of rumen fill on short-term ingestive behavior and circulating concentrations of ghrelin, insulin, and glucose of dairy cows foraging vegetative micro-swards. Journal of Dairy Science 93 2095–2105. (doi:10.3168/jds.2008-1803)
Gueorguiev M, Lecoeur C, Meyre D, Benzinou M, Mein CA, Hinney A, Vatin V, Weill J, Heude B & Hebebrand J et al. 2009 Association studies on ghrelin and ghrelin receptor gene polymorphisms with obesity. Obesity 17 745–754. (doi:10.1038/oby.2008.589)
Hammon HM, Stürmer G, Schneider F, Tuchscherer A, Blum H, Engelhard T, Genzel A, Staufenbiel R & Kanitz W 2009 Performance and metabolic and endocrine changes with emphasis on glucose metabolism in high-yielding dairy cows with high and low fat content in liver after calving. Journal of Dairy Science 92 1554–1566. (doi:10.3168/jds.2008-1634)
Hayashida T, Murakami K, Mogi K, Nishihara M, Nakazato M, Mondal MS, Horii Y, Kojima M, Kangawa K & Murakami N 2001 Ghrelin in domestic animals: distribution in stomach and its possible role. Domestic Animal Endocrinology 21 17–24. (doi:10.1016/S0739-7240(01)00104-7)
Hosoda H, Kojima M, Matsuo H & Kangawa K 2000 Ghrelin and des-acyl ghrelin: two major forms of rat ghrelin peptide in gastrointestinal tissue. Biochemical and Biophysical Research Communications 279 909–913. (doi:10.1006/bbrc.2000.4039)
Hosoda H, Kojima M & Kangawa K 2006 Biological, physiological and pharmacological aspects of ghrelin. Journal of Pharmacological Sciences 100 398–410. (doi:10.1254/jphs.CRJ06002X)
Huda MS, Dovey TM, Wong SP, English PJ, Halford JC, McCulloch P, Cleator J, Martin B, Cashen J & Hayden K et al. 2011 Ghrelin does not orchestrate the metabolic changes seen in fasting but has significant effects on lipid mobilisation and substrate utilization. European Journal of Endocrinology 165 45–55. (doi:10.1530/EJE-10-1122)
Ingvartsen KL 2006 Feeding and management-related diseases in the transition cow. Physiological adaptations around calving and strategies to reduce feeding-related diseases. Animal Feed Science and Technology 126 175–213. (doi:10.1016/j.anifeedsci.2005.08.003)
Inhoff T, Wiedenmann B, Klapp BF, Mönnikes H & Kobelt P 2009 Is desacyl ghrelin a modulator of food intake? Peptides 30 991–994. (doi:10.1016/j.peptides.2009.01.019)
Iqbal J, Kurose Y, Canny B & Clarke IJ 2006 Effects of central infusion of ghrelin on food intake and plasma levels of growth hormone, luteinizing hormone, prolactin, and cortisol secretion in sheep. Endocrinology 147 510–519. (doi:10.1210/en.2005-1048)
Itoh F, Komatsu T, Yonai M, Sugino T, Kojima M, Kangawa K, Hasegawa Y, Terashima Y & Hodate K 2005 GH secretory responses to ghrelin and GHRH in growing and lactating dairy cattle. Domestic Animal Endocrinology 28 34–35. (doi:10.1016/j.domaniend.2004.06.001)
Kirchner H, Gutierrez JA, Solenberg PJ, Pfluger PT, Czyzyk TA, Willency JA, Schurmann A, Joost HG, Jandacek RJ & Hale JE et al. 2009 GOAT links dietary lipids with the endocrine control of energy balance. Nature Medicine 15 741–745. (doi:10.1038/nm.1997)
Komaragiri MV, Casper DP & Erdman RA 1998 Factors affecting body tissue mobilization in early lactation dairy cows. 2. Effect of dietary fat on mobilization of body fat and protein. Journal of Dairy Science 81 169–175. (doi:10.3168/jds.S0022-0302(98)75564-X)
Komatsu M, Fujimori J, Sato Y, Okamura H, Sasaki S, Itoh T, Morita M, Nakamura R, Oe T & Furuta M et al. 2010 Nucleotide polymorphisms and 5′UTR transcriptional analysis of the bovine growth hormone secretagogue receptor 1a (GHSR1a) gene. Animal Science Journal 81 530–550. (doi:10.1111/j.1740-0929.2010.00797.x)
Komatsu M, Itoh T, Fujimori Y, Satoh M, Miyazaki Y, Takahashi H, Shimizu K, Malau-Aduli AE & Morita M 2011 Genetic association between GHSR1a 5′UTR-microsatellite and nt-7(C>A) loci and growth and carcass traits in Japanese Black cattle. Animal Science Journal 82 396–405. (doi:10.1111/j.1740-0929.2010.00852.x)
Komatsu M, Kojima M, Okamura H, Nishio M, Kaneda M, Kojima T, Takeda H, Malau-Aduli AE & Takahashi H 2012 Age-related changes in gene expression of the growth hormone secretagogue and growth hormone-releasing hormone receptors in Holstein-Friesian cattle. Domestic Animal Endocrinology 42 83–93. (doi:10.1016/j.domaniend.2011.09.006)
Kuhla S, Klein M, Renne U, Jentsch W, Rudolph PE & Souffrant WB 2004 Carbon and nitrogen content based estimation of the fat content of animal carcasses in various species. Archives of Animal Nutrition 58 37–46. (doi:10.1080/00039420310001656712)
Kühn C, Edel C, Thaller G & Weikard R 2007 Dominance and parent-of-origin effects of coding and non-coding alleles at the acylCoA-diacylglycerol-acyltransferase (DGAT1): gene on milk production traits in German Holstein cows. BMC Genetics 8 62. (doi:10.1186/1471-2156-8-62)
Kurose Y, Iqbal J, Rao A, Murata Y, Hasegawa Y, Terashima Y, Kojima M, Kangawa K & Clarke IJ 2005 Changes in expression of the genes for the leptin receptor and the growth hormone-releasing peptide/ghrelin receptor in the hypothalamic arcuate nucleus with long-term manipulation of adiposity by dietary means. Journal of Neuroendocrinology 17 331–340. (doi:10.1111/j.1365-2826.2005.01318.x)
Liu J, Prudom CE, Nass R, Pezzoli SS, Oliveri MC, Johnson ML, Veldhuis P, Gordon DA, Howard AD & Witcher DR et al. 2008 Novel ghrelin assays provide evidence for independent regulation of ghrelin acylation and secretion in healthy young men. Journal of Clinical Endocrinology and Metabolism 93 1980–1987. (doi:10.1210/jc.2007-2235)
Miegueu P, St Pierre D, Broglio F & Cianflone K 2011 Effect of desacyl ghrelin, obestatin and related peptides on triglyceride storage, metabolism and GHSR signaling in 3T3-L1 adipocytes. Journal of Cellular Biochemistry 112 704–714. (doi:10.1002/jcb.22983)
Muccioli G, Pons N, Ghe C, Catapano F, Granata R & Ghigo E 2004 Ghrelin and des-acyl ghrelin both inhibit isoproterenol-induced lipolysis in rat adipocytes via a non-type 1a growth hormone secretagogue receptor. European Journal of Pharmacology 798 27–35. (doi:10.1016/j.ejphar.2004.07.066)
Nakazato M, Murakami N, Date Y, Kojima M, Matsuo H, Kangawa K & Matsukura S 2001 A role for ghrelin in the central regulation of feeding. Nature 409 194–198. (doi:10.1038/35051587)
Reist M, Erdin D, vonEuw D, Tschuemperlin K, Leuenberger H, Chillard Y, Hammon HM, Morel C, Philipona C & Zbinden Y et al. 2002 Estimation of energy balance individual and herd level using blood and milk traits in high-yielding dairy cows. Journal of Dairy Science 85 3314–3327. (doi:10.3168/jds.S0022-0302(02)74420-2)
Relling AE, Pate JC, Reynolds CK & Loerch SC 2010 Effect of feed restriction and supplemental dietary fat on gut peptide and hypothalamic neuropeptide messenger ribonucleic acid concentrations in growing wethers. Journal of Animal Science 88 737–748. (doi:10.2527/jas.2009-2316)
Rigault C, Le Borgne F, Georges B & Demarquoy J 2007 Ghrelin reduces hepatic mitochondrial fatty acid beta oxidation. Journal of Endocrinological Investigation 30 RC4–RC8. (doi:10.1016/j.mce.2012.11.007)
Roche JR, Sheahan AJ, Chagas LM, Blache D, Berry DP & Kay JK 2008 Long-term infusions of ghrelin and obestatin in early lactation dairy cows. Journal of Dairy Science 91 4728–4740. (doi:10.3168/jds.2008-1193)
Rodríguez A, Gómez-Ambrosi J, Catalán V, Gil MJ, Becerril S, Sáinz N, Silva C, Salvador J, Colina I & Frühbeck G 2009 Acylated and desacyl ghrelin stimulate lipid accumulation in human visceral adipocytes. International Journal of Obesity 33 541–552. (doi:10.1038/ijo.2009.40)
Sugino T, Hasegawa Y, Kikkawa Y, Yamaura J, Yamagishi M, Kurose Y, Kojima M, Kangawa K & Terashima Y 2002a A transient ghrelin surge occurs just before feeding in a scheduled meal-fed sheep. Biochemical and Biophysical Research Communications 295 255–260. (doi:10.1016/S0006-291X(02)00654-X)
Sugino T, Yamaura J, Yamagishi M, Ogura A, Hayashi R, Kurose Y, Kojima M, Kangawa K, Hasegawa Y & Terashima Y 2002b A transient ghrelin secretion before feeding is modified by different feeding regimens in sheep. Biochemical and Biophysical Research Communications 298 785–788. (doi:10.1016/S0006-291X(02)02572-X)
Thaller G, Kühn C, Winter A, Ewald G, Bellmann O, Wegner J, Zühlke H & Fries R 2003 DGAT1, a new positional and functional candidate gene for intramuscular fat deposition in cattle. Animal Genetics 34 354–357. (doi:10.1046/j.1365-2052.2003.01011.x)
Theander-Carrillo C, Wiedmer P, Cettour-Rose P, Nogueiras R, Perez-Tilve D, Pfluger P, Castaneda TR, Muzzin P, Schurmann A & Szanto I et al. 2006 Ghrelin action in the brain controls adipocyte metabolism. Journal of Clinical Investigation 116 1983–1993. (doi:10.1172/JCI25811)
ThidarMyint H, Yoshida H, Ito T & Kuwayama H 2006 Dose-dependent response of plasma ghrelin and growth hormone concentrations to bovine ghrelin in Holstein heifers. Journal of Endocrinology 189 655–664. (doi:10.1677/joe.1.06746)
Thompson NM, Gill DA, Davies R, Loveridge N, Houston PA, Robinson IC & Wells T 2004 Ghrelin and des-octanoyl ghrelin promote adipogenesis directly in vivo by a mechanism independent of the type 1a growth hormone secretagogue receptor. Endocrinology 145 234–242. (doi:10.1210/en.2003-0899)
Toshinai K, Yamaguchi H, Sun Y, Smith RG, Yamanaka A, Sakurai T, Date Y, Mondal MS, Shimbara T & Kawagoe T et al. 2006 Des-acyl ghrelin induces food intake by a mechanism independent of the growth hormone secretagogue receptor. Endocrinology 147 2306–2314. (doi:10.1210/en.2005-1357)
Tschöp M, Smiley DL & Heiman ML 2000 Ghrelin induces adiposity in rodents. Nature 407 908–913. (doi:10.1038/35038090)
Vernon RG 2005 Lipid metabolism during lactation: a review of adipose tissue-liver interactions and the development of fatty liver. Journal of Dairy Science 72 460–469. (doi:10.1017/S0022029905001299)
Vestergaard ET, Buhl M, Gjedsted J, Madsen M, Jessen N, Nielsen S, Gaylinn BD, Liu J, Thorner MO & Moller N et al. 2011 Acute peripheral metabolic effects of intraarterial ghrelin infusion in healthy young men. Journal of Clinical Endocrinology and Metabolism 96 468–477. (doi:10.1210/jc.2010-1995)
Wertz-Lutz AE, Knight TJ, Pritchard RH, Daniel JA, Clapper JA, Smart AJ, Trenkle A & Beitz DC 2006 Circulating ghrelin concentration fluctuate relative to nutritional status and influences feeding behavior. Journal of Animal Science 84 3285–3300. (doi:10.2527/jas.2006-053)
Wortley KE, Anderson KD, Garcia K, Murray JD, Malinova L, Liu R, Moncrieffe M, Thabet K, Cox HJ & Yancopoulos GD et al. 2004 Genetic deletion of ghrelin does not decrease food intae but influences metabolic fuel preference. PNAS 21 8227–8232. (doi:10.1073/pnas.0402763101)
Yang J, Brown MS, Liang G, Grishin NV & Golstein JL 2008 Identification of the acyltransferase that octanoylates ghrelin, an appetite-stimulating peptide hormone. Cell 132 387–396. (doi:10.1016/j.cell.2008.01.017)
Zhang B, Chen H, Guo Y, Zhang L, Zhao M, Lan X, Zhang C, Pan C, Hu S & Wang J et al. 2009 Associations of polymorphism within the GHSR gene with growth traits in Nanyang cattle. Molecular Biology Reports 36 2259–2263. (doi:10.1007/s11033-008-9442-x)