Next Article in Journal
Detecting the Conspecific: Herbivory-Induced Olfactory Cues in the Fall Armyworm (Lepidoptera: Noctuidae)
Next Article in Special Issue
Correlation between Pre-Ovulatory Follicle Diameter and Follicular Fluid Metabolome Profiles in Lactating Beef Cows
Previous Article in Journal
Increased Fibroblast Growth Factor 21 (FGF21) Concentration in Early Second Trimester Amniotic Fluid and Its Association with Fetal Growth
Previous Article in Special Issue
Uterine Metabolomic Analysis for the Regulation of Eggshell Calcification in Chickens
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 11
Issue 9
10.3390/metabo11090582
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Maternal Undernutrition during Pregnancy Alters Amino Acid Metabolism and Gene Expression Associated with Energy Metabolism and Angiogenesis in Fetal Calf Muscle

by
Susumu Muroya
1,*,
Yi Zhang
2,
Aoi Kinoshita
2,
Kounosuke Otomaru
3,
Kazunaga Oshima
4,
Yuji Gotoh
4,
Ichiro Oshima
2,
Mitsue Sano
5,
Sanggun Roh
6,
Mika Oe
1,
Koichi Ojima
1 and
Takafumi Gotoh
2,*
1
Division of Animal Products Research, NARO Institute of Livestock and Grassland Science (NILGS), Tsukuba 305-0901, Ibaraki, Japan
2
Department of Agricultural Sciences and Natural Resources, Kagoshima University, Korimoto 1-21-24, Kagoshima 890-8580, Japan
3
Joint Faculty of Veterinary Medicine, Kagoshima University, Korimoto 1-21-24, Kagoshima 890-8580, Japan
4
Division of Year-Round Grazing Research, NARO Western Region Agricultural Research Center, 60 Yoshinaga, Ohda 694-0013, Shimane, Japan
5
Faculty of Human Culture, University of Shiga Prefecture, 2500 Hassaka-cho, Hikone 522-8533, Shiga, Japan
6
Graduate School of Agricultural Science, Tohoku University, 468-1 Aoba, Aramaki, Aoba-ku, Sendai 980-8578, Miyagi, Japan
*
Authors to whom correspondence should be addressed.
Metabolites 2021, 11(9), 582; https://doi.org/10.3390/metabo11090582
Submission received: 21 July 2021 / Revised: 26 August 2021 / Accepted: 27 August 2021 / Published: 28 August 2021
(This article belongs to the Special Issue Metabolomic Applications in Animal Science Volume 2)
Browse Figures
Versions Notes

Abstract

:
To elucidate the mechanisms underlying maternal undernutrition (MUN)-induced fetal skeletal muscle growth impairment in cattle, the longissimus thoracis muscle of Japanese Black fetal calves at 8.5 months in utero was analyzed by an integrative approach with metabolomics and transcriptomics. The pregnant cows were fed on 60% (low-nutrition, LN) or 120% (high-nutrition, HN) of their overall nutritional requirement during gestation. MUN markedly decreased the bodyweight and muscle weight of the fetus. The levels of amino acids (AAs) and arginine-related metabolites including glutamine, gamma-aminobutyric acid (GABA), and putrescine were higher in the LN group than those in the HN group. Metabolite set enrichment analysis revealed that the highly different metabolites were associated with the metabolic pathways of pyrimidine, glutathione, and AAs such as arginine and glutamate, suggesting that MUN resulted in AA accumulation rather than protein accumulation. The mRNA expression levels of energy metabolism-associated genes, such as PRKAA1, ANGPTL4, APLNR, CPT1B, NOS2, NOS3, UCP2, and glycolytic genes were lower in the LN group than in the HN group. The gene ontology/pathway analysis revealed that the downregulated genes in the LN group were associated with glucose metabolism, angiogenesis, HIF-1 signaling, PI3K-Akt signaling, pentose phosphate, and insulin signaling pathways. Thus, MUN altered the levels of AAs and expression of genes associated with energy expenditure, glucose homeostasis, and angiogenesis in the fetal muscle.

1. Introduction

Placental insufficiency and various other maternal factors, such as environmental stress, disease, and malnutrition, contribute to intrauterine growth retardation (IUGR) during gestation, which results in impaired fetal development and decreased bodyweight (BW) at birth [1]. IUGR, also called fetal growth restriction (FGR), which leads to stunting of postnatal growth and decreased nutrient utilization efficiency in the offspring, adversely affects the whole-body composition and results in prolonged disruption of homeostasis as it predisposes the fetus to metabolic disorders, such as insulin resistance [1,2]. Maternal undernutrition (MUN) is one of the major factors that contribute to the induction of IUGR. Previous studies have reported that MUN impairs the growth and metabolism of fetal organs, including the skeletal muscle and liver [3,4,5,6,7,8], by altering the secretion and/or sensitivity of glucose, insulin, and insulin-like growth factor-1 (IGF1) [8]. The skeletal muscle accounts for 40–45% of the body mass [1] and 80% uptake of whole-body insulin-stimulated glucose. Thus, MUN-induced impairment of skeletal muscle growth and metabolism adversely affect fetal metabolism by modulating whole-body insulin sensitivity and glucose homeostasis [9].
Recently, we demonstrated that the skeletal muscle weight of fetuses of dams fed with 60% of crude protein (CP), crude fat (CF), and energy requirements was 1.27- to 1.72-fold lower than that of the fetuses of dams fed with 120% of the nutritional requirements [10]. The phenotypic alteration observed in this study promotes understanding the mechanism underlying the effect of MUN on muscle development and tissue metabolism in the fetuses. Previous studies have demonstrated that MUN from the mid-first to mid-second trimester of gestation period in sheep decreases the number of secondary fibers per fasciculi in the fetal longissimus dorsi and semitendinosus muscles [11,12]. In human fetuses diagnosed with IUGR, the DNA content decreased in the skeletal muscle at late gestation but did not show difference in protein/DNA ratios compared with healthy fetuses [13]. IUGR alters the proliferation and differentiation of myoblasts prepared from near-term sheep fetuses [14] and the proportion of type I (oxidative) and II (glycolytic) myofibers, the two major matured muscle cell types in the skeletal muscle tissue [15,16,17]. Fetal muscle growth impairment in meat animals can affect livestock production, especially meat quality parameters such as beef tenderness [18,19]. The intramuscular triglyceride content is upregulated in the nutrition-restricted fetal skeletal muscles [16]. This indicated that the skeletal muscles of nutrition-restricted fetuses are predisposed to insulin resistance by fetal programming, which can lead to intramuscular fat accumulation in beef, an important attribute for high marbled beef, such as Wagyu (Japanese Black) cattle.
In contrast to the effect on the myogenic process, the effect of MUN on metabolism of fetal skeletal muscle tissue remains poorly understood. The fetal skeletal muscle mass was reduced by MUN [10], possibly due to the retardation of whole process of skeletal muscle development. This could be a result of insufficient supply of energy and nutrient components that were required for the myogenic differentiation, muscle fiber maturation, and protein accumulation. Skeletal muscle metabolism in nutrient-restricted fetuses is hypothesized to be forced to prioritize minimal metabolism required for cellular energy and structure maintenance to sustain fetal muscle viability. Hence, metabolic pathways associated with energy metabolism and glucose homeostasis in the fetal skeletal muscle under MUN may be dysregulated with concomitant modulation in gene expression. In pigs, the concentrations of arginine, ornithine, proline, glutamine, and polyamines were markedly downregulated in the skeletal muscle under MUN [20]. The dysregulation of amino acids (AAs) and/or lipids in the muscles of the fetus with IUGR has been reported in the nutrition-restricted ewes [21] and placenta insufficiency model under high temperature and humidity in sheep [22] although the altered metabolites varied depending on the experimental conditions. MUN exerts its effect on the fetus by decreasing insulin sensitivity and promoting hypoxia [23]. Furthermore, MUN altered the gene expression profile in the fetal skeletal muscle tissue [24,25,26,27,28]. These studies have suggested that alteration of fetal skeletal muscle metabolism caused by MUN is crucial for the muscle development, however, the effect of MUN on fetal skeletal muscle metabolism has been poorly elucidated. This could be partly due to lack of experimental system leading to phenotypic alteration of fetal development by MUN, which can be overcome by setting the nutrient levels for dams to 60 and 120% as the low and high levels, respectively [10].
The effect of MUN on the fetal development is different depending on the gestational stage of cows, the extent and period of MUN, and the types of nutrient restricted [29]. In addition, difference in the result of MUN not only between animal species, breeds, or organs, but also between skeletal muscle types can be observed [24,30]. Comparison of experimental conditions among previous studies suggests that the effect of MUN on skeletal muscle gene expression varies between early and middle-to-late pregnancy of dams, between energy and protein restriction, and between levels of nutrition [24,25,26,27,28]. Notably, no significant phenotypic effect was observed in skeletal muscle mass or BW, in cases that the MUN period was limited to early-to-middle gestation [12] or that the nutritional level was compared between 70 and 100% or between 85 and 140% in requirement of energy for maintenance [30,31], even with significant effect on muscle fiber numbers and/or gene expression. Thus, the period of MUN and gestational stage of dams can alter the effect of MUN on the gene expression, metabolism, and phenotype of the fetal skeletal muscle. To elucidate the molecular basis underlying the phenotypic alteration of fetal muscle caused by MUN, an experimental design that causes significant phenotypic alteration is necessary.
This study aimed to elucidate the effect of MUN on fetal skeletal muscle metabolism that links to the retardation of the fetal muscle development in Japanese Black cows. To this end, we used a design in which the pregnant dams were fed on low-nutrition (LN) and high-nutrition (HN) diets (based on protein, fat, and energy contents) during the whole gestation (until month 8.5 post-conception) [10]. To find the phenotypic effect of MUN on fetal development under the limited knowledge regarding the MUN effect in cattle, this design is beneficial due to its marked phenotypic alteration of fetal muscles in the LN group [10]. The LN and HN groups were set to 60 and 120% of the recommended nutritional level, respectively, based on the standard diet model for prepregnant BW in the Japanese Feeding Standard for Beef Cattle (JFSBC, 2008 ed.) [32]. The changes in metabolomic and transcriptomic profiles in the fetal longissimus thoracis (LT) muscle were analyzed using capillary electrophoresis–time-of-flight mass spectrometry (CE-TOFMS) and microarray analysis, respectively. Bioinformatic analysis of the metabolomic and transcriptomic data was performed to understand the impact of MUN on the fetal skeletal muscle metabolism.

2. Results

2.1. Fetal Carcass Traits

To examine the effect of MUN on fetal skeletal muscle growth, the BW, LT muscle weight, and total skeletal muscle weight of half carcass (right side) were measured (Table 1). The BW, LT muscle weight, and total muscle weight of the LN group were lower than those of the HN group (p < 0.01), whereas the percentage of total muscle weight in BW and that of the LT muscle weight in total muscle weight were not different (p > 0.05). The ratios of BW, total muscle weight, and LT muscle weight in the LN group to those in the HN group were 0.72, 0.72, and 0.69, respectively. Thus, MUN markedly decreased the mass of skeletal muscles in the fetuses of the LN group.

2.2. Metabolomic Profile of Fetal Skeletal Muscle

