Skip to main content

Advertisement

Corn dried distillers grains with solubles (cDDGS) in the diet of pigs change the expression of adipose genes that are potential therapeutic targets in metabolic and cardiovascular diseases

Article metrics

Abstract

Background

Corn dried distillers grains with solubles (cDDGS) are a byproduct of biofuel and alcohol production. cDDGS have been used in pig feed for many years, because they are readily available and rich in protein, fiber, unsaturated fatty acids and phytosterols. However, feed mixtures too high in cDDGS result in the worsening of backfat quality. We performed RNA-sequencing analysis of backfat from crossbred pigs fed different diets. The diets were isoenergetic but contained different amounts of cDDGS and various sources of fats. The animals were divided into four dietary groups during the two months of experimentation: group I (control (-cDDGS+rapeseed oil)), group II (+cDDGS+rapeseed oil), group III (+cDDGS+beef tallow), and group IV (+cDDGS+coconut oil). The aim of the present experiment was to evaluate changes in the backfat transcriptome of pigs fed isoenergetic diets that differed in cDDGS presence.

Results

Via DESeq2 software, we identified 93 differentially expressed genes (DEGs) between groups I and II, 13 between groups I and III, and 125 between groups I and IV. DEGs identified between group I (-cDDGS+rapeseed oil) and group II (+cDDGS+rapeseed oil) were highly overrepresented in several KEGG pathways: metabolic pathways (FDR < 1.21e-06), oxidative phosphorylation (FDR < 0.00189), fatty acid biosynthesis (FDR < 0.00577), Huntington’s disease (FDR < 0.00577), fatty acid metabolism (FDR < 0.0112), Parkinson’s disease (FDR < 0.0151), non-alcoholic fatty liver disease (NAFLD) (FDR < 0.016), Alzheimer’s disease (FDR < 0.0211) and complement and coagulation cascades (FDR < 0.02).

Conclusions

We observed that the addition of cDDGS positively affects the expression of several genes that have been recently proposed as potential targets for the treatment of obesity, diabetes, cardiovascular disease, and Alzheimer’s disease (e.g., FASN, AACS, ALAS1, HMGCS1, and VSIG4). Thus, our results support the idea of including cDDGS into the diets of companion animals and humans and encourage research into the bioactive ingredients of cDDGS.

Background

The optimal feeding strategy for farm animals is one of the fundamental factors determining the profitability of meat production. cDDGS (corn dried distillers grains with solubles) has commonly been used as a replacement for soybean meal in animal feed due to its high protein and fiber content, reasonable price and availability. However, backfat from animals fed diets containing cDDGS is soft and prone to oxidation [1]. Fats containing high amounts of saturated fatty acids, such as beef tallow or coconut oil, are added to the feedstuff to counteract this deterioration in backfat quality. Such cDDGS have also been proposed for use in human diets, especially those of diabetic and celiac patients [1]. Interestingly, a decade ago, a beneficial effect of DDGS on ischemic heart disease was postulated and patented by the US Patent and Trademark Office (U.S. Patent No. 2004/0234630) [2]. Currently, with the increasing production of biofuels, the availability of cDDGS is high; however, its usage in human diets remains negligible. Understanding the molecular processes that occur in animal tissues after the consumption of cDDGS and different fats could resolve doubts surrounding the use of these compounds as a feedstuff and as part of companion animal and human nutrition.

To date, several experiments have revealed that nutrition may induce changes in the transcriptome and specific metabolic pathways [3]. Peñagaricano et al. (2014) observed noticeable gene expression differences in adipose tissue between sheep fetuses whose mothers were fed diets with cDDGS and those whose mothers were fed other isoenergetic diets during pregnancy [4]. Many of these genes were involved in adipogenesis, lipogenesis and adipose tissue development. We hypothesize that these observations may be relevant to other mammalian species. Therefore, we aimed to evaluate whether the addition of cDDGS to the diet of pigs changes the backfat transcriptome. Recently, pigs have become a promising alternative to rodent animal models due to similarities to humans with regard to organ size and physiology [5]; therefore, the results of our study may provide information concerning the possibility of using cDDGS as a health-promoting additive in the human diet.

The aims of the present study were to evaluate changes in the backfat transcriptome of pigs fed isoenergetic diets that differed in cDDGS content and sources of fat. Moreover, we intended to assess potential interactions between cDDGS and fatty acid compositions in the diet.

Results

The performance and subcutaneous backfat quality of animals after different diet treatments

Animal performance and backfat quality are described in detail elsewhere for the following dietary treatments: control - group I (-cDDGS+rapeseed oil), group II (+cDDGS+rapeseed oil), group III (+cDDGS+beef tallow), and group IV (+cDDGS+coconut oil) [6, 7]. In brief, the different dietary treatments did not affect weight gain, feed utilization, backfat thickness, or carcass meatiness. However, backfat from animals receiving cDDGS and rapeseed oil in their feed mixture (group II) displayed a lowest ratio of the sum of saturated fatty acids to the sum of unsaturated fatty acids (SFA:UFA), and the highest C18:2 linoleic acid content and iodine value (p < 0.001) (Table 1).

RNA-seq statistics

After RNA-seq analysis, we obtained from 8,659,689 to 25,360,573 filtered reads per sample (Table 2). The percentage of mapped reads was between 79,5-84,4%. Up to 5,668,453 of reads were mapped uniquely to Ensembl Sus scrofa 11.1 with an annotation version 87. The number of identified genes with normalized counts > 1 was between 13,148–14,391 depending on sample.

Differentially expressed genes (DEGs) identified by RNA-seq analysis

In the present study, we classified genes as differentially expressed when the fold change was > ± 1.3 and the adjusted p-value was < 0.05. Initially, we performed a comparison between the –cDDGS+rapeseed oil group (n = 7) and three other groups, namely, the +cDDGS+rapeseed oil (n = 6), +cDDGS+beef tallow (n = 5), and + cDDGS+coconut oil (n = 5) groups (3 analyses in total). The highest number of DEGs was identified when group I (-cDDGS+rapeseed oil) was compared to group IV (+cDDGS+coconut oil) and to group II (+cDDGS+rapeseed oil), which showed 125 DEGs with 32 overexpressed and 93 underexpressed in the -cDDGS+rapeseed oil group and 93 DEGs with 44 underexpressed and 49 overexpressed in the -cDDGS+rapeseed oil group, respectively. The number of DEGs was much lower in group I versus group III, which had 13 DEGs, 5 of which were underexpressed and 8 of which were overexpressed (Additional file 1: Table S1). Next, we combined all groups with cDDGS (group II, III and IV) and compared them to the group without cDDGS (group I). In this comparison, we identified 155 DEGs, with 59 upregulated and 96 downregulated, in the –cDDGS+rapeseed oil group (Additional file 1: Table S1). All expression data were submitted to the GEO database (accession number: GSE101433).

Table 1 Chosen performance and backfat quality parameters
Table 2 RNA-seq statistics

