Skip to main content

Integrative analysis of indirect calorimetry and metabolomics profiling reveals alterations in energy metabolism between fed and fasted pigs

Abstract

Background

Fasting is a simple metabolic strategy that is used to estimate the maintenance energy requirement where the energy supply for basic physiological functions is provided by the mobilization of body reserves. However, the underlying metabolic components of maintenance energy expenditure are not clear. This study investigated the differences in heat production (HP), respiratory quotient (RQ) and plasma metabolites in pigs in the fed and fasted state, using the techniques of indirect calorimetry and metabolomics.

Methods

Nine barrows (45.2 ± 1.7 kg BW) were fed corn-soybean based meal diets and were kept in metabolism crates for a period of 14 d. After 7 d adaptation, pigs were transferred to respiratory chambers to determine HP and RQ based on indirect calorimetry. Pigs were fed the diet at 2,400 kJ ME/(kg BW0.6·d) during d 8 to 12. The last 2 d were divided into 24 h fasting and 48 h fasting treatment, respectively. Plasma samples of each pig were collected from the anterior vena cava during the last 3 d (1 d while pigs were fed and 2 d during which they were fasted). The metabolites of plasma were determined by high-resolution mass spectrometry using a metabolomics approach.

Results

Indirect calorimetry analysis revealed that HP and RQ were no significant difference between 24 h fasting and 48 h fasting, which were lower than those of fed state (P <  0.01). The nitrogen concentration of urine tended to decrease with fasting (P = 0.054). Metabolomics analysis between the fed and fasted state revealed differences in 15 compounds, most of which were not significantly different between 24 h fasting and 48 h fasting. Identified compounds were enriched in metabolic pathways related to linoleic acid metabolism, amino acid metabolism, sphingolipid metabolism, and pantothenate and CoA biosynthesis.

Conclusion

These results suggest that the decreases in HP and RQ of growing pigs under fasting conditions were associated with the alterations of linoleic acid metabolism and amino acid metabolism. The integrative analysis also revealed that growing pigs under a 24-h fasting were more appropriate than a 48-h fasting to investigate the metabolic components of maintenance energy expenditure.

Background

Maintenance corresponds to the basal energy requirements for supporting body function, body temperature, and normal physical activity at a time when there is no net gain or loss of energy [1]. Though it is vital to the organism, little is known of the metabolic components of maintenance energy expenditure [2]. For nearly five decades, the fasting method has been used to estimate maintenance. During fasting, pigs mobilize body reserves resulting in a negative energy balance and, during a short duration of fasting, the fasting heat production may be related to the maintenance energy expenditure in the fed state [3, 4]. Therefore, the fasting metabolism may, to some extent reflect, the metabolic components of maintenance in pigs [5]. A comparison between the fed and the fasted state will contribute to identify specific metabolic indicator of maintenance in pigs [6]. However, the duration of fasting has not been consistent in different studies [7, 8]. A series of studies have demonstrated that fasting duration plays a role in blood biochemistry and transcriptional regulation of metabolic genes. During the short-term fasting, plasma non-esterified fatty acids (NEFA) increased, while leptin and insulin concentrations were reduced [9, 10]. In addition, pyruvate dehydrogenase in rat muscle and fatty acid synthase in rat adipose tissue decreased after 1 d of fasting [11]. During the long-term fasting, plasma glucose level was reduced by 48 h of fasting in fasted pigs and dairy cattle [10]. It was noted that the transcriptions of genes involved in the fatty acid oxidation began to decrease in 3-d fasted rat liver [11]. Thus, it is important to investigate the effect of fasting duration on plasma metabolic profiling to understand deeply the fasting metabolism.

Indirect calorimetry is a noninvasive method to study energy expenditure by respiratory gas exchange analysis [12]. However, changes in heat production (HP) and in the respiratory quotient (RQ) after fasting based on indirect calorimetry are a phenotypic reaction of an organism’s metabolism, as it adapts to a different physiological and nutritional condition [13]. Fasting metabolism is associated with many metabolic changes that occur in the body when the animal has to rely on its body reserves to sustain maintenance [14, 15]. Typically, the NEFA level in plasma increases 12 h after fasting in growing pigs [9]. Also, important increases in the ratio between NEFA and glycerol can be observed after prolonged fasting from 12 h to 36 h [6].

Plasma can be used as an accessible metabolic footprint that provides a picture of the metabolic events in the organisms and may reveal changes in metabolic pathways under various physiological or nutritional conditions [16, 17]. The metabolome is defined as the collection and global analysis of all small molecular metabolites generated in a cell, organ or organism [18]. Rubio-Aliaga et al. [6] used metabolomics of prolonged fasting in humans to reveal new catabolic markers. Thus, metabolomics is an ideal tool to explore the effect of fasting on plasma metabolites that result from synthetic and catabolic processes in growing pigs [19, 20].

However, what the association do the plasma metabolites of pigs have with alterations of components of energy expenditure or substrate oxidation pattern under fasting duration, to our knowledge, is not known clearly. We hypothesized that the changes in HP and RQ of fasted pigs were modulated by potential metabolic pathways related to energy metabolism. Therefore, the objective of present study was to investigate the effect of fasting treatment in growing pigs on the alterations in energy metabolism based on the integrative analysis of indirect calorimetry and metabolomics profiling.

Methods

Animal, diets and experimental procedures

The experimental protocol used in the present study was approved by the Institutional Animal Care and Use Committee at China Agricultural University. Nine growing barrows (Duroc × Landrace × Yorkshire), with an average initial BW of 45.2 ± 1.7 kg, were selected from the Fengning Swine Research Unit of China Agricultural University (Hebei, China). The experiment was conducted at the Laboratory of Animal Metabolism of China Agricultural University (Beijing, China). The basal diet (Table 1) was formulated based on corn and soybean meal.

Table 1 Composition and nutrient analysis of the basal diet, as-fed basis

The experiment was conducted in 3 consecutive periods (3 pigs per period) using the same facilities and similar experimental procedures because only 3 respiration chambers were available for the study. In each period, pigs were housed individually in metabolism cages and adapted to the cages and diet for 7 d prior to the start of the experimental period. Pigs were fed 2,400 kJ ME/(kg BW0.6·d) daily during the adaptation period, which was close to the ad libitum feed intake. The experimental period was 7 d, which was divided into three treatments. The first 5 d were feeding treatment. Pigs were also fed 2,400 kJ ME/(kg BW0.6·d) daily. The last 2 d were divided into 24 h fasting and 48 h fasting treatment, respectively. Pigs were fasted with ad libitum access to water via a low-pressure nipple drinker. During the adaptation and experimental period, all pigs were fed the same basal diet and had free access to water. During the experimental period, pigs were housed in metabolism cages that were placed inside open circuit respiratory chambers to determine HP and RQ. The HP of fed state was calculated by averaged the daily HP of the 5-d feeing treatment. The 8-h HP from 22:00 h (d 13) to 06:00 h (d 14) and from 22:00 h (d 14) to 06:00 h (d 15) were calculated and then were extrapolated to a 24-h period, which were considered as the HP of 24 h fasting treatment and 48 h fasting treatments, respectively. Pigs were fed twice daily an equal amount of meal, at 09:00 and 15:30 h. Feed refusals and spillage were recorded daily.

Pigs were weighed at the beginning of adaptation (d 0), and on d 8, 13, 14 and 15, which corresponded to the start or end of the different feeding regimens or fasting times. Temperature in the chambers was maintained at 22 ± 1 °C during the feeding and fasting periods according to Zhang et al. [21].

Blood collection and animal sampling

The feeding, 24 h fasting and 48 h fasting treatments were terminated around 06:00 h on d 13, 14 and 15, respectively. Three blood samples were collected from the anterior vena cava of each pig at these three time points and represented the blood samples of feeding, 24 h fasting and 48 h fasting treatments. Samples were centrifuged (Heraeus, Hanau, Germany) at 3,000×g for 10 min at 4 °C, then the supernatants were transferred to storage tubes, frozen in liquid nitrogen, and stored at − 80 °C for subsequent metabolomics testing.

Urine was collected during the time that pigs were in the respiration chambers according to the methods described by Liu et al. [4]. In brief, urine was collected each morning for each pig in plastic buckets containing 50 mL of 6 N HCl and sieved with cotton gauze, and filtered into a plastic bottle every day. The total quantity of collected urine was weighed and 5% of the daily urinary excretion was stored at − 20 °C. At the end of the experiment, urine collected during the first 5 d was thawed and mixed separately by individual animal and a sub-sample was saved for chemical analysis. Urine was collected separately during the 24 h fasting and 48 h fasting periods. The feed was ground through a 1-mm screen and mixed thoroughly for chemical analysis.

Chemical analysis and calculation

All chemical analyses of ingredients and diet were conducted in duplicate. Dry matter (DM) was measured by drying the samples in a 135 °C drying oven for 2 h (method 930.15) [22]. The total crude protein (N × 6.25) content of the samples was determined using the Kjeldahl method (method 984.13) [22]. Calcium was measured by titration with 0.1 N KMnO4 (method 927.02) [22]. Total phosphorus was measured colorimetrically using a molybodovanadate reagent (method 965.17) [22]. Ether extract (EE) was determined using the Thiex method [23]. Neutral detergent fibre (NDF) and acid detergent fibre (ADF) were determined using filter bags and fibre analyser equipment (Fibre Analyzer; Ankom Technology, Macedon, NY, USA) following a modification of the procedure of van Soest et al. [24]. The gross energy (GE) in diets was determined using a Parr 6400 bomb calorimeter (Parr Instruments, Moline, IL). The nitrogen in urine were measured according to Liu et al. [4].

Heat production and RQ were calculated daily from CO2 and CH4 production, O2 consumption, and nitrogen excretion in the urine (UN) during the 5 d feeding period and the 2 d fasting period according to the following formulas by Brouwer et al. [25]:

$$ \mathrm{HP}\;\left(\mathrm{kJ}\right)=16.1753\times {\mathrm{O}}_2\left(\mathrm{L}\right)+5.0208\times {\mathrm{CO}}_2\left(\mathrm{L}\right)-2.1673\times {\mathrm{CH}}_4\left(\mathrm{L}\right)-5.9873\times \mathrm{UN}\left(\mathrm{g}\right) $$
$$ \mathrm{RQ}={\mathrm{CO}}_2\left(\mathrm{L}\right)/{\mathrm{O}}_2\left(\mathrm{L}\right) $$

Sample preparation for metabolomics

Six plasma samples selected randomly from each treatment were extracted using 800 μL ice-cold extraction mix (acetonitrile:methanol, 1:1, v:v) at a 1:4 sample:extract solution ratio. After vortexing for 5 min, the samples were centrifuged (Eppendorf, Germany) at 15,000×g for 10 min at 4 °C for deproteinization. The supernatant fractions were then collected and evaporated to dryness using a vacuum concentrator (Concentrator plus, Eppendorf). The resulting dry residues were re-suspended in 200 μL of methanol:water (4:1), vortexed and centrifuged again at 15,000×g for 10 min at 4 °C. Lastly, the supernatant fractions were filtered through a 0.1-μm membrane and transferred to sampler vials to be analyzed on a UPLC-MS system.

UPLC-MS analysis

Plasma samples were analyzed with an UPLC-HRMS system (UPLC, ACQUITYUPLC H-Class Bio, Waters; MS, Q-Exactive, Termo Scientifc) equipped with a heated electrospray ionization (HESI) source. The UPLC separation was operated on a BEH C18 column (2.1 mm× 100 mm, 1.7 μm, Waters). The HPLC grade solvents and additives from Fisher Scientifc (ThermoFisher Scientific, NJ, USA) were used. Mobile phase was comprised of 0.1% formic acid water solution (A) and 0.1% formic acid ACN solution (B). The gradient program was as follows: 95% A at 0 min to 70% A at 5 min, 5% A at 10 min and held for 3 min, then returned to initial condition. The flow rate was 0.3 mL/min. A sample of pooled plasma was re-injected after each six samples for quality control. The column temperature was set at 35 °C and the injection volume was 5 μL.

