Spermatogenesis is a continuous process in which functional sperm are produced through a series of mitotic and meiotic divisions and morphological changes in germ cells. The aberrant development and fate transitions of spermatogenic cells cause hybrid sterility in mammals. Cattle-yak, a hybrid animal between taurine cattle (Bos taurus) and yak (Bos grunniens), exhibits male-specific sterility due to spermatogenic failure. In the present study, we performed single-cell RNA sequencing analysis to identify differences in testicular cell composition and the developmental trajectory of spermatogenic cells between yak and cattle-yak. The composition and molecular signatures of spermatogonial subtypes were dramatically different between these 2 animals, and the expression of genes associated with stem cell maintenance, cell differentiation and meiotic entry was altered in cattle-yak, indicating the impairment of undifferentiated spermatogonial fate decisions. Cell communication analysis revealed that signaling within different spermatogenic cell subpopulations was weakened, and progenitor spermatogonia were unable to or delayed receiving and sending signals for transformation to the next stage in cattle-yak. Simultaneously, the communication between niche cells and germ cells was also abnormal. Collectively, we obtained the expression profiles of transcriptome signatures of different germ cells and testicular somatic cell populations at the single-cell level and identified critical regulators of spermatogonial differentiation and meiosis in yak and sterile cattle-yak. The findings of this study shed light on the genetic mechanisms that lead to hybrid sterility and speciation in bovid species.

Hybrids between closely related but genetically incompatible species are often inviable or sterile. Cattle-yak, an interspecific hybrid of yak and cattle, exhibits male-specific sterility, which limits the fixation of its desired traits and prevents genetic improvement in yak through crossbreeding. Transcriptome profiles of testicular tissues have been generated in cattle, yak, and cattle-yak; however, the genetic variations underlying differential gene expression associated with hybrid sterility have yet to be elucidated. We detected differences in the cellular composition and gene expression of testes from yak and cattle-yak at 3 mo of age, 10 mo of age, and adulthood. Histological analysis revealed that the most advanced germ cells were gonocytes (prospermatogonia) at 3 mo and spermatocytes at 10 mo. Complete spermatogenesis occurred in the seminiferous tubules of adult yak, whereas only spermatogonia and a limited number of spermatocytes were detected in the testis of adult cattle-yak. Transcriptome analysis revealed 180, 6,310, and 6,112 differentially expressed genes (DEG) in yak and cattle-yak at each stage, respectively. Next, we examined the spermatogenic cell types in the backcross generation (BC1) and detected the appearance of round spermatids, indicating the partial recovery of spermatogenesis in these animals. Compared with those in cattle-yak, 272 DEG were identified in the testes of BC1 animals. Notably, we discovered that the expression of X chromosome-linked genes was upregulated in the testis of cattle-yak compared with yak, suggesting a possible abnormality in the process of meiotic sex chromosome inactivation in hybrid animals. We next screened DEG harboring structural variations (SV) and identified a list of SV genes associated with spermatogonial development, meiotic recombination, and double-strand break repair. Furthermore, we found that the SV genes ESCO2 (establishment of sister chromatid cohesion N-acetyltransferase 2) and BRDT (bromodomain testis associated) may be involved in meiotic arrest of cattle-yak spermatocytes. Overall, our research provides a valuable database for identifying structural variant loci that contribute to hybrid sterility.

The goal of this study was to determine the metabolism of multiparous female yaks during the late perinatal period and identify its effects on reproductive recovery in order to explain the low reproduction rate of yaks. Eight multiparous female yaks were randomly selected as the sample, and serum was collected from the yaks every 7 days from the day of delivery until 28 days after the delivery (five time points). The presence of serum metabolic profiles and reproductive hormones was identified using ELISA. The key metabolites were identified using liquid chromatography–mass spectrometry, and a dynamic metabolic network representation was created using bioinformatics analysis. A total of 117 different metabolites were identified by calculating the fold change of the metabolite expression at each time point. The dynamic metabolic network was created to represent the activities of the key metabolites, metabolic indexes and reproductive hormones. The initial efficiency of the glucose metabolism in the late perinatal period was found to be low, but it increased during the final period. The initial efficiencies of the lipid and amino acid metabolisms were high but decreased during the final period. We inferred that there was a postpartum negative energy balance in female yaks and that the synthesis and secretion of estrogen were blocked due to an excessive fatty acid mobilization. As a result, the reproductive hormone synthesis and secretion were maintained at a low level in the late perinatal period, and this was the main reason for the delayed recovery of the reproductive function postpartum. However, the specific mechanism needs to be further verified.

