Age is an important physiological factor that affects the metabolism and immune function of beef cattle. While there have been many studies using the blood transcriptome to study the effects of age on gene expression, few have been reported on beef cattle. To this end, we used the blood transcriptomes of Japanese black cattle at different ages as the study subjects and screened 1055, 345, and 1058 differential expressed genes (DEGs) in the calf vs. adult, adult vs. old, and calf vs. old comparison groups, respectively. The weighted co-expression network consisted of 1731 genes. Finally, blue, brown, and yellow age-specific modules were obtained, in which genes were enriched in signaling pathways related to growth and development and immune metabolic dysfunction, respectively. Protein-protein interaction (PPI) analysis showed gene interactions in each specific module, and 20 of the highest connectivity genes were chosen as potential hub genes. Finally, we identified 495, 244, and 1007 genes by exon-wide selection signature (EWSS) analysis of different comparison groups. Combining the results of hub genes, we found that VWF, PARVB, PRKCA, and TGFB1I1 could be used as candidate genes for growth and development stages of beef cattle. CORO2B and SDK1 could be used as candidate marker genes associated with aging. In conclusion, by comparing the blood transcriptome of calves, adult cattle, and old cattle, the candidate genes related to immunity and metabolism affected by age were identified, and the gene co-expression network of different age stages was constructed. It provides a data basis for exploring the growth, development, and aging of beef cattle.

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.

Maternal parity is an important physiological factor influencing beef cow reproductive performance. However, there are few studies on the influence of different calving periods on early growth and postpartum diseases. Here, we conducted blood transcriptomic analysis on cows of different parities for gene discovery. We used Short Time Series Expression Miner (STEM) analysis to determine gene expression levels in cows of various parities and divided multiple parities into three main periods (nulliparous, primiparous, and multiparous) for subsequent analysis. Furthermore, the top 15,000 genes with the lowest median absolute deviation (MAD) were used to build a co-expression network using weighted correlation network analysis (WGCNA), and six independent modules were identified. Combing with Exon Wide Selection Signature (EWSS) and protein-protein interaction (PPI) analysis revealed that TPCN2, KIF22, MICAL3, RUNX2, PDE4A, TESK2, GPM6A, POLR1A, and KLHL6 involved in early growth and postpartum diseases. The GO and KEGG enrichment showed that the Parathyroid hormone synthesis, secretion, and action pathway and stem cell differentiation function-related pathways were enriched. Collectively, our study revealed candidate genes and gene networks regulating the early growth and postpartum diseases and provided new insights into the potential mechanism of reproduction advantages of different parity selection.

Genomic selection (GS) is a technique that integrates genomic data, pedigree information, and individual phenotypes to enhance genetic improvements of economically important traits in livestock. While it has shown significant effects in dairy cattle, its efficacy in beef cattle is lower due to breed diversity and differences in reproductive structures. Therefore, this study evaluated the impact of heritability levels, marker densities, and assessment methods (such as pedigree-based BLUP, genomic BLUP, and weighted genomic BLUP) on genomic prediction accuracy across multiple beef cattle breeds through simulations. Three beef cattle populations were simulated with heritability levels set at 0.3, 0.5, and 0.7 and marker densities set at 50 k and 770 k. The results showed that the predictive accuracy of PBLUP and GBLUP increased with higher heritability and larger reference populations. Increasing the marker density also improved the accuracy of genomic predictions; even a low marker density (50 k SNP) can significantly enhance the accuracy of genetic evaluation, although the size of the reference population needs to be optimized according to population structure, heritability, and the genetic architecture of the trait. Overall, integrating pedigree, genomic, and weighted SNP information can significantly improve the precision of GEBV prediction and reduce bias. In particular, the wGBLUP method demonstrated an improvement in the prediction accuracy of low-heritability traits in small but high-density marker populations.

The Bactrian camel exemplifies mammalian adaptation to deserts, but the spatial organization of its innate immune system remains uncharacterized. This study integrated transcriptomes from 110 samples across 11 major tissues and organs to resolve tissue-specific gene expression and innate immune modularity. Through differential expression analysis, Tau specificity index (τ > 0.8), and machine learning validation (Random Forest F1-score = 0.86 ± 0.11), we identified 4242 high-confidence tissue-specific genes (e.g., LIPE/PLIN1 in adipose). Weighted gene co-expression network analysis (WGCNA) of 1522 innate immune genes revealed 11 co-expression modules, with six exhibiting significant tissue associations (FDR < 0.01): liver-specific (r = 0.96), spleen-adipose-enriched (r = 0.88), muscle-associated (r = 0.82), and blood-specific (r = 0.80) modules. These networks demonstrated multifunctional coordination of immune pathways—including Pattern Recognition, Cytokine Signaling, and Phagocytosis—rather than isolated functions. Our results establish that camel innate immunity is organized into spatially modular networks tailored to tissue microenvironments, providing the first systems-level framework for understanding immune resilience in desert-adapted mammals and may inform strategies for enhancing livestock resilience in arid regions.

The Bactrian camel exemplifies mammalian adaptation to deserts, but the spatial organization of its innate immune system remains uncharacterized. This study integrated transcriptomes from 110 samples across 11 tissues to resolve tissue-specific gene expression and innate immune modularity. Through differential expression analysis, Tau specificity index (τ > 0.8), and machine learning validation (Random Forest F1-score = 0.86 ± 0.11), we identified 4,242 high-confidence tissue-specific genes (e.g., LIPE/PLIN1 in adipose). Weighted gene co-expression network analysis (WGCNA) of 1,522 innate immune genes revealed 11 co-expression modules, with six exhibiting significant tissue associations (FDR < 0.01): liver-specific (r = 0.96), spleen-adipose-enriched (r = 0.88), muscle-associated (r = 0.82), and blood-specific (r = 0.80) modules. These networks demonstrated multifunctional coordination of immune pathways—including Pattern Recognition, Cytokine Sig-naling, and Phagocytosis—rather than isolated functions. Our results establish that camel innate immunity is organized into spatially modular networks tailored to tissue microenvironments, providing the first systems-level framework for understanding immune resil-ience in desert-adapted mammals and may inform strategies for enhancing livestock re-silience in arid regions.