Comparative study of growth performance, nutrient digestibility, and ruminal and fecal bacterial community between yaks and cattle-yaks raised by stall-feeding

The experiment was conducted to compare the growth performance, rumen fermentation, nutrient digestibility, and ruminal and fecal bacterial community between yaks and cattle-yaks. Ten male yaks (36-month-old) were used as the yak (YAK) group and 10 male cattle-yaks with similar age were selected as the cattle-yak (CAY) group. All the animals were fed same ration and the experiment lasted for 60 days. The results showed that the average daily gain and dry matter intake of CAY group were higher (P < 0.05) than those of YAK group. The ruminal concentrations of total volatile fatty acids, acetate, and butyrate were higher (P < 0.05) in CAY group than those in YAK group. However, the neutral detergent fiber and acid detergent fiber digestibility exhibited an opposite between two groups. In the rumen, the relative abundances of Prevotella 1 and Prevotellaceae UCG-001 were higher (P < 0.05) and Succiniclasticum and Butyrivibrio 2 were lower (P < 0.05) in YAK group compared to CAY group. In the feces, the unclassified Lachnospiraceae, Lachnospiraceae NK4A136 group, and Lachnospiraceae AC2044 group were significantly enriched (P < 0.05) in YAK group, whereas the Ruminococcaceae UCG-010, Ruminococcaceae UCG-013, and Succiniclasticum were significantly enriched (P < 0.05) in CAY group. Overall, under the same diet, the yaks have higher fiber utilization and cattle-yaks have higher energy utilization.

Expression of Rad51 and the histo-morphological evaluation of testis of the sterile male cattle-yak

Meiotic recombination is key to the repair of DNA double-strand break damage, provide a link between homologs for proper chromosome segregation as well as ensure genetic diversity in organisms. Defects in recombination often lead to sterility. The ubiquitously expressed Rad51 and the meiosis-specific DMC1 are two closely related recombinases that catalyze the key strand invasion and exchange step of meiotic recombination. This study cloned and sequenced the coding region of cattle-yak Rad51 and determined its mRNA and protein expression levels, evaluated its molecular and evolutionary relationship as well as evaluated the histo-morphological structure of testes in the yellow cattle, yak and the sterile cattle-yak hybrid. The Rad51 gene was amplified using PCR, cloned and sequenced using testicular cDNA from yak and cattle-yak. Real-time PCR was used to examine the expression levels of Rad51/DMC1 mRNA in the cattle, yak and cattle-yak testis while western blotting, immunofluorescence and immunohistochemistry were used to assess the protein expression and localization of Rad51/DMC1 protein in the testicular tissue sections. The results revealed that the mRNA and protein expression of Rad51 and DMC1 are extremely low in the male cattle-yak testis with a corresponding higher incidence of germ cell apoptosis. There was also thinning of the germinal epithelium possibly due to the depletion of the germ cells leading to the widening of the lumen area of the cattle-yak seminiferous tubule. Our findings provide support for the hypothesis that the low expression of Rad51 and DMC1 may contribute to the male hybrid sterility in the cattle-yak.

Single-cell transcriptome analysis and in vitro differentiation of testicular cells reveal novel insights into male sterility of the interspecific hybrid cattle-yak

BACKGROUND

Interspecific hybridization plays vital roles in enriching animal diversity, while male hybrid sterility (MHS) of the offspring commonly suffered from spermatogenic arrest constitutes the postzygotic reproductive isolation. Cattle-yak, the hybrid offspring of cattle (Bos taurus) and yak (Bos grunniens) can serve as an ideal MHS animal model. Although meiotic arrest was found to contribute to MHS of cattle-yak, yet the cellular characteristics and developmental potentials of male germline cell in pubertal cattle-yak remain to be systematically investigated.

RESULTS

Single-cell RNA-seq analysis of germline and niche cell types in pubertal testis of cattle-yak and yak indicated that dynamic gene expression of developmental germ cells was terminated at late primary spermatocyte (meiotic arrest) and abnormal components of niche cell in pubertal cattle-yak. Further in vitro proliferation and differentially expressed gene (DEG) analysis of specific type of cells revealed that undifferentiated spermatogonia of cattle-yak exhibited defects in viability and proliferation/differentiation potentials.

CONCLUSION

Comparative scRNA-seq and in vitro proliferation analysis of testicular cells indicated that not only meiotic arrest contributed to MHS of cattle-yak. Spermatogenic arrest of cattle-yak may originate from the differentiation stage of undifferentiated spermatogonia and niche cells of cattle-yak may provide an adverse microenvironment for spermatogenesis.

