Risu Na

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.

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.

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.

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.

Follicle-stimulating hormone (FSH) regulates ovarian follicle development through specific gene expression programs. Granulosa cells (GCs) are somatic cells surrounding the oocytes, secreting gonadotropins to regulate ovulation and promote follicular development. By analyzing the effects of different doses of FSH on the proliferation of GCs, we found that adding 10 ng/mL of FSH, as the optimal concentration, could promote the growth of GCs. Furthermore, we have successfully constructed the first CRISPR-Cas9 knockout library targeting the genes on chromosomes 2 and 3 and the X chromosomes of the sheep massively parallel coding gene, as well as an ovarian GCs knockout cell library. For the first time, we have exposed the knockout cell library to a concentration of 10 ng/mL FSH to explore the underlying mechanisms. Through this screening, we have identified 836 positive–negative screening genes that are responsive to FSH, thereby revealing the regulatory mechanisms and screening the functionality of candidate genes. Next, RNA-Seq of control (0 ng/mL), low (10 ng/mL), and high (100 ng/mL) doses of FSH revealed 1708 differentially expressed genes, and combined with 836 genes, we obtained 129 FSH dose-dependent genes with extremely significant differences. This enables us to delve deeper into investigating and identifying the mechanisms by which FSH regulates GCs. More generally, we have discovered new regulatory factors and identified reproductivity-associated major effectors. These findings provide novel research directions for further studies on sheep reproduction.

Objective: Over the past few decades, the development of anti-cancer drugs in China has made outstanding achievements based on the support of national policies. To assess the progress of non-small cell lung cancer (NSCLC) drugs, we conducted a statistical analysis of clinical trials of drugs targeting NSCLC in China from 2005 to 2023.Methods: We downloaded, screened and analysed the data from three official websites, the Centre for Drug Evaluation of China National Medical Products Administration website (NMPA), ClinicalTrials.gov and the Chinese Clinical Trial Registry (ChiCTR).Results: From January 1, 2005 to April 15, 2023, a total of 1,357 drug clinical trials that met the standards were included, and the number of registered drug clinical trials has been increasing year by year, reaching the maximum of 199 in 2021. Among them, the maximum of 462 items (34.05%) in phase II clinical trials, followed by 333 (24.54%) in phase III clinical trials, and 139 (10.24%) in phase IV clinical trials. In all drug clinical trials, industry sponsored trials (ISTs) have 722 items (53.21%), which are higher than investigator-initiated trials (IITs). The clinical trials of chemical drugs have a maximum of 723 items (53.28%), while biopharmaceuticals have grown rapidly in the past 10 years, with a total of 374 (27.56%), and 48.19% of the drug clinical trials of combined medication. In addition, the geographical distribution of the leading units and participating units of Chinese drug clinical trials are uneven, and economic regions such as Beijing, Shanghai, Jiangsu are obviously ahead of other regions.Conclusion: From 2005 to 2023, the clinical trials of registered drugs for the treatment of NSCLC increased rapidly. Among them, due to the development of immunotherapy, the clinical trials of biopharmaceuticals and drugs for combined medication are growing most rapidly, while the exploration of the original drugs is a little far from enough. Our research provides a direction for the future drug clinical trials of NSCLC, laying foundation for further extending the survival rate of patients with NSCLC.

