Search

Abstract In recent years, the rapid proliferation of digital services and resources on Industrial Internet has imposed higher demands on universality and privacy of identity management. Particularly with the advent of the digital economy, prudent users are urged to maintain control over their digital identity credentials. However, traditional identity management methods have failed to meet this requirement and thus have been prone to raise users' concerns about potential financial loss. Specifically, conventional identity management systems(IDMS) have been plagued by imperceptible privacy disclosure, which derives from the flaws in single points of failure, excessive disclosure, correlation analysis, traceability, and revocation. The emerging Self-Sovereign Identity (SSI) architecture aims to tackle these issues and is propelling the evolution of privacy-enhanced distributed identity management. To this end, we proposed a privacy-enhanced distributed identity management scheme with sequential aggregate issuance, threshold traceability and revocability in the setting of multiple credential issuers and regulators. We adopted the Decentralized identifiers(DIDs) and verifiable credentials(VCs) based on the SSI architecture to ensure the hierarchical identity authentication. The security and performance analysis shows that our proposal achieves the desired design goals and is feasible for distributed Industrial Internet.

Abstract 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.

Milk lipids greatly affect the volatile flavor of milk, and the relationship between lipids and volatile flavor in yak milk was explored in this study. The volatile flavor compounds (VFCs), lipids profile, fatty acids in yak ordinary milk and colostrum were detected with HP/SPME-GC–MS, the semiquantitative lipidomics based on LC-MS/MS, GC–MS, respectively. The VFCs differences in yak milk were closely related to 1-((1 s,3ar,4r,7 s,7as)-4-hydroxy-7-isopropyl-4-methyloctahloctahydro-1h-inden-1-Yl)-ethanone,2,6,6-trimethyl-2,4-cycloheptadien-1-one, pentanal, 2-phenylethyl propionate, octanoic acid methyl ester, diphosphoric acid diisooctyl ester, (Z)-3,4,4-trimethyl-5-oxo-2-hexenoic acid and acetic acid. The volatile flavor in yak milk was well correlated with milk lipids, and TG(4:0_12:3_18:1), TG(6:0_8:0_18:1), TG(4:0_12:3_18:1), TG(12:0_18:2_18:3) and TG(16:0e_18:1_22:5) were the crucial lipid molecules affecting volatile flavor. The degeneration of above lipids by hydrolysis produced some fatty acids and alcohol, then these compounds were further derived into other VFCs especially above crucial 8 molecules. This study provided a theoretical basis for improving the volatile flavor by controlling lipids in yak milk.

Lysozyme like 4 (LYZL4), lysozyme like 6 (LYZL6) and proliferating cell nuclear antigen (PCNA) are implicated in the regulation of testicular function, but there was no research reported available on the expression patterns of LYZL4, LYZL6 and PCNA genes at different developmental stages of yak testes. In this study, we used the qRT-PCR, western blotting and immunohistochemistry estimated the LYZL4, LYZL6 and PCNA gene expression and protein lo-calization at different developmental stages of yak testes. The qPCR results showed that the mRNA expression of LYZL4, LYZL6 and PCNA genes significantly increased with age in the testes of yaks. Western blot results showed that the protein abundance of LYZL4, LYZL6 and PCNA in yak testes was significantly higher after puberty than before puberty. Furthermore, the results of immunohistochemistry indicated that LYZL4, LYZL6 and PCNA may be involved in the regulation of spermatogonia proliferation and Leydig cell function in immature testis. In adult yak testes, LYZL4, LYZL6 and PCNA may involve in the development of round spermatids and primary spermatocytes during testicular development. Our results indicated that LYZL4, LYZL6 and PCNA may be involved in the development of Sertoli cells, Leydig cells and gonocytes in yak testes.

Abstract 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

Abstract 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.

Abstract Pathogenic gut microbiota is responsible for a few debilitating gastrointestinal diseases. While the host immune cells do produce extracellular vesicles to counteract some deleterious effects of the microbiota, the extracellular vesicles are of insufficient doses and at unreliable exposure times. Here we use mechanical stimulation of hydrogel-embedded macrophage in a bioelectronic controller that on demand boost production of up to 20 times of therapeutic extracellular vesicles to ameliorate the microbes’ deleterious effects in vivo. Our miniaturized wireless bioelectronic system termed inducible mechanical activation for in-situ and sustainable generating extracellular vesicles (iMASSAGE), leverages on wireless electronics and responsive hydrogel to impose mechanical forces on macrophages to produce extracellular vesicles that rectify gut microbiome dysbiosis and ameliorate colitis. This in vivo controllable extracellular vesicles-produced system holds promise as platform to treat various other diseases.

Background Despite the high risk of eating disorder (ED)-related attitudes and behaviors among female dancers, targeted scientific dietary regimens are currently inadequate. Time-restricted eating (TRE), a popular intermittent fasting protocol, has been shown to be effective in enhancing body composition and exercise performance in athletes. In this study, TRE was employed as a dietary regimen to improve body composition and exercise performance and address ED attitudes and behaviors in DanceSport dancers.Methods Twenty female DanceSport dancers were recruited and divided into two groups: TRE (n = 10) and normal diet (ND) (n = 10). The TRE group consumed their self-selected necessary energy intake exclusively between 11 a.m. and 7 p.m. (utilizing a 16-hour fasting and 8-hour eating window) for 6 weeks, while the ND group maintained their regular dieting patterns. The consumption of water, black tea, or coffee without added sugar or milk was not restricted. Physical activity and calorie intake were systematically recorded during the TRE intervention. Body composition, aerobic and anaerobic performance, and ED attitudes and behaviors were assessed before and after the TRE intervention. The trial was registered in the Chinese Clinical Trial Registry under the identifier ChiCTR2200063780.Results The fixed effects tests (p 0.05). TRE resulted in a reduction in hip circumference (p = 0.039), fat mass (kg) (p = 0.0004), and body fat percentage (p = 0.0005), with no significant decrease in fat-free mass (p > 0.05). No significant improvement was observed in aerobic performance (p > 0.05). The average power (AP) (p = 0.01) and AP/Body weight ratio (p = 0.003) significantly increased. Additionally, the power drop decreased significantly (p = 0.019). Group-by-time interactions were observed for fat mass (kg) (p = 0.01), body fat percentage (p = 0.035), and AP/Body weight (p = 0.020).Conclusion TRE can be considered a feasible nutritional strategy for DanceSport dancers, facilitating improvements in body composition without compromising aerobic and anaerobic exercise performance or exacerbating ED attitudes and behaviors. Moreover, TRE may facilitate more favorable physiological adaptations, potentially contributing to improved exercise performance.

Yaks (Bos grunniens) exhibit exceptional adaptation to the challenging high-altitude environment of the Qinghai-Tibetan plateau, making them the sole bovine species capable of thriving in such exreme conditions. Investigating the cellular and molecular characteristics of yak ovaries across different reproductive states is crucial for gaining insight into their ovarian functions. Herein, the cellular atlases of yak ovaries in different reproductive states were depicted by single-cell RNA-sequencing (scRNA-seq). The cellular atlases of the ovaries were established by identifying specific gene expression patterns of various cell types, including granulosa cells, theca cells, stromal cells, smooth muscle cells, endothelial cells, glial cell, macrophages, natural killer cells, and proliferating cells. The cellular compositions of the ovaries vary among different reproductive states. Furthermore, the granulosa cells comprise six cell subtypes, while theca cells consist of eight cell subtypes. The granulosa cells and theca cells exhibit distinct biological functions throughout different reproductive states. The two cell types were aligned along their respective pseudotime trajectories. Moreover, a cell-to-cell communication network was constructed among distinct cell types within the ovary, spanning the three reproductive states. Notably, during the estrus period, the granulosa cells demonstrated more prominent interactions with other cell types compared to the remaining reproductive states.

