Single-cell RNA-seq and single-cell bisulfite-sequencing reveal insights into yak preimplantation embryogenesis.

Extensive epigenetic reprogramming occurs during preimplantation embryonic development. However, the impact of DNA methylation in plateau yak preimplantation embryos and how epigenetic reprogramming contributes to transcriptional regulatory networks are unclear. In this study, we quantified gene expression and DNA methylation in oocytes and a series of yak embryos at different developmental stages and at single-cell resolution using single-cell bisulfite-sequencing and RNA-seq. We characterized embryonic genome activation and maternal transcript degradation and mapped epigenetic reprogramming events critical for embryonic development. Through cross-species transcriptome analysis, we identified 31 conserved maternal hub genes and 39 conserved zygotic hub genes, including SIN3A, PRC1, HDAC1/2, and HSPD1. Notably, by combining single-cell DNA methylation and transcriptome analysis, we identified 43 candidate methylation driver genes, such as AURKA, NUSAP1, CENPF, and PLK1, that may be associated with embryonic development. Finally, using functional approaches, we further determined that the epigenetic modifications associated with the histone deacetylases HDAC1/2 are essential for embryonic development and that the deubiquitinating enzyme USP7 may affect embryonic development by regulating DNA methylation. Our data represent an extensive resource on the transcriptional dynamics of yak embryonic development and DNA methylation remodeling, and provide new insights into strategies for the conservation of germplasm resources, as well as a better understanding of mammalian early embryonic development that can be applied to investigate the causes of early developmental disorders.

A drug-resistant mutation in plant target of rapamycin validates the specificity of ATP-competitive TOR inhibitors in vivo.

Kinases are major components of cellular signaling pathways, regulating key cellular activities through phosphorylation. Kinase inhibitors are efficient tools for studying kinase targets and functions, however assessing their kinase specificity in vivo is essential. The identification of resistant kinase mutants has been proposed to be the most convincing approach to achieve this goal. Here, we address this issue in plants via a pharmacogenetic screen for mutants resistant to the ATP-competitive TOR inhibitor AZD-8055. The eukaryotic TOR (Target of Rapamycin) kinase is emerging as a major hub controlling growth responses in plants largely thanks to the use of ATP-competitive inhibitors. We identified a dominant mutation in the DFG motif of the Arabidopsis TOR kinase domain that leads to very strong resistance to AZD-8055. This resistance was characterized by measuring root growth, photosystem II (PSII) activity in leaves and phosphorylation of YAK1 (Yet Another Kinase 1) and RPS6 (Ribosomal protein S6), a direct and an indirect target of TOR respectively. Using other ATP-competitive TOR inhibitors, we also show that the dominant mutation is particularly efficient for resistance to drugs structurally related to AZD-8055. Altogether, this proof-of-concept study demonstrates that a pharmacogenetic screen in Arabidopsis can be used to successfully identify the target of a kinase inhibitor in vivo and therefore to demonstrate inhibitor specificity. Thanks to the conservation of kinase families in eukaryotes, and the possibility of creating amino acid substitutions by genome editing, this work has great potential for extending studies on the evolution of signaling pathways in eukaryotes.

Seasonal stability of the rumen microbiome contributes to the adaptation patterns to extreme environmental conditions in grazing yak and cattle.

