Skip to main content
Download PDF

Genome-wide expression profiling in muscle and subcutaneous fat of lambs in response to the intake of concentrate supplemented with vitamin E

BMC Genomics volume 18, Article number: 92 (2017) Cite this article

Abstract

Background

The objective of this study was to acquire a broader, more comprehensive picture of the transcriptional changes in the L. Thoracis muscle (LT) and subcutaneous fat (SF) of lambs supplemented with vitamin E. Furthermore, we aimed to identify novel genes involved in the metabolism of vitamin E that might also be involved in meat quality. In the first treatment, seven lambs were fed a basal concentrate from weaning to slaughter (CON). In the second treatment, seven lambs received basal concentrate from weaning to 4.71 ± 2.62 days and thereafter concentrate supplemented with 500 mg dl-α-tocopheryl acetate/kg (VE) during the last 33.28 ± 1.07 days before slaughter.

Results

The addition of vitamin E to the diet increased the α-tocopherol muscle content and drastically diminished the lipid oxidation of meat. Gene expression profiles for treatments VE and CON were clearly separated from each other in the LT and SF. Vitamin E supplementation had a dramatic effect on subcutaneous fat gene expression, showing general up-regulation of significant genes, compared to CON treatment. In LT, vitamin E supplementation caused down-regulation of genes related to intracellular signaling cascade. Functional analysis of SF showed that vitamin E supplementation caused up-regulation of the lipid biosynthesis process, cholesterol, and sterol and steroid biosynthesis, and it down-regulated genes related to the stress response.

Conclusions

Different gene expression patterns were found between the SF and LT, suggesting tissue specific responses to vitamin E supplementation. Our study enabled us to identify novel genes and metabolic pathways related to vitamin E metabolism that might be implicated in meat quality. Further exploration of these genes and vitamin E could lead to a better understanding of how vitamin E affects the oxidative process that occurs in manufactured meat products.

Background

Vitamin E is a lipid-soluble essential component of human and animal diets due to its powerful antioxidant activity. There are eight different isoforms of vitamin E, and the most active isoform for cell protection appears to be α-tocopherol [1]. Alpha-tocopherol is a cell signaling molecule involved in a broad range of effects on cellular systems [2] and specific regulation of gene expression [2–6]; therefore, α-tocopherol can influence a number of biological functions by regulating cell signaling at both the mRNA and microRNA (miRNA) levels [7]. So far, the impact and beneficial effects of vitamin E in the prevention of chronic diseases have been mainly associated with its antioxidant properties [8].

Vitamin E is extensively used in animal diets as an antioxidant supplement. Feeding strategy is an important tool to manipulate intramuscular α-tocopherol content in animals. An increase in α-tocopherol intake can be achieved through grazing systems [9] or by feeding animals with concentrate supplemented with α-tocopherol, which can be used as an effective method to reduce the oxidative processes reported in meat products [10, 11].

During oxidative processes, the muscle haeminic pigment changes from red oxymyoglobin to brown metmyoglobin, giving the meat an undesirable brownish color [9]. Moreover, lipid oxidation results in the production of free radicals, which are linked to the formation of off-flavors and odors, a reduction in polyunsaturated fatty acids and the production of undesirable compounds, such as potentially toxic peroxides and aldehydes [12]. All of these modifications cause decreases in the freshness and quality of meat and lower consumer acceptance, resulting in considerable economic loss for the meat industry.

So far, efforts to discover the genes and metabolic pathways involved in vitamin E metabolism have mainly focused on rodents. These studies have reported important long-term effects of vitamin E deficiency on liver gene expression, up-regulation of genes involved in cholesterol synthesis and steroidogenesis, lipid uptake, and anti-oxidative mechanisms [3, 13]. In addition, vitamin E in the rat brain affected genes associated with hormones, nerve growth, apoptosis, dopaminergic neurotransmission, and clearance of amyloid-α [14]. Furthermore, genes encoding for muscle structure and extra cellular matrix and those involved in anti-oxidative and anti-inflammatory processes were up-regulated in rat skeletal muscle in response to vitamin E deficiency [15].

Although nutritional science has embraced the tools of genomics, including cDNA arrays in rodent models, few attempts at large-scale or global evaluation of nutritional gene regulation have mainly used rodents as a model. Only a small number of genes are currently known to be related to vitamin E metabolism in muscle or fat. In the last few years, several studies have explored the transcriptomic adaptations of skeletal muscle in response to different nutrition variables in ruminants [16–19]. Data from these studies have allowed scientists to identify the biochemical mechanisms that could be associated with key physiological processes in animals and to define specific markers for evaluating meat quality.

A better understanding of the genes and metabolic pathways associated with vitamin E metabolism is critical for identifying the key physiological processes associated with vitamin E content, oxidative stress and other metabolic pathways associated with meat quality. The main effect of vitamin E on meat is to change the lipid oxidation and thus to increase shelf-life of meat. Because lipid oxidation depends mainly on enzyme activity, post-transcriptional and post-translational changes are processes implicated in the shelf-life of meat. However, the identification of genes related to vitamin E content would be a good starting point for candidate gene selection and single nucleotide polymorphism (SNP) association studies with variations in vitamin E content and therefore implicated in meat quality. Furthermore, non-antioxidant roles of vitamin E, serving as a regulator of gene/protein, have been described. In this sense, transcriptomic analysis of skeletal muscle could identify metabolic pathways modified by vitamin E that could influence meat traits.

Therefore, the main objective of this study was to investigate transcriptional changes in the L. Thoracis muscle (LT) and subcutaneous fat (SF) of lambs supplemented with vitamin E using the Affymetrix Ovine Gene 1.1 ST whole-genome array. Furthermore, we aimed to identify novel genes that could play important roles in the metabolism of vitamin E and that might be associated with meat quality traits.

Results

Alpha-tocopherol muscle content, intramuscular fat, TBARS and metmyoglobin formation

Significant differences in weaning weight and slaughter age (SA), and average daily gain (ADG), from birth to weaning and from birth to slaughter, were found between treatments (Table 1). Animals from the CON group were younger at slaughter (P < 0.05). In addition, these animals had greater ADG from birth to weaning (P < 0.01) and from birth to slaughter (P < 0.05). During the experimental period, the ADG from weaning to slaughter was not significant between treatments. Average daily gain from birth to slaughter and SA were highly correlated between them (r 2 = -0.90; P < 0.01). The results of α-tocopherol content, intramuscular fat (IMF) content, thiobarbituric acid-reactive substances (TBARS) and metmyoglobin (MMb) formation in the LT muscle for each treatment are shown in Table 1. Our results showed that the content of α-tocopherol in the muscle was significantly higher in lambs that received a basal concentrate supplemented with dl-α-tocopheryl acetate (VE), while lambs fed a basal concentrate (CON) showed higher values of TBARS (mg of malonaldehyde per kg of L. Thoracis) (P < 0.05). Furthermore, metmyoglobin formation was significantly lower in VE lambs, compared to the CON group (P < 0.05). The K/S572/525 ratio decreased when the metmyoglobin content increased. Intramuscular fat (IMF) content was not different between the treatments (P > 0.05).

Table 1 Effect of the dietary treatment on slaughter age and weight, growth rate, alpha-tocopherol muscle content, intramuscular fat, TBARS and metmyoglobin formation
Full size table

Microarray gene expression results

Identification and classification of differentially expressed genes by microarray analysis in LT and SF

Significance analysis of microarray (SAM) was performed to compare the VE treatment with CON. SAM identified a total of 29 genes with a false discovery rate (FDR) = 0.008, 26 down- regulated and 3 up-regulated genes (Table 2). The results of SAM regarding LT muscle are shown in Fig. 1.

Table 2 The significant features identified by SAM in VE vs. CON contrast in L. Thoracis muscle
Full size table
Fig. 1
figure 1

Significant features identified by SAM in a VE-CON contrast in LT, and b VE-CON contrast in SF. The green circles represent features that exceed the specified threshold

Full size image