Abstract 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.

Abstract Industrial data trading can considerably enhance the economic and social value of abundant data resources. However, traditional data trading models are plagued by critical flaws in fairness, security, privacy and regulation. To tackle the above issues, we first proposed a distributed industrial data trading architecture based on blockchain and cloud for multiple data owners. Subsequently, we realized implemented distributed identity management by the distributed verifiable credentials scheme that possesses the desirable properties, i.e., selective disclosure, multi-show unlinkability, threshold traceability, and public verifiability. Finally, we presented a fair trading mechanism without trusted third parties based on smart contracts, and we employed blockchain and multi-signature to ensure data integrity during data storage and trading. The security and performance analysis shows that our proposal is feasible for sensitive data trading for multiple data owners and provides a useful exploration for future industrial data trading and management.

Coronary restenosis is an important cause of poor long-term prognosis in patients with coronary heart disease. Here, we show that lysine methyltransferase SMYD2 expression in the nucleus is significantly elevated in serum- and PDGF-BB-induced vascular smooth muscle cells (VSMCs), and in tissues of carotid artery injury-induced neointimal hyperplasia. Smyd2 overexpression in VSMCs (Smyd2-vTg) facilitates, but treatment with its specific inhibitor LLY-507 or SMYD2 knockdown significantly inhibits VSMC phenotypic switching and carotid artery injury-induced neointima formation in mice. Transcriptome sequencing revealed that SMYD2 knockdown represses the expression of serum response factor (SRF) target genes and that SRF overexpression largely reverses the inhibitory effect of SMYD2 knockdown on VSMC proliferation. HDAC3 directly interacts with and deacetylates SRF, which enhances SRF transcriptional activity in VSMCs. Moreover, SMYD2 promotes HDAC3 expression via tri-methylation of H3K36 at its promoter. RGFP966, a specific inhibitor of HDAC3, not only counteracts the pro-proliferation effect of SMYD2 overexpression on VSMCs, but also inhibits carotid artery injury-induced neointima formation in mice. HDAC3 partially abolishes the inhibitory effect of SMYD2 knockdown on VSMC proliferation in a deacetylase activity-dependent manner. Our results reveal that the SMYD2-HDAC3-SRF axis constitutes a novel and critical epigenetic mechanism that regulates VSMC phenotypic switching and neointimal hyperplasia.

This study explored the polymorphism of the leucine aminopeptidase (LAP3) gene and its relationship with milk quality characteristics in Gannan yak. A cohort of 162 Gannan yak was genotyped utilizing the Illumina Yak cGPS 7K BeadChip, and the identified single nucleotide polymorphisms (SNPs) were evaluated for their association with milk protein, casein, lactose, and fat concentrations. The results showed that four SNPs (g.4494G > A, g.5919A > G, g.8033G > C, and g.15,615A > G) in the LAP3 gene exhibited polymorphism with information content values of 0.267, 0.267, 0.293, and 0.114, respectively. All four SNPs were in Hardy–Weinberg equilibrium (p > 0.05). The g.4494G > A and g.5919A > G SNPs were significantly associated with protein content (p < 0.05), with homozygous genotypes showing significantly higher protein content than heterozygous genotypes (p < 0.05). The g.8033G > C SNP was significantly associated with casein content, protein content, non-fat solids, and acidity (p < 0.05), with the CC genotype having significantly higher casein, protein, and non-fat solids content than the GG and GC genotypes (p < 0.05). The g.15,615A > G SNP was significantly associated with average fat globule diameter (p < 0.05). In general, the mutations within the LAP3 gene demonstrated a positive impact on milk quality traits in Gannan yak, with mutated genotypes correlating with enhanced milk quality. These results indicate that the LAP3 gene could be a significant or candidate gene affecting milk quality traits in Gannan yak and offer potential genetic markers for molecular breeding programs in this species.

Abstract Yak has been subject to natural selection, human domestication and interspecific introgression during its evolution. However, genetic variants favored by each of these processes have not been distinguished previously. We constructed a graph-genome for 47 genomes of 7 cross-fertile bovine species. This allowed detection of 57,432 high-resolution structural variants (SVs) within and across the species, which were genotyped in 386 individuals. We distinguished the evolutionary origins of diverse SVs in domestic yaks by phylogenetic analyses. We further identified 334 genes overlapping with SVs in domestic yaks that bore potential signals of selection from wild yaks, plus an additional 686 genes introgressed from cattle. Nearly 90% of the domestic yaks were introgressed by cattle. Introgression of an SV spanning the KIT gene triggered the breeding of white domestic yaks. We validated a significant association of the selected stratified SVs with gene expression, which contributes to phenotypic variations. Our results highlight that SVs of different origins contribute to the phenotypic diversity of domestic yaks.

Energy feed can provide animals with balanced nutrition, thereby enhancing their growth performance. This study aimed to evaluate the effects of dietary energy levels on the growth performance, serum metabolites, and meat quality of Jersey cattle–yaks. A total of 24 male Jersey cattle–yaks were randomly divided into three groups. Each group was fed diets with metabolizable energy levels of 8.21 MJ/kg (LE), 9.50 MJ/kg (ME), and 10.65 MJ/kg (HE), respectively. The HE and ME groups showed significantly higher final body weight, average daily gain (ADG), and feed efficiency compared to the LE group (p < 0.05). The glucose (GLU) and total cholesterol (TC) concentrations were significantly increased in the serum of the ME and HE groups (p < 0.05). The low-density lipoprotein cholesterol (LDL-C) and alanine aminotransferase (ALT) levels were significantly higher in the serum of the HE group than in the ME group (p < 0.05). Blood urea nitrogen (BUN) levels exhibited a significant decrease with increasing metabolizable energy levels in the diet (p < 0.05). Increasing dietary energy levels enhances the eye muscle area and intramuscular fat content of Jersey cattle–yaks (p < 0.05), with no effect on pH45 min, pH24 h, and shear force. In the HE group, the levels of heneicosanoic acid (C21:0), palmitoleic acid (C16:1), elaidic acid (C18:1n9t), and eicosadienoic acid (C20:2n6) were notably elevated (p < 0.05) when compared to the LE group. We concluded that a higher dietary energy level enhanced the growth performance and meat quality traits of male Jersey cattle–yaks.

Male reproductive health is largely determined already in the early development of the testis. Although much work has been carried out to study the mechanisms of testicular development and spermatogenesis, there was previously no information on the differences in the protein composition of yak testicles during early development. In this study, the protein profiles in the testicles of 6- (M6), 18- (M18), and 30-month-old (M30) yaks were comparatively analyzed using TMT proteomics. A total of 5521 proteins were identified, with 13, 1295, and 1397 differentially expressed proteins (DEPs) in 30- vs. 18-, 18- vs. 6-, and 30- vs. 6-month-old testes, respectively. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that DEPs were mainly involved in signaling pathways related to testicular development and spermatogenesis, including the MAPK, PI3K–Akt, Wnt, mTOR, TGF-β, and AMPK signaling pathways. Furthermore, we also identified eight potential proteins (TEX101, PDCL2, SYCP2, SYCP3, COL1A1, COL1A2, ADAM10, and ATF1) that may be related to the testicular development and spermatogenesis of yaks. This study may provide new insights into the molecular mechanisms of the testicular development and spermatogenesis of yaks.