BACKGROUND: The rumen microbiome plays an essential role in maintaining ruminants' growth and performance even under extreme environmental conditions, however, which factors influence rumen microbiome stability when ruminants are reared in such habitats throughout the year is unclear. Hence, the rumen microbiome of yak (less domesticated) and cattle (domesticated) reared on the Qinghai-Tibetan Plateau through the year were assessed to evaluate temporal changes in their composition, function, and stability. RESULTS: Rumen fermentation characteristics and pH significantly shifted across seasons in both cattle and yak, but the patterns differed between the two ruminant species. Ruminal enzyme activity varied with season, and production of xylanase and cellulase was greater in yak compared to cattle in both fall and winter. The rumen bacterial community varied with season in both yak and cattle, with higher alpha diversity and similarity (beta diversity) in yak than cattle. The diversity indices of eukaryotic community did not change with season in both ruminant species, but higher similarity was observed in yak. In addition, the similarity of rumen microbiome functional community was higher in yak than cattle across seasons. Moreover, yak rumen microbiome encoded more genes (GH2 and GH3) related to cellulose and hemicellulose degradation compared to cattle, and a new enzyme family (GH160) gene involved in oligosaccharides was uniquely detected in yak rumen. The season affected microbiome attenuation and buffering values (stability), with higher buffering value in yak rumen microbiome than cattle. Positive correlations between antimicrobial resistance gene (dfrF) and CAZyme family (GH113) and microbiome stability were identified in yak, but such relationship was negatively correlated in cattle. CONCLUSIONS: The findings of the potential of cellulose degradation, the relationship between rumen microbial stability and the abundance of functional genes varied differently across seasons and between yak and cattle provide insight into the mechanisms that may underpin their divergent adaptation patterns to the harsh climate of the Qinghai-Tibetan Plateau. These results lay a solid foundation for developing strategies to maintain and improve rumen microbiome stability and dig out the potential candidates for manufacturing lignocellulolytic enzymes in the yak rumen to enhance ruminants' performance under extreme environmental conditions.

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.

Genetic structure analysis of yak breeds and their response to adaptive evolution.

Yaks are crucial genetic resources in the Tibetan Plateau and surrounding regions. Throughout the long process of domestication, natural and artificial selection pressures have enabled yaks to demonstrate adaptive characteristics to the environment in terms of physiological structure and genetic molecules, but no systematic cell analysis has been carried out on this phenomenon of yaks. Here, the population structure and genetic diversity of yak were studied by WGRS, and the genes related to yak adaptability were excavated. Combined with scRNA-seq method, the transcription map of yak lung tissue and skin tissue was constructed, which provided a new comprehensive insight into yak adaptability. The analysis of yak population structure showed that there was obvious genetic differentiation between TZ _ yak and other seven yak populations, while there was significant genetic exchange between PL _ yak and SB _ yak at high altitude. WGRS and scRNA-seq analysis revealed that the gene HIF1A related to high altitude adaptation was expressed in various cell types, while EPAS1 was predominantly expressed in epithelial and endothelial cells of yak lung tissue. Endothelial cells play a critical role in hypoxia-adapted VEGF signaling, which correlates closely with the high expression of KDR and VEGFA genes in endothelial cells and monocytes. Furthermore, in the selection signal of High _ yak vs Low _ yak, 19.8 % of the genes overlapped with the genes screened by skin scRNA-seq, including genes related to coat color such as RORA, BNC2, and KIT. Notably, BNC2 is a gene associated with melanin deposition and shows high expression levels in HS cells. Additionally, GRN in melanocytes and SORT1 in IRS play an important role in cell communication between melanocytes and IRS. These findings offer new insights into the natural polymorphism of yaks and provide a valuable reference for future research on high-altitude mammals, and potentially even human genetics.

A conjoint analysis of renal structure and omics characteristics reveal new insight to yak high-altitude hypoxia adaptation.

BACKGROUND: Yaks have unique adaptive mechanisms to the hypoxic environment, in which the kidney plays an important role. The aim of this study was to explore the histological changes of yak kidney at different altitudes and the metabolites and genes associated with adaptation to the hypoxic environment. METHODS: We analyzed the tissue structure and transcriptomic metabolomic data of yak kidney tissue at two altitudes, 2600 and 4400 m. We compared and identified the morphological adaptations of the kidney and the metabolites and genes associated with hypoxia adaptation in yaks. Changes in renal morphological adaptations, differential metabolites and genes were compared and identified, combining the two in a joint analysis. RESULTS: High-altitude yak kidneys showed significant adaptive changes: increased mitochondria, increased glomerular thylakoid area, and decreased localized ribosomes. Transcriptomics and metabolomics identified 69 DAMs (Differential metabolites) and 594 DEGs (differential genes). Functional enrichment analysis showed that the DAMs were associated with protein digestion and absorption, ABC transporter, and MTOR signaling pathway; the DEGs were significantly enriched in Cholesterol metabolism and P53 signaling pathway. The joint analysis indicated that metabolites such as lysine and arginine, as well as key genes such as ABCB5 and COL1A2, were particularly affected under hypoxic conditions, whereas changes in mitochondria in the tissue structure may be related to the expression of MFN1 and OPA1, and changes in glomerular thylakoid membranes are related to VEGFA and TGFB3. CONCLUSION: The kidney regulates metabolites and gene expression related to hormone synthesis, protein metabolism, and angiogenesis by adjusting the mitochondrial and glomerular thylakoid membrane structure to support the survival of yaks in high-altitude environments.

