Next Article in Journal
Effects of Intravenous Antimicrobial Drugs on the Equine Fecal Microbiome
Next Article in Special Issue
A High-Energy Diet and Spirulina Supplementation during Pre-Gestation, Gestation, and Lactation do Not Affect the Reproductive and Lactational Performance of Primiparous Sows
Previous Article in Journal
Pluripotent Core in Bovine Embryos: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Maternal Feed Restriction at Different Stages of Gestation on the Proteomic Profile of the Newborn Skeletal Muscle

by
Thaís Correia Costa
1,2,
Luana Lucas Dutra
3,
Tiago Antônio de Oliveira Mendes
3,
Marta Maria dos Santos
1,2,
Renata Veroneze
1,
Mateus Pies Gionbelli
4 and
Marcio de Souza Duarte
5,*
1
Department of Animal Science, Universidade Federal de Viçosa, Viçosa 36570-000, Brazil
2
Muscle Biology and Nutrigenomics Laboratory, Universidade Federal de Viçosa, Viçosa 36570-000, Brazil
3
Department of Biochemistry and Molecular Biology, Universidade Federal de Viçosa, Viçosa 36570-000, Brazil
4
Department of Animal Science, Universidade Federal de Lavras, Lavras 37200-900, Brazil
5
Department of Animal Biosciences, University of Guelph, Guelph, ON N1G 2W1, Canada
*
Author to whom correspondence should be addressed.
Submission received: 24 February 2022 / Revised: 8 April 2022 / Accepted: 11 April 2022 / Published: 13 April 2022

Abstract

:

Simple Summary

Due to forage seasonality, pregnant ruminants raised in pasture are commonly subjected to a period of feed restriction, which may lead to changes in the skeletal muscle metabolism of the offspring. Such alteration may have long-term consequences for the offspring’s performance and carcass composition. In the current study, we evaluated the effects of feed restriction at different stages of gestation on the proteomic profile in the skeletal muscle of the offspring. The results showed that both periods of restriction, whether applied in the first or last half of gestation, influence muscle energy metabolism of the progeny. Specifically, the restriction in the first half of gestation influences muscle glycogen utilization, fatty acid oxidation, and the energy-investment phase of glycolysis, while the restriction in the second half of gestation influences the energy-generation phase of glycolysis and the production of glutamine in the skeletal muscle of the offspring.

Abstract

We aimed to investigate the effects of the maternal plane of nutrition during gestation on the proteome profile of the skeletal muscle of the newborn. Pregnant goats were assigned to the following experimental treatments: restriction maintenance (RM) where pregnant dams were fed at 50% of their maintenance requirements from 8–84 days of gestation, and then feed of 100% of the maintenance requirements was supplied from 85—parturition (n = 6); maintenance restriction (MR) where pregnant dams were fed at 100% of their maintenance requirements from 8–84 days of gestation, and then experienced feed restriction of 50% of the maintenance requirements from 85—parturition (n = 8). At birth, newborns were euthanized and samples of the Longissimus dorsi muscle were collected and used to perform HPLC-MS/MS analysis. The network analyses were performed to identify the biological processes and KEGG pathways of the proteins identified as differentially abundant protein and were deemed significant when the adjusted p-value (FDR) < 0.05. Our results suggest that treatment RM affects the energy metabolism of newborns’ skeletal muscle by changing the energy-investment phase of glycolysis, in addition to utilizing glycogen as a carbon source. Moreover, the RM plane of nutrition may contribute to fatty acid oxidation and increases in the cytosolic α-KG and mitochondrial NADH levels in the skeletal muscle of the newborn. On the other hand, treatment MR likely affects the energy-generation phase of glycolysis, contributing to the accumulation of mitochondrial α-KG and the biosynthesis of glutamine.

1. Introduction

Improvement in the efficiency of livestock production has been largely demanded due to population growth, in addition to cost reduction, food safety, and quality. Among the cost reduction approaches, feed restriction is commonly adopted postnatally as a strategy of compensatory growth, which consists in the application of a period of undernutrition, followed by realimentation, to promote the increase in the animal’s growth rate and directly reduce costs of animal feed [1]. However, feed restriction during certain periods of gestation may negatively impact the embryo and fetus development, leading to long-term consequences that may impact postnatal growth and health [2,3]. Specifically, the early gestational period in utero development is markedly characterized by fetal organ development, as well as the beginning of skeletal muscle development where primary and secondary myogenesis occurs, while the majority of intramuscular adipogenesis and fibrogenesis occur in the last half of gestation [4]. Moreover, it is well-established that muscle hyperplasia and most of the fibro-adipogenic progenitor cells are limited to the prenatal period [4]. Hence, depending on the period in which maternal feed restriction is applied, distinct biological processes may be activated or deactivated by the action of transcripts, proteins, or metabolites, modifying the whole course of cell fate in the skeletal muscle of the progeny.
Commonly, pregnant ruminants raised in pasture are subjected to a period of feed restriction, resulting from the lack of quantity and quality of forages. Hence, under this situation, maternal metabolism prioritizes the delivery of nutrients for the formation of vital tissues instead of secondary tissues, such as skeletal muscle, which is a disadvantage in animals raised for meat production [5,6]. Studies evaluating the impacts of maternal undernutrition on offspring’s development, health, and performance have been reported, using targeted approaches [7] as well as sequencing-based approaches for profiling mRNAs (RNAseq) [3,8]. Although RNAseq provides important information about differentially expressed genes related to biological events comparing different experimental scenarios, the expression of mRNAs does not reflect its translation in a functional protein [9], as the mRNA transcription is only partially associated with the abundance of proteins [10]. Therefore, proteomic approaches have been demonstrated to be a useful tool in livestock studies that provides the identification of differentially abundant functional proteins, allowing the establishment of its interaction and enriched molecular pathways [11].
Given the importance of meat production for human consumption, the mechanisms involved in skeletal muscle development have been widely studied over the past years [12]. Muscle proteins are classified based on their function and solubility and can be divided into sarcoplasmic, myofibrillar, and structural proteins [13]. The sarcoplasmic fraction plays roles in regulating cellular metabolism, where the glycolytic enzymes comprise the majority of sarcoplasmic proteins [14]. Because the utilization of substrate for producing energy for skeletal muscle development is dependent on maternal nutritional status, we hypothesized that maternal feed restriction programs the skeletal muscle proteome of the offspring, based on alterations of sarcoplasmic proteins. Moreover, fractioning the protein extract allows the identification of proteins of interest that are often less expressed and hard to detect [15,16]. Besides the expansive potential of proteomics, studies linking the effects of maternal nutrition at different stages of gestation with alteration of the offspring’s skeletal muscle proteome are still little explored in livestock animals. Thus, in the present study, we aimed to investigate the impact of the maternal plane of nutrition where pregnant dams were feed-restricted at different stages of gestation (first or second half) on the proteomic profile of the skeletal muscle of newborns.

2. Materials and Methods

2.1. Animals and Sampling