To assess the effect of MUN on fetal LT muscle metabolism, the metabolomic difference was analyzed using CE-TOFMS (Table 2). In the CE-TOFMS analysis, four fetuses with the highest BW in the HN group and four fetuses with the lowest BW in the LN group were used for the analysis. In total, 169 uniquely annotated metabolites were detected (Table S1). Of these metabolites, the levels of carnosine, glutamine, glycerol, creatine, N6-methyllysine, phosphorylcholine, phenylalanine, proline, (p < 0.05), alanine, putrescine, creatinine, gamma-aminobutyric acid (GABA), histidine, and glycine (p < 0.10) in the fetal muscles of the LN group were higher than those in the fetal muscles of the HN group (Table 2). The levels of 2-aminoethylphosphonic acid (2-AEP, also called ciliatine), myo-inositol 2-phosphate (also called inositol 5-phosphate; Ins2P) (p = 0.039), and 2-hydroxyvaleric acid (2-HVA) (p = 0.098) in the LN fetal muscles were lower than those in the HN fetal muscles.
Unsupervised hierarchical clustering analysis (HCA) using the top 50 significantly different metabolites between the LN and HN groups classified the fetal muscle samples into LN and HN groups. Figure 1 shows that glutamine, glycerol, phosphorylcholine, carnosine, N6-methyllysine, creatinine, putrescine, creatine, valine, and allantoin were abundant in the LN group, whereas 2-AEP, Ins2P, taurine, S-adenosylmethionine (SAM), UDP-N-acetylgalactosamine/UDP-N-acetylglucosamine, N5-ethylglutamine, terephthalic acid, 2-HVA, guanosine, and NAD+ were abundant in the HN group, compared with the counterpart group. This result indicated the important contribution of these metabolites to the separation of the fetal muscles between the LN and HN groups. Thus, the fetal muscle of the LN group was characterized by the abundance of metabolites such as carnosine, glutamine, glycerol, creatine, N6-methyllysine, phenylalanine, and phosphorylcholine.
To understand biological processes associated with differentially expressed metabolites between the LN and HN groups, metabolite set enrichment analysis (MSEA) was performed using the top 50 differentially expressed metabolites. The differentially expressed metabolites were significantly associated with the metabolic pathways of pyrimidine (p = 0.003), aminoacyl-tRNA biosynthesis (p = 0.007), glycerolipid (p = 0.007), AAs (arginine, alanine, aspartate, glutamate, histidine, and proline) (p < 0.05), and glutathione (p = 0.011) (Table 3). The metabolisms associated with glyoxylate and dicarboxylate (p = 0.013), and phosphonate and phosphinate (p = 0.017), glycerophospholipid (p = 0.018), galactose (p = 0.022), primary bile acid biosynthesis (p = 0.025), purine (p = 0.025), D-glutamine and D-glutamate (p = 0.029), nitrogen (p = 0.029), glycine, serine and threonine (p = 0.036), and β-alanine (p = 0.040) were also significantly extracted (Table 3).

2.3. Effect of MUN on the Gene Expression Profile of the Fetal Skeletal Muscle

The MUN-mediated metabolic alterations in the LT muscles of the LN fetuses could result from disturbed gene expression. To investigate the association of MUN-induced metabolic alterations with gene expression, microarray-based gene expression analysis was performed using the pooled muscle RNA samples. The RNA samples were extracted from four fetuses with the highest BW in the HN group and four fetuses with the lowest BW in the LN group. Of the 12,786 detected unique genes, the expression levels of 219 genes in the fetal muscle varied > two-fold between the LN and HN groups (187 upregulated and 32 downregulated genes in the LN group). The top five downregulated genes in the LN group were retinoid X receptor α(RXRA), angiopoietin-like 4 (ANGPTL4), uncoupling protein 2 (UCP2), and apelin receptor (APLNR, also called APJ). Meanwhile, the top five upregulated genes in the LN group were syndecan-binding protein (SDCBP), olfactomedin 4 (OLFM4), nuclear receptor subfamily 2 group E member 3 (NR2E3), tribbles pseudokinase 3 (TRIB3), and an uncharacterized gene (XM_005208051). The downregulation of ANGPTL4, UCP2, and ALPNR, which are involved in energy homeostasis, in the fetuses of the LN group was validated using qRT-PCR (p < 0.05; Figure 2). However, the expression levels of RXRA, SDCBP, OLFM4, NR2E3, and TRIB3 were not significantly different between the groups in the qRT-PCR analysis (p > 0.10).
The microarray analysis results suggested that the MUN affected energy metabolism in the fetal LT muscle of the LN group. Moreover, MSEA revealed that MUN altered the metabolic pathways associated with AAs, glutathione, nitrogen, and one-carbon cycle. Therefore, this study further focused mainly on energy metabolism-related genes that showed difference between the LN and HN groups in the microarray analysis. The qRT-PCR analysis revealed that the expression levels of carnitine palmitoyltransferase 1B (CPT1B) (p = 0.027), NOS2 (also called inducible NOS, iNOS) (p = 0.024), and pyruvate dehydrogenase kinase 4 (PDK4) (p = 0.094) were downregulated, whereas those of solute carrier family 7A5 (SLC7A5) (p = 0.034), methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2), aldehyde dehydrogenase 1 family member L2 (ALDH1L2), peroxisome proliferator-activated receptor α (PPARA), resistin (RETN), glycoprotein nonmetastatic melanoma protein B (GPNMB), rapamycin-insensitive companion of mammalian target of rapamycin (RICTOR), and forkhead box protein P1 (FOXP1) were upregulated in the LN fetal muscles (p < 0.10) (Figure 2).

2.4. GO Analysis of MUN-Associated Metabolic Pathways

The analysis of differentially expressed genes between the LN and HN groups revealed that MUN modulated gene expression associated with energy production/consumption, NO synthesis, angiogenesis, and the metabolic pathways involved in skeletal muscle development. To understand the pathways in which the differentially expressed genes are enriched, gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes database (KEGG) pathway analyses were performed using 3949 genes whose expression in the LN fetal muscle was >1.25-fold downregulated (Table 4) or upregulated (Table 5) when compared with those in the HN fetal muscles.
The downregulated genes in the LN fetal muscle were significantly enriched in the following KEGG pathways: glycolysis/gluconeogenesis, MAPK signaling, HIF-1 signaling, carbon metabolism, PI3K-Akt signaling, thyroid hormone signaling, pentose phosphate, insulin signaling, and p53 signaling (p < 0.01, Table 4). In addition, the downregulated genes were associated with GO terms, such as response to cytokine, negative regulation of transcription from RNA polymerase II promoter, angiogenesis, negative regulation of lipid catabolic process, positive regulation of blood vessel endothelial cell migration, negative regulation of protein catabolic process, as well as with those related to glycolysis/gluconeogenesis, apoptosis, and release of cytochrome c (biological process) (p < 0.01). Meanwhile, the upregulated genes in the LN fetal muscle were enriched in the following KEGG pathways: biosynthesis of AAs, aminoacyl-tRNA biosynthesis, glycine, serine and threonine metabolism, synthesis and degradation of ketone, and protein digestion and absorption (p < 0.01, Table 5). The upregulated genes were also associated with the following GO terms (biological process): one-carbon (1C) metabolic process and tetrahydrofolate metabolic process (p < 0.01).
To validate the resultant biological processes in the GO/pathway analyses, the expression levels of the representative genes involved in these pathways or metabolisms were examined with qRT-PCR. The expression levels of genes related to glycolysis/gluconeogenesis or carbon metabolism (GPI, ENO3, PFKM, FBP2, and PGAM1), MAPK signaling, HIF-1 signaling, and PI3K-Akt signaling pathways (ANGPT4, NOS2, NOS3, ENO3, MAPK3, NFKB1, PLCG1, PRKAA1, RELA, and/or THBS1), pentose phosphate pathway (GPI, PFKM, and FBP2), and insulin signaling pathway (FBP2 and MAPK3) were downregulated in the LN fetal muscle (Table 4; Figure 3). In addition, the expression level of the cellular energy sensor PRKAA1 (cyclic AMP-activated protein kinase (AMPK) α1 subunit) was downregulated in the LN fetal muscle (p = 0.029), although the AMP/ATP ratio was not significantly different between the LN and HN groups based on the metabolome analysis (p = 0.462). Most of the downregulated genes of interest were also associated with GO terms, such as response to cytokine (NFKB1 and RELA), negative regulation of transcription from RNA polymerase II promoter (FOXP1, NFKB1, and RELA), angiogenesis (GPI, NOS3, ANGPT4, and ANGPTL4), regulation of apoptotic process (PINK1 and RELA), negative regulation of lipid catabolic process (PRKAA1), positive regulation of blood vessel endothelial cell migration (ANGPTL4, PLCG1, and THBS1), and negative regulation of protein catabolic process (NOS2, SIRT2, and RELA) (Table 4).
The qRT-PCR analysis also revealed that the genes associated with synthesis and degradation of ketone bodies (BDH1) and one-carbon metabolic process (MTHFD2 and ALDH1L2) were upregulated in the LN fetal muscle (Table 5). The associated genes validated for the upregulated genes were markedly lower than those for the downregulated genes. These results indicated that the downregulated genes were associated with transcription, angiogenesis, glucose metabolism, AA metabolism, and protein catabolism, and MAPK, HIF-1, and PI3K-Akt signaling pathways.

3. Discussion

3.1. Altered Amino Acid Metabolism Was the Major Response to IUGR

In this study, MUN altered the metabolomic profile of the fetal LT muscle. Among various maternal nutritional restriction conditions, this study examined the effect of restricting global nutrition, including CP, CF, and total digestible nutrients, during whole gestation. This nutrient-restricted condition markedly decreased BW, skeletal muscle weight, and fat content in the fetuses. The weight of liver, kidney, thymus, spleen, heart, lung, rumen, omasum, and large intestine were also lower in fetuses of the LN group than in those of the HN group as reported previously [10].
The retardation of fetal organs such as liver and kidney might affect the muscle growth retardation and disturbance of metabolism. In addition, MUN altered the fetal LT muscle metabolome. The levels of AAs, such as glutamine, phenylalanine, and proline in the LN fetal muscle were 1.2- to 1.5-fold higher than those in the HN fetal muscles. The increased levels of AAs, such as glutamine and alanine, were consistent with the results of MSEA, which revealed that metabolic pathways associated with arginine, alanine, aspartate, glutamate, histidine, and proline in the fetal muscles were significantly different between the LN and HN groups (Table 3). On the other hand, levels of ATP, AMP, and AMP/ATP ratio were not significantly different between the LN and HN groups. This suggests that energy balance was maintained in the LN fetal muscle to the same level as the HN fetal muscle, at the expense of energy consumption in nutrient-restricted condition. Hence, the increased AA levels in the LN fetal muscles can be attributed to the decreased utilization of AAs caused by energy insufficiency, increased AA uptake, and/or enhanced protein degradation (Figure 4).
AAs, such as glutamine and alanine, accumulate likely due to decreased protein synthesis even under optimal AA supply conditions, which can be due to restricted available energy. Hypothetically, AAs that are constitutively taken up from the circulation or generated through intracellular protein degradation accumulate in the muscle cells due to their decreased utilization in protein synthesis. This is supported by the decreased muscle mass and downregulation of genes involved in the PI3K-Akt signaling and insulin signaling pathways that activate protein synthesis [33], in the LN fetal muscles (Table 4). Similar results were also reported in previous studies, which demonstrated that the skeletal muscle mass decreased even though the AA levels were increased [22,34]. The decreased AA utilization could be associated with the downregulated expression of genes involved in the glycolytic and pentose phosphate metabolic pathways (Table 4) and decreased energy production in the LN fetal muscles. This restricted energy-mediated AA accumulation is also supported by decreased glucose uptake, which can be attributed to downregulated expression of glucose transporter 4 (GLUT4) observed previously (in submission). These changes in energy metabolism can force the fetal metabolism to prioritize skeletal muscle maintenance by saving energy rather than protein synthesis for cell growth and proliferation (Figure 4).
In this study, the increased levels of glutamine, phenylalanine, proline, and histidine were accompanied by upregulation of SLC7A5 (a large neutral AA transporter) in the LN fetal muscles. This suggests that the AA uptake/release is activated in the bovine fetuses of the LN group. High temperature-induced IUGR decreased the uptake of BCAAs and aminogenic AAs, and decreased the levels of metabolites, such as alanine, arginine, aspartate, and glutamate in the ovine fetal muscle [22]. These results indicate that IUGR affects the AA uptake in the skeletal muscle, although the differences in the altered AAs between the bovine and ovine models can be attributed to the experimental conditions. In addition, the contribution of protein degradation to the increase in AAs cannot be ruled out. The fetal muscle of the LN group can catabolize its proteins as evidenced by the GO analysis results that the downregulated genes were associated with the negative regulation of the protein catabolic process. In previous studies examining MUN impact on proteolytic system, the contents of calpastatin, an endogenous specific inhibitor of protease calpains, and ubiquitinated proteins were not lower in the fetal muscle under MUN, compared with those in the control fetuses [31,35]. Collectively, protein degradation is considered not to significantly contribute to the increased levels of AA in this study.
The five increased proteinogenic AAs (glutamine, proline, alanine, histidine, and glycine) in the LN fetal muscles are glucogenic, whereas phenylalanine is both glucogenic and ketogenic [36,37]. In particular, glutamine was the most abundant proteinogenic AA (3879–5444 nmol/g tissue in the fetal muscle) and is important for glucose synthesis for energy production [38,39]. Previous studies have reported that gluconeogenesis [40] and urea synthesis [41,42] are activated in the nutrient-restricted ovine fetuses. Therefore, the AAs in the LN fetal muscle may be transported to the liver and utilized for energy production through the glucose-alanine cycle.