The heritability of litter size in sheep is low and controlled by multiple genes, but the research on its related genes is not sufficient. Here, to explore the expression pattern of multi-tissue genes in Chinese native sheep, we selected 10 tissues of the three adult ewes with the highest estimated breeding value in the early study of the prolific Xinggao sheep population. The global gene expression analysis showed that the ovary, uterus, and hypothalamus expressed the most genes. Using the Uniform Manifold Approximation and Projection (UMAP) cluster analysis, these samples were clustered into eight clusters. The functional enrichment analysis showed that the genes expressed in the spleen, uterus, and ovary were significantly enriched in the Ataxia Telangiectasia Mutated Protein (ATM) signaling pathway, and most genes in the liver, spleen, and ovary were enriched in the immune response pathway. Moreover, we focus on the expression genes of the hypothalamic–pituitary–ovarian axis (HPO) and found that 11,016 genes were co-expressed in the three tissues, and different tissues have different functions, but the oxytocin signaling pathway was widely enriched. To further explore the differences in the expression genes (DEGs) of HPO in different sheep breeds, we downloaded the transcriptome data in the public data, and the analysis of DEGs (Xinggao sheep vs. Sunite sheep in Hypothalamus, Xinggao sheep vs. Sunite sheep in Pituitary, and Xinggao sheep vs. Suffolk sheep in Ovary) revealed the neuroactive ligand–receptor interactions. In addition, the gene subsets of the transcription factors (TFs) of DEGs were identified. The results suggest that 51 TF genes and the homeobox TF may play an important role in transcriptional variation across the HPO. Altogether, our study provided the first fundamental resource to investigate the physiological functions and regulation mechanisms in sheep. This important data contributes to improving our understanding of the reproductive biology of sheep and isolating effecting molecular markers that can be used for genetic selection in sheep.

The phenotypic characteristics of existing domestic pigs (DPs) greatly differ from those of wild boar (WB) populations thousands of years ago. After thousands of years of human domestication, WB and DP have diverged greatly in terms of genetics. Theoretically, worldwide local pigs have independent contributions from their local WBs at the beginning of Sus scrofa domestication. The investigation of the vicissitude of the heredity material between domestic populations and their wild ancestors will help in further understanding the domestication history of domestic animals. In the present study, we performed a genome-wide association scan (GWSA) and phylogeny estimation with a total of 1098 public European Illumina 60K single nucleotide polymorphism data, which included 650 local DPs and 448 WBs. The results revealed that the phylogenetic relationship of WBs corresponds to their geographical distribution and carries large divergence with DPs, and all WB breeds (e.g., HRWB, SBWB, and TIWB) presents a closely linkage with the middle WB (e.g., HRWB, and PLWB). In addition, 64 selected candidate genes (e.g., IDH2, PIP5K1B, SMARCA2, KIF5C, and TJP2) were identified from GWSA. A total of 63 known multiple biological functional pathways were annotated by 22 genes, and ubiquinone and other terpenoid-quinone biosynthesis pathways that belong to the metabolism of cofactors and vitamins were significantly enriched (p < 0.05). The most frequent (28.57%) pathways were classified under metabolism. We confirmed that the middle European WB has made an important genetic contribution to the entire European WB populations. A series of selected genes discovered from this study provides the scientific community with a deeper understanding of the heredity performance of metabolism and emotion and the real purpose behind domestication.

This study aimed to explore the genetic basis of muscle development in goats. The transcriptome dataset for differentially expressed lncRNAs (DELs) and differentially expressed genes (DEGs) of goat muscle at different developmental stages were obtained using RNA-Seq. A total of 447,806,481 and 587,559,465 clean reads in the longissimus dorsi muscle of Dazu black goats between 75d embryonic stage and 1d after birth were generated through Illumina paired-end sequencing, and their mapping rates were 89.82 and 90.99%, respectively. Moreover, 4517 DEGs and 648 DELs were identified, and 4784 lncRNA-mRNA targeting relationships were predicted. Gene function annotation results showed that 4101 DEGs were significantly enriched to 1098 GO terms, and 2014 DEGs were significantly enriched to 40 KEGG pathways, including many GO terms and pathways related to muscle development, such as cell differentiation and Wnt signaling pathway. Then, 10 DELs and 20 DEGs were randomly selected for RT-qPCR verification, and the agreement rate between the verification and RNA-Seq results was 90%, indicating the high reliability of the RNA-Seq data analysis. In conclusion, this study obtained several mRNAs and lncRNAs related to the muscle development of Dazu black goats and identified several targeted regulatory pairs of lncRNA–mRNA. This study may serve as a reference to understand the genetic basis and molecular mechanism of muscle development in goats.