Skip to main content
Download PDF
Download ePub

The effects of short-term high-fat feeding on exercise capacity: multi-tissue transcriptome changes by RNA sequencing analysis

Lipids in Health and Disease volume 16, Article number: 28 (2017) Cite this article

Abstract

Background

The effects of short-term high fat diets on physiology are elusive and the molecular changes following fat overconsumption remain largely unknown. In this study, we aimed to evaluate exercise capacity in mice fed with a high fat diet (HFD) for 3 days and investigate the molecular mechanisms in the early response to high-fat feeding.

Methods

Exercise capacity was assessed by weight-loaded swimming test in mice fed a control diet (10 kcal% fat) or a HFD (60 kcal% fat) for 3 days. Global gene expression of ten important tissues (brain, heart, liver, spleen, lung, kidney, stomach, duodenum, skeletal muscle and blood) was analyzed using RNA Sequencing.

Results

A HFD for just 3 days can induce 71% decrease of exercise performance prior to substantial weight gain (P <0.01). Principle component analysis revealed that differential gene expression patterns existed in the ten tissues. Out of which, the brain, spleen and lung were demonstrated to have more pronounced transcriptional changes than other tissues. Biological process analysis for differentially expressed genes in the brain, spleen and lung showed that dysregulation of peripheral and central immune response had been implicated in the early stage of HFD exposure. Neurotransmission related genes and circulatory system process related genes were significantly down-regulated in the brain and lung, respectively.

Conclusions

Our findings provide new insights for the deleterious effects of high-fat feeding, especially revealing that the lung maybe as a new important target attacked by short-term high-fat feeding.

Background

High-fat and high-calorie diets along with a low physical activity lifestyle have contributed to the onset or development of type 2 diabetes, metabolic syndrome and cardiovascular disease [1]. There has been conflicting results on whether consumption of a high fat diet (HFD) is detrimental or beneficial for endurance performance. Studies in rats have demonstrated a beneficial effect of a fat-rich diet on exercise capacity via increasing the ability to oxidize fat and concomitantly sparing glycogen content [2–4]. In contrast, Murray et al. [5] reported that 9 days of high-fat feeding impaired energy production and physical performance associated with respiratory uncoupling in skeletal muscle mitochondria. In the present study, we aimed to evaluate exercise capacity in mice fed with a HFD for 3 days and investigate the molecular mechanisms in the early response to high-fat feeding.

Advances in genomic technologies may help to reveal the early molecular changes by enabling simultaneous analysis of thousands of genes in response to a HFD. The serial analysis of gene expression strategy identified 12 transcripts of hypothalamus which regulated by food intake in mice at 3 h after high-fat meal ingestion [6]. The transcriptomic analysis of duodenum mucosa after high-fat meal ingestion in C57BL/6 J mice found substantial changes of genes related to lipid metabolism [7]. Microarray analysis showed markedly changes of numerous genes involved in various biological processes including morphogenesis, fatty acid catabolism and amino acid metabolism following 3 days of high-fat feeding in the skeletal muscle of C57BL/6 J mice [8]. cDNA microarrays analysis of mRNA expression showed down-regulation of genes related to fatty acid biosynthesis in the liver of one week HFD-fed BALB/c mice [9].

However, no study to date has simultaneously analyzed the systemic gene expression profile of multi-tissues in response to short-term HFD and it remains unknown that which tissue has the most pronounced changes of gene expression profile in the early stage after high-fat feeding. Recently, RNA sequencing (RNA-seq), as an attractive alternative to microarrays for transcriptome analysis, provides major advances in robustness, comparability and richness of expression profiling data [10]. Thus we utilized RNA-seq to investigate gene expression profile of ten tissues (brain, heart, liver, spleen, lung, kidney, stomach, duodenum, skeletal muscle and blood) in C57BL/6 J mice with 3 days of high-fat feeding, which may contribute to the understanding of molecular mechanisms of changes in exercise performance induced by short-term HFD.

Methods

Animals and study protocol

Animal experiments were approved by the Animal Care and Use Committee of Southern Medical University (Approval No.2013027). The methods were carried out in accordance with the approved guidelines. Forty male C57BL/6 J mice at the age of 8 weeks were obtained from Laboratory Animal Center of Southern Medical University (Approval No. SCXK (Yue) 2011–0015). All the animals were maintained in a temperature-controlled room (22–25 °C; 35–55% humidity) with a twelve-hour light/dark cycle. Mice were randomly divided into two groups, where 20 mice were fed a control diet (CD, D12450B, 10 kcal% fat) and 20 mice were fed a high fat diet (HFD, D12492, 60 kcal% fat) for 3 days. Mice were allowed free access to food and water. The changes of body weight were observed after 3 days.

Assessment of exercise capacity

A weight-loaded swimming test has been commonly used for assessment of exercise capacity in murine [11, 12]. After 3 days, 10 mice were taken out from each group for swimming exercise performance test which was conducted as previously described with some modifications [13]. The mice were not fasted and were loaded the constant weight (1.5 g tin wire, attached to the tail). The mice were dropped individually into a swimming pool (30 cm high, 25 cm in diameter) filled with water at 25 ± 1 °C. It was considered that the mice were exhausted when they failed to return to the surface of water within a 10 s period. The swimming time to exhaustion was used as the index of exercise capacity.

Biochemical assays

