Cattleyak are the hybrid offspring between cattle and yak and combine yak hardiness with cattle productivity. Much attempt has been made to examine the mechanisms of male sterility caused by spermatogenic arrest, but yet there is no research systematically and precisely elucidated testis gene expression profiling between cattleyak and yak. To explore the higher resolution comparative transcriptome map between the testes of yak and cattleyak, and further analyze the mRNA expression dynamics of spermatogenic arrest in cattleyak. We characterized the comparative transcriptome profile from the testes of yak and cattleyak using high-throughput sequencing. Then we used quantitative analysis to validate several differentially expressed genes (DEGs) in testicular tissue and spermatogenic cells. Testis transcriptome profiling identified 6477 DEGs (2919 upregulated and 3558 downregulated) between cattleyak and yak. Further analysis revealed that the marker genes and apoptosis regulatory genes for undifferentiated spermatogonia were upregulated, while the genes for differentiation maintenance were downregulated in cattleyak. A majority of DEGs associated with mitotic checkpoint, and cell cycle progression were downregulated in cattleyak during spermatogonial mitosis. Furthermore, almost all DEGs related to synaptonemal complex assembly, and meiotic progression presented no sign of expression in cattleyak. Even worse, dozens of genes involved in acrosome formation, and flagellar development were dominantly downregulated in cattleyak. DEGs indicated that spermatogenic arrest of cattleyak may originate from the differentiation stage of spermatogonial stem cells and be aggravated during spermatogonial mitosis and spermatocyte meiosis, which contributes to the scarcely presented sperms in cattleyak.

Skeletal muscle is a complex heterogeneous tissue and characterizing its cellular heterogeneity and transcriptional and epigenetic signatures are important for understanding the details of its ontogeny. In our study, we applied scRNA‐seq and scATAC‐seq to investigate the cell types, molecular features, transcriptional and epigenetic regulation, and patterns of developing bovine skeletal muscle from gestational, lactational and adult stages. Detailed molecular analyses were used to dissect cellular heterogeneity, and we deduced the differentiation trajectory of myogenic cells and uncovered their dynamic gene expression profiles. SCENIC analysis was performed to demonstrate key regulons during cell fate decisions. We explored the future expression states of these heterogeneous cells by RNA velocity analysis and found extensive networks of intercellular communication using the toolkit CellChat. Moreover, the transcriptomic and chromatin accessibility modalities were confirmed to be highly concordant, and integrative analysis of chromatin accessibility and gene expression revealed key transcriptional regulators acting during myogenesis. In bovine skeletal muscle, by scRNA‐seq and scATAC‐seq analysis, different cell types such as adipocytes, endothelial cells, fibroblasts, lymphocytes, monocytes, pericyte cells and eight skeletal myogenic subpopulations were identified at the three developmental stages. The pseudotime trajectory exhibited a distinct sequential ordering for these myogenic subpopulations and eight distinct gene clusters were observed according to their expression pattern. Moreover, specifically expressed TFs (such as MSC, MYF5, MYOD1, FOXP3, ESRRA, BACH1, SIX2 and ATF4) associated with muscle development were predicted, and likely future transcriptional states of individual cells and the developmental dynamics of differentiation among neighbouring cells were predicted. CellChat analysis on the scRNA‐seq data set then classified many ligand–receptor pairs among these cell clusters, which were further categorized into significant signalling pathways, including BMP, IGF, WNT, MSTN, ANGPTL, TGFB, TNF, VEGF and FGF. Finally, scRNA‐seq and scATAC‐seq results were successfully integrated to reveal a series of specifically expressed TFs that are likely to be candidates for the promotion of cell fate transition during bovine skeletal muscle development. Overall, our results outline a single‐cell dynamic chromatin/transcriptional landscape for normal bovine skeletal muscle development; these provide an important resource for understanding the structure and function of mammalian skeletal muscle, which will promote research into its biology.

The male infertility of cattleyak resulted from spermatogenic arrest has greatly restricted the effective utilization of the heterosis from crossbreeding of cattle and yak. Based on our previous studies, the significant divergences of the transcriptomic and proteomic sequencing between yak and cattleyak prompt us to investigate the critical roles of microRNAs in post-transcriptional regulation of gene expression during spermatogenesis. TUNEL-POD analysis presented sharply decreased spermatogenic cell types and the increased apoptotic spermatogonia in cattleyak. The STA-PUT velocity sedimentation was employed to obtain spermatogonia and spermatocytes from cattle, yak and cattleyak and these spermatogenic cells were verified by the morphological and phenotypic identification. MicroRNA microarray showed that 27 differentially expressed miRNAs were simultaneously identified both in cattleyak vs cattle and in cattleyak vs yak comparisons. Further analysis revealed that the down-regulation of bta-let-7 families, bta-miR-125 and bta-miR-23a might impair the RA-induced differentiation of spermatogonia. Target gene analysis for differentially expressed miRNAs revealed that miRNAs targeted major players involved in vesicle-mediated transport, regulation of protein kinase activity and Pathways in cancer. In addition, spermatogonia transfection analysis revealed that the down-regulation of bta-miR-449a in the cattleyak might block the transition of male germ cells from the mitotic cycle to the meiotic program. The present study provided valuable information for future elucidating the regulatory roles of miRNAs involved in spermatogenic arrest of cattleyak.