In order to investigate the composition and differences in volatile organic compounds (VOCs) in yak and cattle-yak meat and determine the key metabolites and metabolic pathways related to flavor formation. In this study, the VOCs and non-volatile metabolites in Longissimus dorsi muscle of two groups of samples were detected and analyzed by gas chromatography–ion migration spectrometry (GC-IMS) and gas chromatography–mass spectrometry (GC-MS). The results showed that 31 VOCs were identified by GC-IMS, including 5 alcohols, 5 ketones, 5 esters, 3 aldehydes, 2 furans, 2 hydrocarbons, 1 amine, 1 acid, 1 thiazole, 1 pyrazine, and 5 others. Most of them were alcohols, ketones, esters, and aldehydes. A total of 75 non-volatile metabolites with significant differences were obtained by GC-MS screening, among which amino acid contents such as serine, glycine, phenylalanine, and aspartic acid were significantly up-regulated in cattle-yak, and glutamic acid and tyrosine were significantly up-regulated in yak. The non-volatile differential metabolites in the two groups were significantly enriched in the metabolic pathways of arginine biosynthesis and oxidative phosphorylation. By combining GC-IMS and GC-MS, this study comprehensively and intuitively reflected the differences in VOCs between yak and cattle-yak meat, and clarified the metabolomic reasons for the differences in VOCs, so as to provide a theoretical basis for meat quality improvement.

Yak (Bos grunniens) is a unique large ruminant species in the Qinghai-Tibetan Plateau (QTP). Changing the energy levels of their rations can significantly improve their growth performance. Therefore, studying the effects of dietary energy levels on the rumen microflora and metabolites of yak is crucial for enhancing the development of the yak industry. Currently, there is a lack of understanding regarding the impact of feeding energy diets on rumen fermentation parameters, microbial functions, and metabolites. This study was designed to determine the appropriate energy level for feeding yak. Three test diets with metabolizable energy levels of 7.57 MJ/kg, 9.44 MJ/kg, and 11.9 MJ/kg were used and the concentration of volatile fatty acids (VFA) in rumen fluid was measured. The microbial communities, functions, and metabolites in yaks were studied by 16S rRNA sequencing, metagenome, and LC-MS non-targeted metabolomics to investigate the relationships among rumen fermentation parameters, microbial diversity, and metabolites. Ration energy levels significantly affect total VFA, acetate, propionate, butyrate, iso-valerate, valerate, and acetate/propionate (p < 0.05). At the phylum level, the dominant phyla in all three treatment groups were Bacteroidota, Firmicutes, and Actinobacteriota. At the genus level, the abundance of the unclassified_o__Bacteroidales, norank_f_Muribaculaceae, Lachnospiraceae_NK4A136_group, and Family _XIII_AD3011_group showed significant differences (p < 0.05) and were significantly correlated with differential metabolites screened for phosphatidylcholine [PC(16:0/0:0), PC(18:3/0:0)], uridine 3′-monophosphate, and adenosine monophosphate, etc. CAZymes family analysis showed that GHs and CEs differed significantly among the three groups. In addition, differential metabolites were mainly enriched in the pathways of lipid metabolism, nucleotide metabolism, and biosynthesis of other secondary metabolites, and the concentrations of differential metabolites were correlated with microbial abundance. In summary, this study analyzed the effects of ration energy levels on rumen microorganisms and metabolites of yaks and their relationships. The results provided a scientific basis for the selection of dietary energy for yaks in the house feeding period in the future.

The growth and development of skeletal muscle also determine the meat production of yak, ultimately affecting the economic benefits. Hence, improving growth performance is a top priority in the yak industry. Skeletal muscle development is a complex process involving the regulation of several genes, including microRNAs (miRNAs). However, the transcription of miRNAs in yak skeletal muscle during prenatal to postnatal stages is unknown. We used small RNA sequencing (small RNA-Seq) to determine the global miRNAs of longissimus dorsi muscle from yak (the samples were collected from three fetuses and three adults). Totally 264 differently expressed miRNAs (|log2(fold change)|>1 and P-value≤0.05) were detected between the two groups. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that differently expressed miRNAs-targeted genes participated in pathways associated with muscle development, such as MAPK, PI3K-Akt, and Hippo signaling pathways, etc. MiR-652, which was up-regulated in the fetal group, was transfected into C2C12 myoblasts to examine its role. miR-652 promoted (P≤0.05) proliferation and differentiation, but inhibited (P≤0.001) apoptosis at early period. Furthermore, miR-652 reduced (P≤0.001) the proportion of C2C12 myoblasts in the G1 phase while increasing (P≤0.01) the proportion of cells in the S and G2 phases. Dual-luciferase reporter assays indicated that ISL1 served as a target of miR-652. In general, these findings expand our understanding of yak skeletal muscle miRNAs, and suggested that miR-652 probably regulated myogenesis by regulating ISL1.

A single-nucleotide polymorphism (SNP) is a genome-level trait that arises from a variation in a single nucleotide, leading to diversity in DNA sequences. SNP screening is commonly used to provide candidate genes for yak breeding efforts. Integrin Subunit Alpha 9 (ITGA9) is an integrin protein. It plays an important role in cell adhesion, signalling, and other processes. The aim of this study was to discuss the association between genetic polymorphisms in the ITGA9 gene and milk quality traits and to identify potential molecular marker loci for yak breeding quality. We genotyped 162 yaks using an Illumina Yak cGPS 7K liquid chip and identified the presence of polymorphisms at nine SNP loci in the ITGA9 gene of yaks. The results showed that the mutant genotypes in the loci g.285,808T>A, g.306,600T>C, and g.315,413C>T were positively correlated with the contents of casein, protein, total solids (TS), and solid nonfat (SNF) in yak milk. In other loci, heterozygous genotypes had a positive correlation with nutrient content in yak milk. Then, two ITGA9 haplotype blocks were constructed based on linkage disequilibrium, which facilitated a more accurate screening of ITGA9 as a candidate gene for yak milk quality improvement. In conclusion, we identified SNPs and haplotype blocks related to yak milk quality traits and provided genetic resources for marker-assisted selection in yak breeding.