Supplementary feeding of cattle-yak in the cold season alters rumen microbes, volatile fatty acids, and expression of SGLT1 in the rumen epithelium

Cattle-yak, a hybrid offspring of yak (Bos grunniens) and cattle (Bos taurus), inhabit the Qinghai-Tibet Plateau at an altitude of more than 3,000 m and obtain nutrients predominantly through grazing on natural pastures. Severe shortages of pasture in the cold season leads to reductions in the weight and disease resistance of grazing cattle-yak, which then affects their production performance. This study aimed to investigate the effect of supplementary feeding during the cold season on the rumen microbial community of cattle-yak. Six cattle-yak (bulls) were randomly divided into two groups-"grazing + supplementary feeding" (G+S) (n = 3) and grazing (G) (n = 3)-and rumen microbial community structure (based on 16S rRNA sequencing), volatile fatty acids (VFAs), and ruminal epithelial sodium ion-dependent glucose transporter 1 (SGLT1) expression were assessed. There were significant differences in the flora of the two groups at various taxonomic classification levels. For example, Bacteroidetes, Rikenellaceae, and Rikenellaceae_RC9_gut_group were significantly higher in the G+S group than in the G group (P < 0.05), while Firmicutes and Christensenellaceae_R-7_group were significantly lower in the G+S group than in the G group (P < 0.05). Kyoto Encyclopedia of Genes and Genomes (KEGG) and Clusters of Orthologous Groups (COG) analyses revealed that functions related to carbohydrate metabolism and energy production were significantly enriched in the G+S group (P < 0.05). In addition, the concentration of total VFAs, along with concentrations of acetate, propionate, and butyrate, were significantly higher in the G+S group than in the G group (P < 0.05). Furthermore, SGLT1 expression in ruminal epithelial tissue was significantly lower in the G+S group (P < 0.01). Supplementary feeding of cattle-yak after grazing in the cold season altered the microbial community structure and VFA contents in the rumen of the animals, and decreased ruminal epithelial SGLT1 expression. This indicated that supplementary feeding after grazing aids rumen function, improves adaptability of cattle-yak to the harsh environment of the Qinghai-Tibet Plateau, and enhances ability of the animals to overwinter.

Quantitative proteomic analysis of cattle-yak and yak longissimus thoracis provides insights into the differential mechanisms of meat quality

In this study, proteins of cattle-yak longissimus thoracis (CYLT) and yak longissimus thoracis (YLT) were compared using tandem mass tag-labeled quantitative proteomic analysis. A total of 157 proteins were screened as differentially abundant proteins (DAPs) derived from 1551 quantitative proteins. Bioinformatics analysis revealed that the upregulated DAPs in CYLT were mainly involved in energy metabolism, oxidative stress, muscle fiber structure, and extracellular matrix (ECM), while the downregulated DAPs were mainly involved in energy metabolism and ECM function. The upregulated myoglobin, downregulation of NADH dehydrogenase, and upregulation of cytochrome oxidase indicated that CYLT initiates compensatory regulation in response to hypoxic high-altitude environments. Two differentially abundant myosins and five collagens suggested that CYLT and YLT may have distinct differences in the assembly structure of muscle fibers and connective tissue. These differences in energy metabolism and muscle structure will inevitably affect the postmortem physiology of "muscle to meat" and consequently the meat qualities. Therefore, our results will provide important clues to gain insight into the potential causes of meat quality differences between cattle-yak and yak based on high-altitude response.

Full-length transcriptome sequencing of the longissimus dorsi muscle of yak and cattle-yak using nanopore technology

Short-read RNA sequencing has been used to sequence the transcriptome of the skeletal muscle of yak and cattle-yak; however, full-length transcripts cannot be obtained and alternative splicing (AS) events cannot be inferred using this sequencing approach. Here, we used Oxford Nanopore Technologies (ONT) full-length sequencing to sequence the transcriptome of the longissimus dorsi of yak and cattle-yak. A total of 20,323 novel genes and 172,870 novel transcripts were identified, and 159,700 novel transcripts were successfully annotated. A total of 157,812 AS events, 58,073 simple sequence repeats, 57,468 complete open reading frames, 2296 transcription factors, and 20,404 lncRNAs were detected. Differentially expressed transcripts (DETs) in the longissimus dorsi muscle of yak and cattle-yak were involved in the MAPK and JAK-STAT signaling pathways related to muscle development and growth. Protein-protein interaction analysis of DETs suggested that TNNI2 might make a major contribution to differences in muscle growth and meat quality traits between yak and cattle-yak. The results have enriched the transcriptome data of dorsal muscles, providing new ideas for the study of transcriptional regulation processes, and also providing useful information for the production of higher yields of yak meat.

