The growth and development of muscle tissue play a pivotal role in the economic value and quality of meat in agricultural animals, garnering close attention from breeders and researchers. The quality and palatability of muscle tissue directly determine the market competitiveness of meat products and the satisfaction of consumers. Therefore, a profound understanding and management of muscle growth is essential for enhancing the overall economic efficiency and product quality of the meat industry. Despite this, systematic research on muscle development-related genes across different species still needs to be improved. This study addresses this gap through extensive cross-species muscle transcriptome analysis, combined with interpretable machine learning models. Utilizing a comprehensive dataset of 275 publicly available transcriptomes derived from porcine, bovine, and ovine muscle tissues, encompassing samples from ten distinct muscle types such as the semimembranosus and longissimus dorsi, this study analyzes 113 porcine (n = 113), 94 bovine (n = 94), and 68 ovine (n = 68) specimens. We employed nine machine learning models, such as Support Vector Classifier (SVC) and Support Vector Machine (SVM). Applying the SHapley Additive exPlanations (SHAP) method, we analyzed the muscle transcriptome data of cattle, pigs, and sheep. The optimal model, adaptive boosting (AdaBoost), identified key genes potentially influencing muscle growth and development across the three species, termed SHAP genes. Among these, 41 genes (including NANOG, ADAMTS8, LHX3, and TLR9) were consistently expressed in all three species, designated as homologous genes. Specific candidate genes for cattle included SLC47A1, IGSF1, IRF4, EIF3F, CGAS, ZSWIM9, RROB1, and ABHD18; for pigs, DRP2 and COL12A1; and for sheep, only COL10A1. Through the analysis of SHAP genes utilizing Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, relevant pathways such as ether lipid metabolism, cortisol synthesis and secretion, and calcium signaling pathways have been identified, revealing their pivotal roles in muscle growth and development.

The average daily gain (ADG) is a critical index for evaluating growth rates in cattle and is closely linked to the economic benefits of the cattle industry. Heredity is one of the factors affecting the daily gain of cattle. However, the molecular mechanisms regulating ADG remain incompletely understood. This study aimed to systematically unravel the molecular mechanisms underlying the divergence in ADG between high average daily gain (HADG) and low average daily gain (LADG) Angus cattle through integrated multi-omics analyses (microbiome, metabolome, and transcriptome), hypothesizing that the gut microbiota–host gene–metabolism axis is a key regulatory network driving ADG divergence. Thirty Angus cattle were classified according to their HADG and LADG. Fecal and serum samples were collected for 16S, fecal metabolome, and blood transcriptome analysis. The results showed that compared with the LADG group, the abundance of Firmicutes increased in the HADG group, while the abundance of Bacteroidetes and Proteobacteria decreased. Metabolomics and transcriptomic analysis revealed that KEGG pathways associated with differentially expressed genes (DEGs) and differentially abundant metabolites (DAMs) were enriched in bile acid metabolism. Spearman correlation analysis showed that Oscillospira was positively correlated with ZBTB20 and negatively correlated with RADIL. ZBTB20 was negatively correlated with dgA-11_gut_group. This study analyzed the regulatory mechanism of average daily gain of beef cattle from genetic, metabolic, and microbial levels, providing a theoretical basis for analyzing the mechanism of differential daily gain of beef cattle, and has important significance for improving the production performance of beef cattle. The multi-omics network provides biomarker foundations for machine learning-based ADG prediction models, offering potential applications in precision breeding. While these biomarkers show promise for precision breeding, their causal roles require further validation. The conclusions are derived from a single breed (Angus) and gender (castrated males). Future studies should include females and diverse breeds to assess generalizability.

The skeletal muscle of cattle is the main component of their muscular system, responsible for supporting and movement functions. However, there are still many unknown areas regarding the ranking of the importance of different types of cell populations within it. This study conducted in-depth research and made a series of significant findings. First, we trained 15 bovine skeletal muscle models and selected the best-performing model as the initial model. Based on the SHAP (Shapley Additive exPlanations) analysis of this initial model, we obtained the SHAP values of 476 important genes. Using the contributions of these 476 genes, we reconstructed a 476-gene SHAP value matrix, and relying solely on the interactions among these 476 genes, successfully mapped the single-cell atlas of bovine skeletal muscle. After retraining the model and further interpretation, we found that Myofiber cells are the most representative cell type in bovine skeletal muscle, followed by neutrophils. By determining the key genes of each cell type through SHAP values, we conducted analyses on the correlations among key genes and between cells for Myofiber cells, revealing the critical role these genes play in muscle growth and development. Further, by using protein language models, we performed cross-species comparisons between cattle and pigs, deepening our understanding of Myofiber cells as key cells in skeletal muscle, and exploring the common regulatory mechanisms of muscle development across species.

Muscle is a crucial component of cattle, playing a vital role in determining the final quality of beef. This study aimed to identify candidate genes associated with muscle growth and lipid metabolism in beef and dairy cattle by utilizing the public database of the National Center for Biotechnology Information (NCBI) to download bovine muscle transcriptome data. Through differential expression analysis, weighted gene co-expression network analysis (WGCNA), and SHapley Additive exPlanation (SHAP) explains machine learning models, we integrated and screened for relevant genes. The results showed a total of 2588 differentially expressed genes (DEGs), with 933 upregulated and 1655 downregulated in beef cattle compared to dairy cattle. In the WGCNA, the purple, black, green, red, brown, and blue modules were identified as significant modules. Based on the results of five different machine learning models, the Adaptive Boosting (AdaBoost) model demonstrated superior classification performance (accuracy = 0.84) compared to the other four models and was therefore selected as the optimal model. SHAP analysis was then employed to interpret the results, yielding the top 500 SHAP genes. In combination with DEGs and WGCNA, a total of 117 genes were identified. Subsequent functional enrichment analysis of these 117 genes revealed significant enrichment in pathways such as lipoprotein metabolic process, muscle contraction, and cytoskeleton in muscle cells, followed by interaction network analysis of genes and pathways. Ultimately, the APOA1, ACTB, S1PR1, PKLR, and SLC27A6 genes were identified as potential key regulators of lipid metabolism and muscle growth in beef and dairy cattle. In summary, this study provides a feasible method for handling large-scale transcriptome data and lays a foundation for future research on meat quality and improving the economic benefits of Holstein cattle.