The yak is a unique species of livestock found in the Qinghai-Tibet Plateau and its surrounding areas. Due to factors such as late sexual maturity and a low rate of estrus, its reproductive efficiency is relatively low. The process of estrus synchronization in yaks plays a crucial role in enhancing their reproductive success and ensuring the continuation of their species. In order to clarify the characteristics of the serum metabolites of yak estrus synchronization, the yaks with inactive ovaries were compared with the estrus synchronization yaks. In this study, yaks were divided into the inactive ovaries group (IO), gonarelin-induced yak estrus group (GnRH), and chloprostenol sodium-induced yak estrus group (PGF). After the completion of the estrus synchronization treatment, blood samples were collected from the jugular veins of the non-estrus yaks in the control group and the yaks with obvious estrus characteristics in the GnRH and PGF groups. Metabolites were detected by ultra-high performance liquid chromatography-mass spectrometry, and differential metabolites were screened by multivariate statistical analysis. The results showed that a total of 70 significant differential metabolites were screened and identified in the GnRH vs. IO group, and 77 significant differential metabolites were screened and identified in the PGF vs. IO group. Compared with non-estrus yaks, 36 common significant differential metabolites were screened out after the induction of yak estrus by gonarelin (GnRH) and cloprostenol sodium (PGF), which were significantly enriched in signaling pathways such as the beta oxidation of very long chain fatty acids, bile acid biosynthesis, oxidation of branched chain fatty acids, steroidogenesis, steroid biosynthesis, and arginine and proline metabolism. This study analyzed the effects of gonadotropin releasing hormone (GnRH) and prostaglandin F (PGF) on the reproductive performance of yaks treated with estrus synchronization, which provides a theoretical basis for the optimization and application of yak estrus synchronization technology and promotes the healthy development of the yak industry.

Protein kinase D1 (PRKD1) functions primarily in normal mammary cells, and the potassium voltage-gated channel subfamily Q member 3 (KCNQ3) gene plays an important role in controlling membrane potential and neuronal excitability, it has been found that this particular gene is linked to the percentage of milk fat in dairy cows. The purpose of this study was to investigate the relationship between nucleotide polymorphisms (SNPs) of PRKD1 and KCNQ3 genes and the milk quality of Gannan yak and to find molecular marker sites that may be used for milk quality breeding of Gannan yak. Three new SNPs were detected in the PRKD1 (g.283,619T>C, g.283,659C>A) and KCNQ3 gene (g.133,741T>C) of 172 Gannan lactating female yaks by Illumina yak cGPS 7K liquid-phase microarray technology. Milk composition was analyzed using a MilkoScanTM milk composition analyzer. We found that the mutations of these three loci significantly improved the lactose, milk fat, casein, protein, non-fat milk solid (SNF) content and acidity of Gannan yaks. The lactose content of the TC heterozygous genotype population at g.283,619T>C locus was significantly higher than that of the TT wild-type population (p < 0.05); the milk fat content of the CA heterozygous genotype population at g.283,659C>A locus was significantly higher than that of the CC wild-type and AA mutant populations (p < 0.05); the casein, protein and acidity of the CC mutant and TC heterozygous groups at the g.133,741T>C locus were significantly higher than those of the wild type (p < 0.05), and the SNF of the TC heterozygous group was significantly higher than that of the mutant group (p < 0.05). The results showed that PRKD1 and KCNQ3 genes could be used as candidate genes affecting the milk traits of Gannan yak.

Software-defined networking (SDN) faces many of the same security threats as traditional networks. The separation of the SDN control plane and data plane makes the controller more vulnerable to cyber attacks. The conventional “perimeter defense” network security model cannot prevent lateral movement attacks caused by malicious insider users or hardware and software vulnerabilities. The “zero trust architecture” has become a new security network model to protect enterprise network security. In this article, we propose an intelligent zero-trust security framework IZTSDN for the software-defined networking by integrating deep learning and zero-trust architecture, which adopts zero-trust architecture to protect every resource and network connection in the network. IZTSDN uses a traffic anomaly detection mode CALSeq2Seql based on a deep learning algorithm to analyze users’ network behavior in real-time and achieve continuous tracking and analysis of users, restrict malicious users from accessing network resources, and realize the dynamic authorization process. Finally, the Mininet simulation platform is extended to build the simulation platform MiniIZTA supporting zero-trust architecture and the proposed security framework IZTSDN is experimentally analyzed. The experimental results show that the IZTSDN security framework can provide about 80.5% of throughput when the network is attacked. The accuracy of abnormal traffic detection reaches 99.56% on the SDN dataset, which verifies that the reliability and availability of the IZTSDN security framework are verified.

Abstract Background Normal testicular development is highly crucial for male reproduction and is a precondition for spermatogenesis that is the production of spermatozoa in the testes. MiRNAs have been implicated in several testicular biological processes, including cell proliferation, spermatogenesis, hormone secretion, metabolism and reproductive regulation. In the present study, we used deep sequencing data to study the functions of miRNAs during testicular development and spermatogenesis, by analyzing the expression patterns of small RNAs in 6-, 18- and 30-month-old yak testis tissues. Results A total of 737 known and 359 novel miRNAs were obtained from 6-, 18- and 30-month-old yak testes. In all, we obtained 12, 142 and 139 differentially expressed (DE) miRNAs in 30- vs. 18-, 18- vs. 6-, and 30- vs. 6-month-old testes, respectively. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of all DE miRNA target genes revealed BMP2, TGFB2, GDF6, SMAD6, TGFBR2 and other target genes as participants in different biological processes, including TGF-β, GnRH, Wnt, PI3K–Akt, MAPK signaling pathways and several other reproductive pathways. In addition, quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) was used to detect the expression of seven randomly selected miRNAs in 6-, 18- and 30-month-old testes, and the results were consistent with the sequencing data. Conclusions The differential expression of miRNAs in yak testes at different development stages was characterized and investigated using deep sequencing technology. We believe that the results will contribute to further understanding the functions of miRNAs in regulating the development of yak testes and improving the reproductive performance of male yaks.

Abstract Background As the direct organ of villus, hair follicles have obvious seasonal cycles. The hair follicle cycle is orchestrated by multiple cell types that together direct cell renewal and differentiation. But the regulation property of hair follicle cells from anagen to catagen in yak is still unknown. Results In this study, single-cell RNA sequencing was performed on 24,124 single cells of the scapular skin from white yak. Based on tSNE cluster analysis, the cell types of IFE-DC, epidermal cell lines, fibroblasts, keratinocytes, IRS, DS, INFU, and other cells in yak hair follicles during anagen and catagen were successfully identified, and the gene expression profiles were described. The GO enrichment analysis indicated the different cells characteristic genes to be mainly enriched in the epidermal development, epithelial cell differentiation and wound healing pathways. The pseudotime trajectory analysis described the differentiation trajectory of the epidermal lineage and dermal lineage of the hair follicle during anagen and catagen. Moreover, the dynamic changes of the genes like LHX2, KRT25, and KRT71 were found to be highly expressed in HS and IRS, but not in the IFE-DC, INFU, and keratinocyte during differentiation. Conclusions Our results analyzed the time-varying process of gene expression in the dermal cell lineage and epidermal cell lineage of hair follicles during anagen and catagen during fate differentiation was expounded at the single cell level, revealing the law of fate specialization of different types of cells. In addition, based on the enrichment analysis, the transcriptional regulatory factors involved in the different cell fates were also revealed. These results will help to enhance our understanding of yak hair follicle cycle and promote the development and utilization of yak villus.

