Analysis of Chromatin Openness in Testicle Tissue of Yak and Cattle-Yak

Identifying low-density, ancestry-informative SNP markers through whole genome resequencing in Indian, Chinese, and wild yak

The current investigation was undertaken to elucidate the population-stratifying and ancestry-informative markers in Indian, Chinese, and wild yak populations using whole genome resequencing (WGS) analysis while employing various selection strategies (Delta, Pairwise Wright's Fixation Index-FST, and Informativeness of Assignment) and marker densities (5-25 thousand). The study used WGS data on 105 individuals from three separate yak cohorts i.e., Indian yak (n = 29), Chinese yak (n = 61), and wild yak (n = 15). Variant calling in the GATK program with strict quality control resulted in 1,002,970 high-quality and independent (LD-pruned) SNP markers across the yak autosomes. Analysis was undertaken in toolbox for ranking and evaluation of SNPs (TRES) program wherein three different criteria i.e., Delta, Pairwise Wright's Fixation Index-FST, and Informativeness of Assignment were employed to identify population-stratifying and ancestry-informative markers across various datasets. The top-ranked 5,000 (5K), 10,000 (10K), 15,000 (15K), 20,000 (20K), and 25,000 (25K) SNPs were identified from each dataset while their composition and performance was assessed using different criteria. The average genomic breed clustering of Indian, Chinese, and wild yak cohorts with full density dataset (105 individuals with 1,002,970 markers) was 81.74%, 80.02%, and 83.62%, respectively. Informativeness of Assignment criterion with 10K density emerged as the best combination for three yak cohorts with 86.94%, 96.46%, and 98.20% clustering for Indian, Chinese, and wild yak, respectively. There was an average increase of 7.56%, 22.72%, and 30.35% in genomic breed clustering scores of Indian, Chinese, and wild yak cohorts over the estimates of the original dataset. The selected markers showed overlap multiple protein-coding genes within a 10 kb window including ADGRB3, ANK1, CACNG7, CALN1, CHCHD2, CREBBP, GLI3, KHDRBS2, and OSBPL10. This is the first report ever on elucidating low-density SNP marker sets with population-stratifying and ancestry-informative properties in three yak groups using WGS data. The results gain significance for application of genomic selection using cost-effective low-density SNP panels in global yak species.

Whole transcriptome sequencing revealed the gene regulatory network of hypoxic response in yak Sertoli cells

Yaks live in the Qinghai-Tibet Plateau for a long time where oxygen is scarce, but can ensure the smooth development of testis and spermatogenesis. The key lies in the functional regulation of the Sertoli cells under hypoxia. In this study, we sequenced yak Sertoli cells cultured in normal oxygen concentration (Normoxia) and treated with low oxygen concentration (Hypoxia) by whole transcriptomics, and screened out 194 differentially expressed mRNAs (DEmRNAs), 934 differentially expressed LncRNAs (DELncRNAs) and 129 differentially expressed miRNAs (DEmiRNAs). GO and KEGG analysis showed that these differential genes were mainly concentrated in PI3K-AKT, MAPK, RAS, and other signaling pathways, and were associated with glucose metabolism, tight junction, steroid hormone synthesis, cell fusion, and immunity of yak Sertoli cells. We constructed the gene interaction network of yak Sertoli cells in hypoxia and screened out the relationship pairs related to glucose metabolism and tight junction. The results suggested that the changes in energy metabolism, tight junction, and immune regulation of yak Sertoli cells under hypoxia might provide favorable conditions for spermatogenesis. This study provides data for further study on the role of non-coding RNA in testis development and spermatogenesis of yak.

Integrative ATAC-seq and RNA-seq Analysis of the Longissimus Dorsi Muscle of Gannan Yak and Jeryak

Jeryak is the F1 generation of the cross between Gannan yak and Jersey cattle, which has the advantages of fast growth and high adaptability. The growth and development of skeletal muscle is closely linked to meat production and the quality of meat. However, the molecular regulatory mechanisms of muscle growth differences between Gannan yak and Jeryak analyzed from the perspective of chromatin opening have not been reported. In this study, ATAC-seq was used to analyze the difference of chromatin openness in longissimus muscle of Gannan yak and Jeryak. It was found that chromatin accessibility was more enriched in Jeryak compared to Gannan yak, especially in the range of the transcription start site (TSS) ± 2 kb. GO and KEGG enrichment analysis indicate that differential peak-associated genes are involved in the negative regulation of muscle adaptation and the Hippo signaling pathway. Integration analysis of ATAC-seq and RNA-seq revealed overlapping genes were significantly enriched during skeletal muscle cell differentiation and muscle organ morphogenesis. At the same time, we screened FOXO1, ZBED6, CRY2 and CFL2 for possible involvement in skeletal muscle development, constructed a genes and transcription factors network map, and found that some transcription factors (TFs), including YY1, KLF4, KLF5 and Bach1, were involved in skeletal muscle development. Overall, we have gained a comprehensive understanding of the key factors that impact skeletal muscle development in various breeds of cattle, providing new insights for future analysis of the molecular regulatory mechanisms involved in muscle growth and development.

mRNA-Seq of testis and liver tissues reveals a testis-specific gene and alternative splicing associated with hybrid male sterility in dzo

Alternative splicing (AS) plays an important role in the co-transcription and post-transcriptional regulation of gene expression during mammalian spermatogenesis. The dzo is the male F1 offspring of an interspecific hybrid between a domestic bull (Bos taurus ♂) and a yak (Bos grunniens ♀) which exhibits male sterility. This study aimed to identify the testis-specific genes and AS associated with hybrid male sterility in dzo. The iDEP90 program and rMATS software were used to identify the differentially expressed genes (DEG) and differential alternative splicing genes (DSG) based on RNA-seq data from the liver (n = 9) and testis (n = 6) tissues of domestic cattle, yak, and dzo. Splicing factors (SF) were obtained from the AmiGO2 and the NCBI databases, and Pearson correlation analysis was performed on the differentially expressed SFs and DSGs. We focused on the testis-specific DEGs and DSGs between dzo and cattle and yak. Among the top 3,000 genes with the most significant variations between these 15 samples, a large number of genes showed testis-specific expression involved with spermatogenesis. Cluster analysis showed that the expression levels of these testis-specific genes were dysregulated during mitosis with a burst downregulation during the pachynema spermatocyte stage. The occurrence of AS events in the testis was about 2.5 fold greater than in the liver, with exon skipping being the major AS event (81.89% to 82.73%). A total of 74 DSGs were specifically expressed in the testis and were significantly enriched during meiosis I, synapsis, and in the piRNA biosynthesis pathways. Notably, STAG3 and DDX4 were of the exon skipping type, and DMC1 was a mutually exclusive exon. A total of 36 SFs were significantly different in dzo testis, compared with cattle and yak. DDX4, SUGP1, and EFTUD2 were potential SFs leading to abnormal AS of testis-specific genes in dzo. These results show that AS of testis-specific genes can affect synapsis and the piRNA biosynthetic processes in dzo, which may be important factors associated with hybrid male sterility in dzo.

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.

Integration of ATAC-Seq and RNA-Seq Analysis to Identify Key Genes in the Longissimus Dorsi Muscle Development of the Tianzhu White Yak

During the postnatal stages, skeletal muscle development undergoes a series of meticulously regulated alterations in gene expression. However, limited studies have employed chromatin accessibility to unravel the underlying molecular mechanisms governing muscle development in yak species. Therefore, we conducted an analysis of both gene expression levels and chromatin accessibility to comprehensively characterize the dynamic genome-wide chromatin accessibility during muscle growth and development in the Tianzhu white yak, thereby elucidating the features of accessible chromatin regions throughout this process. Initially, we compared the differences in chromatin accessibility between two groups and observed that calves exhibited higher levels of chromatin accessibility compared to adult cattle, particularly within ±2 kb of the transcription start site (TSS). In order to investigate the correlation between alterations in chromatin accessible regions and variations in gene expression levels, we employed a combination of ATAC-seq and RNA-seq techniques, leading to the identification of 18 central transcriptional factors (TFs) and 110 key genes with significant effects. Through further analysis, we successfully identified several TFs, including Sp1, YY1, MyoG, MEF2A and MEF2C, as well as a number of candidate genes (ANKRD2, ANKRD1, BTG2 and LMOD3) which may be closely associated with muscle growth and development. Moreover, we constructed an interactive network program encompassing hub TFs and key genes related to muscle growth and development. This innovative approach provided valuable insights into the molecular mechanism underlying skeletal muscle development in the postnatal stages of Tianzhu white yaks while also establishing a solid theoretical foundation for future research on yak muscle development.

Comprehensive Analysis of miRNA and mRNA Expression Profiles during Muscle Development of the Longissimus Dorsi Muscle in Gannan Yaks and Jeryaks

A hybrid offspring of Gannan yak and Jersey cattle, the Jeryak exhibits apparent hybrid advantages over the Gannan yak in terms of production performance and other factors. The small non-coding RNAs known as miRNAs post-transcriptionally exert a significant regulatory influence on gene expression. However, the regulatory mechanism of miRNA associated with muscle development in Jeryak remains elusive. To elucidate the regulatory role of miRNAs in orchestrating skeletal muscle development in Jeryak, we selected longissimus dorsi muscle tissues from Gannan yak and Jeryak for transcriptome sequencing analysis. A total of 230 (DE) miRNAs were identified in the longissimus dorsi muscle of Gannan yak and Jeryak. The functional enrichment analysis revealed a significant enrichment of target genes from differentially expressed (DE)miRNAs in signaling pathways associated with muscle growth, such as the Ras signaling pathway and the MAPK signaling pathway. The network of interactions between miRNA and mRNA suggest that some (DE)miRNAs, including miR-2478-z, miR-339-x, novel-m0036-3p, and novel-m0037-3p, played a pivotal role in facilitating muscle development. These findings help us to deepen our understanding of the hybrid dominance of Jeryaks and provide a theoretical basis for further research on the regulatory mechanisms of miRNAs associated with Jeryak muscle growth and development.

Long-stranded non-coding RNAs temporal-specific expression profiles reveal longissimus dorsi muscle development and intramuscular fat deposition in Tianzhu white yak

The process of muscle development and intramuscular fat (IMF) deposition is quite complex and controlled by both mRNAs and ncRNAs. Long-stranded non-coding RNAs (LncRNAs) are involved in various biological processes in mammals while also playing a critical role in muscle development and fat deposition. In the present study, RNA-Seq was used to comprehensively study the expression of lncRNAs and mRNAs during muscle development and intramuscular fat deposition in postnatal Tianzhu white yaks at three stages, including 6 mo of age (calve, n = 6), 30 mo of age (young cattle, n = 6) and 54 mo of age (adult cattle, n = 6). The results indicated that a total of 2,101 lncRNAs and 20,855 mRNAs were screened across the three stages, of which the numbers of differential expression (DE) lncRNAs and DE mRNAs were 289 and 1,339, respectively, and DE lncRNAs were divided into eight different expression patterns based on expression trends. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that some DE mRNAs overlapped with target genes of lncRNAs, such as NEDD4L, SCN3B, AGT, HDAC4, DES, MYH14, KLF15 (muscle development), ACACB, PCK2, LIPE, PIK3R1, PNPLA2, and MGLL (intramuscular fat deposition). These DE mRNAs were significantly enriched in critical muscle development and IMF deposition-related pathways and GO terms, such as AMPK signaling pathway, PI3K-Akt signaling pathway, PPAR signaling pathway, etc. In addition, lncRNA-mRNA co-expression network analysis revealed that six lncRNAs (MSTRG.20152.2, MSTRG.20152.3, XR_001351700.1, MSTRG.8190.1, MSTRG.4827.1, and MSTRG.11486.1) may play a major role in Tianzhu white yak muscle development and lipidosis deposition. Therefore, this study enriches the database of yak lncRNAs and could help to further explore the functions and roles of lncRNAs in different stages of muscle development and intramuscular fat deposition in the Tianzhu white yak.