Regarding subcutaneous fat, when VE treatment was compared with the CON group, SAM identified a total of 330 genes with a FDR < 0.001. Among these genes, 295 were up-regulated, and 35 were down-regulated. The results of the top 50 genes identified with SAM for SF are shown in Table 3. In Additional file 1: Table S1 all of the significant genes in SF are ranked according to their fold change (FC). Notably, H1F0 gene was found to be significantly down-regulated in VE lambs in both tissues.

Table 3 Top 50 genes identify with SAM in VE vs. CON contrast in subcutaneous fat
Full size table

Treatment-dependent multivariate analysis results of gene expression

In the LT muscle, principal component analysis (PCA) of the complete set of 32 genes identified by SAM showed that the first 2 PCs covered 39.7% of the observed variance in the sample set (Fig. 2a). The PCA score plot revealed differences corresponding to lambs fed with the two different treatments. The ellipse corresponding to CON was clearly separated from the VE treatment. Partial least squares-discriminate analysis (PLS-DA) showed a clear separation of the two groups (Fig. 2b). In addition, PLS-DA allowed for the identification of the genes that were most important for the separation observed in the score plots. DEF8 gene showed the highest score, followed by ASPN and AKR7A2 (Fig. 2c). Moreover, we investigated trends or patterns in gene expression changes (Fig. 2d). For example, ABCC4, DEPTOR, IGFR1, MYLK2 and ACAT1 were positively correlated with each other in the two treatments, showing a down- and up-regulation in the VE and CON treatments, respectively. In contrast, they were negatively correlated with SAT1, ASPN, or LRRTM2. These genes, which are either positively or negatively correlated with each other, appeared to play an important role in light of the PCA and PLS-DA cluster analysis.

Fig. 2
figure 2

Multivariate analysis based on gene expression profile data in LT muscle. a Principal component analysis (PCA) score plots distinguishing between the LT muscle of lambs fed concentrate (triangle), and supplemented with vitamin E in the treatment, (+). b Partial least squares-discriminant analysis (PLC-DA) based on gene expression profile data. c Important features identified by PLS-DA. The top 15 genes are ranked by VIP scores. The colored boxes on the right indicate the relative expression of the corresponding gene in each group under study. d A bar graph showing the top 25 genes correlating with diet (CON, VE). Variables with the same distance from 0 with similar positions are positively correlated. Those with the opposite direction are negatively correlated. The light blue bars indicate genes showing a negative correlation, and the light pink bars indicate those with a positive correlation with the given pattern of change of treatment

Full size image

Multivariate analysis results from subcutaneous fat were used to cluster the samples based on gene expression profiles of animals fed two different diets: CON and VE. The PCA of the complete set of 330 genes identified by SAM showed that the first 2 PCs covered 68% of the observed variance of the sample set (Fig. 3a). The figure shows the score plot of the two first principal components extracted in this study. PLS-DA showed that the VE and CON groups were clearly separate from each other. (Fig. 3b). In addition, PLS-DA allowed for the identification of the most important genes contributing to the separation observed in the score plots, and ALG11 had the highest score, followed by SRPBR and PPX47.

Fig. 3
figure 3

Multivariate analysis based on gene expression profile data in subcutaneous fat. a Principal component analysis (PCA) score plots distinguishing between the subcutaneous fat of lambs fed concentrate (triangle), and supplemented with vitamin E in the treatment (+). b Partial least squares-discriminant analysis (PLC-DA) based on gene expression profile data. c Important features identified by PLS-DA. The top 15 genes are ranked by VIP scores. The colored boxes on the right indicate the relative expression of the corresponding gene in each group under study. d A bar graph showing the top 25 genes correlating with treatment (CON, VE). Variables with the same distance from 0 with similar positions are positively correlated. Those with the opposite direction are negatively correlated The light pink bars indicate genes with a positive correlation with the given pattern of change of treatment

Full size image

Moreover, we investigated patterns in gene expression changes in subcutaneous fat (Fig. 3d). For example, the gene expressions of ALG11, SRPBR, and PPX47 were positively correlated with each other, being up-regulated in VE treatment.

Hierarchical clustering analysis (HCA)

HCA was performed using the significant genes obtained by SAM for both contrasts. The results of HCA for LT muscle are presented in Fig. 4a. The expression profile of these genes was able to cluster and to classify correctly the samples within their corresponding groups. The heatmap shows the presence of 2 different clusters containing different genes. The responses of each variable to the two different treatments are indicated with changes in the color intensity on the heatmap. The VE and CON groups showed very different gene expression patterns. For instance, LRRTM2, ASPN and SAT1 were up-regulated in the VE group. Furthermore, a second cluster including the remainder of the genes was found to be down-regulated in the VE group. These two clusters distinguished the VE group from the CON group. These genes are involved in different metabolic processes.

Fig. 4
figure 4

Hierarchical clustering analysis of gene expression in a) the L. Thoracis muscle, b) subcutaneous fat tissue of lambs receiving different treatments (CON and VE), using the most significant genes of each contrast. Cells are colored based on the signal intensity measured. Dark brown represents high gene expression levels, blue indicates low signal intensity, and gray indicates the intermediate level (see color scale above the heatmap)

Full size image

Regarding subcutaneous fat, the results of clustering analysis of top 50 genes are presented in Fig 4b. Analyzing the heatmap, it can be observed that the CON group showed very different gene expression patterns from the VE group. There is a general down-regulation of gene expression in the CON group and an up-regulation of gene expression in the VE group. The results of clustering analysis of the total 330 genes in subcutaneous fat were similar to the 50 genes analysis results, but one VE animal was clustered in the CON group (Additional file 2: Figure S1). More specifically, 34 genes were down- regulated in the VE group (SOD3, CLEC3B, METTL7A, IER3, PLK2, CLEC1A, HSPA1A, TNFSF10, MERTK, MLLT3, TPPP2, HMGB2, PECAM1, IFITM3, H1F0, HMG20B, GALK1, NFIL3, CDH5, ERG, CEBPA, ANKRD44, HCST, UACA, IFITM1, RPL10, PID1, LUC7L3, H3F3A, AMOTL2, PSIP1, RPL17, TRA2A, and RPS259).

Functional clustering annotation

In the LT muscle, the results of Database for Annotation, Visualization and Integrated Discovery (DAVID) functional annotation clustering (FAC), including the 29 significant genes in the VE-CON contrast, showed that 4 genes from the “intracellular signaling cascade” (IGF1R, DEF8, AKAP7 and CISH) had the most enrichment and were all down-regulated in the VE treatment (Additional file 3: Table S2). However, the confident enrichment scores were less than 1.3 in both cases.

In subcutaneous fat, the results of DAVID FAC of 330 significant genes in the VE-CON contrast, revealed that the most significantly enriched cluster was “lipid biosynthesis process”, (enrichment score of 5.38) with a total of 20 genes involved (PGAP3, EBP, CRLS1, MVD, CYP51A1, GNE, HMGCS1, DPAGT1, LSS, SIGMAR1, LPCAT3, FDFT1, DOLK, PIGM, SQLE, DHCR7, AGPAT9, LTA4H, PCYT2 and HSD17B7), all of them up-regulated in the VE group compared to the CON group. Many of these genes play roles in the cholesterol, sterol and steroid biosynthesis. In Table 4, the DAVID FAC results of the 330 genes in the VE-CON contrast are shown.

Table 4 DAVID Functional Annotation Clustering of SAM genes in VE vs. CON subcutaneous fat
Full size table

Validation of microarrays results using qPCR

Thirteen genes were selected to validate the microarray results by quantitative real-time PCR (qPCR). The gene set included 4 genes in muscle and 9 genes in subcutaneous fat. In LT muscle, the genes were selected because they were significant differentially expressed between the VE and CON treatments (CISH, ABCC4, ACAT1 and IGF1R). In SF, the significant genes in the VE-CON contrast (LDLR, SOD3, SQLE, SREBF1, MTTL1, HAX1, HMGB2, HSPB8 and IER3) were selected. Although the magnitude of FC obtained by microarray and qPCR was slightly different in some instances, the qPCR results demonstrated a similar trend compared with the microarray results of these 13 genes, demonstrating the reliability of microarray analysis (Pearson’s correlation coefficient 0.99, P = 0.008 for LT muscle and 0.99, P < 0.001 for SF) (Table 5).