The Ethical Committee on Animal Use of the Department of Animal Science at Universidade Federal de Viçosa, Minas Gerais, Brazil (protocol number 09/2017) approved all the procedures prior to the beginning of the experiment.
The detailed description of the experimental management practices was previously described in [17]. Briefly, 14 nulliparous pregnant dairy goats were randomly assigned to one of the following dietary treatments: treatment RM (n = 6) consisted of feed at 50% of maintenance requirement from 8–84 days of gestation followed by feed at 100% of maintenance requirement from day 85 of gestation to parturition (~150 days), while the treatment MR (n = 8) consisted of animals fed at 100% of maintenance requirement from 8–84 days of gestation and then fed at 50% of maintenance requirement from day 85 of gestation to parturition. It is worth mentioning that effects of maternal nutrient restriction compared to feed adhering to the nutrient requirements result in the impairment of the offspring’s skeletal muscle development [2,8,18,19]. Thus, the choice of the experimental design was based on the aim of evaluating what metabolic processes were likely affected by varying the maternal plane of nutrition when the restriction occurs during the gestational period.
The dams were weighed once a week and the week of gestation was taken into consideration for the adjustment of the dry matter intake, in order to meet the objective of the experimental treatments. The same diet was offered to the dams for both treatments and consisted of 111.6 g/kg of crude protein (CP) and 676 g/kg of total digestible nutrients (TDN) on a dry matter (DM) basis, which was composed of corn silage (723 g/kg DM basis), soybean meal (96 g/kg DM basis), ground corn (165 g/kg DM basis), and mineral mixture (16 g/kg DM basis). The amount of dry matter offered to the dams adhered to the nutritional requirements for pregnant dairy goats [20], meeting either 50% or 100% of their requirements according to the experimental treatments.
After birth, male newborn goats were immediately separated from the dams, stunned using a non-penetrating captive bolt pistol, and exsanguinated. In the case of twins (4 in RM and 5 in MR), the heaviest newborn goat was selected to be sampled. Skeletal muscle samples (~1 g) were collected from the Longissimus muscle and stored in liquid nitrogen for further protein extraction.

2.2. Protein Extraction

