Maternal nutrition alters mRNA isoform expression, usage, and splicing dynamics in skeletal muscle of beef cattle offspring

Background: Maternal nutrition during gestation plays a critical role in fetal muscle development and long-term metabolic programming; however, its persistent molecular effects on offspring skeletal muscle remain unclear. Therefore, the main objective of this study was to investigate the influence of prenatal nutrition on long-term differential gene expression (DGE), differential mRNA transcript expression (DTE), and differential transcript usage (DTU) in skeletal muscle of beef cattle.

Methods: A total of 126 pregnant Nellore cows were assigned to three dietary treatments: mineral supplementation only (NP), protein-energy supplementation during late gestation (PP), and protein-energy supplementation throughout gestation (FP). At 676 ± 28 days of age, muscle samples were collected from offspring for RNA sequencing. The DGE and DTE analyses were performed using the edgeR package, while DTU was evaluated with the IsoformSwitchAnalyzeR package. Over-representation analysis was conducted using g:Profiler.

Results: A total of 27,412 genes and 111,185 transcripts, including novel loci and isoforms were identified. Gene-level differences were modest (16 genes), whereas transcript-level analyses revealed stronger effects, with a higher number of significant expression and usage changes across conditions. The FP × NP comparison exhibited the greatest impact on gene expression, with 14 DTEs and 87 DTUs, compared with 12 and 30 in PP × NP, and 10 and three in FP × PP, respectively. Isoform switching was observed in key genes including SLC7A8, SLC25A30, SORBS3, and CDH13 genes, influencing coding potential, functional domains, and mRNA stability, with potential consequences for amino acid transport, cytoskeletal organization, and muscle regeneration. Functional enrichment analyses highlighted significant metabolic pathways related to amino acid and biotin metabolism, intracellular trafficking, and immune regulation.

Conclusion: Overall, prenatal nutrition, particularly protein-energy supplementation throughout gestation in comparison to mineral supplementation, modulates offspring muscle mainly through transcript usage and splicing, suggesting long-term adaptive mechanisms beyond gene-level regulation.

1 Introduction

Beef is a major source of high-quality protein, essential fatty acids, and trace minerals, playing a critical role in human health and nutrition (Murimi, 2022). However, improving beef quality remains a challenge for the cattle sector (da Silveira et al., 2021) and continues to be a central focus in beef cattle breeding and nutritional genomics research. Among the extrinsic factors influencing meat quality, maternal nutrition has been suggested to play an important role in muscle development and meat characteristics in the offspring, potentially exerting long-term effects (Du et al., 2015). This is particularly relevant given that skeletal muscle, the primary tissue contributing to meat yield, starts its development during the intrauterine period; thus, optimizing fetal muscle development may enhance both production efficiency and overall beef quality (Yan et al., 2013; Du et al., 2015; Du et al., 2017; Costa et al., 2021).

During fetal development, the total number of skeletal muscle fibers is determined exclusively within the first two-thirds of gestation and is highly sensitive to maternal nutrition and other uterine environmental factors (Du et al., 2010). As skeletal muscle has lower priority in nutrient partitioning compared to vital organs such as the brain and liver, it is particularly vulnerable to maternal nutritional fluctuations (Long et al., 2009; Meyer et al., 2010). Phenotypic variations related to skeletal muscle tissue, including body morphology, body weight, growth rate, longissimus dorsi area, intramuscular fat (marbling content), and tenderness, observed postnatally originate from developmental alterations that occur during the prenatal period (Summers et al., 2011; Copping et al., 2021; Polizel et al., 2021b; Dean et al., 2024; Pombo et al., 2025). These outcomes are associated with epigenetic modifications established during gestation, which may induce long-lasting changes in gene expression throughout postnatal life and, in some cases, be inherited across generations (Shokrollahi et al., 2025).