Abstract Background Physical activity at pre-older ages (55–64 years) can greatly affect one’s physical fitness, health, physical-activity behaviour, and quality of life at older ages. The objective of this study was to conduct a 24-week walking-exercise programme among sedentary pre-older females and investigate the influence of different walking cadences on cardiorespiratory fitness and associated biomarkers. Methods A total of 78 pre-older sedentary female participants were recruited and randomly assigned to normal (n = 36), paced (n = 15), music-synchronised (n = 15) walking, and no-exercise control (n = 12) groups, respectively. The normal, paced, and music-synchronised walking groups walked at a cadence of 120 steps/min, 125 steps/min, and 120–128 steps/min, respectively, under supervised conditions. Anthropometric characteristics, step length, nutrient intake, blood pressure and composition, and cardiorespiratory fitness were measured at baseline, the 12th week of the programme, the 24th week of the programme (completion), and after a 12-week retention period, which began immediately upon completion of the programme and did not feature any supervised exercises. Results All walking conditions improved high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol, step length, maximum oxygen consumption (VO2max), and oxidative capacity at anaerobic threshold (all P

Abstract Background The hair coat is available for the yak to live in the harsh environment of the plateau. Besides, improving the hair production of yak is necessary for its textile industry development. Hair grows from hair follicles (HFs). The HFs undergo periodic growth after birth and are regulated by the complex gene regulatory network. However, the molecular mechanism of HFs regeneration in the Tianzhu white yak remains unclear. RNA editing is a post-transcriptional mechanism that regulates gene expression and produces new transcripts. Hence, we investigated the influence of the A-to-I RNA editing events on the HFs cycle of the Tianzhu white yak. Results We finally identified 54,707 adenosine-to-inosine (A-to-I) RNA editing sites (RESs) from RNA sequencing data of the HFs cycle in the Tianzhu white yak. Annotation results showed RESs caused missense amino acid changes in 7 known genes. And 202 A-to-I editing sites altered 23 target genes of 140 microRNAs. A total of 1,722 differential RESs were identified during the HFs cycle of Tianzhu white yak. GO and KEGG enrichment analysis revealed several signaling pathways and GO terms involved skin development, hair growth, and HFs cycle. Such as genes with differential RNA editing levels were significantly enriched in the peroxisome, metabolic pathways, Notch signaling pathway, and PPAR signaling pathway. Besides, the editing sites in HFs development-related genes FAS, APCDD1, WWOX, MPZL3, RUNX1, KANK2, DCN, DSC2, LEPR, HEPHL1, and PTK2B were suggested as the potential RESs involving HFs development. Conclusion This study investigated the global A-to-I RNA editing events during the HFs cycle of yak skin tissue and expanded the knowledge of A-to-I RNA editing on the HFs cycle. Furthermore, this study revealed that RNA editing-influenced genes may regulate the HFs cycle by participating in the HFs development-related pathways. The findings might provide new insight into the regulation of RNA editing in hair growth.

The application of the valuable natural product thaxtomin A, a potent bioherbicide from the potato scab pathogenic Streptomyces strains, has been greatly hindered by the low yields from its native producers. Here, we developed an orthogonal transcription system, leveraging extra-cytoplasmic function (ECF) sigma (σ) factor 17 (ECF17) and its cognate promoter Pecf17, to express the thaxtomin gene cluster and improve the production of thaxtomin A. The minimal Pecf17 promoter was determined, and a Pecf17 promoter library with a wide range of strengths was constructed. Furthermore, a cumate inducible system was developed for precise temporal control of the ECF17 transcription system in S. venezuelae ISP5230. Theoretically, the switchable ECF17 transcription system could reduce the unwanted influences from host and alleviate the burdens introduced by overexpression of heterologous genes. The yield of thaxtomin A was significantly improved to 202.1 ± 15.3 μ g/mL using the switchable ECF17 transcription system for heterologous expression of the thaxtomin gene cluster in S. venezuelae ISP5230. Besides, the applicability of this transcription system was also tested in Streptomyces albus J1074, and the titer of thaxtomin A was raised to as high as 239.3 ± 30.6 μg/mL. Therefore, the inducible ECF17 transcription system could serve as a complement of the generally used transcription systems based on strong native constitutive promoters and housekeeping σ factors for the heterologous expression of valuable products in diverse Streptomyces hosts.

A Chinese Compilation of Physical Activities was compiled to estimate the energy costs of physical activities (PAs) using data on adults aged 18–64. Data were obtained from published articles and laboratory measurements. Databases, including PubMed, Embase, Scopus, Ebsco, Web of Science, Chinese National Knowledge Infrastructure, Wan Fang Data, National Science and Technology Report Service, Public Health Scientific Data were searched to collect data from inception to January 2022, on energy expenditure associated with PA in the healthy Chinese population. Two reviewers independently screened the literature and extracted, classified, and summarized data. Data were measured for 36 PAs using indirect calorimetry. Detailed descriptions of specific activities and metabolic equivalent values were provided by summarizing 241 physical activities in 13 categories. The first edition of the Chinese Compilation of PAs in Healthy Adults Aged 18–64(CCPA) was created. It provides valuable resources for people who regularly engage in physical exercise, researchers, educators, fitness professionals, and health or commercial sectors to quickly obtain various PA MET intensities. In the future, the energy expenditure of various PAs of different ages within the Chinese population can be measured based on the CCPA.

Yak meat is nutritionally superior to beef cattle but has a low fat content and is slow-growing. The liver plays a crucial role in lipid metabolism, and in order to determine whether different feeding modes affect lipid metabolism in yaks and how it is regulated, we employed RNA sequencing (RNA-seq) technology to analyze the genome-wide differential gene expression in the liver of yaks maintained under different raising systems. A total of 1663 differentially expressed genes (DEGs) were identified (|log2FC| ≥ 0 and p-value ≤ 0.05), including 698 down-regulated and 965 up-regulated genes. According to gene ontology (GO) and KEGG enrichment analyses, these DEGs were significantly enriched in 13 GO terms and 26 pathways (p < 0.05). Some DEGs were enriched in fatty acid degradation, PPAR, PI3K-Akt, and ECM receptor pathways, which are associated with lipid metabolism. A total of 16 genes are well known to be related to lipid metabolism (e.g., APOA1, FABP1, EHHADH, FADS2, SLC27A5, ACADM, CPT1B, ACOX2, HMGCS2, PLIN5, ACAA1, IGF1, FGFR4, ALDH9A1, ECHS1, LAMA2). A total of 11 of the above genes were significantly enriched in the PPAR signaling pathway. The reliability of the transcriptomic data was verified using qRT-PCR. Our findings provide new insights into the mechanisms regulating yak meat quality. It shows that fattening improves the expression of genes that regulate lipid deposition in yaks and enhances meat quality. This finding will contribute to a better understanding of the various factors that determine yak meat quality and help develop strategies to improve yield and quality.

Jersey-yak is a hybrid offspring of Jersey cattle and yak (Bos grunniens). Changing the feeding system of Jersey-yak can significantly improve its growth performance. In this study, tandem mass tag (TMT) proteomics technology was used to determine the differentially expressed proteins (DEPs) of the longissimus lumborum (LL) muscle of Jersey-yak fed different protein levels of diet. The results showed that compared with the traditional grazing feeding, the growth performance of Jersey-yaks was significantly improved by crude protein supplementation after grazing. A total of 3368 proteins were detected in these muscle samples, of which 3365 were quantified. A total of 434 DEPs were identified. Through analyses, it was found that some pathways related to muscle growth and development were significantly enriched, such as Rap1 signaling pathway, mTOR signaling pathway, and TGF-beta signaling pathway. A number of DEPs enriched in these pathways are related to muscle cell development, differentiation, and muscle development, including integrin subunit alpha 7 (ITGA7), myosin heavy chain 8 (MYH8), and collagen type XII alpha 1 chain (COL12A1). In conclusion, the results of this study provide insights into the proteomics of different feeding patterns of Jersey-yak, providing a stronger basis for further understanding the biological mechanism of hybrid varieties.

The morphogenesis of hair follicle structure is accompanied by the differentiation of skin tissue. Mammalian coats are produced by hair follicles. The formation of hair follicles requires signal transmission between the epidermis and dermis. However, knowledge of the transcriptional regulatory mechanism is still lacking. We used single-cell RNA sequencing to obtain 26,573 single cells from the scapular skin of yaks at hair follicle telogen and anagen stages. With the help of known reference marker genes, 11 main cell types were identified. In addition, we further analyzed the DP cell and dermal fibroblast lineages, drew a single-cell map of the DP cell and dermal fibroblast lineages, and elaborated the key genes, signals, and functions involved in cell fate decision making. The results of this study provide a very valuable resource for the analysis of the heterogeneity of DP cells and dermal fibroblasts in the skin and provide a powerful theoretical reference for further exploring the diversity of hair follicle cell types and hair follicle morphogenesis.

Abstract Background High spontaneous miscarriage rate in yak, especially during late pregnancy, have caused a great economic loss to herdsmen living in the Qinghai-Tibet plateau. However, the mechanism underlying spontaneous miscarriage is still poorly understood. In the present study, placenta protein markers were identified to elucidate the pathological reasons for yak spontaneous miscarriage through isobaric tags for relative and absolute quantification (iTRAQ) proteomic technology and bioinformatic approaches. Results Subsequently, a total of 415 differentially expressed proteins (DEPs) were identified between aborted and normal placentas. The up-regulated DEPs in the aborted placentas were significantly associated with “spinocerebellar ataxia”, “sphingolipid signalling”, “relaxin signalling”, “protein export”, “protein digestion and absorption” and “aldosterone synthesis and secretion” pathway. While the down-regulated DEPs in the aborted placentas mainly participated in “valine, leucine and isoleucine degradation”, “PPAR signalling”, “peroxisome”, “oxidative phosphorylation”, “galactose metabolism”, “fatty acid degradation”, “cysteine and methionine metabolism” and “citrate cycle” pathway. Conclusions The results implied that the identified DEPs could be considered as placental protein markers for yak miscarriage during late pregnancy, and biomacromolecule metabolic abnormality and oxidative damage might be responsible for the high spontaneous miscarriage rate in yak. These findings provide an important theoretical basis for deciphering the pathologic mechanism of late spontaneous miscarriage in yak.

IntroductionThe crude protein level in the diet will affect the fermentation parameters, microflora, and metabolites in the rumen of ruminants. It is of great significance to study the effect of crude protein levels in supplementary diet on microbial community and metabolites for improving animal growth performance. At present, the effects of crude protein level in supplementary diet on rumen fermentation parameters, microbial community, and metabolites of Jersey-Yak (JY) are still unclear.MethodsThe purpose of this experiment was to study the appropriate crude protein level in the diet of JY. The rumen fermentation indexes (volatile fatty acids and pH) were determined by supplementary diets with crude protein levels of 15.16 and 17.90%, respectively, and the microbial community and metabolites of JYs were analyzed by non-target metabonomics and metagenome sequencing technology, and the changes of rumen fermentation parameters, microbial flora, and metabolites in the three groups and their interactions were studied.Results and DiscussionThe crude protein level in the supplementary diet had significant effects on pH, valeric acid, and the ratio of acetic acid to propionic acid (p < 0.05). The protein level had no significant effect on the dominant microflora at the phylum level (p > 0.05), and all three groups were Bacteroides and Firmicutes. The results of metabolite analysis showed that the crude protein level of supplementary diet significantly affected the metabolic pathways such as Bile secretion and styrene degradation (p < 0.05), and there were different metabolites between the LP group and HP group, and these different metabolites were related to the dominant microbial to some extent. To sum up, in this experiment, the effects of crude protein level in supplementary diet on rumen microorganisms and metabolites of JY and their relationship were studied, which provided the theoretical basis for formulating a more scientific and reasonable supplementary diet in the future.

Abstract 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.

In recent years, many companies have chosen to outsource data and other data computation tasks to cloud service providers to reduce costs and increase efficiency. However, there are risks of security and privacy breaches when users outsource data to a cloud environment. Many researchers have proposed schemes based on cryptographic primitives to address these risks under the assumption that the cloud is a semi-honest participant and query users are honest participants. However, in a real-world environment, users’ data privacy and security may be threatened by the presence of malicious participants. Therefore, a novel scheme based on secure multi-party computation is proposed when attackers gain control over both the cloud and a query user in the paper. We prove that our solution can satisfy our goals of security and privacy protection. In addition, our experimental results based on simulated data show feasibility and reliability.

Coiled-coil serine-rich protein 1 (CCSER 1) gene is a regulatory protein gene. This gene has been reported to be associated with various economic traits in large mammals in recent years. The aim of this study was to investigate the association between CCSER1 gene single nucleotide polymorphisms (SNPs) and Gannan yaks and to identify potential molecular marker loci for breeding milk quality in Gannan yaks. We genotyped 172 Gannan yaks using Illumina Yak cGPS 7K liquid microarrays and analyzed the correlation between the three SNPs loci of the CCSER1 gene and the milk qualities of Gannan yaks, including milk fat, protein and casein. It was found that mutations at the g.183,843A>G, g.222,717C>G and g.388,723G>T loci all affected the fat, protein, casein and lactose traits of Gannan yak milk to varying extents, and that the milk quality of individuals with mutant phenotypes was significantly improved. Among them, the milk fat content of AG heterozygous genotype population at g.183,843A>G locus was significantly higher than that of AA and GG genotype populations (p < 0.05); the casein and protein content of mutant GG and CG genotype populations at g.222,717C>G locus was significantly higher than that of wild-type CC genotype population (p < 0.05); and the g.388,723G>T locus of the casein and protein contents of the mutant TT genotype population were significantly higher (p < 0.05) than those of the wild-type GG genotype population. These results provide potential molecular marker sites for Gannan yak breeding.

Studying the mechanism of spermatogenesis is key to exploring the reproductive characteristics of male yaks. Although N6-methyladenosine (m6A) RNA modification has been reported to regulate spermatogenesis and reproductive function in mammals, the molecular mechanism of m6A in yak testis development and spermatogenesis remains largely unknown. Therefore, we collected testicular tissue from juvenile and adult yaks and found that the m6A level significantly increased after sexual maturity in yaks. In MeRIP-seq, 1702 hypermethylated peaks and 724 hypomethylated peaks were identified. The hypermethylated differentially methylated RNAs (DMRs) (CIB2, AK1, FOXJ2, PKDREJ, SLC9A3, and TOPAZ1) mainly regulated spermatogenesis. Functional enrichment analysis showed that DMRs were significantly enriched in the adherens junction, gap junction, and Wnt, PI3K, and mTOR signaling pathways, regulating cell development, spermatogenesis, and testicular endocrine function. The functional analysis of differentially expressed genes showed that they were involved in the biological processes of mitosis, meiosis, and flagellated sperm motility during the sexual maturity of yak testis. We also screened the key regulatory factors of testis development and spermatogenesis by combined analysis, which included BRCA1, CREBBP, STAT3, and SMAD4. This study indexed the m6A characteristics of yak testicles at different developmental stages, providing basic data for further research of m6A modification regulating yak testicular development.

Abstract Background Ovarian follicle fluid (FF) as a microenvironment surrounding oocyte plays critical roles in physio-biochemical processes of follicle development and oocyte maturation. It is hypothesized that proteins in yak FF participate in the physio-biochemical pathways. The primary aims of this study were to find differentially expressed proteins (DEPs) between mature and immature FF, and to elucidating functions of the mature and immature FF in yak. Results The mature and immature FF samples were obtained from three healthy yaks that were nonpregnant, aged from four to five years, and free from any anatomical reproductive disorders. The FF samples were subjected to mass spectrometry with the isobaric tags for relative and absolute quantification (iTRAQ). The FF samples went through correlation analysis, principle component analysis, and expression pattern analysis based on quantification of the identified proteins. Four hundred sixty-three DEPs between mature and immature FF were identified. The DEPs between the mature and immature FF samples underwent gene ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG), and protein-protein interaction (PPI) analysis. The DEPs highly expressed in the mature FF mainly took parts in the complement and coagulation cascades, defense response, acute-phase response, response to other organism pathways to avoid invasion of exogenous microorganisms. The complement activation pathway contains eight DEPs, namely C2, C5, C6, C7, C9, C4BPA, CFH, and MBL2. The three DEPs, CATHL4, CHGA, and PGLYRP1, take parts in defense response pathway to prevent invasion of exogenetic microorganism. The coagulation cascades pathway involves many coagulation factors, such as F7, F13A1, FGA, FGB, FGG, KLKB1, KNG1, MASP1, SERPINA1, and SERPIND1. While the DEPs highly expressed in the immature FF participated in protein translation, peptide biosynthetic process, DNA conformation change, and DNA geometric change pathways to facilitate follicle development. The translation pathway contains many ribosomal proteins, such as RPL3, RPL5, RPS3, RPS6, and other translation factors, such as EIF3J, EIF4G2, ETF1, MOV10, and NARS. The DNA conformation change and DNA geometric change involve nine DEPs, DDX1, G3BP1, HMGB1, HMGB2, HMGB3, MCM3, MCM5, MCM6, and RUVBL2. Furthermore, the expressed levels of the main DEPs, C2 and SERPIND1, were confirmed by western blot. Conclusions The differential proteomics revealed the up-regulated DEPs in mature FF take parts in immunoreaction to prevent invasion of microorganisms and the up-regulated DEPs in immature FF participate in protein synthesis, which may improve our knowledge of the follicular microenvironment and its biological roles for reproductive processes in yak. The DEPs, C2 and SERPIND1, can be considered as protein markers for mature yak follicle.