Table 5 Real-time PCR confirmation of the microarray results
Full size table

Discussion

The aim of the present study was to assess the effects of adding VE to the fattening concentrate, fed between weaning and slaughter, on the transcriptional changes in the LT and SF. The experimental period ended when the lambs reached the target slaughter LW (22–24 kg), according to the specifications of Ternasco de Aragón Protected Geographical Indication (Regulation (EC) No. 1107/96). The light lamb production is based on two periods: suckling and fattening. The suckling period is usually limited by the time, usually lasting approximately 45 d (in the present study, it was 48 d), while the fattening period is limited by the weight (22–24 kg LW). Thus, lambs with high ADG during lactation period spent fewer days in the fattening period than those with low ADG. To be able to evaluate the potential effects of the inclusion of VE on the transcriptome of the L. Thoracis muscle and subcutaneous fat, we planned to feed at least 30 days of VE concentrate. According to Ripoll et al. [11], supplementation with VE for 30 days was found to increase noticeably concentration of the α-tocopherol content in muscle. In the present study, the ADG during the experimental period was not different, however during the suckling period it was lower in VE treatment (p < 0.01; Table 1). This implies that the lambs of the VE treatment during the suckling period presented lower growth compared to the CON group and needed more days to reach the target weight. Moreover, there was a strong correlation between ADG from birth to slaughter and SA (r 2 =−0.9). Therefore, to overcome this unbalance between treatments ADG values were included as covariates in the statistical model used to validate gene expression differences by qPCR, thus avoiding estimation biases. We analyzed LT muscle and SF because meat cuts for human consumption include both intramuscular and subcutaneous fat [20], constituting the total amount of meat fat purchased at retail. The results related to the content of α-tocopherol in the muscle and TBARS confirm what was planned in the methodology. Vitamin E has been successfully used to increase the muscle α-tocopherol and to reduce lipid oxidation in beef and chicken meat [21]. In this study, as expected, there were significant differences between both diets in the α-tocopherol content, lipid oxidation and MMb formation in the LT muscle, as previously found by Ripoll et al. [11]. Moreover, TBARS values were lower in VE lambs which showed that the lipid oxidation process was slower in animals fed this type of diet.

In the present study, we used microarray technology to study changes in gene expression profiles in LT muscle and SF in response to vitamin E supplementation in lambs. Our results, showed that vitamin E supplementation caused different responses in gene expression in LT muscle and SF, suggesting a specific response of tissue to vitamin E supplementation. In our study, we did not measure the concentration of VE in the SF, however the gene expression results might be related to the greater α-tocopherol accumulation in adipose tissues than in skeletal muscles [22, 23].

It has been reported that α-tocopherol can influence a number of biological functions by regulating cell signaling at both the mRNA and miRNA levels [7]. Indeed, our results showed that vitamin E supplementation affected the expression of 29 genes in the LT muscle. The results of functional analysis showed that genes related to the intracellular signaling cascade (CISH, IGF1R, DEF8, and AKAP7) and metabolic processes (ZNF79, MAFB, MYLK2, ACACB, ACAT1, CISH, IGF1R, PGLS, DUSP26, AKR7A2, FBXL4, AKAP7, and RSC1A1) were down-regulated. The enrichment scores were less than 1.3, likely because the number of significant genes in this contrast was low. The most down-regulated gene by vitamin E in LT was CISH (Cytokine Inducible SH2 containing protein). Chen et al. [24] reported that CISH activates protein kinase C (PKC) activity by G-protein coupled receptor protein, which is important for the activation of both the activator protein 1 (AP-1) and nuclear factor- κB (NF-κB) pathways. In a previous study Boscoboinik et al. [25], demonstrated that PKC is inhibited by α-tocopherol. NF-κB proteins are a family of transcription factors and are of central importance in inflammation, immunity and apoptosis [26]. Evidence suggests a role for reactive oxidative intermediates (ROIs) as a common and critical intermediate for various NF-κB-activating signals, based on inhibition of NF-κB activation by a variety of antioxidants [27, 28]. There is a bidirectional relationship between cytokines and oxidative stress. Exposure of myotubes to reactive oxygen species (ROS)-producing agents resulted in an increase in interleukin 6 (IL-6) release through the activation of the redox-sensitive transcription factor, NF-κB [28]. In this sense, IL-6 induces marked increases in expression of CISH, SOCS-1, SOCS-2, and SOCS-3 in tissues, which in turn result in inhibition of the signaling of wide range of cytokines [29]. Thus, CISH expression could be related to oxidative status in the muscle. In this sense, low levels of ROS and ROIs in muscle caused by α-tocopherol could be associated with low expression of CISH. In this sense, lipid oxidation in the LT muscle and metmyoglobin formation were lower when lambs were supplemented with VE. In the same manner, RSC1A, which inhibits a dynamin and PKC-dependent exocytotic pathway of SLC5A1 gene, was down-regulated in the VE group. Interestingly enough CISH and MYLK, also down-regulated in the VE treatment, are involved in inflammation mediated by chemokine and cytokines signaling pathway, which could suggest another possible role for VE that of decreased inflammation.

In our study, alterations of the Ras homolog gene family, member A (RhoA) and actin cytoskeleton signaling were identified (IGF1R and MYLK2 genes). Both IGF1R and MYLK2 were down-regulated in VE treatment. Insulin-like growth hormone 1 (IGF-1) is a protein structurally similar to insulin, and it regulates tissue growth and development in several vertebrates [30]. As a main receptor of IGFs, IGF1R mediates the transduction of metabolic signals of cell proliferation, bone growth, and protein synthesis in the GH/IGF pathway [31]. In our study, IGF1R was down-regulated in VE treatment, and because of IGF1R polymorphisms have been associated to growth traits [32–36], we thought that this effect could be due to the higher ADG in CON animals. However, ADG values were included as a covariate in the statistical model used to validate expression differences by qPCR, thus avoiding estimation biases related to ADG. Therefore these results suggest that supplementation with VE causes down-regulation of IGF1R in LT muscle. Our findings are in agreement with Araujo et al. [37], who showed that VE supplementation reduces IGF1R expression by 17% in hyperthyroid of Wistar rats. In addition, Chuang et al. [38] also found that VE alone significantly and dose-dependently reduced the cell surface expression of IGFIR in HL-60 cells. Moreover, Holzenberger et al. [39] reported that Igf1r +/−mice display greater resistance to oxidative stress.

In addition, vitamin E down-regulated two genes related to lipid metabolisms (ACAT1 and ACACB) in the LT muscle. ACAT1 (acetyl-coenzyme A acetyltransferase 1) encodes a mitochondrial localized enzyme that catalyzes the reversible formation of acetoacetyl- CoA from two molecules of acetyl-CoA. ACAT1 is responsible for cholesterol homeostasis and maintain appropriate cholesterol availability in cell membranes, whereas ACACB is the key regulator of the fatty acid oxidation pathways [40]. It controls fatty acid oxidation by means of the ability of malonyl-CoA to inhibit carnitine-palmitoyl-CoA transferase I (CPT1B). As in our work, Shige et al. [41] showed that vitamin E reduced the uptake of modified low-density lipoprotein (LDL) and suppressed ACAT activity, resulting in less cholesterol esterification in macrophages. Interestingly enough ABCC4, an ATP-binding cassette (ABC) transporter, AKR7A2 which is involved in the detoxification of aldehydes and ketones, and finally RSC1A1 which transport carbohydrate across the plasma membrane, were all down- regulated with VE treatment. In hepatocytes, ABCC4 was shown to be induced by oxidative stress through binding of the oxidative sensor nuclear factor E2-related factor 2 (Nrf2) to antioxidant-responsive element sequences in the promoter of ABCC4 [42]. Our results showed for the first time that vitamin E down-regulates ABCC4 expression. Because of significant differences in ADG between treatments, we validated the expression results of the CISH, ABCC4, ACAT1 and IGFR1 genes by qPCR, including in the statistical model ADG as a covariate. ADG was not significant, with a similar FC between the microarray and qPCR results (r 2 = 0.99; P = 0.008). Therefore, treatment was the main effect over gene expression.