After 3 days of feeding, the remaining 10 mice in each group were anesthetized with sodium pentobarbital (75 mg/kg, ip) following a 6 h fasting period. The blood samples were collected by removing the left eyeball of the mice and rapidly centrifuged at 1000 g at 4 °C for 10 min. Plasma levels of glucose, triglycerides, total cholesterol, low-density lipoprotein-cholesterol (LDL-C), high-density lipoprotein-cholesterol (HDL-C), free fatty acids (FFAs), apolipoprotein E (ApoE), C-reactive protein (CRP), superoxide dismutase (SOD), homocysteine (HCY), alanine aminotransferase (ALT), aspartate aminotransferase (AST), alkaline phosphatase (ALP), total protein (TP), albumin (ALB), globulin (GLB), ALB/GLB, total bilirubin (TBIL), direct bilirubin (DBIL), indirect bilirubin (IBIL), total bile acid (TBA), uric acid (UA), Creatinine (Cr), urea, Cystatin C (CysC), creatine kinase (CK), lactate dehydrogenase (LDH), A-hydroxybutyric acid dehydrogenase (HBDH), potassium (K), sodium (Na), chlorine (Cl) and calcium (Ca) were measured using a multifunctional biochemistry analyzer (Olympus AU2700, Tokyo, Japan). Statistical analyses for biochemical assays and assessment of exercise capacity were conducted using SPSS (version 19.0) for Windows. The data are reported as mean ± standard error of the mean (SEM). Differences between the compared groups were analyzed by Student’s t test. A P value less than 0.05 was considered to be statistically significant.

Tissue Processing and RNA Isolation

Five mice of each group were selected randomly from the mice which did not perform weight-loaded swimming test for RNA sequencing. Tissue samples included the whole brain and heart, liver, spleen, lung, kidney, stomach, duodenum, skeletal muscle and blood. The samples were dissected and immediately immersed in RNA later solution (Ambion, California, USA). All samples were stored at – 80 °C before processing. Total RNA was extracted from all samples using Trizol reagent (Invitrogen, Carlsbad, CA). The RNA concentration was quantified using a spectrophotometer (NanoDrop echnologies, Wilmington, DE) and the integrity was evaluated by the Agilent Bioanalyzer 2100 (Agilent, Santa Clara, CA).

RNA sequencing and gene expression analysis

In the CD and HFD group, fixed quantities of RNA of five samples from the same kind of tissue were combined into a single sample. The cDNA library was conducted by Illumina Tru-Seq RNA Sample Prep Kits (Illumina, San Diego, CA) with Ribosomal RNA depletion following manufacturer’s instructions. Samples were sequenced for 50 bp single read using the HiSeq2000 platform. Before alignment, reads with a low quality and adapters were screened by FastQC and removed. The remaining reads were mapped to the mice reference genome (UCSC mm10) with TopHat v2.0.9. The maximum number allowed for mismatch mapping was 2. Reads Per Kilobase of exon model per Million mapped reads (RPKM) was calculated to express the mRNA abundances. Analysis of differential expression was performed using edgeR, which could be used even with the most minimal levels of replication [14]. The read counts per gene were normalized to counts per million (CPM). CPM values were utilized for differential expression analysis, whereas RPKM values were used for principle component analysis (PCA) with the GeneSpring Gx 12.0 software (Agilent Technologies, Palo Alto, CA).

Biological process analysis

Molecule annotation system (MAS) is a set of web tools for function annotation based on integration of various public resources such as Gene Ontology, KEGG, BioCarta, GenMapp, UniGene, OMIM and more [15]. Biological process analysis for differentially expressed genes (DEGs) was performed using the CapitalBio MAS 3.0 software (CapitalBio Corporation, Beijing, China). Absolute fold change >2 with P < 0.05 was considered statistically significant in the RNA-seq analysis.

Real-time quantitative RT-PCR verification of RNA-seq data

To further confirm the findings from the RNA-seq analysis, we selectively examined 22 genes expression (8 genes in the brain, 9 genes in the spleen and 5 genes in the lung) using real-time quantitative RT-PCR (qRT-PCR) method. Five samples from the same kind of tissue of purified RNA in each group were used for qRT-PCR. Total RNA from the samples was first reverse- transcribed into cDNA templates with the PrimeScriptTM RT reagent Kit (TaKaRa, Otsu, Japan) according to the manufacturer’s instruction. PCR was run on a ABI 7500 Real-Time PCR System (Applied Biosystems, Inc., Foster City, CA, USA) using the SYBR Premix Ex TaqTM II (Otsu-Shi, Shiga, Japan). The reaction volume was 20 μL and the PCR conditions were as follows: 30 s. at 95 °C, 40 cycles of 5 s. at 95 °C and 34 s. at 60 °C, followed by a melting curve analysis step. Every sample was measured in duplicate, and relative quantification was determined by the comparative Ct method (2-ΔΔCT). β-actin was used as a housekeeping gene to normalize the expression data. The primers used for gene validation are listed in Additional file 1: Table S1.

Results

Body weight and blood plasma metabolites

As shown in Table 1, initial and final body weights showed no significantly differences among the groups. HFD feeding tended to increase weight gain, but this effect failed to reach statistical significance (P = 0.094). Plasma glucose levels were 67% higher in HFD-fed mice than CD-fed mice (P < 0.001). Although no differences were seen in plasma LDL-C, FFA and ApoE, the cholesterol level was significantly greater in the HFD-fed mice (P < 0.001). The increase in HDL-C (P < 0.001) and decrease in plasma triglycerides (P = 0.005) were also observed in HFD-fed mice. Plasma levels of CRP, SOD, HCY, ALT, AST, ALP, TP, ALB, GLB, ALB/GLB, TBIL, DBIL, IBIL, TBA, UA, Cr, Urea, CysC, LDH, HBDH, K and Ca were unchanged between the CD and HFD groups except for the CK, Na and Cl (Additional file 2: Table S2).

Table 1 Body weight and plasma biochemical parameters of CD-fed and HFD-fed mice
Full size table

Exercise capacity in a weight-loaded swimming test

The swimming time indicated the exercise capacity. Both groups of mice swam the same mean time at baseline (CD = 687.3 ± 93.9 s, HFD = 669.0 ± 87.8 s). After 3 days of feeding, CD-fed mice maintained a similar swimming time of 748.9 ± 77.4 s, whereas HFD-fed mice swam 213.7 ± 44.4 s on average, 71% less far than the CD-fed mice (P <0.001) (Fig. 1).

Fig. 1
figure 1

Effects of high-fat diet feeding on exercise capacity in C57BL/6 J mice (CD, control diet; HFD, high fat diet; *P < 0.01 vs. CD-fed mice). Values are expressed as means ± SEM

Full size image

Summary of sequencing data and global gene expression profiles

A range of 28.1 to 67.4 million raw reads were generated among samples. After removing reads with a low quality, an average of 32.3 million clean reads per sample was obtained (range, 21.3 to 51.5 million reads). Approximately 98.64% of clean reads per sample were mapped to the mice reference genome among samples (Table 2). Totally 33151 unique genes among all samples was detected. To assess the effect of sequencing depth on RNA-seq data, we conducted sequencing saturation analysis. In the beginning of the RNA-seq, with increase of the counts of reads, the number of identified genes in each tissue was increasing. However, when the counts of reads rose to approximately 30 million, the growth rate of identified genes flattened which indicated that the number of identified genes tended to saturation.

Table 2 Summary of sequence statistics
Full size table

To investigate the global gene expression profiles of ten tissues and identify the tissue with the most pronounced transcriptional changes after short-term high-fat feeding, we performed PCA on all samples (Fig. 2). The results showed differential gene expression patterns in the ten tissues. Each sphere represented an individual sample. The sphere representing liver, heart, kidney, skeletal muscle and blood in CD and HFD group overlapped, indicating that gene expression patterns of HFD-fed mice in the five tissues were almost not changed. The sphere representing stomach and duodenum in CD and HFD group were closely to each other, indicating that gene expression patterns of CD-fed and HFD-fed mice in the two tissues were nearly similar. Gene expression patterns of brain, spleen and lung in HFD group showed substantial differences as compared with CD group.

Fig. 2
figure 2

Principle component analysis of ten tissues in CD-fed and HFD-fed mice (CD, control diet; HFD, high fat diet). PCA analysis was conducted using the GeneSpring Gx 12.0 software. Each sphere represents an individual sample. The sphere representing liver, heart, kidney, skeletal muscle and blood in the two groups overlapped

Full size image

Genes and the related biological processes altered in the brain of HFD-fed mice

According to the results of PCA, we found that the brain, spleen and lung had more pronounced transcriptional changes than other tissues following 3 days of HFD intervention. Consequently, we focused on the analysis of the genes and related biological processes altered in the brain, spleen and lung of HFD-fed mice. We found 145 DEGs in the brain, of which less than half of the genes were annotated with known function from the Ensembl database (Table 3). To gain insight into the possible biologic functions of the genes affected by high-fat feeding, enrichment analysis of Gene Ontology for the DEGs was conducted. After 3 days of HFD exposure, the overrepresented biological processes in the brain were mainly enriched in neurological system process and immune response (Fig. 3). In the neurological system process related group, LIM homeobox transcription factor 1 beta (Lmx1b), and NK2 homeobox 1(Nkx2-1), which involved in neuron migration and development, were down-regulated to 8.11-fold and 10.56-fold respectively in the HFD-fed mice. Genes related to central nervous system morphogenesis were significantly down-regulated, such as homeobox D11 (Hoxd11) (12.21-fold) and UNC homeobox (Uncx) (5.39-fold). Inflammatory/immune related processes were altered as well. The mRNA levels of chemokine (C-C motif) receptor 1(Ccr1) was up-regulated to 13.18-fold. The rest immune-related genes including chemokine (C-C motif) receptor 4 (Ccr4), CD200 receptor 3 (Cd200r3), CD274 molecule (Cd274), CD300 antigen like family member G (Cd300lg) and transcription factor AP-2, alpha (Tfap2a) were down-regulated to 8.28-,3.92-, 3.48-, 11.96- and 4.79-fold respectively.

Table 3 Differentially expressed genes in the brain of HFD-fed mice
Full size table
Fig. 3
figure 3

Enriched biological process significantly altered in HFD-fed mice compared to CD-fed mice in the brain, spleen and lung (CD, control diet; HFD, high fat diet)

Full size image

Genes and the related biological processes altered in the spleen of HFD-fed mice

61 genes were markedly changed in the spleen, however, half of which were largely unknown. The overrepresented biological processes in the spleen were mainly related to acute-phase response and immune system, with a significant change in the expression of immune-related genes (Fig. 3). As shown in Table 4, immunoglobulin kappa joining 4 (Igkj4) and T cell receptor alpha joining 37(Traj37) were increased to 8.34- and 70.52-fold respectively, while Fc receptor-like S, scavenger receptor (Fcrls), immunoglobulin heavy variable 1–84 (Ighv1-84), immunoglobulin kappa joining 1(Igkj1), immunoglobulin kappa variable 5–39 (Igkv5-39), regenerating islet-derived 3 alpha (Reg3a) and regenerating islet-derived 3 beta (Reg3b) in various immunological pathways were decreased to 74.54-, 11.63-, 54.57-, 3.60-, 16.34- and 8.11- fold respectively. Other overrepresented biological processes included carbohydrate metabolism, ubiquitin-dependent protein catabolism, G2/M transition of mitotic cell cycle and neurotransmitter uptake.

Table 4 Differentially expressed genes in the spleen of HFD-fed mice
Full size table

Genes and the related biological processes altered in the lung of HFD-fed mice

In the lung, 83 genes were significantly altered. The overrepresented biological processes were mainly enriched in immune-related processes, including T cell and B cell mediated immune response (Fig. 3). As shown in Table 5, Immunoglobulin heavy constant epsilon (Ighe) and immunoglobulin heavy constant gamma 1 (G1m marker) (Ighg1) were increased to 52.71- and 6.02-fold respectively, while the mRNA levels of chemokine (C-C motif) receptor 9 (Ccr9), CD8 antigen, alpha chain (Cd8a), recombination activating gene 1(Rag1), recombination activating gene 2 (Rag2) and suppression inducing transmembrane adaptor 1 (Sit1) were decreased to 6.32-, 3.94-, 512-, 19.70- and 5.10- fold respectively. Circulatory system process was also significantly changed in the lung of HFD-fed mice, with the down-regulation of natriuretic peptide precursor A (Nppa) (7.78-fold) and natriuretic peptide precursor B (Nppb) (21.26-fold). Leptin (Lep) involved in the regulation of cholesterol absorption were increased to 5.70-fold and CART prepropeptide (Cartpt) related to cell glucose homeostasis were decreased to 62.25-fold. Other overrepresented biological processes included circadian regulation of gene expression and positive regulation of transmission of nerve impulse.

Table 5 Differentially expressed genes in the lung of HFD-fed mice
Full size table

Verification of RNA-seq data

qRT-PCR was used to validate the expression levels measured by RNA-seq for 22 selected genes (8 genes in the brain, 9 genes in the spleen and 5 genes in the lung) from the list of differently expressed genes. As demonstrated in Fig. 4, qRT-PCR showed significant alterations in the expression of the 17 genes in correspondence with the findings from the RNA-sequencing analysis, while no obvious differential expression was detected for the 5 genes of Trap2a, Cd274, Ighv1-84, Igkv5-39 and Traj37 by qRT-PCR.

Fig. 4
figure 4

Confirmation of differential gene expression via qRT-PCR analysis (CD, control diet; HFD, high fat diet; *P < 0.05 vs. CD group). The expression value was normalized to the β-actin expression level. Values are expressed as average mRNA expression ± SEM bars

Full size image

Discussion

Consistent with previous reports [6, 16], our results showed that the consumption of a high-fat diet for 3 days significantly increased plasma glucose level. The high-fat diet also increased plasma cholesterol and HDL cholesterol concentrations. Plasma triglycerides concentration decreased significantly after short-term high-fat feeding. Indeed, decreased triglycerides level was previously reported as early as 3 days after beginning a high-fat diet in a study that involved mice fed a chronic high-fat diet [17]. A randomized, double-blind, crossover study in 12 healthy subjects reported that plasma triglycerides concentration was significantly lower after a 3-d high-fat diet [18]. The low plasma triglycerides could be due to increased liver triglyceride content, possibly resulting in triglycerides being stored in the liver [19].

To our knowledge, we first reported that 3 days of high-fat feeding can induce exercise performance decrease in mice prior to substantial weight gain. In contrast to previous studies which conducted microarray analysis in a single tissue [6–9], we used RNA-seq to investigate the global gene expression profiles of ten tissues in the early response to fat intake. Interestingly, our study showed that differential gene expression patterns existed in the ten tissues. Previous researches on the effect of HFD mainly focused on the liver, skeletal muscle, intestine and heart, which were thought to have significant responses to consumption of high fat diets. However, in our results, the brain, spleen and lung were demonstrated to have more pronounced transcriptional changes than other tissues following 3 days of high-fat feeding. The simultaneous analysis of multi-tissues by RNA sequencing yielded information which had not been revealed by previous analyses of a single tissue.