Spermatogenesis is a highly orchestrated germ cell differentiation process involving the dynamic regulation of cell fate transitions. Dissecting the molecular landscapes of spermatogenic cell types is crucial for identifying fertility-related problems and improving the reproductive performance of farm animals. Here, we conducted transcriptomic and chromosome spreading across meiotic stages of testicular cells from taurine cattle (Bos taurus), yak (Bos grunniens), and their hybrid progenies to describe the transcriptional landscape of normal spermatogenesis and identify potential regulators that are involved in hybrid sterility. The results revealed seven types of spermatogonia, ten spermatocytes and ten types of spermatids in the cattle or yak testes. In sharp contrast, the testes of the cattle–yak hybrids contained only seven spermatogonial subtypes and six types of spermatocytes. Notably, the arrest of spermatocytes at the diplotene-to-diakinesis transition was accompanied by defects in double-strand break repair. In the testes of backcrossed offspring, spermatogenic arrest was partially rescued, and round spermatozoa were produced. By performing joint analysis, we identified 115 genes that exhibited differential protein abundance in spermatocytes of cattle-yak. Among them, 24 genes carrying genomic structural variations were differentially expressed in spermatocytes of cattle-yak but recovered in those of backcrossed offspring. This work provides important insights into spermatogenesis in large animals and serves as a valuable resource for identifying the factors determining reproductive isolation.

Meiotic recombination is initiated by the formation of programmed DNA double-strand breaks during spermatogenesis. PRDM9 determines the localization of recombination hotspots by interacting with several protein complexes in mammals. The function of PRDM9 is not well understood during spermatogenesis in mice or yaks. In this study, we applied yeast two-hybrid assays combined with next-generation sequencing techniques to screen the complete set of PRDM9-interacting proteins and explore its novel functions in yak spermatogenesis. Our results showed that 267 PRDM9-interacting proteins were identified. The gene ontology (GO) analysis of the interacting proteins revealed that the GO terms were primarily associated with spermatogenesis, positive regulation of double-strand break repair via homologous recombination, RNA splicing, the ubiquitin-dependent ERAD pathway, and other biological processes. MKX and PDCD5 were verified to be strongly interacting with PRDM9 and expressed in prophase I of meiosis in both mouse and yak testes. The localizations of RNA splicing genes including THOC5, DDX5, and XRCC6 were expressed in spermatocytes. Cattleyak is the hybrid offspring of a yak and a domestic cow, and the male offspring are sterile. The gene expression of the interacting proteins was also examined in the sterile male hybrid of yak and cattle. Among the 58 detected genes, 55 were downregulated in cattleyak. In conclusion, we established a complete PRDM9 interaction network, and a novel function of PRDM9 was identified, which will further promote our understanding of spermatogenesis. It also provides new insights for the study of hybrid male sterility.

To investigate prenatal muscle satellite cell (MuSC) development and the associated epigenetic modifications in yak. Here, we conducted morphological and protein co-localization analyses of fetal longissimus dorsi muscle at various developmental stages using histology and immunofluorescence staining methods. Our study observed that primary muscle fibers began forming at 40 days of gestation, fully developed by 11 weeks, and secondary muscle fibers were predominantly formed by around 105 days. Throughout development, MuSCs were mainly located between the muscle fiber membrane and the basement membrane, acting as a reserve for the stem cell pool. MuSCs appeared within myotubes only during critical phases of primary and secondary muscle fiber formation. The proliferation of MuSCs gradually decreases until birth. MuSCs with 5mC modification show a trend of increasing first and then decreasing. MuSCs with 5hmC modification also present a dynamic change trend. The 41st day and 11th week are the critical periods for the changes of both. From the 11th week to around the 110th day of gestation, the modification effect of histone H3K4me3 is crucial for MuSCs during the development of the fetal longissimus dorsi muscle. Combined, our data identify key time points for yak fetal skeletal muscle growth and development and demonstrate that DNA methylation and histone modifications in MuSCs are closely related to this process, offering a valuable basis for future research into the molecular mechanisms underlying yak muscle development.