Background: Qaidam cattle are local breeds that habitats in northwest China. It has many excellent characteristics, such as high cold and roughage tolerance, low oxygen adaptability, and tender meat quality. Copy number variation (CNV) can induce phenotypic changes in animals by a variety of effects, and thus affects the biological functions of the animals. To explore the molecular mechanism of its adaptation to extreme cold weather and muscle fat development, the CNV variations in the genome of three Qaidam cattle were detected by whole-genome sequencing, in this study. Results: A total of 16,743 CNVs and 9498 copy number variable regions (CNVRs) were obtained after the screening, which accounts for 2.18% of the bovine genome. The CNVR length detected ranged from 0.3 KB to 10.77 KB, with a total length of 58.17 MB and an average length of 6.12 KB/ CNVR. Through functional enrichment of CNVR related genes, LDHB, and ME1 genes were screened as the key genes for Qaidam cattle to adapt to the cold and low oxygen environments, whereas KIT and FGF18 genes might be related to the coat color and growth. In the CNVR overlapped with QTLs, variation in CAPN1 and CAST genes might be closely related to the tender meat quality of Qaidam cattle. Conclusions: Therefore, this study provides new genetic insights on the environmental adaptability and important economic traits of Qaidam cattle.How to cite: Guo S, Wu X, Pei J, et al. Genome-wide CNV analysis reveals variants associated with high-altitude adaptation and meat traits in Qaidam cattle. Electron J Biotechnol 2021;54. https://doi.org/10.1016/j.ejbt.2021.07.006