3.2. Association of Altered AAs with Metabolism of Nitrogen, β-Alanine, and Glycerophospholipids

Among the increased AAs, glutamine plays important physiological roles in the synthesis of other AAs and the related metabolites [38,39]. Glutamine is a major precursor of arginine [43] that can be converted to metabolites, including NO, proline, glutamate, creatine, and polyamines, such as putrescine [38], although the content of arginine was not altered in this study. The levels of these metabolites and glutamine were increased in the LN fetal muscle, which suggested that the increased level of glutamine raised the production of these metabolites through the arginine metabolic pathways. This is consistent with the results of MSEA, which revealed that the altered metabolites in the fetuses of the LN group were mainly associated with arginine metabolism (Table 3).
Putrescine, which was one of the increased metabolites, can be generated from glutamine via ornithine. In addition, putrescine is a precursor of bioactive polyamines, spermine and spermidine. These polyamines, with their roles in the metabolism of pyrimidine and purine, are essential for cell growth, proliferation, and differentiation [44]. As MSEA revealed the enrichment of the metabolism of pyrimidine and purine (Table 3), the increased putrescine levels in the LN group can be attributed to decreased nucleotide synthesis. Similarly, the accumulation of catabolic products (proline and GABA) of arginine and glutamine can be attributed to decreased synthesis of further reaction products due to restricted energy (Figure 4).
The level of creatine, a product of arginine catabolism, was also increased. However, the increased glutamine levels may not simply have contributed to increased creatine levels in the LN fetal muscle. This is because even though creatine is abundant in the skeletal muscle, it is generated in the liver and kidney [45]. The enhanced production of creatine can be attributed to creatinine kinase-catalyzed ATP production, in which creatine is generated by transfer of phosphate from creatine phosphate to ADP [46]. Creatine catabolism may also be impaired in the fetuses of the LN group due to the shift of reaction equilibrium toward enhanced ATP production. The increased levels of creatinine, a product of muscle creatine catabolism, can result from creatine accumulation in the fetuses of the LN group (Figure 4). This is consistent with an observation of increased creatine and creatinine levels in the urine of human neonates with IUGR [47].
Carnosine (β-alanyl-L-histidine) was previously reported to be decreased in the ovine fetal muscle under restricted maternal nutrition conditions, while anserine is increased at the late gestation stage [21], which may be associated with altered histidine metabolism. The levels of carnosine and histidine were increased in the fetuses of the LN group. As the LN fetuses showed altered histidine metabolism (Table 3) and increased levels of histidine transporter SLC7A5, the enhanced histidine uptake might contribute to histidine and carnosine accumulation.

3.3. Other Metabolites Altered in the Fetal Muscles of the LN Group

In this study, glycerol was increased in the LN fetal muscles, which may be derived from the intake of glycerol from nutrient-restricted maternal circulation through the placenta during late gestation [48,49]. The effect of MUN on glycerol levels in the fetal muscle has not been examined. Thus, the reason for increased glycerol levels in the fetuses of the LN group is unclear. Nevertheless, pregnant ewes can produce glycerol from the triglycerides in the adipose tissues during late gestation, especially during the fasting period [48,49]. The increased glycerol level in the LN group may indicate that MUN promotes fetal muscle glycerol uptake from the pregnant dams for energy production.
Inositol phosphoglycans (IPGs), including Ins2P and myo-inositol (myoIns), are involved in the insulin signaling pathway [50]. The administration of IPGs to non-diabetic or diabetic rats dose-dependently alleviates hyperglycemia and promotes muscular glycogenesis [51]. Moreover, IPGs released in response to insulin receptor activation promote oxidative glucose metabolism and the tricarboxylic acid cycle [52]. Thus, Ins2P may be involved in the dysregulation of glucose metabolism in the fetuses of the LN group in this study. The fetal plasma inositol concentrations are increased in sheep and pigs with IUGR [53,54,55]. Hence, Ins2P can be a biomarker for IUGR [47,56].
The correlation between IUGR and 2-AEP or 2-HVA has not been elucidated. Rats and other higher vertebrates cannot generate 2-AEP that contains a carbon–phosphorus bond [57,58]. The ruminant protozoa can synthesize protein-bound and lipid-bound forms of 2-AEP, as well as a free form of 2-AEP [59]. Previous studies have reported that 2-AEP is detected in the goat liver, bovine brain and liver, and major tissues of rats, chicken, and humans, including skeletal muscle [59,60]. In the rat liver, 2-AEP is incorporated in phosphonolipids, such as diacylglyceryl-AEP [61], which suggests that 2-AEP is a membrane component. Thus, 2-AEP in the bovine fetus may be taken up from the pregnant cows through the placenta and utilized as a membrane component. Further studies are needed to examine the biological role of dysregulated 2-AEP in the LN fetal muscle.
Phosphorylcholine level was increased in the LN fetal muscle. Free phosphorylcholine can be generated from the catabolism of sphingomyelins as evidenced by the changes in ceramide, sphingolipid, and phospholipid metabolites in the offspring of ewes fed on a nutrient-restricted diet [21]. The roles of ceramide and sphingomyelins in the muscle are poorly understood. However, the altered phosphorylcholine level may indicate that undernutrition modulates sphingolipid signaling by promoting the production of sphingolipid-1-phosphate in the fetal skeletal muscle, which can affect not only the activation of muscle satellite cells from quiescence for postnatal muscle growth and repair [62] and also neuronal development.
N6-methyllysine is a product of protein degradation. The interaction between modified histones and altered metabolites, including methylated lysine, in fetuses with IUGR has not been completely understood. IUGR-induced metabolic alteration can be associated with alteration of histone modification. The decreased N6-acetyllysine levels in the offspring of over-fed ewes at birth suggested epigenetic changes [21]. Thus, N6-methyllysine may be released from specific lysine residues of histones in relation to altered gene expression in the fetuses of the LN group.

3.4. Potential Mechanisms of NOS Downregulation and Its Impact in the Fetuses of the LN Group

MUN downregulated the expression of NOS2, PRKAA1, ANGPTL4, PDK4, APLNR, UCP2, and CPT1B, all of which are associated with energy metabolism. NO synthase (NOS) including NOS2 product (iNOS), plays an essential role in energy production, autoregulation of blood flow, myocyte differentiation, respiration, and glucose homeostasis of skeletal muscle [63], as well as angiogenesis of the tissue [64]. NOS3 (endothelial NOS gene) was also downregulated in the fetuses of the LN group. To date, the effect of IUGR on NOS gene expression in the fetal tissues has not been demonstrated. On the other hand, MUN impairs NO-dependent vasodilation and increases arterial blood pressure in the ovine fetus [65]. Collectively, the present results suggest that downregulation of NOS genes is involved in impairment of the LN fetal muscle development through dysregulated blood flow and/or nutrient supply.
In this study, MUN downregulated genes participating in the HIF-1 signaling pathway (Table 4). The HIF-1 signaling pathway can be promoted by hypoxia-induced factor-α (HIF-1α) in response to inflammatory mediators, as well as hypoxia, and regulates angiogenesis and immunity through activation of genes including NOS2 [66]. HIF-1α expression was not significantly different between the LN and HN groups (p = 0.142), suggesting that downregulation of NOS2 and NOS3 was not mediated by HIF-1α. The mechanism underlying the lower expression of HIF-1 signaling pathway genes in the LN fetuses is unclear; however, the LN fetal muscle may have reduced inflammatory activity as observed in the lower expression associated with inflammatory response (GO:0006954, p = 0.011, data not shown).
The downregulation of NOS genes is assumed to decrease NO synthesis and consequently can affect energy homeostasis through a decrease in the activity of AMPK, as indicated in the role for NO in the GLUT4 regulation via regulation of AMPK [67]. AMPK is a cellular energy sensor that can sense changes in energy metabolism, including glucose uptake and oxidation [68]. In addition, AMPK regulates whole-body energy homeostasis by integrating hormonal and nutritional signals from the cellular environment and the whole organism [69]. As MUN did not alter the AMP/ATP ratio, factors other than cellular energy levels could be involved in the downregulation of AMPK α1 subunit, which may be mediated by NOSs. The NOS3 expression affects glucose uptake via GLUT4 activation in response to insulin [70]. Hence, the downregulation of NOS3 may decrease glucose uptake, which may affect energy homeostasis (Figure 4), in coordination with the downregulated expression of glycolytic genes, such as GPI, ENO3, and PGAM1 (Table 4).
The downregulation of NOS may also affect the muscle angiogenesis in the LN fetuses, as NO is essential for placental angiogenesis and fetal development [71]. This can be coordinated with downregulation of ANGPTL4, a multifunctional factor involved in metabolism and vascular homeostasis [72]. ANGPTL4 promoted NO production through the integrin/JAK/STAT3-mediated upregulation of NOS2 expression in wound epithelia [73]. This is consistent with the results of GO analysis, which revealed that the downregulated genes were associated with angiogenesis and the relevant biological events (Table 4). Furthermore, ANGPTL4 is involved in the vascularization of the tissues in angiogenesis-related disorders [74]. The downregulated expression levels of NOS3 and ANGPTL4 may affect fetal muscle angiogenesis and thereby impair skeletal muscle development.

3.5. Downregulated Genes in the Fetuses of the LN Group Are Associated with Energy Metabolism or Angiogenesis

ANGPTL4 is involved in lipid partitioning and maintenance of lipid homeostasis [75,76], and also has a role in insulin sensitivity and glucose homeostasis [77]. Although ANGPTL4 is considered a cytokine secreted from the adipose tissue and liver, the expression of ANGPTL4 was altered in the biceps femoris muscle of the ovine fetuses after in utero exposure to excess maternal cortisol during late gestation [78]. This suggested that ANGPTL4 regulates the fetal muscle energy metabolism, as well as lipid metabolism, by targeting the distant adipose tissue in response to the nutritional environment. The altered ANGPTL4 expression in the ovine study was accompanied by altered PDK4 expression [78]. PDK4, a major isoform of the pyruvate dehydrogenase kinases in the skeletal and cardiac muscles, is transcriptionally regulated by both glucocorticoids and insulin [79]. Thus, MUN dysregulates energy metabolism and glucose homeostasis partly through the downregulation of ANGPTL4 and PDK4.
Previous studies have demonstrated that IUGR modulates glucose homeostasis in the adipose tissues by regulating APLN (apelin) and its receptor APLNR [80,81] but not in the skeletal muscles so far. Although the expression levels of APLN and APLNR vary depending on the tissues and developmental stages [80], they are broadly expressed in almost all adult tissues, including the skeletal muscle [82]. The apelinergic system, APLN and APLNR, are involved in regulation of glucose uptake [83] and insulin sensitivity in the soleus muscle [84], which can be mediated by AMPK [84]. Thus, the downregulated APLNR expression may contribute to decreased glucose uptake in the LN fetal muscles in this study. However, further studies are needed to elucidate the roles of APLNR in MUN-mediated impaired fetal muscle development.
MUN downregulated the fetal muscle expression of UCP2 (an isoform of UCPs; known as mitochondrial transporters) that is involved in thermogenesis and energy dissipation [85]. In the human skeletal muscle, the downregulated expression of UCP2 was associated with decreased weight loss under calorie restriction conditions [86]. This suggested that the expression of UCP2 is downregulated under energy-depleted conditions. The fetal muscle metabolism in the LN group may shift toward energy-saving pathways by suppressing energy expenditure, such as thermogenesis due to UCP2 activity (Figure 4).
CPT1B, which is involved in the rate-limiting step of mitochondrial β-oxidation, regulates the mitochondrial uptake of long-chain acyl-CoAs [87]. MUN downregulated the CPT-1 activity in the skeletal muscle of offspring of the nutrition-restricted ewes, which may be associated with the increased fat accumulation in the offspring [16]. The findings of this study suggest that MUN downregulated the levels of CPT1B, which may decrease fatty acid oxidation metabolism in the fetuses of the LN group (Figure 4) due to decreased abundance of fatty acids.
RETN and GPNMB were upregulated in the fetuses of the LN group. Previous studies have reported that RETN, an adipokine associated with insulin resistance, is downregulated in the skeletal muscle of rats fed on a high-fat and high-protein diet [88]. The upregulated expression of GPNMB in the blood is correlated with AMPK activation in the skeletal muscle of Zucker diabetic fatty rats [89]. Moreover, GPNMB is involved not only in hypoxic adaptation by regulating blood vessel development and energy metabolism [90] but also in adipose tissue lipogenesis [91]. These genes may contribute to altered energy metabolism in the fetuses of the LN group.