In the brain of HFD-fed mice, dopamine neurons differentiation related genes including Lmx1b, Nkx2-1 and Uncx were found to be significantly down-regulated. Lmx1b, a key transcription factor for the specification of dopaminergic cell fate, has been reported to increase midbrain size and allocation of dopamine progenitors by promoting Wnt1/Wnt signaling [20]. Deficiency of Nkx2-1 in mice would lead to a remarkable abnormality in the trajectory of the ascending dopamine pathway [21]. Uncx, also known as Uncx4.1, is involved in the development of midbrain dopaminergic neurons [22]. Therefore, the down-regulation of these genes may contribute to impaired exercise capacity, which is supported by the observation of a link between dopamine and levels of physical exercise [23]. Moreover, Lmx1b and Nkx2-1 were demonstrated to regulate the migration of the superficial dorsal horn neurons and interneurons to the striatum or cortex, respectively [24, 25]. The decreased expression of the two genes suggested disturbed neuronal migration in HFD-fed mice, which may result in abnormal development of the nervous system. Inflammation is considered to be one of the important factors for deterioration of physical performance [26]. Similarly, our results also suggested that many immune-related genes were significantly altered in the brain of HFD-fed mice. Ccr1, which involved in the host response to pathogens and several inflammatory conditions [27], were significantly up-regulated. On the other hand, Ccr4 and Cd200r3 showed decreased expression in HFD-fed mice. Ccr4 were found to be functionally expressed on peripheral blood CD4+CD25+ regulatory T (Treg) cells [28]. CD200 imparts an immunoregulatory signal through the receptor for CD200, leading to the suppression of T-cell–mediated immune responses [29]. The up-regulation of Ccr1, coupled with the down-regulation of Ccr4 and Cd200r3, suggested the up-regulation of an inflammatory response toward the high fat diet in the brain. Consistently, it was reported that consumption of a HFD with 1 to 3 days induced hypothalamic inflammatory in both rats and mice [30]. However, the genes related to hypothalamic inflammatory signaling in the above study were not found to be significantly altered in our results. This discrepancy may reflect a difference between specific brain region and the whole brain.

The immune system has been considered to be affected by HFD exposure over a period of several weeks [31]. However, in the present study, we found that 3 days of HFD feeding had induced disturbed immune response in spleen, which is a major organ involved in B-cell maturation. Fcrls, which belongs to Fc receptor-like family possessing inhibitory and/or activating signaling motifs in B cell differentiation [32], were significantly down-regulated in response to HFD feeding. Similarly, microarray profiling carried out by Cui et al. in the spleen of C57BL/6 mice fed on a HFD also showed a decreased expression of Fc receptor [33]. Meanwhile, we observed abnormal expression of a few genes involved in immunoglobulin/B cell receptor signaling. Igkj4 showed significant up-regulation in HFD-fed mice, whereas the expression of Igkj1 were significantly decreased. Additionally, Reg3a and Reg3b, which play critical roles in acute-phase response [34, 35], were significantly down-regulated. The findings suggested that immune dysfunction was implicated in the spleen in response to 3 days of HFD feeding. Interestingly, neurotransmitter uptake related genes were significantly changed, such as solute carrier family 17 (sodium-dependent inorganic phosphate cotransporter), member 6(Slc17a6) which plays a key role in the transport of glutamate into synaptic vesicles before exocytotic release and the regulation of glutamate signaling [36].

Evidence has increasingly shown that a HFD regarded as a primary cause for numerous diseases including diabetes, hypertension, and steatohepatitis. However, few studies have been carried out to examine the effect of a HFD on the lung. To our surprise, the PCA results showed pronounced transcriptional changes in the lung and this is the first report to investigate the lung transcriptome profile after 3 days of HFD feeding. As the lung is a major site of immune regulation, our results revealed that many immune-related genes were significantly altered. Chemokine receptor Ccr9, which have proved to be important in the Treg cells mediated self-tolerance [37], was markedly down-regulated in HFD-fed mice. Sit1, a critical negative regulator of TCR-mediated signaling [38], showed significant down-regulation as well. The decreased expression of Ccr9 and Sit1 suggested an activated inflammatory response affected by HFD in the lung, which was in agreement with the previous study of the involvement of a HFD on lung inflammation [39]. In addition, we discovered decreased expression of CD8a in HFD-fed mice, which is important in cell-mediated immune defense and T-cell development [40]. Moreover, genes involved in T cell receptor signaling and immunoglobulin/ B cell receptor signaling were found to be significant altered. This study also revealed that circulatory system related genes were down-regulated by HFD in the lung. Nppa and Nppb are the precursor of atrial natriuretic peptide (ANP) and b-type natriuretic peptide (BNP), which have important physiological functions in the regulation of vascular tone and plasma volume [41]. ANP exhibits a protective role in the lung function in acute lung injury apart from its vasodilatory and natriuretic effects [42]. Intake of a HFD had been proved to make a slower pulmonary O2 uptake kinetics and attenuate microvascular blood flow and O2 delivery during the transition to moderate intensity exercise [43]. Therefore, we speculated that the decreased expression of Nppa and Nppb may contribute to the impaired exercise capacity in HFD-fed mice.

Conclusions

The simultaneous analysis of ten tissues following 3 days of high-fat feeding by RNA-seq technology revealed that the brain, spleen and lung had more pronounced transcriptional changes than other tissues. Dysregulation of peripheral and central immune response had been implicated in the early stage of the response to HFD exposure. Neurotransmission-related genes and circulatory system process related genes were markedly down-regulated in the brain and lung, respectively. These findings provide new insights for the deleterious effects of a HFD and contribute to the understanding of molecular mechanisms of exercise performance decrease induced by short-term high-fat feeding.

Abbreviations

ALB:

Albumin

ALP:

Alkaline phosphatase

ALT:

Alanine aminotransferase

ANP:

Atrial natriuretic peptide

ApoE:

Apolipoprotein E

AST:

Aspartate aminotransferase