Researchers have linked microgravity in space to the significant imbalance between bone formation and bone resorption that induces persistent bone loss in load-bearing bones. However, the underlying molecular mechanisms are still unclear, which hinders the development of therapeutic measures. The aim of this study was to identify hub genes and explore novel molecular mechanisms underlying microgravity-induced bone loss using transcriptome datasets obtained from the GEO and SRA databases. In summary, comparative RNA expression pattern studies that differ in species (Homo or Mus), models (in vitro or in vivo), microgravity conditions (real microgravity or ground-based simulators) and microgravity duration showed that it is difficult to reach a consistent conclusion about the pathogenesis of microgravity-induced bone loss across these studies. Even so, we identified 11 hub genes and some miRNA-mRNA interactions mainly based on the GSE100930 dataset. Also, the expression of CCL2, ICAM1, IGF1, miR-101-3p and miR-451a markedly changed under clinorotation-microgravity condition. Remarkedly, ICAM1 and miR-451a were key mediators of the osteogenesis of hMSCs under clinorotation-microgravity condition. These findings provide novel insights into the molecular mechanisms of bone loss during microgravity and could indicate potential targets for further countermeasures against this condition.

Mammalian spermatozoa are not mature after ejaculation and must undergo additional functional and structural changes within female reproductive tracts to achieve subsequent fertilization, including both capacitation and acrosome reaction (AR), which are dominated by post-translational modifications (PTMs), especially phosphorylation. However, the mechanism of protein phosphorylation during frozen-thawed sperm capacitation and AR has not been well studied. In this study, the phosphoproteomics approach was employed based on tandem mass tag (TMT) labeling combined with liquid chromatography-tandem mass spectrometry (LC-MS/MS) strategy to analyze frozen-thawed sperm in Ashidan yak under three sequential conditions (density gradient centrifugation-based purification, incubation in the capacitation medium and induction of AR processes by the calcium ionophore A23187 treatment). The identification of 1,377 proteins with 5,509 phosphorylation sites revealed changes in phosphorylation levels of sperm-specific proteins involved in regulation of spermatogenesis, sperm motility, energy metabolism, cilium movement, capacitation and AR. Some phosphorylated proteins, such as AKAP3, AKAP4, SPA17, PDMD11, CABYR, PRKAR1A, and PRKAR2A were found to regulate yak sperm capacitation and AR though the cAMP/PKA signaling pathway cascades. Notably, the phosphorylation level of SPA17 at Y156 increased in capacitated sperm, suggesting that it is also a novel functional protein besides AKAPs during sperm capacitation. Furthermore, the results of this study suggested that the phosphorylation of PRKAR1A and PRKAR2A, and the dephosphorylation of CABYR both play key regulatory role in yak sperm AR process. Protein-protein interaction analysis revealed that differentially phosphorylated proteins (AKAP3, AKAP4, FSIP2, PSMD11, CABYR, and TPPP2) related to capacitation and AR process played a key role in protein kinase A binding, sperm motility, reproductive process, cytoskeleton and sperm flagella function. Taken together, these data provide not only a solid foundation for further exploring phosphoproteome of sperm in yak, but an efficient way to identify sperm fertility-related marker phosphorylated proteins.