3.6. Altered Expression of Genes Associated with One-Carbon Cycle and Uptake of Large Neutral AAs

Among the LN-induced upregulated genes involved in the target GO and KEGG metabolisms/pathways, the expression levels of 29 genes were examined using qRT-PCR. Of these 29 genes, the qRT-PCR analysis revealed that the expression levels of 26 genes were not significantly different between the LN and HN groups (p > 0.10, data not shown); only three genes showing upregulation in the microarray analysis were validated in qRT-PCR analysis (BDH1, MTHFD2, and ALDH1L2). Thus, the upregulation of genes varied among fetuses of the LN group compared with those of the HN group.
MTHFD2 and ALDH1L2, which are involved in folate-mediated 1C metabolism, were highly upregulated in the fetuses of the LN group (Figure 2). The 1C cycle is the crucial metabolic pathway for purine and thymidine biosynthesis, AA homeostasis, and epigenetic regulation [92]. MSEA revealed that purine metabolism and AA metabolism were activated. Glycine, an important byproduct of 1C metabolism, functions as a precursor for the biosynthesis of metabolites, including glutathione, purine, and heme [92]. The upregulated levels of glycine and MTHFD2 support the activation of AA metabolism, which may be utilized to sustain the LN fetal muscle viability.
SLC7A5, also known as large amino acid transporter 1 (LAT1), was upregulated in the fetuses of the LN group. The uptake of large neutral AAs, such as phenylalanine, leucine, and histidine, is mediated by SLC7A5 [93]. Although the mechanism of SLC7A5 upregulation is unclear, glucose deprivation and low insulin concentration are reported to induce SLC7A5 upregulation in the retina [94] and skeletal muscle [95], respectively. SLC7A5 also functions as an AA sensor and stimulates the mTORC1 pathway [96]. Hence, a global nutrient-restricted environment may upregulate SLC7A5 and consequently upregulate the levels of glutamine, phenylalanine, and histidine in the LN fetal muscle. Further studies are needed to examine the mechanism underlying SLC7A5 upregulation under nutrient-restricted conditions.

3.7. Levels of Nutrients for Dams to Be Compared in the Present Study

The present study lacks animals fed on 100% of the nutrient requirement. In a previous study investigating MUN effect on fetuses in cattle, no significant difference in fetal carcass was observed between the dams fed on 70 and 100% of nutritional requirements during day 30–125 of gestation [31]. This result suggested that a relatively small restriction in maternal nutrition level may not cause a significant phenotypic effect on the fetal carcass. In another study, the fetal whole-body, heart, liver, kidney, and placenta masses were not different between 85 and 140% of the required metabolic energy in the feeding for the dams [30]. No significant phenotypic effect of MUN in the latter study can be explained by the milder nutrient restriction [30]. According to their results, the nutrient level within the range of 85–140% may not be significantly different in terms of MUN effect on fetal muscle. To understand the molecular basis underlying the phenotypic effect of MUN on fetal muscle development, it was necessary to use an experimental design that causes significant phenotypic alteration. Furthermore, a broader range of nutrient levels to be compared was was needed, employing the lowest level as much as possible. The 60% was considered as a minimum level to avoid accidents of pregnant cows [32]. During the last two months of gestation, JFSBC recommends 125% and 141% levels of total digestible nutrients (TDN) and CP, respectively, when the nutrient levels required for maintenance are considered as 100%, for pregnant Japanese Black cows [32,97]. Taking these previous studies into account, we considered that 120% of the requirement was an appropriate option as a counterpart treatment to the low level of 60%. In this experimental condition, significant reduction of the masses of fetal skeletal muscles in the LN group was observed as mentioned above [10]. This enabled us to further investigate molecular mechanism of the MUN effect on fetal muscle development in the present study. Thus, we focused on 60 and 120% of the nutrient requirement to be compared; however, 100% of the requirement might be the best nutritionally standard level and is important as the control level to evaluate the nutritional effect. The MUN effect on fetal muscle development in a case using 60 and 100% levels might be different from that in a case using 60 and 120% levels. The effect of MUN, particularly on fetal skeletal muscle phenotype, needs to be further examined using levels of 60 and 100% of the nutrient requirement.

4. Materials and Methods

4.1. Animals and Dietary Treatments