Background: Against the backdrop of the large-scale and intensive development of the livestock industry, enhancing the reproductive efficiency of cattle has become a crucial factor in industrial development. Holstein cows, as the most predominant dairy cattle breed globally, are characterized by high milk yield and excellent milk quality. However, their reproductive efficiency is comprehensively influenced by a variety of complex factors, and improving their reproductive performance faces numerous challenges. The ovary, as the core organ of the female reproductive system, plays a decisive role in embryonic development and pregnancy maintenance. It is not only the site where eggs are produced and developed but it also regulates the cow’s estrous cycle, ovulation process, and the establishment and maintenance of pregnancy by secreting various hormones. The normal functioning of the ovary is crucial for the smooth development of the embryo and the successful maintenance of pregnancy. Methods: Currently, traditional sequencing technologies have obvious limitations in deciphering ovarian function and reproductive regulatory mechanisms. To overcome the bottlenecks of traditional sequencing technologies, this study selected Holstein cows as the research subjects. Ovarian samples were collected from one pregnant and one non-pregnant Holstein cow, and single-nucleus transcriptome sequencing technology was used to conduct an in-depth study on the ovarian cells of Holstein cows. Results: By constructing a cell type-specific molecular atlas of the ovaries, nine different cell types were successfully identified. This study compared the proportions of ovarian cell types under different physiological states and found that the proportion of endothelial cells decreased during pregnancy, while the proportions of granulosa cells and luteal cells increased significantly. In terms of functional enrichment analysis, oocytes during both pregnancy and non-pregnancy play roles in the “cell cycle” and “homologous recombination” pathways. However, non-pregnant oocytes are also involved in the “progesterone-mediated oocyte maturation” pathway. Luteal cells during pregnancy mainly function in the “cortisol synthesis and secretion” and “ovarian steroidogenesis” pathways; non-pregnant luteal cells are mainly enriched in pathway processes such as the “AMPK signaling pathway”, “pyrimidine metabolism”, and “nucleotide metabolism”. Cell communication analysis reveals that there are 51 signaling pathways involved in the pregnant ovary, with endothelial cells, granulosa cells, and luteal cells serving as the core communication hubs. In the non-pregnant ovary, there are 48 pathways, and the interaction between endothelial cells and stromal cells is the dominant mode. Conclusions: This study provides new insights into the regulatory mechanisms of reproductive efficiency in Holstein cows. The differences in the proportions of ovarian cell types, functional pathways, and cell communication patterns under different physiological states, especially the increase in the proportions of granulosa cells and luteal cells during pregnancy and the specificity of related functional pathways, indicate that these cells play a crucial role in the reproductive process of cows. These findings also highlight the importance of ovarian cells in pathways such as “cell cycle”, “homologous recombination”, and “progesterone-mediated oocyte maturation”, as well as the cell communication mechanisms in regulating ovarian function and reproductive performance.

Objective: Intramuscular fat (IMF) is key for meat quality. Although fibro/adipogenic progenitor cells (FAPs) participate in adipogenesis, roles of other stromal cells remain unclear. This study characterized cellular heterogeneity in cattle IMF deposition and identified essential cell types and mechanisms regulating intramuscular adipogenesis.Methods: We performed single-nucleus RNA sequencing on longissimus dorsi muscle samples from three cattle breeds (Angus, Hereford, and Simmental) with divergent IMF content. Cell populations were annotated using multiple reference-based methods. We applied Scissor to identify phenotype-associated cells, CellChat to infer intercellular communication, high-dimensional weighted gene co-expression network analysis to detect co-expression modules, and single-cell flux estimation analysis to evaluate metabolic activity. Virtual gene perturbation was evaluated using scTenifoldKnk, and drug-responsive cell subpopulations were assessed with scRank. Additional in vitro validation experiments were conducted to confirm the regulatory effect of BMPR1A.Results: Ten cell types were identified. Angus (high marbling) showed elevated proportions of adipocytes, FAPs, and vascular endothelial cells (VEndoCs), while Hereford were enriched in type II fibers. VEndoCs were consistently linked to high IMF, acting as central signaling hubs that interact extensively with adipocytes and FAPs, and exhibited high glucose metabolic activity. Co-expression analysis within VEndoCs revealed 40 candidate regulators including nine BMP signaling pathway genes. Virtual knockdown identified BMPR1A as top regulator, altering expression of 13 lipid metabolism-related genes. In vitro validation confirmed that BMPR1A inhibition downregulated PPARγ mRNA, validating its positive regulatory role in intramuscular adipogenesis. Neurocytes showed highest sensitivity to resveratrol (via ALDH2), whereas FAPs were least responsive.Conclusion: Our findings indicate VEndoCs actively regulate intramuscular adipogenesis via a BMPR1A-mediated signaling pathway, supported by in vitro validation. These results provide insights into the cellular and molecular basis of marbling in cattle, highlighting BMPR1A as a potential target for genetic or nutritional strategies to improve meat quality.