Transcriptomics and metabolomics of blood, urine and ovarian follicular fluid of yak at induced estrus stage.

To gain a deeper understanding of the metabolic differences within and outside the body, as well as changes in transcription levels following estrus in yaks, we conducted transcriptome and metabolome analyses on female yaks in both estrus and non-estrus states. The metabolome analysis identified 114, 13, and 91 distinct metabolites in urine, blood, and follicular fluid, respectively. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis highlighted an enrichment of pathways related to amino acid and lipid metabolism across all three body fluids. Our transcriptome analysis revealed 122 differentially expressed genes within microRNA (miRNA) and 640 within long non-coding RNA (lncRNA). Functional enrichment analysis of lncRNA and miRNA indicated their involvement in cell signaling, disease resistance, and immunity pathways. We constructed a regulatory network composed of 10 lncRNAs, 4 miRNAs, and 30 mRNAs, based on the targeted regulation relationships of the differentially expressed genes. In conclusion, the accumulation of metabolites such as amino acids, steroids, and organic acids, along with the expression changes of key genes like miR-129 during yak estrus, provide initial insights into the estrus mechanism in yaks.

Whole-transcriptome sequencing analysis to identify key circRNAs, miRNAs, and mRNAs in the development of yak testes.

BACKGROUND: The Testis is an important reproductive organ in male mammals and the site for spermatogenesis, androgen synthesis, and secretion. Non-coding RNAs (ncRNAs) play an important regulatory role in various biological processes. However, the regulatory role of ncRNAs in the development of yak testes and spermatogenesis remains largely unclear. RESULT: In this study, we compared the expression profiles of circular RNAs (circRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs) in yak testicular tissue samples collected at 6 months (Y6M), 18 months (Y18M), and 4 years (Y4Y). Using RNA sequencing (RNA-Seq), we observed a significant difference in the expression patterns of ncRNAs in the samples collected at different testicular development stages. Twenty-two differentially expressed (DE) circRNAs, 69 DE miRNAs, and 64 DE mRNAs were detected in Y6M, Y18M, and Y4Y testicular samples, respectively. The results of gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that the source genes of DE circRNAs, predicted target genes of DE miRNAs, and DE mRNAs were specifically associated with signaling pathways and GO terms that were related to sperm synthesis, sperm vitality, and testicular development, such as cell cycle, Wnt signaling pathway, MAPK signaling pathway, GnRH signaling pathway, and spermatogenesis. The analysis of the circRNA-miRNA-mRNA network revealed that some DE ncRNAs, including miR-574, miR-449a, CDC42, and CYP11A1, among others, may be involved in testicular spermatogenesis. Concurrently, various circRNA-miRNA interaction pairs were observed. CONCLUSION: Our findings provide a database of circRNAs, miRNAs, and mRNAs expression profiles in testicular tissue of yaks at different developmental stages and a detailed understanding of the regulatory network of ncRNAs in yak testicular development and provide data that can help elucidate the molecular mechanisms underlying yak testicular development.

Genetic diversity, phylogeography, and maternal origin of yak (Bos grunniens).

BACKGROUND: There is no consensus as to the origin of the domestic yak (Bos grunniens). Previous studies on yak mitochondria mainly focused on mitochondrial displacement loop (D-loop), a region with low phylogenetic resolution. Here, we analyzed the entire mitochondrial genomes of 509 yaks to obtain greater phylogenetic resolution and a comprehensive picture of geographical diversity. RESULTS: A total of 278 haplotypes were defined in 509 yaks from 21 yak breeds. Among them, 28 haplotypes were shared by different varieties, and 250 haplotypes were unique to specific varieties. The overall haplotype diversity and nucleotide diversity of yak were 0.979 ± 0.0039 and 0.00237 ± 0.00076, respectively. Phylogenetic tree and network analysis showed that yak had three highly differentiated genetic branches with high support rate. The differentiation time of clades I and II were about 0.4328 Ma, and the differentiation time of clades (I and II) and III were 0.5654 Ma. Yushu yak is shared by all haplogroups. Most (94.70%) of the genetic variation occurred within populations, and only 5.30% of the genetic variation occurred between populations. The classification showed that yaks and wild yaks were first clustered together, and yaks were clustered with American bison as a whole. Altitude had the highest impact on the distribution of yaks. CONCLUSIONS: Yaks have high genetic diversity and yak populations have experienced population expansion and lack obvious phylogeographic structure. During the glacial period, yaks had at least three or more glacial refugia.