This study was performed using 11 multiparous Japanese Black cattle (initial BW 488 ± 9.6 kg) at the Iriki farm of Kagoshima University and the farm of Western Region Agricultural Research Center, NARO. The animals were maintained according to the Guide for the Care and Use of Experimental Animals [98]. The experimental protocol was approved by the Kagoshima University Animal Care and Use Committee (#A18007). The management and feeding of animals were performed as described previously [10]. Briefly, the individual diets were designed for pregnant Japanese Black cows to meet 60 or 120% of energy requirement and other nutrients based on the standard diet model for prepregnant BW in JFSBC [32]. Water was supplied sufficiently. The diet comprised the formula feed, total mixed ration, and rice straw as designed previously [10]. According to the JFSC, the recommended contents of dry matter (DM), neutral detergent fiber, acidic detergent fiber, ash, CP, Ca, and p in diet are 68.0, 56.1, 36.0, 11.1, 8.0, 0.6, and 0.3%, respectively (all values are expressed as a percentage of DM). The cows randomly assigned to the LN (n = 5) and HN (n = 6) diet groups were fed on the respective diet during gestation including prepregnant period. Metabolizable energy of the mixture feed for 100% requirement was estimated as 8.56 MJ/kg DM. The cows were subjected to artificial insemination (AI) using male-sorted semen of an identical sire. The pregnancy was diagnosed using trans-rectal ultrasonography (HS-1500; Honda Electronics Co., Ltd., Toyohashi, Japan) at approximately day 40 post-AI.

4.2. Sample Collection

In total, 11 fetuses were obtained from the cows through cesarean section at the Kagoshima University Veterinary Teaching Hospital. The fetuses were euthanized by exsanguination at day 260 ± 8.3 of gestation, after injecting 200 mg of lidocaine (AstraZeneca, Osaka, Japan) into the jugular vein. Fetal BW and body length were recorded. The fetuses were dissected to collect the skeletal muscles and adipose tissues for subsequent analysis. The LT muscle was collected from the right side of the carcass, frozen with liquid nitrogen, or soaked in RNAlater® (Thermo Fisher Scientific, Tokyo, Japan), and stored at −80 °C for further analyses.

4.3. Sample Preparation for CE-TOFMS

Among the 11 fetuses, the fetuses with the lowest BW in the LN group (n = 4) and those with the highest BW in the HN group (n = 4) were selected. The LT muscle samples were subjected to metabolomic analysis using CE-TOFMS. The frozen muscle pieces (46.3–90.0 mg) were immediately immersed into a solution of 50% acetonitrile and 10 μM internal standard solution 1 (Human Metabolome Technologies, Tsuruoka, Japan) at 0 °C and homogenized twice at 1500 rpm for 120 s. The samples were centrifuged at 2300× g and 4 °C for 5 min. The upper layer solution was filtered through a Millipore 5-kDa cutoff membrane using centrifugation. The filtrate was lyophilized, suspended in Milli-Q water, and analyzed using CE-TOFMS.

4.4. Instrumentation and Conditions of CE-TOFMS

CE-TOFMS was performed using an Agilent capillary electrophoresis system equipped with an Agilent 6210 time-of-flight mass spectrometer, Agilent 1100 isocratic high-performance liquid chromatography pump, Agilent G1603A CE-MS adapter kit, and Agilent G1607A CE-ESI-MS sprayer kit (Agilent Technologies, Waldbronn, Germany). The analytic conditions were identical to those used in a previous study [99]. The spectrometer was scanned from m/z 50 to 1000.

4.5. Data Analysis of CE-TOFMS Results

The raw data of CE-TOFMS were processed with MasterHands as described in a previous study [99]. Among the detected compounds, the compounds annotated in the Human Metabolome Database (ver. 4.0, http://www.hmdb.ca/, accessed on 17 March 2021) or KEGG database (http://www.genome.jp/kegg/, accessed on 17 March 2021) were further analyzed. The relative contents of the annotated compounds were determined by comparing the peaks of compounds with the same MS properties in the analysis. To compare the relative contents of the compounds between the LN and HN groups, the peak areas were normalized to those of the internal standards (methionine sulfone for cations, (+)-camphor-10-sulfonic acid for anions) as well as to sample weight. The contents of major metabolites, such as glycolytic products, AAs, and ATP degradation products were determined using commercially available standards. The abundance of each compound used for comparative analysis was set as 0 when the level of the compound was not detected. File conversion of raw MS data, peak picking, noise reduction, and alignment of data for multiple samples were conducted as previously described [99].

4.6. RNA Preparation and Complementary DNA (cDNA) Synthesis

Total RNA was extracted from the LT muscle using the mirVanaTM microRNA isolation kit (Thermo Fisher Scientific, Tokyo, Japan), following the manufacturer’s instructions. The quantity and quality of the RNA were determined using an Agilent Bioanalyzer 2100 with an RNA 6000 Pico kit (Agilent Technologies, Santa Clara, CA, USA). For polymerase chain reaction (PCR) analysis, total RNA was prepared using ISOGEN II (Nippon Gene, Toyama, Japan). The cDNAs were synthesized from 1000 ng of total RNA using the ReverTra Ace qPCR RT kit (Toyobo), following the manufacturer’s instructions.

4.7. Microarray Analysis

Similar to metabolomic analysis, the microarray analysis was performed using the muscle samples from fetuses with the lowest BW in the LN group (n = 4) and those with the highest BW in the HN group (n = 4). The total RNA samples of four steers from the HN and LN groups were pooled and applied to a Bovine (v2) Gene Expression 4x44K Microarray (Agilent). The signals of the hybridized probes were detected using an Agilent Microarray Scanner (Agilent). The results were normalized by quantile method using GeneSpring GX (Agilent). Array data were deposited in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database, and are accessible through GEO Series accession number GSE176377 (http://www.ncbi.nlm.nih.gov/geo, accessed on 8 June 2021).

4.8. Quantitative Real-Time PCR (qRT-PCR) Analysis

The cDNA samples of the LN (n = 5) and HN (n = 6) groups were subjected to qRT-PCR analysis. The qRT-PCR analysis was performed using CFX96 thermal cycler (Bio-Rad, Hercules, CA, USA) with the QuantiTect SYBR Green PCR kit (Qiagen, Tokyo, Japan) and the primers listed in Table S2. Ribosomal protein L7 (RPL7) was used as an internal control. Melting curve analysis was performed to confirm the specificity of the amplification reactions.

4.9. Functional Annotation of Target Genes

To classify the genes of interest according to their functional annotation, GO and pathway analyses were performed on the differentially expressed genes, which were determined from qRT-PCR analysis, between the HN and LN groups. The genes of interest were analyzed using Database for Annotation, Visualization, and Integrated Discovery (version 6.7, http://david.abcc.ncifcrf.gov, accessed on 2 February 2021) with the setting of Bos taurus as the background species to enrich GO terms and characterize the KEGG pathway terms defined by KEGG (http://www.genome.jp/kegg/, accessed on 17 March 2021) for the respective biological process associated with the effect of MUN. The terms were considered significant at p < 0.05.

4.10. Statistical Analyses

The effects of LN and HN diets on the metabolite and gene expression profiles were determined. The data were analyzed using the two-sided Student’s t-test for metabolomics data or the one-sided Student’s t-test for PCR data based on the trend of gene expression in microarray analysis. The differences were considered significant at p < 0.05 or a trend at p < 0.10. HCA, and MSEA using KEGG as a metabolite set library were performed using MetaboAnalyst 5.0 (https://www.metaboanalyst.ca/MetaboAnalyst/faces/home.xhtml, accessed on 26 May 2021).

5. Conclusions

This study examined the effect of MUN on the fetal LT muscle development at the late gestational stage, the levels of AAs and other metabolites, and the expression of genes associated with energy metabolism and angiogenesis. The levels of glutamine and its related metabolites were increased, which was associated with decreased muscle mass in the fetuses of the LN group. This suggested that these AAs accumulate due to the downregulation of protein synthesis rather than the activation of protein degradation. AA accumulation may be due to restricted energy caused by the downregulation of energy homeostasis-associated genes, including ANGPTL4, APLNR, NOS2, and glycolytic genes, at the expense of activation of metabolisms such as fatty acids and thermogenesis. MUN also downregulated the angiogenesis-associated genes, such as ANGPT4, ANGPTL4, and NOS3 that are important for fetal muscle development through regulation of local blood flow. The findings of this study indicate that the metabolism of AAs and the expression of genes associated with energy expenditure, glucose homeostasis, and angiogenesis were altered in the MUN fetal muscle for prioritizing cell survival over protein accumulation and muscle growth.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/metabo11090582/s1, Table S1: Relative levels of detected metabolites in LN and HN fetal LT muscle by CE-TOFMS metabolomics, Table S2: PCR primers used in this study.

Author Contributions

S.M. and T.G. conceived and designed the research; Y.Z., K.O. (Kounosuke Otomaru), K.O. (Kazunaga Oshima), Y.G., I.O. and T.G. performed the trials of maternal nutrition restriction; Y.Z., K.O. (Koichi Ojima), M.S., A.K., S.R. and T.G. performed collection and analysis of carcass data; S.M., M.O. and K.O. (Koichi Ojima) conducted metabolomics, microarray, and PCR analyses; S.M. analyzed the metabolism and transcriptomics data; S.M. drafted the manuscript; S.M. and T.G. edited and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Japan Society for the Promotion of Science (JSPS KAKENHI JP18K05960, JP26310312, and JP19KT0013), and the Canon Fund (R15-0089).

Institutional Review Board Statement

All experimental procedures followed the Guide for the Care and Use of Experimental Animals (the Kagoshima University Animal Care and Use Committee), which was approved by the committee (#A18007). There was no use of human participants, data, or tissues.

Informed Consent Statement

Not applicable.

Data Availability Statement

Array data were deposited in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database, and are accessible through GEO Series accession number GSE176377 (http://www.ncbi.nlm.nih.gov/geo, accessed on 8 June 2021).

Acknowledgments

We thank the personnel at the farm staff of the Kagoshima University and Western Region Agricultural Research Center, NARO for their great support to help manage and feed cattle. We especially thank Okamura, R. Saneshima, Y. Nagao (Kagoshima Univ.), M. Futohashi (Tohoku Univ.), C. Shindo, M. Ichimura, and Y. Eguchi (NILGS) for their technical supports.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Sandoval, C.; Lambo, C.A.; Beason, K.; Dunlap, K.A.; Satterfield, M.C. Effect of maternal nutrient restriction on skeletal muscle mass and associated molecular pathways in SGA and Non-SGA sheep fetuses. Domest. Anim. Endocrinol. 2020, 72, 106443. [Google Scholar] [CrossRef]
  2. Yates, D.T.; Macko, A.R.; Nearing, M.; Chen, X.; Rhoads, R.P.; Limesand, S.W. Developmental programming in response to intrauterine growth restriction impairs myoblast function and skeletal muscle metabolism. J. Pregnancy 2012, 2012, 631038. [Google Scholar] [CrossRef]
  3. Osgerby, J.C.; Wathes, D.C.; Howard, D.; Gadd, T.S. The effect of maternal undernutrition on ovine fetal growth. J. Endocrinol. 2002, 173, 131–141. [Google Scholar] [CrossRef] [Green Version]
  4. Vonnahme, K.A.; Hess, B.W.; Hansen, T.R.; McCormick, R.J.; Rule, D.C.; Moss, G.E.; Murdoch, W.J.; Nijland, M.J.; Skinner, D.C.; Nathanielsz, P.W.; et al. Maternal undernutrition from early- to mid-gestation leads to growth retardation, cardiac ventricular hypertrophy, and increased liver weight in the fetal sheep. Biol. Reprod. 2003, 69, 133–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. George, L.A.; Zhang, L.; Tuersunjiang, N.; Ma, Y.; Long, N.M.; Uthlaut, A.B.; Smith, D.T.; Nathanielsz, P.W.; Ford, S.P. Early maternal undernutrition programs increased feed intake, altered glucose metabolism and insulin secretion, and liver function in aged female offspring. Am. J. Physiol. Regul. Integr. Comp. Physiol. 2012, 302, R795–R804. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Lloyd, L.J.; Foster, T.; Rhodes, P.; Rhind, S.M.; Gardner, D.S. Protein-energy malnutrition during early gestation in sheep blunts fetal renal vascular and nephron development and compromises adult renal function. J. Physiol. 2012, 590, 377–393. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Shukla, P.; Ghatta, S.; Dubey, N.; Lemley, C.O.; Johnson, M.L.; Modgil, A.; Vonnahme, K.; Caton, J.S.; Reynolds, L.P.; Sun, C.; et al. Maternal nutrient restriction during pregnancy impairs an endothelium-derived hyperpolarizing factor-like pathway in sheep fetal coronary arteries. Am. J. Physiol. Heart Circ. Physiol. 2014, 307, H134–H142. [Google Scholar] [CrossRef] [Green Version]
  8. Devaskar, S.U.; Chu, A. Intrauterine Growth Restriction: Hungry for an Answer. Physiology 2016, 31, 131–146. [Google Scholar] [CrossRef] [Green Version]
  9. Brown, L.D. Endocrine regulation of fetal skeletal muscle growth: Impact on future metabolic health. J. Endocrinol. 2014, 221, R13–R29. [Google Scholar] [CrossRef] [Green Version]
  10. Zhang, Y.; Otomaru, K.; Oshima, K.; Goto, Y.; Oshima, I.; Muroya, S.; Sano, M.; Saneshima, R.; Nagao, Y.; Kinoshita, A.; et al. Effects of low and high levels of maternal nutrition consumed for the entirety of gestation on the development of muscle, adipose tissue, bone, and the organs of Wagyu cattle fetuses. Anim. Sci. J. 2021, 92, e13600. [Google Scholar] [CrossRef]
  11. Zhu, M.J.; Ford, S.P.; Nathanielsz, P.W.; Du, M. Effect of maternal nutrient restriction in sheep on the development of fetal skeletal muscle. Biol. Reprod. 2004, 71, 1968–1973. [Google Scholar] [CrossRef] [Green Version]
  12. Quigley, S.P.; Kleemann, D.O.; Kakar, M.A.; Owens, J.A.; Nattrass, G.S.; Maddocks, S.; Walker, S.K. Myogenesis in sheep is altered by maternal feed intake during the peri-conception period. Anim. Reprod. Sci. 2005, 87, 241–251. [Google Scholar] [CrossRef] [PubMed]
  13. Widdowson, E.M.; Crabb, D.E.; Milner, R.D. Cellular development of some human organs before birth. Arch. Dis. Child. 1972, 47, 652–655. [Google Scholar] [CrossRef] [Green Version]
  14. Yates, D.T.; Clarke, D.S.; Macko, A.R.; Anderson, M.J.; Shelton, L.A.; Nearing, M.; Allen, R.E.; Rhoads, R.P.; Limesand, S.W. Myoblasts from intrauterine growth-restricted sheep fetuses exhibit intrinsic deficiencies in proliferation that contribute to smaller semitendinosus myofibres. J. Physiol. 2014, 592, 3113–3125. [Google Scholar] [CrossRef]
  15. Maier, A.; McEwan, J.C.; Dodds, K.G.; Fischman, D.A.; Fitzsimons, R.B.; Harris, A.J. Myosin heavy chain composition of single fibres and their origins and distribution in developing fascicles of sheep tibialis cranialis muscles. J. Muscle Res. Cell Motil. 1992, 13, 551–572. [Google Scholar] [CrossRef]
  16. Zhu, M.J.; Ford, S.P.; Means, W.J.; Hess, B.W.; Nathanielsz, P.W.; Du, M. Maternal nutrient restriction affects properties of skeletal muscle in offspring. J. Physiol. 2006, 575, 241–250. [Google Scholar] [CrossRef] [PubMed]
  17. Yates, D.T.; Cadaret, C.N.; Beede, K.A.; Riley, H.E.; Macko, A.R.; Anderson, M.J.; Camacho, L.E.; Limesand, S.W. Intrauterine growth-restricted sheep fetuses exhibit smaller hindlimb muscle fibers and lower proportions of insulin-sensitive Type I fibers near term. Am. J. Physiol. Regul. Integr. Comp. Physiol. 2016, 310, R1020–R1029. [Google Scholar] [CrossRef] [Green Version]
  18. Maresca, S.; Valiente, S.L.; Rodriguez, A.M.; Testa, L.M.; Long, N.M.; Quintans, G.I.; Pavan, E. The influence of protein restriction during mid- to late gestation on beef offspring growth, carcass characteristic and meat quality. Meat Sci. 2019, 153, 103–108. [Google Scholar] [CrossRef]
  19. Ramírez, M.; Testa, L.M.; López Valiente, S.; Latorre, M.E.; Long, N.M.; Rodriguez, A.M.; Pavan, E.; Maresca, S. Maternal energy status during late gestation: Effects on growth performance, carcass characteristics and meat quality of steers progeny. Meat Sci. 2020, 164, 108095. [Google Scholar] [CrossRef]
  20. Wu, G.; Bazer, F.W.; Burghardt, R.C.; Johnson, G.A.; Kim, S.W.; Li, X.L.; Satterfield, M.C.; Spencer, T.E. Impacts of amino acid nutrition on pregnancy outcome in pigs: Mechanisms and implications for swine production. J. Anim. Sci. 2010, 88, E195–E204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Martin, D.E.; Jones, A.K.; Pillai, S.M.; Hoffman, M.L.; McFadden, K.K.; Zinn, S.A.; Govoni, K.E.; Reed, S.A. Maternal Restricted- and Over-Feeding During Gestation Result in Distinct Lipid and Amino Acid Metabolite Profiles in the Longissimus Muscle of the Offspring. Front. Physiol. 2019, 10, 515. [Google Scholar] [CrossRef] [Green Version]
  22. Chang, E.I.; Wesolowski, S.R.; Gilje, E.A.; Baker, P.R., 2nd; Reisz, J.A.; D’Alessandro, A.; Hay, W.W., Jr.; Rozance, P.J.; Brown, L.D. Skeletal muscle amino acid uptake is lower and alanine production is greater in late gestation intrauterine growth-restricted fetal sheep hindlimb. Am. J. Physiol. Regul. Integr. Comp. Physiol. 2019, 317, R615–R629. [Google Scholar] [CrossRef]
  23. Dunlop, K.; Cedrone, M.; Staples, J.F.; Regnault, T.R. Altered fetal skeletal muscle nutrient metabolism following an adverse in utero environment and the modulation of later life insulin sensitivity. Nutrients 2015, 7, 1202–1216. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. 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, skz388. [Google Scholar] [CrossRef] [PubMed]
  25. He, J.; He, Y.; Yu, B.; Wang, X.; Chen, D. Transcriptome Characterization of Repressed Embryonic Myogenesis Due to Maternal Calorie Restriction. Front. Cell Dev. Biol. 2020, 8, 527. [Google Scholar] [CrossRef] [PubMed]
  26. Sanglard, L.P.; Nascimento, M.; Moriel, P.; Sommer, J.; Ashwell, M.; Poore, M.H.; Duarte, M.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, 702. [Google Scholar] [CrossRef]
  27. Sohel, M.M.H.; Akyuz, B.; Konca, Y.; Arslan, K.; Gurbulak, K.; Abay, M.; Kaliber, M.; Cinar, M.U. Differential protein input in the maternal diet alters the skeletal muscle transcriptome in fetal sheep. Mamm. Genome 2020, 31, 309–324. [Google Scholar] [CrossRef]
  28. Diniz, W.J.S.; Crouse, M.S.; Cushman, R.A.; McLean, K.J.; Caton, J.S.; Dahlen, C.R.; Reynolds, L.P.; Ward, A.K. Cerebrum, liver, and muscle regulatory networks uncover maternal nutrition effects in developmental programming of beef cattle during early pregnancy. Sci. Rep. 2021, 11, 2771. [Google Scholar] [CrossRef]
  29. Sandoval, C.; Wu, G.; Smith, S.B.; Dunlap, K.A.; Satterfield, M.C. Maternal Nutrient Restriction and Skeletal Muscle Development: Consequences for Postnatal Health. Adv. Exp. Med. Biol. 2020, 1265, 153–165. [Google Scholar] [CrossRef]
  30. 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, 632. [Google Scholar] [CrossRef] [PubMed]
  31. 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]
  32. National Agriculture and Food Research Organization. Japanese Feeding Standard for Beef Cattle 2008 Edn; Japan Livestock Industry Association: Tokyo, Japan, 2009. (In Japanese)
  33. Wu, G.; Bazer, F.W.; Wallace, J.M.; Spencer, T.E. Board-invited review: Intrauterine growth retardation: Implications for the animal sciences. J. Anim. Sci. 2006, 84, 2316–2337. [Google Scholar] [CrossRef]
  34. Rozance, P.J.; Zastoupil, L.; Wesolowski, S.R.; Goldstrohm, D.A.; Strahan, B.; Cree-Green, M.; Sheffield-Moore, M.; Meschia, G.; Hay, W.W., Jr.; Wilkening, R.B.; et al. Skeletal muscle protein accretion rates and hindlimb growth are reduced in late gestation intrauterine growth-restricted fetal sheep. J. Physiol. 2018, 596, 67–82. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Du, M.; Zhu, M.J.; Means, W.J.; Hess, B.W.; Ford, S.P. Effect of nutrient restriction on calpain and calpastatin content of skeletal muscle from cows and fetuses1. J. Anim. Sci. 2004, 82, 2541–2547. [Google Scholar] [CrossRef] [PubMed]
  36. Brosnan, J.T. Interorgan amino acid transport and its regulation. J. Nutr. 2003, 133, 2068s–2072s. [Google Scholar] [CrossRef] [Green Version]
  37. Aschenbach, J.R.; Kristensen, N.B.; Donkin, S.S.; Hammon, H.M.; Penner, G.B. Gluconeogenesis in dairy cows: The secret of making sweet milk from sour dough. IUBMB Life 2010, 62, 869–877. [Google Scholar] [CrossRef]
  38. Morris, C.R.; Hamilton-Reeves, J.; Martindale, R.G.; Sarav, M.; Ochoa Gautier, J.B. Acquired Amino Acid Deficiencies: A Focus on Arginine and Glutamine. Nutr. Clin. Pract. 2017, 32, 30S–47S. [Google Scholar] [CrossRef]
  39. Yoo, H.C.; Yu, Y.C.; Sung, Y.; Han, J.M. Glutamine reliance in cell metabolism. Exp. Mol. Med. 2020, 52, 1496–1516. [Google Scholar] [CrossRef] [PubMed]
  40. Leury, B.J.; Bird, A.R.; Chandler, K.D.; Bell, A.W. Glucose partitioning in the pregnant ewe: Effects of undernutrition and exercise. Br. J. Nutr. 1990, 64, 449–462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Lemons, J.A.; Schreiner, R.L. Amino acid metabolism in the ovine fetus. Am. J. Physiol. 1983, 244, E459–E466. [Google Scholar] [CrossRef] [PubMed]
  42. Faichney, G.J.; White, G.A. Effects of maternal nutritional status on fetal and placental growth and on fetal urea synthesis in sheep. Aust. J. Biol. Sci. 1987, 40, 365–377. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Morris, S.M., Jr. Arginine Metabolism Revisited. J. Nutr. 2016, 146, 2579S–2586S. [Google Scholar] [CrossRef] [PubMed]
  44. Cervelli, M.; Leonetti, A.; Duranti, G.; Sabatini, S.; Ceci, R.; Mariottini, P. Skeletal Muscle Pathophysiology: The Emerging Role of Spermine Oxidase and Spermidine. Med. Sci. 2018, 6, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Hosten, A.O. BUN and Creatinine. In Clinical Methods: The History, Physical, and Laboratory Examinations; Walker, H.K., Hall, W.D., Hurst, J.W., Eds.; Butterworths: Boston, MA, USA, 1990; pp. 874–878. [Google Scholar]
  46. Bessman, S.P.; Geiger, P.J. Transport of energy in muscle: The phosphorylcreatine shuttle. Science 1981, 211, 448–452. [Google Scholar] [CrossRef] [PubMed]
  47. Dessì, A.; Atzori, L.; Noto, A.; Visser, G.H.; Gazzolo, D.; Zanardo, V.; Barberini, L.; Puddu, M.; Ottonello, G.; Atzei, A.; et al. Metabolomics in newborns with intrauterine growth retardation (IUGR): Urine reveals markers of metabolic syndrome. J. Matern. Fetal Neonatal Med. 2011, 24 (Suppl. 2), 35–39. [Google Scholar] [CrossRef]
  48. Zorzano, A.; Lasunción, M.A.; Herrera, E. Role of the availability of substrates on hepatic and renal gluconeogenesis in the fasted late pregnant rat. Metabolism 1986, 35, 297–303. [Google Scholar] [CrossRef]
  49. Zeng, Z.; Liu, F.; Li, S. Metabolic Adaptations in Pregnancy: A Review. Ann. Nutr. Metab. 2017, 70, 59–65. [Google Scholar] [CrossRef]
  50. Bevilacqua, A.; Bizzarri, M. Inositols in Insulin Signaling and Glucose Metabolism. Int. J. Endocrinol. 2018, 2018, 1968450. [Google Scholar] [CrossRef] [Green Version]
  51. Huang, L.C.; Fonteles, M.C.; Houston, D.B.; Zhang, C.; Larner, J. Chiroinositol deficiency and insulin resistance. III. Acute glycogenic and hypoglycemic effects of two inositol phosphoglycan insulin mediators in normal and streptozotocin-diabetic rats in vivo. Endocrinology 1993, 132, 652–657. [Google Scholar] [CrossRef]
  52. Larner, J. D-chiro-inositol-its functional role in insulin action and its deficit in insulin resistance. Int. J. Exp. Diabetes Res. 2002, 3, 47–60. [Google Scholar] [CrossRef]
  53. Regnault, T.R.; Teng, C.; de Vrijer, B.; Galan, H.L.; Wilkening, R.B.; Battaglia, F.C. The tissue and plasma concentration of polyols and sugars in sheep intrauterine growth retardation. Exp. Biol. Med. 2010, 235, 999–1006. [Google Scholar] [CrossRef] [PubMed]
  54. Nissen, P.M.; Nebel, C.; Oksbjerg, N.; Bertram, H.C. Metabolomics reveals relationship between plasma inositols and birth weight: Possible markers for fetal programming of type 2 diabetes. J. Biomed. Biotechnol. 2011, 2011, 378268. [Google Scholar] [CrossRef] [PubMed]
  55. Metges, C.C.; Görs, S.; Lang, I.S.; Hammon, H.M.; Brüssow, K.P.; Weitzel, J.M.; Nürnberg, G.; Rehfeldt, C.; Otten, W. Low and high dietary protein: Carbohydrate ratios during pregnancy affect materno-fetal glucose metabolism in pigs. J. Nutr. 2014, 144, 155–163. [Google Scholar] [CrossRef]
  56. Priante, E.; Verlato, G.; Giordano, G.; Stocchero, M.; Visentin, S.; Mardegan, V.; Baraldi, E. Intrauterine Growth Restriction: New Insight from the Metabolomic Approach. Metabolites 2019, 9, 267. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Shimizu, H.; Kakimoto, Y.; Nakajima, T.; Kanazawa, A.; Sano, I. Isolation and identification of 2-aminoethylphosphonic acid from bovine brain. Nature 1965, 207, 1197–1198. [Google Scholar] [CrossRef] [PubMed]
  58. Alhadeff, J.A.; Van Bruggen, J.T.; Doyle Daves, G. Biosynthetic studies on 2-aminoethylphosphonic acid in a mammalian (rat) system. Biochim. Biophys. Acta (BBA)—Gen. Subj. 1972, 286, 103–106. [Google Scholar] [CrossRef]
  59. Hasegawa, S.; Tamari, M.; Kametaka, M. Isolation of diacylglyceryl-2-aminoethylphosphonate from bovine liver. J. Biochem. 1976, 80, 531–535. [Google Scholar] [CrossRef]
  60. Tan, S.A.; Tan, L.G. Distribution of ciliatine (2-aminoethylphosphonic acid) and phosphonoalanine (2-amino-3-phosphonopropionic acid) in human tissues. Clin. Physiol. Biochem. 1989, 7, 303–309. [Google Scholar] [PubMed]
  61. Joseph, J.C.; Henderson, T.O. 2-Aminoethylphosphonic acid metabolism in the rat. Lipids 1977, 12, 75–84. [Google Scholar] [CrossRef]
  62. Nagata, Y.; Partridge, T.A.; Matsuda, R.; Zammit, P.S. Entry of muscle satellite cells into the cell cycle requires sphingolipid signaling. J. Cell Biol. 2006, 174, 245–253. [Google Scholar] [CrossRef] [Green Version]
  63. Stamler, J.S.; Meissner, G. Physiology of nitric oxide in skeletal muscle. Physiol. Rev. 2001, 81, 209–237. [Google Scholar] [CrossRef] [PubMed]
  64. Williams, J.L.; Cartland, D.; Hussain, A.; Egginton, S. A differential role for nitric oxide in two forms of physiological angiogenesis in mouse. J. Physiol. 2006, 570, 445–454. [Google Scholar] [CrossRef] [PubMed]
  65. Wu, G.; Bazer, F.W.; Cudd, T.A.; Meininger, C.J.; Spencer, T.E. Maternal nutrition and fetal development. J. Nutr. 2004, 134, 2169–2172. [Google Scholar] [CrossRef]
  66. Hellwig-Bürgel, T.; Stiehl, D.P.; Wagner, A.E.; Metzen, E.; Jelkmann, W. Review: Hypoxia-inducible factor-1 (HIF-1): A novel transcription factor in immune reactions. J. Interferon Cytokine Res. 2005, 25, 297–310. [Google Scholar] [CrossRef]
  67. Lira, V.A.; Soltow, Q.A.; Long, J.H.; Betters, J.L.; Sellman, J.E.; Criswell, D.S. Nitric oxide increases GLUT4 expression and regulates AMPK signaling in skeletal muscle. Am. J. Physiol. Endocrinol. Metab. 2007, 293, E1062–E1068. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Jobgen, W.S.; Fried, S.K.; Fu, W.J.; Meininger, C.J.; Wu, G. Regulatory role for the arginine-nitric oxide pathway in metabolism of energy substrates. J. Nutr. Biochem. 2006, 17, 571–588. [Google Scholar] [CrossRef] [PubMed]
  69. Kjøbsted, R.; Hingst, J.R.; Fentz, J.; Foretz, M.; Sanz, M.-N.; Pehmøller, C.; Shum, M.; Marette, A.; Mounier, R.; Treebak, J.T.; et al. AMPK in skeletal muscle function and metabolism. FASEB J. 2018, 32, 1741–1777. [Google Scholar] [CrossRef] [Green Version]
  70. Bahadoran, Z.; Mirmiran, P.; Ghasemi, A. Role of Nitric Oxide in Insulin Secretion and Glucose Metabolism. Trends Endocrinol. Metab. TEM 2020, 31, 118–130. [Google Scholar] [CrossRef] [PubMed]
  71. Wu, G.; Bazer, F.W.; Satterfield, M.C.; Li, X.; Wang, X.; Johnson, G.A.; Burghardt, R.C.; Dai, Z.; Wang, J.; Wu, Z. Impacts of arginine nutrition on embryonic and fetal development in mammals. Amino Acids 2013, 45, 241–256. [Google Scholar] [CrossRef] [PubMed]
  72. Fernández-Hernando, C.; Suárez, Y. ANGPTL4: A multifunctional protein involved in metabolism and vascular homeostasis. Curr. Opin. Hematol. 2020, 27, 206–213. [Google Scholar] [CrossRef]
  73. Chong, H.C.; Chan, J.S.; Goh, C.Q.; Gounko, N.V.; Luo, B.; Wang, X.; Foo, S.; Wong, M.T.; Choong, C.; Kersten, S.; et al. Angiopoietin-like 4 stimulates STAT3-mediated iNOS expression and enhances angiogenesis to accelerate wound healing in diabetic mice. Mol. Ther. 2014, 22, 1593–1604. [Google Scholar] [CrossRef] [Green Version]
  74. Trayhurn, P.; Wood, I.S. Adipokines: Inflammation and the pleiotropic role of white adipose tissue. Br. J. Nutr. 2004, 92, 347–355. [Google Scholar] [CrossRef] [Green Version]
  75. Aryal, B.; Price, N.L.; Suarez, Y.; Fernández-Hernando, C. ANGPTL4 in Metabolic and Cardiovascular Disease. Trends Mol. Med. 2019, 25, 723–734. [Google Scholar] [CrossRef]
  76. Muroya, S.; Ogasawara, H.; Nohara, K.; Oe, M.; Ojima, K.; Hojito, M. Coordinated alteration of mRNA-microRNA transcriptomes associated with exosomes and fatty acid metabolism in adipose tissue and skeletal muscle in grazing cattle. Asian-Australas. J. Anim. Sci. 2020, 33, 1824–1836. [Google Scholar] [CrossRef] [PubMed]
  77. La Paglia, L.; Listì, A.; Caruso, S.; Amodeo, V.; Passiglia, F.; Bazan, V.; Fanale, D. Potential Role of ANGPTL4 in the Cross Talk between Metabolism and Cancer through PPAR Signaling Pathway. PPAR Res. 2017, 2017, 8187235. [Google Scholar] [CrossRef]
  78. Joseph, S.; Alava, B.; Antolic, A.; Richards, E.M.; Wood, C.E.; Keller-Wood, M. Fetal ovine skeletal and cardiac muscle transcriptomics are differentially altered by increased maternal cortisol during gestation. Physiol. Genom. 2020, 52, 178–190. [Google Scholar] [CrossRef] [PubMed]
  79. Zhang, S.; Hulver, M.W.; McMillan, R.P.; Cline, M.A.; Gilbert, E.R. The pivotal role of pyruvate dehydrogenase kinases in metabolic flexibility. Nutr. Metab. 2014, 11, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Mayeur, S.; Wattez, J.S.; Lukaszewski, M.A.; Lecoutre, S.; Butruille, L.; Drougard, A.; Eberlé, D.; Bastide, B.; Laborie, C.; Storme, L.; et al. Apelin Controls Fetal and Neonatal Glucose Homeostasis and Is Altered by Maternal Undernutrition. Diabetes 2016, 65, 554–560. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Lecoutre, S.; Marousez, L.; Drougard, A.; Knauf, C.; Guinez, C.; Eberlé, D.; Laborie, C.; Vieau, D.; Lesage, J.; Breton, C. Maternal undernutrition programs the apelinergic system of adipose tissue in adult male rat offspring. J. Dev. Orig. Health Dis. 2017, 8, 3–7. [Google Scholar] [CrossRef]
  82. Shin, K.; Kenward, C.; Rainey, J.K. Apelinergic System Structure and Function. Compr. Physiol. 2017, 8, 407–450. [Google Scholar] [CrossRef]
  83. Dray, C.; Knauf, C.; Daviaud, D.; Waget, A.; Boucher, J.; Buléon, M.; Cani, P.D.; Attané, C.; Guigné, C.; Carpéné, C.; et al. Apelin stimulates glucose utilization in normal and obese insulin-resistant mice. Cell Metab. 2008, 8, 437–445. [Google Scholar] [CrossRef] [PubMed]
  84. Yue, P.; Jin, H.; Aillaud, M.; Deng, A.C.; Azuma, J.; Asagami, T.; Kundu, R.K.; Reaven, G.M.; Quertermous, T.; Tsao, P.S. Apelin is necessary for the maintenance of insulin sensitivity. Am. J. Physiol. Endocrinol. Metab. 2010, 298, E59–E67. [Google Scholar] [CrossRef] [Green Version]
  85. Pratt, J.; Boreham, C.; Ennis, S.; Ryan, A.W.; De Vito, G. Genetic Associations with Aging Muscle: A Systematic Review. Cells 2019, 9, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Heinitz, S.; Piaggi, P.; Yang, S.; Bonfiglio, S.; Steel, J.; Krakoff, J.; Votruba, S.B. Response of skeletal muscle UCP2-expression during metabolic adaptation to caloric restriction. Int. J. Obes. 2018, 42, 974–984. [Google Scholar] [CrossRef]
  87. Bonnefont, J.P.; Djouadi, F.; Prip-Buus, C.; Gobin, S.; Munnich, A.; Bastin, J. Carnitine palmitoyltransferases 1 and 2: Biochemical, molecular and medical aspects. Mol. Asp. Med. 2004, 25, 495–520. [Google Scholar] [CrossRef] [PubMed]
  88. Nowacka-Woszuk, J.; Pruszynska-Oszmalek, E.; Szydlowski, M.; Sadkowski, S.; Szczerbal, I. Diet-induced variability of the resistin gene (Retn) transcript level and methylation profile in rats. BMC Genet. 2015, 16, 113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Grempler, R.; Wolff, M.; Simon, E.; Schmid, R.; Eisele, C.; Rieber, K.; Fischer, E.; Mettel, S.; Gabrielyan, O.; Delic, D.; et al. Discovery and translation of a target engagement marker for AMP-activated protein kinase (AMPK). PLoS ONE 2018, 13, e0197849. [Google Scholar] [CrossRef] [PubMed]
  90. Zhang, Y.; Zhang, H.; Zhang, B.; Ling, Y.; Zhang, H. Identification of key HIF-1α target genes that regulate adaptation to hypoxic conditions in Tibetan chicken embryos. Gene 2020, 729, 144321. [Google Scholar] [CrossRef]
  91. Gong, X.M.; Li, Y.F.; Luo, J.; Wang, J.Q.; Wei, J.; Wang, J.Q.; Xiao, T.; Xie, C.; Hong, J.; Ning, G.; et al. Gpnmb secreted from liver promotes lipogenesis in white adipose tissue and aggravates obesity and insulin resistance. Nat. Metab. 2019, 1, 570–583. [Google Scholar] [CrossRef]
  92. Ducker, G.S.; Rabinowitz, J.D. One-Carbon Metabolism in Health and Disease. Cell Metab. 2017, 25, 27–42. [Google Scholar] [CrossRef] [Green Version]
  93. Scalise, M.; Galluccio, M.; Console, L.; Pochini, L.; Indiveri, C. The Human SLC7A5 (LAT1): The Intriguing Histidine/Large Neutral Amino Acid Transporter and Its Relevance to Human Health. Front. Chem. 2018, 6, 243. [Google Scholar] [CrossRef] [PubMed]
  94. Matsuyama, R.; Tomi, M.; Akanuma, S.-i.; Tabuchi, A.; Kubo, Y.; Tachikawa, M.; Hosoya, K.-i. Up-regulation of l-Type Amino Acid Transporter 1 (LAT1) in Cultured Rat Retinal Capillary Endothelial Cells in Response to Glucose Deprivation. Drug Metab. Pharmacokinet. 2012, 27, 317–324. [Google Scholar] [CrossRef]
  95. Walker, D.K.; Drummond, M.J.; Dickinson, J.M.; Borack, M.S.; Jennings, K.; Volpi, E.; Rasmussen, B.B. Insulin increases mRNA abundance of the amino acid transporter SLC7A5/LAT1 via an mTORC1-dependent mechanism in skeletal muscle cells. Physiol. Rep. 2014, 2, e00238. [Google Scholar] [CrossRef]
  96. Taylor, P.M. Role of amino acid transporters in amino acid sensing. Am. J. Clin. Nutr. 2014, 99, 223s–230s. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Takimoto, E.; Fukushima, N.; Kanetani, T.; Nishimura, S. The Optimum Method for Additional Daily Nutrient Requirements for Pregnant Japanese Black Cattle in Later Pregnancy. Bull. Okayama Prefect. Cent. Anim. Husb. Res. 2019, 9, 11–16. (In Japanese) [Google Scholar]
  98. Institute of Laboratory Animal Science Research Support Center Kagoshima University. The Guide for the Care and Use of Experimental Animals. (In Japanese). Available online: https://www.kufm.kagoshima-u.ac.jp/~animal/tebiki.htm (accessed on 9 August 2021).
  99. Muroya, S.; Oe, M.; Nakajima, I.; Ojima, K.; Chikuni, K. CE-TOF MS-based metabolomic profiling revealed characteristic metabolic pathways in postmortem porcine fast and slow type muscles. Meat Sci. 2014, 98, 726–735. [Google Scholar] [CrossRef]
Metabolites 11 00582 g001 550
Figure 1. HCA using the top 50 statistically different metabolites between LN (green) and HN (red) fetuses. The row displays the metabolite and the column represents the sample. Metabolites with relatively low contents are displayed in blue, while metabolites with relatively high contents are displayed in brown. The brightness of each color corresponds to the magnitude of the difference when compared with the average value. UDP-GlcNAc: uridine diphosphate N-acetylglucosamine/uridine diphosphate N-acetylgalactosamine, 8-OHdG: 8-hydroxy-2′-deoxyguanosine, PEP: phosphoenolpyruvic acid, Frc 1,6-diP: fructose 1,6-diphosphate, DHAP: dihydroxyacetone phosphate.
Figure 1. HCA using the top 50 statistically different metabolites between LN (green) and HN (red) fetuses. The row displays the metabolite and the column represents the sample. Metabolites with relatively low contents are displayed in blue, while metabolites with relatively high contents are displayed in brown. The brightness of each color corresponds to the magnitude of the difference when compared with the average value. UDP-GlcNAc: uridine diphosphate N-acetylglucosamine/uridine diphosphate N-acetylgalactosamine, 8-OHdG: 8-hydroxy-2′-deoxyguanosine, PEP: phosphoenolpyruvic acid, Frc 1,6-diP: fructose 1,6-diphosphate, DHAP: dihydroxyacetone phosphate.
Metabolites 11 00582 g001
Metabolites 11 00582 g002 550
Figure 2. qPCR results of highly down- and upregulated fetal skeletal muscle genes in microarray analysis. The ratios of the normalized gene expression of the LN (green) to HN (red) fetuses are shown as columns. RPL7 was used as the internal control. Error bars indicate SEM. * and + indicate differences between the LN and HN fetuses at p < 0.05 and p < 0.10, respectively. PDK4: pyruvate dehydrogenase kinase 4, NOS2: nitric oxide synthase 2, CPT1B: carnitine palmitoyltransferase 1B, MTHFD2: methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 2, methenyltetrahydrofolate cyclohydrolase, RICTOR: RPTOR independent companion of MTOR complex 2, FOXP1: forkhead box P1, PPARA: peroxisome proliferator activated receptor alpha, RETN: resistin, SLC7A5: solute carrier family 7 member 5, GPNMB: glycoprotein nmb, ALDH1L2: aldehyde dehydrogenase 1 family member L2.
Figure 2. qPCR results of highly down- and upregulated fetal skeletal muscle genes in microarray analysis. The ratios of the normalized gene expression of the LN (green) to HN (red) fetuses are shown as columns. RPL7 was used as the internal control. Error bars indicate SEM. * and + indicate differences between the LN and HN fetuses at p < 0.05 and p < 0.10, respectively. PDK4: pyruvate dehydrogenase kinase 4, NOS2: nitric oxide synthase 2, CPT1B: carnitine palmitoyltransferase 1B, MTHFD2: methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 2, methenyltetrahydrofolate cyclohydrolase, RICTOR: RPTOR independent companion of MTOR complex 2, FOXP1: forkhead box P1, PPARA: peroxisome proliferator activated receptor alpha, RETN: resistin, SLC7A5: solute carrier family 7 member 5, GPNMB: glycoprotein nmb, ALDH1L2: aldehyde dehydrogenase 1 family member L2.
Metabolites 11 00582 g002
Metabolites 11 00582 g003 550
Figure 3. Expression of fetal skeletal muscle genes associated with significant metabolisms/pathways extracted in GO and KEGG analyses. The ratios of the normalized gene expression of the LN (green) to HN (red) fetuses are shown as columns. RPL7 was used as the internal control. Error bars indicate SEM. * and + indicate differences between the LN and HN fetuses at p < 0.05 and p < 0.10, respectively.
Figure 3. Expression of fetal skeletal muscle genes associated with significant metabolisms/pathways extracted in GO and KEGG analyses. The ratios of the normalized gene expression of the LN (green) to HN (red) fetuses are shown as columns. RPL7 was used as the internal control. Error bars indicate SEM. * and + indicate differences between the LN and HN fetuses at p < 0.05 and p < 0.10, respectively.
Metabolites 11 00582 g003
Metabolites 11 00582 g004 550
Figure 4. Hypothetic scheme of influences of maternal nutrient restriction on the fetal skeletal muscle molecular pathways and metabolisms. Genes and metabolites are shown in blue and green, respectively. Increase and decrease are indicated in small black and red arrows, respectively. Underlined metabolisms/pathways are significantly extracted in GO/KEGG pathway analysis. Grayed events and metabolites were not determined in this study, but some of these were confirmed in previous studies.
Figure 4. Hypothetic scheme of influences of maternal nutrient restriction on the fetal skeletal muscle molecular pathways and metabolisms. Genes and metabolites are shown in blue and green, respectively. Increase and decrease are indicated in small black and red arrows, respectively. Underlined metabolisms/pathways are significantly extracted in GO/KEGG pathway analysis. Grayed events and metabolites were not determined in this study, but some of these were confirmed in previous studies.
Metabolites 11 00582 g004
Table
Table 1. Phenotypic effect of maternal nutrient restriction on fetal carcass.
Table 1. Phenotypic effect of maternal nutrient restriction on fetal carcass.
HN (n = 5) LN (n = 6)
MeanSEMean *SE *p-Value
Age (d)261.61.5260.71.60.684
BW (kg)32.50.523.42.40.005
Total muscle (g)4889.4128.73496.8401.30.009
% of total BW15.00.414.90.40.733
LT muscle (g)272.07.4187.024.10.008
% of total muscle5.60.15.30.20.183
* The values of the LN group were reported previously [10].
Table
Table 2. Top 20 differently expressed metabolites in fetal LT muscle between LN and HN groups *.
Table 2. Top 20 differently expressed metabolites in fetal LT muscle between LN and HN groups *.
CompoundLNHNRatio (LN/HN)p-Value
MeanSEMeanSE
Carnosine0.01970.00110.01480.00071.30.006
Glutamine0.02150.00090.01680.00081.30.008
Glycerol0.13670.01050.10380.00721.30.022
Creatine0.05500.00260.04470.00221.20.023
N6-Methyllysine0.00370.00030.00260.00021.50.028
Phosphorylcholine0.01250.00110.00890.00071.40.034
Phenylalanine0.00970.00090.00660.00071.50.039
myo-Inositol 2-phosphate0.00070.00000.00090.00010.80.039
2-Aminoethylphosphonic acid0.00020.00000.00030.00000.60.041
Proline0.03860.00330.02920.00161.30.043
Alanine0.15650.00440.12760.01101.20.059
Putrescine0.00220.00020.00150.00021.40.065
Creatinine0.01890.00110.01510.00161.20.070
Gamma-amino butyric acid0.00710.00050.00550.00051.30.077
Histidine0.00340.00060.00190.00031.80.081
Glycine0.07010.00650.05060.00731.40.090
2-Hydroxyvaleric acid0.00020.00000.00020.00010.60.098
Taurine0.02490.00190.03070.00260.80.111
S-Adenosylmethionine0.00060.00010.00080.00010.70.113
Dihydroxyacetone phosphate0.00820.00160.00470.00111.70.118
* Values in table are relative content levels (arbitrary unit). LN and HN: Low and high nutrition treatment, respectively.
Table
Table 3. Top 20 fetal muscle metabolisms different between LN and HN fetuses *.
Table 3. Top 20 fetal muscle metabolisms different between LN and HN fetuses *.
Metabolism/PathwayHits/Total Metabolitesp-ValueIncreased in LNDecreased in LN
Pyrimidine metabolism2/390.003Glutamine, Deoxyuridine
Aminoacyl-tRNA
biosynthesis		13/480.007Glutamine, Phenylalanine, Proline, Alanine, Valine, Histidine, Glycine,
Tryptophan, Serine,
Tyrosine, Glutamate,
Arginine		Methionine
Glycerolipid metabolism2/160.007Glycerol, DHAP
Arginine biosynthesis4/140.007Glutamine, Ornithine,
Glutamate, Arginine
Alanine, aspartate and
glutamate metabolism		5/280.010Glutamine, Alanine, GABA, Pyruvate,
Glutamate
Glutathione metabolism6/280.011Putrescine, Glycine,
Ornithine, Spermidine, Spermine, Glutamate
Histidine metabolism3/160.012Carnosine, Histidine,
Glutamate
Arginine and proline
metabolism		11/380.012Creatine, Proline,
Putrescine, GABA,
Ornithine, Spermidine,
Pyruvate, Spermine,
Glutamate, Arginine		SAM
Glyoxylate and dicarboxylate metabolism5/320.013Glutamine, Glycine,
Serine, Pyruvate,
Glutamate
Phosphonate and
phosphinate metabolism		2/60.017Phosphorylcholine2-AEP
Glycerophospholipid
metabolism		2/360.018Phosphorylcholine, DHAP
Galactose metabolism1/270.022Glycerol
Primary bile acid
biosynthesis		2/460.025GlycineTaurine
Purine metabolism2/650.025Glutamine, Allantoin
D-Glutamine and
D-glutamate metabolism		4/60.029Glutamate, Glutamine
Nitrogen metabolism2/60.029Glutamate, Glutamine
Glycine, serine and threonine metabolism5/330.036Creatine, Glycine, Serine, PyruvateCystathionine
β-Alanine metabolism4/210.040Carnosine, Histidine,
Spermidine, Spermine
Glycolysis/Gluconeogenesis3/260.058Pyruvate, PEP, DHAP
Selenocompound
metabolism		1/200.059Alanine
* Metabolisms and pathways that were extracted in MSEA using top 50 different metabolites are listed above.
Table
Table 4. Fetal growth and energy associated metabolisms and pathways extracted from downregulated genes *.
Table 4. Fetal growth and energy associated metabolisms and pathways extracted from downregulated genes *.
CategoryTermp-ValueFold EnrichmentValidated Genes
KEGG Pathway
bta04940:Type I diabetes mellitus<0.0013.8050
bta00010:Glycolysis/Gluconeogenesis<0.0013.1407GPI, ENO3, PFKM, PGAM1, FBP2
bta04010:MAPK signaling pathway<0.0011.9194MAPK3, NFKB1, RELA
bta04066:HIF-1 signaling pathway<0.0012.4950ANGPT4, NOS2, NOS3, ENO3, MAPK3, NFKB1, PLCG1, RELA
bta01200:Carbon metabolism<0.0012.1974GPI, PGAM1, ENO3, PFKM, FBP2
bta04151:PI3K-Akt signaling pathway<0.0011.5906ANGPT4, NOS3, MAPK3, NFKB1, PRKAA1, RELA, THBS1
bta04919:Thyroid hormone signaling pathway<0.0012.1196MAPK3, PLCG1
bta00030:Pentose phosphate pathway0.0033.4713GPI, PFKM, FBP2
bta04910:Insulin signaling pathway0.0051.8377FBP2, MAPK3, PRKAA1
bta04115:p53 signaling pathway0.0072.2001THBS1
GO: Biological Process
GO:0006096~glycolytic process<0.0014.3200GPI, ENO3, PGAM1
GO:0034097~response to cytokine<0.0014.0909NFKB1, RELA
GO:0006094~gluconeogenesis<0.0014.2078GPI, FBP2, PGAM1
GO:0000122~negative regulation of transcription from RNA polymerase II promoter0.0021.4727FOXP1, NFKB1, RELA
GO:0001525~angiogenesis0.0021.9050GPI, NOS3, ANGPT4, ANGPTL4
GO:0042981~regulation of apoptotic process0.0021.9474PINK1, RELA
GO:0050995~negative regulation of lipid catabolic process0.0095.4545PRKAA1
GO:0090200~positive regulation of release of cytochrome c from mitochondria0.0093.6172PINK1
GO:0043536~positive regulation of blood vessel endothelial cell migration0.0104.2078ANGPTL4, PLCG1, THBS1
GO:0042177~negative regulation of protein catabolic process0.0103.1418NOS2, SIRT2, RELA
* Top 20 metabolisms/pathways extracted at p < 0.05 are listed. GPI: glucose-6-phosphate isomerase, ENO3: enolase 3, PFKM: phosphofructokinase, muscle, FBP2: fructose-bisphosphatase 2, PGAM1: phosphoglycerate mutase 1, ANGPT4: angiopoietin 4, NOS2: nitric oxide synthase 2, NOS3: nitric oxide synthase 3, MAPK3: mitogen-activated protein kinase 3, NFKB1: nuclear factor kappa B subunit 1, PLCG1: phospholipase C gamma 1, PRKAA1: protein kinase AMP-activated catalytic subunit alpha 1, RELA: RELA proto-oncogene, NF-kB subunit, THBS1: thrombospondin 1, FOXP1: forkhead box P1, SIRT2: sirtuin 2.
Table
Table 5. Fetal growth and energy associated metabolisms and pathways extracted from upregulated genes *.
Table 5. Fetal growth and energy associated metabolisms and pathways extracted from upregulated genes *.
CategoryTermp-ValueFold EnrichmentValidated Genes
KEGG Pathway
bta01230:Biosynthesis of amino acids<0.0012.8335
bta00970:Aminoacyl-tRNA
biosynthesis		<0.0013.4187
bta00260:Glycine, serine and
threonine metabolism		<0.0013.4636
bta00072:Synthesis and degradation of ketone bodies0.0016.4548BDH1
bta04512:ECM-receptor interaction0.0022.3124
bta00670:One carbon pool by folate0.0034.6021MTHFD2, ALDH1L2
bta04550:Signaling pathways
regulating pluripotency of stem cells		0.0031.9581
bta04974:Protein digestion and
absorption		0.0042.2541
GO: Biological Process
GO:0006730~one-carbon metabolic process<0.0014.4172MTHFD2, ALDH1L2
GO:0046653~tetrahydrofolate
metabolic process		<0.00110.6749MTHFD2
* Top 20 metabolisms/pathways extracted at p < 0.05 are listed. BDH1: 3-hydroxybutyrate dehydrogenase 1, MTHFD2: methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 2, methenyltetrahydrofolate cyclohydrolase, ALDH1L2: aldehyde dehydrogenase 1 family member L2.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Muroya, S.; Zhang, Y.; Kinoshita, A.; Otomaru, K.; Oshima, K.; Gotoh, Y.; Oshima, I.; Sano, M.; Roh, S.; Oe, M.; et al. Maternal Undernutrition during Pregnancy Alters Amino Acid Metabolism and Gene Expression Associated with Energy Metabolism and Angiogenesis in Fetal Calf Muscle. Metabolites 2021, 11, 582. https://doi.org/10.3390/metabo11090582

AMA Style

Muroya S, Zhang Y, Kinoshita A, Otomaru K, Oshima K, Gotoh Y, Oshima I, Sano M, Roh S, Oe M, et al. Maternal Undernutrition during Pregnancy Alters Amino Acid Metabolism and Gene Expression Associated with Energy Metabolism and Angiogenesis in Fetal Calf Muscle. Metabolites. 2021; 11(9):582. https://doi.org/10.3390/metabo11090582

Chicago/Turabian Style

Muroya, Susumu, Yi Zhang, Aoi Kinoshita, Kounosuke Otomaru, Kazunaga Oshima, Yuji Gotoh, Ichiro Oshima, Mitsue Sano, Sanggun Roh, Mika Oe, and et al. 2021. "Maternal Undernutrition during Pregnancy Alters Amino Acid Metabolism and Gene Expression Associated with Energy Metabolism and Angiogenesis in Fetal Calf Muscle" Metabolites 11, no. 9: 582. https://doi.org/10.3390/metabo11090582

APA Style

Muroya, S., Zhang, Y., Kinoshita, A., Otomaru, K., Oshima, K., Gotoh, Y., Oshima, I., Sano, M., Roh, S., Oe, M., Ojima, K., & Gotoh, T. (2021). Maternal Undernutrition during Pregnancy Alters Amino Acid Metabolism and Gene Expression Associated with Energy Metabolism and Angiogenesis in Fetal Calf Muscle. Metabolites, 11(9), 582. https://doi.org/10.3390/metabo11090582

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