In addition, we found that the transcription factors FHL3, ZNF777 and MAFB were down-regulated. Considering all together, our results showed that supplementation with VE increased the content of α-tocopherol in the LT muscle and decreased metmyoglobin formation and lipid oxidation. We speculate that α-tocopherol in the LT muscle causes a decrease in the catalytic activity of enzymes involved in cellular transport of fatty acids and carbohydrates, and fatty acid oxidation in the mitochondria. Decreased IGF1R expression might be related to the lower lipid oxidation levels in VE animals. We also found that IGFR1, ABCC4 ACAT1, CISH, ACACB, MYLK2, ZNF777 and MAFB were positively correlated with each other and with the diet, which suggest co-expression processes of these genes. We speculate that transcription factors ZNF777, MAFB and FHL3 could be important players in mediating the effects of VE in regulating the expression of these genes.

To establish whether IGFR1, ABCC4 and ACAT1, CISH are markers of meat oxidation or indirect markers of meat quality, further studies with greater numbers of animals are necessary.

A most dramatic effect of VE was observed on SF gene expression. ALG11, SRPBR, DDX47, SEC23ID and TTC37 were among the most important genes in discriminate fed treatments. These genes were up-regulated and positively correlated with each other and with the diet, which might suggest co-expression processes. Four of the up-regulated genes in the VE group were related to heat shock proteins (HSPs) or chaperonin activity (TTC37, DNAJC16, HSB8, AHSA1). Some HSPs are characterized by various specific functions such as anti-apoptotic or anti-inflammatory effects [43]. In our study, these genes were up-regulated and positively correlated, suggesting putative increased stress protection in VE lambs.

Surprisingly VE treatment showed general up-regulation of almost all significant genes, compared to CON treatment. Lipid biosynthetic processes were among the most enriched functional clusters with major biological significance and importance (PGAP3, EBP, CRLS1, MVD, CYP51A1, GNE, HMGCS1, DPAGT1, LSS, SIGMAR1, LPCAT3, FDFT1, DOLK, PIGM, SQLE, DHCR7, AGPAT9, LTA4H, PCYT2, and HSD17B7). Moreover, genes implicated in sterol, steroid and cholesterol biosynthesis processes (SREBF1, EBP, LDLR, MVD, CYP51A1, SQLE, DHCR7, INSIG1, HMGCS1, LSS, and FDFT1) were up-regulated in the SF of VE animals, compared to the CON group. It has been previously reported that tocopherols inhibited de novo cholesterol synthesis within enterocytes [44], and cause repression of genes (DHCR7 and HMGCS1) involved in the de novo synthesis of cholesterol in testes and adrenal glands [45]. The differences reported in previous studies and ours might be due to the different tissues analyzed. This supports even more our idea that there is a tissue specific response in response to VE supplementation. On the other hand, Wang et al. [46] found oxysterol-specific repressive effects in the CYP51A1 and FDFT1 genes, mediated via direct binding of the ligands to liver X receptor (LXR). Because of α–tocopherol and other antioxidants can inhibit the oxidation of cholesterol [47], the VE group could have a lower quantity of oxysterols and then have inhibited the repression of cholesterol biosynthesis via LXR and oxysterols and up-regulated the sterol, steroid and cholesterol biosynthesis processes. In this sense, several genes implicated in “the stress response” (SOD3, IER3, HMGB2, UACA, LUC7L3) were down-regulated in the VE group, compared to the CON group (Table 3 and Supplementary Table S1). Among them, SOD3 showed higher values of FC. This gene encodes a member of the superoxide dismutase (SOD) protein family that protects the extracellular space from the toxic effects of reactive oxygen intermediates by converting superoxide radicals into hydrogen peroxide and oxygen. This finding contradicts the results of previous studies, which suggested that grass-based treatments elevated the activity of antioxidants, such as glutathione and superoxide dismutase (SOD), compared to grain-feeding [1]. However, Kumar et al. [48] in the myocardium muscle and Strobel et al. [49] in the skeletal muscle of exercise-trained and sedentary rats, found that antioxidant supplements reduced the endogenous antioxidants SOD2 gene and protein and the glutathione peroxidase (GPx) gene and enzyme activity.

Another down-regulated gene in VE treatment was IER3. IER3-deficient NCM460 cells exhibited reduced reactive oxygen species levels, indicating increased antioxidative protection [50]. In our study, IER3 was down regulated in the VE group suggesting an increase in antioxidative protection. The role of HMGB1 in recognizing aberrant or damaged DNA has been shown in multiple in vitro experiments. A recent study directly showed the accumulation of HMGB1 at sites of oxidative DNA damage in live cells, thus defining HMGB1 as a component of an early DNA damage response [51]. A similar function has been attributed to HMGB2 [52]. These authors hypothesized that HMGB1/2 proteins act as a sensor of DNA modification, and their interaction with chemically altered DNA changes the chromatin structure, thus inducing DNA damage responses.

Finally, a cluster related to tRNA metabolic processes was also significant (TRNT1, METTL1, ADAT2, NSUN2, IARS, GARS, EPRS and MARS). Aminoacyl-tRNA synthetases perform an essential function in protein synthesis by catalyzing the esterification of an amino acid to its cognate tRNA (IARS, GARS, EPRS and MARS). Considering all together, we speculate that in the SF, vitamin E exerts anti-inflammatory effect and stress protection by increasing heat shock protein expression. Vitamin E reduces stress response probably as results of reduced reactive oxygen species levels. In addition animals supplemented with VE might have inhibited cholesterol oxidation in SF and enhanced sterol, steroid and cholesterol biosynthesis processes. And finally, we speculate that regulation of these gene expressions it could be mediated through tRNAs IARS, GARS, EPRS and MARS.

Regarding LT, we included ADG values as a covariate in the statistical model used to validate the expression differences by qPCR to avoid estimation biases. In this case, ADG was not significant, and a similar FC between microarray and qPCR results was found (r 2 = 0.99; p < 0.0001). Thus, the treatment was also the main effect over gene expression of SF.

Although we found differences in mRNA activity, it did not necessarily cause differences in metabolic processes. An increase in gene expression is not necessarily correlated with an increase in protein concentrations or enzyme activities. There are many processes between transcription and translation, including post-transcriptional, translational and protein degradation regulation, in controlling steady-state protein abundances [53]. Moreover, DNA microarray has been the technology of choice for transcriptome analysis in recent years. Nonetheless, array technology has several limitations which include: using microarray technology limits the researcher to detecting transcripts that correspond to existing genomic sequencing information; and background hybridization limits the accuracy of expression measurements, particularly for transcripts present in low abundance.

Conclusion

This study demonstrated the beneficial effects of vitamin E supplementation during fattening period in lambs by increasing α-tocopherol content in the LT muscle and diminishing drastically the lipid oxidation of the meat. We observed a tissue-specific response to vitamin E supplementation. The gene expression profiles for VE and CON treatments were different in both LT and SF. Vitamin E supplementation had a dramatic effect on subcutaneous fat gene expression, showing a general up-regulation of genes, compared to CON treatment. Our study enabled us to identify novel genes (for example, IGF1R, ACAT1, ABCC4, ACACB, SOD3, and IER3) and metabolic pathways related to vitamin E metabolism that might be implicated in meat quality. To the best of our knowledge, this study was the first to report the effect of vitamin E supplementation on gene expression in the LT muscle and SF of lambs. Future exploration of these genes is necessary for a better understanding of how vitamin E affects the oxidative processes that occur in meat products.

Methods

Ethics statement

All experimental procedures including the care of animals and euthanasia were performed in accordance with the guidelines of the European Union and with Spanish regulations for the use and care of animals in research and were approved by the Animal Welfare Committee of the Centro de Investigación y Tecnología agroalimentaria (CITA) (protocol number 2009-01_MJT). In all cases, euthanasia was performed by penetrating captive bolt followed by immediate exsanguination.

Animals and sample collection

