Background and aim Mitochondrial genome has aseries of characteristics such as simple structure, no recombination, maternalinheritance, stable structure, fast evolution rate, and high copy number. Moreover, it is easy to be sequenced,contains high-resolution phylogenetic information, and exists in a wide rangeof taxa. Therefore, it is widely used in the study of biological phylogeny. Atpresent, phylogenetic studies focus mainly on D-loop region, cytochrome b gene,and protein-coding sequence. Phylogenetic studies using the mitochondrialcomplete sequence are rarely reported in yak. Therefore, the present studyaimed to construct phylogenetic tree using yak mitochondrial complete sequenceand compare the subsequent results with previous findings obtained usingpartial sequences. Results Complete mitochondrial sequences of five yakpopulations from Qinghai and Xinjiang were obtained. The mitotype diversity ofthe five populations was Xueduo yak (0.992 ± 0.015), Pamir yak (0.990 ± 0.014),Yushu yak (0.963 ± 0.033), Qilian yak (0.948 ± 0.036), and Huanhu yak (0.905 ±0.048), which showed a higher mitotype diversity compared with other breeds fromthe previous reports, including Jiulong yak, Maiwa yak, Zhongdian yak, andTianzhu yak. A total of 78 mitotypes were obtained from 111 individuals. Amongthese, Yushu yak, Huanhu yak, Xueduo yak, and Qilian yak all shared mitotypes,but the Pamir yak did not share mitotypes with these four populations.Phylogenetic analysis showed that yak populations were separable into threedistinct branches. The analysis identified a new phylogenetic branch containingboth wild and domestic yaks. The 155 mitotypes found in 206 individuals weredivided into 3 haplogroups by mitotype clustering. Thehaplogroup was not associated with the geographical distribution of yaks. Theyaks in the same population or the same ecological environment were distributedin different haplogroups. Among the threehaplogroups, haplogroup A and haplogroup B showed a star-shaped distribution ofmitotypes. The central mitotypes were widely distributed and had a highfrequency. Conclusions Thegenetic diversity of yaks in Qinghai was high. Both domestic and wild yaks clusteredinto three branches.

The potential regulatory role of N6-methyladenosine (m6A), the most prominent mRNA modification in eukaryotes, has recently been identified in mammals, plants, and yeast. However, whether and how m6A methylation is involved in sexual maturation in mammals remains largely unexplored. In this study, testicular tissue was obtained from yaks before and after sexual maturation, and m6A maps were generated via preliminary experiments and methylated RNA immunoprecipitation sequencing. Only spermatogonial cells and a few primary spermatocytes were observed in the testicular tissue of yaks before sexual maturation, while spermatogenic cells at different stages of maturity could observed after sexual maturation. Experiments examining the expression of methylation-related enzymes and overall methylation levels showed that the methylation levels in yak testes increased after sexual maturation. Overall, 1,438 methylation peaks were differentially expressed before and after sexual maturation; 1,226 showed significant up-regulation and 212 showed significant down-regulation after sexual maturation. Annotation analysis showed that the differential methylation peaks were most commonly concentrated in the exon region, followed by the 3′UTR and finally the 5′UTR region. KEGG pathway analysis demonstrated that homologous recombination, the Notch signaling pathway, growth hormone synthesis, and other signaling pathways may be involved in testicular development and maturation in yaks. Levels of most m6A modifications were positively correlated with mRNA abundance, suggesting that m6A plays a regulatory role in mammalian sexual maturation. To our knowledge, this is the first report of an m6A transcriptional map of the yak testes, and our study lays the foundation for elucidating the function of m6A in the development of yak testes.

No abstract

The ovary as one of the most dynamic organs produces steroids to orchestrate female secondary sexual characteristics, harbors ovarian reserve for oocytes, releases mature oocytes for fertilization, and maintains pregnancy. Yak (Bos grunniens) is the only bovid animal that can adapt to the harsh climatic conditions on the Qinghai‐Tibetan Plateau (altitudes of over 3000 m above sea level). However, the cellular atlas is composed of oocytes and other somatic cells, and their individual molecular characteristics remain to be elucidated in the yak ovary. Here, single‐cell RNA sequencing (scRNA‐seq) was performed to delineate the molecular signature of various cell types in the yak ovarian cortex. A cellular atlas of yak ovarian cortex was constructed successfully on the basis of the differentially expressed genes (DEGs) from the distinct cell types and their functional enrichment analysis, comprising endothelial cells, nature kill cells, stromal cells, smooth muscle cells, oocytes, macrophages, epithelial cells, and granulosa cells. Meanwhile, the signature genes were determined based on their expression specificity in each cell type. A cell‐to‐cell communication network was built in light of the differentially overexpressed ligand and receptor genes from each cell type. Further, the oocytes were subdivided into four subtypes based on their individual DEGs and the functional enrichment of the DEGs. FST and TOP2A were identified as maker genes for oocytes by immunostaining in the yak ovarian cortex. The cellular atlas reveals the biological characteristics of the ovarian cortex at the cellular molecular level and provides insights into female reproductive biology via cellular communications in the yak.

