Skeletal muscle proteome analysis provides insights on high altitude adaptation of yaks

Skeletal muscle transcriptomics of sheep acclimated to cold desert and tropical regions identifies genes and pathways accentuating their diversity

The current study attempts to investigate the differences in gene expression in longissimus thoracis muscles between sheep breeds acclimated to diverse environments. Changthangi sheep inhabits the cold arid plateau of Ladakh, at an altitude above 3000 m with prevalence of rarefied atmosphere. Muzzafarnagri sheep, on the other hand is found in the sub-tropical hot and humid plains at an altitude of about 250 m. Comparative transcriptomics was used to provide a molecular perspective of the differential adaptation of the two breeds. RNA sequencing data was generated from four biological replicates of the longissimus thoracis muscles from both breeds. The common genes expressed in both breeds were involved in muscle contraction and muscle fibre organization. The most significant pathways enriched in Changthangi muscles were glycogen metabolism, reduction of cytosolic Ca++ levels and NFE2L2 regulating anti-oxidant, while those in Muzzafarnagri were extracellular matrix organization and collagen formation. The hub genes identified in Changthangi were involved in hematopoiesis and HIF signaling pathway, suggesting the molecular acclimatization of Changthangi to the high altitude cold desert of Ladakh. The nodal genes discovered in Muzzafarnagri sheep were associated with the extracellular matrix which accentuates its significance in the development, growth and repair of muscles. The observed transcriptomic differences underscore the morphological and adaptive disparity between the two breeds. The candidate genes and pathways identified in this study will form the basis for future research on adaptation to high altitude and body size in small ruminants.

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

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

Quantitative N-Glycoproteomic Analysis of Cattle-Yak and Yak Longissimus Thoracis

In this study, the N-glycosylated protein profiles of cattle-yak longissimus thoracis (CYLT) and yak longissimus thoracis (YLT) were comparatively analyzed using quantitative proteomics techniques. A total of 76 differential N-glycosylated proteins (DGPs) were screened from 181 quantified N-glycoproteins, indicating that differences in N-glycosylation levels are key to the differences between CYLT and YLT. In particular, a variety of N-glycoproteins involved in the extracellular matrix were differentially N-glycosylated between CYLT and YLT, mainly including fibrillin-1, fibromodulin, collagen, and laminins. In addition, the N-glycosylation levels of several lysosomal-related proteolytic enzymes (cathepsin D, dipeptidyl peptidase 1, legumain, and aminopeptidases, etc.) were significantly higher in CYLT. These results indicated that the N-glycosylation of CYLT and YLT proteins plays a crucial role in the regulation of extracellular matrix organization (muscle fiber structure) and lysosomal activity (postmortem meat tenderness). The results remind us that posttranslation modifications, especially N-glycosylation, are still icebergs beneath the surface.

Whole transcriptome analyses and comparison reveal the metabolic differences between oxidative and glycolytic skeletal muscles of yak

Mammalian skeletal muscle is composed of various muscle fibers that exhibit different physiological and metabolic features. Muscle fiber type composition has significant influences on the meat quality of livestock. In this study, we comprehensively analyzed the whole transcriptome profiles of the oxidative muscle biceps femoris (BF) and the glycolytic muscle obliquus externus abdominis (OEA) of yak. A total of 1436 mRNAs, 1172 lncRNAs, and 218 circRNAs were differentially expressed in the oxidative muscles compared with the glycolytic muscles. KEGG annotation showed that differentially expressed mRNAs regulated by lncRNA and circRNA were mainly involved in PPAR signaling pathway, citrate cycle (TCA cycle), and PI3K-Akt signaling pathway, which reflect the different metabolic properties between oxidative and glycolytic muscles. In addition, regulatory networks associated with muscle fiber type conversion and mitochondria energy metabolism in muscles were constructed. Our study provides new evidence for a better understanding of the molecular mechanisms underlying skeletal muscle fiber determination and meat quality traits of yak.

Comparative proteomic analysis of spleen reveals key immune-related proteins in the yak (Bos grunniens) at different growth stages

Spleen plays an indispensable role in the immune system as the largest lymphatic organ in the body. The spleens of yaks at three developmental stages (1 day fetal yak, 15 months juvenile yak and 5 years old adult yak) were sampled and the Tandem mass tag (TMT) quantification method was employed in spleen proteomic analysis. The results showed that 6576 proteins and 529 differentially expressed proteins (DEPs) were identified in the yak spleens at three growth stages. Gene ontology (GO) analysis of DEPs indicated that DEPs were enriched in Oxygen transport, Actin filament movement, DNA replication, Cell cycle process, and Cell macromolecule biosynthesis process, which was conducive to high altitude breathing, protein synthesis and organ growth in yaks. These were indispensable for yak spleen growth and cell metabolism, high altitude adaptation. Those DEPs were further analyzed based on Kyoto encyclopedia of genes and genomes (KEGG) pathways, which principally participated in Th1 and Th2 cell differentiation, NF-kappa B signaling pathway, Phagosome, and Glutathione metabolism. Those pathways were associated with some animal life activities in defense against microbial antigens, indicating that with age, the immune function of the yak's spleen continued to increase. Hemoglobin, Tumor necrosis factor receptor associated factor 1 (TRAF1), T cell receptor (TCR), Macrophage receptor, Fc receptors (FcR), and Gamma-glutamyl transferase (GGT) of DEPs played roles in immune function in yak spleen directly or indirectly. The dynamic changes of Toll like receptor 2 (TLR2), TRAF1 and Heat shock protein 27 (HSP27 or HSPB1) detected by Immunohistochemistry were consistent with those obtained from TMT proteomic. In conclusion, this study provides extensive and functional analyses of the spleen proteome at three developmental stages and will offer a new insight into key proteins involved in the immune function of yak spleen.