Comprehensive analysis of differentially expressed mRNAs, circRNAs, and miRNAs and their ceRNA network in the testis of cattle-yak, yak, and cattle

Cattle-yak is a hybrid offspring resulting from the crossbreeding of yak and cattle, and it exhibits substantial heterosis in production performance. However, male sterility in cattle-yak remains a concern. Reports suggest that noncoding RNAs are involved in the regulation of spermatogenesis. Therefore, in this study, we comprehensively compared testicular transcription profiles among cattle, yak, and cattle-yak. Numerous differentially expressed genes (DEGs), differentially expressed circRNAs (DECs), and differentially expressed miRNAs (DEMs) were identified in the intersection of two comparison groups, namely cattle versus cattle-yak and yak versus cattle-yak, with the number of DEGs, DECs, and DEMs being 4968, 360, and 59, respectively. The DEGs in cattle-yaks, cattle, and yaks were mainly associated with spermatogenesis, male gamete generation, and sexual reproduction. Concurrently, GO and KEGG analyses indicated that DEC host genes and DEM source genes were involved in the regulation of spermatogenesis. The construction of a potential competing endogenous RNA network revealed that some differentially expressed noncoding RNAs may be involved in regulating the expression of genes related to testicular spermatogenesis, including miR-423-5p, miR-449b, miR-34b/c, and miR-15b, as well as previously unreported miR-6123 and miR-1306, along with various miRNA-circRNA interaction pairs. This study serves as a valuable reference for further investigations into the mechanisms underlying male sterility in cattle-yaks.

Whole-transcriptome analysis of longissimus dorsi muscle in cattle-yaks reveals the regulatory functions of ADAMTS6 gene in myoblasts

Cattle-yak, which is the hybrid F1 generation of cattle and yak, demonstrates better production performance compared to yak. However, there is limited research on the molecular mechanisms responsible for the muscle development of cattle-yak. To address this knowledge gap, a comprehensive transcriptomic survey of the longissimus dorsi muscle in cattle-yak was conducted. Three transcript types, namely lncRNAs, miRNAs, and circRNAs, along with protein-coding genes were characterized at two developmental stages (6 m, 18 m) of cattle-yak. The results revealed significant enrichment of these transcripts into pathways related to myoblast differentiation and muscle development signaling. Additionally, the study identified the TCONS00024465/circHIPK3-bta-miR-499-ADAMTS6 regulatory network, which may play a crucial role in the muscle development of cattle-yak by combining the transcriptome data of yak and constructing the ceRNA co-expression network. HEK 293 T cells were used to validate that TCONS00024465 and circHIPK3 are located upstream of bta-miR-499, and can competitively bind to bta-miR-499 as ceRNA. The study also verified that ADAMTS6 regulates skeletal muscle development by inhibiting myoblast proliferation, promoting myoblast differentiation, and positively regulating the apoptosis of myoblasts. Taken together, this study provides new insights into the advantages of cattle-yak production performance and offers a molecular basis for further research on muscle development.

Comparative RNA-Seq analysis of differentially expressed genes in the sertoli cells of yak and cattle-yak

BACKGROUND

To study the problem of male sterility of cattle-yak and improve the yak crossbreeding, this study obtained the testicular Sertoli cells of yak and cattle-yak and compared the differences in transcriptome levels between the two bovine species. The testicular tissues of 3 healthy male cattle-yaks and 3 F1 generation male yaks were collected at the age of 24 months. The Sertoli cells were isolated after enzymatic digestion, differential adhesion and starvation treatment. DATA-4 and SOX9 immunofluorescence staining were used to identify the cell type. Sertoli cells were subjected to transcriptome sequencing, GO analysis, KEGG analysis and differentially expressed gene were validated by RT-qPCR and Western blotting.

RESULTS

The study successfully isolated and purified Sertoli cells of yak and cattle-yak. The transcriptome sequencing data were compared, analyzed and annotated. Compared to yak Sertoli cells, 6592 differentially expressed genes were identified, with 3007 genes upregulated and 3585 genes downregulated in cattle-yak Sertoli cells. GO analysis suggested that the upregulated genes might be mainly involved in processes such as translation, peptide biosynthetic process, amide biosynthetic process, peptide metabolic process, ribosome, cytoplasmic part, structural constituent of ribosome, structural molecule activity, endomembrane system, protein kinase activity, and phosphotransferase activity. The downregulated genes appeared to be primarily involved in protein phosphorylation, phosphorylation, endomembrane system, protein kinase activity, and phosphotransferase activity. KEGG analysis compared differential genes across 316 pathways, with 8 pathways showing significant enrichment. The upregulated pathways were potentially enriched in cattle-yak Sertoli cells, including ribosome, thermogenesis, and oxidative phosphorylation, while the downregulated pathways seemed to be significantly enriched in adherens junction, mTOR signaling pathway, AMPK signaling pathway, FoxO signaling pathway, and focal adhesion. Compared with yak Sertoli cells, ISOC2, RPL27A and FISI were highly expressed in cattle-yak, as confirmed by RT-qPCR analysis. PDPN, SORBS2, TF, PLSCR1, TJP2, KIF2C, ITGA3, SMTNL2, DSP, ADGRG1, DDR1, GSK3A, RBBP6, ZC3H15 and Claudin 11 showed low expression levels in cattle-yak.