Sarcoplasmic protein fraction of Longissimus muscle was extracted from 0.1 g of tissue in 1 mL lysis buffer (20 mM Tris HCl pH 8, 5 mM EDTA, 1% 2-mercaptoethanol, and 1% protease inhibitor cocktail (Sigma-Aldrich®®, San Luis, MO, USA), homogenized using a shaft-type homogenizer (LabGEN 125, Cole-Parmer, Bunker Hill, IL, USA). The homogenate was centrifuged at 20,200× g for 20 min at 4 °C. The supernatant was collected, aliquoted, and stored at −80 °C. Sarcoplasmic protein content was estimated by Bradford Protein Assay (Bio-Rad, Hercules, CA, USA).

2.3. Protein Identification and Data Processing

Proteins were digested using a solution containing 50 mM ammonium bicarbonate and 20 ug trypsin (Promega, Madison, WI, USA) overnight at 37 °C. Tryptic peptides were then dried using a SpeedVac centrifuge (AG-22331, Eppendorf, Hamburg, Germany), resuspended with trifluoroacetic acid (TFA), and desalted using Zip-Tip, according to the manufacturer’s protocol. Samples were again dried in SpeedVac centrifuge (AG-22331, Eppendorf, Germany) and sent to Central Analítica of IQ-USP (São Paulo, Brazil) for protein identification and quantification performed in a NanoAquity high-performance liquid chromatographer (HPLC) coupled with a maXis 3G high-resolution Q-TOF mass spectrometer (Bruker Daltonics, Bremen, Germany). The raw data were processed with MaxQuant (v. 1.6.3.3) software with parameters set to default values, considering the protein amino terminal acetylation, methionine oxidation as variable modification, and the fixed modification as carbamidomethylation of cysteine. The trypsin specificity was kept as the digestion mode and the instrument selected was Bruker-QTOF, set to default, including the parameters of first (20 ppm) and main (10 ppm) search peptide tolerance. For the calculation of the label-free quantification (LFQ) protein intensity, the LFQ mode was added and at least two unique peptide ratios (min LFQ ratio count = 2) were considered. The goat reference proteome used was obtained from UniProt database (ID: UP000291000) available in (www.uniprot.org) (accessed on 27 September 2019).

2.4. Network Analyses

The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment analysis was built with String version 11.0 (string-db.org) (accessed on 1 January 2020). The interaction network of the exclusive proteins from the treatments MR and RM was obtained using the available interaction map from the closest species (Ovis aries), with the default option (medium confidence given by score of 0.4) and connections defined according to statistical evidence. KEGG pathways were considered significantly enriched based on adjusted FDR [21] using p-value < 0.05. For the differentially abundant proteins (DAPs), the network analysis was built in ClueGo 2.5.3 [22], a Cytoscape 3.7.2 application. The biological processes shared by the DAPs were identified based on a unilateral hypergeometric test also adjusted by FDR [23,24]. Only the biological processes identified with the adjusted p-value < 0.05 were considered for posterior analysis involving protein networks. In this context, the node colors represent the functional group; whereas the node size represents the term enrichment significance, being the additional inferences on edges connection based on the Kappa statistic (Kappa score = 0.4) [25,26].

2.5. Statistical Analysis

Statistical analysis was performed based on a peptide-based model through the MSqRob software (version 0.7.6) [27,28]. The peptide intensities were log2-transformed and normalized based on the robust linear regression (RLR) method. Prior to analysis, proteins identified by site, reverse, and potential contaminants were filtered and removed. The statistical model considered the fixed effects of treatments (RM and MR), the fixed effect of the number of offspring (1 or 2), and the random effect of sequence. The random effect was included to provide fixed effects estimation independently of the peptide effect. Since we aimed to investigate the treatment effects, the contrasts were made between treatments comparing RM vs. MR, with significance based on ANOVA and posterior differentially abundant proteins (DAPs) adopting p-value (FDR) < 0.05.

3. Results

3.1. Differentially Abundant Proteins

Initially, a total of 415 proteins were identified (Table S1) in the skeletal muscle of newborn goats. However, after removing the proteins only identified by site, potential contaminants, and reverse sequences, a total of 386 proteins were maintained. We found 181 and 46 exclusive proteins in the treatments RM and MR, respectively (Figure 1). For DAPs, 159 proteins corresponding to common proteins present in both treatments (intercept) were tested (Figure 1).
Controlling the p-value adjusted (FDR) < 0.05, a total of 13 proteins were differentially abundant between treatments (Figure 2, Table 1). Among them, ATP-dependent 6-phosphofructokinase (PFKM), dihydrolipoyl dehydrogenase (DLD), lipoyl-binding domain-containing protein (DLST), S-formylglutathione hydrolase (ESD), IF rod domain-containing protein (DES), triosephosphate isomerase (TPI1), and 5 uncharacterized proteins corresponding to the genes FLNC, HSPA9, SELEMBP1, and MYOM1 were more abundant in RM compared with MR (FRD < 0.05; Table 1). Calponin homology (CH) domain-containing protein (SMTNL1) and pyruvate kinase (PKM) were less abundant in RM compared to MR (FDR < 0.05; Table 1).

3.2. Protein-Protein Interaction Network and KEGG Pathways of the Exclusive Proteins

The interaction network of the exclusive proteins from treatment RM was significant (p-value (FDR) < 1 c 10−16), while in the treatment MR, the p-value was higher than RM (p-value (FDR) < 1.8 × 10−5); however, it was still significant, indicating that the proteins are at least partially biologically connected.
KEGG analysis from the exclusive proteins of treatment RM indicate pathways of interest related to tight junction, carbon metabolism, valine, leucine, and isoleucine degradation, mRNA surveillance pathway, fatty acid degradation, RNA transport, long-term potentiation, tryptophan metabolism, regulation of actin cytoskeleton, and citrate cycle (TCA cycle) (Table 2).
Although the proteins exclusive to MR treatment showed fewer interactions between them, interesting pathways were enriched, such as glycolysis/gluconeogenesis, citrate cycle (TCA cycle), biosynthesis of amino acids, carbon metabolism, metabolic pathways, HIF-1 signaling pathway, pyruvate metabolism, arginine and proline metabolism, and phagosome and 2-oxocarboxylic acid metabolism (Table 2).

3.3. Interaction Network and Functional Enrichment of Differentially Abundant Proteins

Due to the presence of uncharacterized proteins identified between the DAPs, we built a network based on gene names. Moreover, all DAPs (more or less abundant proteins) were combined into the same network. The DAPs network revealed enriched biological processes and sub-biological processes related to negative regulation of hematopoietic stem cell differentiation (HSP9A), protein succinylation (DLD, DLST), ribose phosphate metabolic process (PFKM, PKM, TFI1, DLD, DLST), generation of precursor metabolites and energy (AGL, TFL1, PKM, PFKM), formaldehyde catabolic process (ESD), pyruvate kinase activity (PKM), and various sub-processes of interest, such as, positive regulation of protein secretion (MYOM1), glycolytic process (TFL1, PFKM, PKM), and lysine metabolic process (DLST), among others (Figure 3).

4. Discussion

Through the differences previously reported for the average daily gain of maternal tissues between treatments [17], the effectiveness and the objective of the experimental treatments in causing feed restriction at different time-points of gestation were reached. As such, we were successfully able to assess the consequences of the plane of nutrition where restriction occurred only in the first or last half of gestation for the proteome profile of the newborn skeletal muscle.
Maternal feed restriction impairs the offspring’s skeletal muscle development, causing long-lasting impact throughout the animals’ productive life. The mechanisms of placenta adaptation in the face of a scarce intrauterine environment ensure fetus survival, and may cause permanent structural and functional changes by programming its metabolism through changes in transcripts, proteins, and metabolite profile [29]. Given the importance of vital tissues for fetuses (i.e., organs and viscera), the lack of nutrients during gestation prioritizes nutrient delivery for this purpose, while skeletal muscle development may be impaired. Thus, with the employment of proteomic approaches in the current study, we have identified the differentially abundant proteins (DAPs) and exclusively expressed proteins in the skeletal muscle of the offspring, resulting from maternal feed restriction during the first or last half of gestation. Specifically, the current study identified 13 DAPs involved in general biological processes mostly related to muscle energy metabolism. From the exclusive expressed proteins (181 in RM and 46 in MR), we performed a network analysis and the identification of KEGG pathways revealed enriched pathways related to fatty acid degradation, glycolysis/gluconeogenesis, citrate cycle (TCA cycle), and biosynthesis of amino acids, among others. Although the enriched biological processes and pathways were similar between treatments, the protein abundance inherent to each treatment affects different steps of the pathway and thus influences the muscle energy balance.
The skeletal muscle characteristics of plasticity allow its metabolism to change and adapt according to the environment, such as calorie and nutrient intake [30]. In terms of metabolic properties, skeletal muscle is the primary site of glucose uptake and storage [31]. Glucose is oxidized to generate ATP through two major pathways: the oxidative (aerobic) and the glycolytic (anaerobic) pathway. Through the aerobic pathway, glucose undergoes glycolysis to generate pyruvate, followed by its conversion into acetyl-coA, which enters the citrate cycle, producing substrates for the electron transport chain which results in ATP synthesis.
The proteins found to be differentially abundant and exclusive in our experimental treatments appear in biological processes (BP) and signaling pathways (SP) associated with energy metabolism, such as glycolytic process (BP), pyruvate metabolism (SP), glycolysis/gluconeogenesis (SP), and citrate cycle (TCA cycle; SP) (Figure 4).

4.1. Glycolysis

Glycolysis is commonly divided into the energy-investment and energy-generation phases [32,33,34]. Among the more abundant proteins in RM compared to MR, PFKM and TPI1 play roles in the energy-investment phase of glycolysis. The enzyme PFKM uses ATP to catalyze the irreversible conversion of fructose-6-phosphate into fructose-1,6-bisphosphate (FBP). FBP is then cleaved into glyceraldehyde 3-phosphate (GAP) and dihydroxyacetone phosphate (DHAP). DHAP is a precursor of lipid synthesis, while GAP is utilized in the next steps of glycolysis. Moreover, DHAP and GAP can be reversibly converted by the reaction involving the enzyme TPI1 [35]. A recent study showed that DHAP may function by signaling the availability of glucose to activate the complex of proteins mTORC1 [36], which is an important regulator of cellular growth, protein synthesis, autophagy, and lipogenesis [37]. Studies evaluating post-mortem glycolysis and meat quality parameters observed that TPI1 abundance was associated with beef tenderness [38,39] and was positively related to the rate of pH decline, which may partially be associated with meat color [40,41,42,43]. Therefore, since TPI1 may influence the abundance of both metabolites (DHAP and GAP), and consequently, the synthesis of lipids and ATP, the greater abundance of TPI1 in the skeletal muscle of animals born from dams that were feed-restricted in the first half of gestation compared to those born from dams feed-restricted in the last half of gestation may not be attributed to a specific pathway (lipid or ATP synthesis). However, we speculate that the greater abundance of TPI1 in the skeletal muscle of offspring from RM treatment would negatively impact the conversion from muscle into meat upon slaughter later in life, with consequences for meat color and tenderness parameters.
According to Akram [32], when glycolysis is initiated from the glycogen breakdown (glycogenolysis), fewer ATPs are consumed in the first phase of glycolysis, because in this case, the first step of glucose activation is dispensable. In this context, the protein AGL, which is an important protein in the process of glycogenolysis, was more abundant in the treatment RM compared to MR. Such an observation may indicate that, despite the ATP utilization in the reaction catalyzed by the PFKM, the skeletal muscle of the offspring from treatment RM may use glycogen as its main carbon source, and consequently, more ATP would be spared from the first step of glycolysis.
The second phase of glycolysis, named the energy-generation phase, is characterized by the oxidative conversion of GAP to pyruvate and the formation of ATP and NADH [34]. Three enzymes that participate in this phase of glycolysis were exclusively expressed in the treatment MR or were less abundant in RM compared to MR. GAPDH, which was exclusively detected in treatment MR, catalyzes the reversible conversion of GAP into 1,3-biphosphoglycerate, in addition to generating NADH from this reaction [32]. An important glycolytic enzyme, ENO2, was also exclusively expressed in the treatment MR. Although ENO2 (gamma-enolase) is one of the enolase isoforms found to be predominantly located in the neuron and neuroendocrine tissues [44], ENO2 was found to be expressed in human and rat cultured muscle cells [45]. The reversible reaction that converts 2-phosphoglycerate into the energy-rich phosphoenolpyruvate (PEP) is catalyzed by ENO2 [32,46]. After that, PEP is dephosphorylated in a reaction catalyzed by the enzyme PKM, downregulated in RM compared to MR, leading to the final glycolytic step where pyruvate and ATP are formed [46]. Thus, the results related to glycolysis indicate that the treatment RM may cause an enhancement in the energy-investment phase of glycolysis, while treatment MR may have a major function favoring the energy-generation phase of glycolysis.
Once formed, pyruvate may follow different pathways. In our current study, we identified proteins in both treatments that participate in the same destiny as pyruvate. The oxidative decarboxylation of pyruvate to form acetyl-CoA, CO2, and NADH is catalyzed by a well-orchestrated complex called pyruvate dehydrogenase complex (PDC). The PDC is composed of three catalytic subunits: pyruvate dehydrogenase (PDH), dihydrolipoamide acetyltransferase (DLAT), and dihydrolipoamide dehydrogenase (DLD) [47,48]. Intriguingly, PDHA1, which is a component of the PDH subunit, was exclusively expressed in the treatment MR. On the other hand, the proteins that are components of the subunit DLAT were exclusively expressed in treatment RM, while proteins that are components of DLD were more abundant in RM compared to the MR group. Because our experimental treatment was based on feed restriction, it is possible that some vitamins may have been differentially metabolized by the dams, and consequently, their availability to the fetuses was different between treatments. For instance, vitamin B1 is a cofactor of the subunit PDH [49], while vitamin B5 is the key precursor for the biosynthesis of CoA [49], utilized as a cofactor in subunit DLAT, and the precursors of FAD and NAD are vitamin B2 [49] and B3 [50], respectively, which serve as cofactors in the subunit DLD.

4.2. Citrate Cycle (TCA)

Citrate cycle (TCA) initiates catabolizing of acetyl-CoA molecules through the steps involving many enzymes that produce reducing equivalents NADH and FADH, which will be further driven to ATP production in the electron-transport chain [46]. The first step of TCA produces citrate, and like pyruvate, citrate may follow a variety of pathways both inside and outside the mitochondria matrix, in addition to regulating intermediate metabolites. When cellular iron levels are high, the enzyme ACO1, which was found to be exclusively expressed in treatment RM, catalyzes the reversible conversion of the cytosolic citrate into isocitrate [51,52]. Cytosolic isocitrate is then metabolized by the enzyme IDH1, also found to be exclusively expressed in treatment MR, generating α-ketoglutarate (a-KG) and NADPH [52], which is an important cofactor involved in lipid and cholesterol metabolism [53,54]. Moreover, findings from Moreno et al. [55] showed that ACO1 was positively associated with adipogenic markers, linking the iron metabolism with the adipogenic potential of the adipose tissue. Thus, our results may suggest that, compared to treatment MR, maternal feed restriction during the first half of gestation may influence the iron metabolism in the skeletal muscle of the offspring, which may result in enhancement of the protein abundance of ACO1 and IDH1 and likely increase the biosynthesis of lipids.
Treatment RM may also influence the pathways downstream a-KG, since the components of the a-KG dehydrogenase (α-KGDH) complex were found to be exclusively expressed in treatment RM (OGDH) or were more abundant in RM compared to MR (DLAT and DLST). The a-KGDH complex function catalyzed the conversion of a-ketoglutarate to succinyl-CoA and produced NADH [56]. Thus, based on these results, RM treatment, not only influences the a-KG production in cytosol but also its downstream reaction on the mitochondria, leading to the enhancement of NADH synthesis.

4.3. Glutamine

Regarding the cytosolic formation of α-KG, in contrast to our findings in treatment RM, treatment MR may have influenced the increase in mitochondrial α-KG, since the enzyme responsible for the conversion of the mitochondrial isocitrate into α-KG (IDH3A) was exclusively expressed in treatment MR. The amount of α-KG in mitochondria is dependent on the state of oxidation reduction (redox). For example, high levels of NAD over NADH lead to the conversion of α-KG into succinyl-CoA, while high levels of NADH over NAD lead to α-KG accumulation and consequently its participation in several other biological processes, such as the biosynthesis of amino acids. To prevent the excess of α-KG in the cells, reactions associated with the production of glutamine may occur [57]. The protein GLUL, exclusively expressed in the treatment MR, uses ammonia to catalyze the conversion of glutamate into glutamine [58,59]. When required, skeletal muscle releases glutamine into the blood, supplying the energy and protein requirements of other cells and tissues [60]. Moreover, glutamine oxidation is crucial for energy production and survival of pluripotent stem cells, in addition to donating nitrogen in de novo nucleotide synthesis [61]. Taken together, the exclusive expression of the proteins IDH3A and GLUL in the skeletal muscle of the offspring from treatment MR may indicate the accumulation of mitochondrial α-KG and its fate of glutamine synthesis, increasing the availability of this amino acid, which can be further utilized in several biological processes.

4.4. Fatty Acid Degradation

Under the aerobic conditions, ATP production is obtained through three key biological pathways, including the citrate cycle and fatty acid degradation. As previously discussed, our experimental treatments may have affected the overall skeletal muscle energy metabolism. The proteins ACAT1 and ECHS1 found to be exclusively expressed in treatment RM are involved in the fatty acid beta-oxidation signaling pathways. ECHS1 catalyzes the second reaction, while ACAT1 catalyzes the last reaction, ending in acetyl-CoA synthesis. Our results were contrary to the study which reported the downregulation of the corresponding genes in bulls subjected to compensatory growth conditions (restriction followed by realimentation) [62], reinforcing that the results may depend on the period of feed restriction (prenatal vs. postnatal) or post-transcriptional factors. Interestingly, ACAT1 can acetylate the lysine residues of PDHA1, causing inhibition of the PDC in some conditions and favoring cell proliferation [63]. Hence, the exclusive expression of PDHA1 in treatment MR in particular may be related to the exclusive expression of ACAT1 in RM, which is able to inhibit this protein. Taken together, these findings may result in increased fatty acid degradation in the skeletal muscle of the offspring from treatment RM.

4.5. Pentose Phosphate Pathways

Interestingly, the nucleotide metabolism may have been affected by our experimental treatment. The protein PGD was exclusively expressed in treatment RM, and it links the cytosolic carbohydrate metabolism with protein secretion [64]. PGD is a NADP-dependent enzyme that catalyzes the decarboxylation of 6-phosphogluconate to ribulose 5-phosphate (Ru5P), producing NADPH in the aerobic stage of the pentose phosphate pathways [46]. NADPH is mainly used in the processes of cellular redox balance and antioxidant defense [65], while the Ru5P is used for the synthesis of nucleotides and nucleic acids [46]. Therefore, the greater production of NADPH generated by the reaction involving PGD may be required in the skeletal muscle of the offspring, resulting from maternal feed restriction at the first half of gestation.

5. Conclusions

In conclusion, our data indicate that the maternal plane of nutrition where feed restriction occurs at different stages of gestation affects the overall proteins related to energy metabolism in the skeletal muscle of the newborn. Specifically, our results showed that maternal feed restriction during the first half of gestation followed by non-restriction regulates proteins related to the energy-investment phase of glycolysis. On the other hand, the energy-generation phase of glycolysis is steeper in the skeletal muscle of the offspring resulting from dams that were not feed-restricted at the first half of gestation and were feed-restricted at the second half of gestation. Moreover, carbon sources may be largely provided by glycogen in the skeletal muscle of the offspring born from dams feed-restricted at the first half of gestation. Indeed, skeletal muscle of newborns from dams feed-restricted at the first half of gestation had an increase in cytosolic α-KG; however, it appears that α-KG is used to produce NADH in mitochondria. On the other hand, in the newborns from dams feed-restricted at the second half of gestation, the α-KG accumulation in mitochondria appears to have deviated for the biosynthesis of glutamine. Furthermore, proteins involved in fatty acid degradation were enhanced by feed restriction at the first half of gestation, in addition to influencing the NADPH production by regulating an enzyme that participates in the pentose phosphate pathways. However, it is possible that offspring from a dam undergoing feed restriction during the first half of gestation likely experience extensive use of glycogen and fatty acid storage. Although the long-term effects of maternal nutrition on skeletal muscle development and metabolism have been exhibited in previous studies, if the changes observed in our trial persist in the skeletal muscle throughout the progeny’s life, they need to be further investigated.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ani12081011/s1, Table S1: Sheet containing all the proteins identified, its respective molecular weight, and the LFQ intensities of each replicate.

Author Contributions

Conceptualization, M.d.S.D. and M.P.G.; methodology, T.C.C., T.A.d.O.M. and L.L.D.; software, T.C.C., T.A.d.O.M. and L.L.D.; validation, T.A.d.O.M. and T.C.C.; formal analysis, T.C.C. and T.A.d.O.M.; investigation, T.C.C. and T.A.d.O.M.; resources, M.d.S.D., T.A.d.O.M. and M.P.G.; data curation, T.C.C. and T.A.d.O.M.; writing—original draft preparation, T.C.C.; writing—review and editing, T.C.C., T.A.d.O.M. and M.d.S.D.; visualization, T.C.C., L.L.D., T.A.d.O.M., M.M.d.S., R.V., M.P.G. and M.d.S.D.; supervision, M.d.S.D. and T.A.d.O.M.; project administration, T.C.C.; funding acquisition, M.d.S.D. and M.P.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by CNPq—Conselho Nacional de Desenvolvimento Científico e Tecnológico (Grant#311883/2018–4), CAPES—Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (Grant#001), FAPEMIG—Fundação de Amparo à Pesquisa do Estado de Minas Gerais (Grant#APQ-02496–17), and INCT-CA—Instituto Nacional de Ciência e Tecnologia de Ciência Animal (Grant#465377/2014-9).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and was approved by the Animal Care and Use Committee of the Department of Animal Science at the Universidade Federal de Viçosa, Viçosa, Minas Gerais, Brazil (protocol 09/2017).

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data are provided in the supplementary table. Any additional data will be made available upon reasonable request to the corresponding author.

Acknowledgments

The authors thank Fabyano Fonseca e Silva, from the Department of Animal Sciences of Federal University of Viçosa—Brazil, who contributed deeply to this study, and sadly passed away while the manuscript was being assessed. The authors dedicate this manuscript to Fonseca e Silva, a researcher and more importantly, a friend, who certainly was one of the greateststatisticians and geneticists in the Animal Science field.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hornick, J.L.; Van Eenaeme, C.; Gérard, O.; Dufrasne, I.; Istasse, L. Mechanisms of Reduced and Compensatory Growth. Domest. Anim. Endocrinol. 2000, 19, 121–132. [Google Scholar] [CrossRef]
  2. Paradis, F.; Wood, K.M.; Swanson, K.C.; Miller, S.P.; McBride, B.W.; Fitzsimmons, C. Maternal Nutrient Restriction in Mid-to-Late Gestation Influences Fetal MRNA Expression in Muscle Tissues in Beef Cattle. BMC Genom. 2017, 18, 1–14. [Google Scholar] [CrossRef] [PubMed]
  3. Gauvin, M.C.; Pillai, S.M.; Reed, S.A.; Stevens, J.R.; Hoffman, M.L.; Jones, A.K.; Zinn, S.A.; Govoni, K.E. Poor Maternal Nutrition during Gestation in Sheep Alters Prenatal Muscle Growth and Development in Offspring. J. Anim. Sci. 2020, 98, 1–15. [Google Scholar] [CrossRef] [PubMed]
  4. Du, M.; Tong, J.; Zhao, J.; Underwood, K.R.; Zhu, M.; Ford, S.P.; Nathanielsz, P.W. Fetal Programming of Skeletal Muscle Development in Ruminant Animals. J. Anim. Sci. 2010, 88, E51–E60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Long, N.M.; Vonnahme, K.A.; Hess, B.W.; Nathanielsz, P.W.; Ford, S.P. Effects of Early Gestational Undernutrition on Fetal Growth, Organ Development, and Placentomal Composition in the Bovine. J. Anim. Sci. 2009, 87, 1950–1959. [Google Scholar] [CrossRef] [Green Version]
  6. Meyer, A.M.; Reed, J.J.; Vonnahme, K.A.; Soto-Navarro, S.A.; Reynolds, L.P.; Ford, S.P.; Hess, B.W.; Caton, J.S. Effects of Stage of Gestation and Nutrient Restriction during Early to Mid-Gestation on Maternal and Fetal Visceral Organ Mass and Indices of Jejunal Growth and Vascularity in Beef Cows. J. Anim. Sci. 2010, 88, 2410–2424. [Google Scholar] [CrossRef] [Green Version]
  7. Gonzalez, J.M.; Camacho, L.E.; Ebarb, S.M.; Swanson, K.C.; Vonnahme, K.A.; Stelzleni, A.M.; Johnson, S.E. Realimentation of Nutrient Restricted Pregnant Beef Cows Supports Compensatory Fetal Muscle Growth. J. Anim. Sci. 2013, 91, 4797–4806. [Google Scholar] [CrossRef]
  8. Sanglard, L.P.; Nascimento, M.; Moriel, P.; Sommer, J.; Ashwell, M.; Poore, M.H.; Duarte, M. de S.; Serão, N.V.L. Impact of Energy Restriction during Late Gestation on the Muscle and Blood Transcriptome of Beef Calves after Preconditioning. BMC Genom. 2018, 19, 1–18. [Google Scholar] [CrossRef]
  9. Vogel, C.; Marcotte, E.M. Insights into the Regulation of Protein Abundance from Proteomic and Transcriptomic Analyses. Nat. Rev. Genet. 2012, 13, 227–232. [Google Scholar] [CrossRef]
  10. Abreu, R.S.; Penalva, L.O.; Marcotte, E.M.; Vogel, C. Global Signatures of Protein and MRNA Expression Levels. Mol. Biosyst. 2009, 5, 1512–1526. [Google Scholar]
  11. Mohammadabadi, M.; Bordbar, F.; Jensen, J.; Du, M.; Guo, W. Key Genes Regulating Skeletal Muscle Development and Growth in Farm Animals. Animals 2021, 11, 835. [Google Scholar] [CrossRef]
  12. Du, M.; Ford, S.P.; Zhu, M.-J. Optimizing Livestock Production Efficiency through Maternal Nutritional Management and Fetal Developmental Programming. Anim. Front. 2017, 7, 5–11. [Google Scholar] [CrossRef]
  13. Feiner, G. Meat Products Handbook: Practical Science and Technology; CRC Press: Boca Raton, FL, USA, 2006; ISBN 1845691725. [Google Scholar]
  14. Scopes, R.K.; Stoter, A. [79] Purification of All Glycolytic Enzymes from One Muscle Extract. Methods Enzymol. 1982, 90, 479–490. [Google Scholar] [PubMed]
  15. Jia, X.; Hollung, K.; Therkildsen, M.; Hildrum, K.I.; Bendixen, E. Proteome Analysis of Early Post-Mortem Changes in Two Bovine Muscle Types:M. Longissimus Dorsi and M. Semitendinosis. Proteomics 2006, 6, 936–944. [Google Scholar] [CrossRef] [PubMed]
  16. Sayd, T.; Morzel, M.; Chambon, C.; Franck, M.; Figwer, P.; Larzul, C.; Le Roy, P.; Monin, G.; Chérel, P.; Laville, E. Proteome Analysis of the Sarcoplasmic Fraction of Pig Semimembranosus Muscle: Implications on Meat Color Development. J. Agric. Food Chem. 2006, 54, 2732–2737. [Google Scholar] [CrossRef] [PubMed]
  17. Costa, T.C.; Moura, F.H.; Souza, R.O.; Lopes, M.M.; Fontes, M.M.S.; Serão, N.V.L.; Sanglard, L.P.; Du, M.; Gionbelli, M.P.; Duarte, M.S. Effect of Maternal Feed Restriction in Dairy Goats at Different Stages of Gestation on Skeletal Muscle Development and Energy Metabolism of Kids at the Time of Births. Anim. Reprod. Sci. 2019, 206, 46–59. [Google Scholar] [CrossRef]
  18. Du, M.; Zhu, M.J.; Means, W.J.; Hess, B.W.; Ford, S.P. Nutrient Restriction Differentially Modulates the Mammalian Target of Rapamycin Signaling and the Ubiquitin-Proteasome System in Skeletal Muscle of Cows and Their Fetuses. J. Anim. Sci. 2005, 83, 117–123. [Google Scholar] [CrossRef]
  19. Aragão, R.S.; Guzmán, O.Q.; Pérez, G.G.; Manhães, R.C.; Bolaños, F.J. Maternal Protein Restriction Impairs the Transcriptional Metabolic Flexibility of Skeletal Muscle in Adult Rat Offspring. Br. J. Nutr. 2014, 112, 328–337. [Google Scholar] [CrossRef] [Green Version]
  20. NRC. Nutrient Requirements of Small Ruminants: Sheep, Goats, Cervids, and New World Camelids; National Academy of Science: Washington, DC, USA, 2007; ISBN 978-0-309-10213-1. [Google Scholar]
  21. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  22. Bindea, G.; Mlecnik, B.; Hackl, H.; Charoentong, P.; Tosolini, M.; Kirilovsky, A.; Fridman, W.H.; Pagès, F.; Trajanoski, Z.; Galon, J. ClueGO: A Cytoscape Plug-in to Decipher Functionally Grouped Gene Ontology and Pathway Annotation Networks. Bioinformatics 2009, 25, 1091–1093. [Google Scholar] [CrossRef] [Green Version]
  23. Bei, Y.; Hong, P. A Novel Approach to Minimize False Discovery Rate in Genome-Wide Data Analysis. BMC Syst. Biol. 2013, 7, S1. [Google Scholar] [CrossRef] [Green Version]
  24. Lewin, A.M.; Grieve, I.C. Grouping Gene Ontology Terms to Improve the Assessment of Gene Set Enrichment in Microarray Data. BMC Bioinform. 2006, 7, 1–9. [Google Scholar] [CrossRef] [Green Version]
  25. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Bioinformatics Enrichment Tools: Paths toward the Comprehensive Functional Analysis of Large Gene Lists. Nucleic Acids Res. 2009, 37, 1–13. [Google Scholar] [CrossRef] [Green Version]
  26. McHugh, M.L. Interrater Reliability: The Kappa Statistic. Biochem. Med. 2012, 22, 276–282. [Google Scholar] [CrossRef]
  27. Goeminne, L.J.E.; Gevaert, K.; Clement, L. Experimental Design and Data-Analysis in Label-Free Quantitative LC/MS Proteomics: A Tutorial with MSqRob. J. Proteom. 2018, 171, 23–36. [Google Scholar] [CrossRef]
  28. Goeminne, L.J.E.; Gevaert, K.; Clement, L. Peptide-Level Robust Ridge Regression Improves Estimation, Sensitivity, and Specificity in Data-Dependent Quantitative Label-Free Shotgun Proteomics. Mol. Cell. Proteom. 2016, 15, 657–668. [Google Scholar] [CrossRef] [Green Version]
  29. Zhang, S.; Regnault, T.R.H.; Barker, P.L.; Botting, K.J.; McMillen, I.C.; McMillan, C.M.; Roberts, C.T.; Morrison, J.L. Placental Adaptations in Growth Restriction. Nutrients 2015, 7, 360–389. [Google Scholar] [CrossRef] [Green Version]
  30. Argilés, J.M.; Campos, N.; Lopez-Pedrosa, J.M.; Rueda, R.; Rodriguez-Mañas, L. Skeletal Muscle Regulates Metabolism via Interorgan Crosstalk: Roles in Health and Disease. J. Am. Med. Dir. Assoc. 2016, 17, 789–796. [Google Scholar] [CrossRef] [Green Version]
  31. Meyer, C.; Dostou, J.M.; Welle, S.L.; Gerich, J.E. Role of Human Liver, Kidney, and Skeletal Muscle in Postprandial Glucose Homeostasis. Am. J. Physiol.-Endocrinol. Metab. 2002, 282, E419–E427. [Google Scholar] [CrossRef] [Green Version]
  32. Akram, M. Mini-Review on Glycolysis and Cancer. J. Cancer Educ. 2013, 28, 454–457. [Google Scholar] [CrossRef]
  33. Ferguson, D.M.; Gerrard, D.E. Regulation of Post-Mortem Glycolysis in Ruminant Muscle. Anim. Prod. Sci. 2014, 54, 464–481. [Google Scholar] [CrossRef] [Green Version]
  34. Cox, M.M.; Nelson, D.L. Lehninger Principles of Biochemistry; Wh Freeman: New York, NY, USA, 2008; Volume 5. [Google Scholar]
  35. Roland, B.P.; Stuchul, K.A.; Larsen, S.B.; Amrich, C.G.; VanDemark, A.P.; Celotto, A.M.; Palladino, M.J. Evidence of a Triosephosphate Isomerase Non-Catalytic Function Crucial to Behavior and Longevity. J. Cell Sci. 2013, 126, 3151–3158. [Google Scholar] [PubMed] [Green Version]
  36. Orozco, J.M.; Krawczyk, P.A.; Scaria, S.M.; Cangelosi, A.L.; Chan, S.H.; Kunchok, T.; Lewis, C.A.; Sabatini, D.M. Dihydroxyacetone Phosphate Signals Glucose Availability to MTORC1. Nat. Metab. 2020, 2, 893–901. [Google Scholar] [CrossRef] [PubMed]
  37. Peterson, T.R.; Sengupta, S.S.; Harris, T.E.; Carmack, A.E.; Kang, S.A.; Balderas, E.; Guertin, D.A.; Madden, K.L.; Carpenter, A.E.; Finck, B.N.; et al. MTOR Complex 1 Regulates Lipin 1 Localization to Control the Srebp Pathway. Cell 2011, 146, 408–420. [Google Scholar] [CrossRef] [Green Version]
  38. Grabež, V.; Kathri, M.; Phung, V.; Moe, K.M.; Slinde, E.; Skaugen, M.; Saarem, K.; Egelandsdal, B. Protein Expression and Oxygen Consumption Rate of Early Postmortem Mitochondria Relate to Meat Tenderness. J. Anim. Sci. 2015, 93, 1967–1979. [Google Scholar] [CrossRef] [Green Version]
  39. Picard, B.; Gagaoua, M. Meta-Proteomics for the Discovery of Protein Biomarkers of Beef Tenderness: An Overview of Integrated Studies. Food Res. Int. 2020, 127, 108739. [Google Scholar] [CrossRef]
  40. Gagaoua, M.; Bonnet, M.; De Koning, L.; Picard, B. Reverse Phase Protein Array for the Quantification and Validation of Protein Biomarkers of Beef Qualities: The Case of Meat Color from Charolais Breed. Meat Sci. 2018, 145, 308–319. [Google Scholar] [CrossRef]
  41. Huang, H.; Larsen, M.R.; Karlsson, A.H.; Pomponio, L.; Costa, L.N.; Lametsch, R. Gel-Based Phosphoproteomics Analysis of Sarcoplasmic Proteins in Postmortem Porcine Muscle with PH Decline Rate and Time Differences. Proteomics 2011, 11, 4063–4076. [Google Scholar] [CrossRef] [Green Version]
  42. Nair, M.N.; Li, S.; Beach, C.; Rentfrow, G.; Suman, S.P. Intramuscular Variations in Color and Sarcoplasmic Proteome of Beef during Postmortem Aging. Meat Muscle Biol. 2018, 2, 92–101. [Google Scholar] [CrossRef]
  43. Wu, W.; Dai, R.T.; Bendixen, E. Comparing SRM and SWATH Methods for Quantitation of Bovine Muscle Proteomes. J. Agric. Food Chem. 2019, 67, 1608–1618. [Google Scholar] [CrossRef]
  44. Gagaoua, M.; Hughes, J.; Terlouw, E.M.C.; Warner, R.D.; Purslow, P.P.; Lorenzo, J.M.; Picard, B. Proteomic Biomarkers of Beef Colour. Trends Food Sci. Technol. 2020, 101, 234–252. [Google Scholar] [CrossRef]
  45. Merkulova, T.; Dehaupas, M.; Nevers, M.C.; Créminon, C.; Alameddine, H.; Keller, A. Differential Modulation of α, β and γ Enolase Isoforms in Regenerating Mouse Skeletal Muscle. Eur. J. Biochem. 2000, 267, 3735–3743. [Google Scholar] [CrossRef] [PubMed]
  46. Werner, C.; Doenst, T.; Schwarzer, M. Metabolic Pathways and Cycles. In The Scientist’s Guide to Cardiac Metabolism; Academic Press: San Diego, CA, USA, 2016; pp. 39–55. ISBN 9780128023945. [Google Scholar]
  47. Guevara, E.L.; Yang, L.; Birkaya, B.; Zhou, J.; Nemeria, N.S.; Patel, M.S.; Jordan, F. Global View of Cognate Kinase Activation by the Human Pyruvate Dehydrogenase Complex. Sci. Rep. 2017, 7, 1–11. [Google Scholar]
  48. Biensø, R.S.; Knudsen, J.G.; Brandt, N.; Pedersen, P.A.; Pilegaard, H. Effects of IL-6 on Pyruvate Dehydrogenase Regulation in Mouse Skeletal Muscle. Pflugers Arch. Eur. J. Physiol. 2014, 466, 1647–1657. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Shin, M.; Umezawa, C.; Shin, T. Metabolic Pathways: Metabolism of Minerals and Vitamins. In Encyclopedia of Food Microbiology: Second Edition; Academic Press: Burlington, MA, USA, 2014; pp. 535–543. ISBN 9780123847331. [Google Scholar]
  50. Poljsak, B.; Milisav, I. Vitamin B3 Forms as Precursors to NAD+: Are They Safe? Trends Food Sci. Technol. 2018, 79, 198–203. [Google Scholar] [CrossRef]
  51. Volz, K. The Functional Duality of Iron Regulatory Protein 1. Curr. Opin. Struct. Biol. 2008, 18, 106–111. [Google Scholar] [CrossRef] [Green Version]
  52. McGahan, M.C.; Harned, J.; Mukunnemkeril, M.; Goralska, M.; Fleisher, L.; Ferrell, J.B. Iron Alters Glutamate Secretion by Regulating Cytosolic Aconitase Activity. Am. J. Physiol.-Cell Physiol. 2005, 288, C1117–C1124. [Google Scholar] [CrossRef] [Green Version]
  53. Koh, H.J.; Lee, S.M.; Son, B.G.; Lee, S.H.; Ryoo, Z.Y.; Chang, K.T.; Park, J.W.; Park, D.C.; Song, B.J.; Veech, R.L.; et al. Cytosolic NADP+-Dependent Isocitrate Dehydrogenase Plays a Key Role in Lipid Metabolism. J. Biol. Chem. 2004, 279, 39968–39974. [Google Scholar] [CrossRef] [Green Version]
  54. Tong, W.H.; Rouault, T.A. Metabolic Regulation of Citrate and Iron by Aconitases: Role of Iron-Sulfur Cluster Biogenesis. BioMetals 2007, 20, 549–564. [Google Scholar] [CrossRef]
  55. Moreno, M.; Ortega, F.; Xifra, G.; Ricart, W.; Fernández-Real, J.M.; Moreno-Navarrete, J.M. Cytosolic Aconitase Activity Sustains Adipogenic Capacity of Adipose Tissue Connecting Iron Metabolism and Adipogenesis. FASEB J. 2015, 29, 1529–1539. [Google Scholar] [CrossRef] [Green Version]
  56. Tretter, L.; Adam-Vizi, V. Alpha-Ketoglutarate Dehydrogenase: A Target and Generator of Oxidative Stress. Philos. Trans. R. Soc. B Biol. Sci. 2005, 360, 2335–2345. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Zdzisińska, B.; Żurek, A.; Kandefer-Szerszeń, M. Alpha-Ketoglutarate as a Molecule with Pleiotropic Activity: Well-Known and Novel Possibilities of Therapeutic Use. Arch. Immunol. Ther. Exp. 2017, 65, 21–36. [Google Scholar] [CrossRef] [Green Version]
  58. Krebs, H.A. Metabolism of Amino-Acids. Biochem. J. 1935, 29, 1951–1969. [Google Scholar] [CrossRef] [PubMed]
  59. Owen, O.E.; Kalhan, S.C.; Hanson, R.W. The Key Role of Anaplerosis and Cataplerosis for Citric Acid Cycle Function. J. Biol. Chem. 2002, 277, 30409–30412. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Mittendorfer, B.; Volpi, E.; Wolfe, R.R. Whole Body and Skeletal Muscle Glutamine Metabolism in Healthy Subjects. Am. J. Physiol.-Endocrinol. Metab. 2001, 280, E323–E333. [Google Scholar] [CrossRef] [PubMed]
  61. Cruzat, V.F. Glutamine and Skeletal Muscle. In Nutrition and Skeletal Muscle; Academic Press: San Diego, CA, USA, 2019; pp. 299–313. ISBN 9780128104224. [Google Scholar]
  62. Keogh, K.; Kenny, D.A.; Cormican, P.; McCabe, M.S.; Kelly, A.K.; Waters, S.M. Effect of Dietary Restriction and Subsequent Re-Alimentation on the Transcriptional Profile of Bovine Skeletal Muscle. PLoS ONE 2016, 11, e0149373. [Google Scholar]
  63. Goudarzi, A. The Recent Insights into the Function of ACAT1: A Possible Anti-Cancer Therapeutic Target. Life Sci. 2019, 232, 116592. [Google Scholar] [CrossRef]
  64. Li, H.; Ericsson, M.; Rabasha, B.; Budnik, B.; Chan, S.H.; Freinkman, E.; Lewis, C.A.; Doench, J.G.; Wagner, B.K.; Garraway, L.A.; et al. 6-Phosphogluconate Dehydrogenase Links Cytosolic Carbohydrate Metabolism to Protein Secretion via Modulation of Glutathione Levels. Cell Chem. Biol. 2019, 26, 1306–1314.e5. [Google Scholar] [CrossRef]
  65. Ying, W. NAD+/NADH and NADP+/NADPH in Cellular Functions and Cell Death: Regulation and Biological Consequences. Antioxid. Redox Signal. 2008, 10, 179–206. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Venn diagram of the proteins identified. Number of proteins identified in each treatment (exclusive), the intercept containing the number of proteins common in both treatments. RM = restriction maintenance; MR = maintenance restriction.
Figure 1. Venn diagram of the proteins identified. Number of proteins identified in each treatment (exclusive), the intercept containing the number of proteins common in both treatments. RM = restriction maintenance; MR = maintenance restriction.
Animals 12 01011 g001
Figure 2. Volcano plot comparing protein abundance fold changes (FC) between treatments (ratio RM/MR). Differentially abundant proteins (DAPs; FDR < 0.05) are highlighted in red and green. DAPs highlighted in green are upregulated and DAPs highlighted in red are downregulated. The black squares represent the non-significant proteins (FDR > 0.05).
Figure 2. Volcano plot comparing protein abundance fold changes (FC) between treatments (ratio RM/MR). Differentially abundant proteins (DAPs; FDR < 0.05) are highlighted in red and green. DAPs highlighted in green are upregulated and DAPs highlighted in red are downregulated. The black squares represent the non-significant proteins (FDR > 0.05).
Animals 12 01011 g002
Figure 3. Analyses of the enriched biological process from the differentially abundant proteins (DAPs). The triangles represent the protein names, while the circles represent the biological processes. Different node color means distinct functional group, and node size means the significance. The most significant (FDR < 0.05) term is labeled with bold letters.
Figure 3. Analyses of the enriched biological process from the differentially abundant proteins (DAPs). The triangles represent the protein names, while the circles represent the biological processes. Different node color means distinct functional group, and node size means the significance. The most significant (FDR < 0.05) term is labeled with bold letters.
Animals 12 01011 g003
Figure 4. Summarized metabolic pathways influenced by the experimental treatments. Glycolysis is represented in gray, glycogenolysis in purple, pentose phosphate pathway in pink, oxidative decarboxylation of pyruvate to form acetyl-CoA in green, citrate cycle in yellow, cytoplasmic conversion of citrate into α-ketoglutarate in blue, glutamine synthesis in orange, and fatty acid degradation in red. Proteins inside orange squares were exclusively expressed in MR or were more abundant in MR compared to RM. Proteins inside blue squares were exclusively expressed in RM or were more abundant in RM compared to MR. α-KGDH: α- ketoglutarate dehydrogenase complex; ACAT1: acetoacetyl-CoA thiolase; ACO1: aconitase 1; AGL: glycogen debranching enzyme; DLAT: dihydrolipoyl transacetylase; DLD: dihydrolipoamide dehydrogenase; DLST: dihydrolipoyllysine-residue succinyltransferase; ECHS1: enoyl CoA hydratase, short chain 1; ENO2: enolase 2; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; GLUL: glutamine synthetase; IDH1: isocitrate dehydrogenase 1; IDH3A: isocitrate dehydrogenase 3 catalytic subunit alpha; OGDH: 2-oxoglutarate dehydrogenase; PDC: pyruvate dehydrogenase complex; PDHA1: pyruvate dehydrogenase E1 subunit alpha 1; PFKM: phosphofructokinase; PGD: 6-phosphogluconate dehydrogenase; PKM: pyruvate kinase; TPI1: triosephosphate isomerase.
Figure 4. Summarized metabolic pathways influenced by the experimental treatments. Glycolysis is represented in gray, glycogenolysis in purple, pentose phosphate pathway in pink, oxidative decarboxylation of pyruvate to form acetyl-CoA in green, citrate cycle in yellow, cytoplasmic conversion of citrate into α-ketoglutarate in blue, glutamine synthesis in orange, and fatty acid degradation in red. Proteins inside orange squares were exclusively expressed in MR or were more abundant in MR compared to RM. Proteins inside blue squares were exclusively expressed in RM or were more abundant in RM compared to MR. α-KGDH: α- ketoglutarate dehydrogenase complex; ACAT1: acetoacetyl-CoA thiolase; ACO1: aconitase 1; AGL: glycogen debranching enzyme; DLAT: dihydrolipoyl transacetylase; DLD: dihydrolipoamide dehydrogenase; DLST: dihydrolipoyllysine-residue succinyltransferase; ECHS1: enoyl CoA hydratase, short chain 1; ENO2: enolase 2; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; GLUL: glutamine synthetase; IDH1: isocitrate dehydrogenase 1; IDH3A: isocitrate dehydrogenase 3 catalytic subunit alpha; OGDH: 2-oxoglutarate dehydrogenase; PDC: pyruvate dehydrogenase complex; PDHA1: pyruvate dehydrogenase E1 subunit alpha 1; PFKM: phosphofructokinase; PGD: 6-phosphogluconate dehydrogenase; PKM: pyruvate kinase; TPI1: triosephosphate isomerase.
Animals 12 01011 g004
Table 1. Differentially abundant proteins (DAPs) in the skeletal muscle of the offspring.
Table 1. Differentially abundant proteins (DAPs) in the skeletal muscle of the offspring.
AccessionProtein NameGene NameFDR 1Fold Change 2
(RM vs. MR)
A0A452DPE6Calponin homology (CH) domain-containing proteinSMTNL10.0398564−0.90
A0A452ET82Pyruvate kinasePKM0.0398564−0.90
A0A452FHP7Uncharacterized proteinFLNC0.00021620.25
A0A452FBM2Uncharacterized proteinHSPA90.01389940.33
A0A452G4K3ATP-dependent 6-phosphofructokinasePFKM0.03985640.36
A0A452FWD9; A0A452FWC6Dihydrolipoyl dehydrogenaseDLD0.03985640.42
A0A452FX48I; A0A452FWP3Uncharacterized proteinAGL0.00021620.49
A0A452ESM1; A0A452ERW8; A0A452ERN3Lipoyl-binding domain-containing proteinDLST0.04461520.53
A0A452EYA9; A0A452EYB9; A0A452EY59Uncharacterized proteinSELENBP10.00069710.66
A0A452FIG7S-formylglutathione hydrolaseESD0.03985640.66
A0A452EJS4Uncharacterized proteinMYOM10.02197530.69
A0A452DWL1IF rod domain-containing proteinDES0.00000030.99
A0A452ET55Triosephosphate isomeraseTPI10.03213581.18
RM = restriction maintenance; MR = maintenance restriction; 1 FDR= false discovery rate; 2 negative and positive fold change indicates the less and more abundant proteins in the treatment RM compared to MR.
Table 2. Enriched metabolic pathways of the exclusive proteins from treatment RM and MR.
Table 2. Enriched metabolic pathways of the exclusive proteins from treatment RM and MR.
KEGG IDDescriptionFDR 1Protein Names 2
Treatment RM
oas04530Tight junction9.44 × 10−8ACTN1, ACTN4, EZR, LOC443340, MSN, MYH1, MYH13, MYH4, MYH8, MYL6, OMYHC2A, PPP2R1A, PPP2R1B, RDX, RHOA, YWHAQ
oas01200Carbon metabolism7.73 × 10−7ACAT1, ACO1, ADH5, DLAT, ECHS1, HADHA, IDH1, OGDH, OGDHL, PGD
oas00280Valine, leucine, and isoleucine degradation2.98 × 10−5ACAT1, ALDH2, ECHS1, HADHA, HIBADH, HSD17B10, OXCT1
oas03015mRNA surveillance pathway5.52 × 10−5EIF4A3, PABPC1, PABPC4, PPP1CA, PPP1CB, PPP1CC, PPP2R1A, PPP2R1B
oas00071Fatty acid degradation6.03 × 10−5ACAT1, ADH5, ALDH2, ECHS1, ECI1, HADHA
oas03013RNA transport6.03 × 10−5EEF1A2, EIF4A1, EIF4A2, EIF4A3, EIF4EBP1, PABPC1, PABPC4, SEC13, SUMO2, SUMO3
oas04720Long-term potentiation6.03 × 10−5PPP1CA, PPP1CB, PPP1CC, PPP1R1A, PPP3CA, PPP3CB, PPP3CC
oas00380Tryptophan metabolism9.75 × 10−5ACAT1, ALDH2, ECHS1, HADHA, OGDH, OGDHL
oas04810Regulation of actin cytoskeleton1.30 × 10−4ACTN1, ACTN4, EZR, LOC443340, MSN, PPP1CA, PPP1CB, PPP1CC, RDX, RHOA
oas00020Citrate cycle (TCA cycle)1.80 × 10−4ACO1, DLAT, IDH1, OGDH, OGDHL
Treatment MR
oas00010Glycolysis/Gluconeogenesis2.68 × 10−5AKR1A1, ALDH7A1, ENO2, PDHA1, PDHA2, GAPDH
oas00020Citrate cycle (TCA cycle)2.68 × 10−5IDH3A, IDH3G, PDHA1, PDHA2
oas01230Biosynthesis of amino acids2.68 × 10−5ARG1, ENO2, GLUL, IDH3A, IDH3G, GAPDH
oas01200Carbon metabolism9.08 × 10−5ENO2, IDH3A, IDH3G, PDHA1, PDHA2, GAPDH
oas01100Metabolic pathways0.00064AKR1A1, ALDH7A1, ARG1, ENO2, GLUL, IDH3A, IDH3G, LAP3, MTHFD1, PDHA1, PDHA2, GAPDH
oas04066HIF-1 signaling pathway0.00074EIF4E, ENO2, PDHA1, PDHA2, GAPDH
oas00620Pyruvate metabolism0.0012ALDH7A1, PDHA1, PDHA2
oas00330Arginine and proline metabolism0.0019ALDH7A1, ARG1, LAP3
oas04145Phagosome0.0033RAC1, TUBB, TUBB1, TUBB2A
oas012102-Oxocarboxylic acid metabolism0.0055IDH3A, IDH3G
RM = restriction maintenance; MR = maintenance restriction; 1 FDR= false discovery rate; 2 bold font are the proteins that are present in each treatment that are connected into the protein-protein interaction network.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Costa, T.C.; Dutra, L.L.; Mendes, T.A.d.O.; dos Santos, M.M.; Veroneze, R.; Gionbelli, M.P.; Duarte, M.d.S. Impact of Maternal Feed Restriction at Different Stages of Gestation on the Proteomic Profile of the Newborn Skeletal Muscle. Animals 2022, 12, 1011. https://0-doi-org.brum.beds.ac.uk/10.3390/ani12081011

AMA Style

Costa TC, Dutra LL, Mendes TAdO, dos Santos MM, Veroneze R, Gionbelli MP, Duarte MdS. Impact of Maternal Feed Restriction at Different Stages of Gestation on the Proteomic Profile of the Newborn Skeletal Muscle. Animals. 2022; 12(8):1011. https://0-doi-org.brum.beds.ac.uk/10.3390/ani12081011

Chicago/Turabian Style

Costa, Thaís Correia, Luana Lucas Dutra, Tiago Antônio de Oliveira Mendes, Marta Maria dos Santos, Renata Veroneze, Mateus Pies Gionbelli, and Marcio de Souza Duarte. 2022. "Impact of Maternal Feed Restriction at Different Stages of Gestation on the Proteomic Profile of the Newborn Skeletal Muscle" Animals 12, no. 8: 1011. https://0-doi-org.brum.beds.ac.uk/10.3390/ani12081011

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop