Next Article in Journal
Molecular and Biological Investigation of Isolated Marine Fungal Metabolites as Anticancer Agents: A Multi-Target Approach
Next Article in Special Issue
Impacts of Different Prenatal Supplementation Strategies on the Plasma Metabolome of Bulls in the Rearing and Finishing Phase
Previous Article in Journal
Assessing Metabolic Markers in Glioblastoma Using Machine Learning: A Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Cow-Calf Supplementation on Gene Expression, Processes, and Pathways Related to Adipogenesis and Lipogenesis in Longissimus thoracis Muscle of F1 Angus × Nellore Cattle at Weaning

by
Germán Darío Ramírez-Zamudio
1,
Maria Júlia Generoso Ganga
2,
Guilherme Luis Pereira
2,3,
Ricardo Perecin Nociti
1,
Marcos Roberto Chiaratti
4,
Reinaldo Fernandes Cooke
5,
Luis Artur Loyola Chardulo
2,3,
Welder Angelo Baldassini
2,3,
Otávio Rodrigues Machado-Neto
2,3 and
Rogério Abdallah Curi
2,3,*
1
College of Animal Science and Food Engineering, São Paulo University (USP), Pirassununga 13635-900, SP, Brazil
2
School of Agriculture and Veterinary Sciences (FCAV), São Paulo State University (UNESP), Jaboticabal 14884-900, SP, Brazil
3
School of Veterinary Medicine and Animal Science (FMVZ), São Paulo State University (UNESP), Botucatu 18618-681, SP, Brazil
4
Department of Genetics and Evolution, Federal University of São Carlos (UFSCAR), São Carlos 13565-905, SP, Brazil
5
Department of Animal Science, Texas A&M University, College Station, TX 77843, USA
*
Author to whom correspondence should be addressed.
Metabolites 2023, 13(2), 160; https://doi.org/10.3390/metabo13020160
Submission received: 19 December 2022 / Revised: 12 January 2023 / Accepted: 18 January 2023 / Published: 21 January 2023
(This article belongs to the Special Issue Feedlot Ruminant Nutrition and Metabolism)

Abstract

:
The aim of this study was to identify differentially expressed genes, biological processes, and metabolic pathways related to adipogenesis and lipogenesis in calves receiving different diets during the cow-calf phase. Forty-eight uncastrated F1 Angus × Nellore males were randomly assigned to two treatments from thirty days of age to weaning: no creep feeding (G1) or creep feeding (G2). The creep feed offered contained ground corn (44.8%), soybean meal (40.4%), and mineral core (14.8%), with 22% crude protein and 65% total digestible nutrients in dry matter. After weaning, the animals were feedlot finished for 180 days and fed a single diet containing 12.6% forage and 87.4% corn-based concentrate. Longissimus thoracis muscle samples were collected by biopsy at weaning for transcriptome analysis and at slaughter for the measurement of intramuscular fat content (IMF) and marbling score (MS). Animals of G2 had 17.2% and 14.0% higher IMF and MS, respectively (p < 0.05). We identified 947 differentially expressed genes (log2 fold change 0.5, FDR 5%); of these, 504 were upregulated and 443 were downregulated in G2. Part of the genes upregulated in G2 were related to PPAR signaling (PPARA, SLC27A1, FABP3, and DBI), unsaturated fatty acid synthesis (FADS1, FADS2, SCD, and SCD5), and fatty acid metabolism (FASN, FADS1, FADS2, SCD, and SCD5). Regarding biological processes, the genes upregulated in G2 were related to cholesterol biosynthesis (EBP, CYP51A1, DHCR24, and LSS), unsaturated fatty acid biosynthesis (FADS2, SCD, SCD5, and FADS1), and insulin sensitivity (INSIG1 and LPIN2). Cow-calf supplementation G2 positively affected energy metabolism and lipid biosynthesis, and thus favored the deposition of marbling fat during the postweaning period, which was shown here in an unprecedented way, by analyzing the transcriptome, genes, pathways, and enriched processes due to the use of creep feeding.

Graphical Abstract

1. Introduction

Different strategies have been used in beef cattle to increase intramuscular fat (IMF) deposition. Early weaning and supplementation during the cow-calf phase are nutritional management practices that can induce different metabolic adaptations when compared to conventionally weaned animals [1]. The approach known as creep feeding consists of supplementation during the cow-calf (lactation) phase with grains or forage to obtain heavier individuals at weaning, to reduce the time to carcass finishing for slaughter, and to allow the dam to rest [2,3].
It is well-known that animals can respond to different environmental/nutritional factors by exhibiting phenotypic plasticity as a result of changes in gene expression patterns [4]. Thus, animal characteristics can be modified by nutritional modulations, which means that dietary exposures can have consequences for growth and health [5]. New evidence is constantly emerging that alterations in gene expression and in the phenotype of individuals, as a result of nutritional stimuli, may be due to epigenetic factors, including DNA methylation [6,7,8]. These modifications in gene expression can occur at specific loci or on a genomic scale [9,10,11,12,13].
In Bos taurus, it was identified that glucose is the best precursor for the synthesis of fatty acids in the intramuscular adipose tissue and the only carbonic compound with a higher concentration of its carbons in the intramuscular adipose tissue compared to the subcutaneous tissue, and in addition, lower glucose incorporation in post-pubertal animals was also noted [14]. Thus, the adequate intake of concentrate during the lactation period, especially the one rich in starch (grains), favors the proliferation and growth of adipocytes. This is due to the increase in glucose metabolism, which can lead to increased marbling without, however, increasing fat in the viscera, which contributes to overall meat quality and higher carcass yield [15]. Furthermore, Lopez et al. [16] proposed that higher levels of crude protein in the diet increase the digestion and intestinal absorption of starch, and promote the elevation of plasma insulin and glucose, which, consequently, increases the deposition of IMF, by favoring precursors of fatty acids.
Brazil, the world’s largest exporter of unprocessed beef for almost 20 years [17], has been witnessing the growing use of crossbred animals (Bos taurus × Bos indicus), especially Angus × Nellore. This strategy is used for the production of more tender meat with more IMF in an attempt to add value to the final product and to meet the demands of consumers who are willing to pay more for quality. IMF deposition in beef can positively influence sensory attributes such as flavor, juiciness, and tenderness [18]. In contrast, less IMF (marbling), which is observed in animals with a predominance of the Bos indicus genotype, compromises the sensory attributes of meat. This fact can be explained by the stimulation of fibrogenesis with declining intramuscular adipogenesis, with a consequent increase in connective tissue content [19]. However, the molecular mechanisms that control the growth of these tissues have not yet been completely elucidated.
Considering the lack of information about the changes in gene expression that occur in crossbred Bos taurus × Bos indicus calves supplemented during the cow-calf phase, this study aimed to evaluate the effect of creep feeding on gene expression, biological processes, and metabolic pathways related to adipogenesis and lipogenesis in Longissimus thoracis (LT) muscle biopsies collected at weaning from F1 Angus × Nellore cattle. Such data are of great economic importance for cattle production systems that aim to manipulate marbling and to obtain better quality meat to attend more demanding markets that pay for quality.

2. Materials and Methods

2.1. Animals

Forty-eight uncastrated (intact) F1 Angus × Nellore males born to the same Aberdeen Angus (Bos taurus) sire (half-siblings) were used.
The experiment (cow-calf phase) was conducted on a commercial farm located in the city of Anhembi, State of São Paulo, Brazil. The animals were submitted to two different treatments during most of the cow-calf phase (from 30 days of age to weaning—approximately 210 days), with 24 animals per treatment: group 1 (G1/T1)—no creep feeding (conventional weaning), and group 2 (G2/T2)—creep feeding. In the creep feeding system, the animals were kept with their mothers in the same paddock but had exclusive and free access to supplement corresponding to approximately 1% of body weight during all the cow-calf phase. The creep feed offered to animals of G2 (from 30 to 210 days of life) consisted of dry matter (22% crude protein and 65% total digestible nutrients) containing ground corn (44.8%), soybean meal (40.4%), and mineral core (14.8%).
After weaning (mean of 210 days), the animals of the 2 treatments were transferred to an experimental feedlot (Botucatu, São Paulo, Brazil), where they were housed in covered collective pens (three animals/pen with 10 m2 per animal) for approximately 180 days. The two groups received the same diet, containing 12.6% forage and 87.4% corn-based concentrate. The diet was formulated with the RLM 3.3 software (Ração de Lucro Máximo–Maximum Profit Ration, Piracicaba, São Paulo, Brazil) [20] using the NRC Tropicalizado ESALQ system and consisted of corn, soybean meal, Tifton hay, sugarcane bagasse, urea, and vitamin–mineral supplement (Table S1). The diet was offered ad libitum, twice a day at 8 am and 4 pm.
The animals were weighed after a 16 h fast at the beginning of the cow-calf phase (initial weight—BWi), at the end of weaning (weaning weight—WW), and at the end of the feedlot period (final weight—BWf). The average daily weight gain 1 (ADG1: beginning of cow-calf phase to weaning) was calculated from BWi and WW, and ADG2 (weaning to the end of the feedlot period) from BWf and WW.

2.2. Collection of Muscle Tissue at Weaning and Slaughter

At weaning (210 days), fragments of the LT muscle were collected by biopsy from 12 animals randomly selected from each treatment (n = 24). For biopsy, the lumbar region was shaved, and a local anesthetic was subcutaneously administered. The biopsies were performed at the height of the 13th rib. After cleaning the biopsy site, a 1 cm incision was made with a scalpel and a sterile Bergstrom biopsy needle was used to obtain 1 g of muscle tissue. The sample was immediately transferred to liquid nitrogen and stored in an ultra-freezer at −80 °C.
After the feedlot period, the 48 animals were slaughtered using the cerebral concussion technique and sectioning of the jugular vein. The carcasses were identified, washed, and divided into two halves. The half-carcasses were weighed individually to obtain the hot carcass weight (HCW) and kept in a cold room for approximately 24 h at 1 °C.
After cooling, the carcasses were removed from the cold room and weighed. After weighing, the LT muscle of the left half-carcass was separated and backfat thickness (BFT) and rib eye area (REA) were evaluated between the 12th and 13th ribs before deboning. Beef samples (sirloin steaks) were collected between the 12th and 13th ribs (cranial direction) of the left half-carcass and used for the laboratory analysis of physicochemical quality attributes.

2.3. Analysis of Meat Quality

The following physicochemical meat quality attributes were analyzed in the 48 slaughtered animals (n = 24/treatment): marbling score (MS), total lipids/intramuscular fat percentage (IMF), and Warner–Blatzler shear force (WBSF). The MS was determined by a single, previously trained evaluator following the Brazil Beef Quality reference standards (https://www.brazilbeefquality.com/) (accessed on 12 June 2021). These standards are numbered from 1 to 11 and provide a point scale ranging from 100 to 1100 (adapted from AUS-MEAT [21]), where the closer to 1100, the more marbled the meat, and the closer to 100, the less marbled the meat. The IMF was determined by infrared spectroscopy using a FoodScanTM (Foss NIRSystems, Laurel, MD, USA). The procedure standardized by Shackelford et al. [22] was adopted to measure WBSF at 7 and 14 days of aging (WBSF7 and WBSF14).

2.4. Statistical Analysis of Weight, Weight Gain, Carcass, and Meat Data

The weight and weight gain data and carcass and meat traits (BWi, WW, ADG1, BWf, ADG2, HCW, REA, BFT, MS, IMF, WBSF7, and WBSF14) of the 48 samples (n = 24/group) were analyzed regarding the presence of outliers, homogeneity of variance, and normality of residuals. The data were expressed as means and their respective standard errors. Transformations of the variables to provide an approximation of the normal distribution were not necessary. The means of the groups were compared by the t-test to detect significant differences (p < 0.05). These analyses were performed using the R software v.4.2.1 (Vienna, Austria) [23].

2.5. Analysis of Differential Gene Expression

2.5.1. RNA Extraction and Sequencing

After the extraction of total RNA, 24 genomic libraries were prepared, which consisted of 12 samples from each group/treatment (G1 and G2) collected at weaning. Total RNA was extracted individually from 100 mg of LT muscle using TRIzol® (Life Technologies, Carlsbad, CA, USA), according to the manufacturer’s instructions, and its quality was analyzed in a Bioanalyzer 2100® (Agilent, Santa Clara, CA, USA). A minimum RNA integrity number (RIN) ≥ 7 was adopted to ensure adequate quality of the total RNA. RNAs with a poly-A tail, mainly the mRNA, were purified from total RNA using oligo-dT beads. It is well-known that many eukaryotic non-coding RNAs can have a poly-A tail in intermediary or in mature forms [24] and, in this way, eventually be captured by oligo-dT bead columns.
The cDNA libraries of each sample were prepared and multiplexed using the TruSeq RNA Sample Preparation kit (Illumina, San Diego, CA, USA) from 2 µg of total RNA, according to the TruSeq RNA Sample Preparation kit v2 guide (Illumina). The Bioanalyzer 2100® (Agilent) was used to estimate the average size of the libraries and quantitative PCR (RT-qPCR) using the KAPA Library Quantification kit (KAPA Biosystems, Wilmington, MA, USA) to quantify them. Clustering and sequencing were performed in one lane using the HiSeq2500 kit v4 2 × 100 bp (Illumina) to produce 100 bp paired-end (PE) reads. The HiSeq 2500® sequencer (Illumina) was used to reach a minimum coverage of 16 million reads per sample.

2.5.2. Mapping of Sequences to the Reference Genome and Identification of Differentially Expressed Genes

The sequence data generated by the HiSeq System Illumina platform were converted to FastQ format and separated by library (multiplexed data) using the Casava 1.8.2 software (Illumina). The FastQC v.0.11.9 software [25] was used to analyze the quality of raw reads. Adapter sequences and low-quality sequences were removed using Fastp v.0.23.1 [26]. After this step, the quality of the reads was reassessed by combined visualization of all FastQC outputs using the MultiQC v.1.13 program [27] to confirm improvements in quality. Next, the reads were mapped to the bovine reference genome (Bos taurus—ARS-UCD1.3), available at: http://www.ensembl.org/Bos_taurus/Info/Index/ (accessed on 19 August 2021), using the STAR v.2.7.20 program [28]. Mapping was performed independently for each sample. For each library, a file with the .bam extension was generated, which contained the alignment of the fragments to the reference genome. Mapped reads were counted using featurecounts v.2.0.3 [29] and only PE reads mapped to a single position of the genome (uniquely mapped PE reads) and to known chromosomes were used for differential gene expression analysis.
Differential gene expression was compared between the different groups/treatments at weaning (i.e., G1 vs. G2). First, graphical principal component analysis (PCA) of the read counts normalized to counts per million (CPM) was performed using the factoextra package [30] of the R software [23] to divide the samples based on gene expression patterns, examining the level of similarity/dissimilarity between groups. Differentially expressed genes (DEGs) were identified using the method implemented in edgeR v.3.40.0 [31] of the R software [23] for parameter estimation by the maximum likelihood method. Generalized linear models were used, assuming a negative binomial distribution of the count data. For this purpose, size factors were computed by the trimmed mean of M-values (TMM) for each pair of samples and overdispersion parameters were estimated for each gene by the Cox–Reid method [32]. The expression of each gene was calculated as the mean expression in all samples of each group and is reported as the mean of the logarithmic function of CPM. The fold change was calculated as the logarithmic function of the ratio between the CPM of G2 and G1 for each gene. The p-value associated with the difference in gene expression between groups was obtained by the likelihood ratio test. The Benjamini–Hochberg procedure [33] was used to control the false discovery rate (FDR). A log2 fold change of 0.5 and significance adjusted to FDR < 0.05 were adopted to identify DEGs. The expression profiles of 13 constitutive genes were analyzed to evaluate the quality of sequencing. These genes were: actin beta—ACTB; beta-2-microglobulin—B2M; glyceraldehyde-3-phosphate dehydrogenase—GAPDH; glucuronidase beta—GUSB; hydroxymethylbilane synthase—HMBS; hypoxanthine phosphoribosyltransferase 1—HPRT1; phosphoglycerate kinase 1—PGK1; peptidylprolyl isomerase A—PPIA; ribosomal protein L13a—RPL13A; ribosomal protein lateral stalk subunit P0—RPLP0; succinate dehydrogenase complex flavoprotein subunit A—SDHA; TATA-box binding protein—TBP; transferrin receptor—TFRC.

2.5.3. Functional Analysis of Differentially Expressed Genes

To understand the functional role of the genes identified as DEGs between groups at weaning, lists of up- and down-regulated genes were used in over-representation analysis (ORA) of gene ontology terms (GO terms: biological processes) and metabolic pathways (Kyoto Encyclopedia of Genes and Genomes (KEGG) database) using the online Database for Annotation, Visualization, and Integrated Discovery (DAVID v.6.8; https://david.ncifcrf.gov/home.jsp) (accessed on 12 October 2021) [34]. Biological processes and metabolic pathways were defined as enriched in the presence of at least three genes in each pathway or process and a p-value < 0.05.
The functional annotation clustering function of DAVID and the ClueGO package [35] of Cytoscape v.3.9.1 were used to identify the relationships between enriched processes and pathways. Due to their relationships with each other and their potential relationships with intramuscular fat content in cattle, processes and pathways were highlighted. To identify/visualize the participation of DEGs in the different enriched and highlighted biological processes and pathways, a HeatMap was generated using the ComplexHeatMap package [36] of the R software. The ClueGO and CluePedia packages [37] of Cytoscape were used to generate networks of shared DEGs between KEGG pathways and biological processes (FDR < 5%).
The online STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) database (version 11.5), which collects and integrates information on functional interactions between genes/proteins for a large number of organisms [38], was used to reveal and visualize functional interactions between the DEGs of biological processes and KEGG pathways identified as enriched by DAVID and highlighted due to their relationships visualized in ClueGO and potential relationships with intramuscular adipogenesis and lipogenesis. Protein–protein interactions (PPIs) with a confidence score > 0.4 (a commonly used threshold) and an FDR < 0.05 were considered and are shown in the graph. The Markov cluster algorithm (MCL), an unsupervised clustering algorithm for graphs based on simulation of stochastic flow, was used for the clustering of genes (nodes) using the default inflation parameter 3 of STRING.
The Bos taurus database was used for the analyses that employed the DAVID, ClueGO/CluePedia, and STRING programs/packages.

3. Results

3.1. Pre- and Post-Weaning Performance, Carcass and Meat Quality

Table 1 shows the body weight of the animals at the different time points, weight gain, and carcass-related information. Significant differences (p < 0.05) between groups/treatments were observed for WW and ADG1. Regarding WW, animals of G1 (no supplementation) were lighter and those of G2 (creep feeding) were heavier. The same situation was observed for ADG1, i.e., G2 > G1. Animals of G1 and G2 did not differ in terms of BWf, and likewise for ADG2 during the feedlot period. No significant differences in BWi or HCW were observed between groups.
Table 2 shows the quality variables, which were measured or estimated objectively or subjectively in the carcasses or meat samples of animals by laboratory analysis. No significant differences (p > 0.05) between groups were observed for REA, WBSF7, or WBSF14. On the other hand, there were significant differences (p < 0.05) in BFT, IMF, and MS. Animals of G2 (creep feeding) exhibited a higher mean BFT than G1 animals (no supplementation). Regarding IMF, an objective measure for the assessment of marbling, G2 animals had a higher fat content/percentage than G1 animals. Visual differences in relation to the fat content between LT muscle samples collected from animals of G1 and G2 are shown in Figure S1. The MI, a subjective measure of marbling, was higher in G2 animals compared to G1.

3.2. Analysis of Differential Gene Expression

3.2.1. Concentration and Integrity of Total RNA

The mean total RNA concentration extracted from the 24 samples was 260.27 ng/µL. The 260/280 nm (nucleic acid/protein) and 260/230 nm (nucleic acid/extraction contaminants) ratios were approximately 1.9, a value considered to be adequate. The mean contamination with genomic DNA was 1.04% (range: 0.73% to 1.11%). The mean RIN was 7.6 (range: 7.0 to 8.0). Thus, the samples were intact (all RINs > 7) and free of contaminants.

3.2.2. RNA Sequencing and Mapping of Reads to the Reference Genome

A total of 241.2 million PE reads (2 × 100 bp) were obtained; of these, 230.5 million were uniquely mapped PE reads. The coverage achieved by sequencing was 37X (coverage for all transcripts of all samples). An average of 9.6 million uniquely mapped PE reads were obtained per sample, corresponding to 95.56% of all PE reads generated.
The reads were mapped to 27,607 genes (protein coding and non-coding). However, considering a count of uniquely mapped PE reads ≥ 3 in at least 12 samples, the total number of expressed genes, and thus used in the differential gene expression analyzed, was 16,604. The number of genes detected by functional category after application of a filter/threshold that excludes genes with a low count is presented in Figure S2. Following protein-coding genes (15,672), the largest number of sequenced genes were those encoding transcription factors (656) and non-coding genes such as long non-coding RNA (lncRNA), short nuclear RNA (snRNA), and micro-RNA (miRNA). The latter are related to the control of gene expression and processing of messenger RNA.
The number and percentage of transcripts/fragments aligned to the bovine reference genome identified in samples collected from non-creep-fed (G1) and creep-fed (G2) animals, which showed significant differences between means for marbling measurements—IMF and MS (G1 < G2, p < 0.05), are shown in Table 3 and Table 4, respectively.
The boxplots of the read counts normalized by the size factor showed that the distribution of quartiles was consistent between the samples of the two groups, indicating good quality of the sequencing data (Figure S3).
In relation to the expression profile of the constitutive genes, the expression was similar between the experimental groups (Figure 1).
PCA showed that the first two principal components explained more than 20% of the variation among samples (Figure 2). In addition, the formation of clearly distinct groups of samples was observed at weaning, indicating an evident difference in the expression of genes between treatments. This fact illustrates the effect of cow-calf supplementation (creep feeding) on gene expression at weaning.

3.2.3. Identification of Differentially Expressed Genes

Fragments of LT muscle obtained by biopsy were used to identify differences in global gene expression between G1 and G2 at weaning. A total of 947 DEGs were identified (log2 fold change < −0.5 or >0.5, FDR < 5%) between groups at weaning; of these, 443 were downregulated and 504 were upregulated in G2 (creep feeding).
Figure 3 shows a Volcano plot that illustrates genes differentially expressed in G1 × G2 at weaning. The top 30 DEGs for comparison, with an adjusted p-value (FDR) < 0.05, are shown in Table 5. The complete list of the differentially expressed genes, with log2 fold change, p-value, and adjusted p-value of down- and up-regulated genes for comparison between groups at weaning, is provided in Table S2.

3.2.4. Functional Enrichment Analysis of Differentially Expressed Genes

Regarding the genes upregulated in G2 at weaning (n = 504), ORA identified 22 enriched metabolic KEGG pathways (p < 0.05) (Table S3). Among the pathways identified, five are highlighted due to their relationships with each other (identified by DAVID and ClueGO) and their potential relationships with adipose cell proliferation/adipogenesis and synthesis and degradation of intramuscular fat (lipolysis/lipogenesis) in cattle: PPAR signaling pathway (bta03320), steroid biosynthesis (bta00100), biosynthesis of unsaturated fatty acids (bta01040), apelin signaling pathway (bta04371), and fatty acid metabolism (bta01212). The upregulated DEGs in the enriched pathways and additional information are provided in Table 6. Regarding GO terms, 34 enriched biological processes were identified (p < 0.05) (Table S4). Four of these processes are highlighted due to their relationships with each other and potential relationships with intramuscular fat content in cattle: cholesterol biosynthetic process (GO:0006695), unsaturated fatty acid biosynthetic process (GO:0006636), sterol biosynthetic process (GO:0016126), and cellular response to insulin stimulus (GO:0032869). Table 7 shows the upregulated DEGs related to these processes.
Considering the genes downregulated in G2 at weaning (n = 443), 46 enriched pathways were identified (p < 0.05) (Table S5), and 3 of these pathways are highlighted due to the aforementioned relationships: AMPK signaling pathway (bta04152), glucagon signaling pathway (bta04922), and PPAR signaling pathway (bta03320). The downregulated DEGs participating in these pathways and additional information are provided in Table 6. Regarding GO terms, 52 enriched biological processes were identified (p < 0.05) (Table S6), and 3 are highlighted due to the aforementioned relationships: positive regulation of lipid storage (GO:0010884), fatty acid metabolic process (GO:0006631), and negative regulation of glycolytic process (GO:0045820). Table 7 shows the downregulated DEGs related to these processes.
Figure 4 shows the HeatMap that illustrates the relationship between each of the 52 up- and down-regulated DEGs identified in G2 at weaning and the 14 biological processes and metabolic pathways enriched in ORA, which are highlighted due to their relationships with each other and potential relationships with intramuscular fat content in cattle. This approach permitted to observe the participation of DEGs in one or more processes or pathways, as well as the magnitude of differences in gene expression between treatments.
Networks of shared up- and down-regulated DEGs in G2 at weaning between biological processes (GO_BP) and KEGG pathways, which were highlighted in ORA and due to the potential relationship with intramuscular fat/marbling in Bos taurus, are shown in Figure 5 and Figure 6, respectively.
Figure 7 shows the network of PPI of 29 DEGs upregulated in G2 that participate in the biological processes (GO_BP) and metabolic pathways (KEGG) identified as enriched and highlighted in the G1 × G2 comparison at weaning due to their relationships with each other and potential relationships with intramuscular adipogenesis and lipogenesis. Sixty-nine edges (significant interactions) and six clusters consisting of eight to two genes were identified. The red and yellow clusters were the largest, with eight and five genes, respectively. Three clusters (dark green, light blue, and dark blue) contained only two genes. The mean clustering coefficient was 0.512. Six DEGs were disconnected from the network.
Figure 8 shows the network of PPI of 23 DEGs downregulated in G2 that participate in the biological processes (GO_BP) and metabolic pathways (KEGG) identified as enriched and highlighted in the G1 × G2 comparison at weaning due to their relationships with each other and potential relationships with intramuscular adipogenesis and lipogenesis. Twenty-three edges (significant interactions) and four clusters consisting of eight to two genes were identified. The red, yellow, and dark green clusters were the largest, with eight, three, and three genes, respectively. One of the clusters (dark blue) contained only two genes. The mean clustering coefficient was 0.467. Seven DEGs were disconnected from the network.

4. Discussion

4.1. Pre- and Post-Weaning Performance, Carcass and Meat Quality

In tropical breeding systems, milk production of Nellore cows, which form the basis of the Brazilian cattle herd, no longer meets the calf’s requirements for expression of its growth potential after the third month of lactation [39]. In F1 Angus × Nellore cattle, this interval is believed to be even shorter. Thus, obtaining the nutrients necessary to meet the growth requirements during lactation increasingly depends on the forage consumed by the animal. An important fact of the pasture-based systems in Brazil is that, in addition to the decline in the lactation curve of cows, the mass and nutritive value of pastures decrease due to seasonality, while the calf’s nutritional requirements increase with the progression of growth [40]. Within this context, creep feeding is a strategy used to compensate the nutrient deficits of milk and forage, in addition to stimulating greater muscle development in the animal, which results in a shorter time to slaughter and improvement in carcass quality [41]. However, due to issues related to the offer of supplement, high supplies can affect fiber digestion from pasture, with negative consequences for feed efficiency [42,43] and calf performance [44]. This fact was not observed in the present study; on the contrary, calves that received 1% of their body weight as a supplement exhibited an additional daily weight gain of approximately 100 g and an additional total gain of 14.6 kg at weaning when compared to non-supplemented animals. However, in a meta-analysis conducted by Carvalho et al. [38], the additional daily weight gain was about 200 g and the total gain at weaning was 30 kg in male calves supplemented between 3 and 8 months of age when compared to non-supplemented animals. In the present study, supplementation during the lactation period did not exert any long-term effect on the performance of animals during the finishing period (BWf, ADG2, and HCW; Table 1). This finding suggests that the lack of difference in the efficiency of animals during the postweaning period may be related to gain composition [45], i.e., an early increase in the fat deposition rate may reduce the rate of lean tissue deposition [46]; consequently, weight gain in subsequent phases may be moderate. Similarly, calves supplemented during the lactation period showed a higher degree of carcass finishing (BFT) and greater fat deposition in beef, while no differences in weight gain were observed (Table 2).

4.2. Differentially Expressed Genes and Alterations in Biological Processes and Metabolic Pathways

Intramuscular fat is a desirable characteristic in some niche markets because of its positive effect on flavor, juiciness, and greater consumer perception of meat tenderness [47]. Although this quality trait is attributed to intense cell proliferation during the fetal period, there is a period of postnatal life within the first 250 days, known as the “marbling window”, during which a high-grain diet can specifically lead to the recruitment of intramuscular adipocytes (adipocyte hyperplasia) that provide sites for fat deposition during finishing [48].
Studies indicate that the PPAR signaling pathway is important for the regulation of cell differentiation, energy balance, and lipid metabolism [49,50]. During lipogenesis, activation of this pathway upregulates the expression of FABP, FASN, and SCD [51,52]. The last two genes are important for de novo fatty acid synthesis [53] and were upregulated in creep-fed animals in the present study. Furthermore, Ward et al. [54] reported that, in addition to the upregulation of FASN, an increase in the level of marbling in beef is associated with the expression of Δ-6 and Δ-5 desaturases, enzymes that are regulated by PPARA [55]. This confirms that PPARA also acts on lipogenesis in skeletal muscle [56]. The Δ-5 and Δ-6 desaturases encoded by the FADS1 and FADS2 genes, respectively, belong to the fatty acid desaturase family, key proteins in the first desaturation reaction for endogenous formation of polyunsaturated fatty acids from dietary essential fatty acids (linoleic and linolenic) [57]. Linoleic and linolenic acids are unsaturated fatty acids found in oilseeds such as soybeans [58]. Although these fatty acids may undergo high ruminal biohydrogenation, part of them may escape into the duodenum and be absorbed [59]. In the present study, consumption of approximately 375 g of soybean meal per day is estimated in animals supplemented with 1% of body weight. Thus, part of the unsaturated fatty acids consumed from the supplement may have reached the skeletal muscle, with consequent upregulation of FADS1 and FADS2.
Although calf supplementation during the lactation period had positive effects on the PPAR signaling pathway in lipid metabolism, downregulation of this pathway can downregulate the APOA1 and ANGPTL4 genes [60]. The APOA1 gene encodes the structural and functional protein component of high-density lipoprotein (HDL), which promotes the reverse flow of cholesterol from tissues to excretion in the liver [61]. In the present study, downregulation of APOA1 in supplemented calves may be a long-term effect of unsaturated fatty acid intake from soybean meal. Some studies have shown that the consumption of unsaturated fatty acids promotes a decrease in the expression of APOA1 [62,63]. On the other hand, angiopoietin-like 4 encoded by the ANGPTL4 gene mediates the inactivation of lipoprotein lipase involved in lipid metabolism [64], thereby inhibiting the uptake of triglycerides in adipocytes [65]. Downregulation of APOA1 and ANGPTL4 expression in supplemented calves suggests the accumulation of fatty acids, probably controlled by lipoprotein lipase that reduces the efflux of triglycerides from adipose tissue.
Calf supplementation at 1% of body weight during the lactation period positively influenced the lipogenic program by regulating genes related to the biosynthesis (FASN, SCD, SCD5, FADS1, and FADS2) and uptake of fatty acids (FABP3), while reducing the expression of transcription factors related to β-oxidation (CPT1A, CPT1B, and UCP3). A similar result was reported in a study on yaks, herbivores of the genus Bos, in which supplementation during the growth period stimulated de novo synthesis of fatty acids by upregulating SREBF1, ACACA, FASN, and SCD1, as well as the transcription factor H-FABP (FABP3), and downregulating CPT1 [66]. Studies with rodents concluded that increased intracellular glucose availability inhibits CPT1 due to elevated malonyl-CoA concentrations which, in turn, reduce fatty acid oxidation [67,68]. In the present study, calves fed via creep feeding ingested high amounts of starch (>500 g/day), which can be transformed into glucose. Therefore, a high amount of glucose reaching the muscle tissue of these animals may have reduced the expression of genes related to fatty acid oxidation (CPT1A, CPT1B, and UCP3).
We also observed the regulation of genes that participate in the AMPK pathway of energy metabolism. However, the upregulation of the PRKAA2 and PRKAB2 genes, which encode the catalytic subunits necessary for AMPK activation during fatty acid oxidation [69], observed in the present study in supplemented calves was contradictory. A comparison between Angus and crossbred Angus × Simmental cattle, with the former exhibiting greater marbling potential, showed lower expression of these AMPK regulatory genes [70]. However, a study using cardiac muscle cells of mice indicated that AMPK acts as a downstream signaling molecule of apelin [71]. The latter plays an important role in energy metabolism by improving insulin sensitivity [72]. Another important enzyme in energy metabolism and fatty acid metabolism regulated by AMPK is malonyl-CoA decarboxylase [73]. This mitochondrial enzyme encoded by the MLYCD gene, which is responsible for increasing fatty acid oxidation by converting malonyl-CoA to acetyl-CoA [74], was downregulated in supplemented calves (Table 5). Studies on humans [74] and mice [75] have shown that silencing of this gene in muscle tissue increases the concentration of malonyl-CoA, thus increasing the utilization of glucose while reducing fatty acid oxidation. This fact may explain the downregulation of genes related to β-oxidation (CPT1A, CPT1B, and UCP3) in supplemented calves and may represent a protection mechanism against the development of diet-induced insulin resistance [76].
Higher expression of the complement 3 (C3) gene has been demonstrated in individuals with insulin resistance or hyperinsulinemia [77,78]. In the present study, calves supplemented at 1% of body weight showed downregulation of C3 despite higher consumption of rapidly fermentable carbohydrates (starch). This finding may suggest that the insulin signaling cascade in muscle tissue is not dysregulated by high-carbohydrate intake. In addition to this reasoning, the activation of the apelin signaling pathway by upregulation of PRKAA2 and PRKAB2 and downregulation of MLYCD in supplemented calves suggests that these animals utilize more glucose as an energy source instead of fatty acids, a fact that may be closely related to changes in muscle fiber metabolism.
Skeletal muscle exhibits a certain metabolic plasticity, which allows to change substrate utilization (fat or glucose) for ATP production depending on the growth pattern of the animal [79] or nutritional stimuli [80]. During periods of accelerated muscle growth, energy expenditure for intramuscular fat and protein deposition is expected to increase [70]. Supplementation with high amounts of carbohydrates (starch) during the suckling period of calves probably stimulates glycolytic metabolism in skeletal muscle. Thus, a reduction in the fatty acid oxidation pathway for ATP production may increase the accumulation of lipids in adipose tissue. Skeletal muscle consists of a heterogeneous group of fibers that contain different myosin heavy-chain isoforms used to identify contractile and metabolic activity of muscle cells [81]. Although these isoforms were not identified in the present study, the downregulation of CPT1A, CPT1B, and UCP3 in supplemented animals suggests a lower proportion of oxidative (type I) fibers in skeletal muscle since these transcripts are abundant in mitochondria [82]. According to Gagaoua and Picard [79], glycolytic (type II β) fibers contain a smaller number of mitochondria. These fibers are more abundant in muscles of animals with accelerated growth [83].
Fatty acid transport proteins (FATPs) are a family of six isoforms (FATP1–6). The isoform encoded by the FATP1 gene (also known as SLC27A1) is highly expressed in muscle fibers, adipocytes, and hepatocytes due to the high absorption and accelerated metabolism of fatty acids in these cells [84]. In the present study, upregulation of SLC27A1 was observed in skeletal muscle collected at weaning from animals supplemented during the cow-calf phase, and these animals thus exhibited higher IMF and MS (Table 2). These results agree with other studies on cattle that showed upregulation of SLC27A1 in animals with higher intramuscular fat deposition [70,85]. Despite evidence that FATP1 (SLC27A1) expression in muscle tissue is significantly associated with lipid accumulation [70,85,86,87], the results are contradictory since some studies have shown the opposite [88,89]. This effect on lipid metabolism (oxidation or esterification) can be explained by the localization of FATP1 that varies according to cell type and is tissue-specific [84]. For example, in muscle cells, FATP1 is most abundant in mitochondria since it is the key protein for supplying energy from fatty acid oxidation [89]. On the other hand, the cytoplasm is the site of highest abundance of FATP1 in adipocytes. When stimulated with insulin, these transporters are translocated to the plasma membrane for uptake, esterification, and accumulation of lipids [90]. An increase in FATP1 in adipose tissue increases the clearance of triglycerides from the vascular bed of muscle, thus reducing the availability of this substrate to be broken down into fatty acids [91]. Considering the negative effect of fatty acids on insulin-mediated glucose metabolism, a reduction in the flow of this substrate into muscle fibers may improve insulin sensitivity [91]. Since skeletal muscle is a heterogeneous tissue composed of muscle fibers, adipocytes, and fibroblasts, among other cell groups, the upregulated SLC27A1 (FATP1) in the supplemented animals of the present study may be derived from adipocytes and not from muscle fibers. Furthermore, lower expression of SLC27A1 in muscle fibers has been shown to be related to downregulation of CPT1A. Consequently, a lower rate of fatty acid oxidation in muscle fibers leads to a greater accumulation of intramuscular fat [88], if upregulation of SLC27A1 (FATP1) in supplemented calves of the present study occurred in intramuscular adipose tissue. The supply of energy produced by the uptake and oxidation of fatty acids to muscle fibers may be reduced, a fact that partly explains the downregulation of pro-oxidative genes (CPT1A, CPT1B, and UCP3) in muscle of these animals.
In addition to stimulating lipogenesis, FATP1 (SLC27A1) participates in the differentiation of pre-adipocytes regulated by PPAR [86,92]. During adipogenesis, FATP1 increases the release of triglycerides into plasma and promotes the uptake of fatty acids by pre-adipocytes [93]. Once inside the cells, fatty acids bind to the ligand-binding domain that modifies the structure of PPAR, forming a heterodimer with the retinoid X receptor (RXR) [56]. Finally, the PPAR/RXR complex binds to the specific promoter regions of the target genes, inducing or inhibiting their expression [94]. In the present study, animals supplemented by creep feeding consumed approximately 375 g of soybean meal (Table 1). Diets containing soybean meal are rich in oleic, linoleic, and linolenic acids [58], which are potent regulators of PPARA [56]. It is therefore possible that, during the marbling window (250 days of age), the upregulation of PPARA and RXRG in supplemented animals is an indicator of greater recruitment of intramuscular adipocytes (adipogenesis).
The FOXO1 gene encodes a transcription factor of the forkhead box class O family and plays an important role in the regulation of glucose metabolism through insulin signaling, in fatty acid metabolism (oxidation), and in the recruitment of pre-adipocytes in adipogenesis [95,96]. The role of FOXO1 in adipogenesis may be either a promotor or an inhibitor depending on the stage of differentiation. Thus, in the post-mitotic stage, upregulation of FOXO1 inhibits adipogenesis by activating the transcription of p21, a cell cycle inhibitor [97]. On the other hand, during the final differentiation and maturation of adipocytes (lipid metabolism), the inhibition of FOXO1 impairs its binding to PPARγ, which leads to the formation of the PPARγ–RXRα–DNA complex in the transcription program [97]. In the present study, downregulation of FOXO1 in calves supplemented during the lactation period suggests that the intramuscular adipose tissue collected at weaning may be in the final stage of differentiation and that adipocytes were ready to start fat accumulation during the finishing phase. Within this context, a larger number of differentiated intramuscular adipocytes in supplemented animals may also support the upregulation of genes related to the biosynthesis of fatty acids and cholesterol discussed in the previous paragraphs.
The INSIG1 gene encodes an endoplasmic reticulum membrane protein that regulates glucose homeostasis and provides a negative feedback mechanism for cholesterol biosynthesis and lipogenesis [70,98]. Murine models [99] indicated that INSIG1 expression is activated relatively late by the regulation of PPARG and SREBF1 during adipogenesis and that its role in adipose tissue is to block the release of SREBP1 into the endoplasmic reticulum [100]. SREBP1 is a key protein in the control of ACACA, FASN, and SCD gene expression during de novo fatty acid synthesis [101]. However, the present findings do not support this mechanism since, despite the upregulation of INSIG1 in supplemented calves, higher abundance of FASN, SCD, and SCD5 was also observed, which would explain the higher IMF and MS (Table 2). Similarly, in a study on Angus and crossbred Angus × Simmental cattle, upregulation of INSIG1 was observed in Angus animals with higher IMF and MS [70]. This fact suggests that overexpression of INSIG1 in ruminants acts as a pro-lipogenic factor [70] and not as an anti-lipogenic gene, as demonstrated in mouse models [99].
Lipin is a protein encoded by a family of genes (LPIN1, LPIN2, and LPIN3) that play a key role in the lipogenesis and energy metabolism of adipose and muscle tissue [102,103]. Studies on transgenic mice with lipin deficiency or exclusive overexpression in muscle have shown that the upregulation of lipin reduces the utilization of fatty acids as energy, increasing the expression of lipogenic genes in adipocytes. Thus, the expression of lipin reduces muscle tissue energy expenditure and fat oxidation, slightly increasing the obesity of individuals [104,105]. In addition to the effects on fat deposition, altered lipin expression in muscle and adipose tissue affects insulin sensitivity. Lipin deficiency promotes insulin resistance, probably as a consequence of low leptin and adiponectin levels and impaired glucose absorption [106,107]. In contrast, lipin overexpression in transgenic mice was found to increase insulin sensitivity, although adipose tissue mass was doubled compared to normal mice [104]. This finding demonstrates that lipin improves the efficiency of fatty acid storage in adipocytes and predicts ectopic deposition of lipids in muscle tissue and consequent insulin impairment [105], in agreement with the results of the present study. Another study on cattle reported that starch-rich diets increased the expression of LPIN2 and that this expression was positively correlated with insulin sensitivity and therefore associated with greater adipogenesis [70].

5. Conclusions

Supplementation of F1 Angus × Nellore calves during the cow-calf (lactation) phase prepares intramuscular adipose tissue for fat deposition during the postweaning period. The contribution of fatty acids such as linoleic and linolenic acids and high amounts of starch through ingredients such as soybean meal and ground corn included in supplements via creep feeding largely affected gene expression in the LT muscle at weaning, including adipogenic and lipogenic genes. Therefore, the upregulation of these genes that participate in critical metabolic pathways and biological processes related to PPAR signaling in adipogenesis and the insulin response, biosynthesis of fatty acids, and cholesterol in lipogenesis, contributes to the long-term intramuscular fat deposition during the finishing phase.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/metabo13020160/s1. Table S1: Ingredients and chemical-bromatological composition of the feedlot diet. Table S2: Complete list of differentially expressed genes (DEGs) between groups (G1 × G2) at weaning. Table S3: Enriched KEGG pathways for upregulated genes in group 2 at weaning. Table S4: Enriched biological processes (GO terms) for upregulated genes in group 2 at weaning. Table S5: Enriched KEGG pathways for downregulated genes in group 2 at weaning. Table S6: Enriched biological processes (GO terms) for downregulated genes in group 2 at weaning. Figure S1: Visual differences in relation to the fat content between Longissimus thoracis muscle samples collected from animals of G1 and G2. Figure S2: Number of genes by functional category detected by RNA sequencing considering the application of a filter that excludes genes with a low relative read count after normalization. Figure S3: Log2 boxplot of the read count normalized by the size factor per sample in G1 (control, no creep feeding) and G2 (creep feeding).

Author Contributions

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

Funding

This research was funded by the Research Foundation of the State of São Paulo (Fundação de Amparo à Pesquisa do Estado de São Paulo—FAPESP; grant number 2019/22621-3). This study was also financed in part by CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior), Finance code 001.

Institutional Review Board Statement

The study was conducted in accordance with the Animal Use Ethics Committee of the School of Agriculture and Veterinary Sciences (FCAV), UNESP, Jaboticabal, São Paulo, Brazil (certificate number 013689/19).

Informed Consent Statement

Not applicable.

Data Availability Statement

All relevant data are presented within the paper. RNA-Seq data may be made available by contacting the corresponding author. Data is not publicly due to privacy.

Acknowledgments

The authors thank the São Paulo State University (Universidade Estadual Paulista, UNESP) for providing the infrastructure for all the stages of this research and the company Trouw Nutrition Brasil for providing part of the animals’ diet.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Scheffler, J.; McCann, M.; Greiner, S.; Jiang, H.; Hanigan, M.; Bridges, G.; Lake, S.; Gerrard, D. Early metabolic imprinting events increase marbling scores in fed cattle. J. Anim. Sci. 2014, 92, 320–324. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Dantas, C.C.O.; de Mattos Negrão, F.; Geron, L.J.V.; Mexia, A.A. O uso da técnica do Creep-feeding na suplementação de bezerros. PUBVET 2010, 4, 899–904. [Google Scholar]
  3. Stewart, L. Creep feeding beef calves. UGA Coop. Ext. Bull. 2017, 1315, 1–8. [Google Scholar]
  4. Zhang, N. Epigenetic modulation of DNA methylation by nutrition and its mechanisms in animals. Anim. Nutr. 2015, 1, 144–151. [Google Scholar] [CrossRef]
  5. McKay, J.A.; Mathers, J.C. Diet induced epigenetic changes and their implications for health. Acta Physiol. 2011, 202, 103–118. [Google Scholar] [CrossRef]
  6. Farias, N.; Ho, N.; Butler, S.; Delaney, L.; Morrison, J.; Shahrzad, S.; Coomber, B.L. The effects of folic acid on global DNA methylation and colonosphere formation in colon cancer cell lines. J. Nutr. Biochem. 2015, 26, 818–826. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Farkas, S.A.; Befekadu, R.; Hahn-Strömberg, V.; Nilsson, T.K. DNA methylation and expression of the folate transporter genes in colorectal cancer. Tumor Biol. 2015, 36, 5581–5590. [Google Scholar] [CrossRef] [PubMed]
  8. Day, J.J.; Kennedy, A.J.; Sweatt, J.D. DNA methylation and its implications and accessibility for neuropsychiatric therapeutics. Annu. Rev. Pharmacol. Toxicol. 2015, 55, 591. [Google Scholar] [CrossRef] [Green Version]
  9. Vucetic, Z.; Kimmel, J.; Totoki, K.; Hollenbeck, E.; Reyes, T.M. Maternal high-fat diet alters methylation and gene expression of dopamine and opioid-related genes. Endocrinology 2010, 151, 4756–4764. [Google Scholar] [CrossRef] [Green Version]
  10. Bogdarina, I.; Haase, A.; Langley-Evans, S.; Clark, A.J. Glucocorticoid effects on the programming of AT1b angiotensin receptor gene methylation and expression in the rat. PLoS ONE 2010, 5, e9237. [Google Scholar] [CrossRef] [Green Version]
  11. Jousse, C.; Parry, L.; Lambert-Langlais, S.; Maurin, A.C.; Averous, J.; Bruhat, A.; Carraro, V.; Tost, J.; Letteron, P.; Chen, P. Perinatal undernutrition affects the methylation and expression of the leptin gene in adults: Implication for the understanding of metabolic syndrome. FASEB J. 2011, 25, 3271–3278. [Google Scholar] [CrossRef] [PubMed]
  12. Dudley, K.J.; Sloboda, D.M.; Connor, K.L.; Beltrand, J.; Vickers, M.H. Offspring of mothers fed a high fat diet display hepatic cell cycle inhibition and associated changes in gene expression and DNA methylation. PLoS ONE 2011, 6, e21662. [Google Scholar] [CrossRef] [PubMed]
  13. Altmann, S.; Murani, E.; Schwerin, M.; Metges, C.C.; Wimmers, K.; Ponsuksili, S. Somatic cytochrome c (CYCS) gene expression and promoter-specific DNA methylation in a porcine model of prenatal exposure to maternal dietary protein excess and restriction. Br. J. Nutr. 2012, 107, 791–799. [Google Scholar] [CrossRef] [Green Version]
  14. Smith, S.B.; Crouse, J.D. Relative contributions of acetate, lactate and glucose to lipogenesis in bovine intramuscular and subcutaneous adipose tissue. J. Nutr. 1984, 114, 792–800. [Google Scholar] [CrossRef]
  15. Baldassini, W.A.; Ferreira, M.S.; Santiago, B.M.; Chardulo, L.A.L.; Curi, R.A.; Lanna, D.P.D.; Ribeiro, R.V.; Martins, C.L.; Arrigoni, M.; Neto, O.R.M. Intake, performance, meat quality and fatty acid profile of Nellore bulls finished in feedlot with diets containing dry corn gluten feed. Livest. Sci. 2021, 253, 104715. [Google Scholar] [CrossRef]
  16. Lopez, R.; Krehbiel, C.R.; Thomas, M.G.; Hallford, D.M.; Keisler, D.H.; Obeidat, B.S.; Hernandez, J.A.; Bryant, W.D.; Garcia, M.; Flores, R. Effect of increasing level of dietary protein on serum concentrations of metabolic hormones and mammary development in Holstein heifers consuming a moderate-energy diet. J. Dairy Sci. 2001, 84, 161. [Google Scholar]
  17. ABIEC. Beef Report. 2022. Available online: https://www.abiec.com.br/publicacoes/beef-report-2022/ (accessed on 19 September 2022).
  18. Hocquette, J.-F.; Gondret, F.; Baeza, E.; Medale, F.; Jurie, C.; Pethick, D. Intramuscular fat content in meat-producing animals: Development, genetic and nutritional control, and identification of putative markers. Animal 2010, 4, 303–319. [Google Scholar] [CrossRef] [Green Version]
  19. Du, M.; Ford, S.P.; Zhu, M.-J. Optimizing livestock production efficiency through maternal nutritional management and fetal developmental programming. Anim. Front. 2017, 7, 5–11. [Google Scholar] [CrossRef] [Green Version]
  20. Lanna, D.; Barioni, L.; Nepomuceno, N.; Almeida, R.; Tedeschi, L. Ração de Lucro Máximo—RLM, Version 3.2.1; RLM: Piracicaba, Brazil, 2011.
  21. AUS-MEAT Limited. Handbook of Australian Beef Processing; AUS-MEAT Limited: Murarrie, Australia, 2018. [Google Scholar]
  22. Shackelford, S.; Wheeler, T.; Koohmaraie, M. Evaluation of slice shear force as an objective method of assessing beef longissimus tenderness. J. Anim. Sci. 1999, 77, 2693–2699. [Google Scholar] [CrossRef] [Green Version]
  23. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria; Available online: https://www.R-project.org/ (accessed on 2 June 2022).
  24. Saini, H.K.; Griffiths-Jones, S.; Enright, A.J. Genomic analysis of human microRNA transcripts. Proc. Natl. Acad. Sci. USA 2007, 104, 17719–17724. [Google Scholar] [CrossRef] [Green Version]
  25. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 2 June 2022).
  26. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. Fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef] [PubMed]
  27. Ewels, P.; Magnusson, M.; Lundin, S.; Käller, M. MultiQC: Summarize analysis results for multiple tools and samples in a single report. Bioinformatics 2016, 32, 3047–3048. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  29. Liao, Y.; Smyth, G.K.; Shi, W. FeatureCounts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014, 30, 923–930. [Google Scholar] [CrossRef] [Green Version]
  30. Kassambara, A.; Mundt, F. Factoextra: Extract and Visualize the Results of Multivariate Data Analyses. Available online: https://CRAN.R-project.org/package=factoextra (accessed on 2 June 2022).
  31. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. EdgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
  32. Robinson, M.D.; Oshlack, A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010, 11, R25. [Google Scholar] [CrossRef] [Green Version]
  33. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  34. Huang, D.W.; Sherman, B.T.; Tan, Q.; Kir, J.; Liu, D.; Bryant, D.; Guo, Y.; Stephens, R.; Baseler, M.W.; Lane, H.C. DAVID Bioinformatics Resources: Expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res. 2007, 35, W169–W175. [Google Scholar] [CrossRef] [Green Version]
  35. Bindea, G.; Mlecnik, B.; Hackl, H.; Charoentong, P.; Tosolini, M.; Kirilovsky, A.; Fridman, W.-H.; Pagès, F.; Trajanoski, Z.; Galon, J. ClueGO: A Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 2009, 25, 1091–1093. [Google Scholar] [CrossRef] [Green Version]
  36. Gu, Z.; Eils, R.; Schlesner, M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics 2016, 32, 2847–2849. [Google Scholar] [CrossRef] [Green Version]
  37. Bindea, G.; Galon, J.; Mlecnik, B. CluePedia Cytoscape plugin: Pathway insights using integrated experimental and in silico data. Bioinformatics 2013, 29, 661–663. [Google Scholar] [CrossRef] [PubMed]
  38. Szklarczyk, D.; Morris, J.H.; Cook, H.; Kuhn, M.; Wyder, S.; Simonovic, M.; Santos, A.; Doncheva, N.T.; Roth, A.; Bork, P. The STRING database in 2017: Quality-controlled protein–protein association networks, made broadly accessible. Nucleic Acids Res. 2016, 45, gkw937. [Google Scholar] [CrossRef]
  39. Costa e Silva, L.F.; Engle, T.E.; de Valadares Filho, S.; Rotta, P.P.; Villadiego, F.A.C.; Silva, F.A.S.; Martins, E.C.; Silva, L.H.R.; Paulino, M.F. Nellore cows and their calves during the lactation period: Performance, intake, milk composition, and total apparent digestibility. Trop. Anim. Health Prod. 2015, 47, 735–741. [Google Scholar] [CrossRef] [PubMed]
  40. Fonseca, M.A.; Valadares Filho, S.D.C.; Henriques, L.T.; Paulino, P.V.R.; Detmann, E.; Fonseca, E.A.; Benedeti, P.D.B.; Silva, L.D.D. Nutritional requirements of nursing Nellore calves. Rev. Bras. Zootec. 2012, 41, 1212–1221. [Google Scholar] [CrossRef] [Green Version]
  41. Carvalho, V.V.; Paulino, M.F.; Detmann, E.; Valadares Filho, S.C.; Lopes, S.A.; Rennó, L.N.; Sampaio, C.B.; Silva, A.G. A meta-analysis of the effects of creep-feeding supplementation on performance and nutritional characteristics by beef calves grazing on tropical pastures. Livest. Sci. 2019, 227, 175–182. [Google Scholar] [CrossRef]
  42. Cremin, J., Jr.; Faulkner, D.; Merchen, N.; Fahey, G., Jr.; Fernando, R.; Willms, C. Digestion criteria in nursing beef calves supplemented with limited levels of protein and energy. J. Anim. Sci. 1991, 69, 1322–1331. [Google Scholar] [CrossRef] [Green Version]
  43. Faulkner, D.; Hummel, D.; Buskirk, D.; Berger, L.; Parrett, D.; Cmarik, G. Performance and nutrient metabolism by nursing calves supplemented with limited or unlimited corn or soyhulls. J. Anim. Sci. 1994, 72, 470–477. [Google Scholar] [CrossRef] [Green Version]
  44. Lopes, S.A.; Paulino, M.F.; Detmann, E.; Valente, É.E.L.; Rennó, L.N.; Valadares, R.F.D.; Cardenas, J.E.G.; de Almeida, D.M.; de Moura, F.H.; Oliveira, C.A.S. Evaluation of supplementation plans for suckling beef calves managed on tropical pasture. Semin. Ciênc. Agrár. 2017, 38, 1027–1039. [Google Scholar] [CrossRef] [Green Version]
  45. Asher, A.; Shabtay, A.; Cohen-Zinder, M.; Aharoni, Y.; Miron, J.; Agmon, R.; Halachmi, I.; Orlov, A.; Haim, A.; Tedeschi, L. Consistency of feed efficiency ranking and mechanisms associated with inter-animal variation among growing calves. J. Anim. Sci. 2018, 96, 990–1009. [Google Scholar] [CrossRef]
  46. Gerrits, W.J.; France, J.; Dijkstra, J.; Bosch, M.W.; Tolman, G.H.; Tamminga, S. Evaluation of a model integrating protein and energy metabolism in preruminant calves. J. Nutr. 1997, 127, 1243–1252. [Google Scholar] [CrossRef] [Green Version]
  47. Hunt, M.; Garmyn, A.; O’Quinn, T.; Corbin, C.; Legako, J.; Rathmann, R.; Brooks, J.; Miller, M. Consumer assessment of beef palatability from four beef muscles from USDA Choice and Select graded carcasses. Meat Sci. 2014, 98, 1–8. [Google Scholar] [CrossRef]
  48. Du, M.; Huang, Y.; Das, A.; Yang, Q.; Duarte, M.; Dodson, M.; Zhu, M.-J. Meat Science and Muscle Biology Symposium: Manipulating mesenchymal progenitor cell differentiation to optimize performance and carcass value of beef cattle. J. Anim. Sci. 2013, 91, 1419–1427. [Google Scholar] [CrossRef] [PubMed]
  49. Gross, B.; Pawlak, M.; Lefebvre, P.; Staels, B. PPARs in obesity-induced T2DM, dyslipidaemia and NAFLD. Nat. Rev. Endocrinol. 2017, 13, 36–49. [Google Scholar] [CrossRef] [PubMed]
  50. Chen, H.; Peng, T.; Shang, H.; Shang, X.; Zhao, X.; Qu, M.; Song, X. RNA-Seq Analysis Reveals the Potential Molecular Mechanisms of Puerarin on Intramuscular Fat Deposition in Heat-Stressed Beef Cattle. Front. Nutr. 2022, 9, 817557. [Google Scholar] [CrossRef]
  51. Brown, J.D.; Plutzky, J. Peroxisome proliferator–activated receptors as transcriptional nodal points and therapeutic targets. Circulation 2007, 115, 518–533. [Google Scholar] [CrossRef] [PubMed]
  52. Rangwala, S.M.; Lazar, M.A. Peroxisome proliferator-activated receptor γ in diabetes and metabolism. Trends Pharmacol. Sci. 2004, 25, 331–336. [Google Scholar] [CrossRef]
  53. Ladeira, M.; Schoonmaker, J.; Gionbelli, M.; Dias, J.; Gionbelli, T.; Carvalho, J.R.; Teixeira, P. Nutrigenomics and beef quality: A review about lipogenesis. Int. J. Mol. Sci. 2016, 17, 918. [Google Scholar] [CrossRef] [Green Version]
  54. Ward, R.; Woodward, B.; Otter, N.; Doran, O. Relationship between the expression of key lipogenic enzymes, fatty acid composition, and intramuscular fat content of Limousin and Aberdeen Angus cattle. Livest. Sci. 2010, 127, 22–29. [Google Scholar] [CrossRef]
  55. Nakamura, M.T.; Nara, T.Y. Structure, function, and dietary regulation of (delta) 6,(delta) 5, and (delta) 9 desaturases. Annu. Rev. Nutr. 2004, 24, 345. [Google Scholar] [CrossRef]
  56. Bionaz, M.; Chen, S.; Khan, M.J.; Loor, J.J. Functional Role of PPARs in Ruminants: Potential Targets for Fine-Tuning Metabolism during Growth and Lactation. PPAR Res. 2013, 2013, 28. [Google Scholar] [CrossRef] [Green Version]
  57. Ibeagha-Awemu, E.M.; Akwanji, K.A.; Beaudoin, F.; Zhao, X. Associations between variants of FADS genes and omega-3 and omega-6 milk fatty acids of Canadian Holstein cows. BMC Genet. 2014, 15, 25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Ladeira, M.M.; de Oliveira, D.M.; Schoonmaker, J.P.; Chizzotti, M.L.; Barreto, H.G.; Paiva, L.V.; Coelho, T.C.; Neto, O.R.M.; Gionbelli, M.P.; Chalfun-Junior, A. Expression of lipogenic genes in the muscle of beef cattle fed oilseeds and vitamin E. Agri Gene 2020, 15, 100097. [Google Scholar] [CrossRef]
  59. Duckett, S.; Gillis, M. Effects of oil source and fish oil addition on ruminal biohydrogenation of fatty acids and conjugated linoleic acid formation in beef steers fed finishing diets. J. Anim. Sci. 2010, 88, 2684–2691. [Google Scholar] [CrossRef] [PubMed]
  60. Schroyen, M.; Li, B.; Arevalo Sureda, E.; Zhang, Y.; Leblois, J.; Deforce, D.; Van Nieuwerburgh, F.; Wavreille, J.; Everaert, N. Pre-Weaning Inulin Supplementation Alters the Ileal Transcriptome in Pigs Regarding Lipid Metabolism. Vet. Sci. 2021, 8, 207. [Google Scholar] [CrossRef] [PubMed]
  61. Mangaraj, M.; Nanda, R.; Panda, S. Apolipoprotein AI: A molecule of diverse function. Indian J. Clin. Biochem. 2016, 31, 253–259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Poudyal, H.; Panchal, S.K.; Diwan, V.; Brown, L. Omega-3 fatty acids and metabolic syndrome: Effects and emerging mechanisms of action. Prog. Lipid Res. 2011, 50, 372–387. [Google Scholar] [CrossRef]
  63. Fox, J.; McGill Jr, H.; Carey, K.; Getz, G. In vivo regulation of hepatic LDL receptor mRNA in the baboon. Differential effects of saturated and unsaturated fat. J. Biol. Chem. 1987, 262, 7014–7020. [Google Scholar] [CrossRef]
  64. Blædel, T.; Holm, J.B.; Sundekilde, U.K.; Schmedes, M.S.; Hess, A.L.; Lorenzen, J.K.; Kristiansen, K.; Dalsgaard, T.K.; Astrup, A.; Larsen, L.H. A randomised, controlled, crossover study of the effect of diet on angiopoietin-like protein 4 (ANGPTL4) through modification of the gut microbiome. J. Nutr. Sci. 2016, 5, E45. [Google Scholar] [CrossRef] [Green Version]
  65. Yoo, J.Y.; Kim, S.S. Probiotics and prebiotics: Present status and future perspectives on metabolic disorders. Nutrients 2016, 8, 173. [Google Scholar] [CrossRef] [Green Version]
  66. Zhang, H.B.; Wang, Z.S.; Peng, Q.H.; Tan, C.; Zou, H.W. Effects of different levels of protein supplementary diet on gene expressions related to intramuscular deposition in early-weaned yaks. Anim. Sci. J. 2014, 85, 411–419. [Google Scholar] [CrossRef]
  67. Turcotte, L.P.; Swenberger, J.R.; Yee, A.J. High carbohydrate availability increases LCFA uptake and decreases LCFA oxidation in perfused muscle. Am. J. Physiol. Endocrinol. Metab. 2002, 282, E177–E183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Yee, A.J.; Turcotte, L.P. Insulin fails to alter plasma LCFA metabolism in muscle perfused at similar glucose uptake. Am. J. Physiol. Endocrinol. Metab. 2002, 283, E73–E77. [Google Scholar] [CrossRef] [Green Version]
  69. Zhang, M.; Zheng, D.; Peng, Z.; Zhu, Y.; Li, R.; Wu, Q.; Li, Y.; Li, H.; Xu, W.; Zhang, M. Identification of differentially expressed genes and lipid metabolism signaling pathways between muscle and fat tissues in broiler chickens. J. Poult. Sci. 2021, 58, 131–137. [Google Scholar] [CrossRef] [PubMed]
  70. Graugnard, D.E.; Piantoni, P.; Bionaz, M.; Berger, L.L.; Faulkner, D.B.; Loor, J.J. Adipogenic and energy metabolism gene networks in longissimus lumborum during rapid post-weaning growth in Angus and Angus× Simmental cattle fed high-starch or low-starch diets. BMC Genom. 2009, 10, 142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Simpkin, J.C.; Yellon, D.M.; Davidson, S.M.; Lim, S.Y.; Wynne, A.M.; Smith, C.C. Apelin-13 and apelin-36 exhibit direct cardioprotective activity against ischemiareperfusion injury. Basic Res. Cardiol. 2007, 102, 518–528. [Google Scholar] [CrossRef]
  72. Castan-Laurell, I.; Dray, C.; Knauf, C.; Kunduzova, O.; Valet, P. Apelin, a promising target for type 2 diabetes treatment? Trends Endocrinol. Metab. 2012, 23, 234–241. [Google Scholar] [CrossRef]
  73. Ruderman, N.; Prentki, M. AMP kinase and malonyl-CoA: Targets for therapy of the metabolic syndrome. Nat. Rev. Drug Discov. 2004, 3, 340–351. [Google Scholar] [CrossRef]
  74. Bouzakri, K.; Austin, R.; Rune, A.; Lassman, M.E.; Garcia-Roves, P.M.; Berger, J.P.; Krook, A.; Chibalin, A.V.; Zhang, B.B.; Zierath, J.R. Malonyl CoenzymeA decarboxylase regulates lipid and glucose metabolism in human skeletal muscle. Diabetes 2008, 57, 1508–1516. [Google Scholar] [CrossRef] [Green Version]
  75. Dyck, J.R.; Hopkins, T.A.; Bonnet, S.; Michelakis, E.D.; Young, M.E.; Watanabe, M.; Kawase, Y.; Jishage, K.-I.; Lopaschuk, G.D. Absence of malonyl coenzyme A decarboxylase in mice increases cardiac glucose oxidation and protects the heart from ischemic injury. Circulation 2006, 114, 1721–1728. [Google Scholar] [CrossRef] [Green Version]
  76. Koves, T.R.; Ussher, J.R.; Noland, R.C.; Slentz, D.; Mosedale, M.; Ilkayeva, O.; Bain, J.; Stevens, R.; Dyck, J.R.; Newgard, C.B. Mitochondrial overload and incomplete fatty acid oxidation contribute to skeletal muscle insulin resistance. Cell Metab. 2008, 7, 45–56. [Google Scholar] [CrossRef] [Green Version]
  77. Weyer, C.; Tataranni, P.A.; Pratley, R.E. Insulin action and insulinemia are closely related to the fasting complement C3, but not acylation stimulating protein concentration. Diabetes Care 2000, 23, 779–785. [Google Scholar] [CrossRef] [Green Version]
  78. Barbu, A.; Hamad, O.A.; Lind, L.; Ekdahl, K.N.; Nilsson, B. The role of complement factor C3 in lipid metabolism. Mol. Immunol. 2015, 67, 101–107. [Google Scholar] [CrossRef] [PubMed]
  79. Picard, B.; Gagaoua, M. Muscle fiber properties in cattle and their relationships with meat qualities: An overview. J. Agric. Food Chem. 2020, 68, 6021–6039. [Google Scholar] [CrossRef] [PubMed]
  80. Aragão, R.D.S.; Guzmán-Quevedo, O.; Pérez-García, G.; Manhaes-de-Castro, R.; Bolanos-Jiménez, F. Maternal protein restriction impairs the transcriptional metabolic flexibility of skeletal muscle in adult rat offspring. Br. J. Nutr. 2014, 112, 328–337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Schiaffino, S.; Reggiani, C. Fiber types in mammalian skeletal muscles. Physiol. Rev. 2011, 91, 1447–1531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Turner, N.; Hariharan, K.; TidAng, J.; Frangioudakis, G.; Beale, S.M.; Wright, L.E.; Zeng, X.Y.; Leslie, S.J.; Li, J.-Y.; Kraegen, E.W. Enhancement of muscle mitochondrial oxidative capacity and alterations in insulin action are lipid species dependent: Potent tissue-specific effects of medium-chain fatty acids. Diabetes 2009, 58, 2547–2554. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Brandstetter, A.M.; Picard, B.; Geay, Y. Muscle fibre characteristics in four muscles of growing bulls: I. Postnatal differentiation. Livest. Prod. Sci. 1998, 53, 15–23. [Google Scholar] [CrossRef]
  84. Huang, J.; Zhu, R.; Shi, D. The role of FATP1 in lipid accumulation: A review. Mol. Cell. Biochem. 2021, 476, 1897–1903. [Google Scholar] [CrossRef]
  85. Jeong, J.; Kwon, E.G.; Im, S.K.; Seo, K.S.; Baik, M. Expression of fat deposition and fat removal genes is associated with intramuscular fat content in longissimus dorsi muscle of Korean cattle steers. J. Anim. Sci. 2012, 90, 2044–2053. [Google Scholar] [CrossRef]
  86. Chen, X.; Luo, Y.; Wang, R.; Zhou, B.; Huang, Z.; Jia, G.; Zhao, H.; Liu, G. Effects of fatty acid transport protein 1 on proliferation and differentiation of porcine intramuscular preadipocytes. Anim. Sci. J. 2017, 88, 731–738. [Google Scholar] [CrossRef]
  87. Qi, R.; Long, D.; Wang, J.; Wang, Q.; Huang, X.; Cao, C.; Gao, G.; Huang, J. MicroRNA-199a targets the fatty acid transport protein 1 gene and inhibits the adipogenic trans-differentiation of C2C12 myoblasts. Cell. Physiol. Biochem. 2016, 39, 1087–1097. [Google Scholar] [CrossRef] [PubMed]
  88. Qiu, F.; Xie, L.; Ma, J.-E.; Luo, W.; Zhang, L.; Chao, Z.; Chen, S.; Nie, Q.; Lin, Z.; Zhang, X. Lower expression of SLC27A1 enhances intramuscular fat deposition in chicken via down-regulated fatty acid oxidation mediated by CPT1A. Front. Physiol. 2017, 8, 449. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Guitart, M.; Osorio-Conles, Ó.; Pentinat, T.; Cebrià, J.; García-Villoria, J.; Sala, D.; Sebastián, D.; Zorzano, A.; Ribes, A.; Jiménez-Chillarón, J.C. Fatty acid transport protein 1 (FATP1) localizes in mitochondria in mouse skeletal muscle and regulates lipid and ketone body disposal. PLoS ONE 2014, 9, e98109. [Google Scholar] [CrossRef]
  90. Lobo, S.; Wiczer, B.M.; Smith, A.J.; Hall, A.M.; Bernlohr, D.A. Fatty acid metabolism in adipocytes: Functional analysis of fatty acid transport proteins 1 and 4. J. Lipid Res. 2007, 48, 609–620. [Google Scholar] [CrossRef] [Green Version]
  91. Martin, G.; Schoonjans, K.; Lefebvre, A.-M.; Staels, B.; Auwerx, J. Coordinate regulation of the expression of the fatty acid transport protein and acyl-CoA synthetase genes by PPARα and PPARγ activators. J. Biol. Chem. 1997, 272, 28210–28217. [Google Scholar] [CrossRef] [Green Version]
  92. Qi, R.; Feng, M.; Tan, X.; Gan, L.; Yan, G.; Sun, C. FATP1 silence inhibits the differentiation and induces the apoptosis in chicken preadipocytes. Mol. Biol. Rep. 2013, 40, 2907–2914. [Google Scholar] [CrossRef]
  93. Liu, X.; Li, S.; Wang, L.; Zhang, W.; Wang, Y.; Gui, L.; Zan, L.; Zhao, C. The effect of FATP1 on adipocyte differentiation in Qinchuan beef cattle. Animals 2021, 11, 2789. [Google Scholar] [CrossRef]
  94. La Cour Poulsen, L.; Siersbæk, M.; Mandrup, S. PPARs: Fatty acid sensors controlling metabolism. Semin. Cell. Dev. Biol. 2012, 23, 631–639. [Google Scholar] [CrossRef] [PubMed]
  95. Nakae, J.; Kitamura, T.; Kitamura, Y.; Biggs III, W.H.; Arden, K.C.; Accili, D. The forkhead transcription factor Foxo1 regulates adipocyte differentiation. Dev. Cell 2003, 4, 119–129. [Google Scholar] [CrossRef] [Green Version]
  96. Ioannilli, L.; Ciccarone, F.; Ciriolo, M.R. Adipose tissue and FoxO1: Bridging physiology and mechanisms. Cells 2020, 9, 849. [Google Scholar] [CrossRef] [Green Version]
  97. Dowell, P.; Otto, T.C.; Adi, S.; Lane, M.D. Convergence of peroxisome proliferator-activated receptor γ and Foxo1 signaling pathways. J. Biol. Chem. 2003, 278, 45485–45491. [Google Scholar] [CrossRef] [Green Version]
  98. Attie, A.D. Insig: A significant integrator of nutrient and hormonal signals. J. Clin. Investig. 2004, 113, 1112–1114. [Google Scholar] [CrossRef] [PubMed]
  99. Kast-Woelbern, H.R.; Dana, S.L.; Cesario, R.M.; Sun, L.; de Grandpre, L.Y.; Brooks, M.E.; Osburn, D.L.; Reifel-Miller, A.; Klausing, K.; Leibowitz, M.D. Rosiglitazone induction of Insig-1 in white adipose tissue reveals a novel interplay of peroxisome proliferator-activated receptor γ and sterol regulatory element-binding protein in the regulation of adipogenesis. J. Biol. Chem. 2004, 279, 23908–23915. [Google Scholar] [CrossRef] [PubMed]
  100. Li, J.; Takaishi, K.; Cook, W.; McCorkle, S.K.; Unger, R.H. Insig-1 “brakes” lipogenesis in adipocytes and inhibits differentiation of preadipocytes. Proc. Natl. Acad. Sci. USA 2003, 100, 9476–9481. [Google Scholar] [CrossRef] [Green Version]
  101. Teixeira, P.D.; Oliveira, D.M.; Chizzotti, M.L.; Chalfun-Junior, A.; Coelho, T.C.; Gionbelli, M.; Paiva, L.V.; Carvalho, J.R.R.; Ladeira, M.M. Subspecies and diet affect the expression of genes involved in lipid metabolism and chemical composition of muscle in beef cattle. Meat Sci. 2017, 133, 110–118. [Google Scholar] [CrossRef] [PubMed]
  102. Reue, K. The lipin family: Mutations and metabolism. Curr. Opin. Lipidol. 2009, 20, 165. [Google Scholar] [CrossRef]
  103. Zhang, P.; Reue, K. Lipin proteins and glycerolipid metabolism: Roles at the ER membrane and beyond. Biochim. Biophys. Acta BBA Biomembr. 2017, 1859, 1583–1595. [Google Scholar] [CrossRef]
  104. Phan, J.; Reue, K. Lipin, a lipodystrophy and obesity gene. Cell Metab. 2005, 1, 73–83. [Google Scholar] [CrossRef] [Green Version]
  105. Reue, K.; Donkor, J. Lipin: A determinant of adiposity, insulin sensitivity and energy balance. Future Lipidol. 2006, 1, 91–101. [Google Scholar] [CrossRef]
  106. Péterfy, M.; Phan, J.; Xu, P.; Reue, K. Lipodystrophy in the fld mouse results from mutation of a new gene encoding a nuclear protein, lipin. Nat. Genet. 2001, 27, 121–124. [Google Scholar] [CrossRef]
  107. Reue, K.; Xu, P.; Wang, X.-P.; Slavin, B.G. Adipose tissue deficiency, glucose intolerance, and increased atherosclerosis result from mutation in the mouse fatty liver dystrophy (fld) gene. J. Lipid Res. 2000, 41, 1067–1076. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Expression profile of reference genes in the control group (G1, no creep feeding) and the group submitted to creep feeding (G2).
Figure 1. Expression profile of reference genes in the control group (G1, no creep feeding) and the group submitted to creep feeding (G2).
Metabolites 13 00160 g001
Figure 2. Principal component (PC) analysis performed based on normalized count data of gene expression for samples collected at weaning from G1 (control, no creep feeding) and G2 (creep feeding).
Figure 2. Principal component (PC) analysis performed based on normalized count data of gene expression for samples collected at weaning from G1 (control, no creep feeding) and G2 (creep feeding).
Metabolites 13 00160 g002
Figure 3. Volcano plot of log2 fold change (x-axis) versus -log10 p-value (FDR, y-axis) indicating differences in DEGs identified by the edgeR method (downregulated: log2 fold change < −0.5 and FDR < 0.05; upregulated: log2 fold change > 0.5 and FDR < 0.05) between G1 (control) × G2 (creep feeding) at weaning.
Figure 3. Volcano plot of log2 fold change (x-axis) versus -log10 p-value (FDR, y-axis) indicating differences in DEGs identified by the edgeR method (downregulated: log2 fold change < −0.5 and FDR < 0.05; upregulated: log2 fold change > 0.5 and FDR < 0.05) between G1 (control) × G2 (creep feeding) at weaning.
Metabolites 13 00160 g003
Figure 4. Relationship between up- and down-regulated DEGs identified in G2 at weaning (x-axis) and enriched biological processes and pathways, which are highlighted due to their relationships with each other and possible relationships with intramuscular fat content in cattle (y-axis). The color scale indicating the log2 fold change is given on the right side of the figure.
Figure 4. Relationship between up- and down-regulated DEGs identified in G2 at weaning (x-axis) and enriched biological processes and pathways, which are highlighted due to their relationships with each other and possible relationships with intramuscular fat content in cattle (y-axis). The color scale indicating the log2 fold change is given on the right side of the figure.
Metabolites 13 00160 g004
Figure 5. Circular network of shared upregulated DEGs in G2 at weaning between biological processes and KEGG pathways, which were highlighted in ORA and due to their potential relationship with intramuscular fat in cattle.
Figure 5. Circular network of shared upregulated DEGs in G2 at weaning between biological processes and KEGG pathways, which were highlighted in ORA and due to their potential relationship with intramuscular fat in cattle.
Metabolites 13 00160 g005
Figure 6. Circular network of shared downregulated DEGs in G2 at weaning between biological processes and KEGG pathways, which were highlighted in ORA and due to their potential relationship with intramuscular fat in cattle.
Figure 6. Circular network of shared downregulated DEGs in G2 at weaning between biological processes and KEGG pathways, which were highlighted in ORA and due to their potential relationship with intramuscular fat in cattle.
Metabolites 13 00160 g006
Figure 7. Protein–protein interaction network of DEGs upregulated in G2 that participate in the biological processes and KEGG pathways identified as enriched and highlighted in the G1 × G2 comparison at weaning. Type of interaction (edge color) between nodes (DEGs): light green—text mining; dark green—neighborhood; light blue—curated databases; dark blue—co-occurrence; pink—experiments; purple—protein homology; black—co-expression. A larger number of edges indicate stronger evidence/greater strength of the interaction between two nodes. The nodes that make up the clusters formed (joined by dotted edges) are illustrated with different colors.
Figure 7. Protein–protein interaction network of DEGs upregulated in G2 that participate in the biological processes and KEGG pathways identified as enriched and highlighted in the G1 × G2 comparison at weaning. Type of interaction (edge color) between nodes (DEGs): light green—text mining; dark green—neighborhood; light blue—curated databases; dark blue—co-occurrence; pink—experiments; purple—protein homology; black—co-expression. A larger number of edges indicate stronger evidence/greater strength of the interaction between two nodes. The nodes that make up the clusters formed (joined by dotted edges) are illustrated with different colors.
Metabolites 13 00160 g007
Figure 8. Protein–protein interaction network of DEGs downregulated in G2 that participate in the biological processes and KEGG pathways identified as enriched and highlighted in the G1 × G2 comparison at weaning. Type of interaction (edge color) between nodes (DEGs): light green—text mining; dark green—neighborhood; light blue—curated databases; dark blue—co-occurrence; pink—experiments; purple—protein homology; black—co-expression. A larger number of edges indicate stronger evidence/greater strength of the interaction between two nodes. The nodes that make up the clusters formed (joined by dotted edges) are illustrated with different colors.
Figure 8. Protein–protein interaction network of DEGs downregulated in G2 that participate in the biological processes and KEGG pathways identified as enriched and highlighted in the G1 × G2 comparison at weaning. Type of interaction (edge color) between nodes (DEGs): light green—text mining; dark green—neighborhood; light blue—curated databases; dark blue—co-occurrence; pink—experiments; purple—protein homology; black—co-expression. A larger number of edges indicate stronger evidence/greater strength of the interaction between two nodes. The nodes that make up the clusters formed (joined by dotted edges) are illustrated with different colors.
Metabolites 13 00160 g008
Table 1. Means and respective standard errors of body weight at the beginning of the cow-calf phase, body weight at weaning, body weight at the end of the feedlot period, pre- and post-weaning average daily gain, and hot carcass yield obtained for the two treatments.
Table 1. Means and respective standard errors of body weight at the beginning of the cow-calf phase, body weight at weaning, body weight at the end of the feedlot period, pre- and post-weaning average daily gain, and hot carcass yield obtained for the two treatments.
BWi (kg)WW (kg)ADG1 (kg)BWf (kg)ADG2 (kg)HCW (kg)
G1 161.29 ± 2.41228.92 ± 5.07 b0.93 ± 0.02 b484.64 ± 5.961.36 ± 0.02269.22 ± 8.23
G2 157.55 ± 2.61243.57 ± 5.70 a1.03 ± 0.03 a491.85 ± 7.851.32 ± 0.03273.62 ± 9.28
1 G1: no creep feeding, G2: creep feeding, BWi: initial body weight, WW: weaning weight, ADG1: average daily gain between WW and BWi, BWf: final body weight, ADG2: average daily gain between BWf and WW, HCW: hot carcass weight. a,b Means followed by different superscript letters significantly differ (p < 0.05).
Table 2. Means and respective standard errors of carcass backfat thickness and rib eye area, fat percentage, marbling index, and shear force at 7 and 14 days of aging obtained for the two treatments.
Table 2. Means and respective standard errors of carcass backfat thickness and rib eye area, fat percentage, marbling index, and shear force at 7 and 14 days of aging obtained for the two treatments.
BFT (mm)IMF (%)MSREA (cm2)WBSF7 (kg)WBSF14 (kg)
G1 110.61 ± 0.42 b4.95 ± 0.20 b321.50 ± 13.65 b67.94 ± 1.164.52 ± 0.113.45 ± 0.11
G2 112.96 ± 0.83 a5.80 ± 0.23 a366.11 ± 12.39 a65.50 ± 0.934.28 ± 0.123.42 ± 0.09
1 G1: group 1, G2: group 2, BFT: backfat thickness, IMF: fat percentage, MS: marbling score, REA: rib eye area, WBSF7 and WBSF14: Warner–Bratzler shear force at 7 and 14 days post-mortem, respectively. a,b Means followed by different superscript letters significantly differ (p < 0.05).
Table 3. Total number of generated PE reads aligned to the reference genome, and total number and percentage of uniquely mapped PE reads in samples collected from G1 animals (no creep feeding) at weaning.
Table 3. Total number of generated PE reads aligned to the reference genome, and total number and percentage of uniquely mapped PE reads in samples collected from G1 animals (no creep feeding) at weaning.
Animal/SampleNo. of Generated PE Reads No. of Mapped PE Reads No. of Uniquely Mapped PE Reads% of Uniquely Mapped PE Reads
1/RC111,679,44011,506,60611,203,82095.93
2/RC210,814,46810,674,01310,406,09596.22
3/RC39,853,0629,675,8639,402,09695.42
4/RC410,773,91210,497,72710,223,07094.89
5/RC510,407,98410,276,90110,003,08396.11
6/RC69,827,8479,684,3049,411,73195.77
7/RC710,068,9259,867,6429,628,73895.63
8/RC89,866,3009,636,1039,384,36195.12
9/RC99,846,3989,631,7699,327,98794.74
10/RC109,854,6649,722,2899,476,47096.16
11/RC119,880,0309,722,3999,479,34495.94
12/RC1212,482,58612,229,11511,918,10195.48
Mean10,446,30110,260,3949,988,74195.62
Table 4. Total number of generated PE reads aligned to the reference genome, and total number and percentage of uniquely mapped PE reads in samples collected from G2 animals (creep feeding) at weaning.
Table 4. Total number of generated PE reads aligned to the reference genome, and total number and percentage of uniquely mapped PE reads in samples collected from G2 animals (creep feeding) at weaning.
Animal/SampleNo. of Generated PE ReadsNo. of Mapped PE ReadsNo. of Uniquely Mapped PE Reads% of Uniquely Mapped PE Reads
13/RC139,479,6459,250,1469,033,77495.30
14/RC1410,640,50810,341,88910,068,52394.62
15/RC159,553,3139,373,2669,141,01795.68
16/RC169,333,5929,207,4128,986,54396.28
17/RC179,129,5209,022,7138,783,63696.21
18/RC1810,213,50210,046,5199,797,82995.93
19/RC1910,152,41510,000,4839,757,87496.11
20/RC2010,134,1399,981,2769,736,55596.08
21/RC219,651,4819,260,1509,018,95193.45
22/RC228,606,9668,434,8958,233,20195.66
23/RC239,538,1579,359,8969,132,29795.74
24/RC249,435,2529,215,5918,995,31795.34
Mean9,655,7089,457,8539,223,79395.53
Table 5. Top 30 genes identified as differentially expressed between G1 (no creep feeding) × G2 (creep feeding) at weaning.
Table 5. Top 30 genes identified as differentially expressed between G1 (no creep feeding) × G2 (creep feeding) at weaning.
Gene ID EnsemblGene SymbolRegulated 1log2 FC 2FDR 3
ENSBTAG00000017280C3Down−2.882.35 × 10−19
ENSBTAG00000046307CEBPDDown−2.562.31 × 10−16
ENSBTAG00000021672RGS1Down−2.512.70 × 10−16
ENSBTAG00000048501ST8SIA2Up2.378.32 × 10−15
ENSBTAG00000011121CLCN4Up1.021.52 × 10−14
ENSBTAG00000014069PDK4Down−3.341.52 × 10−14
ENSBTAG00000004248MLYCDDown−1.892.63 × 10−14
ENSBTAG00000013242MYMKUp1.472.67 × 10−14
ENSBTAG00000002834CCDC69Up1.136.93 × 10−14
ENSBTAG00000022989FAM174BUp1.128.62 × 10−14
ENSBTAG00000013631GLULDown−1.541.08 × 10−13
ENSBTAG00000017956KCTD15Up1.021.27 × 10−13
ENSBTAG00000013860GADD45ADown−1.961.75 × 10−13
ENSBTAG00000002783PCYOX1Up0.836.15 × 10−13
ENSBTAG00000007890SEC14L5Up2.076.85 × 10−13
ENSBTAG00000021999CPT1ADown−1.707.95 × 10−13
ENSBTAG00000004118ALAS1Up0.741.13 × 10−12
ENSBTAG00000046548ST6GALNAC4Up0.868.42 × 10−12
ENSBTAG00000007578SHTN1Down−1.692.05 × 10−11
ENSBTAG00000050158----------Up0.862.05 × 10−11
ENSBTAG00000011909ACVR1Up0.623.64 × 10−11
ENSBTAG00000012314LDLRUp2.714.10 × 10−11
ENSBTAG00000015942DNAJA4Up1.334.11 × 10−11
ENSBTAG00000014265SREBF2Up0.904.15 × 10−11
ENSBTAG00000050852CXCL9Down−2.344.23 × 10−11
ENSBTAG00000000163DDIT4Down−1.165.15 × 10−11
ENSBTAG00000011437----------Down−1.546.69 × 10−11
ENSBTAG00000016819FABP3Up1.417.36 × 10−11
ENSBTAG00000048728----------Up1.367.55 × 10−11
ENSBTAG00000017280PMEPA1Up0.901.19 × 10−10
1 Downregulated or upregulated in G2; 2 log2 fold change; 3 adjusted p-value (FDR) obtained by edgeR used to rank the differentially expressed genes presented in the table.
Table 6. Enriched metabolic pathways (p < 0.05) highlighted in over-representation analysis of up- and down-regulated differentially expressed genes identified in G2 at weaning.
Table 6. Enriched metabolic pathways (p < 0.05) highlighted in over-representation analysis of up- and down-regulated differentially expressed genes identified in G2 at weaning.
Term (KEGG)Up/DownNo. of Genesp-ValueGenes
PPAR signaling pathwayUp9<0.001FADS2, FABP3, SLC27A1, SCD, SCD5, AQP7, DBI, PPARA, RXRG
Steroid biosynthesisUp5<0.001SQLE, EBP, CYP51A1, DHCR24, LSS
Biosynthesis of unsaturated fatty acidsUp40.029FADS2, SCD, SCD5, FADS1
Apelin signaling pathwayUp80.039PRKAB2, SMAD3, PRKAA2, CCND1, MYL2, MYL3, APLNR, CALM3
Fatty acid metabolismUp50.041FADS2, SCD, FASN, SCD5, FADS1
AMPK signaling pathwayDown80.023PFKFB4, CPT1A, PFKFB3, EIF4EBP1, MLYCD, CPT1B, FBP1, FOXO1
Glucagon signaling pathwayDown70.031CPT1A, GCGR, SIK1, CPT1B, FBP1, PLCB2, FOXO1
PPAR signaling pathwayDown60.044CPT1A, APOA1, ME3, ANGPTL4, CPT1B, PLIN5
Table 7. Enriched biological processes (p < 0.05) highlighted in over-representation analysis of up- and down-regulated differentially expressed genes identified in G2 at weaning.
Table 7. Enriched biological processes (p < 0.05) highlighted in over-representation analysis of up- and down-regulated differentially expressed genes identified in G2 at weaning.
Term (GO_BP)Up/DownNo. of Genesp-ValueGenes
Cholesterol biosynthetic processUp50.002EBP, INSIG1, CYP51A1, DHCR24, LSS
Unsaturated fatty acid biosynthetic processUp40.003FADS2, SCD, SCD5, FADS1
Sterol biosynthetic processUp30.021SQLE, EBP, INSIG1
Cellular response to insulin stimulusUp50.032GOT1, INSIG1, INHBB, LPIN2, HDAC9
Positive regulation of lipid storageDown4<0.001C3, IKBKE, FAM71F2, PLIN5
Fatty acid metabolic processDown60.005C3, CPT1A, ACOT7, UCP3, HACL1, CPT1B
Negative regulation of glycolytic processDown30.022DDIT4, NUPR1, FBP1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ramírez-Zamudio, G.D.; Ganga, M.J.G.; Pereira, G.L.; Nociti, R.P.; Chiaratti, M.R.; Cooke, R.F.; Chardulo, L.A.L.; Baldassini, W.A.; Machado-Neto, O.R.; Curi, R.A. Effect of Cow-Calf Supplementation on Gene Expression, Processes, and Pathways Related to Adipogenesis and Lipogenesis in Longissimus thoracis Muscle of F1 Angus × Nellore Cattle at Weaning. Metabolites 2023, 13, 160. https://doi.org/10.3390/metabo13020160

AMA Style

Ramírez-Zamudio GD, Ganga MJG, Pereira GL, Nociti RP, Chiaratti MR, Cooke RF, Chardulo LAL, Baldassini WA, Machado-Neto OR, Curi RA. Effect of Cow-Calf Supplementation on Gene Expression, Processes, and Pathways Related to Adipogenesis and Lipogenesis in Longissimus thoracis Muscle of F1 Angus × Nellore Cattle at Weaning. Metabolites. 2023; 13(2):160. https://doi.org/10.3390/metabo13020160

Chicago/Turabian Style

Ramírez-Zamudio, Germán Darío, Maria Júlia Generoso Ganga, Guilherme Luis Pereira, Ricardo Perecin Nociti, Marcos Roberto Chiaratti, Reinaldo Fernandes Cooke, Luis Artur Loyola Chardulo, Welder Angelo Baldassini, Otávio Rodrigues Machado-Neto, and Rogério Abdallah Curi. 2023. "Effect of Cow-Calf Supplementation on Gene Expression, Processes, and Pathways Related to Adipogenesis and Lipogenesis in Longissimus thoracis Muscle of F1 Angus × Nellore Cattle at Weaning" Metabolites 13, no. 2: 160. https://doi.org/10.3390/metabo13020160

APA Style

Ramírez-Zamudio, G. D., Ganga, M. J. G., Pereira, G. L., Nociti, R. P., Chiaratti, M. R., Cooke, R. F., Chardulo, L. A. L., Baldassini, W. A., Machado-Neto, O. R., & Curi, R. A. (2023). Effect of Cow-Calf Supplementation on Gene Expression, Processes, and Pathways Related to Adipogenesis and Lipogenesis in Longissimus thoracis Muscle of F1 Angus × Nellore Cattle at Weaning. Metabolites, 13(2), 160. https://doi.org/10.3390/metabo13020160

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