Integrated proteomic, phosphoproteomic, and N-glycoproteomic analyses of the longissimus thoracis of yaks

Yaks (Bos mutus) live in the Qinghai–Tibet plateau. The quality of yak meat is unique due to its genetic and physiological characteristics. Identification of the proteome of yak muscle could help to reveal its meat-quality properties. The common proteome, phosphoproteome, and N-glycoproteome of yak longissimus thoracis (YLT) were analyzed by liquid chromatography-tandem mass spectrometry-based shotgun analysis. A total of 1812 common proteins, 1303 phosphoproteins (3918 phosphorylation sites), and 204 N-glycoproteins (285 N-glycosylation sites) were identified in YLT. The common proteins in YLT were involved mainly in myofibril structure and energy metabolism; phosphoproteins were associated primarily with myofibril organization, regulation of energy metabolism, and signaling; N-glycoproteins were engaged mainly in extracellular-matrix organization, cellular immunity, and organismal homeostasis. We reported, for the first time, the “panorama” of the YLT proteome, specifically the N-glycoproteome of YLT. Our results provide essential information for understanding post mortem physiology (rigor mortis and aging) and the quality of yak meat.

•

A total of 2650 proteins were identified in yak longissimus thoracis.

•

Common proteins were involved mainly in myofibril structure and energy metabolism.

•

Phosphoproteins were associated with myofibrils, energy metabolism, and signaling.

•

N-glycoproteins were engaged mainly in ECM organization, immunity, and homeostasis.

Differences in proteomic profiles between yak and three cattle strains provide insights into molecular mechanisms underlying high-altitude adaptation

Yaks display unique properties of the lung and heart, enabling their adaptation to high-altitude environments, but the underlying molecular mechanisms are still largely unknown. In the present study, the proteome differences in lung and heart tissues were compared between yak (Bos grunniens) and three cattle strains (Bos taurus, Holstein, Sanjiang and Tibetan cattle) using the sequential window acquisition of all theoretical mass spectra/data-independent acquisition (SWATH/DIA) proteomic method. In total, 51,755 peptides and 7215 proteins were identified. In the lung tissue, there were 162, 310 and 118 differential abundance proteins (DAPs) in Tibetan, Holstein and Sanjiang cattle compared to yak respectively. In the heart tissue, there were 71, 57 and 78 DAPs in Tibetan, Holstein and Sanjiang cattle compared to yak respectively. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the DAPs were enriched for the retinol metabolism and toll-like receptor categories in lung tissue. The changes in these two pathways may regulate hypoxia-induced factor and immune function in yaks. Moreover, DAPs in heart tissues were enriched for cardiac muscle contraction, Huntington's disease, chemical carcinogenesis and drug metabolism-cytochrome P450. Further exploration indicated that yaks may alter cardiac function through regulation of type 2 ryanodine receptor (RyR2) and Ca²⁺ -release channels. The present results are useful to further develop an understanding of the mechanisms underlying adaptation of animals to high-altitude conditions.

Signature of high altitude adaptation in the gluteus proteome of the yak

Yak is the unique Bovidae species in the Qinghai-Tibetan Plateau. A previous proteomic study has compared the yak muscle tissue to one cattle strain using the isobaric tags for relative and absolute quantification approach. In this study, to further investigate the molecular mechanisms underlying yak adaptation, the proteomic profiles of gluteus were compared between yak and one moderate-altitude cattle strain (Tibetan cattle) and two low-altitude cattle strains (Holstein and Sanjiang cattle) using a label-free quantitative method. The comparisons identified 20, 364, 143 upregulated proteins and 4, 6, 37 downregulated proteins in yak, compared with Tibetan, Holstein, and Sanjiang cattle, respectively. Protein-protein interaction analysis indicated that these differentially expressed proteins were mainly related to "oxidative phosphorylation" and "electron transport chain." Further analysis revealed that NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 11, NADH dehydrogenase [ubiquinone] 1 beta subcomplex subunit 4, cytochrome C oxidase subunit 6A2, mitochondrial and cytochrome c oxidase subunit NDUFA4 were all increased in the yak, suggesting that yak might increase mitochondrial capacity to sustain metabolic rates under high altitude conditions, which might be a long-term adaptive mechanism underlying adaptation to high altitude environments. Yak increased the level of thioredoxin reductase 2 to protect themselves from oxidative damages. Moreover, the increased expression levels of phosphatidylinositol 4,5-bisphosphate 3-kinase catalytic subunit alpha isoform and caveolin-1 in yak suggested that yaks promoted glucose uptake for adaptation to high altitude. These results provided more information to better understand the molecular mechanisms underlying yak adaption.