N6-methyladenosine (m6A) is the most common form of eukaryotic mRNA modification, and it has been shown to exhibit broad regulatory activity in yeast, plants, and mammals. The specific role of m6A methylation as a regulator of spermatogenesis, however, has yet to be established. In this experiment, through a series of preliminary studies and methylated RNA immunoprecipitation sequencing, the m6A map of cattle-yak testicular tissue was established as a means of exploring how m6A modification affects cattle-yak male infertility. Cattle-yak testis tissues used in this study were found to contain sertoli cells and spermatogonia. Relative to sexually mature yak samples, those isolated from cattle-yak testis exhibited slightly reduced levels of overall methylation, although these levels were significantly higher than those in samples from pre-sexually mature yaks. Annotation analyses revealed that differentially methylated peaks were most concentrated in exonic regions, with progressively lower levels of concentration in the 3'-untranslated region (UTR) and 5'-UTR regions. To further explore the role of such m6A modification, enrichment analyses were performed on differentially methylated and differentially expressed genes in these samples. For the cattle-yaks vs. 18-months-old yaks group comparisons, differentially methylated genes were found to be associated with spermatogenesis-related GO terms related to the cytoskeleton and actin-binding, as well as with KEGG terms related to the regulation of the actin cytoskeleton and the MAPK signaling pathway. Similarly, enrichment analyses performed for the cattle-yaks vs. 5-years-old yaks comparison revealed differentially methylated genes to be associated with GO terms related to protein ubiquitination, ubiquitin ligase complexes, ubiquitin-dependent protein catabolism, and endocytotic activity, as well as with KEGG terms related to apoptosis and the Fanconi anemia pathway. Overall, enrichment analyses for the cattle-yaks vs. 18-months-old yaks comparison were primarily associated with spermatogenesis, whereas those for the cattle-yaks vs. 5-years-old yaks comparison were primarily associated with apoptosis.

Yaks (Bos grunniens) are the only bovine species that adapt well to the harsh high-altitude environment in the Qinghai-Tibetan plateau. However, the reproductive adaptation to the climate of the high elevation remains to be elucidated. Cell composition and molecular characteristics are the foundation of normal ovary function which determines reproductive performance. So, delineating ovarian characteristics at a cellular molecular level is conducive to elucidating the mechanism underlying the reproductive adaption of yaks. Here, the single-cell RNA-sequencing (scRNA-seq) was employed to depict an atlas containing different cell types with specific molecular signatures in the yak ovary. The cell types were identified on the basis of their specifically expressed genes and biological functions. As a result, a cellular atlas of yak ovary was established successfully containing theca cells, stromal cells, endothelial cells, smooth muscle cells, natural killer cells, macrophages, and proliferating cells. A cell-to-cell communication network between the distinct cell types was constructed. The theca cells were clustered into five subtypes based on their biological functions. Further, CYP11A1 was confirmed as a marker gene for the theca cells by immunofluorescence staining. Our work reveals an ovarian atlas at the cellular molecular level and contributes to providing insights into reproductive adaption in yaks.

Male-derived sterility in cattle-yaks, a hybrid deriving from yak and cattle, is a challenging problem. This study compared and analyzed the histomorphological differences in testis between sexually mature yak and cattle-yak, and examined the transcriptome differences employing RNA-seq. The study found that yak seminiferous tubules contained spermatogenic cells at all levels, while cattle-yak seminiferous tubules had reduced spermatogonia (SPG) and primary spermatocyte (Pri-SPC), fewer secondary spermatocytes (Sec-SPC), an absence of round spermatids (R-ST) and sperms (S), and possessed large vacuoles. All of these conditions could have significantly reduced the volume and weight of cattle-yak testis compared to that of yak. RNA-seq analysis identified 8473 differentially expressed genes (DEGs; 3580 upregulated and 4893 downregulated). GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment evaluations for DEGs found their relation mostly to spermatogenesis and apoptosis. Among the DEGs, spermatogonia stem cell (SSCs) marker genes (Gfra1, CD9, SOHLH1, SALL4, ID4, and FOXO1) and genes involved in apoptosis (Fas, caspase3, caspase6, caspase7, caspase8, CTSK, CTSB and CTSC) were significantly upregulated, while differentiation spermatogenic cell marker genes (Ccna1, PIWIL1, TNP1, and TXNDC2) and meiosis-related genes (TEX14, TEX15, MEIOB, STAG3 and M1AP) were significantly downregulated in cattle-yak. Furthermore, the alternative splicing events in cattle-yak were substantially decreased than in yak, suggesting that the lack of protein subtypes could be another reason for spermatogenic arrest in cattle-yak testis.

Background and aim Yak estrus is a seasonal phenomenon, probably involving epigenetic regulation of synthesis and secretion of sex hormones as well as growth and development of follicles. N6-methyladenosine (m6A) is the most common internal modification of the eukaryotic mRNA. However, there are no detailed reports on the m6A transcriptome map of yak ovary. Therefore, this study aimed to collected the yak ovarian tissues at three different states of anestrus (YO-A), estrus (YO-F), and pregnancy (YO-P), and obtained the full transcriptome m6A map in yak by MeRIP-seq. Results The HE staining revealed that the number of growing follicles and mature follicles in the ovary during the estrus period was relatively higher than those in the anestrus period and the pregnancy period. The RT-qPCR showed that the expression of METTL3, METTL14, FTO, YTHDC1 were significantly different across different periods in the ovaries, which suggests that m6A may play a regulatory role in ovarian activity. Next, we identified 20,174, 19,747 and 13,523 m6A peaks in the three ovarian samples of YO-A, YO-F and YO-P using the methylated RNA immunoprecipitation sequencing (MeRIP-seq). The m6A peaks are highly enriched in the coding sequence (CDS) region and 3′untranslated region (3′UTR) as well as the conserved sequence of “RRACH.” The GO, KEGG and GSEA analysis revealed the involvement of m6A in many physiological activities of the yak’s ovary during reproductive cycle. The association analysis found that some genes such as BNC1, HOMER1, BMP15, BMP6, GPX3, and WNT11 were related to ovarian functions. Conclusions The comparison of the distribution patterns of methylation peaks in the ovarian tissues across different periods further explored the m6A markers related to the regulation of ovarian ovulation and follicular development in the yak ovary. This comprehensive map provides a solid foundation for revealing the potential function of the mRNA m6A modification in the yak ovary.