Fourteen single reared male lambs of the Rasa Aragonesa breed were weaned at 48.28 ± 0.85 days of age. Seven lambs were fed a basal concentrate from weaning to slaughter (CON treatment, 24 ± 2.62 days). The remaining 7 lambs received for 4.71 ± 2.62 days the same basal concentrate as the CON group, and thereafter until slaughter, they received a similar concentrate with the same characteristics but enriched with 500 mg dl-α-tocopheryl acetate per kg of feed (VE treatment) for 33.28 ± 1.07 days. The ingredients and chemical composition of the feedstuffs are shown in Table 6. Prior to weaning, the CON and VE lambs suckled their mothers and had free access to their concentrates. The average concentrate intake of the CON and VE lambs during the experimental period was 24.3 and 25.1 kg per lamb, respectively. The experimental procedures, management of the animals and sample details for each group are described in detail in Ripoll et al. [11].

Table 6 Ingredients and chemical composition of the feedstuffs used in the experiment
Full size table

All of the lambs were slaughtered when they attained 22–24 kg of slaughter weight (SW), according to the specifications of Ternasco de Aragón Protected Geographical Indication (Regulation (EC) No. 1107/96), which stipulates that lambs must be less than 90 days old with a SW between 22 and 24 kg. The lambs were slaughtered using EU laws in the same commercial abattoir, and the carcasses were hung by the Achilles tendon and chilled for 24 h at 4 °C in total darkness. The slaughter age, slaughter weight, and growth rate of the 2 management strategies are presented in Table 1.

Immediately after slaughter, a piece of LT muscle from the 12th thoracic vertebra and a piece of SF between the atlas and axis cervical vertebrae were cut, frozen in liquid nitrogen and stored at−80 °C until RNA isolation.

Analysis of α-tocopherol, intramuscular fat content, TBARS and metmyoglobin formation

After chilling, a piece of the LT muscle between the 4th and the 6th lumbar vertebrae was vacuum-packed and kept at−20 °C in darkness until the α-tocopherol analysis. The α-tocopherol concentration was determined by liquid extraction as described in González-Calvo et al. [54]. A portion of the loin between the 7th and the 13th thoracic vertebrae was used to measure the LT muscle oxidative processes. The color (metmyoglobin content, MMb) and lipid oxidation analysis (thiobarbituric acid-reactive substance, TBARS) were quantified after 7 days of being maintained in darkness at 4 °C. The LT muscle color and LT intramuscular fat oxidation were measured as described in González-Calvo et al. [54]. Briefly, the relative content of metmyoglobin (MMb) was estimated by the K/S572/525 ratio [55]. This ratio decreases when the MMb content increases. Intramuscular fat oxidation of the M. Longissimus thoracis was determined using the procedure reported by Pfalzgraf et al. [56]. The TBARS values are expressed as milligrams of malonaldehyde (MDA) kilogram−1 of muscle.

RNA isolation and assessment of RNA integrity

Total RNA was extracted from approximately 500 mg of LT muscle or SF using RNeasy Tissue mini kits (QIAGEN, Madrid, Spain), following the manufacturer’s protocol. Prior to microarray analysis, RNA integrity and quality were assessed by a RNA 6000 Nano LabChip on the Agilent 2100 Bioanalyzer (Agilent, Madrid, Spain) and were quantified using a nanophotometric spectrophotometer (Implen, Madrid, Spain). All RNA integrity number (RIN) values were greater than 8.

Microarray hybridization and data processing

RNA samples (n = 14, 7 samples from each treatment) were analyzed using Ovine Gene 1.1 ST Array Strip (Affymetrix, High Wycombe, UK). The Ovine gene 1.1 ST Array Strip allows for processing of four samples in parallel. This array contains a collection of 508538 probes that interrogate up to 26 unique sequences of each transcript (median probes/transcript = 23). These probes correspond to 22047 ovine genes. This study was part of a larger study in which we analyzed 84 samples corresponding to 3 treatments and 4 tissues (7 animals per treatment). In each strip, we hybridized the four tissues, selecting the animal tissue for each strip by random sampling. Microarray hybridization and scanning were performed at the Functional Genomics Core facility (Institute for Research in Biomedicine, IRB Barcelona, Spain), following the recommendations of the manufacturer. Scanned images (DAT files) were transformed into intensities (CEL files) by Affymetrix GeneChip Operating Software (GCOS). Overall array intensity was normalized between arrays to correct for systematic bias in the data and to remove the impact of non-biological influences on biological data. The imported data were analyzed at the gene-level, with exons summarized to genes, using the mean expression of all of the exons of a gene. Normalization was performed with the Robust Multi-Array Average (RMA) algorithm using quantile normalization, median polish probe summarization, and log2 probe transformation. The data sets supporting the results and discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus repository [57] and are accessible through GEO Series accession number GSE63774 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE63774? acc = GSE63774).

Validation of microarray data by qPCR