Several recent studies have investigated the impact of prenatal nutrition on the muscle transcriptome (Diniz et al., 2021; Palmer et al., 2021; Carvalho et al., 2022; Shokrollahi et al., 2024; Keogh et al., 2025). However, most of these studies have focused primarily on gene-level expression analyses, which while informative, may overlook subtle yet functionally important differences in transcript-level dynamics (i.e., the expression levels of individual mRNA isoforms and in the relative usage of alternative splice transcripts (Muniz et al., 2021; Asselstine et al., 2022; Singh and Ahi, 2022). Changes in transcript usage (i.e., the relative proportion of mRNA isoforms expressed from a single gene) and mRNA isoform switching (i.e., a change in which isoform is predominantly expressed) can result in substantial downstream effects on protein structure, cellular metabolism, and ultimately, phenotypic outcomes (Wang et al., 2008; Pal et al., 2011; Jiang et al., 2023).

Here, we hypothesized that prenatal nutrition modulates the muscle transcriptome of beef cattle offspring at the isoform level, primarily through alternative splicing (AS) mechanisms. Therefore, the overarching objective of this study was to investigate the long-term effects of maternal nutrition during gestation on muscle transcript expression and transcript usage in both novel and known genes and transcripts. The specific objectives were to: (1) identify novel beef muscle genes and transcripts; (2) evaluate differential gene- and transcript-level expression profiles among offspring from contrasting maternal nutrition groups; (3) characterize transcript usage patterns, mRNA isoform switches, and AS events associated with maternal nutrition; and, (4) perform functional enrichment analyses to elucidate the molecular pathways and biological processes underlying the observed transcriptomic differences.

2 Materials and methods

2.1 Experimental design

This study is part of a broader research initiative to evaluate the long-term molecular effects of maternal nutrition on the development and performance of beef cattle offspring. The experimental animals were managed at the Faculty of Animal Science and Food Engineering (FZEA-USP) campus (Pirassununga, SP, Brazil). The Research Ethics Committee of FZEA-USP approved the study under protocol number 1843241117, in accordance with the ethical standards and guidelines of the National Council for the Control of Animal Experimentation.

The nutritional interventions applied to 126 pregnant Nellore cows during gestation were previously described by Schalch Junior et al. (2022). Briefly, cows were artificially inseminated (AI) using semen from four sires, and pregnancy was confirmed 30 days post-AI. The dams were then stratified into three experimental groups (n = 42 per group), balanced for age, body weight (BW), and body condition score. Each group was maintained in grazing paddocks with Urochloa brizantha cv. Marandu, which were equipped with troughs for water and feed supplementation. The cows were assigned to one of three prenatal nutritional treatments: (1) NP (Not Programmed; control) – mineral supplementation only (0.3 g/kg BW per day) throughout gestation; (2) PP (Partially Programmed) – protein–energy supplementation (3 g/kg BW per day) provided only during the third trimester; and, (3) FP (Fully Programmed) – the same protein–energy supplementation (3 g/kg BW per day) provided continuously from pregnancy confirmation until calving. All groups received mineral supplementation at 0.3 g/kg of BW per day; however, in the PP and FP groups, the minerals were incorporated into the protein–energy supplement formulation (Table 1). The quality of the grazing areas was monitored through periodic bromatological analyses to ensure comparable nutrient availability among treatments. The total digestible nutrient (TDN) content averaged 63.07%, 64.1%, and 61.43% for the NP, PP, and FP groups, respectively, while crude protein (CP) levels were 7.38%, 7.82%, and 7.40%, respectively. Further methodological details and the effects of maternal nutritional treatments on dam performance and metabolism are provided in Schalch Junior et al. (2022).

Table 1

Table 1. Composition of the different supplements offered to the Nellore cows during the gestational period.

After calving, protein–energy supplementation for the PP and FP dams was discontinued. Thereafter, all groups, including NP, received 0.3 g/kg BW of mineral supplement daily and were managed together in a single grazing system with Urochloa brizantha cv. Marandu. Standardized health and feeding protocols were followed for cow–calf management until weaning at 240 ± 28 days. Post-weaning, calves were separated by sex and reared under uniform management conditions until reaching 570 ± 28 days of age. The young bulls received two types of supplements depending on the season: an energy supplement during the dry season and a protein supplement during the wet season, while continuing to graze on the same pasture. The finishing phase began at 570 ± 28 days of age and included 63 bulls, lasting until slaughter at 676 ± 28 days of age. During this period, the animals received three sequential diets: a 15-day adaptation diet, a subsequent diet for 35 days, and a final diet for 56 days (Polizel et al., 2022a). At the end of the experimental period, all animals were slaughtered at the FZEA-USP slaughterhouse, following the regulations established by the Brazilian Ministry of Agriculture, Livestock, and Food Supply (MAPA; Brasilia, DF, Brazil). Figure 1 illustrates a schematic overview of the experimental design.

Figure 1

Figure 1. Schematic overview of the experimental design and nutritional treatments (NP, PP, and FP) applied during gestation and postnatal development in Nellore cattle.

2.2 Collection of muscle samples

Immediately after slaughter, samples of the longissimus thoracis muscle (located between the 12th and 13th ribs) were collected from all bulls for transcriptomic analyses. Tissue samples were obtained within 15 min postmortem, rapidly snap-frozen in liquid nitrogen, and stored at −80 °C until RNA extraction. From the 63 bulls, a subset of 15 animals (n = 5 per maternal nutrition group), all sired by the same bull (to minimize paternal genetic variation), was randomly selected for RNA sequencing (RNA-Seq) analyses.

2.3 RNA isolation, library preparation, and sequencing

Total RNA was extracted from 100 mg of longissimus thoracis muscle using TRIzol reagent (Thermo Fisher Scientific, Waltham, MA, United States), following the manufacturer’s instructions. RNA concentration and purity were assessed with a DS-11 spectrophotometer (DeNovix, Wilmington, DE, United States), and RNA integrity was evalauted with the Bioanalyzer 2,100 (Agilent Technologies, Santa Clara, CA, United States). Only samples with a RNA Integrity Number (RIN) > 7.0 were retained for subsequent analyses.

For library construction, 0.1–1 µg of high-quality RNA was used following the TruSeq Stranded mRNA Sample Preparation protocol (Illumina, San Diego, CA, United States). Library quantification was performed via quantitative PCR (qPCR) using the KAPA Library Quantification kit (Roche, Basel, Switzerland), and average fragment size was determined with the Bioanalyzer. Sequencing was performed on a single flow-cell lane using the TruSeq PE Cluster kit v3-cBot-HS and a paired-end sequencing configuration on the Illumina HiSeq 2,500 platform. All sequencing procedures were conducted by NGS Soluções Genômicas (Piracicaba, São Paulo, Brazil).

2.4 Mapping, assembly, classification, and quantification of novel and known transcripts

Quality control of the raw RNA-Seq reads was performed using FastQC (v0.12.1), and adapter sequences as well as low-complexity reads were removed with SeqyClean (v1.10.09; Zhbannikov et al., 2017). The cleaned reads were then aligned to the Bos taurus reference genome (ARS-UCD1.3) using the STAR aligner (version 2.7.11b; Dobin et al., 2013) in two-pass mode (Veeneman et al., 2016) based on the Ensembl genome annotation (release 113) database. The resulting coordinate-sorted BAM files (“Aligned.sortedByCoord.out.bam”) were assembled into transcripts using StringTie (v.2.1.7; Pertea et al., 2015) with default parameters. Subsequently, the transcriptomes from all the samples were merged using the StringTie “merge” function to produce a unified transcriptome assembly. By default, StringTie assigns the prefix “MSTRG” to novel genes and transcripts.

To classify known and novel transcripts, the unified GTF file was compared to the Bos taurus reference annotation (ARS-UCD1.3.113) using GffCompare (v.0.12.10; Pertea and Pertea, 2020). Predicted transcripts were categorized according to GffCompare class codes as follows: known isoforms of a reference gene (class code “=”), novel isoforms of a reference gene (class codes “c”, “k”, “j”, “m”, “n”, or “o”), novel loci (class codes “i”, “u”, “y”, or “x”), and, potential artifacts (class codes “e”, “s”, or “p”) (Halstead et al., 2021).

The GffRead tool (v.0.12.7; Pertea and Pertea, 2020) was used to generate a transcriptome FASTA file from the StringTie-assembled GTF file and the Bos taurus reference genome FASTA. This custom transcriptome assembly was then indexed and used for transcript quantification with Salmon (v1.10.1; Patro et al., 2017) operating in quasi-mapping mode to estimate transcript-level abundance.

2.5 Differential gene expression (DGE) and differential transcript expression (DTE) analyses

Transcript-level counts from Salmon were imported into R (v4.4.1) using the tximport package (v.1.32.0; Soneson et al., 2016). A DGEList object was then created using edgeR (v.4.2.2; Robinson et al., 2010). Lowly expressed transcripts were filtered out using the “filterByExpr” function and 39,280 transcripts were kept for subsequent analyses. Data normalization was performed using the Trimmed Mean of M-values (TMM) method based on library sizes. Differential expression analyses were conducted in edgeR using the exact test based on the negative binomial distribution. The analysis included the prenatal nutritional group as the only factor of interest. Statistical significance was assessed using the “exactTest” function, and p-values were adjusted for multiple testing using the false discovery rate (FDR) approach.

The normalized counts were used for differential expression (DE) analysis, whereas log₂-transformed counts per million (CPM) values were used for performing principal component analysis (PCA) using the “prcomp” R function. For gene-level analyses, the same computational workflow was employed, but gene expression matrices were obtained using the “extractGeneExpression” function from the IsoformSwitchAnalyzeR package (v.2.2.0; Vitting-Seerup et al., 2019). After filtering, 15,832 genes were retained for downstream analyses. Transcripts and genes were considered significantly differentially expressed when the FDR-adjusted p-value was lower than 0.05 and the absolute log₂ fold change (|log₂FC|) exceeded 1.

2.6 Differential transcript usage (DTU)

The pipeline of IsoformSwitchAnalyzeR (Vitting-Seerup et al., 2019) was used to investigate differential transcript usage (DTU). Only genes with two or more isoforms (“removeSingleIsoformGenes” = TRUE), an average expression >10 transcripts per million (TPM; “geneExpressionCutoff” = 10) across experimental groups, and at least one isoform expressed >3 TPM (“isoformExpressionCutoff” = 3) in at least one group were retained. After filtering, 7,973 isoforms remained for downstream analyses.

Isoform switching in each contrast (FP × NP, FP × PP, and NP × PP) was assessed using the DEXSeq method (Anders et al., 2012), which implements a generalized linear model (GLM) framework based on the negative binomial distribution to test for differences in relative transcript (or exon) usage between conditions. In this approach, count data from transcript regions are modeled so that changes in relative usage across groups can be detected while accounting for biological variability. DEXSeq fits models with and without condition–region (exon) interactions and performs likelihood ratio tests to identify significant differential usage, with p-values adjusted for multiple testing using the false discovery rate (FDR) method. The exon- or region-level differential usage results are subsequently integrated within the IsoformSwitchAnalyzeR framework to infer transcript-level (isoform) usage patterns. The difference in isoform fraction (dIF) was calculated as IF₂ – IF₁, where IF₁ and IF₂ represent the isoform fractions of the first and second groups in each contrast, respectively (e.g., FP in FP × NP or FP × PP). The dIF value serves as an estimate of effect size, analogous to fold change in differential expression analysis. Isoform switches were considered as significant when |dIF| > 0.1 and FDR <0.05.

Further analyses of the switched isoforms were performed to evaluate their potential consequences on gene function, including intron retention (IR), changes in the amino acid sequences of open reading frames (ORFs), susceptibility to nonsense-mediated decay (NMD), modifications in protein domain structures, and variations in coding potential. A comprehensive examination of notable isoform switches was conducted using multiple bioinformatics tools: IUPred2A web server (Mészáros et al., 2018) for predicting intrinsically disordered regions; SignalP (v.5.0; Teufel et al., 2022); for signal peptide prediction; CPC2 (v.1.0.1; Kang et al., 2017) for coding potential evaluation; Pfam (v.37.0; Finn et al., 2014); for identifying conserved protein domains; DeepTMHMM (v.1.0.44; Hallgren et al., 2022); for protein topology prediction; and DeepLoc2 (v.2.0; Thumuluri et al., 2022); for subcellular location inference. In addition, AS events associated with isoform switching were investigated, including alternative 3′ and 5′ splice sites (A3, A5), alternative transcription start and termination sites (ATSS, ATTS), exon skipping (ES), and intron retention (IR). These events can reveal potential modifications in gene function mediated by altered exon junctions and transcript structural variation.

2.7 Over-representation analysis (ORA) and overlapping transcripts and processes

As the primary objective of this study was to evaluate the effects of maternal nutrition at the transcript level, rather than through conventional gene-level analyses, the over-representation analysis (ORA) was performed using transcript-based data. Genes corresponding to significant transcripts identified in differential transcript expression (DTE) or DTU analyses were used as input for enrichment testing. ORA was conducted separately for each contrast (FP × PP, PP × NP, and FP × NP) and for both analysis type (DTE and DTU). Functional enrichment was carried out using Gene Ontology (GO) analysis including the three main GO categories: Biological Process (BP), Molecular Function (MF), and Cellular Component (CC). Furthermore, metabolic pathways analyses were performed using the Reactome (REAC) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases via the g:Profiler platform (Kolberg et al., 2023). Enriched processes were considered significant when FDR was lower than 0.05. The overlap of significant transcripts between DTE and DTU analyses, as well as the intersection of enriched biological processes, was visualized using the R packages “eulerr” (v7.0.2) and “ggvenn” (v0.1.10). A schematic overview of the main bioinformatics workflow is presented in Figure 2.

Figure 2

Figure 2. Bioinformatics workflow to identify and characterize genes, transcripts and processes relevant to prenatal nutrition in offspring muscle.

3 Results

3.1 Classification of muscle genes and transcripts

A total of 27,412 genes were identified in the assembled muscle transcriptome, including 26,080 known genes and 1,332 novel genes (Figure 3A). At the transcript level, 111,185 transcripts were detected, of which 75,284 transcripts were classified as known and 35,901 as novel (Figure 3B). Among the novel transcripts, 32,705 originated from reference genes, 3,196 from novel genes, and 177 were categorized as potential artifacts (Figure 3C). Comprehensive details for each gene and transcript are provided in the Supplementary Material: Additional File 1 (known transcripts), Additional File 2 (novel loci), Additional File 3 (novel transcripts from reference genes), and Additional File 4 (potential artifacts).

Figure 3

Figure 3. Classification of the transcripts and genes identified. (A) Proportion of known and novel genes. (B) Proportion of known and novel transcripts. (C) Categories of novel transcripts identified.

3.2 Principal component analysis (PCA) and differential expression of genes and transcripts (DGE and DTE)

The PCA of both gene-level (Figure 4A; PC1 = 14% and PC2 = 12.4%) and transcript-level (Figure 4B; PC1 = 10.7% and PC2 = 9.3%) expression data showed similar distribution patterns, with considerable overlap among the experimental groups. This overlap indicates the lack of a clear global distinction between the maternal nutrition treatments at either gene or transcript level.

Figure 4

Figure 4. Principal component analysis (PCA) plot and volcano plot illustrating muscle gene and transcript profiles according to the maternal nutritional groups (NP, PP, and FP). (A) PCA of gene expression levels from different prenatal groups (NP, PP, and FP). (B) PCA of transcript expression levels from different prenatal groups (NP, PP, and FP). (C) Volcano plot of the differentially expressed genes across the group comparisons (FP × NP, FP × PP, and PP × NP). (D) Volcano plot of the differentially expressed transcripts across the group comparisons (FP × NP, FP × PP, and PP × NP).

Differential expression analysis enabled the identification of 13 DEGs between the PP and NP groups (FDR <0.05 and |log₂FC| > 1), two genes between FP × NP, and one gene between FP × PP (Figure 4C; Additional File 5). At the transcript level (Figure 4D; Additional File 6), 14 transcripts were DE in the FP × NP contrast, 12 in PP × NP, and 10 in FP × PP. Notably, among the 13 DEGs, five corresponded to novel genes (MSTRG.19105, MSTRG.2378, MSTRG.10953, MSTRG.16026, and MSTRG.13630), while at the transcript level, 10 unique novel transcripts were identified (Figure 4).

3.3 Differentially transcript usage (DTU), functional consequences and splicing events

A total of 120 transcripts exhibited DTU across the three maternal nutrition contrasts: NP × PP, FP × NP, and FP × PP (Figure 5A; Additional File 7). Among these, 30 transcripts were detected in NP × PP comparison, 87 in FP × NP, and three in FP × PP. We further assessed the functional impact of isoform switching for each significant isoform detected in the DTU analysis (Figure 6A; Additional File 8). Although some isoform switching events occurred in the FP × PP comparison, none resulted in detectable functional consequences (Figure 6A). In contrast, 21 switch consequences were detected in the NP × PP comparison, involving different domains, coding potential, ORF similarity, intron retention, and other features across 11 genes. For the FP × NP contrast, 93 switch consequences were detected, affecting 38 unique genes. Representative examples of switch consequences in key genes are illustrated in Figure 5B (SORBS3 gene), Figure 5C (SLC7A8 gene), and Figure 5D (SLC25A30 gene). These genes were selected as illustrative cases because they showed robust and biologically meaningful isoform switching events (high dIF values and significant FDR), and because they are functionally linked to metabolic, transport, and cellular adaptation processes previously reported to be sensitive to maternal nutritional status. As such, they provide representative examples of how prenatal nutrition may influence transcript usage with potential functional consequences in offspring muscle. Figure 6A provides an overview of the number of genes exhibiting isoform switching with associated functional consequences across all pairwise comparisons.

Figure 5

Figure 5. Differential transcript usage (DTU) analysis and main isoforms illustrated. (A) Volcano plot of DTU analysis across group comparisons (FP × NP, FP × PP, and NP × PP). (B) SORBS3 gene in the NP × PP contrast showing the isoforms identified and results from DTE, DGE, and DTU analyses. (C) SLC7A8 gene in the FP × NP contrast showing the isoforms identified and results from DTE, DGE, and DTU analyses. (D) SLC25A30 gene in the FP × NP contrast showing the isoforms identified and results from DTE, DGE, and DTU analyses.

Figure 6

Figure 6. Number of genes in each contrast (FP × NP, FP × PP, and NP × PP) with (A) associated switch isoform consequences and (B) alternative splicing events.

AS events associated with isoform switching are shown in Figure 6B, with detailed information for each affected gene and isoform provided in Additional File 9. In the FP × NP contrast, the most frequent AS event was ATTS in the less-used isoform. For the FP × PP comparison, only two events AS were detected: ATSS in the more-used isoform and ES in the less-used isoform. Finally, in the NP × PP contrast, the predominant AS event was ATSS in the more-used isoform.

3.4 Over-representation analyses (ORA)

In the FP × NP comparison, ORA of the genes identified through the DTE analysis revealed significant enrichment of pathways related to biotin metabolism (GO:MF:0004077 [FDR = 0.0196], GO:MF:0009374 [FDR = 0.0196], and KEGG:00780 [FDR = 0.0226]). Additional enriched pathways were associated with transport of small molecules (REAC:R-BTA-382551 [FDR = 0.0260]), and with cell adhesion and cadherin regulation (REAC:R-BTA-9759476 [FDR = 0.0260], REAC:R-BTA-9764260 [FDR = 0.0260], REAC:R-BTA-9762292 [FDR = 0.0260], and REAC:R-BTA-9759475 [FDR = 0.0260]. The endosomal/vacuolar metabolic pathway was also significantly enriched (REAC:R-BTA-1236977, FDR = 0.0260).

The DTU analyses further revealed significant enrichment of pathways related to intracellular trafficking and vesicle-mediated transport (REAC:R-BTA-1236977 [FDR = 0.00047] and REAC:R-BTA-1236974 [FDR = 0.0021]), as well as cellular structural organization (GO:CC:0043229 [FDR = 0.0024], GO:CC:0031982 [FDR = 0.0024], GO:CC:0043226 [FDR = 0.0024], GO:CC:0110165 [FDR = 0.0035], GO:CC:0005622 [FDR = 0.0038], and GO:CC:0005575 [FDR = 0.0038]). In addition, immune-related pathways were significantly enriched (REAC:R-BTA-983169 [FDR = 0.0033] and REAC:R-BTA-983170 [FDR = 0.0033]).

In the FP × PP contrast, over-represented processes identified from the DTE analysis were predominantly associated with amino acid metabolism and transmembrane transport (GO:BP:0035524 [FDR = 0.0191], GO:BP:1903801 [FDR = 0.0191], GO:BP:1904273 [FDR = 0.0191], GO:BP:1904557 [FDR = 0.0191], GO:CC:0015827 [FDR = 0.0191], and GO:CC:0015829 [FDR = 0.0191]). Additional enriched processes were associated with inositol metabolism and biosynthesis (GO:MF:0004512 [FDR = 0.0109], GO:BP:0006021 [FDR = 0.0191], GO:BP:0006020 [FDR = 0.0191], and GO:MF:0016872 [FDR = 0.0109]).

In the PP × NP contrast, genes identified through DTE analysis were predominantly enriched in pathways related to biotin metabolism (GO:MF:0004077 [FDR = 0.0050], GO:BP:0009374 [FDR = 0.0050], and KEGG:00780 [FDR = 0.0155]). Additional enriched biological processes were associated with inositol metabolism and biosynthesis (GO:MF:0004512 [FDR = 0.0050], GO:MF:0016872 [FDR = 0.0050], and GO:BP:0006021 [FDR = 0.0331]). Regarding DTU-based ORA, significant enrichment was observed for molecular functions related to adiponectin binding and small nuclear ribonucleoprotein (snRNP) binding (GO:MF:0055100 [FDR = 0.028] and GO:MF:0070990 [FDR = 0.028], respectively). A summary of all over-represented biological processes and molecular functions across the three nutritional contrasts and both analytical approaches (DTE and DTU) is presented in Figure 7. Complete details of the ORA are provided in Additional File 10.

Figure 7

Figure 7. Over-representation analysis for each contrast (FP × NP, FP × PP, and PP × NP) and analysis type (DTE and DTU), highlighting significantly enriched metabolic pathways and GO terms associated with the significant transcripts. The y-axis shows enriched terms grouped by source (GO:BP, GO:MF, GO:CC, KEGG, and REAC), and the x-axis indicates the enrichment significance, expressed as–log10(FDR).

3.5 Overlapping transcripts and processes across comparisons and analyses

Among the significant isoforms identified in both the DTE and DTU analyses, only MSTRG.1066.2, an isoform of the SLC7A8 gene, was found to overlap between them (Figure 8A). The ORA also revealed processes shared across comparisons and analysis types (Figure 8B), with detailed information provided in Table 2. In the FP × NP comparison, the DTU and DTE analyses shared 15 over-represented processes, indicating convergence in transcript-level and isoform-usage responses to maternal nutrition. In contrast, no overlap was observed between DTU and DTE in the PP × NP comparison. As described in Section 3.4, no over-represented processes were identified for the DTU analysis in the FP × PP contrast; therefore, this comparison was excluded from the overlap assessment.

Figure 8

Figure 8. Exclusive and overlapping transcripts between analyses (DTE and DTU; (A)), and exclusive and overlapping enriched metabolic pathways and GO terms from comparisons with significant processes in both analyses (FP × NP and PP × NP; (B)).

Table 2

Table 2. Overlapping enriched terms across prenatal nutrition group comparisons with significant processes in both analyses (DTE and DTU).

4 Discussion

Previous studies investigated how different maternal nutrition strategies affect differential gene expression across tissues and their functional implications for Nellore cattle offspring (Polizel et al., 2022b; Polizel et al., 2025b; Cracco et al., 2024). However, identifying AS events with biologically relevant regulatory effects remains poorly explored in cattle. To our knowledge, this is the first study in Nellore cattle to investigate the impact of maternal nutrition on differential transcript regulation and its functional consequences. Splicing can generate distinct isoforms of the same gene, which may differ in regulatory elements, coding potential, and functional domains. These alterations can modify protein structure, stability, subcellular localization, and interaction networks (Hallegger et al., 2010). This study was designed to detect long-term transcriptomic changes in AS patterns leading to DTU and DTE in response to maternal nutrition strategies during gestation, with a focus on the skeletal muscle of offspring at the end of the finishing phase. Although PCA of both gene- and transcript-level expression did not reveal a clear global separation among maternal nutrition treatments, several genes and transcripts were differentially expressed. This pattern is consistent with a fine-tuning of transcriptional regulation rather than broad transcriptomic reprogramming, which would be more likely under drastic nutritional conditions such as severe overnutrition or undernutrition.

Overall, the FP treatment exhibited the strongest transcriptomic response relative to the control (NP), as reflected by the higher number of differentially expressed (n = 14) and differentially used transcripts (n = 87). Among these, a single novel transcript (MSTRG.1066.2), an isoform of the SLC7A8 gene overlapped between the two analyses. This transcript was completely absent in the FP group (IF = 0) while representing 32.6% of the total SLC7A8 expression in the NP group. The DTE analysis confirmed this pattern, showing consistent downregulation in the same isoform in the FP group. Importantly, MSTRG.1066.2 is annotated as nonsense-mediated decay (NMD)-sensitive, suggesting that it is normally targeted for degradation through the NMD pathway (Peltz et al., 1993). This finding implies that prenatal nutrition strategies can influence both isoform production and overall mRNA stability, as changes in AS may regulate the balance between functional and degradable transcripts (Kishor et al., 2019). Consequently, the suppression of this NMD-sensitive isoform in the FP group could contribute to stabilized gene expression by reducing the abundance of transcripts prone to degradation. In line with this, the main coding transcript of the SLC7A8 gene (ENSBTAT00000109260) showed increased usage in the FP group relative to NP, supporting the hypothesis that the suppression of the nonfunctional isoform favored the upregulation of the coding isoform.

These findings highlight the importance of AS and transcript-specific regulation as additional layers of post-transcriptional control that extend beyond conventional gene-level expression, potentially impacting protein synthesis, cellular function, and metabolic efficiency. Among the genes showing notable isoform switching with predicted functional consequences, SLC7A8 and SLC25A30 were selected for further discussion due to their roles in nutrient transport and energy metabolism. The SLC7A8 gene belongs to the solute carrier (SLC) superfamily, which includes more than 458 transport proteins responsible for mediating the transmembrane movement of diverse substrates, including amino acids, ions, and metabolites (Pizzagalli et al., 2021). Specifically, SLC7A8 encodes a neutral amino acid transporter that facilitates the uptake of leucine (Leu), valine (Val), glycine (Gly), alanine (Ala), serine (Ser), glutamate (Glu), and cysteine (Cys) (Meixner et al., 2020). Experimental evidence in mice demonstrates that the deletion of SLC7A8 confers protection against diet-induced obesity (DIO) and improves glucose metabolism, thereby reducing lipid accumulation (Pitere et al., 2022). In parallel, SLC25A30, a mitochondrial solute carrier, has been reported to be associated with fatty acid metabolism, adipogenesis, and lipid metabolic pathways (Gunawan et al., 2021). In the context of beef cattle, the increased usage of the coding isoforms of SLC7A8 (ENSBTAT00000109260) and SLC25A30 (ENSBTAT00000090651) observed in the FP group reflects transcript-level modulation of genes with established roles in amino acid transport and mitochondrial metabolism. SLC7A8 encodes a neutral amino acid transporter involved in the uptake of branched-chain and other neutral amino acids that contribute to metabolic homeostasis and nutrient availability (Fotiadis et al., 2013; Pizzagalli et al., 2021), whereas SLC25A30 belongs to the mitochondrial solute carrier family, whose members are central to lipid and energy metabolism (Palmieri, 2013). Although the present study does not include protein-level or functional validation, the observed transcript usage patterns are consistent with differences in carcass subcutaneous fat previously reported for the same experimental population (Fernandes et al., 2023), supporting the biological relevance of these transcriptomic changes while warranting further functional investigation.

Several overlapping biological processes found in the FP × NP comparison were associated with amino acid transport and other molecular transport activities, functions closely related to SLC family genes (Hediger et al., 2013). Similar observations were reported by Palmer et al. (2021), who demonstrated that prenatal nutrition influences amino acid transport pathways in the beef muscle transcriptome. Such pathways likely contribute to muscle development by regulating amino acid availability for protein synthesis, modulating energy metabolism, and ultimately affecting growth efficiency and meat quality in beef cattle.

In the PP × NP comparison, no overlapping transcripts or processes were identified between the DTE and DTU analyses; however, several functional implications were noteworthy. For instance, the novel transcript MSTRG.17793.12, originating from the SORBS3 gene was affected by maternal supplementation during the final third of gestation and exhibited higher usage in the PP group compared with NP. Although this transcript is protein-coding, it differs from the downregulated isoform (MSTRG.17793.7) in two key aspects: it encodes a longer ORF and shows a domain gain. Specifically, the upregulated isoform contains the Sorbin (SoHo) domain, whereas the downregulated variant lacks it, suggesting a domain-dependent functional divergence. The SORBS3 gene codes for the adaptor protein vinexin and has been demonstrated to be involved in signal transduction induced by growth factors as well as in the structure of the cytoskeleton (Kioka et al., 2002). Alterations in SoHo and SH3 domain containing 3 (SORBS3) were previously linked to epigenetic regulation in skeletal muscle (Day et al., 2017). Obesity increased promoter methylation and suppressed SORBS3 expression, whereas weight loss following Roux-en-Y gastric bypass reversed these changes and improved metabolic outcomes, including glucose homeostasis (Day et al., 2017). Thus, the differential isoform usage observed may reflect domain-dependent regulatory shifts, in which the absence of the SoHo domain in the downregulated isoform could diminish cytoskeletal signaling and cellular responsiveness, whereas the upregulated isoform may restore or even modulate these functions through distinct molecular interactions.

Another important gene with functional relevance in the PP × NP comparison was CDH13. Among the four isoforms associated with this gene, the transcript ENSBTAT00000079761 exhibited decreased usage in the PP group. Compared with the upregulated isoform (ENSBTAT00000106688), this transcript showed a loss of protein domains, suggesting a potential shift in functional capacity. The CDH13 gene encodes the T-cadherin protein, a member of the cadherin superfamily involved in cell–cell adhesion, actin cytoskeleton organization, tissue morphogenesis, homeostasis, and angiogenesis (Hulpiau and van Roy, 2009; Rubina et al., 2021). Many of these functions are mediated through ligand binding, as T-cadherin serves as a receptor for adiponectin and LDL (Rubina et al., 2021). In skeletal muscle, adiponectin–T-cadherin signaling enhances muscle regeneration through the specific binding of adiponectin to T-cadherin (Tanaka et al., 2019). Therefore, the observed differences in protein domains among CDH13 isoforms identified in the DTU analysis may directly affect muscle regenerative capacity and vascularization. Notably, one of the two enriched GO terms for DTU was adiponectin binding (GO:0055100), reinforcing this interpretation and suggesting potential molecular mechanisms through which late-gestation maternal supplementation may influence splicing regulation and exert long-term functional effects in offspring muscle.

Additional enriched terms associated with biotin (e.g., biotin binding, biotin metabolism, biotin [biotin carboxyl-carrier protein] ligase activity, and response to biotin) were affected by late gestation supplementation (PP × NP) and whole gestation supplementation (FP × NP). Biotin is a B vitamin (B7), which is an essential cofactor necessary for the activation of various carboxylase enzymes (Zempleni et al., 2008). These enzymes play key roles in fatty acids synthesis and catabolism, leucine metabolism, and gluconeogenesis (Zempleni and Kuroishi, 2012). One of the genes responsible for the differences observed in biotin metabolism is HLCS (ENSBTAG00000033679), which encodes holocarboxylase synthetase, the enzyme that catalyzes the attachment of biotin to carboxylases (León-Del-Rió et al., 2017). In our DTE analysis, HLCS transcript MSTRG.802.7 was upregulated in the NP group, linked to both PP and FP contrasts, suggesting that maternal supplementation influences biotin-dependent enzyme activity in offspring muscle. The activity of these enzymes may influence metabolic efficiency as well as energy and protein metabolism in animals that receive protein-energy supplementation during gestation (PP and FP) compared with the control (NP). In our previous study investigating the long-term effects of prenatal nutrition on metabolic pathways through integrated metabolomics and metagenomics analyses (Polizel et al., 2025a), biotin-related pathways were functionally enriched in rumen fluid. These pathways exhibited distinct regulation patterns across maternal nutrition groups, as revealed by Weighted Gene Co-expression Network Analysis (WGCNA), with some modules showing opposite correlations, indicating differential regulation of biotin metabolism (Polizel et al., 2025a). This suggests that systemic shifts in biotin metabolism may occur due to differences linked to maternal nutrition.

The differences observed between PP and FP were subtle, with no functional isoform switching detected by the DTU analyses. Notably, the PP group showed upregulation of the novel transcript MSTRG.1066.2, within the SLC7A8 gene, which was absent in the FP group, and may influence amino acid transport, as also indicated by ORA. Overall, however, the timing of protein-energy supplementation during gestation had only minor effects on muscle differential gene expression (DGE), DTE, DTU, and splicing dynamics.

5 Conclusions

Collectively, this study provides novel evidence that maternal nutrition during gestation may modulate the offspring muscle transcriptome primarily through AS and isoform regulation. Prenatal nutrition in beef cattle had a limited impact on gene expression level as revealed by differential gene expression analysis; however, more detailed transcript-level assessments (DTU and DTE) showed a greater impact of whole gestation protein–energy supplementation (FP) in comparison to the control group (NP). These transcriptomic shifts highlight the role of AS events and functional consequences mainly associated with amino acid transporters and members of the SLC superfamily. These findings indicate that transcript usage and isoform regulation provide an additional layer of molecular adaptation beyond differential gene expression, with potential long-term effects on muscle metabolism and growth efficiency. Although the overall number of significant events was modest, the consistent enrichment of processes related to amino acid and biotin metabolism suggests that maternal nutrition strategies can shape long-term offspring physiology through subtle but functionally meaningful transcript-level regulation. By applying more complex bioinformatics analyses, such as DTU and DTE, our study advances beyond traditional gene-level approaches and helps to address important gaps in the literature, offering new perspectives on how maternal nutrition can shape long-term molecular regulation in offspring muscle.

Data availability statement

The datasets presented in this study can be found in online repositories. The transcriptome dataset analyzed during the current study is available in the European Nucleotide Archive (ENA) repository (EMBL-EBI), under accession number PRJEB84398 (https://www.ebi.ac.uk/ena/browser/view/PRJEB84398).

Ethics statement

The animal study was approved by the Research Ethics Committee of the Faculty of Animal Science and Food Engineering, University of São Paulo, under protocol number 1843241117, in compliance with the National Council for the Control of Animal Experimentation recommendations. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

GP: Methodology, Writing – original draft, Formal Analysis, Investigation, Funding acquisition. ÁnC: Writing – review and editing, Investigation, Supervision, Methodology. GR-Z: Writing – review and editing. AlC: Writing – review and editing. WD: Methodology, Writing – review and editing. JF: Writing – review and editing. HF: Writing – review and editing. LB: Writing – review and editing. FC: Writing – review and editing. MS: Project administration, Funding acquisition, Investigation, Supervision, Conceptualization, Writing – review and editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This study was partially funded by the São Paulo Research Foundation (FAPESP, Sao Paulo, SP, Brazil; Process Numbers: 17/12105-2, 21/03265-1, and 23/09113-4), and the National Council for Scientific and Technological Development (CNPq, Brasilia, DF, Brazil; grant number 304968/2024-2].

Acknowledgements

The authors acknowledge the São Paulo Research Foundation (FAPESP; Sao Paulo, SP, Brazil), the Faculty of Animal Science and Food Engineering (FZEA-USP; Pirassununga, SP, Brazil) and the University of Guelph (Guelph, ON, Canada) for all their support.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The authors AnC, AlC, WZ, JF, LB, MS declared that they were an editorial board member of Frontiers at the time of submission. This had no impact on the peer review process and the final decision.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2026.1734234/full#supplementary-material

Additional File S1 | Annotation file with known transcripts identified.

Additional File S2 | Annotation file with novel genes and transcripts identified.

Additional File S3 | Annotation file with novel transcripts from reference genes identified.

Additional File S4 | Annotation file with potential artifacts found.

References

Anders, S., Reyes, A., and Huber, W. (2012). Detecting differential usage of exons from RNA-seq data. Nat. Preced. 22, 1. doi:10.1038/npre.2012.6837.2

PubMed Abstract | CrossRef Full Text | Google Scholar

Asselstine, V., Medrano, J. F., and Cánovas, A. (2022). Identification of novel alternative splicing associated with mastitis disease in Holstein dairy cows using large gap read mapping. BMC Genomics 23, 222. doi:10.1186/s12864-022-08430-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Carvalho, E. B., Costa, T. C., Sanglard, L. P., Nascimento, K. B., Meneses, J. A. M., Galvão, M. C., et al. (2022). Transcriptome profile in the skeletal muscle of cattle progeny as a function of maternal protein supplementation during mid-gestation. Livest. Sci. 263, 104995. doi:10.1016/J.LIVSCI.2022.104995

CrossRef Full Text | Google Scholar

Copping, K. J., Callaghan, M. J., Geesink, G. H., Gugusheff, J. R., McMillen, I. C., Rodgers, R. J., et al. (2021). Periconception and first trimester diet modifies appetite, hypothalamic gene expression, and carcass traits in bulls. Front. Genet. 12, 1591. doi:10.3389/fgene.2021.720242

PubMed Abstract | CrossRef Full Text | Google Scholar

Costa, T. C., Gionbelli, M. P., and Duarte, M. de S. (2021). Fetal programming in ruminant animals: understanding the skeletal muscle development to improve meat quality. Anim. Front. 11, 66–73. doi:10.1093/AF/VFAB061

PubMed Abstract | CrossRef Full Text | Google Scholar

Cracco, R. C., Alexandre, P. A., Polizel, G. H. G., Fernandes, A. C., and de Almeida Santana, M. H. (2024). Evaluation of muscle long non-coding RNA profile during rearing and finishing phase of bulls subjected to different prenatal nutritional strategies. Animals 14, 652. doi:10.3390/ANI14040652/S1

PubMed Abstract | CrossRef Full Text | Google Scholar

da Silveira, D. D., Schmidt, P. I., Campos, G. S., de Vargas, L., de Souza, F. R. P., Roso, V. M., et al. (2021). Genetic analysis of growth, visual scores, height, and carcass traits in nelore cattle. Animal Sci. J. 92, e13611. doi:10.1111/asj.13611

PubMed Abstract | CrossRef Full Text | Google Scholar

Day, S. E., Garcia, L. A., Coletta, R. L., Campbell, L. E., Benjamin, T. R., De Filippis, E. A., et al. (2017). Alterations of sorbin and SH3 domain containing 3 (SORBS3) in human skeletal muscle following Roux-en-Y gastric bypass surgery. Clin. Epigenetics 9, 1–8. doi:10.1186/S13148-017-0396-5/FIGURES/4

PubMed Abstract | CrossRef Full Text | Google Scholar

Dean, S., Gomes, M., Silva, W., Steele, M., Wood, K., Du, M., et al. (2024). Vitamin A–Enriched diet at late gestation affects intramuscular fat deposition in beef offspring. Meat Muscle Biol. 8, 1–12. doi:10.22175/MMB.17646

CrossRef Full Text | Google Scholar

Diniz, W. J. S., Bobe, G., Ward, A. K., Hall, J. A., Klopfenstein, J. J., Gultekin, Y., et al. (2021). Supranutritional maternal organic selenium supplementation during different trimesters of pregnancy affects the muscle gene transcriptome of newborn beef calves in a time-dependent manner. Genes (Basel) 12, 1884. doi:10.3390/GENES12121884/S1

PubMed Abstract | CrossRef Full Text | Google Scholar

Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., et al. (2013). STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21. doi:10.1093/bioinformatics/bts635

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, M., Tong, J., Zhao, J., Underwood, K. R., Zhu, M., Ford, S. P., et al. (2010). Fetal programming of skeletal muscle development in ruminant animals. J. Anim. Sci. 88, E51–E60. doi:10.2527/jas.2009-2311

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, M., Wang, B., Fu, X., Yang, Q., and Zhu, M. J. (2015). Fetal programming in meat production. Meat Sci. 109, 40–47. doi:10.1016/J.MEATSCI.2015.04.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, M., Ford, S. P., and Zhu, M. J. (2017). Optimizing livestock production efficiency through maternal nutritional management and fetal developmental programming. Anim. Front. 7, 5–11. doi:10.2527/af.2017-0122

CrossRef Full Text | Google Scholar

Fernandes, A. C., Beline, M., Polizel, G. H. G., Cracco, R. C., Dias, E. F. F., Furlan, É., et al. (2023). Fetal programming and its effects on meat quality of nellore bulls. Veterinary Sci. 10, 672. doi:10.3390/VETSCI10120672

CrossRef Full Text | Google Scholar

Finn, R. D., Bateman, A., Clements, J., Coggill, P., Eberhardt, R. Y., Eddy, S. R., et al. (2014). Pfam: the protein families database. Nucleic Acids Res. 42, D222–D230. doi:10.1093/NAR/GKT1223

PubMed Abstract | CrossRef Full Text | Google Scholar

Fotiadis, D., Kanai, Y., and Palacín, M. (2013). The SLC3 and SLC7 families of amino acid transporters. Mol. Asp. Med. 34, 139–158. doi:10.1016/J.MAM.2012.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Gunawan, A., Listyarini, K., Harahap, R. S., Roosita, K., Sumantri, C., Llorian, M., et al. (2021). Hepatic transcriptome analysis identifies genes, polymorphisms and pathways involved in the fatty acids metabolism in sheep. PLoS One 16, e0260514. doi:10.1371/JOURNAL.PONE.0260514

PubMed Abstract | CrossRef Full Text | Google Scholar

Hallegger, M., Llorian, M., and Smith, C. W. J. (2010). Alternative splicing: global insights. FEBS J. 277, 856–866. doi:10.1111/J.1742-4658.2009.07521.X

PubMed Abstract | CrossRef Full Text | Google Scholar

Hallgren, J., Tsirigos, K. D., Damgaard Pedersen, M., Juan, J., Armenteros, A., Marcatili, P., et al. (2022). DeepTMHMM predicts alpha and beta transmembrane proteins using deep neural networks. bioRxiv 2022. doi:10.1101/2022.04.08.487609

CrossRef Full Text | Google Scholar

Halstead, M. M., Islas-Trejo, A., Goszczynski, D. E., Medrano, J. F., Zhou, H., and Ross, P. J. (2021). Large-scale multiplexing permits full-length transcriptome annotation of 32 bovine tissues from a single nanopore flow cell. Front. Genet. 12, 664260. doi:10.3389/FGENE.2021.664260/BIBTEX

PubMed Abstract | CrossRef Full Text | Google Scholar

Hediger, M. A., Clémençon, B., Burrier, R. E., and Bruford, E. A. (2013). The ABCs of membrane transporters in health and disease (SLC series): introduction. Mol. Asp. Med. 34, 95–107. doi:10.1016/J.MAM.2012.12.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Hulpiau, P., and van Roy, F. (2009). Molecular evolution of the cadherin superfamily. Int. J. Biochem. Cell Biol. 41, 349–369. doi:10.1016/J.BIOCEL.2008.09.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, M., Zhang, S., Yin, H., Zhuo, Z., and Meng, G. (2023). A comprehensive benchmarking of differential splicing tools for RNA-seq analysis at the event level. Brief. Bioinform 24, 1–13. doi:10.1093/BIB/BBAD121

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, Y. J., Yang, D. C., Kong, L., Hou, M., Meng, Y. Q., Wei, L., et al. (2017). CPC2: a fast and accurate coding potential calculator based on sequence intrinsic features. Nucleic Acids Res. 45, W12–W16. doi:10.1093/NAR/GKX428

PubMed Abstract | CrossRef Full Text | Google Scholar

Keogh, K., Kelly, A. K., and Kenny, D. A. (2025). Effect of enhanced early life nutrition on the M. longissimus thoracis et lumborum transcriptome of crossbred beef heifer calves. Sci. Rep. 15, 1–10. doi:10.1038/S41598-025-07587-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kioka, N., Ueda, K., and Amachi, T. (2002). Vinexin, CAP/ponsin, ArgBP2: a novel adaptor protein family regulating cytoskeletal organization and signal transduction. Cell Struct. Funct. 27, 1–7. doi:10.1247/CSF.27.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kishor, A., Fritz, S. E., and Hogg, J. R. (2019). Nonsense-mediated mRNA decay: the challenge of telling right from wrong in a complex transcriptome. Wiley Interdiscip. Rev. RNA 10, e1548. doi:10.1002/WRNA.1548

PubMed Abstract | CrossRef Full Text | Google Scholar

Kolberg, L., Raudvere, U., Kuzmin, I., Adler, P., Vilo, J., and Peterson, H. (2023). g:Profiler—interoperable web service for functional enrichment analysis and gene identifier mapping (2023 update). Nucleic Acids Res. 51, W207–W212. doi:10.1093/NAR/GKAD347

PubMed Abstract | CrossRef Full Text | Google Scholar

León-Del-Rió, A., Valadez-Graham, V., and Gravel, R. A. (2017). Holocarboxylase synthetase: a moonlighting transcriptional coregulator of gene expression and a cytosolic regulator of biotin utilization. Annu. Rev. Nutr. 37, 207–223. doi:10.1146/ANNUREV-NUTR-042617-104653/CITE/REFWORKS

PubMed Abstract | CrossRef Full Text | Google Scholar

Long, N. M., Vonnahme, K. A., Hess, B. W., Nathanielsz, P. W., and Ford, S. P. (2009). Effects of early gestational undernutrition on fetal growth, organ development, and placentomal composition in the bovine. J. Anim. Sci. 87, 1950–1959. doi:10.2527/jas.2008-1672

PubMed Abstract | CrossRef Full Text | Google Scholar

Meixner, E., Goldmann, U., Sedlyarov, V., Scorzoni, S., Rebsamen, M., Girardi, E., et al. (2020). A substrate-based ontology for human solute carriers. Mol. Syst. Biol. 16, 9652. doi:10.15252/MSB.20209652

PubMed Abstract | CrossRef Full Text | Google Scholar

Mészáros, B., Erdös, G., and Dosztányi, Z. (2018). IUPred2A: context-dependent prediction of protein disorder as a function of redox state and protein binding. Nucleic Acids Res. 46, W329–W337. doi:10.1093/NAR/GKY384

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyer, A. M., Reed, J. J., Vonnahme, K. A., Soto-Navarro, S. A., Reynolds, L. P., Ford, S. P., et al. (2010). Effects of stage of gestation and nutrient restriction during early to mid-gestation on maternal and fetal visceral organ mass and indices of jejunal growth and vascularity in beef cows. J. Anim. Sci. 88, 2410–2424. doi:10.2527/jas.2009-2220

PubMed Abstract | CrossRef Full Text | Google Scholar

Muniz, M. M. M., Fonseca, L. F. S., dos Santos Silva, D. B., de Oliveira, H. R., Baldi, F., Chardulo, A. L., et al. (2021). Identification of novel mRNA isoforms associated with meat tenderness using RNA sequencing data in beef cattle. Meat Sci. 173, 108378. doi:10.1016/J.MEATSCI.2020.108378

PubMed Abstract | CrossRef Full Text | Google Scholar

Murimi, M. (2022). What is the place of beef in a healthful sustainable Homo sapiens diet? Am. J. Health Promot. 36, 898–900. doi:10.1177/08901171221088661c

PubMed Abstract | CrossRef Full Text | Google Scholar

Pal, S., Gupta, R., Kim, H., Wickramasinghe, P., Baubet, V., Showe, L. C., et al. (2011). Alternative transcription exceeds alternative splicing in generating the transcriptome diversity of cerebellar development. Genome Res. 21, 1260–1272. doi:10.1101/GR.120535.111

PubMed Abstract | CrossRef Full Text | Google Scholar

Palmer, E. A., Peñagaricano, F., Vedovatto, M., Oliveira, R. A., Field, S. L., Laporta, J., et al. (2021). Effects of maternal gestational diet, with or without methionine, on muscle transcriptome of Bos indicus-influenced beef calves following a vaccine-induced immunological challenge. PLoS One 16, e0253810. doi:10.1371/JOURNAL.PONE.0253810

PubMed Abstract | CrossRef Full Text | Google Scholar

Palmieri, F. (2013). The mitochondrial transporter family SLC25: identification, properties and physiopathology. Mol. Asp. Med. 34, 465–484. doi:10.1016/J.MAM.2012.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Patro, R., Duggal, G., Love, M. I., Irizarry, R. A., and Kingsford, C. (2017). Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 14, 417–419. doi:10.1038/NMETH.4197

PubMed Abstract | CrossRef Full Text | Google Scholar

Peltz, S. W., Brown, A. H., and Jacobson, A. (1993). mRNA destabilization triggered by premature translational termination depends on at least three cis-acting sequence elements and one trans-acting factor. Genes Dev. 7, 1737–1754. doi:10.1101/GAD.7.9.1737

PubMed Abstract | CrossRef Full Text | Google Scholar

Pertea, G., and Pertea, M. (2020). GFF utilities: Gffread and GffCompare. F1000Res 9, ISCB Comm. J-304 9. doi:10.12688/F1000RESEARCH.23297.2

PubMed Abstract | CrossRef Full Text | Google Scholar

Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T. C., Mendell, J. T., and Salzberg, S. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290–295. doi:10.1038/NBT.3122

PubMed Abstract | CrossRef Full Text | Google Scholar

Pitere, R. R., van Heerden, M. B., Pepper, M. S., and Ambele, M. A. (2022). Slc7a8 deletion is protective against diet-induced obesity and attenuates lipid accumulation in multiple organs. Biol. (Basel) 11, 311. doi:10.3390/BIOLOGY11020311/S1

PubMed Abstract | CrossRef Full Text | Google Scholar

Pizzagalli, M. D., Bensimon, A., and Superti-Furga, G. (2021). A guide to plasma membrane solute carrier proteins. FEBS J. 288, 2784–2835. doi:10.1111/febs.15531

PubMed Abstract | CrossRef Full Text | Google Scholar

Polizel, G. H. G., Fantinato-Neto, P., Rangel, R. B., Grigoletto, L., Bussiman, F. de O., Cracco, R. C., et al. (2021a). Evaluation of reproductive traits and the effect of nutrigenetics on bulls submitted to fetal programming. Livest. Sci. 247, 104487. doi:10.1016/j.livsci.2021.104487

CrossRef Full Text | Google Scholar

Polizel, G. H. G., Strefezzi, R. de F., Cracco, R. C., Fernandes, A. C., Zuca, C. B., Castellar, H. H., et al. (2021b). Effects of different maternal nutrition approaches on weight gain and on adipose and muscle tissue development of young bulls in the rearing phase. Trop. Anim. Health Prod. 53, 1–9. doi:10.1007/S11250-021-02982-Y

CrossRef Full Text | Google Scholar

Polizel, G. H. G., Cançado, F. A. C. Q., Dias, E. F. F., Fernandes, A. C., Cracco, R. C., Carmona, B. T., et al. (2022a). Effects of different prenatal nutrition strategies on the liver metabolome of bulls and its correlation with body and liver weight. Metabolites 12, 441. doi:10.3390/metabo12050441

PubMed Abstract | CrossRef Full Text | Google Scholar

Polizel, G. H. G., Cesar, A. S. M., Cracco, R. C., Fernandes, A. C., Reginato, G. M., Xavier, P. L. P., et al. (2022b). Identification of eQTLs and differential gene expression associated with fetal programming in beef cattle. J. Appl. Genet. 63, 717–728. doi:10.1007/s13353-022-00711-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Polizel, G. H. G., Diniz, W. J. S., Cesar, A. S. M., Ramírez-Zamudio, G. D., Cánovas, A., Dias, E. F. F., et al. (2025a). Impacts of prenatal nutrition on metabolic pathways in beef cattle: an integrative approach using metabolomics and metagenomics. BMC Genomics 26, 359. doi:10.1186/s12864-025-11545-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Polizel, G. H. G., Santos, M. E. P., Cesar, A. S. M., Diniz, W. J. S., Ramírez-Zamudio, G. D., Fantinato-Neto, P., et al. (2025b). Late gestation maternal nutrition has a stronger impact on offspring liver transcriptome than full-gestation supplementation in beef cattle. Vet. Sci. 12, 406. doi:10.3390/VETSCI12050406

PubMed Abstract | CrossRef Full Text | Google Scholar

Pombo, G. do V., Polizel, G. H. G., Fernandes, A. C., Furlan, É., Pratti, B. C. T., Schalch Junior, F. J., et al. (2025). The influence of maternal nutrition on the lifetime performance of nellore cattle offspring. Ruminants 5, 18. doi:10.3390/RUMINANTS5020018

CrossRef Full Text | Google Scholar

Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi:10.1093/bioinformatics/btp616

PubMed Abstract | CrossRef Full Text | Google Scholar

Rubina, K. A., Semina, E. V., Kalinina, N. I., Sysoeva, V. Y., Balatskiy, A. V., and Tkachuk, V. A. (2021). Revisiting the multiple roles of T-cadherin in health and disease. Eur. J. Cell Biol. 100, 151183. doi:10.1016/J.EJCB.2021.151183

PubMed Abstract | CrossRef Full Text | Google Scholar

Schalch Junior, F. J., Polizel, G. H. G., Cançado, F. A. C. Q., Fernandes, A. C., Mortari, I., Pires, P. R. L., et al. (2022). Prenatal supplementation in beef cattle and its effects on plasma metabolome of dams and calves. Metabolites 12, 347. doi:10.3390/metabo12040347

PubMed Abstract | CrossRef Full Text | Google Scholar

Shokrollahi, B., Park, M., Baek, Y.-C., Jin, S., Jang, G.-S., Moon, S.-J., et al. (2024). Differential gene expression in neonatal calf muscle tissues from Hanwoo cows overfed during mid to late pregnancy period. Sci. Rep. 14, 23298. doi:10.1038/s41598-024-74976-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Shokrollahi, B., Park, M., Jang, G.-S., Jin, S., Moon, S.-J., Um, K.-H., et al. (2025). Maternal overnutrition in beef cattle: effects on fetal programming, metabolic health, and postnatal outcomes. Biology 14, 645. doi:10.3390/BIOLOGY14060645

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, P., and Ahi, E. P. (2022). The importance of alternative splicing in adaptive evolution. Mol. Ecol. 31, 1928–1938. doi:10.1111/mec.16377

PubMed Abstract | CrossRef Full Text | Google Scholar

Soneson, C., Love, M. I., and Robinson, M. D. (2016). Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. F1000Res 4, 1521. doi:10.12688/F1000RESEARCH.7563.2

PubMed Abstract | CrossRef Full Text | Google Scholar

Summers, A. F., Ramsay, K. H., and Funston, R. N. (2011). CASE STUDY: the effects of maternal nutrition on steer progeny performance. Prof. Anim. Sci. 27, 251–256. doi:10.15232/S1080-7446(15)30481-2

CrossRef Full Text | Google Scholar

Tanaka, Y., Kita, S., Nishizawa, H., Fukuda, S., Fujishima, Y., Obata, Y., et al. (2019). Adiponectin promotes muscle regeneration through binding to T-cadherin. Sci. Rep. 9 (1), 1–12. doi:10.1038/s41598-018-37115-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Teufel, F., Almagro Armenteros, J. J., Johansen, A. R., Gíslason, M. H., Pihl, S. I., Tsirigos, K. D., et al. (2022). SignalP 6.0 predicts all five types of signal peptides using protein language models. Nat. Biotechnol. 40, 1023–1025. doi:10.1038/S41587-021-01156-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Thumuluri, V., Almagro Armenteros, J. J., Johansen, A. R., Nielsen, H., and Winther, O. (2022). DeepLoc 2.0: multi-label subcellular localization prediction using protein language models. Nucleic Acids Res. 50, W228–W234. doi:10.1093/NAR/GKAC278

PubMed Abstract | CrossRef Full Text | Google Scholar

Veeneman, B. A., Shukla, S., Dhanasekaran, S. M., Chinnaiyan, A. M., and Nesvizhskii, A. I. (2016). Two-pass alignment improves novel splice junction quantification. Bioinformatics 32, 43–49. doi:10.1093/BIOINFORMATICS/BTV642

PubMed Abstract | CrossRef Full Text | Google Scholar

Vitting-Seerup, K., Sandelin, A., and Berger, B. (2019). IsoformSwitchAnalyzeR: analysis of changes in genome-wide patterns of alternative splicing and its functional consequences. Bioinformatics 35, 4469–4471. doi:10.1093/BIOINFORMATICS/BTZ247

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, E. T., Sandberg, R., Luo, S., Khrebtukova, I., Zhang, L., Mayr, C., et al. (2008). Alternative isoform regulation in human tissue transcriptomes. Nature 456, 470–476. doi:10.1038/nature07509

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, X., Zhu, M.-J., Dodson, M. V., and Du, M. (2013). Developmental programming of fetal skeletal muscle and adipose tissue development. J. Genomics 1, 29–38. doi:10.7150/jgen.3930

PubMed Abstract | CrossRef Full Text | Google Scholar

Zempleni, J., and Kuroishi, T. (2012). Biotin. Adv. Nutr. 3, 213–214. doi:10.3945/AN.111.001305

PubMed Abstract | CrossRef Full Text | Google Scholar

Zempleni, J., Hassan, Y. I., and Wijeratne, S. S. K. (2008). Biotin and biotinidase deficiency. Expert Rev. Endocrinol. Metab. 3, 715–724. doi:10.1586/17446651.3.6.715

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhbannikov, I. Y., Hunter, S. S., Foster, J. A., and Settles, M. L. (2017). “Seqyclean: a pipeline for high-throughput sequence data preprocessing,” in ACM-BCB 2017 - Proceedings of the 8th ACM International Conference on Bioinformatics, Computational Biology, and Health Informatics 17, 407–416. doi:10.1145/3107411.3107446

CrossRef Full Text | Google Scholar

Glossary

A3 Alternative 3′ splice site

A5 Alternative 5′ splice site

AI Artificial insemination

Ala Alanine

AS Alternative splicing

ATSS Alternative transcription start site

ATTS Alternative transcription termination site

BP Biological processes

BW Body weight

CC Cellular component

CP Crude protein

CPM Counts per million

Cys Cysteine

DGE Differential gene expression

dIF Difference in isoform fraction

DIO Diet-induced obesity

DTE Differential transcript expression

DTU Differential transcript usage

FDR False discovery rate

FP Full programming

Glu Glutamate

Gly Glycine

GO Gene ontology

IF Isoform fraction

IR Intron retention

KEGG Kyoto Encyclopedia of Genes and Genomes

Leu Leucine

MF Molecular functions

NMD Nonsense-mediated decay

NP Not-programmed

ORA Over-representation analysis

ORFs Open reading frames

PCA Principal component analysis

PP Partial programming

REAC Reactome

RIN RNA integrity number

Ser Serine

SLC Solute Carrier superfamily

SLC16A6 Solute Carrier Family 16 Member 6

SLC25A30 Solute Carrier Family 25 Member 30

SLC7A8 Solute Carrier Family 7 Member 8

SoHo Sorbin domain

SORBS3 Sorbin and SH3 Domain Containing 3

TDN Total digestible nutrients

TMM Trimmed Mean of M-values

TPM Transcripts per million

Val Valine