Yak genome database: a multi-omics analysis platform.

BACKGROUND: The yak (Bos grunniens) is a large ruminant species that lives in high-altitude regions and exhibits excellent adaptation to the plateau environments. To further understand the genetic characteristics and adaptive mechanisms of yak, we have developed a multi-omics database of yak including genome, transcriptome, proteome, and DNA methylation data. DESCRIPTION: The Yak Genome Database ( http://yakgenomics.com/ ) integrates the research results of genome, transcriptome, proteome, and DNA methylation, and provides an integrated platform for researchers to share and exchange omics data. The database contains 26,518 genes, 62 transcriptomes, 144,309 proteome spectra, and 22,478 methylation sites of yak. The genome module provides access to yak genome sequences, gene annotations and variant information. The transcriptome module offers transcriptome data from various tissues of yak and cattle strains at different developmental stages. The proteome module presents protein profiles from diverse yak organs. Additionally, the DNA methylation module shows the DNA methylation information at each base of the whole genome. Functions of data downloading and browsing, functional gene exploration, and experimental practice were available for the database. CONCLUSION: This comprehensive database provides a valuable resource for further investigations on development, molecular mechanisms underlying high-altitude adaptation, and molecular breeding of yak.

Metabolomics and proteomics insights into subacute ruminal acidosis etiology and inhibition of proliferation of yak rumen epithelial cells in vitro.

BACKGROUND: Untargeted metabolomics and proteomics were employed to investigate the intracellular response of yak rumen epithelial cells (YRECs) to conditions mimicking subacute rumen acidosis (SARA) etiology, including exposure to short-chain fatty acids (SCFA), low pH5.5 (Acid), and lipopolysaccharide (LPS) exposure for 24 h. RESULTS: These treatments significantly altered the cellular morphology of YRECs. Metabolomic analysis identified significant perturbations with SCFA, Acid and LPS treatment affecting 259, 245 and 196 metabolites (VIP > 1, P < 0.05, and fold change (FC) ≥ 1.5 or FC ≤ 0.667). Proteomic analysis revealed that treatment with SCFA, Acid, and LPS resulted in differential expression of 1251, 1396, and 242 proteins, respectively (FC ≥ 1.2 or ≤ 0.83, P < 0.05, FDR < 1%). Treatment with SCFA induced elevated levels of metabolites involved in purine metabolism, glutathione metabolism, and arginine biosynthesis, and dysregulated proteins associated with actin cytoskeleton organization and ribosome pathways. Furthermore, SCFA reduced the number, morphology, and functionality of mitochondria, leading to oxidative damage and inhibition of cell survival. Gene expression analysis revealed a decrease the genes expression of the cytoskeleton and cell cycle, while the genes expression associated with inflammation and autophagy increased (P < 0.05). Acid exposure altered metabolites related to purine metabolism, and affected proteins associated with complement and coagulation cascades and RNA degradation. Acid also leads to mitochondrial dysfunction, alterations in mitochondrial integrity, and reduced ATP generation. It also causes actin filaments to change from filamentous to punctate, affecting cellular cytoskeletal function, and increases inflammation-related molecules, indicating the promotion of inflammatory responses and cellular damage (P < 0.05). LPS treatment induced differential expression of proteins involved in the TNF signaling pathway and cytokine-cytokine receptor interaction, accompanied by alterations in metabolites associated with arachidonic acid metabolism and MAPK signaling (P < 0.05). The inflammatory response and activation of signaling pathways induced by LPS treatment were also confirmed through protein interaction network analysis. The integrated analysis reveals co-enrichment of proteins and metabolites in cellular signaling and metabolic pathways. CONCLUSIONS: In summary, this study contributes to a comprehensive understanding of the detrimental effects of SARA-associated factors on YRECs, elucidating their molecular mechanisms and providing potential therapeutic targets for mitigating SARA.

Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differential isoform in skin tissues of different hair length Yak.

BACKGROUND: The hair follicle development process is regulated by sophisticated genes and signaling networks, and the hair grows from the hair follicle. The Tianzhu white yak population exhibits differences in hair length, especially on the forehead and shoulder region. However, the genetic mechanism is still unclear. Isoform sequencing (Iso-seq) technology with advantages in long reads sequencing. Hence, we combined the Iso-seq and RNA-seq methods to investigate the transcript complexity and difference between long-haired yak (LHY) and normal-haired yak (NHY). RESULTS: The hair length measurement result showed a significant difference between LHY and NHY on the forehead and the shoulder (P-value < 0.001). The skin samples from the forehead and the shoulder of LHY and NHY were pooled for isoform sequencing (Iso-seq). We obtained numerous long transcripts, including novel isoforms, long non-coding RNA, alternative splicing events, and alternative polyadenylation events. Combined with RNA-seq data, we performed differential isoforms (DEIs) analysis between LHY and NHY. We found that some hair follicle and skin development-related DEIs, like BMP4, KRT2, IGF2R, and COL1A2 in the forehead skin; BMP1, KRT1, FGF5, COL2A1, and IGFBP5 in the shoulder skin. Enrichment analysis revealed that DEIs in both two comparable groups significantly participated in skin and hair follicle development-related pathways, like ECM-receptor interaction, focal adhesion, and PI3K-Akt signaling pathways. The results indicated that the hair follicle development of Tianzhu white yak may influence the hair length difference. Besides, the protein-protein interaction (PPI) network of DEIs showed COL2A1 and COL3A1 exhibited a high degree of centrality, and these two genes were suggested as potential candidates for the hair length growth of Tianzhu white yak. CONCLUSIONS: The results provided a comprehensive analysis of the transcriptome complexity and identified differential transcripts that enhance our understanding of the molecular mechanisms underlying the variation in hair length growth in Tianzhu white yak.

Metagenomic analysis for exploring the potential of Lactobacillus yoelii FYL1 to mitigate bacterial diarrhea and changes in the gut microbiota of juvenile yaks.

Diarrhea in calves is a common disease that results in poor nutrient absorption, poor growth and early death which leads to productivity and economic losses. Therefore, it is important to explore the methods to reduce diarrhea in yak's calves. Efficacy of lactic acid bacteria (LAB) for improvement of bacterial diarrhea is well recognized. For this purpose, two different doses (107 CFU, 1011 CFU) of Lactobacillus yoelii FYL1 isolated from yaks were fed to juvenile yaks exposed to E. coli O78. After a trial period of ten days fresh feces and intestinal contents of the experimental yaks were collected and metagenomics sequencing was performed. It was found that feeding a high dose of Lactobacillus yoelii FYL1 decreased abundance of phylum Firmicutes in the E. coli O78 infected group whereas, it was high in animals fed low dose of Lactobacillu yoelii FYL1. Results also revealed that counts of bacteria from the family Oscillospiraceae, genus Synergistes and Megasphaera were higher in control group whereas, order Bifidobacteriales and family Bifidobacteriaceae were higher in infected group. It was observed that bacterial counts for Pseudoruminococcus were significantly (P < 0.05) higher in animals of group that were given high dose of Lactobacillus yoelii FYL1 (HLAB). Compared to infected group multiple beneficial bacterial genera such as Deinococus and Clostridium were found higher in the animals that were given a low dose of Lactobacillus yoelii FYL1 (LLAB). The abundance of pathogenic bacterial genera that included Parascardovia, Bacteroides and Methanobrevibacter was decreased (P < 0.05) in the lower dose treated group. The results of functional analysis revealed that animals of LLAB had a higher metabolism of terpenoids and polyketides compared to animals of infected group. Virus annotation also presented a significant inhibitory effect of LLAB on some viruses (P < 0.05). It was concluded that L. yoelii FYL1 had an improved effect on gut microbiota of young yaks infected with E. coli O78. This experiment contributes to establish the positive effects of LAB supplementation while treating diarrhea.

Exploring the induction of endometrial epithelial cell apoptosis in clinical-type endometritis in yaks through the cyt-c/caspase-3 signaling axis.

Endometritis is a significant contributor to reduced productivity in yaks in Tibet, China. The Cyt-c/Caspase-3 signaling axis plays a crucial role in the mitochondrial pathway that triggers cell apoptosis due to endogenous factors. In this study, we examined the endometrial epithelial tissue of yaks with endometritis using pathological examination, immunohistochemical analysis, TUNEL staining, qRT-PCR, and Western blot. The results indicated significant changes in the apoptotic factors of the Cyt-c/Caspase-3 signaling axis. The expression levels of Bak1, Bax, Cyt-c, Apaf-1, Caspase-9, and Caspase-3 were significantly increased (P < 0.05), while the expression level of Bcl-2 was significantly decreased. Immunohistochemistry results revealed significant increase in Bak1, Bax, Cyt-c, Apaf-1, Caspase-9, and Caspase-3 expression in the cytoplasm compared to the healthy group, except for Bcl-2, which showed a significant decrease. Pathological section analysis demonstrated that clinical endometritis in yaks led to structural damage, bleeding, congestion, and inflammatory cell infiltration in the endometrial epithelium. Our study findings indicated that clinical endometritis in yaks can modulate apoptosis of endometrial epithelial cells via the Cyt-c/Caspase-3 signaling pathway, resulting in different levels of damage. This research is pioneering in exploring cell apoptosis induced by clinical endometritis in yaks, offering novel insights and potential strategies for the future prevention and treatment of endometritis in yaks.

A single-cell transcriptomic atlas characterizes cell types and their molecular features in yak ovarian cortex.

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.

DEFB114 protein enhances host resistance to fungal infection through the NOD1/2-ATG16L1-NF-κB signaling pathway.

The overuse of antibiotics has led to the enhanced resistance of many pathogenic bacteria, posing a threat to human health. Therefore, there is a need to develop green and safe alternatives to antibiotics. Beta-defensins play a crucial role in host defense against pathogens and have multifunctional properties, exerting key roles in innate and adaptive immunity, as well as non-immune processes. In this study, a 210 bp long cDNA sequence of yak DEFB114 gene was amplified and successfully expressed in a prokaryotic system. The DEFB114 protein exhibited significant inhibitory effects on the growth of Aspergillus fumigatus in vitro. When co-cultured with yak macrophages, DEFB114 protein enhanced macrophage phagocytic activity and increased nucleic acid fluorescence intensity (P < 0.05). DEFB114 protein also enhanced the activity of yak macrophages stimulated by inactivated Aspergillus fumigatus spores, increased the release of nitric oxide (NO), and promoted the expression of genes such as γ-actin, Lgals, Man2b, and Capg (P < 0.05). In mice experiments, DEFB114 protein promoted resistance against Aspergillus fumigatus infection, by regulating the NOD1/2-ATG16L1-NF-κB pathway to modulate the host immune response and exert its anti-infective effects. In summary, the yak DEFB114 protein could inhibit the growth of Aspergillus fumigatus and enhance the animal's resistance to pathogenic microorganisms, thereby having significant implications in the treatment and prevention of fungal infections.

Identification of circRNA-associated ceRNA networks in the longissimus dorsi of yak under different feeding systems.

BACKGROUND: Yaks (Bos grunniens), prized for their ability to thrive in high-altitude environments, are indispensable livestock in the plateau region. Modifying their feeding systems holds significant promise for improving their growth and meat quality. Tenderness, a key determinant of yak meat quality and consumer appeal, is demonstrably influenced by dietary regimen. Indoor feeding regimes have been shown to enhance tenderness by lowering shear stress and optimizing pH values. CircRNAs, well-known modulators of circulatory function, also play a crucial role in skeletal muscle development across various animal species. However, their functional significance in yak skeletal muscle remains largely unexplored. RESULTS: In this study, we identified a total of 5,534 circRNAs within the longissimus dorsi muscle, and we found 51 differentially expressed circRNAs (20 up-regulated and 31 down-regulated) between the two feeding groups. Constructing a comprehensive ceRNA network illuminated intricate regulatory mechanisms, with PGP and circRNA_0617 converging on bta-miR-2285q, mirrored by KLF15/circRNA_0345/bta-miR-20b and CTSF/circRNA_0348/bta-miR-146a. These findings shed light on the potential of circRNAs to influence yak muscle development and meat quality, offering valuable insights for future research. CONCLUSIONS: This investigation unraveled a complex interaction network between circRNAs、mRNAs and miRNAs in yak skeletal muscle. We further elucidated the target genes regulated by these target genes within the network, offering valuable insights into the potential regulatory mechanisms governing muscle development and meat quality-related traits in yaks.

Correlation between microbial characteristics and reproductive status of the yak uterus based on macrogenomic analysis.

INTRODUCTION: This study aimed to investigate the microbial characteristics of yak uteri collected using intrauterine cotton swabs (CS) during different reproductive stages and the correlation of these microbial characteristics with reproductive status. METHODS: We used a macrogenomic approach to analyze the functional aspects of different microorganisms in samples collected during the pre-estrus, estrus, late estrus, and diestrus stages. RESULTS: The results revealed the presence of 1293 microbial genera and 3401 microbial species in the uteri of yaks at different reproductive stages. The dominant bacterial species varied across the different periods, with Micrococcus and Proteus being dominant during pre-estrus; Pseudomonas, Clostridium, Flavobacterium, Bacillus, and Staphylococcus during estrus; Acinetobacter, Bacillus and Proteus during late estrus; and Pseudomonas, Escherichia coli, and Proteus during diestrus. DISCUSSION: The primary functions of these bacteria are enriched in various metabolic pathways, including carbohydrate and amino acid metabolism, intracellular transport and secretion, post-translational protein modification, and drug resistance. These findings suggest that the microbial diversity in the uterus of yaks plays a crucial role in reproductive regulation and can help prevent reproductive tract-related diseases.

Alterations in the diversity and composition of the fecal microbiota of domestic yaks (Bos grunniens) with pasture alteration-induced diarrhea.

Diarrhea is a common issue in domestic yaks (Bos grunniens) that can occur with pasture alterations and significantly impacts growth performance. Previous research has examined the microbiota of diarrhetic yaks; however, the structural changes in gut bacterial community and microbial interactions in yaks with grassland alteration-induced diarrhea remain poorly understood. To explore variations in gut microbiota homeostasis among yaks suffering from diarrhea, fecal microbiota diversity and composition were analyzed using 16 S rRNA amplicon sequencing. Gut fecal microbiota diversity was lower in diarrhetic yaks than in non-diarrhetic yaks. Furthermore, the bacterial community composition (including that of Proteobacteria and Actinobacteria) in the feces of diarrhetic yaks displayed significant alterations. Co-occurrence network analysis further underscored the compromised intestinal flora stability in yaks with diarrhea relative to that in non-diarrhetic yaks. Interestingly, the abundance of beneficial bacteria, such as Lachnospiraceae_AC2044_group and Lachnospiraceae_NK4A136_group, were decreased in yaks with diarrhea, and the reductions were negatively correlated with the fecal water content. Collectively, these findings indicate that diminished microbial stability and increased abundance of certain bacteria in the gut may contribute to diarrhea occurrence in yaks.

Whole genome resequencing revealed genomic variants and functional pathways related to adaptation in Indian yak populations.

The present study aims to identify genomic variants through a whole genome sequencing (WGS) approach and uncover biological pathways associated with adaptation and fitness in Indian yak populations. A total of 30 samples (10 from each population) were included from Arunachali, Himachali and Ladakhi yak populations. WGS analysis revealed a total of 32171644, 27260825, and 32632460 SNPs and 4865254, 4429941, and 4847513 Indels in the Arunachali, Himachali, and Ladakhi yaks, respectively. Genes such as RYR2, SYNE2, BOLA, HF1, and the novel transcript ENSBGRG00000011079 were found to have the maximum number of high impact variants in all three yak populations, and might play a major role in local adaptation. Functional enrichment analysis of genes harboring high impact SNPs revealed overrepresented pathways related to response to stress, immune system regulation, and high-altitude adaptation. This study provides comprehensive information about genomic variants and their annotation in Indian yak populations, thus would serve as a data resource for researchers working on the yaks. Furthermore, it could be well exploited for better yak conservation strategies by estimating population genetics parameters viz., effective population size, inbreeding, and observed and expected heterozygosity.