Tibet Gaoshan Yak (Bos grunniens) is a local Yak breed that mainly produces meat in Tibet Autonomous Region, China. In this study, the complete mitochondrial genome of Tibet Gaoshan Yak was sequenced. The total length of the mitochondrial genome is 16,323 bp, and the base composition is 33.71% for A, 13.21% for G, 27.27% for T, and 25.81% for C. The genome includes 13 protein-coding genes (ND1-ND6, ND4L, COX1-COX3, ATP6, ATP8, and CYTB), two rRNA genes (12S rRNA and 16SrRNA), 22 tRNA genes, and a noncoding control region (D-loop). Phylogenetic analysis showed that Tibet Gaoshan Yak has the closest relationship with Polled Yak. The sequence analysis provided in this study will be helpful to the management of Yak breeds, the origin, and evolution of Yak, and the protection and utilization of genetic resources.

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.

Skeletal muscle growth and development is a complicated process that is regulated at multiple steps and by numerous myogenesis genes. RNA editing represents one of the events at the post-transcriptional level, which contributes to the diversity of transcriptome and proteome by altering the nucleotides of RNAs. However, RNA editing events in the skeletal muscle of yaks are still not well defined. This study conducted whole-genome RNA-editing identification in skeletal muscle of yaks at embryonic stage (ES) and adult stage (AS). We found a total of 11,168 unique RNA editing sites, most of which were detected in the intergenic region. After annotation, we totally identified 2,718 editing sites within coding regions, among which 858 were missense changes. Moreover, totally 322 editing sites in the 3′ untranslated regions (UTR) were also predicted to alter the set of miRNA target sites, indicating that RNA editing may be involved in translational repression or mRNA degradation. We found 838 RNA editing sites (involving 244 common genes) that are edited differentially in ES as compared to AS. According to the KEGG enrichment analysis, these differentially edited genes were mainly involved in pathways highly related to skeletal muscle development and myogenesis, including MAPK, AMPK, Wnt, and PI3K-Akt signaling pathways. Altogether, our work presents the first characterization of RNA editing sites within yak skeletal muscles on a genome-wide scale and enhances our understanding of the mechanism of skeletal muscle development and myogenesis.

Abstract 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.

Background: Transmembrane protein 95 (TMEM95) plays a role in male fertility. Previous studies showed that genes with a significant impact on reproductive traits can also affect the growth traits of livestock. Thus, we speculated that the genetic variation of TMEM95 gene may have effects on growth traits of cattle. Results: Two SNPs were genotyped. The rs136174626 and rs41904693 were in the intron 4 and 3′-untranslated region, respectively. The linkage disequilibrium analysis illustrated that these two loci were not linked. The rs136174626 was associated with six growth traits of Nanyang cattle, four traits of Luxi cattle, and three traits of Ji’an cattle. For rs41904693 locus, the GG individuals had greater body height and abdominal girth in Ji’ an cattle than TT and TG individuals. In Jinnan cattle, GG and TT individuals had greater body height, height at hip cross, body length, and heart girth than TG individuals. The potential splice site prediction results suggest that the rs136174626 may influence the splicing efficiency of TMEM95, and the miRNA binding site prediction results showed that the rs41904693 may influence the expression of TMEM95 by affecting the binding efficiency of Bta-miR-1584 and TMEM95 3′-UTR. Conclusions: The findings of the study suggested that the two SNPs in TMEM95 could be a reliable basis for molecular breeding in cattle.How to cite: Guo X, Zhang S, Yang H, et al. Bovine TMEM95 gene: Polymorphisms detecting in five Chinese indigenous cattle breeds and their association with growth traits. Electron J Biotechnol 2021;51. https://doi.org/10.1016/j.ejbt.2021.03.004

Abstract Although persistent efforts have identified and characterized a few candidate genes and related biological processes with potential functions in the adaptation of many species to extreme environments, few works have been conducted to determine the genomic basis of adaptation in endangered livestock breeds that have been living in extreme conditions for more than thousands of years. To fill this gap, we sequenced the whole genomes of nine individuals from three Chinese native endangered cattle breeds that are living in high‐altitude or arid environments. Phylogenetic and evolutionary history analyses of these three and other six breeds showed that the genetic structure of the cattle populations is primarily related to geographic location. Interestingly, we identified pervasive introgression from the yak to Zhangmu cattle (ZMC) that cover several genes (e.g., NOS2, EGLN1 and EPAS1) involved in the hypoxia response and previously identified as positive selection genes in other species, which suggested that the adaptive introgression from yak may have contributed to the adaptation of ZMC to high‐altitude environments. In addition, by contrasting the breeds in opposite living conditions, we revealed a set of candidate genes with various functions from hypoxia response, water metabolism, immune response and body shape change to embryo development and skeletal system development, etc., that may be related to high‐altitude or arid adaptation. Our research provides new insights into the recovery and adaptation of endangered native cattle and other species in extreme environments and valuable resources for future research on animal husbandry to cope with climate change.

This study aimed to investigate the nutritional properties of yak milk in various areas of Gannan. The milk composition analyzer, automatic amino acid analyzer, and flavor analyzer were used to detect the conventional nutrients, amino acids, and volatile flavor substances of 249 yak milks in Meiren grassland, Xiahe grassland, and Maqu grassland (hereinafter referred to as Meiren yak, Xiahe yak, and Maqu yak) in the Gannan area. The results showed that the fat content of Meiren yak milk was significantly higher than that of Maqu yak and Xiahe yak (p < 0.05). The protein content of Meiren yak milk was significantly higher than that of Xiahe yak (p < 0.05), but not significantly different from that of Maqu yak (p > 0.05). The casein content in the milk of Maqu yak was significantly higher than that of Meiren yak and Xiahe yak (p < 0.05). There was no significant difference in the lactose content of yak milk in the three regions (p > 0.05). The content of glutamic acid in the milk of Meiren yak, Xiahe yak, and Maqu yak was noticeably high, which was 1.03 g/100 g, 1.07 g/100 g, and 1.10 g/100 g, respectively. The total amino acid (TAA) content was 4.78 g/100 g, 4.87 g/100 g, and 5.0 g/100 g, respectively. The ratios of essential amino acids (EAA) and total amino acids (TAA) in the milk of Meiren yak, Xiahe yak, and Maqu yak were 42.26%, 41.27%, and 41.39%, respectively, and the ratios of essential amino acids (EAA) and nonessential amino acids (NEAA) were 73.19%, 70.28%, and 70.61%, respectively. In the yak milk samples collected from three different regions, a total of 34 volatile flavor compounds were detected, including 10 aldehydes, five esters, six ketones, four alcohols, two acids, and seven others. The main flavor substances qualitatively obtained from Meiren yak milk were ethyl acetate, n-valeraldehyde, acetic acid, heptanal, and n-hexanal. Xiahe yak milk mainly contains ethyl acetate, isoamyl alcohol, n-valeraldehyde, heptanal, and ethyl butyrate. Maqu yak milk mainly contains ethyl acetate, n-valeraldehyde, isoamyl alcohol, heptanal, ethyl butyrate, and n-hexanal. Principal component analysis showed that the flavor difference between Xiahe yak and Maqu yak was small, while the flavor difference between Xiahe yak, Maqu yak, and Meiren yak was large. The findings of this research can serve as a foundation for the future advancement and application of yak milk.