BNP:

b-type natriuretic peptide

Ca:

Calcium

CD:

Control diet

CK:

Creatine kinase

Cl:

Chlorine

CPM:

Counts per million

Cr:

Creatinine

CRP:

C-reactive protein

CysC:

Cystatin C

DBIL:

Direct bilirubin

DEGs:

Differentially expressed genes

HBDH:

A-hydroxybutyric acid dehydrogenase

HCY:

Homocysteine

HDL-C:

High-density lipoprotein-cholesterol

HFD:

High fat diet

IBIL:

Indirect bilirubin

K:

Potassium

LDH:

Lactate dehydrogenas

LDL-C:

Low-density lipoprotein-cholesterol

MAS:

Molecule annotation system

Na:

Sodium

PCA:

Principle component analysis

RNA-seq:

RNA sequencing

RPKM:

Reads Per Kilobase of exon model per Million mapped reads

SOD:

Superoxide dismutase

TBA:

Total bile acid

TBIL:

Total bilirubin

TP:

Total protein

Treg:

Regulatory T

UA:

Uric acid

References

  1. Riccardi G, Giacco R, Rivellese AA. Dietary fat, insulin sensitivity and the metabolic syndrome. Clin Nutr. 2004;23:447–56.

    Article  CAS  PubMed  Google Scholar 

  2. Miller WC, Bryce GR, Conlee RK. Adaptations to a high-fat diet that increase exercise endurance in male rats. J Appl Physiol Respir Environ Exerc Physiol. 1984;56:78–83.

    CAS  PubMed  Google Scholar 

  3. Simi B, Sempore B, Mayet MH, Favier RJ. Additive effects of training and high-fat diet on energy metabolism during exercise. J Appl Physiol. 1991;71:197–203.

    CAS  PubMed  Google Scholar 

  4. Lapachet RA, Miller WC, Arnall DA. Body fat and exercise endurance in trained rats adapted to a high-fat and/or high-carbohydrate diet. J Appl Physiol. 1996;80:1173–9.

    CAS  PubMed  Google Scholar 

  5. Murray AJ, Knight NS, Cochlin LE, McAleese S, Deacon RM, Rawlins JN, et al. Deterioration of physical performance and cognitive function in rats with short-term high-fat feeding. FASEB J. 2009;23:4353–60.

    Article  CAS  PubMed  Google Scholar 

  6. De Giorgio MR, Yoshioka M, St-Amand J. Feeding induced changes in the hypothalamic transcriptome. Clin Chim Acta. 2009;406:103–7.

    Article  PubMed  Google Scholar 

  7. Yoshioka M, Bolduc C, Raymond V, St-Amand J. High-fat meal-induced changes in the duodenum mucosa transcriptome. Obesity (Silver Spring). 2008;16:2302–7.

    Article  CAS  Google Scholar 

  8. de Wilde J, Mohren R, van den Berg S, Boekschoten M, Dijk KW, de Groot P, et al. Short-term high fat-feeding results in morphological and metabolic adaptations in the skeletal muscle of C57BL/6 J mice. Physiol Genomics. 2008;32:360–9.

    Article  PubMed  Google Scholar 

  9. Nishikawa S, Sugimoto J, Okada M, Sakairi T, Takagi S. Gene expression in livers of BALB/C and C57BL/6 J mice fed a high-fat diet. Toxicol Pathol. 2012;40:71–82.

    Article  CAS  PubMed  Google Scholar 

  10. HP T, Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RH, de Menezes RX, et al. Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008;36:e141.

    Article  Google Scholar 

  11. Wang SY, Huang WC, Liu CC, Wang MF, Ho CS, Huang WP, et al. Pumpkin (Cucurbita moschata) fruit extract improves physical fatigue and exercise performance in mice. Molecules. 2012;17:11864–76.

    Article  CAS  PubMed  Google Scholar 

  12. Wu RM, Sun YY, Zhou TT, Zhu ZY, Zhuang JJ, Tang X, et al. Arctigenin enhances swimming endurance of sedentary rats partially by regulation of antioxidant pathways. Acta Pharmacol Sin. 2014;35:1274–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Tan W, Yu KQ, Liu YY, Ouyang MZ, Yan MH, Luo R, et al. Anti-fatigue activity of polysaccharides extract from Radix Rehmanniae Preparata. Int J Biol Macromol. 2012;50:59–62.

    Article  CAS  PubMed  Google Scholar 

  14. Robinson MD, McCarthy DJ, Smyth GK. edger: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.

    Article  CAS  PubMed  Google Scholar 

  15. Wang L, Sun L, Huang J, Jiang M. Cyclin-dependent kinase inhibitor 3 (CDKN3) novel cell cycle computational network between human non-malignancy associated hepatitis/cirrhosis and hepatocellular carcinoma (HCC) transformation. Cell Prolif. 2011;44:291–9.

    Article  CAS  PubMed  Google Scholar 

  16. Lee YS, Li P, Huh JY, Hwang IJ, Lu M, Kim JI, et al. Inflammation is necessary for long-term but not short-term high-fat diet-induced insulin resistance. Diabetes. 2011;60:2474–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Petit V, Arnould L, Martin P, Monnot MC, Pineau T, Besnard P, et al. Chronic high-fat diet affects intestinal fat absorption and postprandial triglyceride levels in the mouse. J Lipid Res. 2007;48:278–87.

    Article  CAS  PubMed  Google Scholar 

  18. Tremblay AJ, Lamarche B, Guay V, Charest A, Lemelin V, Couture P. Short-term, high-fat diet increases the expression of key intestinal genes involved in lipoprotein metabolism in healthy men. Am J Clin Nutr. 2013;98:32–41.

    Article  CAS  PubMed  Google Scholar 

  19. Jin ES, Beddow SA, Malloy CR, Samuel VT. Hepatic glucose production pathways after three days of a high-fat diet. Metabolism. 2013;62:152–62.

    Article  CAS  PubMed  Google Scholar 

  20. Anderegg A, Lin HP, Chen JA, Caronia-Brown G, Cherepanova N, Yun B, et al. An Lmx1b-miR135a2 regulatory circuit modulates Wnt1/Wnt signaling and determines the size of the midbrain dopaminergic progenitor pool. PLoS Genet. 2013;9:e1003973.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Kawano H, Horie M, Honma S, Kawamura K, Takeuchi K, Kimura S. Aberrant trajectory of ascending dopaminergic pathway in mice lacking Nkx2.1. Exp Neurol. 2003;182:103–12.

    Article  CAS  PubMed  Google Scholar 

  22. Rabe TI, Griesel G, Blanke S, Kispert A, Leitges M, van der Zwaag B, et al. The transcription factor Uncx4.1 acts in a short window of midbrain dopaminergic neuron differentiation. Neural Dev. 2012;7:39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Mathes WF, Nehrenberg DL, Gordon R, Hua K, Garland TJ, Pomp D, et al. Dopaminergic dysregulation in mice selectively bred for excessive exercise or obesity. Behav Brain Res. 2010;210:155–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Ding YQ, Yin J, Kania A, Zhao ZQ, Johnson RL, Chen ZF, et al. Lmx1b controls the differentiation and migration of the superficial dorsal horn neurons of the spinal cord. Development. 2004;131:3693–703.

    Article  CAS  PubMed  Google Scholar 

  25. Elias LA, Potter GB, Kriegstein AR. A time and a place for nkx2-1 in interneuron specification and migration. Neuron. 2008;59:679–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Katafuchi T, Kondo T, Take S, Yoshimura M. Enhanced expression of brain interferon-alpha and serotonin transporter in immunologically induced fatigue in rats. Eur J Neurosci. 2005;22:2817–26.

    Article  PubMed  Google Scholar 

  27. Rossi D, Zlotnik A. The biology of chemokines and their receptors. Annu Rev Immunol. 2000;18:217–42.

    Article  CAS  PubMed  Google Scholar 

  28. Iellem A, Mariani M, Lang R, Recalde H, Panina-Bordignon P, Sinigaglia F, et al. Unique chemotactic response profile and specific expression of chemokine receptors CCR4 and CCR8 by CD4(+)CD25(+) regulatory T cells. J Exp Med. 2001;194:847–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Rosenblum MD, Olasz E, Woodliff JE, Johnson BD, Konkol MC, Gerber KA, et al. CD200 is a novel p53-target gene involved in apoptosis-associated immune tolerance. Blood. 2004;103:2691–8.

    Article  CAS  PubMed  Google Scholar 

  30. Thaler JP, Yi CX, Schur EA, Guyenet SJ, Hwang BH, Dietrich MO, et al. Obesity is associated with hypothalamic injury in rodents and humans. J Clin Invest. 2012;122:153–62.

    Article  CAS  PubMed  Google Scholar 

  31. Miyazaki Y, Iwabuchi K, Iwata D, Miyazaki A, Kon Y, Niino M, et al. Effect of high fat diet on NKT cell function and NKT cell-mediated regulation of Th1 responses. Scand J Immunol. 2008;67:230–7.

    Article  CAS  PubMed  Google Scholar 

  32. Davis RS, Stephan RP, Chen CC, Dennis GJ, Cooper MD. Differential B cell expression of mouse Fc receptor homologs. Int Immunol. 2004;16:1343–53.

    Article  CAS  PubMed  Google Scholar 

  33. Cui J, Xiao Y, Shi YH, Wang B, Le GW. Lipoic acid attenuates high-fat-diet-induced oxidative stress and B-cell-related immune depression. Nutrition. 2012;28:275–80.

    Article  CAS  PubMed  Google Scholar 

  34. Wang L, Huang J, Jiang M, Lin H. Tissue-specific transplantation antigen P35B (TSTA3) immune response-mediated metabolism coupling cell cycle to postreplication repair network in no-tumor hepatitis/cirrhotic tissues (HBV or HCV infection) by biocomputation. Immunol Res. 2012;52:258–68.

    Article  CAS  PubMed  Google Scholar 

  35. van Ampting MT, Loonen LM, Schonewille AJ, Konings I, Vink C, Iovanna J, Chamaillard M, et al. Intestinally secreted C-type lectin Reg3b attenuates salmonellosis but not listeriosis in mice. Infect Immun. 2012;80:1115–20.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Brumovsky P, Watanabe M, Hokfelt T. Expression of the vesicular glutamate transporters-1 and −2 in adult mouse dorsal root ganglia and spinal cord and their regulation by nerve injury. Neuroscience. 2007;147:469–90.

    Article  CAS  PubMed  Google Scholar 

  37. Wermers JD, McNamee EN, Wurbel MA. The chemokine receptor CCR9 is required for the T-cell-mediated regulation of chronic ileitis in mice. Gastroenterology. 2011;140:1526–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Arndt B, Krieger T, Kalinski T, Thielitz A, Reinhold D, Roessner A, et al. The transmembrane adaptor protein SIT inhibits TCR-mediated signaling. PLoS One. 2011;6:e23761.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Naura AS, Hans CP, Zerfaoui M, Errami Y, Ju J, Kim H, et al. High-fat diet induces lung remodeling in ApoE-deficient mice: an association with an increase in circulatory and lung inflammatory factors. Lab Invest. 2009;89:1243–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Xu SW, Wu JY, Hu KS, Ping HL, Duan ZG, Zhang HF. Molecular cloning and expression of orange-spotted grouper (Epinephelus coioides) CD8alpha and CD8 beta genes. Fish Shellfish Immunol. 2011;30:600–8.

    Article  CAS  PubMed  Google Scholar 

  41. Baxter GF. Natriuretic peptides and myocardial ischaemia. Basic Res Cardiol. 2004;99:90–3.

    Article  CAS  PubMed  Google Scholar 

  42. Mitaka C, Hirata Y, Nagura T, Tsunoda Y, Amaha K. Beneficial effect of atrial natriuretic peptide on pulmonary gas exchange in patients with acute lung injury. Chest. 1998;114:223–8.

    Article  CAS  PubMed  Google Scholar 

  43. Raper JA, Love LK, Paterson DH, Peters SJ, Heigenhauser GJ, Kowalchuk JM. Effect of high-fat and high-carbohydrate diets on pulmonary O2 uptake kinetics during the transition to moderate-intensity exercise. J Appl Physiol. 2014;117:1371–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Science Foundation of China (Nos. 81373707, 81403447 and 81603520), the Natural Science Foundation of Guangdong Province, China (Nos. 2014A030313292, 2014A030310072 and 2016A030310084), the Science & Technical Plan of Guangzhou, Guangdong, China (No. 2014Y2-00504), the Administration of Traditional Medicine of Guangdong province (No. 20161063), the Fundamental Research Funds for the Central Universities (No. 21616315) and the Special Funds for the Cultivation of Guangdong College Students Scientific and Technological Innovation (No. pdjh2016b0093).