CONCLUSIONS

Compared with yak Sertoli cells, the expression of genes related to protein activation, cell function, and membranous organelle composition in cattle-yak Sertoli cells appeared to be abnormal. The potential defects in cattle-yak Sertoli cells may hinder the creation of a suitable environment for spermatogenesis, which could be one of the factors contributing to male cattle-yak sterility. Claudin-11 might be a potentially important gene for further research into cattle-yak male sterility.

Complete characterization of the yak testicular development using accurate full-length transcriptome sequencing

Alternative splicing is a prevalent phenomenon in testicular tissues. Due to the low assembly accuracy of short-read RNA sequencing technology in analyzing post-transcriptional regulatory events, full-length (FL) transcript sequencing is highly demanded to accurately determine FL splicing variants. In this study, we performed FL transcriptome sequencing of testicular tissues from 0.5, 1.5, 2.5, and 4-year-old yaks and 4-year-old cattle-yaks using Oxford Nanopore Technologies. The obtained sequencing data were predicted to have 47,185 open reading frames (ORFs), including 26,630 complete ORFs, detected 7645 fusion transcripts, 15,355 alternative splicing events, 25,798 simple sequence repeats, 7628 transcription factors, and 35,503 long non-coding RNAs. A total of 40,038 novel transcripts were obtained from the sequencing data, and the proportion was almost close to the number of known transcripts identified. Structural analysis and functional annotation of these novel transcripts resulted in the successful annotation of 9568 transcripts, with the highest and lowest annotation numbers in the Nr and KOG databases, respectively. Weighted gene co-expression network analysis revealed the key regulatory pathways and hub genes at various stages of yak testicular development. Our findings enhance our comprehension of transcriptome complexity, contribute to genome annotation refinement, and provide foundational data for further investigations into male sterility in cattle-yaks.

Single-cell RNA sequencing and UPHLC-MS/MS targeted metabolomics offer new insights into the etiological basis for male cattle-yak sterility

The variety of species can be efficiently increased by interspecific hybridization. However, because the males in the hybrid progeny are usually sterile, this heterosis cannot be employed when other cattle and yaks are hybridized. While some system-level studies have sought to explore the etiological basis for male cattle-yak sterility, no systematic cellular analyses of this phenomenon have yet been performed. Here, single-cell RNA sequencing and UPHLC-MS/MS targeted metabolomics methods were used to study the differences in testicular tissue between 4-year-old male yak and 4-year-old male cattle-yak, providing new and comprehensive insights into the causes of male cattle-yak sterility. Cattle-yak testes samples detected 6 somatic cell types and one mixed germ cell type. Comparisons of these cell types revealed the more significant differences in Sertoli cells (SCs) and [Leydig cells and myoid cells (LCs_MCs)] between yak and cattle-yak samples compared to other somatic cell clusters. Even though the LCs and MCs from yaks and cattle-yaks were derived from the differentiation of the same progenitor cells, a high degree of overlap between LCs and MCs was observed in yak samples. Still, only a small overlap between LCs and MCs was observed in cattle-yak samples. Functional enrichment analyses revealed that genes down-regulated in cattle-yak SCs were primarily enriched in biological activity, whereas up-regulated genes in these cells were enriched for apoptotic activity. Furthermore, the genes of up-regulated in LCs_MCs of cattle-yak were significantly enriched in enzyme inhibitor and molecular function inhibitor activity. On the other hand, the genes of down-regulated in these cells were enriched for signal receptor binding, molecular function regulation, positive regulation of biological processes, and regulation of cell communication activity. The most significant annotated differences between yak and cattle-yak LCs_MCs were associated with cell-to-cell communication. While yak LCs_MCs regulated spermatogenic cells at spermatogonia, spermatocyte, and spermatid levels, no such relationships were found between cattle-yak LCs_MCs and germ cells. This may suggest that the somatic niche in male cattle-yak testes is a microenvironment that is ultimately not favorable for spermatogenesis.