The MS analysis was performed in an electrospray ionization positive mode. Key parameters of the HESI source were as follows: spray voltage at 3 kV, capillary temperature at 320 °C, sheath gas flow rate at 30 arb.units, aux gas flow rate at 10 arb. Units, sweep gas flow rate at 5 arb. Units, heater temperature at 350 °C, and s-lens RF level at 55. Full scan data was acquired with a resolution of 70,000 in the mass range of m/z 67.7–1,000. For MS/MS analysis, normalized collision was performed at an energy of 35 V. In addition, an isolation window of 0.8 m/z and a mass resolution of 35,000 were selected.

Data mining and processing

SIEVE 2.1 software (Thermo Scientific, NJ, USA) was used for metabolomics data processing. This software achieved background subtraction, peak alignment and component extraction of the raw data. Component extraction was performed according to the rules of retention time from 0.5 to 16 min, intensity threshold at 500,000, minimum scan at 9 and signal to noise ratio at 10. Principal components analysis (PCA) was carried out using SIMCA-P 13 software (Umetrics, Umea, Sweden) after data were scaled to Pareto variance. The compounds with P <  0.05, fold change > 1.5 and CV < 30% were picked by EXCEL for further identification.

Identification of differential compounds was performed in compound database of METLIN (https://metlin.scripps.edu/landing_page.php?pgcontent=mainPage) and Human Metabolome Database (http://www.hmdb.ca) using accurate mass of molecular ions. The MS/MS spectra database was used to match fragment ion spectra of the candidate compounds. The MS/MS spectra were also compared with theoretical fragmentation patterns with mass accuracy at 5 ppm using Xcalibur™ (Thermo Scientific, NJ, USA).

Statistical analysis

Data generated in the present experiment were analyzed using the MIXED procedure of SAS (SAS Inst. Inc., Cary, NC) and repeated measurements were considered when analyzed the effects of fasting duration. The statistical model included the main effects of fasting duration. Pig was treated as the experimental unit. Treatment means were separated using the LSMEANS statement and Tukey’s test was used for adjustment in multiple comparison. Results were considered significant at P <  0.05 and considered as trends at 0.05 < P <  0.10. Boxplot analysis of identified differential compounds were achieved using the R software package (R Development Core Team, 2017,) version 3.4.1). The relative concentrations of differential compounds were imported into Metaboanalyst 3.0 (http://www.metaboanalyst.ca/faces/upload/PathUploadView.xhtml) to generate the metabolome view map, which integrates enrichment analysis and pathway topology analysis. A metabolic pathway which pathway impact value is higher than 0.1 is characterized as the significantly relevant pathways.

Results

Indirect calorimetry

Compared to the fed state, the O2 consumption decreased at 24 h fasting and continued to decrease at 48 h fasting (P <  0.01) (Table 2). The CO2 and CH4 production of fasted pigs were significant lower than those of fed pigs, while those values were not significantly different between 24 h and 48 h fasting (P <  0.01). For HP and RQ, a significant decrease was observed after fasting treatment, but there were not significantly different between 24 h and 48 h fasting (P <  0.01). In addition, the UN production tended to decrease with fasting (P = 0.054).

Table 2 The effect of fasting treatment on heat production, respiratory quotient and urine nitrogen in growing pigs

Plasma metabolic profiling based on UPLC-HRMS

Principal components analysis was performed to reveal intrinsic differences within the signals (Fig. 1). The PCA score plot for plasma of pigs in the fed state, and after 24 h and 48 h of fasting showed clear clustering. PC 1 and PC 2 explained 67% of the total variances within the data. PC1 that described 52% of the variance between the fed and fasted state indicated that there were major differences in the metabolome between these two states, more than between 24 h and 48 h of fasting. The fasted samples differed along PC2 in the PCA score plots. After prolonged fasting, samples moved from the second quadrant to the third quadrant on the score plot, indicating that prolonged fasting alters the composition of the plasma metabolome.

Fig. 1
figure 1

PCA models demonstrating the separation of plasma samples of pigs under feeding, 24 h fasting and 48 h fasting. FD: Feeding; FS_24: 24 h fasting; FS_48: 48 h fasting. Each triangle represents an individual plasma sample

Based on the accurate mass measurement of molecular ions and fragment ions with high resolution, 15 compounds were identified (Table 3). These metabolites belong to different metabolic classes: fatty acid, phospholipid, amino acids, and choline. Fold change was calculated by dividing the mean of normalized intensity of each plasma metabolite in the former by the mean of normalized intensity of each plasma in the latter. A fold change > 1 indicates that the metabolite was down-regulated, whereas a fold change < 1 indicates the metabolite was up-regulated. To illustrate the direction of changes at different time points, data were visualized in the form of box-plots.

Table 3 Metabolites with significant differences among feeding, 24 h of fasting and 48 h of fastinga

Compared to the fed state, plasma 12,13-dihydroxyoctadecenoic acid (12,13-DHOME) and sphinganine levels significantly increased after fasting, slightly at 48 h of fasting and moderately at of 24 h fasting (Fig. 2). In accordance with the feeding stage, linoleic acid, stearidonic acid, oleic acid, palmitoleic acid, and tyrosine levels were low before fasting but increased and showed much larger variability after 24 h or 48 h of fasting. Also, pantothenic acid, lysophosphatidylcholine 18:0 (LysoPC (18:0)) and 5-aminopentanoic acid levels were higher at 48 h of fasting compared to those at 24 h of fasting. Not all metabolites were up-regulated during fasting. Levels of glycerophosphocholine (GPC), ornithine, and glutamine were significantly higher in the fed state compared to those at 24 h or 48 h of fasting. Plasma aminoadipic acid and betaine levels were similar in the fed state and at 48 h of fasting, but were lower at 24 h of fasting showed.

Fig. 2
figure 2

Identified compounds that change during feeding, 24 h fasting and 48 h fasting. Relative concentrations of identified compounds are presented on the Y-axis. Time points of sampling are presented on the X-axis and are defined as follows: FD: Feeding; FS24: 24 h fasting; FS48: 48 h fasting. 12,13-DHOME: 12,13-dihydroxyoctadecenoic acid; LysoPC(18:0): lysophosphatidylcholine 18:0

Pathway enrichment and pathway topology analysis were performed using MetaboAnalyst 3.0, which based on high-quality KEGG metabolic pathways as the backend knowledgebase (Fig. 3). Based on the identified metabolites and changes in their concentrations, five metabolic pathway had a pathway impact value higher than 0.1, which is the cutoff value for relevance. The five significantly relevant metabolic pathways that indicated the growing pigs’ responses to fasting treatment included: linoleic acid metabolism, pantothenate and CoA biosynthesis, arginine and proline metabolism, alanine, aspartate, and glutamate metabolism, and sphingolipid metabolism.

Fig. 3
figure 3

Topology analysis of metabolic pathways identified among the feeding, 24 h fasting and 48 h fasting comparisons. The X-axis represents the pathway impact, and Y-axis represents the pathway enrichment. Larger sizes and darker colors represent greater pathway enrichment and higher pathway impact values, respectively. I: Linoleic acid metabolism; II: Arginine and proline metabolism; III: Sphingolipid metabolism; IV: Alanine, aspartate and glutamate metabolism; V: Pantothenate and CoA biosynthesis

Discussion

The term “indirect calorimetry” refers to the fact that heat production is determined by measuring O2 consumption, and CO2 and CH4 production, which are indicative for the oxidation of metabolites [26]. During fasting, these substrates are mainly glucose, lipids, and proteins mobilized from body reserves [27]. In the present study, heat production during fasting was significantly lower than that in the fed state, which indicates that the oxidation of substrate is down-regulated. Heat production at 24 and 48 h of fasting was 823 and 733 kJ/(kg BW0.6·d), respectively. These values are within the range of fasting HP from 700 to 955 kJ/(kg BW0.6·d), which were determined in previous studies [28,29,30,31]. The variation might be partly explained by differences in the duration of fasting. The RQ is indicative for macronutrients that are being oxidized, as oxidation of glucose, lipids, and proteins yield different RQ values [32]. In the present study, the RQ at 24 and 48 h of fasting were 0.84 and 0.81, respectively. A value of 0.7 indicates that lipids are being catabolized, 0.8 for proteins, and 1.0 for glucose [33]. The only occurrence of glucose oxidation in the body can be reflected in RQ values equal to 1.0 [34]. An RQ below 1.0 indicates that lipids and protein are being catabolized [35]. Therefore, the RQ of 0.84 and 0.81 indicates that mainly lipids and protein were oxidized during fasting [36, 37]. In addition, indirect calorimetry analysis revealed that the HP, RQ and UN of growing pigs were no significant difference between 24 h fasting and 48 h fasting. These results can be explained that the energy metabolism of growing pigs under a 2-d fasting treatment was relatively stable, which were in agreement with the reports of Liu et al. [4] who also measured HP of fasted pigs using indirect calorimetry. Based on these results, a 24-h fasting treatment was more appropriate than a 48-h fasting treatment to determine the effect of fasting on the energy metabolism of growing pigs.

The body’s principal lipid classes are mainly triglycerides, phospholipids, and steroids but triglycerides are quantitatively the most important lipids [38, 39]. Humans and animals will break down lipids to meet their energy requirements when the energy supply is restricted. This is a strategy to save glucose and protein that are crucial fuels for some important organs [40]. A series of studies have reported the effect of feed restriction and fasting on lipid metabolism [41,42,43]. Compared to biochemical traits in the blood of lipid metabolism, metabolomics provides new and a more in-depth information of the global metabolite profile of plasma [44, 45].

As mentioned above, five metabolites of unsaturated fatty acids (i.e., 12,13-DHOME, linoleic acid, stearidonic acid, oleic acid and palmitoleic acid were significantly up-regulated during fasting. Similar results have been reported in other studies in animals [13]. In the fed state, unsaturated fatty acids, saturated fatty acids, and glycerol are used to synthesize triglycerides [46]. During fasting, fatty acids are released from triacylglycerol stored in adipocytes of growing pigs, resulting in increased levels in plasma [47]. Wood et al. [48] concluded that more than half of the fatty acids in animal lipids were unsaturated fatty acids, and unsaturated fatty acids were oxidized more rapidly than saturated fatty acids [49, 50]. The least oxidized of the saturated fatty acids having the greatest retention in the carcass [51]. In addition, multiple transcripts involved in the pathway of saturated fatty acids synthesis were inhibited in rat adipose tissue after short-term fasting [11]. Thus, the unsaturated fatty acids were more efficient energy substrate than saturated fatty acids, which may be explained that unsaturated fatty acids were primarily identified in the current study. Compared with the fed state, the fold change in the level of unsaturated fatty acids at 24 h and 48 h of fasting was close to 1, which indicates that the level of unsaturated fatty acids level was not significantly different at 24 h and 48 h of fasting (Table 3). The pathway of linoleic acid metabolism was identified as the most significant pathway through pathway topology analysis (Fig. 3). On the whole, the change of unsaturated fatty acid primarily contribute to fasting metabolism in our study and suggested that 24 h of fasting was more appropriate than 48 h of fasting for indicators of metabolic components of fasting metabolism.

In contrast, the GPC concentration was found to be significantly down-regulated during fasting. The GPC is a choline derivative that functions as a substrate in many bio-metabolic pathways [52]. As we know, the GPC is formed during the breakdown of phosphatidylcholine and is part of a phospholipid pathway that is active in many body tissues [53]. In the present study, fasting induced a significant decrease in the relative concentration of GPC, while the GPC level was not affected by the duration of fasting. Some studies have shown that the GPC level can be indicative for the ability to break down phospholipids as a fatty acid source to meet the energy requirements [54, 55]. Low GPC values can be accompanied by a high level of ketone bodies, which are produced by the liver from fatty acids and serve as an energy source for tissues during starvation or prolonged exercise [56, 57]. A recent study by Klein et al. [55] suggested that GPC could be used as a prognostic method for the risk of ketosis. Therefore, a decline of GPC during fasting reflects a switch of energy metabolism from phospholipids to fatty acid oxidation when the available energy is limited [55]. In addition, the almost constant levels of GPC during 24 and 48 h of fasting indicates that phospholipids metabolism of growing pigs was relatively stable during fasting.

In addition to triglycerides, phospholipids and sphingolipids are also involved in lipid metabolism. Among the identified compounds of the current study, pantothenic acid and LysoPC (18:0) participated in phospholipid metabolism. Most of these metabolites showed a significant increase during the 48 h of fasting [58, 59]. Pantothenic acid is an essential vitamin and required precursor for the biosynthesis of coenzyme A (CoA) in mammalian tissue [59]. Reibel et al. [60] reported that fasting resulted in higher tissue concentrations of Pantothenic acid, increased incorporation of Pantothenic acid into CoA, and elevated tissue concentrations of CoA in the liver. The CoA may act as an acyl group carrier in all living organisms, where it diverse cellular functions as an indispensable cofactor in energy metabolism, including the oxidation of fatty acids, carbohydrates, pyruvate, ketone bodies, and amino acids [61, 62]. The increase of Pantothenic acid can be explained that the interconversion of Pantothenic acid and CoA in the tissue was accelerated in vivo by prolonged fasting to meet the energy requirement [63].

Interestingly, we found a significant decrease in the level of betaine at 24 h of fasting. Betaine has a main physiological function as a methyl donor in the (re)formation of methionine from homocysteine [64]. A recent study of Rubio-Aliaga et al. [6] showed that methionine levels declined over 50% during fasting. The significant decline of betaine at 24 h fasting may be explained that betaine actively participates in the methionine cycle, which is a methyl donor taking part in several highly important methylation reactions [65]. However, the unexpected increase of betaine at 48 h of fasting may be explained by the observation that betaine can be maintained endogenously from choline, a process that occur mainly in the mitochondria of liver and kidney cells [66, 67].

Ornithine is a non-proteinogenic amino acid that plays an indispensable role in the urea cycle [68, 69]. In the current study, fasting induced a significant decrease in the relative concentration of ornithine, whereas ornithine levels were not affected by prolonged fasting. The change in urinary nitrogen concentration coincided with the changes in ornithine level, which indicated that ornithine was intimately associated with urinary nitrogen excretion as an important intermediate of the urea cycle in mammals [70]. Similarly, glutamine concentrations were also decreased during fasting. Rubio-Aliaga et al. [6] reported that glutamine showed negative correlations with prototypical markers of fasting such as NEFA. Glutamine is an abundant amino acid in the plasma where it functions as a non-toxic nitrogen vehicle and a respiratory fuel [71]. The decline of glutamine may partly be explained by its role as a glucose precursor during fasting by providing carbon for gluconeogenesis [16, 72]. The pathway topology analysis identified that ornithine and glutamine were enriched in the arginine and proline metabolic pathway. Recent studies reported by Pang et al. [69] indicated that proline catabolism was associated with lipid utilization by transcription factor SKN-1 during fasting. We speculated that there was a potential interaction between amino acid and lipid metabolism through those identified metabolites.

Tyrosine is not only a conditionally essential amino acid, but also a potent ketogenic amino acid [73]. Ketone bodies are formed in the liver and contributed as a fuel during fasting [74]. The level of tyrosine was significantly up-regulated during fasting in the current study. The increase of tyrosine may partly be explained by tyrosine acting as a ketogenic amino acid to meet the energy requirements of growing pigs during fasting [75, 76]. In addition, compared with amino acids level in the fed state, the fold change of tyrosine, ornithine and glutamine level at 24 h and 48 h of fasting were close to 1, which indicated that amino acids level were not significantly different at 24 h and 48 h of fasting.

Conclusions

In conclusion, a differential compound library containing 15 metabolites was identified that contributed to the differences in HP and RQ between the fed and the fasted state in growing pigs. Integrative analysis of indirect calorimetry and metabolomics profiling revealed that the decreases in HP and RQ of growing pigs under fasting conditions were associated with the alterations of linoleic acid metabolism and amino acid metabolism. The integrative analysis also indicated that growing pigs under a 24-h fasting were more appropriate than a 48-h fasting to investigate the metabolic components of maintenance energy expenditure. Our findings help to improve knowledge regarding potential mechanisms responsible for metabolic components of maintenance energy expenditure and provide possible important implications for the design of effective strategies to study fasting mechanisms.

Abbreviations

ADF:

Acid detergent fibre

DHOME:

Dihydroxyoctadecenoic acid

DM:

Dry matter

EE:

Ether extract

GE:

Gross energy

GPC:

Glycerophosphocholine

HESI:

Heated electrospray ionization

HP:

Heat production

LysoPC:

Lysophosphatidylcholine

NDF:

Neutral detergent fibre

NEFA:

Nonesterified fatty acid

PCA:

Principal component analysis

RQ:

Respiratory quotient

UN:

Nitrogen excretion in urine

References

  1. Noblet J, Fortune H, Shi XS, Dubois S. Prediction of net energy value of feeds for growing pigs. J Anim Sci. 1994;72:344–54.

    Article  PubMed  CAS  Google Scholar 

  2. Milligan L, Summers M. The biological basis of maintenance and its relevance to assessing responses to nutrients. Proc Nutr Soc. 1986;45:185–93.

    Article  PubMed  CAS  Google Scholar 

  3. de Lange K, van Milgen J, Noblet J, Dubois S, Birkett S. Previous feeding level influences plateau heat production following a 24 h fast in growing pigs. Br J Nutr. 2006;95:1082–7.

    Article  PubMed  CAS  Google Scholar 

  4. Liu DW, Jaworski NW, Zhang GF, Li ZC, Li DF, Wang FL. Effect of experimental methodology on fasting heat production and the net energy content of corn and soybean meal fed to growing pigs. Arch Anim Nutr. 2014;68:281–95.

    Article  PubMed  CAS  Google Scholar 

  5. Noblet J, Labussière E, Dubois S, de Lange C, Barea R, Lasnier J, et al. Fasting heat production and metabolic body size in non-ruminant growing farm animals. Netherlands: Wageningen Academic Publisgers; 2013. p. 313–4.

    Google Scholar 

  6. Rubio-Aliaga I, de Roos B, Duthie SJ, Crosley LK, Mayer C, Horgan G, et al. Metabolomics of prolonged fasting in humans reveals new catabolic markers. Metabolomics. 2011;7:375–87.

    Article  CAS  Google Scholar 

  7. Li ZC, Li P, Liu DW, Li DF, Wang FL, Su YB, et al. Determination of the energy value of corn distillers dried grains with solubles containing different oil levels when fed to growing pigs. J Anim Physiol Anim Nutr. 2017;101:339–48.

    Article  Google Scholar 

  8. Velayudhan D, Heo J, Nyachoti C. Net energy content of dry extruded-expelled soybean meal fed with or without enzyme supplementation to growing pigs as determined by indirect calorimetry. J Anim Sci. 2015;93:3402–9.

    Article  PubMed  CAS  Google Scholar 

  9. Inoue H, Watanuki M, Myint HT, Ito T, Kuwayama H, Hidari H. Effects of fasting and refeeding on plasma concentrations of leptin, ghrelin, insulin, growth hormone and metabolites in swine. Anim Sci J. 2005;76:367–74.

    Article  CAS  Google Scholar 

  10. Chelikani PK, Ambrose JD, Keisler DH, Kennelly JJ. Effect of short-term fasting on plasma concentrations of leptin and other hormones and metabolites in dairy cattle. Domest Anim Endocrinol. 2004;26:33–48.

    Article  PubMed  CAS  Google Scholar 

  11. Li RY, Zhang QH, Liu Z, Qiao J, Zhao SX, Shao L, et al. Effect of short-term and long-term fasting on transcriptional regulation of metabolic genes in rat tissues. Biochem Biophys Res Commun. 2006;344:562–70.

    Article  PubMed  CAS  Google Scholar 

  12. Schadewaldt P, Nowotny B, Straßburger K, Kotzka J, Roden M. Indirect calorimetry in humans: a postcalorimetric evaluation procedure for correction of metabolic monitor variability. Am J Clin Nutr. 2013;97:763–73.

    Article  PubMed  CAS  Google Scholar 

  13. Krueger R, Derno M, Goers S, Metzler-Zebeli BU, Nuernberg G, Martens K, et al. Higher body fatness in intrauterine growth retarded juvenile pigs is associated with lower fat and higher carbohydrate oxidation during ad libitum and restricted feeding. Eur J Nutr. 2013;53:583–97.

  14. Moro T, Tinsley G, Bianco A, Marcolin G, Pacelli QF, Battaglia G, et al. Effects of eight weeks of time-restricted feeding (16/8) on basal metabolism, maximal strength, body composition, inflammation, and cardiovascular risk factors in resistance-trained males. J Transl Med. 2016;14:290.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Robertson DG, Ruepp SU, Stryker SA, Hnatyshyn SY, Shipkova PA, Aranibar N, et al. Metabolomic and transcriptomic changes induced by overnight (16 h) fasting in male and female Sprague-Dawley rats. Chem Res Toxicol. 2011;24:481–7.

    Article  PubMed  CAS  Google Scholar 

  16. Hedemann MS, Damgaard BM. Metabolomic study of plasma from female mink (Neovison vison) with low and high residual feed intake during restrictive and ad libitum feeding. Comp Biochem Physiol Part D Genomics Proteomics. 2012;7:322–7.

    Article  PubMed  CAS  Google Scholar 

  17. Dunn WB. Current trends and future requirements for the mass spectrometric investigation of microbial, mammalian and plant metabolomes. Phys Biol. 2008;5:011001.

    Article  PubMed  CAS  Google Scholar 

  18. Wishart DS. Metabolomics: applications to food science and nutrition research. Trends Food Sci Technol. 2008;19:482–93.

    Article  CAS  Google Scholar 

  19. Lin G, Liu C, Feng C, Fan Z, Dai Z, Lai C, et al. Metabolomic analysis reveals differences in umbilical vein plasma metabolites between normal and growth-restricted fetal pigs during late gestation. J Nutr. 2012;142:990–8.

    Article  PubMed  CAS  Google Scholar 

  20. Bovo S, Mazzoni G, Calo DG, Galimberti G, Fanelli F, Mezzullo M, et al. Deconstructing the pig sex metabolome: targeted metabolomics in heavy pigs revealed sexual dimorphisms in plasma biomarkers and metabolic pathways. J Anim Sci. 2015;93:5681–93.

    Article  PubMed  CAS  Google Scholar 

  21. Zhang GF, Liu DW, Wang FL, Li DF. Estimation of the net energy requirements for maintenance in growing and finishing pigs. J Anim Sci. 2014;92:2987–95.

    Article  PubMed  CAS  Google Scholar 

  22. AOAC. Official methods of analysis. 18th ed. Arlington: Association of Official Chemists; 2007.

    Google Scholar 

  23. Thiex NJ, Anderson S, Gildemeister B. Crude fat, diethyl ether extraction, in feed, cereal grain, and forage (Randall/Soxtec/submersion method): collaborative study. J AOAC Int. 2003;86:888–98.

    PubMed  CAS  Google Scholar 

  24. Van Soest PV, Robertson J, Lewis B. Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition. J Dairy Sci. 1991;74:3583–97.

    Article  PubMed  CAS  Google Scholar 

  25. Brouwer E. Report of sub-committee on constants and factors. Proceedings of the 3rd EAAP Symposium on Energy Metabolism; Troonn, Publ. 11. London: Academic; 1965. p. 441–3.

  26. Simonson DC, DeFronzo RA. Indirect calorimetry: methodological and interpretative problems. Am J Physiol Endocrinol Metab. 1990;258:E399–412.

    Article  CAS  Google Scholar 

  27. Chwalibog A, Thorbek G. Quantitative partition of protein, carbohydrate and fat pools in growing pigs. Arch Anim Nutr. 1995;48:53–61.

    CAS  Google Scholar 

  28. Le Bellego L, Van Milgen J, Dubois S, Noblet J. Energy utilization of low-protein diets in growing pigs. J Anim Sci. 2001;79:1259–71.

    Article  PubMed  CAS  Google Scholar 

  29. Van Milgen J, Noblet J, Dubois S. Energetic efficiency of starch, protein and lipid utilization in growing pigs. J Nutr. 2001;131:1309–18.

    Article  PubMed  CAS  Google Scholar 

  30. Le Goff G, Dubois S, Van Milgen J, Noblet J. Influence of dietary fibre level on digestive and metabolic utilisation of energy in growing and finishing pigs. Anim Res. 2002;51:245–59.

    Article  CAS  Google Scholar 

  31. Labussiere E, van Milgen J, de Lange CF, Noblet J. Maintenance energy requirements of growing pigs and calves are influenced by feeding level. J Nutr. 2011;141:1855–61.

    Article  PubMed  CAS  Google Scholar 

  32. Widmaier EP, Raff H, Strang KT. Vander’s human physiology: the mechanisms of human body function. New York: McGraw-Hill; 2006.

    Google Scholar 

  33. Chwalibog A, Tauson AH, Thorbek G. Energy metabolism and substrate oxidation in pigs during feeding, starvation and re-feeding. J Anim Physiol Anim Nutr (Berl). 2004;88:101–12.

    Article  CAS  Google Scholar 

  34. Ferrannini E. The theoretical bases of indirect calorimetry: a review. Metabolism. 1988;37:287–301.

    Article  PubMed  CAS  Google Scholar 

  35. Chwalibog A, Jakobsen K, Henckel S, Thorbek G. Estimation of quantitative oxidation and fat retention from carbohydrate, protein and fat in growing pigs. J Anim Physiol Anim Nutr (Berl). 1992;68:123–35.

    Article  CAS  Google Scholar 

  36. Chwalibog A, Tauson AH, Thorbek G. Diurnal rhythm in heat production and oxidation of carbohydrate and fat in pigs during feeding, starvation and re-feeding. J Anim Physiol Anim Nutr. 2004;88:266–74.

    Article  CAS  Google Scholar 

  37. Chwalibog A, Jakobsen K, Tauson A-H, Thorbek G. Energy metabolism and nutrient oxidation in young pigs and rats during feeding, starvation and re-feeding. Comp Biochem Physiol A Mol Integr Physiol. 2005;140:299–307.

    Article  PubMed  CAS  Google Scholar 

  38. Sinclair AJ, Murphy K, Li D. Marine lipids: overview" news insights and lipid composition of Lyprinol". Allerg Immunol (Paris). 2000;32:261–71.

    CAS  Google Scholar 

  39. Godsland IF. Biology: risk factor modification by OCs and HRT lipids and lipoproteins. Maturitas. 2004;47:299–303.

    Article  PubMed  CAS  Google Scholar 

  40. Lomb DJ, Laurent G, Haigis MC. Sirtuins regulate key aspects of lipid metabolism. Biochim Biophys Acta. 2010;1804:1652–7.

    Article  PubMed  CAS  Google Scholar 

  41. Taskinen MR, Nikkilä EA. Effects of caloric restriction on lipid metabolism in man changes of tissue lipoprotein lipase activities and of serum lipoproteins. Atherosclerosis. 1979;32:289–99.

    Article  PubMed  CAS  Google Scholar 

  42. Rodgers JT, Puigserver P. Fasting-dependent glucose and lipid metabolic response through hepatic sirtuin 1. Proc Natl Acad Sci U S A. 2007;104:12861–6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Zemel MB, Thompson W, Milstead A, Morris K, Campbell P. Calcium and dairy acceleration of weight and fat loss during energy restriction in obese adults. Obesity. 2004;12:582–90.

    Article  CAS  Google Scholar 

  44. Nørskov NP, Hedemann MS, Lærke HN, Knudsen KEB. Multicompartmental nontargeted LC–MS metabolomics: explorative study on the metabolic responses of rye fiber versus refined wheat fiber intake in plasma and urine of hypercholesterolemic pigs. J Proteome Res. 2013;12:2818–32.

    Article  PubMed  CAS  Google Scholar 

  45. Nørskov NP, Hedemann MS, Theil PK, Knudsen KEB. Oxylipins discriminate between whole grain wheat and wheat aleurone intake: a metabolomics study on pig plasma. Metabolomics. 2013;9:464–79.

    Article  CAS  Google Scholar 

  46. Miyazaki M, Kim YC, Ntambi JM. A lipogenic diet in mice with a disruption of the stearoyl-CoA desaturase 1 gene reveals a stringent requirement of endogenous monounsaturated fatty acids for triglyceride synthesis. J Lipid Res. 2001;42:1018–24.

    PubMed  CAS  Google Scholar 

  47. Jensen MD, Ekberg K, Landau BR. Lipid metabolism during fasting. Am J Physiol Endocrinol Metab. 2001;281:E789–E93.

    Article  PubMed  CAS  Google Scholar 

  48. Wood J, Richardson R, Nute G, Fisher A, Campo M, Kasapidou E, et al. Effects of fatty acids on meat quality: a review. Meat Sci. 2004;66:21–32.

    Article  PubMed  CAS  Google Scholar 

  49. DeLany JP, Windhauser MM, Champagne CM, Bray GA. Differential oxidation of individual dietary fatty acids in humans. Am J Clin Nutr. 2000;72:905–11.

    Article  PubMed  CAS  Google Scholar 

  50. Jones PJH, Schoeller DA. Polyunsaturated-saturated ratio of diet fat influences energy substrate utilization in the human. Metabolism. 1988;37:145–51.

    Article  PubMed  CAS  Google Scholar 

  51. Leyton J, Drury PJ, Crawford MA. Differential oxidation of saturated and unsaturated fatty-acids invivo in the rat. Br J Nutr. 1987;57:383–93.

    Article  PubMed  CAS  Google Scholar 

  52. Loening NM, Chamberlin AM, Zepeda AG, Gonzalez RG, Cheng LL. Quantification of phosphocholine and glycerophosphocholine with 31P edited 1H NMR spectroscopy. NMR Biomed. 2005;18:413–20.

    Article  PubMed  CAS  Google Scholar 

  53. Iorio E, Mezzanzanica D, Alberti P, Spadaro F, Ramoni C, D'Ascenzo S, et al. Alterations of choline phospholipid metabolism in ovarian tumor progression. Cancer Res. 2005;65:9369–76.

    Article  PubMed  CAS  Google Scholar 

  54. Fang F, He X, Deng H, Chen Q, Lu J, Spraul M, et al. Discrimination of metabolic profiles of pancreatic cancer from chronic pancreatitis by high-resolution magic angle spinning 1H nuclear magnetic resonance and principal components analysis. Cancer Sci. 2007;98:1678–82.

    Article  PubMed  CAS  Google Scholar 

  55. Klein MS, Buttchereit N, Miemczyk SP, Immervoll AK, Louis C, Wiedemann S, et al. NMR metabolomic analysis of dairy cows reveals milk glycerophosphocholine to phosphocholine ratio as prognostic biomarker for risk of ketosis. J Proteome Res. 2011;11:1373–81.

    Article  PubMed  CAS  Google Scholar 

  56. Enjalbert F, Nicot M, Bayourthe C, Moncoulon R. Ketone bodies in milk and blood of dairy cows: relationship between concentrations and utilization for detection of subclinical ketosis. J Dairy Sci. 2001;84:583–9.

    Article  PubMed  CAS  Google Scholar 

  57. Newman JC, Verdin E. Ketone bodies as signaling metabolites. Trends Endocrinol Metab. 2014;25:42–52.

    Article  PubMed  CAS  Google Scholar 

  58. Elshorbagy AK, Nurk E, Gjesdal CG, Tell GS, Ueland PM, Nygard O, et al. Homocysteine, cysteine, and body composition in the Hordaland homocysteine study: does cysteine link amino acid and lipid metabolism? Am J Clin Nutr. 2008;88:738–46.

    Article  PubMed  CAS  Google Scholar 

  59. Smith CM, Song WO. Comparative nutrition of pantothenic acid. J Nutr Biochem. 1996;7:312–21.

    Article  CAS  Google Scholar 

  60. Reibel DK, Wyse BW, Berkich DA, Palko WM, Neely JR. Effects of diabetes and fasting on pantothenic-acid metabolism in rats. Am J Phys. 1981;240:E597–601.

    CAS  Google Scholar 

  61. Tahiliani AG, Beinlich CJ. Pantothenic acid in health and disease. Vitam Horm. 1991;46:165–228.

    Article  PubMed  CAS  Google Scholar 

  62. Martinez DL, Tsuchiya Y, Gout I. Coenzyme a biosynthetic machinery in mammalian cells. Biochem Soc Trans. 2014;42:1112–7.

    Article  PubMed  CAS  Google Scholar 

  63. Smith CM, Narrow CM, Kendrick ZV, Steffen C. The effect of pantothenate deficiency in mice on their metabolic response to fast and exercise. Metabolism. 1987;36:115–21.

    Article  PubMed  CAS  Google Scholar 

  64. Lever M, Slow S. The clinical significance of betaine, an osmolyte with a key role in methyl group metabolism. Clin Biochem. 2010;43:732–44.

    Article  PubMed  CAS  Google Scholar 

  65. Tan Y, Ko J, Liu X, Lu C, Li J, Xiao C, et al. Serum metabolomics reveals betaine and phosphatidylcholine as potential biomarkers for the toxic responses of processed Aconitum carmichaelii Debx. Mol BioSyst. 2014;10:2305–16.

    Article  PubMed  CAS  Google Scholar 

  66. Pekkinen J, Olli K, Huotari A, Tiihonen K, Keski-Rahkonen P, Lehtonen M, et al. Betaine supplementation causes increase in carnitine metabolites in the muscle and liver of mice fed a high-fat diet as studied by nontargeted LC-MS metabolomics approach. Mol Nutr Food Res. 2013;57:1959–68.

    Article  PubMed  CAS  Google Scholar 

  67. Craig SA. Betaine in human nutrition. Am J Clin Nutr. 2004;80:539–49.

    Article  PubMed  CAS  Google Scholar 

  68. Hallows WC, Yu W, Smith BC, Devires MK, Ellinger JJ, Someya S, et al. Sirt3 promotes the urea cycle and fatty acid oxidation during dietary restriction. Mol Cell. 2011;41:139–49.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  69. Pang S, Lynn DA, Lo JY, Paek J, Curran SP. SKN-1 and Nrf2 couples proline catabolism with lipid metabolism during nutrient deprivation. Nat Commun. 2014;5:5048.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  70. Weber AL, Miller SL. Reasons for the occurrence of the twenty coded protein amino acids. Mol Evol. 1981;17:273–84.

    Article  CAS  Google Scholar 

  71. Aledo JC. Glutamine breakdown in rapidly dividing cells: waste or investment? BioEssays. 2004;26:778–85.

    Article  PubMed  CAS  Google Scholar 

  72. Hankard RG, Haymond MW, Darmaun D. Role of glutamine as a glucose precursor in fasting humans. Diabetes. 1997;46:1535–41.

    Article  PubMed  CAS  Google Scholar 

  73. Tanaka K, Ichihara A. Control of ketogenesis from amino-acids .3. Invitro and invivo studies on ketone-body formation, lipogenesis and oxidation of tyrosine by rats. Biochim Biophys Acta. 1975;399:302–12.

    Article  PubMed  CAS  Google Scholar 

  74. Balasse EO, Féry F. Ketone body production and disposal: effects of fasting, diabetes, and exercise. Diabetes Metab Res Rev. 1989;5:247–70.

    Article  CAS  Google Scholar 

  75. Weinhouse S, Millington RH. Ketone body formation from tyrosine. J Biol Chem. 1949;181:645–53.

    PubMed  CAS  Google Scholar 

  76. Meister A. Biochemistry of the amino acids. New York: Academic Press;1957.

  77. NRC. Nutrient requirements of swine. 11th ed. Washington: Nat Acad Press; 2012.

    Google Scholar 

Download references

Acknowledgements

The authors would like to thank L. Johnston for excellent assistance in editing the manuscript.

Funding

This study was financially supported by the National Natural Science Foundation of China (31372317), Developing Key Equipment for Digital Management and Monitoring Environment in Animal Production (2013AA10230602), Prevention and Control of Nutritional Metabolism and Toxic Diseases in Livestock and Poultry (2016YFD0501204) and the 111 Project (B16044).

Availability of data and materials

All the data were presented in the main manuscript and available to readers.

Author’s contributions

HL carried out the animal trial, performed the statistics and drafted the manuscript. YFC, DXM, and JW participated in the experiments. ZL, XM, JJW, JVM and FLW critically evaluated the manuscript. All authors read and approved the final manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Fenglai Wang.

Ethics declarations

Ethics approval

All procedures used in this study were performed according to the guidelines for the ethical treatment of animal by the Institutional Animal Care and Use Committee of China Agricultural University (Beijing, China).

Competing interests

The authors declare that they have no competing interests.

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

Liu, H., Chen, Y., Ming, D. et al. Integrative analysis of indirect calorimetry and metabolomics profiling reveals alterations in energy metabolism between fed and fasted pigs. J Animal Sci Biotechnol 9, 41 (2018). https://doi.org/10.1186/s40104-018-0257-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40104-018-0257-x

Keywords