The comparison of DEG lists with the Venny integrative tool (http://bioinfogp.cnb.csic.es/tools/venny/index.html) revealed 39 (13,3%) common genes between (group I vs group II) and (group I vs II + III + IV), 48 (16,4%) common genes between (group I vs group IV) and (group I vs group II + III + IV) and only 3 (1%) common genes between (group I vs group III) and (group I vs group II + III + IV) (Fig. 1).

Fig. 1
figure1

The Venn diagram shows common DEGs among four comparisons: (group I vs group II), (group I vs group III), (group I vs group IV) and (group I vs II + III + IV)

Functional analysis of identified DEGs

To draw proper conclusions about the effect of each dietary factor on biological processes, molecular functions and pathways we analyzed each dataset of DEGs separately and did not analyze DEGs from (group I vs group II + III + IV) comparison. We performed several analyses to identify overrepresented DEGs using the PANTHER classification system. The most interesting results were obtained when the DEGs overexpressed in group I (-cDDGS+rapeseed oil) were compared to those in group II (+cDDGS+rapeseed oil) (47 mapped genes of the 49 DEGs). The DEGs were mainly involved in the following pathways: metabolic processes, fatty acid biosynthetic process, lipid metabolic process, and coenzyme metabolic process (Table 3). Many of these genes encoded enzymes that catalyze metabolic processes; therefore, catalytic activity and ligase activity were overrepresented among the molecular functions. Overexpressed genes were overrepresented in several Reactome pathways, of which the most significant was ChREBP, which activates metabolic gene expression (FDR < 3.94E-06), metabolism (FDR < 1.51E-10), metabolism of lipids and lipoproteins (FDR < 3.10E-05), fatty acyl-CoA biosynthesis (FDR < 4.57E-05), the citric acid (TCA) cycle and respiratory electron transport (FDR < 1.61E-02) (Table 3). No significant enrichments were observed for DEGs in the group I vs group III or the group I vs group IV comparisons.

Next, we employed STRING Software (https://string-db.org) to find networks of DEGs identified in our experiment. The network of DEGs identified in group I (–cDDGS+rapeseed oil) vs group II (+cDDGS+rapeseed oil) had significantly more interactions than expected (PPI enrichment p < 1.01E-06), while for all other comparisons, this p-value was greater than 0.05, suggesting the absence of strong interactions. However, we found that in the comparison of group I vs group IV, genes of tight junction proteins (FDR < 0.016) as well as genes of organelles and vesicles (FDR < 0.0012, FDR < 0.014) were overrepresented. We performed three analyses of the network of genes from the group I vs group II comparison: an analysis of all differentially expressed genes (90 mapped genes of 93 DEGs), an analysis of downregulated genes (48 mapped genes of 49 DEGs) (Fig. 2a), and an analysis of upregulated genes (42 mapped genes of 44 DEGs) (Fig. 2b). In all these analyses, we observed highly significant networks (p < 1.01E-06, p < 3.6E-14, p < 0.000975, respectively). The most overrepresented pathways were metabolic pathways, fatty acid biosynthesis, oxidative phosphorylation, and complement and coagulation cascades. A list of all functional enrichments is presented in Additional file 2: Table S2.

Fig. 2
figure2

Network of interacting proteins obtained after String 10.0 bioinformatics software analysis (http://string-db.org) of DEGs in (–cDDGS+rapeseed oil) vs (+cDDGS+rapeseed oil) comparison. a DEGs downregulated in (+cDDGS+rapeseed oil), b DEGs upregulated in (+cDDGS+rapeseed oil genes)

To further analyze the function of genes regulated by cDDGS, we proposed a set of high-confidence genes that were differentially expressed in both DESeq2 analyses (group I vs group II and group I vs group II + III + IV) (Table 4) and noticed that according to the literature, many are potential therapeutic targets in obesity, diabetes, and cardiovascular and neurodegenerative diseases.

Validation of RNA-seq analysis by qPCR

We performed qPCR analysis of eight genes, five of which were downregulated (ACACA, ACLY, FASN, FITM2, and ALAS) and three of which were upregulated (VSIG4, C5AR1, and MS4A2) in group II (+cDDGS+rapeseed oil) (Additional file 1: Table S1). The results of qPCR analysis for each dietary group are presented in Fig. 3. Statistical analysis with GLM procedure revealed significance of diet for ACACA, FITM2, ALAS1, MS4A2 and VSIG4 expression, while for ACLY, FASN and C5AR1 trends were observed (p < 0.074, 0.063, 0.058 respectively). Sex was a significant factor only for VSIG4 gene expression (Fig. 3). To compare fold changes obtained by RNA-seq and qPCR we used analysis with combined groups (group I) vs (group II + II + IV)(Additional file 3: Figure S1). All expression patterns obtained by qPCR were in the same direction as the patterns obtained by RNA-seq (Additional file 3: Figure S1). The Pearson correlation (R2 = 0.99) between the qPCR and RNA-seq results was highly significant, at p < 0000001.

Fig. 3
figure3

Results of qPCR analysis of selected genes in each dietary groups: group I (-cDDGS+rapeseed oil), group II (+cDDGS+rapeseed oil), group III (+cDDGS+beef tallow), group IV (+cDDGS+coconut oil) after GLM procedure with diet and sex as fixed factors. *p-value < 0.05, ** f3:2 p-value < 0.01, ns - p > 0.1

Discussion

Using cDDGS in the human diet was proposed more than a decade ago due to its beneficial impact on cardiovascular diseases and possible use as a low-cost flavor enhancer and sodium reduction enabler (US Patents: 2004/0234630 and US20080160132A1). Recently, the idea of including DDGS as an ingredient in food has been reconsidered [8, 9]. Moreover, a new method for the mitigation of mycotoxins, which are the main hazardous contaminant of cDDGS, has been developed, increasing the safety of this byproduct [10]. Nevertheless, the putative molecular mechanisms by which cDDGS prevent cardiovascular diseases have not been identified to date. Several experiments conducted to investigate gene expression changes after the inclusion of cDDGS in the diet of animals [11,12,13,14] revealed that cDDGS influence the expression of nutrient receptor genes as well as immune-associated and lipid metabolic processes. It was shown, for example, that a diet based on cDDGS increased the number of goblet cells and the intestinal expression of the MUC2, PEPT1 and CASR genes in pigs. The PEPT1 and CASR genes are responsible, for the immune response and barrier function in piglets [11], an observation that corresponds well with our results demonstrating the upregulation of genes connected to immunity in adipose tissue. Interestingly, the authors also noted a reduction in the intestinal expression of the umami taste receptor TAS1R3, which has been suggested to precipitate in metabolic regulation. Our results show that cDDGS significantly change the expression of genes associated with metabolism (FDR < 1.51E-10) (Table 3) in adipose tissue. It is tempting to speculate that these changes may be the result of processes occurring in the gut after cDDGS consumption. Moreover, many metabolism-related genes have previously been reported as potential therapeutic targets in metabolic disorders as they are associated with several putative common pathomechanisms suggested in obesity, diabetes, cancer, cardiovascular and neurodegenerative diseases (Table 4).

Table 3 Panther Classification System analysis of enrichment of Biological processes. Molecular functions and reactome pathways among overexpressed genes in the animals fed control diet (-DDGS + rapeseed oil) when compared to animals fed (+DDGS+rapeseed oil) diet (FDR) - False Discovery Rate
Table 4 Common genes differentially expressed in two comparisons: group I vs group II and group I vs group II + III + IV, related to human disease: NAFLD-Non Alcoholic Fatty Liver Disease, NASH –Non-Alcoholic SteatoHepatitis, OB-obesity, CAD-Cardiovascular Disease, T1D – diabetes Mellitus type I, T2DM – Diabetes Melitus type II, ND-Neurodegenarative Disease, ID – Immunity Disease

In our experiment, we observed the highest number of DDGS in two comparisons: -DDGS+rapeseed oil vs + cDDGS +rapeseed oil and -cDDGS +rapeseed oil vs + cDDGS+coconut oil (Additional file 1: Table S1). However, only the genes from the first comparison were highly connected to each other, as only this dataset created a significant network of interactions in STRING analysis (Fig. 3). The reason for the lack of significant interactions in group I vs group III and group I vs group IV may be the presence of two different dietary factors (cDDGS and fat), which could make the gene dataset too heterogeneous. We performed additional analysis of -cDDGS vs + cDDGS by combining all three groups with cDDGS to take advantage of a larger number of samples and to verify how many DEGs would still be significant with the additional dietary factors (beef tallow and coconut oil). Interestingly, many of the DEGs common in both analyses are connected to metabolic diseases and are considered potential therapeutic targets (e.g., FASN, RTN3) (Table 4). In contrast, the comparison of group I with group III led us to identify only 13 DEGs, which suggests that cDDGS and beef tallow may have an opposite effect on these biological pathways. Previously, it was suggested that dietary fat and fiber have antagonistic effects on adipose inflammation [15]. In our experiment, the diets with cDDGS had only 3.2–4.7% more fat than the diets without cDDGS, but the diets differed considerably in the SFA:UFA ratio and PUFA n-6 content (Table 1), which suggests that the source and fatty acid composition of dietary fat may also be important in modulating metabolic and inflammatory pathways. The comparison of -cDDGS+rapessed oil with +cDDGS+coconut oil resulted in the identification of 125 DEGs mainly engaged in metabolic processes and catalytic activity (data not shown), suggesting that coconut oil and cDDGS may have a synergistic effect on some metabolism-related genes. Nonetheless, the genes identified in this comparison did not create significant networks. In addition, there was a strong prevalence of overexpressed to underexpressed genes in the +cDDGS+coconut oil group. Therefore, it would be interesting to further investigate the transcriptomic effects of the combination of cDDGS and coconut oil.

Fatty acid biosynthesis is downregulated by cDDGS

Several genes involved in fatty acid biosynthesis were downregulated by the presence of cDDGS in the diet (FASN, ACLY, and SCD) (Additional file 1: Table S1, Additional file 2: Table S2) (Fig. 3). A key enzyme of fatty acid biosynthesis, FASN, was downregulated when group I was compared to group II and to the combined three +cDDGS groups. In the other comparisons, (group I vs group IV and group I vs group III), the differences were not statistically significant for RNA-seq, but qPCR analysis suggested that all three groups with cDDGS had decreased expression of FASN (Fig. 3). The same result was noted for the ACLY, FITM2 and ACACA genes. Simultaneously, KAT8, a gene responsible for the acetylation and degradation of FASN [16], and ZEB1, a transcription factor which is a well-known adipogenesis repressor [17, 18], were upregulated in group II and in the combined groups with cDDGS (II + III + IV) when compared to group I. Recently, it has been proposed that FASN inhibitors have therapeutic utility in diseases associated with elevated lipogenesis, such as obesity, type 2 diabetes, and NAFLD. The first FASN inhibitor successfully advanced through the drug development process and is under clinical evaluation in oncology [19]. Moreover, it was shown that certain anti-diabetic foods, like baker’s yeast glucan (BYG) or tartary buckwheat extracts, downregulate genes responsible for fatty acid biosynthesis (ACLY, ACC, and FASN.) [19,20,21,22]. Our results show that cDDGS may have similar properties. In fact, as a byproduct of alcohol production, cDDGS contains a substantial amount of yeast debris [13]. It also contains high amounts of polyunsaturated fatty acids, which are well-known inhibitors of lipogenesis through the downregulation of SCD and FASN [23, 24]. In addition, cDDGS contain high amounts of fiber, which has strong anti-lipogenic activity [9]. A significant decrease in lipogenic gene expression (ACACA, ACLY, FASN, and SCD) in chickens and rats fed a fiber-enriched diet has been recently demonstrated [25, 26]. The potential mechanism by which dietary fiber decreases lipogenesis may be connected to altered gut microbiota and is currently being intensively investigated [27].

Oxidative phosphorylation is downregulated by cDDGS

The second crucial process downregulated by cDDGS and connected to metabolic diseases was oxidative phosphorylation (Fig. 3b and Additional file 2: Table S2). In our study, we observed the upregulation of mitochondrial genes in animals fed a standard diet with no cDDGS. In particular, genes (NDUFS5, NDUFA10, NDUFB8, COX7A2, and UQCRFS1), (COX8A) and (NDUFA10, NDUFB8, COX1, and UQCRQ) were upregulated when group I was compared to groups II, IV and (II + III + IV), respectively (Additional file 1: Table S1 and Fig. 2a). Since the standard diet for pigs may be considered obesity-inducing, we speculate that the addition of cDDGS prevents diet-induced mitochondrial hyperactivation and possible mitochondrial dysfunction in the long term. In obesity and neurodegenerative diseases, impaired mitochondrial activity and biogenesis were observed [28, 29]. It has been recently demonstrated that mitochondrial alterations in the adipose tissue of mice are biphasic after diet-induced oxidative stress [30]. After two months of a high fat - high sucrose diet, mitochondrial biogenesis increased, while after 12 months, a substantial decrease was observed, accompanied by mitochondrial ultrastructure changes. Similarly, it was shown that in humans, yeast, worms, flies, mice, and monkeys that middle-age is accompanied by increased mitochondrial activity that subsequently declines at advanced ages [31]. A decrease in mitochondrial activity has been proposed as an anti-senescence strategy, since mitochondrial dysfunction and a hyperactive TCA cycle are important contributors to senescence [32, 31]. Interestingly, it has been recently shown that the KAT8 gene, which was upregulated by the addition of cDDGS in our experiment, regulates oxidative phosphorylation by controlling the expression of respiratory genes [33]. This finding suggests that deregulation of this gene may be responsible for the alterations in respiratory electron transport and fatty acid biosynthesis observed in our study.

There were several downregulated genes (ALAS1, CYB5B, PDHX, and ACLY) (Fig. 3b) that were linked to both of the previously mentioned processes: fatty acid metabolism and cellular respiration. ACLY catalyzes the synthesis of acetylcoenzyme A, which is the main substrate in fatty acid biosynthesis as well as in the TCA cycle, while ALAS 1 catalyzes heme biosynthesis for cytochrome P450, the terminal oxidase enzyme in the electron transfer chain. It has been recently proposed that heme biosynthesis is linked to adipogenesis through mitochondrial respiratory activity [34]. Researchers hypothesized that heme biosynthesis is required to achieve optimal adipocyte differentiation by sustaining mitochondrial function.

Complement and coagulation cascade are upregulated by cDDGS

Previous studies indicate that cDDGS may have a stimulatory effect on the immunity of animals [12,13,14, 35]. The results of our experiment support this hypothesis since we observed that cDDGS significantly affect adipose expression of genes connected to immunity. STRING analysis of DEGs in group I (-cDDGS+rapeseed oil) vs group II (+cDDGS+rapeseed oil) shows that the addition of cDDGS results in differential expression of complement and coagulation cascade genes (C1Q, C5AR1, F3, and F13A1) (Fig. 3a) as well as IgE receptor activity and Fc epsilon receptor complex genes (FCER1G and MS4A2) (Additional file 2: Table S2). These genes (except for FCER1G) were also downregulated in the group with no cDDGS in the diet (group I) compared to the combined +cDDGS groups (II + III + IV) (Additional file 1: Table S1). Two additional genes related to the complement system (C1QTNF6 and C1QTNF9) were differentially expressed in group I compared to the(+cDDGS +coconut oil) and (+cDDGS +beef tallow groups), respectively. Activation of complement cascades enhances the removal of bacteria and apoptotic cells and regulates inflammation. Both excessive stimulation and inhibition of this system negatively affect organisms. Chronic low-grade inflammation is considered one of the main factors contributing to the pathogenesis of obesity and neurodegenerative and cardiovascular diseases. Inflammation is initiated by polarization of M0 macrophages into pro-inflammatory M1 macrophages, while polarization to M2 macrophages has anti-inflammatory effects and is responsible for the phagocytosis of apoptotic cells. Recently, it has been shown that the VSIG4 gene, which was one of the most upregulated genes by cDDGS in our study, (Additional file 1: Table S1) inhibits macrophage-mediated inflammation, high-fat diet-induced obesity, insulin resistance, severe fibrinogen formation in the liver, mitochondrial oxidation and ROS formation [36]. Researchers demonstrated that VSIG4 reprograms mitochondrial pyruvate metabolism by increasing PDK2 expression in macrophages and the phosphorylation of pyruvate dehydrogenase (PDH) complex PDH-E. In our study, we observed that the upregulation of the VSIG4 gene in the +cDDGS groups (Additional file 1: Table S1) was accompanied by the downregulation of subunit E3 of the pyruvate dehydrogenase-PDHX complex component X, suggesting that not only the phosphorylation of PDH but also the downregulation of the PDHX-E3 binding protein subunit of the pyruvate dehydrogenase complex may contribute to the decrease in PDH activity. In fact, mutation of the PDHX gene was shown to result in pyruvate dehydrogenase deficiency [37] in Moroccan patients, supporting our hypothesis.

There were several other anti-inflammatory murine markers among the genes upregulated by cDDGS (MRC1 and FOLR1) [38] (Additional file 1: Table S1), while some pro-inflammatory markers were downregulated in this group (IL1RAP) (Additional file 1: Table S1), suggesting that cDDGS may have anti-inflammatory properties. On the other hand, C5AR1, which was also upregulated after the addition of cDDGS, is connected to inflammation. Moreover, the inhibition of C5AR1 is proposed as a therapeutic strategy in cardiovascular diseases; nonetheless, the role of C5AR1 in obesity and inflammation is still not fully understood. It was recently shown that C5AR1 plays an important role in cardiac regeneration [39] and B-1 cell homeostasis [40]. Additionally, some researchers observed downregulation of C5AR1 expression in adipose tissue with increasing body mass index in women [41], while other researchers reported upregulation of C5AR1 expression in acquired adiposity in monozygotic twins [42]. Our results suggest a strong connection between diet and macrophage polarization and adipose inflammation, but we are not able to clearly state its consequences for health. Therefore, further in vitro investigations are needed to fully evaluate the effects of the bioactive ingredients of cDDGS on these processes in adipose tissue. The same finding applies to IgE receptor activation (FCER1G, MS4A2), which is mainly triggered by the presence of parasites or bacteria and connected to allergy response.

In the case of cDDGS, the factor that may activate the complement system may be beta-glucan, which is an element of dietary fiber of grain origin that can be present also in yeast cell walls. Upregulation of inflammation and immunity genes was observed in pigs fed a high-fat high-fiber diet [15]. Researchers concluded that dietary fat and fiber may have an antagonistic effect on proinflammatory and anti-inflammatory signaling pathways. However, contrary to our results, those researchers observed immune gene modulation only in perirenal fat but not in subcutaneous fat. It is possible that cDDGS may have much stronger effects than fiber from straw, leading to gene expression modulation also in subcutaneous fat. On the other hand, such a reaction may be the result of the presence of mycotoxins in cDDGS. Biochemical analysis revealed no such contamination in our feed; however, due to the wide range of feed contaminants, it cannot be completely ruled out.

The limitation of our study is that our observation has not been proven using in vitro cultured porcine adipocytes; however, this genome-wide experiment confirmed several previous reports on in vitro and in silico interactions between regulatory genes of fatty acid biosynthesis and genes coding for enzymes engaged in this process in human and murine cells [16, 17]. Including both sexes in the experiment could also be considered a limitation; however, to assess possible bias, we performed DESeq2 analysis of 6 samples from group I vs 6 samples from group II with balanced sexes and found that the differences are negligible. The key identified genes are among statistically significant genes in both analysis (Additional file 4: Table S4). Moreover functional analysis with String software showed that genes engaged in discussed pathways (fatty acid biosynthesis, oxidative phosphorylation and complement and coagulation cascade) pathways are overrepresented in this dataset (Additional file 4: Table S4). Some doubts may arise from the fact that the diets containing cDDGS had a slightly higher fat content than the other diets, but the difference was only ~ 4%, which is not enough to individually induce such changes in the transcriptome.

Conclusions

Currently, obesity, cardiovascular and neurodegenerative diseases are some of the main causes of death in highly developed countries. This problem has also begun to include domestic animals such as dogs or cats, and cDDGS seems to be a product that could counteract these disorders. The results of our research show that the addition of cDDGS to the diet causes a reduction in the expression of genes involved in lipogenesis and cellular respiratory processes while stimulating the genes of the immune system. The use of this product in livestock diets for many years shows the relative safety of cDDGS; however, further studies are required on laboratory animals and in in vitro cultures to fully understand the molecular mechanisms that are activated after the addition of cDDGS. Moreover, isolating biologically active components of cDDGS and exploring their activity would be useful in assessing the potential of cDDGS as a health-promoting component of the diet.

Materials and methods

Animals and diets

All procedures included in this study relating to the use of live animals were in agreement with the local Ethics Committee for Experiments with Animals in Cracow (Resolution No 912, of the day 26.04.2012).

The experiment was performed at the pig farm of the Experimental Station of the National Research Institute of Animal Production on the pigs born there. National Research Institute of Animal Production was the owner of all animals used in the study. Experimental procedures regarding the nutrition, housing and management of the animals are described elsewhere [6, 7]. Briefly, 32 crossbred fatteners originating from sows (Polish Landrace × White Large Polish) mated with a boar (Duroc × Pietrain) were kept in individual straw-bedded pens in uniform conditions. The animals were healthy and as similar in regard to weight as possible. The animals were divided into four dietary groups, eight pigs in each group (four gilts and four barrows). The diets of all groups were isonitrogenous and isoenergetic and were formulated to cover the nutritional requirements of the pigs. The ingredient composition and nutritive value of the diets, as well as the fatty acid compositions of the fat sources and feed mixtures, are presented elsewhere [6]. Briefly, the diets differed among each other in terms of the presence of 20% cDDGS and a source of fat: group I = -cDDGS+ 3% rapeseed oil, group II = +cDDGS+ 3% rapeseed oil, group III = +cDDGS+ 3% beef tallow, and group IV = +cDDGS+ 3% coconut oil. The control diet in group I was characterized by a slightly lower amount of crude fat (13.7% vs 16.5–18.4%), a lower amount of fiber and a slightly higher amount of starch (45.4% vs 38.9–39.5%) compared to the diet in the other groups, while the dietary content of crude protein was similar among all groups (Table 5). Experimental fattening lasted from 60 to 118 kg of live animal weight (approximately two months).. At the end of the experiment, the pigs were transferred to a professional slaughterhouse, where they were slaughtered in accordance with the legal regulations in force at the commercial slaughterhouse. No other pharmacological anesthesia and euthanasia methods were used. The experiment was planned for the sequencing of 24 libraries, and the 7 samples most uniform in regard to weight from each dietary group of backfat were collected for transcriptome analysis (one extra for each group) and stored in a freezer (− 85 °C). The meat from the experimental animals was intended for consumption, since the dietary treatments did not contain any harmful components.

Table 5 Nutritional values and fatty acids composition (g/100 g of detected fatty acids) of feed mixtures used in the experiment

RNA isolation, RNA-seq and qPCR procedures

Total RNA was isolated using the Direct-zol RNA kit (Zymo Research, Irvine, CA) according to the manufacturer’s protocol. The quality and concentration of the obtained RNA were evaluated on a TapeStation Instrument (Agilent Technologies, Inc., Santa Clara, CA). Since we could not obtain undegraded RNA from our samples (probably due to the problems with the freezer) we decided to employ the SMARTer Stranded Total RNA Sample Prep Kit - HI Mammalian (Takara, Clontech, Mountain View, CA) dedicated for library construction from total RNA with a RIN (RNA integrity number) between 3 and 10. We obtained satisfactory RNA from 23 samples (Additional file 5: Table S3). A 260 /280 ratio was between 1.98 and 2.06. The entire protocol included several steps: ribosomal RNA depletion with RiboGone technology, first-strand cDNA synthesis and purification, Illumina-specific library amplification by PCR, and RNA-Seq library purification. The quality and concentration of the libraries were evaluated on a TapeStation 2000 Instrument (Agilent Technologies, Inc., Santa Clara, CA) and a Qubit fluorometer (Thermo Fisher Scientific, Waltham, MA). RNA libraries (n = 23 [5–7 per group]) were prepared for sequencing using standard Illumina protocols: the libraries were diluted to a final concentration of 10 nM, and 5 or 6 samples were pooled. Cluster generation was performed on a cBot instrument (Illumina, Inc., San Diego, CA) using the TruSeq SR Cluster Kit v3-cBot-HS. Sequencing (single-end) was performed in two replicates on an Illumina HiScanSQ 2000 in one flow cell with the TruSeq SBS Kit v3 - HS (50 cycles).

Validation of the RNA-seq results was performed for 8 genes (ACACA, ACLY, FASN, ALAS1, FITM2, C5AR1, VSIG4, MS4A2) by quantitative real-time PCR (qPCR). The cDNA was synthesized with the cDNA Archive Kit (Thermo Fisher Scientific, Waltham, MA). The qPCR was performed in duplicate on a QuantStudio 7 Flex instrument (Thermo Fisher Scientific, Waltham, MA) under the fast thermal profile. The reaction mix contained 1 μl of cDNA, 5 μl of GoTaq® qPCR Master Mix (Promega Corporation, Madison, WI), 0.1 μl of CXR dye, 3.23 μl of water, 0.17 μl of 60 × TaqMan assay for OAZ1 (endogenous control) amplification (Assay ID: Ss03397505_u1) and 0.5 μl of 20 × TaqMan gene expression assay for amplification of the target gene (ACACA assay ID: Ss03389963_m1, ACLY assay ID: Ss03389566_m1, FASN assay ID: Ss03386194_u1, ALAS1 assay ID: Ss04652684_m1, FITM2 assay ID: Ss03267236_m1, MS4A2 assay ID: Ss03394007_m1, VSIG4 assay ID: Ss04328828_g1, and C5AR1 assay ID: Ss03375530_u1). The relative quantitation (RQ) of each sample was calculated basing on the ΔΔCt method using QuantStudio Real-Time PCR software.

RNA-Seq data processing and statistical methods

The demultiplexing of the RNA-seq samples was performed with the bcl2fastq conversion software v1.8.4 (Illumina). The quality check, trimming of reads and mapping of reads were conducted with FastQC, FLEXBAR and TopHat software, respectively. The mapping statistics and read counts were generated with SAMStat, RSeQC, RSEM and HTSeq software. Differential expression analysis was performed using DEseq2 software. Genes with an adjusted p of < 0.05 and a fold change > ± 1.3 were regarded as differentially expressed. We decided to apply a relatively low fold-change threshold since we expected rather subtle differences in gene expression. Classification of differentially expressed genes (DEGs) was performed with the PANTHER Classification System (Bonferroni’s correction was applied, http://pantherdb.org/), and further analyses were conducted with the Reactome database (http://www.reactome.org/). To compare the RNA-seq results with the qPCR results, the Pearson correlation between the fold change obtained after RNA-seq and qPCR was calculated using SAS software. Statistical differences between phenotypic traits between the groups were evaluated using ANOVA (Duncan’s multiple range test), and for qPCR data using nonparametric ANOVA (Kruskal-Wallis test) and GLM with covariates of sex and diet to assess the significance of each factor. (SAS Enterprise guide 7.1).

Abbreviations

AACS :

Acetoacetyl-CoA synthetase

ACACA :

Acetyl-CoA carboxylase alpha

ACC :

Acetyl-CoA carboxylase

ACLY :

ATP citrate lyase

ALAS1 :

5′-Aminolevulinate synthase 1

ARVCF :

ARVCF, delta catenin family member

C1QTNF6 :

C1q and tumor necrosis factor related protein 6

C1QTNF9 :

C1q and tumor necrosis factor related protein 9

C5AR1 :

Complement component 5a receptor 1

CASR :

Calcium sensing receptor

CEP128 :

Centrosomal protein 128

ChREBP :

Carbohydrate-responsive element-binding protein

COX7A2 :

Cytochrome c oxidase subunit VIIa

COX8A :

Cytochrome c oxidase subunit VIIIA

CYB5B :

Cytochrome B5 type B (outer mitochondrial membrane)

CYGB :

Cytoglobin

CYP2U1 :

Cytochrome P450 family 2 subfamily U member 1

DSPP :

Dentin sialophosphoprotein

F13A1 :

Coagulation factor XIII, A1 polypeptide

F3 :

Coagulation Factor III (thromboplastin, tissue factor)

FAM105A :

Family with sequence similarity 105, member A

FASN :

Fatty acid synthase

FCER1G :

Fc fragment of IgE, high affinity I, receptor for; gamma polypeptide

FITM2 :

Fat storage-inducing transmembrane protein 2

GFM1 :

G elongation factor mitochondrial 1

GIMAP4 :

GTPase, IMAP family member 4

HIP1R :

Huntingtin interacting protein 1 related

HMGCS1 :

3-hydroxy-3-methylglutaryl-CoA synthase 1

KAT8 :

Lysine acetyltransferase 8

MS4A2 :

Membrane spanning 4-domains A2

MUC2 :

Mucin 2

MYO1E :

Myosin IE

NCOA7 :

Nuclear receptor coactivator 7

NDUFA10 :

NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 10

NDUFB8 :

NADH:ubiquinone oxidoreductase subunit B8

NDUFS5 :

NADH:ubiquinone oxidoreductase subunit S5

PCDH12 :

Protocadherin 12

PDHX :

Pyruvate dehydrogenase complex component X

PEPT1 :

Peptide transporter 1

RGS1 :

Regulator of G protein signaling 1

RTN3 :

Reticulon 3

SCD :

stearoyl-CoA desaturase (delta-9-desaturase)

SLC25A1 :

Solute carrier family 25 member 1

SP110 :

SP110 nuclear body protein

SPARCL1 :

SPARC like 1

TAF4B :

TATA-box binding protein associated factor 4b

TAS1R3 :

Taste 1 receptor member 3

TCHP :

Trichoplein keratin filament binding

UMPS :

Uridine monophosphate synthetase

UQCRFS1 :

Ubiquinol-cytochrome c reductase, rieske iron-sulfur polypeptide 1

UQCRQ :

Ubiquinol-cytochrome c reductase, complex III subunit VII

VSIG4 :

V-set and immunoglobulin domain containing 4

ZEB1 :

zinc finger E-box binding homeobox 1

ZMAT1 :

Zinc finger matrin-type 1

References

  1. 1.

    Liu K, Rosentrater K. Distillers Grains: Production, Properties, and Utilization CRC Press; 2016.

  2. 2.

    Sethi R. Beneficial effects of DDG in ischemic heart disease. FASEB J. 2007;21:A735–LB127.

  3. 3.

    Fenech M, El-Sohemy A, Cahill L, Ferguson LR, French TA, Tai ES, Milner J, Koh WP, Xie L, Zucker M, Buckley M, Cosgrove L, Lockett T, Fung KY, Head R. Nutrigenetics and nutrigenomics: viewpoints on the current status and applications in nutrition research and practice. J Nutrigenet Nutrigenomics. 2011;4:69–89.

  4. 4.

    Peñagaricano F, Wang X, Rosa GJ, Radunz AE, Khatib H. Maternal nutrition induces gene expression changes in fetal muscle and adipose tissues in sheep. BMC Genomics. 2014;15:1034.

  5. 5.

    Mentzel CMJ, Cardoso TF, Lex AMJ, Sørensen DB, Fredholm M, Cirera S. Fat and carbohydrate content in the diet induces drastic changes in gene expression in young Göttingen minipigs. Mamm Genome. 2017;28:166–75.

  6. 6.

    Świątkiewicz M, Oczkowicz M, Ropka-Molik K, Hanczakowska E. The effect of dietary fatty acids composition on adipose tissue quality and expression of genes related to lipid metabolism in porcine livers. Anim Feed Sci Technol. 2016;216:204–15.

  7. 7.

    Oczkowicz M, Świątkiewicz M, Ropka-Molik K, Gurgul A, Żukowski K. Effects of different sources of fat in the diet of pigs on the liver transcriptome estimated by RNA-seq. Ann Anim Sci. 2016;16(4):1073–90.

  8. 8.

    Singha P, Muthukumarappan K, Krishnan P. Influence of processing conditions on apparent viscosity and system parameters during extrusion of distiller’s dried grains-based snacks. Food Sci Nutr. 2018;6:101–10.

  9. 9.

    Pourafshar S, Rosentrater KA, Krishnan PG. Production of Barbari bread (traditional Iranian bread) using different levels of distillers dried grains with Solubles (DDGS) and sodium Stearoyl lactate (SSL). Foods. 2018;7(3):pii: E3.

  10. 10.

    Hassan YI, Zhou T. Promising detoxification strategies to mitigate mycotoxins in food and feed. Toxins. 2018;10(3):116. https://doi.org/10.3390/toxins10030116.

  11. 11.

    Saqui-Salces M, Huang Z, Vila MF, Li J, Mielke JA, Urriola PE, Shurson GC. Modulation of intestinal cell differentiation in growing pigs is dependent on the fiber source in the diet. J Anim Sci. 2017;95(3):1179–90.

  12. 12.

    Alizadeh M, Rodriguez-Lecompte JC, Echeverry H, Crow GH, Slominski BA. Effect of yeast-derived products and distillers dried grains with solubles (DDGS) on antibody-mediated immune response and gene expression of pattern recognition receptors and cytokines in broiler chickens immunized with T-cell dependent antigens. Poult Sci. 2016;95(4):823–33.

  13. 13.

    Alizadeh M, Rodriguez-Lecompte JC, Yitbarek A, Sharif S, Crow G, Slominski BA. Effect of yeast-derived products on systemic innate immune response of broiler chickens following a lipopolysaccharide challenge. Poult Sci. 2016;95(10):2266–73.

  14. 14.

    Alizadeh M, Rodriguez-Lecompte JC, Rogiewicz A, Patterson R, Slominski BA. Effect of yeast-derived products and distillers dried grains with solubles (DDGS) on growth performance, gut morphology, and gene expression of pattern recognition receptors and cytokines in broiler chickens. Poult Sci. 2016;95(3):507–17.

  15. 15.

    Gondret F, Vincent A, Houée-Bigot M, Siegel A, Lagarrigue S, Louveau I, Causeur D. Molecular alterations induced by a high-fat high-fiber diet in porcine adipose tissues: variations according to the anatomical fat location. BMC Genomics. 2016;17:120.

  16. 16.

    Lin HP, Cheng ZL, He RY, Song L, Tian MX, Zhou LS, Groh BS, Liu WR, Ji MB, Ding C, Shi YH, Guan KL, Ye D, Xiong Y. Destabilization of fatty acid synthase by acetylation inhibits De novo lipogenesis and tumor cell growth. Cancer Res. 2016;76(23):6924–36.

  17. 17.

    Gubelmann C, Schwalie PC, Raghav SK, Röder E, Delessa T, Kiehlmann E, Waszak SM, Corsinotti A, Udin G, Holcombe W, Rudofsky G, Trono D, Wolfrum C, Deplancke B. Identification of the transcription factor ZEB1 as a central component of the adipogenic gene regulatory network. elife. 2014;3:e03346.

  18. 18.

    Saykally JN, Dogan S, Cleary MP, Sanders MM. The ZEB1 transcription factor is a novel repressor of adiposity in female mice. PLoS One. 2009;4(12):e8460.

  19. 19.

    Angeles TS, Hudkins RL. Recent advances in targeting the fatty acid biosynthetic pathway using fatty acid synthase inhibitors. Expert Opin Drug Discov. 2016;11(12):1187–99.

  20. 20.

    Zhang JS, Lei JP, Wei GQ, Chen H, Ma CY, Jiang HZ. Natural fatty acid synthase inhibitors as potent therapeutic agents for cancers: a review. Pharm Biol. 2016;54(9):1919–25.

  21. 21.

    Lee MS, Shin Y, Jung S, Kim SY, Jo YH, Kim CT, Yun MK, Lee SJ, Sohn J, Yu HJ, Kim Y. The inhibitory effect of Tartary buckwheat extracts on Adipogenesis and inflammatory response. Molecules. 2017;22(7):pii: E1160.

  22. 22.

    Cao Y, Sun Y, Zou S, Li M, Xu X. Orally administered Baker's yeast β-glucan promotes glucose and lipid homeostasis in the livers of obesity and diabetes model mice. J Agric Food Chem. 2017;65(44):9665–74.

  23. 23.

    Zulkifli RM, Parr T, Salter AM, Brameld JM. Regulation of ovine and porcine stearoyl coenzyme a desaturase gene promoters by fatty acids and sterols. J Anim Sci. 2010;88:2565–75.

  24. 24.

    Kellner TA, Gabler NK, Patience JF. The composition of dietary fat alters the transcriptional profile of pathways associated with lipid metabolism in the liver and adipose tissue in the pig. J Anim Sci. 2017;95(8):3609–19. https://doi.org/10.2527/jas.2017.1658.

  25. 25.

    Alves CC, Waitzberg DL, de Andrade LS, Dos Santos Aguiar L, Reis MB, Guanabara CC, Júnior OA, Ribeiro DA, Sala P. Prebiotic and Synbiotic modifications of Beta oxidation and Lipogenic gene expression after experimental hypercholesterolemia in rat liver. Front Microbiol. 2017;8:2010.

  26. 26.

    Desert C, Baéza E, Aite M, Boutin M, Le Cam A, Montfort J, Houee-Bigot M, Blum Y, Roux PF, Hennequet-Antier C, Berri C, Metayer-Coustard S, Collin A, Allais S, Le Bihan E, Causeur D, Gondret F, Duclos MJ, Lagarrigue S. Multi-tissue transcriptomic study reveals the main role of liver in the chicken adaptive response to a switch in dietary energy source through the transcriptional regulation of lipogenesis. BMC Genomics. 2018;19(1):187.

  27. 27.

    Boulangé CL, Neves AL, Chilloux J, Nicholson JK, Dumas ME. Impact of the gut microbiota on inflammation, obesity, and metabolic disease. Genome Med. 2016;8(1):42.

  28. 28.

    Cameron RB, Beeson CC, Schnellmann RG. Development of therapeutics that induce mitochondrial biogenesis for the treatment of acute and chronic degenerative diseases. J Med Chem. 2016;59(23):10411–34.

  29. 29.

    Moreno-Viedma V, Amor M, Sarabi A, Bilban M, Staffler G, Zeyda M, Stulnig TM. Common dysregulated pathways in obese adipose tissue and atherosclerosis. Cardiovasc Diabetol. 2016;15(1):120.

  30. 30.

    Wang PW, Kuo HM, Huang HT, Chang AY, Weng SW, Tai MH, Chuang JH, Chen IY, Huang SC, Lin TK, Liou CW. Biphasic response of mitochondrial biogenesis to oxidative stress in visceral fat of diet-induced obesity mice. Antioxid Redox Signal. 2014;20(16):2572–88.

  31. 31.

    Baker DJ, Peleg S. Biphasic modeling of mitochondrial metabolism dysregulation during aging. Trends Biochem Sci. 2017;42(9):702–11.

  32. 32.

    Theurey P, Pizzo P. The aging mitochondria. Genes (Basel). 2018;9(1):pii: E22.

  33. 33.

    Chatterjee A, Seyfferth J, Lucci J, Gilsbach R, Preissl S, Böttinger L, Mårtensson CU, Panhale A, Stehle T, Kretz O, Sahyoun AH, Avilov S, Eimer S, Hein L, Pfanner N, Becker T, Akhtar A. MOF Acetyl Transferase Regulates Transcription and Respiration in Mitochondria. Cell. 2016;167(3):722–738.e23.

  34. 34.

    Moreno-Navarrete JM, Rodríguez A, Ortega F, Becerril S, Girones J, Sabater-Masdeu M, Latorre J, Ricart W, Frühbeck G, Fernández-Real JM. Heme biosynthetic pathway is functionally linked to Adipogenesis via mitochondrial respiratory activity. Obesity. 2017;25(10):1723–33.

  35. 35.

    Ruan D, Jiang SQ, Hu YJ, Ding FY, Fan QL, Chen F, Lin XJ, Li L, Wang Y. Effects of corn distillers dried grains with solubles on performance, oxidative status, intestinal immunity and meat quality of Chinese yellow broilers. J Anim Physiol Anim Nutr. 2017;101(6):1185–93.

  36. 36.

    Li J, Diao B, Guo S, Huang X, Yang C, Feng Z, Yan W, Ning Q, Zheng L, Chen Y, Wu Y. VSIG4 inhibits proinflammatory macrophage activation by reprogramming mitochondrial pyruvate metabolism. Nat Commun. 2017;8(1):1322.

  37. 37.

    Tajir M, Arnoux JB, Boutron A, Elalaoui SC, De Lonlay P, Sefiani A, Brivet M. Pyruvate dehydrogenase deficiency caused by a new mutation of PDHX gene in two Moroccan patients. Eur J Med Genet. 2012;55(10):535–40.

  38. 38.

    Jablonski KA, Amici SA, Webb LM, Ruiz-Rosado Jde D, Popovich PG, Partida-Sanchez S, Guerau-de-Arellano M. Novel markers to delineate murine M1 and M2 macrophages. PLoS One. 2015;10(12):e0145342.

  39. 39.

    Natarajan N, Abbas Y, Bryant DM, Gonzalez-Rosa JM, Sharpe M, Uygur A, Cocco-Delgado LH, Ho NN, Gerard NP, Gerard CJ, MacRae CA, Burns CE, Burns CG, Whited JL, Lee RT. Complement receptor C5AR1 plays an evolutionarily conserved role in successful cardiac regeneration. Circulation. 2018;137(20):2152–65.

  40. 40.

    Bröker K, Figge J, Magnusen AF, Manz RA, Köhl J, Karsten CM. A novel role for C5a in B-1 cell homeostasis. Front Immunol. 2018;19:9–258.

  41. 41.

    Poursharifi P, Rezvani R, Gupta A, Lapointe M, Marceau P, Tchernof A, Cianflone K. Association of immune and metabolic receptors C5aR and C5L2 with adiposity in women. Mediat Inflamm. 2014;2014:413921.

  42. 42.

    Kaye S, Lokki AI, Hanttu A, Nissilä E, Heinonen S, Hakkarainen A, Lundbom J, Lundbom N, Saarinen L, Tynninen O, Muniandy M, Rissanen A, Kaprio J, Meri S, Pietiläinen KH. Upregulation of early and downregulation of terminal pathway complement genes in subcutaneous adipose tissue and adipocytes in acquired obesity. Front Immunol. 2017;8:545.

  43. 43.

    de Conti A, Dreval K, Tryndyak V, Orisakwe OE, Ross SA, Beland FA, Pogribny IP. Inhibition of the cell death pathway in nonalcoholic Steatohepatitis (NASH)-related Hepatocarcinogenesis is associated with histone H4 lysine 16 deacetylation. Mol Cancer Res. 2017;15(9):1163–72.

  44. 44.

    Chen QY, Costa M, Sun H. Structure and function of histone acetyltransferase MOF. AIMS biophysics. 2015;2(4):555–69.

  45. 45.

    Doğan A, Demirci S, Kıratlı B, Şahin F. Cytoglobin: a potential marker for adipogenic differentiation in preadipocytes in vitro. Cytotechnology. 2017;69(1):157–65.

  46. 46.

    Lv Y, Wang Q, Diao Y, Xu R. Cytoglobin: a novel potential gene medicine for fibrosis and cancer therapy. Curr Gene Ther. 2008;8(4):287–94.

  47. 47.

    Heinonen MT, Laine AP, Söderhäll C, Gruzieva O, Rautio S, Melén E, Pershagen G, Lähdesmäki HJ, Knip M, Ilonen J, Henttinen TA, Kere J, Lahesmaa R. Finnish Pediatric Diabetes Registry. GIMAP GTPase family genes: potential modifiers in autoimmune diabetes, asthma, and allergy. J Immunol. 2015;194(12):5885–94.

  48. 48.

    Qu Z, D'Mello SR. Proteomic analysis identifies NPTX1 and HIP1R as potential targets of histone deacetylase-3-mediated neurodegeneration. Exp Biol Med (Maywood). 2018;243(7):627–38.

  49. 49.

    Pedersen HK, Gudmundsdottir V, Brunak S. Pancreatic islet protein complexes and their dysregulation in type 2 diabetes. Front Genet. 2017;8:43.

  50. 50.

    Philibert C, Bouillot S, Huber P, Faury G. Protocadherin-12 deficiency leads to modifications in the structure and function of arteries in mice. Pathol Biol (Paris). 2012;60(1):34–40.

  51. 51.

    Hawksworth OA, Li XX, Coulthard LG, Wolvetang EJ, Woodruff TM. New concepts on the therapeutic control of complement anaphylatoxin receptors. Mol Immunol. 2017;89:36–43.

  52. 52.

    Small AG, Al-Baghdadi M, Quach A, Hii C, Ferrante A. Complement receptor immunoglobulin: a control point in infection and immunity, inflammation and cancer. Swiss Med Wkly. 2016;146:w14301.

  53. 53.

    Pan H, Yan BS, Rojas M, Shebzukhov YV, Zhou H, Kobzik L, Higgins DE, Daly MJ, Bloom BR, Kramnik I. Ipr1 gene mediates innate immunity to tuberculosis. Nature. 2005;434(7034):767–72.

  54. 54.

    Gagliardi F, Narayanan A, Mortini P. SPARCL1 a novel player in cancer biology. Crit Rev Oncol Hematol. 2017;109:63–8.

  55. 55.

    Timofeeva AV, Goriunova LE, Khaspekov GL, Il'inskaia OP, Sirotkin VN, Andreeva ER, Tararak EM, Bulkina OS, Buza VV, Britareva VV, Karpov IA, Bibilashvili RS. Comparative transcriptome analysis of human aorta atherosclerotic lesions and peripheral blood leukocytes from essential hypertension patients. Kardiologiia. 2009;49(9):27–38.

  56. 56.

    Vafadar-Isfahani B, Ball G, Coveney C, Lemetre C, Boocock D, Minthon L, Hansson O, Miles AK, Janciauskiene SM, Warden D, Smith AD, Wilcock G, Kalsheker N, Rees R, Matharoo-Ball B, Morgan K. Identification of SPARC-like 1 protein as part of a biomarker panel for Alzheimer's disease in cerebrospinal fluid. J Alzheimers Dis. 2012;28(3):625–36.

  57. 57.

    Rangaraju S, Dammer EB, Raza SA, Rathakrishnan P, Xiao H, Gao T, Duong DM, Pennington MW, Lah JJ, Seyfried NT, Levey AI. Identification and therapeutic modulation of a pro-inflammatory subset of disease-associated-microglia in Alzheimer's disease. Mol Neurodegener. 2018;13(1):24.

  58. 58.

    Patel J, Chuaiphichai S, Douglas G, Gorvin CM, Channon KM. Vascular wall regulator of G-protein signalling-1 (RGS-1) is required for angiotensin II-mediated blood pressure control. Vasc Pharmacol. 2018;S1537-1891(17):30111–8.

  59. 59.

    Cai C, Yu H, Huang G, Du X, Yu X, Zhou Y, Shen W. Histone modifications in fatty acid synthase modulated by carbohydrate responsive element binding protein are associated with non-alcoholic fatty liver disease. Int J Mol Med. 2018;42(3):1215–228.

  60. 60.

    Catalina-Rodriguez O, Kolukula VK, Tomita Y, Preet A, Palmieri F, Wellstein A, Byers S, Giaccia AJ, Glasgow E, Albanese C, Avantaggiati ML. The mitochondrial citrate transporter, CIC, is essential for mitochondrial homeostasis. Oncotarget. 2012;3(10):1220–35.

  61. 61.

    Ding J, Reynolds LM, Zeller T, Müller C, Lohman K, Nicklas BJ, Kritchevsky SB, Huang Z, de la Fuente A, Soranzo N, Settlage RE, Chuang CC, Howard T, Xu N, Goodarzi MO, Chen YD, Rotter JI, Siscovick DS, Parks JS, Murphy S, Jacobs DR Jr, Post W, Tracy RP, Wild PS, Blankenberg S, Hoeschele I, Herrington D, McCall CE, Liu Y. Alterations of a cellular cholesterol metabolism network are a molecular feature of obesity-related type 2 diabetes and cardiovascular disease. Diabetes. 2015;64(10):3464–74.

  62. 62.

    Fader KA, Nault R, Kirby MP, Markous G, Matthews J, Zacharewski TR. Convergence of hepcidin deficiency, systemic iron overloading, heme accumulation, and REV-ERBα/β activation in aryl hydrocarbon receptor-elicited hepatotoxicity. Toxicol Appl Pharmacol. 2017;321:1–17.

  63. 63.

    Steidle-Kloc E, Schönfelder M, Müller E, Sixt S, Schuler G, Patsch W, Niebauer J. Does exercise training impact clock genes in patients with coronary artery disease and type 2 diabetes mellitus? Eur J Prev Cardiol. 2016;23(13):1375–82.

  64. 64.

    Xiang R, Fan LL, Huang H, Chen YQ, He W, Guo S, Li JJ, Jin JY, Du R, Yan R, Xia K. Increased RTN3 Leads to Obesity and Hypertriglyceridemia by Interacting with HSPA5. Circulation. 2018; pii: CIRCULATIONAHA;117:030718.

  65. 65.

    Chen R, Jin R, Wu L, Ye X, Yang Y, Luo K, Wang W, Wu D, Ye X, Huang L, Huang T, Xiao G. Reticulon 3 attenuates the clearance of cytosolic prion aggregates via inhibiting autophagy. Autophagy. 2011;7(2):205–16.

  66. 66.

    Tsai MT, Chen YJ, Chen CY, Tsai MH, Han CL, Chen YJ, Mersmann HJ, Ding ST. Identification of potential plasma biomarkers for nonalcoholic fatty liver disease by integrating transcriptomics and proteomics in laying hens. J Nutr. 2017;147(3):293–303.

  67. 67.

    Perkins M, Wolf AB, Chavira B, Shonebarger D, Meckel JP, Leung L, Ballina L, Ly S, Saini A, Jones TB, Vallejo J, Jentarra G, Valla J. Altered energy metabolism pathways in the posterior cingulate in young adult apolipoprotein E ɛ4 carriers. J Alzheimers Dis. 2016;53(1):95–106.

  68. 68.

    Wei P, Milbauer LC, Enenstein J, Nguyen J, Pan W, Hebbel RP. Differential endothelial cell gene expression by African Americans versus Caucasian Americans: a possible contribution to health disparity in vascular disease and cancer. BMC Med. 2011;9:2.

  69. 69.

    McReynolds MR, Wang W, Holleran LM, Hanna-Rose W. Uridine monophosphate synthetase enables eukaryotic de novo NAD+ biosynthesis from quinolinic acid. J Biol Chem. 2017;292(27):11147–53.

  70. 70.

    Bardeleben C, Sharma S, Reeve JR, Bassilian S, Frost P, Hoang B, Shi Y, Lichtenstein A. Metabolomics identifies pyrimidine starvation as the mechanism of 5-aminoimidazole-4-carboxamide-1-β-riboside-induced apoptosis in multiple myeloma cells. Mol Cancer Ther. 2013;12(7):1310–21.

  71. 71.

    Griffith M, Mwenifumbo JC, Cheung PY, Paul JE, Pugh TJ, Tang MJ, Chittaranjan S, Morin RD, Asano JK, Ally AA, Miao L, Lee A, Chan SY, Taylor G, Severson T, Hou YC, Griffith OL, Cheng GS, Novik K, Moore R, Luk M, Owen D, Brown CJ, Morin GB, Gill S, Tai IT, Marra MA. Novel mRNA isoforms and mutations of uridine monophosphate synthetase and 5-fluorouracil resistance in colorectal cancer. Pharmacogenomics J. 2013;13(2):148–58.

  72. 72.

    Sridhar GR, Rao AA, Srinivas K, Nirmala G, Lakshmi G, Suryanarayna D, Rao PV, Kaladhar DG, Kumar SV, Devi TU, Nitesh T, Hanuman T. Butyrylcholinesterase in metabolic syndrome. Med Hypotheses. 2010;75(6):648–51.

  73. 73.

    Miranda DA, Kim JH, Nguyen LN, Cheng W, Tan BC, Goh VJ, Tan JS, Yaligar J, Kn BP, Velan SS, Wang H, Silver DL. Fat storage-inducing transmembrane protein 2 is required for normal fat storage in adipose tissue. J Biol Chem. 2014;289(14):9560–72.

Download references

Acknowledgments

Not applicable.

Funding

This research was financed by the National Science Center Poland (grant no 2014/13/B/NZ9/02134).

Availability of data and materials

All expression data were submitted to the GEO database: https://www.ncbi.nlm.nih.gov/geo/ (accession number: GSE101433).

Author information

MO, MŚ, KP, AG performed experiments; TS, MO, AG, analyzed data; MO, TZ interpreted results of experiments; MO prepared figures; MO drafted manuscript; MO, MŚ, TS, KP, TZ, AG edited and revised manuscript; MO, MŚ, TS, KP, TZ, AG approved final version of manuscript; MŚ, MO, conceived and designed research. All authors read and approved the final manuscript.

Correspondence to Maria Oczkowicz.

Ethics declarations

Ethics approval and consent to participate

All procedures included in this study relating to the use of live animals were in agreement with the local Ethics Committee for Experiments with Animals in Cracow (Resolution No 912, of the day 26.04.2012).

Consent for publication

Not applicable.

Competing interests

Author Artur Gurgul is a member of the editorial board (Associate Editor) of BMC Genomics Journal, however, was not involved in any stage of editorial manuscript assessment and did not influence the review process. Other authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Table S1. Results of the identification of DEGs by DESeq 2 software of all comparisons. (XLSX 98 kb)

Additional file 2:

Table S2. List of all enrichments identified by String software among DEGs in (-cDDGS+rapessed oil) vs (+cDDGS+rapeseed oil) comparison. (XLSX 58 kb)

Additional file 3:

Figure S1. Validation of RNA-seq results by qPCR.– comparison of fold changes between group I (n = 7) vs group II + III + IV (n = 16) obtained after RNA-seq and qPCR. (JPG 32 kb)

Additional file 4:

Table S4. Results of the identification of DEGs between groups I and II by DESeq 2 software and functional analysis with String software of identified DEGs after excluding one sample to obtain sex balance among groups. (XLSX 6539 kb)

Additional file 5:

Table S3. Characteristics of samples used in the study. (XLSX 10 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

Verify currency and authenticity via CrossMark

Cite this article

Oczkowicz, M., Szmatoła, T., Świątkiewicz, M. et al. Corn dried distillers grains with solubles (cDDGS) in the diet of pigs change the expression of adipose genes that are potential therapeutic targets in metabolic and cardiovascular diseases. BMC Genomics 19, 864 (2018) doi:10.1186/s12864-018-5265-x

Download citation

Keywords

  • Obesity
  • Pigs
  • cDDGS
  • RNA-seq
  • Backfat
  • Metabolic diseases
  • CVD
  • Alzheimer’s