One microgram of RNA from each sample were treated with DNAse (Invitrogen, Carlsbad, CA, USA), and single-stranded cDNA was synthesized using the SuperScript®III Reverse Transcriptase kit (Invitrogen, Carlsbad, CA, USA), following the manufacturer’s recommendations. Specific exon-spanning primers for genes were generated and confirmed for specificity using BLAST (National Center for Biotechnology Information: http://www.ncbi.nlm.nih.gov/BLAST/). Before performing the qPCR reactions,conventional PCR was performed for all of the genes to test the primers and to verify the amplified products. The PCR products were sequenced to confirm gene identity using an ABI Prism 3700 (Applied Biosystems) with standard protocols. Homology searches were performed with BLAST (National Center for Biotechnology Information: https://blast.ncbi.nlm.nih.gov/Blast.cgi) to confirm the identity of the amplified fragments. Quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) was performed using SYBR Green Master Mix: SYBR Premix Ex Taq II (Tli RNase H Plus) and an ABI Prism 7500 platform (Applied Biosystem, Madrid, Spain). Standard curves for each gene were generated to calculate the amplification efficiency through 4-fold serial dilution of cDNA pooled from the LT muscle and SF. The efficiency of PCR amplification for each gene was calculated using the standard curve method (E = 10(−1/slope)). Two “connector samples” were replicated in all of the plates to remove technical variation from this source of variability. The annealing temperatures, primer concentrations, and primer sequences for CISH, ABCC4, ACAT1, IGF1R, LDLR, SQLE, SREBF1, SOD3, HAX1, HMGB2, METTL1, HSPB8, and IER3 (genes of interest: GOIs) and the reference genes (GUSB and YWHAZ) are described in Table 7. These reference genes were chosen because in previous studies they have been shown to be the most stable in these tissues [5].

Table 7 Primers forward and reverse used in RT-PCR
Full size table

Statistical analysis

Statistical analysis of performance, α-tocopherol and lipid oxidation in the LT muscle

Statistical analysis was performed using a general lineal model (GLM). Treatment (CON and VE) was included as the fixed factor for weaning age (WA), slaughter age (SA), slaughter weight, intramuscular fat content (IMF) and average daily gain (ADG). For muscle α-tocopherol, myoglobin and TBARS, the model also included intramuscular fat content and slaughter age as covariates. The results are expressed as least square means ± the standard error (SE) values and the differences were tested at a level of significance of 0.05 with the t statistic.

Statistical analysis for the identification of differentially expressed genes by microarray analysis in LT and SF

Normalized data were further analyzed using Babelomics (http://babelomics.bioinfo.cipf.es/) and MetaboAnalyst software [58]. Genes showing a statistically significant value of the Limma-test; (P < 0.01) were screened out as differentially expressed between the treatments. Significant genes were annotated based on similarity scores in blastn comparisons of Affymetrix Transcript cluster sequences against ovine sequences in GenBank. In addition, a second method, significance analysis of microarray (SAM), was used to identify and reconfirm differentially expressed genes in VE-CON contrasts. Furthermore, SAM was used to detect false positive significant genes from Limma-testing. SAM is a well-established statistical method for the identification of differentially expressed genes in microarray data analysis [58]. SAM is designed to address the FDR when running multiple tests on high-dimensional microarray data. First, it assigns a significance score to each variable based on its relative change from the standard deviation of repeated measurements. Then, it chooses variables with scores greater than an adjustable threshold Δ and compares their relative differences from the distribution estimated by random permutations of the class labels. For each Δ, a certain proportion of the variables in the permutation set will be found to be significant by chance. This number is used to calculate the FDR. Underestimation of the variability inflates the value of the statistic and can result in an increased number of false positives.

Multivariate analysis of gene expression

Multivariate analysis was performed using MetaboAnalyst, according to Xia et al. [58]. Principal component analysis (PCA), partial least squares discriminate analysis (PLS-DA), and variable importance of projection (VIP) were used to cluster the samples based on the selected gene expression profiles. Principal component analysis was used to reduce the large set of variables (genes) into 2 principal components (PCA 1 and 2). PLS-DA was used to enhance the separation between the groups by summarizing the data into a few latent variables that maximized covariance between the response and the predictors. The corresponding loading plot was used to determine the genes most responsible for separation in the PLS-DA score plot. Based on the PLS-DA results, genes were plotted according to their importance in separating the dietary groups and each gene received a value called the variable importance in the projection. Variable importance in the projection values >1 suggests that the variable is significantly involved in the separation of groups [59]. Variables with the highest VIP values were the most powerful group of discriminators. We also investigated trends or patterns in gene expression changes.

Hierarchical clustering analysis (HCA)

Cluster analysis was performed using MetaboAnalyst. In (agglomerative) hierarchical cluster analysis, each sample begins as a separate cluster and the algorithm proceeds to combine them until all of the samples belong to one cluster. Two parameters must be considered when performing hierarchical clustering. The first is the similarity measure, and the other is the clustering algorithm. For distance measurements we used the Euclidian and Ward algorithm for clustering. The results are shown as a heatmap.

Statistical analysis of gene expression validated by qPCR

The corresponding mRNA levels were measured and analyzed by their quantification cycles (Cq). The statistical methods to analyze differences in the expression rate were performed following the method proposed by Steibel et al. [60]. The mixed model fitted was:

$$ {{\mathrm{y}}_{\mathrm{rigkm}}}_{=}\mathrm{T}{\mathrm{G}}_{\mathrm{gi}} + \mathrm{b}{\left(\mathrm{I}\mathrm{M}\mathrm{F}\right)}_{\mathrm{m}} + \mathrm{b}\left(\mathrm{A}\mathrm{D}\mathrm{G}\right)\mathrm{m} + {\mathrm{P}}_{\mathrm{k}} + {\mathrm{A}}_{\mathrm{m}} + {\mathrm{e}}_{\mathrm{rigkm}} $$

where, y rigkm is the C q value (transformed data taking into account E < 2) of the gth gene (GOIs and housekeeping) from the rth well in the kth plate collected from the mth animal subjected to the ith treatment (CON, VE and ALF); TGgi is the fixed interaction among the ith treatment and the gth gene; IMFm (only used in L. Thoracis muscle tissue gene expression), and ADGm are the effects of intramuscular fat, and the average daily gain of the mth animal included as covariates; Pk is the fixed effect of the kth plate; Am is the random effect of the mth animal from which samples were collected (Am ~ (0,σ2 A)); and erigkm is the random residual. Gene specific residual variance (heterogeneous residual) was fitted to the gene by treatment effect (e rigkm  ~ N (0, σ2 egi).

To test differences (diff GOI ) in the expression rates of the target genes between treatments in terms of fold changes (FCs), the approach suggested by Steibel et al. [60] was used. The significance of the diff GOI estimates was determined with the t statistic.

Functional annotation analyses

The Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7b [61], was used to determine the pathways and processes of major biological significance and importance through the Functional Annotation Cluster (FAC) tool based on the Gene Ontology (GO) annotation function. DAVID FAC analysis was performed with the gene lists obtained after SAM analysis. Medium stringency EASE score parameters were selected to indicate confident enrichment scores of functional significance and the importance of the given pathways and processes investigated. An enrichment score of 1.3 was employed as a threshold for cluster significance [61].

References

  1. Daley CA, Abbott A, Doyle PS, Nader GA, Larson S. A review of fatty acid profiles and antioxidant content in grass-fed and grain-fed beef. Nutr J. 2010;9:10.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Azzi A, Gysin R, Kempná P, Munteanu A, Negis Y, Villacorta L, Visarius T, Zingg J. Vitamin E mediates cell signaling and regulation of gene expression. Ann N Y Acad Sci. 2004;1031(1):86–95.

    Article  CAS  PubMed  Google Scholar 

  3. Rimbach G, Moehring J, Huebbe P, Lodge JK. Gene-Regulatory Activity of alpha-Tocopherol. Molecules. 2010;15(3):1746–61.

    Article  CAS  PubMed  Google Scholar 

  4. Nakamura YK, Omaye ST. Vitamin E-modulated gene expression associated with ROS generation. J Funct Foods. 2009;1(3):241–52.

    Article  CAS  Google Scholar 

  5. González-Calvo L, Joy M, Alberti C, Ripoll G, Molino F, Serrano M, Calvo J. Effect of finishing period length with α-tocopherol supplementation on the expression of vitamin E-related genes in the muscle and subcutaneous fat of light lambs. Gene. 2014;552:225–33.

    Article  PubMed  Google Scholar 

  6. González-Calvo L, Joy M, Blanco M, Dervishi E, Molino F, Sarto P, Ripoll G, Serrano M, Calvo J. Effect of vitamin E supplementation or alfalfa grazing on fatty acid composition and expression of genes related to lipid metabolism in lambs. J Anim Sci. 2015;93:3044–54.

    Article  PubMed  Google Scholar 

  7. Gaedicke S, Zhang X, Schmelzer C, Lou Y, Doering F, Frank J, Rimbach G. Vitamin E dependent microRNA regulation in rat liver. FEBS Lett. 2008;582(23):3542–6.

    Article  CAS  PubMed  Google Scholar 

  8. Tucker J, Townsend D. Alpha-tocopherol: roles in prevention and therapy of human disease. Biomed Pharmacother. 2005;59(7):380–7.

    Article  CAS  PubMed  Google Scholar 

  9. Jose C, Pethick D, Gardner G, Jacob R. The colour stability of aged lamb benefits from Vitamin E supplementation. Proceedings 54th international congress of meat science and technology. 2008;2008:1–3.

    Google Scholar 

  10. Jensen C, Lauridsen C, Bertelsen G. Dietary vitamin E: quality and storage stability of pork and poultry. Trends Food Sci Tech. 1998;9(2):62–72.

    Article  CAS  Google Scholar 

  11. Ripoll G, González-Calvo L, Molino F, Calvo JH, Joy M. Effects of finishing period length with vitamin E supplementation and alfalfa grazing on carcass color and the evolution of meat color and the lipid oxidation of light lambs. Meat Sci. 2013;93(4):906–13.

    Article  CAS  PubMed  Google Scholar 

  12. Morrissey PA, Buckley DJ, Sheehy P, Monahan F. Vitamin E and meat quality. Proc Nutr Soc. 1994;53(02):289–95.

    Article  CAS  PubMed  Google Scholar 

  13. Barella L, Muller PY, Schlachter M, Hunziker W, Stocklin E, Spitzer V, Meier N, De Pascual-Teresa S, Minihane A-M, Rimbach G. Identification of hepatic molecular mechanisms of action of alpha-tocopherol using global gene expression profile analysis in rats. BBA Mol Basis Dis. 2004;1689(1):66–74.

    Article  CAS  Google Scholar 

  14. Rota C, Rimbach G, Minihane A-M, Stoecklin E, Barella L. Dietary vitamin E modulates differential gene expression in the rat hippocampus: potential implications for its neuroprotective properties. Nutr Neurosci. 2005;8(1):21–9.

    Article  CAS  PubMed  Google Scholar 

  15. Nier B, Weinberg PD, Rimbach G, Stöcklin E, Barella L. Differential gene expression in skeletal muscle of rats with vitamin E deficiency. IUBMB life. 2006;58(9):540–8.

    Article  CAS  PubMed  Google Scholar 

  16. Dervishi E, Joy M, Alvarez-Rodriguez J, Serrano M, Calvo J. The forage type (grazing versus hay pasture) fed to ewes and the lamb sex affect fatty acid profile and lipogenic gene expression in the longissimus muscle of suckling lambs. J Anim Sci. 2012;90(1):54–66.

    Article  CAS  PubMed  Google Scholar 

  17. Hebbachi AM, Knight BL, Wiggins D, Patel DD, Gibbons GF. Peroxisome Proliferator-activated Receptor alpha deficiency abolishes the response of lipogenic gene expression to re-feeding resoration of the normal response by activation of liver X receptor alpha. J Biol Chem. 2008;283(8):4866–76.

    Article  CAS  PubMed  Google Scholar 

  18. Lehnert S, Byrne K, Reverter A, Nattrass G, Greenwood P, Wang Y, Hudson N, Harper G. Gene expression profiling of bovine skeletal muscle in response to and during recovery from chronic and severe undernutrition. J Anim Sci. 2006;84(12):3239–50.

    Article  CAS  PubMed  Google Scholar 

  19. Reverter A, Byrne K, Bruce H, Wang Y, Dalrymple B, Lehnert S. A mixture model-based cluster analysis of DNA microarray gene expression data on Brahman and Brahman composite steers fed high-, medium-, and low-quality diets. J Anim Sci. 2003;81(8):1900–10.

    Article  CAS  PubMed  Google Scholar 

  20. Culioli J, Berri C, Mourot J. Muscle foods: consumption, composition and quality. Sci Aliment. 2003;23(1):13–34.

    Article  CAS  Google Scholar 

  21. Bou R, Grimpa S, Baucells MD, Codony R, Guardiola F. Dose and duration effect of alpha-tocopheryl acetate supplementation on chicken meat fatty acid composition, tocopherol content, and oxidative status. J Agric Food Chem. 2006;54(14):5020–6.

    Article  CAS  PubMed  Google Scholar 

  22. Lee JH, Waller JC, Melton SL, Saxton AM. Incorporation of alpha-tocopherol and linoleic acid in fresh lambs by feeding chemically treated dietary supplements containing dl-alpha-tocopheryl acetate and sunflower oil. J Agric Food Chem. 2006;54:568–73.

    Article  CAS  PubMed  Google Scholar 

  23. Liu Q, Lanari M, Schaefer D. A review of dietary vitamin e supplementation for improvement of beef quality. J Anim Sci. 1995;73:3131–40.

    Article  CAS  PubMed  Google Scholar 

  24. Chen S, Anderson PO, Li L, Sjögren H-O, Wang P, Li S-L. Functional association of cytokine-induced SH2 protein and protein kinase C in activated T cells. Int Imm. 2003;15(3):403–9.

    Article  CAS  Google Scholar 

  25. Boscoboinik D, Szewczyk A, Hensey C, Azzi A. Inhibition of cell proliferation by alpha-tocopherol: Role of protein kinase C. J Biol Chem. 1991;266(10):6188–94.

    CAS  PubMed  Google Scholar 

  26. Hayden MS, Ghosh S. Shared principles in NF-kappaB signaling. Cell. 2008;132:344–62.

    Article  CAS  PubMed  Google Scholar 

  27. Schreck R, Albermann K, Baeuerle PA. Nuclear factor kappa B: an oxidative stress-responsive transcription factor of eukaryotic cells (a review). Free Rad Res Commun. 1992;17:221–37.

    Article  CAS  Google Scholar 

  28. Kosmidou I, Vassilakopoulos T, Xagorari A, Zakynthinos S, Papapetropoulos A, Roussos C. Production of interleukin-6 by skeletal myotubes: Role of reactive oxygen species. Am J Respir Cell Mol Biol. 2002;26:587–93.

    Article  CAS  PubMed  Google Scholar 

  29. Starr R, Willson TA, Viney EM, Murray LJ, Rayner JR, Jenkins BJ, Gonda TJ, Alexander WS, Metcalf D, Nicola NA, Hilton DJ. A family of cytokine-inducible inhibitors of signalling. Nature. 1997;387:917–21.

    Article  CAS  PubMed  Google Scholar 

  30. Jones JI, Clemmons DR. Insulin-Like Growth Factors and Their Binding Proteins: Biological Actions. Endocr Rev. 1995;16(1):3–34.

    CAS  PubMed  Google Scholar 

  31. Delafontaine P, Song Y-H, Li Y. Expression, regulation, and function of IGF-1, IGF-1R, and IGF-1 binding proteins in blood vessels. Arterioscler Thromb Vasc Biol. 2004;24(3):435–44.

    Article  CAS  PubMed  Google Scholar 

  32. Lei M, Peng X, Zhou M, Luo C, Nie Q, Zhang X. Polymorphisms of the IGF1R gene and their genetic effects on chicken early growth and carcass traits. BMC Genet. 2008;9(1):70.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Curi R, De Oliveira H, Silveira A, Lopes C. Association between IGF-I, IGF-IR and GHRH gene polymorphisms and growth and carcass traits in beef cattle. Livest Prod Sci. 2005;94(3):159–67.

    Article  Google Scholar 

  34. Zhang R, Li X. Association between IGF1R, m-calpain and UCP-3 gene polymorphisms and growth traits in nanyang cattle. Mol Biol Rep. 2011;38:2179–84.

    Article  CAS  PubMed  Google Scholar 

  35. Szewczuk M, Zych S, Czerniawska-Piatkowska E, Wajcik J. Association between IGF1R/i16/TaqI and IGF1/SnaBI polymorphisms and milk production traits in Polish Holstein-Friesian cows. Anim Sci Pap Rep. 2012;30(1):13–24.

    CAS  Google Scholar 

  36. El-Magd M, Abbas H, El-kattawy A, Mokhbatly A. Novel polymorphisms of the IGF1R gene and their association with average daily gain in Egyptian buffalo (Bubalus bubalis). Domest Anim Endocrinol. 2013;45(2):105–10.

    Article  CAS  PubMed  Google Scholar 

  37. Araujo AS, Enzveiler AT, Schenkel P, Fernandes TR, Ribeiro MF, Partata WA, Llesuy S, Belló-Klein A. Oxidative stress activates insulin-like growth factor I receptor protein expression, mediating cardiac hypertrophy induced by thyroxine Mol. Cell Biochem. 2007;303:89–95.

    Article  CAS  Google Scholar 

  38. Chuang CH, Huang CS, Hu ML. Vitamin E and rutin synergistically inhibit expression of vascular endothelial growth factor through down-regulation of binding activity of activator protein-1 in human promyelocytic leukemia (HL-60) cells. Chem-Bioll Interact. 2010;183(3):434–41.

    Article  CAS  Google Scholar 

  39. Holzenberger M, Dupont J, Ducos B, Leneuve P, Geloen A, Even PC, Cervera P, Le Bouc Y. IGF-1 receptor regulates lifespan and resistance to oxidative stress in mice. Nature. 2002;421(6919):182–7.

    Article  PubMed  Google Scholar 

  40. Abu-Elheiga L, Matzuk MM, Abo-Hashema KA, Wakil SJ. Continuous fatty acid oxidation and reduced fat storage in mice lacking acetyl-CoA carboxylase 2. Science. 2001;291(5513):2613–6.

    Article  CAS  PubMed  Google Scholar 

  41. Shige H, Ishikawa T, Suzukawa M, Nishiwaki M, Yamashita T, Nakajima K, Ito T, Higashi K, Ayaori M, Yonemura A, Nestel P, Nakamura H. Vitamin E reduces cholesterol esterification and uptake of acetylated low density lipoprotein in macrophages. Lipids. 1998;33(12):1169–75.

    Article  CAS  PubMed  Google Scholar 

  42. Maher JM, Dieter MZ, Aleksunes LM, Slitt AL, Guo G, Tanaka Y, Scheffer GL, Chan JY, Manautou JE, Chen Y. Oxidative and electrophilic stress induces multidrug resistance-associated protein transporters via the nuclear factor-2 related factor-2 transcriptional pathway. Hepatology. 2007;46(5):1597–610.

    Article  CAS  PubMed  Google Scholar 

  43. Vafiadaki E, Lam C, Zhao W, Arvanitis D, Sanoudou D, Kontrogianni-Konstantopoulos A, Kranias E. HAX-1: a mitochondrial anti-apoptotic protein with emerging roles in cardiac muscle. Heart. 2011;97(24):e8.

    Article  Google Scholar 

  44. Landrier J-F, Gouranton E, Reboul E, Cardinault N, El Yazidi C, Malezet-Desmoulins C, André M, Nowicki M, Souidi M, Borel P. Vitamin E decreases endogenous cholesterol synthesis and apo-AI-mediated cholesterol secretion in Caco-2 cells. J Nutr Biochem. 2010;21(12):1207–13.

    Article  CAS  PubMed  Google Scholar 

  45. Barella L, Rora C, Stöcklin E, Rimbach G. Alpha-Tocopherol Affects Androgen Metabolism in Male Rat. Ann N Y Acad Sci. 2004;1031(1):334–6.

    Article  CAS  PubMed  Google Scholar 

  46. Wang Y, Rogers PM, Su C, Varga G, Stayrook KR, Burris TP. Regulation of cholesterologenesis by the oxysterol receptor, LXRalpha. J Biol Chem. 2008;283(39):26332–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Csallany AS, Hee-Lee J, Shoeman DW. Protection of superoxide-induced cholesterol oxidation by antioxidants in protic conditions. Int J Food Sci Nut. 2002;53:403–9.

    Article  CAS  Google Scholar 

  48. Kumar CT, Reddy VK, Prasad M, Thyagaraju K, Reddanna P. Dietary supplementation of vitamin E protects heart tissue from exercise-induced oxidant stress. Mol Cell Biochem. 1992;111(1-2):109–15.

    Article  CAS  PubMed  Google Scholar 

  49. Strobel NA, Peake JM, Matsumoto A, Marsh SA, Coombes JS, Wadley GD. Antioxidant supplementation reduces skeletal muscle mitochondrial biogenesis. Med Sci Sports Exerc. 2010;43(6):1017–24.

    Article  Google Scholar 

  50. Stachel I, Geismann C, Aden K, Deisinger F, Rosenstiel P, Schreiber S, Sebens S, Arlt A, Schafer H. Modulation of nuclear factor E2-related factor-2 (Nrf2) activation by the stress response gene immediate early response-3 (IER3) in colonic epithelial cells. J Biol Chem. 2014;289(4):1917–29.

    Article  CAS  PubMed  Google Scholar 

  51. Prasad R, Liu Y, Deterding LJ, Poltoratsky VP, Kedar PS, Horton JK, Kanno S-I, Asagoshi K, Hou EW, Khodyreva SN. HMGB1 is a cofactor in mammalian base excision repair. Mol Cell. 2007;27(5):829–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Krynetskaia NF, Phadke MS, Jadhav SH, Krynetskiy EY. Chromatin-associated proteins HMGB1/2 and PDIA3 trigger cellular response to chemotherapy-induced DNA damage. Mol Cancer Ther. 2009;8(4):864–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Vogel C, Marcotte EM. Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nat Rev Genet. 2012;13(4):227–32.

    CAS  PubMed  PubMed Central  Google Scholar 

  54. González-Calvo L, Ripoll G, Molino F, Calvo JH, Joy M. The relationship between muscle α-tocopherol concentration and meat oxidation in light lambs fed vitamin E supplements prior to slaughter. J Sci Food Agri. 2015;95:103–10.

    Article  Google Scholar 

  55. Hunt M. Meat color measurements. Reciprocal Meat Conf Proc. 1980;33:41–6.

    Google Scholar 

  56. Pfalzgraf A, Frigg M, Steinhart H. Alpha-tocopherol contents and lipid oxidation in pork Muscle and adipose tissue during storage. J Agr Food Chem. 1995;43(5):1339–42.

    Article  CAS  Google Scholar 

  57. Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M. NCBI GEO: archive for functional genomics data sets-update. Nucleic Acids Res. 2013;41(D1):D991–5.

    Article  CAS  PubMed  Google Scholar 

  58. Xia J, Psychogios N, Young N, Wishart DS. MetaboAnalyst: a web server for metabolomic data analysis and interpretation. Nucleic Acids Res. 2009;37 suppl 2:W652–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Szeto SS, Reinke SN, Sykes BD, Lemire BD. Mutations in the Saccharomyces cerevisiae succinate dehydrogenase result in distinct metabolic phenotypes revealed through 1H NMR-based metabolic footprinting. J Proteome Res. 2010;9(12):6729–39.

    Article  CAS  PubMed  Google Scholar 

  60. Steibel JP, Poletto R, Coussens PM, Rosa GJM. A powerful and flexible linear mixed model framework for the analysis of relative quantification RT-PCR data. Genomics. 2009;94(2):146–52.

    Article  CAS  PubMed  Google Scholar 

  61. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2008;4(1):44–57.

    Article  Google Scholar 

Download references

Acknowledgments

The authors wish to thank the staff of “CITA de Aragón” for its assistance in sample collection and analysis.

Funding

This study was funded by the Ministry of Education and Science of Spain and the European Union Regional Development funds (INIA-RTA 2012-0041 and RZP2013-0001) and the Research Group Funds of the Aragón Government (A49). L. González-Calvo was supported by a doctoral grant from INIA.

Availability of data and materials

The data sets supporting the results and discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus repository and are accessible through GEO Series accession number GSE63774.

Authors’ contributions

LGC and PS performed the experiments. MJ, MS and JHC designed the research and obtained funding for this research. LGC, ED, MS, JMO and JHC wrote the paper. LGC, MS, ED, MJ, RMH, JMO and JHC analyzed the data. MJ provided the animals. JHC had primary responsibility for the final content. All of the authors contributed to the manuscript discussion. All of the authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Ethics approval and consent to participate

All of the experimental procedures, including the care of animals and euthanasia, were performed in accordance with the guidelines of the European Union and with Spanish regulations for the use and care of animals in research and were approved by the Animal Welfare Committee of the Centro de Investigación y Tecnología Agroalimentaria (CITA) (protocol number 2009-01_MJT). In all cases, euthanasia was performed by penetrating captive bolt followed by immediate exsanguination.

Author information

Authors and Affiliations

  1. Unidad de Tecnología en Producción Animal, CITA, 59059, Zaragoza, Spain

    Laura González-Calvo, Margalida Joy, Pilar Sarto & Jorge H. Calvo

  2. University of Alberta, 116 St and 85 Ave, Edmonton, AB, T6G 2R3, Canada

    Elda Dervishi

  3. IMDEA-Alimentación, 28049, Madrid, Spain

    Roberto Martin-Hernandez

  4. Departamento de Mejora Genética Animal, INIA, 28040, Madrid, Spain

    Magdalena Serrano

  5. Jean Mayer-USDA Human Nutrition Research Center on Aging at Tufts University, Boston, MA, USA

    Jose M. Ordovás

  6. ARAID, 50004, Zaragoza, Spain

    Jorge H. Calvo

Authors
  1. Laura González-Calvo
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Elda Dervishi
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Margalida Joy
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Pilar Sarto
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Roberto Martin-Hernandez
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Magdalena Serrano
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Jose M. Ordovás
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Jorge H. Calvo
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Jorge H. Calvo.

Additional files

Additional file 1: Table S1.

The significant genes identified by SAM in VE vs. CON contrast. (DOCX 62 kb)

Additional file 2: Figure S1.

Hierarchical clustering analysis in subcutaneous fat using 330 SAM genes. (DOCX 290 kb)

Additional file 3: Table S2.

DAVID Functional Annotation Clustering of SAM genes in VE vs. CON muscle L. Thoracis. Only is shown the 2 most enrichment cluster. (DOCX 16 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

González-Calvo, L., Dervishi, E., Joy, M. et al. Genome-wide expression profiling in muscle and subcutaneous fat of lambs in response to the intake of concentrate supplemented with vitamin E. BMC Genomics 18, 92 (2017). https://doi.org/10.1186/s12864-016-3405-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-016-3405-8

Keywords

  • Vitamin E
  • Genome-wide expression profiling
  • Microarray
  • Sheep
  • Musclem
  • Subcutaneous fat
  • Meat quality

BMC Genomics

ISSN: 1471-2164

Contact us