Availability of data and materials

All data generated or analysed during this study are included in this published article and its Additional file 1: Table S1 and Additional file 2: Table S2.

Authors’ contributions

XZ and YL conceived and designed the experiments. YX, WW, JC and PJ acquired the data. XF, XN and LC analyzed and interpreted the data. YX, WW and LC drafted the manuscript. XZ and HK revised the manuscript for important intellectual content. XZ and YL supervised the study. All authors were involved in the formulation of the research questions. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval

Animal experiments were approved by the Animal Care and Use Committee of Southern Medical University (Approval No.2013027).

Author information

Author notes

    Authors and Affiliations

    1. Department of Traditional Chinese Medicine, School of Medicine, Jinan University, Guangzhou, China

      Ya Xiao & Liguo Chen

    2. School of Traditional Chinese Medicine, Southern Medical University, Guangzhou, China

      Ya Xiao, Jieyu Chen, Pingping Jiang, Xiaoli Nie, Yanyan Liu & Xiaoshan Zhao

    3. Experimental Animal Center, Southern Medical University, Guangzhou, China

      Wanshan Wang

    4. School of Chinese Medicine, Hong Kong Baptist University, Hong Kong, China

      Xiuqiong Fu & Hiuyee Kwan

    Authors
    1. Ya Xiao
      View author publications

      You can also search for this author in PubMed Google Scholar

    2. Wanshan Wang
      View author publications

      You can also search for this author in PubMed Google Scholar

    3. Liguo Chen
      View author publications

      You can also search for this author in PubMed Google Scholar

    4. Jieyu Chen
      View author publications

      You can also search for this author in PubMed Google Scholar

    5. Pingping Jiang
      View author publications

      You can also search for this author in PubMed Google Scholar

    6. Xiuqiong Fu
      View author publications

      You can also search for this author in PubMed Google Scholar

    7. Xiaoli Nie
      View author publications

      You can also search for this author in PubMed Google Scholar

    8. Hiuyee Kwan
      View author publications

      You can also search for this author in PubMed Google Scholar

    9. Yanyan Liu
      View author publications

      You can also search for this author in PubMed Google Scholar

    10. Xiaoshan Zhao
      View author publications

      You can also search for this author in PubMed Google Scholar

    Corresponding authors

    Correspondence to Yanyan Liu or Xiaoshan Zhao.

    Additional files

    Additional file 1: Table S1.

    List of Primer sequences used for qRT-PCR. (DOC 57 kb)

    Additional file 2: Table S2.

    Plasma biochemical parameters of CD-fed and HFD-fed mice. (DOCX 21 kb)

    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

    Xiao, Y., Wang, W., Chen, L. et al. The effects of short-term high-fat feeding on exercise capacity: multi-tissue transcriptome changes by RNA sequencing analysis. Lipids Health Dis 16, 28 (2017). https://doi.org/10.1186/s12944-017-0424-7

    Download citation

    • Received:

    • Accepted:

    • Published:

    • DOI: https://doi.org/10.1186/s12944-017-0424-7

    Keywords

    Lipids in Health and Disease

    ISSN: 1476-511X